Boundary force

Liquid between the surface of two solid bodies gives rise to boundary forces. A pressure difference arises and is known as the capillary pressure (Pc). This can be calculated from Laplace s equation. 

It is noted that the virtual body force Fp depends not only on the unsteady fluid velocity, but also on the velocity and location of the particle surface, which is also a function of time. There are several ways to specify this boundary force, such as the feedback forcing scheme (Goldstein et al., 1993) and direct forcing scheme (Fadlun et al., 2000). In 3-D simulation, the direct forcing scheme can give higher stability and efficiency of calculation. In this scheme, the discretized momentum equation for the computational volume on the boundary is given as... 

Since the computational grids are generally not coincident with the location of the particle surface, a velocity interpolation procedure needs to be carried out in order to calculate the boundary force and apply this force to the control volumes close to the immersed particle surface (Fadlun et al., 2000). 

Which of the various immersed or embedded boundary methods is best— generally or for a particular case—is still an open question. Thornock and Smith (2005) introduced a Cell Adjusted Boundary Force Method for a stirred vessel. All methods proposed so far have their own pros and cons. Immersed boundary methods are also exploited in LB techniques (e.g., Derksen and Van den Akker, 1999). Rohde (2004) investigated the use of triangular facets for representing a spherical particle. 

From the strict physical point of view, the forces at the boundary of y acting on the material system it contains are due to a second system the material boundary. Let H be the Hamiltonian of the system within y, not including boundary forces, and H that of the material boundary, with a similar exclusion. Together they form a joint system with a Hamiltonian H + H + W where W is the mutual interaction potential between the two component systems. We may write W = W(X, X, t), where X represents the totality of spacial coordinates of the first system, and X , that of the second. Clearly W is the instrumentality of the action of the boundary on our system within y. 

Thus the boundary forces involve a random function (or random choice among an infinitude of functions) a stochastic process. 

We are thus led to the fundamental question How can general results regarding the macroscopic motion be obtained in spite of such a gulf in the knowledge of the boundary forces ... 

Now if the boundary forces are thought of as drawn at random from a mixture of a proportion 9 of w(x, t), 9 of w x, t), etc., simple probability reasoning shows that, with regard to this body of knowledge, the actual... 

In this section, the case of a semiinfinite solid with a concentrated force acting on the boundary is introduced. This case was originally solved by Boussinesq (1885). It should be noted that the only difference between this case and the case of a point force in an infinite solid medium is the boundary conditions. Shear stresses vanish on the boundary of the semiinfinite solid. In the following, the concept of a center of compression is introduced. The stress field in a semiinfinite solid with a boundary force can be obtained by superimposing the stress fields from a point force and a series of centers of compression. A center of compression is defined as the combination of three perpendicular pair forces. 

Figure 6. Hypothetical coupling of length scales in the attachment of two hydroxylated oxide nanoparticles in aqueous solution. The crystal cores are represented by continuum finite elements with elastic moduli Xy and dielectric tensor sy. The far-field continuum solvent has viscosity p, dielectric constant s, and exerts random boundary forces Fstoch on the fluid inside the large sphere modeled using particle methods. SPC is a simple point charge model for water.

The most common boundary representation is periodic boundary conditions which assumes that the system consists of a periodic array (or a crystal ) of identical systems [1], Another common method, developed for the simulation of biomacromolecules, is the stochastic boundary approach, in which the influences of the atoms outside the boundary are replaced by a simple boundary force [78, 79, 80], Warshel uses a Langevin dipoles model in which the solvent is explicitly replaced by a grid of polarizable dipoles. The energy is calculated in a similar way to the polarization energy in a molecular mechanics force field (see above) [15]. 

An adhesion force component A y caused by the boundary force at the solid-liquid-gaseous contact line, which is determined by the surface tension of the liquid ot ... 

In Eq. 12, FB is the boundary force at r0, F(r0 — rT) is the force of interaction between a particle at rT in RR and a particle at r0 in RZ, and dTTpTg( o rr) is the probability of the pair (0, T) having a separation r<> — tT. The boundary force may be written as the gradient of a potential, the boundary potential. The boundary potential for the oxygen atom of ST2 water12 in an 11 A reaction zone is plotted in Fig. 9. In the calculation of this potential only the van der Waals part of the ST2-ST2 interaction was included in Eq. 12. A methodology that consistently incorporates electrostatic forces into the boundary potential is under development.110 In its present simplified form, the model has proven successful in the simulation of localized regions of pure... 

The second force component, the dielectric boundary pressure, results from the tendency of a high dielectric medium to displace a low dielectric constant medium if an electrical field is present. This force is always directed along the gradient of the dielectric constant, which means that it constitutes a pure pressure at the solute-solvent interface. The dielectric boundary pressure is the force component that balances the reaction field component of the qE force. For example, in the case of an isolated charge inside a low dielectric cavity, the dielectric boundary pressure provides the equal and opposite force to the qE force urging the charge toward the solvent. It is therefore evident that the dielectric boundary force is quantitatively as important as the reaction field force and that its neglect will lead to a violation of Newton s third law of motion. 

If Eq. [35] is to be used in molecular simulations, methods must be devised to compute the contribution of each of the three terms to the forces on atoms. This problem is not particularly difficult for qE forces. 2-54 However, to date, computations of atomic forces based on the PB equation either have neglected dielectric and ionic boundary forces or have dealt with them by methods whose accuracy has not yet been examined.-52.54 As recently demonstrated, ionic boundary forces will usually be negligibly weak in comparison with dielectric boundary forces.5° However, dielearic boundary forces are strong,5o and the use of inaccurate methods to compute them can lead to rather large errors in the electrostatic forces on atoms. 

It has recently become possible to compute rather accurate dielectric boundary forces from finite difference solutions of the PB equations.50 The method is based on the finite difference discretization of the second term of the force density expression (Eq. [35]), together with the use of dielectric boundary smoothing. A full description of the method and accuracy testing is beyond the scope of the present chapter. However, the results of a sample calculation may be instructive. 

Adding the two adhesion force components, one caused by the negative capillary pressure in the bridge and the other by the boundary force at the solid/liquid/gas... 

The correspondence principle states that for problems of a statically determinate nature involving bodies of viscoelastic materials subjected to boundary forces and moments, which are applied initially and then held constant, the distribution of stresses in the body can be obtained from corresponding linear elastic solutions for the same body subjected to the same sets of boundary forces and moments. This is because the equations of equilibrium and compatibility that are satisfied by the linear elastic solution subject to the same force and moment boundary conditions of the viscoelastic body will also be satisfied by the linear viscoelastic body. Then the displacement field and the strains derivable from the stresses in the linear elastic body would correspond to the velocity field and strain rates in the linear viscoelastic body derivable from the same stresses. The actual displacements and strains in the linear viscoelastic body at any given time after the application of the forces and moments can then be obtained through the use of the shift properties of the relaxation moduli of the viscoelastic body. Below we furnish a simple example. 

To overcome problems arising fi-om the finite system size used in MC or MD simulation, boundary conditions are imposed using periodic-stochastic approximations or continuum models. In particular, in stochastic boundary conditions the finite system is not duplicated but a boundary force is applied to interact with atoms of the system. This force is set as to reproduce the solvent regions that have been neglected. Anyway, in general any of the methods used to impose boundary conditions in MC or MD can be used in the QM/MM approach. 

In some variants of the immersed boundary method dealing with the boundary between a fluid and a solid body, the fictitious fluid is introduced inside the solid body, and a boundary force, F, in place of Af in Eq. 4, is distributed into the augmented flow field in the same way as described for the fluid-fluid interface (e.g., Goldstein et al. [6]). To impose the no-slip and no-penetration boundary conditions at the interface, F is modeled as... 

---