Boundary conditions macroscopic

Surface waves at an interface between two innniscible fluids involve effects due to gravity (g) and surface tension (a) forces. (In this section, o denotes surface tension and a denotes the stress tensor. The two should not be coiifiised with one another.) In a hydrodynamic approach, the interface is treated as a sharp boundary and the two bulk phases as incompressible. The Navier-Stokes equations for the two bulk phases (balance of macroscopic forces is the mgredient) along with the boundary condition at the interface (surface tension o enters here) are solved for possible hamionic oscillations of the interface of the fomi, exp [-(iu + s)t + i V-.r], where m is the frequency, is the damping coefficient, s tlie 2-d wavevector of the periodic oscillation and. ra 2-d vector parallel to the surface. For a liquid-vapour interface which we consider, away from the critical point, the vapour density is negligible compared to the liquid density and one obtains the hydrodynamic dispersion relation for surface waves + s>tf. The temi gq in the dispersion relation arises from... 

Isolated gas ph ase molecules are th e sim plest to treat com pii tation -ally. Much, if not most, ch emistry lakes place in the liq iiid or solid state, however. To treat these condensed phases, you must simulate continnons, constant density, macroscopic conditions. The usual approach is to invoke periodic boundary conditions. These simulate a large system (order of 10" inoleeti les) as a contiruiotis replication in all direction s of a sm nII box, On ly th e m olceti Ics in the single small box are simulated and the other boxes arc just copies of the single box. 

In continuum boundary conditions the protein or other macromolecule is treated as a macroscopic body surrounded by a featureless continuum representing the solvent. The internal forces of the protein are described by using the standard force field including the Coulombic interactions in Eq. (6), whereas the forces due to the presence of the continuum solvent are described by solvation tenns derived from macroscopic electrostatics and fluid dynamics. 

There are cases where non-regular lattices may be of advantage [36,37]. The computational effort, however, is substantially larger, which makes the models less flexible concerning changes of boundary conditions or topological constraints. Another direction, which may be promising in the future, is the use of hybrid models, where for example local attachment kinetics are treated on a microscopic atomistic scale, while the transport properties are treated by macroscopic partial differential equations [5,6]. 

The constant Q is evaluated under the boundary conditions at the pipe wall r = a, u = Us. Then, the slip velocity Ug is determined (Goldstein 1965) from a macroscopic point of view ... 

Since MPC dynamics yields the hydrodynamic equations on long distance and time scales, it provides a mesoscopic simulation algorithm for investigation of fluid flow that complements other mesoscopic methods. Since it is a particle-based scheme it incorporates fluctuations, which are essential in many applications. For macroscopic fluid flow averaging is required to obtain the deterministic flow fields. In spite of the additional averaging that is required the method has the advantage that it is numerically stable, does not suffer from lattice artifacts in the structure of the Navier-Stokes equations, and boundary conditions are easily implemented. 

If the Brownian particles were macroscopic in size, the solvent could be treated as a viscous continuum, and the particles would couple to the continuum solvent through appropriate boundary conditions. Then the two-particle friction may be calculated by solving the Navier-Stokes equations in the presence of the two fixed particles. The simplest approximation for hydrodynamic interactions is through the Oseen tensor [54],... 

The Gibbs Ensemble MC simulation methodology [17-19] enables direct simulations of phase equilibria in fluids. A schematic diagram of the technique is shown in Fig. 10.1. Let us consider a macroscopic system with two phases coexisting at equilibrium. Gibbs ensemble simulations are performed in two separate microscopic regions, each within periodic boundary conditions (denoted by the dashed lines in Fig. 10.1). The thermodynamic requirements for phase coexistence are that each... 

Mehdizadeh et al. exploited the separability of current distribution on different scales to model the macroscopic current distribution on patterns made up of lines or points distributed over a large workpeice [136], They solved the secondary distribution of the superficial current density sup using a boundary condition which captures the density of small features but not their geometry. The boundary condition is based on a smoothly varying parameter representing the Faradaically active fraction of surface area. 

Molecular calculations provide approaches to supramolecular structure and to the dynamics of self-assembly by extending atomic-molecular physics. Alternatively, the tools of finite element analysis can be used to approach the simulation of self-assembled film properties. The voxel4 size in finite element analysis needs be small compared to significant variation in structure-property relationships for self-assembled structures, this implies use of voxels of nanometer dimensions. However, the continuum constitutive relationships utilized for macroscopic-system calculations will be difficult to extend at this scale because nanostructure properties are expected to differ from microstructural properties. In addition, in structures with a high density of boundaries (such as thin multilayer films), poorly understood boundary conditions may contribute to inaccuracies. 

Many other interesting examples of spontaneous reflection symmetry breaking in macroscopic domains, driven by boundary conditions, have been described in LC systems. For example, it is well known that in polymer disperse LCs, where the LC sample is confined in small spherical droplets, chiral director structures are often observed, driven by minimization of surface and bulk elastic free energies.24 We have reported chiral domain structures, and indeed chiral electro-optic behavior, in cylindrical nematic domains surrounded by isotropic liquid (the molecules were achiral).25... 

The molecular modelling of systems consisting of many molecules is the field of statistical mechanics, sometimes called statistical thermodynamics [28,29], Basically, the idea is to go from a molecular model to partition functions, and then, from these, to predict thermodynamic observables and dynamic and structural quantities. As in classical thermodynamics, in statistical mechanics it is essential to define which state variables are fixed and which quantities are allowed to fluctuate, i.e. it is essential to specify the macroscopic boundary conditions. In the present context, there are a few types of molecular systems of interest, which are linked to so-called ensembles. 

There are several commercial packages that realise the above strategy for molecularly realistic systems. It is useful to discuss some of the limitations. Ideally, one would like to do simulations on macroscopic systems. However, it is impossible to use a computer to deal with numbers of degrees of freedom on the order of /Vav. In lipid systems, where the computations of all the interactions in the system are expensive, a typical system can contain of the order of tens of thousands of particles. Recently, massive systems with up to a million particles have been considered [33], Even for these large simulations, this still means that the system size is limited to the order of 10 nm. Because of this small size, one refers to this volume as a box, although the system boundaries are typically not box-like. Usually the box has periodic boundary conditions. This implies that molecules that move out of the box on one side will enter the box on the opposite side. In such a way, finite size effects are minimised. In sophisticated simulations, i.e. (N, p, y, Tj-ensembles, there are rules defined which allow the box size and shape to vary in such a way that the intensive parameters (p, y) can assume a preset value. 

Tieleman, D. P. and Berendsen, H. J. C. (1996). Molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine bilayer with different macroscopic boundary conditions and parameters, J. Chem. Phys., 105, 4871 —4880. 

Recent times have seen much discussion of the choice of hydrodynamic boundary conditions that can be employed in a description of the solid-liquid interface. For some time, the no-slip approximation was deemed acceptable and has constituted something of a dogma in many fields concerned with fluid mechanics. This assumption is based on observations made at a macroscopic level, where the mean free path of the hquid being considered is much smaller... 

In most problems involving boundary conditions, the boundary is assigned a specific empirical or deterministic behavior, such as the no-slip case or an empirically determined slip value. The condition is defined based on an averaged value that assumes a mean flow profile. This is convenient and simple for a macroscopic system, where random fluctuations in the interfacial properties are small enough so as to produce little noise in the system. However, random fluctuations in the interfacial conditions of microscopic systems may not be so simple to average out, due to the size of the fluctuations with respect to the size of the signal itself. To address this problem, we consider the use of stochastic boundary conditions that account for random fluctuations and focus on the statistical variability of the system. Also, this may allow for better predictions of interfacial properties and boundary conditions. 

Instabilities and dissipative structures play a basic role in our understanding of chemical kinetics in far from equilibrium situations. Through dissipative structures the characteristics of the system become explicitly dependent on macroscopic parameters (dimensions, boundary conditions) that relate it to its environment. 

Surface Tension Effect. As another aspect, when K is very small, we may well expect that the surface tension a at the gel-solvent interface comes into play in the macroscopic instability [31]. Since the gel radius is R + u(R, t), the modified boundary condition Eq. (Z42) is... 

Sensitivity to Boundary Conditions. The macroscopic instability crucially depends on the boundary condition and becomes nonexistent when all the gel surfaces are clamped, while the bulk instability occurs in the interior region and is insensitive to the boundary condition. For example, we set u(R) — 0 in Eq. (4.20) for spherical gels to obtain tan s/s = 1, from which Sq 0.11. The slowest relaxation rate is thus fi0 S 0.11(K + ft) lR2 and goes to zero at K = — ffi in the damped spherical case. 

In closing this section, we note that the stochastic master equation, Eq. (16), can be used to study the effect of boundary conditions on transport equations. If a(x, y, t) is sufficiently peaked as a function of x — y, that is if transitions occur from y to states in the near neighborhood of y, only, then the master equation can be approximated by a Fokker-Planck equation. The effects of the boundary on the master equation all appear in the properties of a(x, y, t). However, in the transition to the Fokker-Planck differential equation, these boundary effects appear as boundary conditions on the differential equation.7 These effects are prototypes for the study of how molecular boundary conditions imposed on the Liouville equation are reflected in the macroscopic boundary conditions imposed on the hydrodynamic equations. 

To find a principle which so narrows down the statistics (e.g., to a subclass of probability densities S ) that 9 relations in addition to the 5 transport equations will hold between our macroscopic quantities, so that the initial value problem for (p, ua, E), with the basic boundary conditions, will become determinate. In actuality, this will amount to finding a principle from which the various phenomenological relations of two paragraphs back are (approximate) deductions. 

Within broad limits it is true that if basic macroscopic boundary conditions are imposed, as in Section II, the boundary system and motions can be such as to bring them into effect. But even with this restriction, the essentially stochastic element in w(X, t) remains. 

---