Closure approximation analytic

Analytic expressions for the Gaussian spread parameters cr may be derived from the Navier-Stokes equation, through the implementation of second-order closure approximations (Sykes et al., 2004). 

The numerical solution of Eqs.(108)-(111) suffers from the fact that the evaluation of requires a suuunation over the complete positive energy spectrum even within the no-pair approximation. In the nonrelativistic case the semi-analytical KLI approximation [72] to the 0PM proved to be highly accurate for atoms [64], molecules [73], and solids [138]. The KLI approximation can be most easily extended into the relativistic domain by use of a closure approximation for Gt [54],... 

Specific results require two additional pieces of information. First, the single diain polymer structure as contained in the species-dependent structure factors must Ik known. In principle, these correlation fimctions should be self-consis-tently computed along with the intermolecular pair correlations. This most rigorous approach is now feasible and is briefly discussed in Sect. 10. However, all published PRISM work on blends has employed the Flory ideality ansatz discussed in Sect. 2, i.e., the required functions are presumed given and effectively independent of blend composition and proximity to the spinodal. To date, published analytical and numerical work for alloys has employed coarsegrained models sudi as the Gaussian, freely-jointed, and Kmi-fiexible chain. Second, a closure approximation is required which relates the intermolecular pair and direct site-site correlation functions outside the hard core diameter. Both traditional atomic and novel moleculaF dosures have been studied. 

The apparent SANS chi-parameter is also easily determined analytically for stiffness asymmetric Berthelot thread model with the R-MMSA or R-MPY/ HTA closure approximations. For algebraic simplicity we consider the neutron data analysis approach which leads to Eq. (6.IS). In the effectively incompressible regime, defined here as Cmm- > IPHmm- in Eq. (8.11), one easily obtains the result [67]... 

For the analytically tractable thread polymer model, and the R-MMSA or R-MPY/HTA closure approximations, k = 0 values of the direct correlation functions are precisely the same in the long chain limit as found for polymer blends in Sect. 8. In particular, for the symmetric block coprdymer, the R-MMSA closure yields [67,86] the mean field result Xinc = 3Co- Thus, within the symmetric thread idealization and the incompres fe approximation ctf Eq. (9.5), PRISM/R-MMSA theory reduces to Leibler theory for all compositions and block architectures [67,86]. 

Liquid-vapor equilibrium of chain molecule fluids. Both analytic and numerical work has been recently done by Schweizer and co-workers. The compressibility route predictions of PRISM for this problem are extremely sensitive to closure approximation since the relevant fluid densities are very low and large-scale density fluctuations are present. The atomiclike MSA closure leads to qualitatively incorrect results as does the R-MMSA closure. However, the R-MPY/HTA approximation appears to be in excellent accord with the computer simulation studies of n-alkanes and model chain polymers, including a critical density that decreases weakly with N and a critical temperature that increases approximately logarithmically with N. 

Solution of the PRISM equation for a particular system requires determining the intramolecular correlation functions. Day (k), and specifying a closure approximation. For rigid molecules, Qay (k) can be determined analytically. For example, the... 

The EMSA requires the degree of dimerization A as an input parameter. This is quite disappointing. However, it ehminates the deficiency of the Percus-Yevick approximation, Eq. (38). The EMSA represents a simpHfied version, to obtain an analytic solution, of a more sophisticated site-site extended mean spherical approximation (SSEMSA) [67-69]. The results of the aforementioned closures can be used as an input for subsequent calculations of the structure of nonuniform associating fluids. 

The comparisons presented in Figures 7 and 8, show that the RY closure relation works better than HNC, PYand RMSA, for the two kinds of systems with repulsive interactions here considered. Before the introduction of the RY approximation, the picture was that HNC and PY were better approximations to describe the static structure of systems with repulsive long-range and hard-sphere interactions, respectively. The RMSA approximation, however, has been used extensively in the comparison with experimental data for the static structure of aqueous suspensions of polystyrene spheres, mainly because it has an analytical solution even for mixtures [36]. 

The RISM integral equations in the KH approximation lead to closed analytical expressions for the free energy and its derivatives [29-31]. Likewise, the KHM approximation (7) possesses an exact differential of the free energy. Note that the solvation chemical potential for the MSA or PY closures is not available in a closed form and depends on a path of the thermodynamic integration. With the analytical expressions for the chemical potential and the pressure, the phase coexistence envelope of molecular fluid can be localized directly by solving the mechanical and chemical equilibrium conditions. 

This implies that higher order moments are introduced, thus the system of PDEs cannot be closed analytically. It is possible to show that similar effects will occur for the other source terms as well. This problem limits the application of the exact method of moments to the particular case where we have constant kernels only. In other cases one has to introduce approximate closures in order to eliminate the higher order moments ensuring that the transport equations for the moments of the particle size distribution can be expressed in terms of the lower order moments only (i.e., a modeling process very similar to turbulence modeling). 

The LHNC and QHNC approximations have not been solved analytically, but numerical solutions can be obtained by iteration. This is also true of the MSA except for the previously discussed dipolar hard-sphere system solved by Wertheim. The details of the numerical solution are described in Refs. 30, 38, 58, and 59. Essentially, (3.11) and the appropriate closure relations are written in terms of c" " and -q "" and iterated until a solution is obtained. This means that all equations defining a particular approximation are simultaneously satisfied. The present problem is very similar to that... 

The closed analytical expressions (4.14) or (4.15) for the excess chemical potential of solvation are no more valid in the SC-3D-RISM approach, since the orientational averaging (4.56) breaks the symmetry of the 3D-RISM integral equation with respect to the solvent indices. Nevertheless, the solvation chemical potential obtained from the SC-3D-RISM/HNC equation (4.59) does take the HNC form (4.14) within the additive approximation (4.64), (see Appendix). The use of the 3D-KH closure leads with account of (4.64) to the solvation chemical potential (4.15). 

Curro and Schweizer have carried out numerical [60,61] and analytical [23] studies of the symmetric blend using the Mean Spherical Approximation (MSA) closure successfully employed for atomic, colloidal, and small molecule fluids [5,6]. This closure corresponds to the approximation ... 

Extensive analytic results for the symmetric thread blend have also been derived [68,70b]. In the thread-polymer limit the hard core condition becomes irrelevant for the molecular closure relations. In particular, for the R-MMSA and R-MPY/HTA approximations the MM (k) functions are fully specified by the closure relations, and their k = 0 values are given in general by... 

The full solution of the binding MSA for dimer association was discussed elsewhere (BIMSA)[48, 39]. Imposing an exponential closure reminiscent of Bjerrum s approximation [39] for the contact pair distribution function results in simple analytic expressions for the excess thermodynamic quantities. 

This method is employed for problems where analytical or differential equations approaches are not feasible in view of high dimensionality, or only lead to approximate solutions or averages for one or more dimensions. Examples of such limitations have been discussed above. The Galerkin-FEM method in the pseudodistribution mode was confronted with closure problems in the case of PVAc with more than one TDB per chain. The pgf-method is applicable only to systems with identical statistics of the constituting elements. In principle, Monte Carlo simulation does not suffer from such limitations. Here, we will introduce classical MC, but subsequently mainly discuss some successful applications of advanced MC methods. 

To understand the processing and dispensing behavior of adhesives we often need to construct flow models of the process. While all models are approximations, they still have to satisfy continuity of mass and the balance of momentum and thermal energy. The constitutive equation, plus the appropriate initial and/or boundary conditions of the problem at hand, provide closure to these balance laws. While any realistic solution to a particular processing problem will generally involve numerical computations, several generic problems of interest may be amenable to analytical development. For instance, when a pressure-sensitive adhesive is pressed down unto a surface, we have an example of squeeze flow. Similarly, if a paste is spread onto a surface via a knife, we have an example of wedge flow. These two flows will be discussed for the PLF. 

Based upon the effective-mode construction, a systematic approximation procedure for the environment can be formulated in terms of a series of coarse-grained spectral densities [32,33]. These spectral densities are generated from successive orders of a truncated chain model with Markovian closure. Analytical expressions can be given in terms of Mori type continued fractions. Assuming that an - a priori arbitrarily complicated - reference spectral density can be obtained independently, e.g., from experiments or classical simulations, one can thus (1) extract those features of the spectral density that determine the interaction with the subsystem... 

Integral Equation and Eield-Theoretic Approaches In addition to theories based on the direct analytical extension of the PB or DH equation, PB results are often compared with statistical-mechanical approaches based on integral equation or density functional methods. We mention only a few of the most recent theoretical developments. Among the more popular are the mean spherical approximation (MSA) and the hyper-netted chain (HNC) equation. Kjellander and Marcelja have developed an anisotropic HNC approximation that treats the double layer near a flat charged surface as a series of discrete layers.Attard, Mitchell and Ninham have used a Debye-Hiickel closure for the direct correlation function to obtain an analytical extension (in terms of elliptic integrals) to the PB equation for the planar double layer. Both of these approaches, which do not include finite volume corrections, treat the fluctuation potential in a manner similar to the MPB theory of Outhwaite. 

---