Calculating the equation of state parameters and predicting the spinodal curve of isotactic polypropylene/poly(ethylene-co-octene) blend by molecular dynamics simulations combined with Sanchez-Lacombe lattice fluid theory

Molecular dynamics simulations are adopted to calculate the equation of state characteristic parameters P*, rho*, and T* of isotactic polypropylene (iPP) and poly(ethylene-co-octene) (PEOC), which can be further used in the Sanchez-Lacombe lattice fluid theory (SLLFT) to describe the respective physical properties. The calculated T* is a function of the temperature, which was also found in the literature. To solve this problem, we propose a Boltzmann fitting of the data and obtain T* at the high-temperature limit. With these characteristic parameters, the pressure-volume-temperature (PVT) data of iPP and PEOC are predicted by the SLLFT equation of state. To justify the correctness of our results, we also obtain the PVT data for iPP and PEOC by experiments. Good agreement is found between the two sets of data. By integrating the Euler-Lagrange equation and the Cahn-Hilliard relation, we predict the density profiles and the surface tensions for iPP and PEOC, respectively. Furthermore, a recursive method is proposed to obtain the characteristic interaction energy parameter between iPP and PEOC. This method, which does not require fitting to the experimental phase equilibrium data, suggests an alternative way to predict the phase diagrams that are not easily obtained in experiments.

Effect of chain-length dependence of interaction parameter on spinodals for polydisperse polymer blends

The chain-length dependence of the Flory-Huggins (FH) interaction parameter is introduced into the FH lattice theory for polydisperse polymer-blend systems. The spinodals are calculated for the model polymer blends with different chain lengths and distributions. It is found that all the related variables r(n), r(w), r(z), and chain-length distribution, have effects on the spinodals for polydisperse polymer blends.

Estimation thermal expansion coefficient from lattice energy for inorganic crystals

By using the study of the lattice energy and the structural parameters of binary inorganic crystals, a new parameter reflecting the thermal expansion property has been found, the relation between the linear expansion coefficient and new parameter has been established. A semiempirical method for evaluation of linear expansion coefficient from the lattice energy is presented, and developed to the complex crystals. The estimated values of the linear expansion coefficients of both simple and complex crystals are in good agreement with the experimental values.

Pressure effects on the thermodynamics of trans-decahydronaphthalene/polystyrene polymer solutions: Application of the Sanchez-Lacombe lattice fluid theory

The cloud-point temperatures (T-cl's) of trans-decahydronaphthalene(TD)/polystyrene (PS, (M) over bar (w) = 270 000) solutions were determined by light scattering measurements over a range of temperatures (1-16degreesC), pressures (100-900 bar), and compositions (4.2-21.6 vol.-% polymer). The system phase separates upon cooling and T-cl was found to increase with rising pressure for constant composition. In the absence of special effects, this finding indicates positive excess volume for the mixing. Special attention was paid to the demixing temperatures as a function of pressure for different polymer solutions and the plots in the T-phi plane (where phi signifies volume fractions). The cloud-point curves of polymer solutions under different pressures were observed for different compositions, which demonstrated that pressure has a greater effect on the TD/PS solutions when far from the critical point as opposed to near the critical point. The Sanchez-Lacombe lattice fluid theory (SLLFT) was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of mixing, and the volume changes of mixing. The calculated results show that modified PS scaling parameters can describe the thermodynamics of the TD/PS system well. Moreover the SLLFT describes the experimental results well.

Theoretical estimation of thermodynamic properties of the system PS/PPO on the basis of modified combining rule of Sanchez-Lacombe lattice fluid model

The three scaling parameters described in Sanchez-Lacombe lattice fluid theory (SLLFT), T*, P* and rho* of pure polystyrene (PS), pure poly(2,6-dimethyl-1,4-phenylene oxide) (PPO) and their mixtures are obtained by fitting corresponding experimental pressure volume-temperature data with equation-of-state of SLLFT. A modified combining rule in SLLFT used to match the volume per mer, v* of the PS/PPO mixtures was advanced and the enthalpy of mixing and Flory-Huggins (FH) interaction parameter were calculated using the new rule. It is found that the difference between the new rule and the old one presented by Sanchez and Lacombe is quite small in the calculation of the enthalpy of mixing and FH interaction parameter and the effect of volume-combining rule on the calculation of thermodynamic properties is much smaller than that of energy-combining rule. But the relative value of interaction parameter changes much due to the new volume-based combining rule. This effect can affect the position of phase diagram very much, which is reported elsewhere [Macromolecules 34 (2001) 6291]

Thermodynamics of PMMA/SAN blends: Application of the Sanchez-Lacombe lattice fluid theory

Polymer blends of poly(methyl methacrylate) (PMMA) and poly(styrene-co-acrylonitrile) (SAN) with an acrylonitrile content of about 30 wt % were prepared by means of solution-casting and characterized by virtue of pressure-volume-temperature (PVT) dilatometry. The Sanchez-Lacombe (SL) lattice fluid theory was used to calculate the spinodals, the binodals, the Flory-Huggins (FH) interaction parameter, the enthalpy of the mixing, the volume change of the mixing, and the combinatorial and vacancy entropies of the mixing for the PMMA/SAN system. A new volume-combining rule was used to evaluate the close-packed volume per mer, upsilon*, of the PMMA/SAN blends. The calculated results showed that the new and the original volume-combining rules had a slight influence on the FH interaction parameter, the enthalpy of the mixing, and the combinatorial entropy of the mixing. Moreover, the spinodals and the binodals calculated with the SL theory by means of the new volume-combining rule could coincide with the measured data for the PMMA/SAN system with a lower critical solution temperature, whereas those obtained by means of the original one could not.

Phase-separation behavior of the system PES/phenoxy: An application of the Sanchez-Lacombe lattice fluid theory

Cloud-point curves reported for the system polyethersulfone (PES)/phenoxy were calculated by means of the Sanchez-Lacombe (SL) lattice fluid theory. The one adjustable parameter epsilon(12)*/k (quantifying the interaction energy between mers of the different components) can be evaluated by comparison of the theoretical and experimental phase diagrams. The Flory-Huggins (FH) interaction parameters are computed based on the evaluated epsilon(12)*/k and are approximately a linear function of volume fraction and of inverse temperature. The calculated enthalpies of mixing of PES/phenoxy blends for different compositions are consistent with the experimental values obtained previously by Singh and Walsh [1].

Application of the Sanchez-Lacombe lattice fluid theory to the system pvme/ps and model calculations

Cloud point curves reported in the literature for five representatives of the system poly(vinyl methyl ether)/polystyrene were evaluated theoretically by means of the Sanchez-Lacombe lattice fluid theory. The measured phase separation behavior can be described within experimental error using only one adjustable parameter (quantifying the interaction between the unlike mers). The Flory-Huggins interaction parameters calculated from this theoretical description depend in good approximation linearly on composition (volume fractions) and on the inverse temperature. An evaluation of these data yields a maximum heat effect which is almost one order of magnitude less (ca. -0.25 J/cm(3)) than obtained via Hess's cycle (dissolution of the components and of the blend) from calorimetric measurements. Model calculations on the basis of the present theory demonstrate that the critical points shift to a different extent upon a certain relative change in the molar mass of the blend components. The sensitivity of the calculated phase diagrams against changes in the scaling parameter decreases in the following order: interaction energies between unlike mers, differences in the scaling temperatures, pressures and densities.

STATISTICAL THERMODYNAMICS IN THE FRAMEWORK OF THE LATTICE FLUID MODEL .1. POLYDISPERSE POLYMERS OF SPECIAL DISTRIBUTION

The Gibbs free energies and equations of state of polymers with special molar mass distributions, e.g., Flory distribution, uniform distribution and Schulz distribution, are derived based on a lattice fluid model. The influence of the polydispersity (or t

Simple lattice Boltzmann model for traffic flows

A lattice Boltzmann model with 5-bit lattice for traffic flows is proposed. Using the Chapman-Enskog expansion and multi-scale technique, we obtain the higher-order moments of equilibrium distribution function. A simple traffic light problem is simulated by using the present lattice Boltzmann model, and the result agrees well with analytical solution.

Rossby Vortex Simulation on a Paraboloidal Coordinate System Using the Lattice Boltzmann Method

In this paper, we apply our compressible lattice Boltzmann model to a rotating parabolic coordinate system to simulate Rossby vortices emerging in a layer of shallow water flowing zonally in a rotating paraboloidal vessel. By introducing a scaling factor, nonuniform curvilinear mesh can be mapped to a flat uniform mesh and then normal lattice Boltzmann method works. Since the mass per unit area on the two-dimensional (2D) surface varies with the thickness of the water layer, the 2D flow seems to be "compressible" and our compressible model is applied. Simulation solutions meet with the experimental observations qualitatively. Based on this research, quantitative solutions and many natural phenomena simulations in planetary atmospheres, oceans, and magnetized plasma, such as the famous Jovian Giant Red Spot, the Galactic Spiral-vortex, the Gulf Stream, and the Kuroshio Current, etc,, can be expected.

Lattice Boltzmann方法模拟多孔介质Klinkenberg效应

格子Boltzmann数值模拟方法是研究复杂的多孔介质结构特别是Klinkenberg效应的有效方法之一，对处理复杂边值问题尤其有效。用格子Boltzmann方法研究了气流穿越多孔介质问题，并将数值计算结果与实验结果进行了比较，结果表明格子Boltzmann方法是数值模拟气流穿越多孔介质问题的有效方法之一。

A lattice Boltzmann method for KDV equation

We prepose a 5-bit lattice Boltzmann model for KdV equation. Using Chapman-Enskog expansion and multiscale technique, we obtained high order moments of equilibrium distribution function, and the 3rd dispersion coefficient and 4th order viscosity. The parameters of this scheme can be determined by analysing the energy dissipation.

Lattice Boltzmann method for compressible flows with high Mach numbers

In this paper we present a lattice Boltzmann model to simulate compressible flows by introducing an attractive force. This scheme has two main advantages: one is to soften sound speed effectively, which greatly raises the Mach number (up to 5); another is its relative simple procedure. Simulations of the March cone and the comparison between theoretical expectations and simulations demonstrate that the scheme is effective in the simulation of compressible flows with high Mach numbers, which would create many new applications.

A Continuation Method of Parameter Inversion for Non-Equilibrium Convection-Dispersion Equation

Based on the homotopy mapping, a globally convergent method of parameter inversion for non-equilibrium convection-dispersion equations (CDEs) is developed. Moreover, in order to further improve the computational efficiency of the algorithm, a properly smooth function, which is derived from the sigmoid function, is employed to update the homotopy parameter during iteration. Numerical results show the feature of global convergence and high performance of this method. In addition, even the measurement quantities are heavily contaminated by noises, and a good solution can be found.