Distinctive features of random local motions in critical and supercritical fluids

Equations for large-scale local fluctuations in fluids, from an ideal gas to an incompressible fluid, including the critical and supercritical state are derived for the first time based on the first principles. The modern phenomenological representation of the critical state of fluids is confirmed and essentially refined; in particular, it is demonstrated that that local density fluctuations in a compressible fluid are accompanied by nonthermodynamic fluctuations in the collective velocity and temperature of the fluid. Distinctive features of the development of these fluctuations near the critical point determine the specific behavior of fluids in the critical and supercritical states.     

热噪声的相干态和压缩态中的高阶涨落

利用热场动力学的方法研究了具有热噪声的相干态和压缩态中场正交分量的高阶涨落和高阶压缩.利用测量相位算符讨论了这些热化态中相位的高阶涨落.从而得到在有限温度下这些态的低阶及高阶涨落与温度的关系和压缩特性与温度的关系.由于实际的场态总是在一有限的温度下，所以讨论温度对涨落和压缩的影响是很有意义的. 关键词：     

The phase behaviour of a fluid in a slitlike pore with permeable walls

The confinement of a lattice fluid in a set of slitlike pores separated by semipermeable walls with a finite width has been studied. The walls are modelled by a square-well repulsive potential with a finite height. The thermodynamic properties and the phase behaviour of the system are evaluated by means of Monte Carlo simulations. For some states theoretical calculations have been made using a mean-field-type theory. These investigations confirm previous findings for confined Lennard-Jones fluids, obtained from a density functional approach. For intermediate and low potential barriers that separate the pores, the isotherms exhibit two hysteresis loops and the liquid-vapour coexistence curve divides into two branches describing condensation inside the pore and inside the permeable wall. These two branches are separated by a triple point. At temperatures lower than the triple point temperature, the condensation takes place instantaneously in both the pore and inside the permeable wall. It was found that when the temperature is scaled by the bulk critical temperature, the phase diagram emerging from this simple mean-field treatment is close to the phase diagram obtained from simulation.     

Properties of the quantum universe in quasistationary states and cosmological puzzles

Old and new puzzles of cosmology are reexamined from the point of view of the quantum theory of the universe developed here. It is shown that in the proposed approach the difficulties of the standard cosmology do not arise. The theory predicts the observed dimensions of the non-homogeneities of matter density and the amplitude of fluctuations of the cosmic background radiation temperature in the Universe and points to a new quantum mechanism of their origin. The large-scale structure in the Universe is explained by the growth of non-homogeneities which arise from primordial quantum fluctuations due to the finite width of the quasistationary states. The theory allows one to obtain the value of the deceleration parameter, which is in good agreement with the recent SNe Ia measurements. It explains the large value of the entropy of the Universe and describes other parameters. Received: 30 July 2001 / Published online: 8 February 2002     

The effects of interaction range,porosity and molecular association on the phase equilibrium of a fluid confined in a disordered porous media

Grand-canonical transition-matrix Monte Carlo (GC-TMMC) is employed to analyse the effects of range of interaction, packing fraction and molecular association on phase coexistence properties of square-well (SW) based fluids in disordered pores. The nature of the phase equilibria were studied inside a repulsive disordered porous media with packing fractions, η _m = 0.05 and 0.10. Three values of the SW attractive well range parameter were studied: λ = 1.5, 1.75, and 2.0. Coexistence number probability distribution reflects the signature of the disordered structure of the porous matrix. Yet, no multiple fluid–fluid transition was observed. The effect of strength of molecular association on coexistence densities, density profile, saturation pressure, and monomer fraction for the SW based dimerizing fluids inside a repulsive disordered media is reported. Association is found to increase as the packing fraction of the matrix increase. Critical properties of these confined fluids are calculated via a rectilinear diameter approach. Fractional shift in the critical temperature linearly decreases with the increase in the attractive well width for non-associating fluids. The rate of decrease in the critical temperature shift increases with the increase in packing fraction. Associating sites are found to suppress the shift in the critical temperature.     

Correlation functions for simple fluids in a finite system under nonequilibrium constraints

The Landau-Lifschitz fluctuating hydrodynamics formalism is applied to study the statistical properties of simple fluids in a finite system under nonequilibrium constraints. The boundary conditions are explicitly taken into account so that the results can be compared with particle simulations. Two scenarios are investigated: a fluid subjected to a constant shear and a fluid subjected to a constant temperature gradient. By considering a fluid with vanishing thermal expansivity, exact results are obtained for the static and dynamic correlation functions.     

具有有限温度的热克尔态

利用热场动力学方法，提出了具有有限温度的热克尔态，研究了这些态的量子统计性质及其与温度的关系，讨论了这些态中场振幅的涨落与压缩性质．对于热化的场，还提出了热压缩的定义，并对热克尔态，讨论了温度对热压缩的影响     

Long-range order and dynamic structure factor of a nematic under a thermal gradient

We use a fluctuating hydrodynamic approach to calculate the orientation fluctuations correlation functions of a thermotropic nematic liquid crystal in a nonequilibrium state induced by a stationary heat flux. Since in this nonequilibrium stationary state the hydrodynamic fluctuations evolve on three widely separated times scales, we use a time-scale perturbation procedure in order to partially diagonalize the hydrodynamic matrix. The wave number and frequency dependence of these orientation correlation functions is evaluated and their explicit functional form on position is also calculated analytically in and out of equilibrium. We show that for both states these correlations are long-ranged. This result shows that indeed, even in equilibrium there is long-range orientational order in the nematic, consistently with the well known properties of these systems.We also calculate the dynamic structure of the fluid in both states for a geometry consistent with light scattering experiments. We find that as with isotropic simple fluids, the external temperature gradient introduces an asymmetry in the spectrum shifting its maximum by an amount proportional to the magnitude of the gradient. This effect may be of the order of 7 per cent. Also, the width at half height may decrease by a factor of about 10 per cent. Since to our knowledge there are no experimental results available in the literature to compare with, the predictions of our model calculation remains to be assessed.     

有限温度下一维Gaudin-Yang模型的热力学性质

本文通过数值求解有限温度下一维均匀费米Gaudin-Yang模型的热力学Bethe-ansatz方程, 研究了此模型的基本性质,得到了在给定的温度或给定的相互作用下, 化学势、相互作用、粒子密度和熵的相互变化图像. 对结果分析发现, 在给定温度和相互作用下, 熵随着化学势的变化有一个量子临界区域.     

任意不可压多介质流的粒子有限元方法模拟

基于粒子有限元方法(particle finite element method,PFEM),利用细分混合单元的界面识别思想,模拟种类任意多的不可压多介质流问题.对分步算法采用基于有限增量微积分理论的稳定措施,以适应流体特性差异;将混合单元细分为代表单一流体的小单元,进而得到流体间的边界;通过加密边界、控制粒子速度、自动检查穿透来防止粒子穿透外部边界.瑞利-泰勒不稳定性和水柱在空气中倒塌的模拟与已有结果的对比验证了PFEM及界面识别方法的可靠性和准确性.七种流体混合的模拟结果表明PFEM可有效处理任意多种类不相溶流体的混合流动问题.     

Fluctuations in fluids out of thermal equilibrium

After a brief review of dynamic correlations in equilibrium fluids, we consider the long-range correlations between the fluctuations in a fluid subjected to a large stationary temperature gradient. These long-range correlations enhance and modify the Rayleigh spectrum of the fluid. We elucidate that the modifications of the Rayleigh line are determined by the coupling of the entropy fluctuations to the transverse velocity fluctuations. Recent attempts to test the theoretical predictions with the aid of light-scattering experiments are discussed.     

Stability of fluid flow past a membrane

The stability of the flow of a fluid past a solid membrane of infinitesimal thickness is investigated using a linear stability analysis. The system consists of two fluids of thicknesses R and H R and bounded by rigid walls moving with velocities and , and separated by a membrane of infinitesimal thickness which is flat in the unperturbed state. The fluids are described by the Navier-Stokes equations, while the constitutive equation for the membrane incorporates the surface tension, and the effect of curvature elasticity is also examined for a membrane with no surface tension. The stability of the system depends on the dimensionless strain rates and in the two fluids, which are defined as and for a membrane with surface tension , and and for a membrane with zero surface tension and curvature elasticity K. In the absence of fluid inertia, the perturbations are always stable. In the limit , the decay rate of the perturbations is O(k ³ ) smaller than the frequency of the fluctuations. The effect of fluid inertia in this limit is incorporated using a small wave number asymptotic analysis, and it is found that there is a correction of smaller than the leading order frequency due to inertial effects. This correction causes long wave fluctuations to be unstable for certain values of the ratio of strain rates and ratio of thicknesses H. The stability of the system at finite Reynolds number was calculated using numerical techniques for the case where the strain rate in one of the fluids is zero. The stability depends on the Reynolds number for the fluid with the non-zero strain rate, and the parameter , where is the surface tension of the membrane. It is found that the Reynolds number for the transition from stable to unstable modes, , first increases with , undergoes a turning point and a further increase in the results in a decrease in . This indicates that there are unstable perturbations only in a finite domain in the plane, and perturbations are always stable outside this domain. Received: 29 May 1997 / Revised: 9 October 1997 / Accepted: 26 November 1997     

Multi-layer flows of immiscible fractional Maxwell fluids with generalized thermal flux

Unsteady laminar flows and heat transfer of n-immiscible fractional Maxwell fluids in a channel are investigated under influence of time-dependent pressure gradient. The isothermal channel walls have translational motions in their planes with time-dependent velocities. Governing equations of the mathematical model are based on the generalized constitutive equations for shear stress and thermal flux described by the time-fractional Caputo derivative. Analytical and semi-analytical solutions for velocity, shear stress, and temperature fields are obtained by using finite sine-Fourier and Laplace transforms. In the case of semi-analytical solutions, the inverse Laplace transforms are obtained numerically by employing the Talbots algorithms. Using the software Mathcad, numerical calculations have carried out and results are presented in graphical illustrations in order to analyze the memory effects on the fluid temperature and motion. It is found that in fluids with thermal memory the heat transfer is slower compared with the ordinary fluid, while the fractional velocity parameters act as braking/accelerating factors of the fluids.     

The Alpha Effect: The Connection between Cyclonic Events and Current Helicity

In a magnetized plasma, resistive diffusion of large-scale magnetic fields can be suppressed or even overcome by a turbulently generated electromotive force. For a plasma in which the turbulence is homogeneous and isotropic this EMF is characterized by the ensemble average = ?B0, where ?v and ?b represent the turbulent fields and B0 defines the large-scale field. Determination of the statistical properties of the turbulence that are required to generate a finite alpha effect, as it has become known, is one of the central subjects of dynamo theory. Parker has shown that helical velocity fluctuations possessing a net amount of kinetic helicity are capable of dynamo action. These "cyclonic events" produce electromagnetic fluctuations characterized by their own statistical properties. Within the context of "mean-field electrodynamics" we show that these fluctuations possess a net amount of current helicity, and find that a necessary condition for dynamo action is that the turbulent current helicity and the current helicity in the large-scale field be of opposite sign.     

Fokker-Planck equation approach to fluctuations about nonequilibrium steady states

We examine the properties of steady states in systems which interact at the boundary with a nonequilibrium environment. The examination is based on a nonlinear Fokker-Planck equation, the structure of which is determined by the fact that it also governs the time evolution of the equilibrium fluctuations of the system. The nonlinearities in the Fokker-Planck equation may have two origins: thermodynamic nonlinearities which arise if the thermodynamic potential is not a bilinear function of the state variables, and nonlinear mode coupling which arises if the transport coefficients depend on the state. While these nonlinearities have only a small effect on the equilibrium fluctuations of a system away from critical points, they are shown to be important for the determination of fluctuations about nonequilibrium steady states. In particular the state dependence of the transport coefficients may lead to deviations from local equilibrium and to a breakdown of detail balance. An explicit formula for the time correlations of fluctuations about the nonequilibrium steady state is obtained. The formula leads to long-range correlations in fluids in the presence of a temperature gradient. The result is compared with earlier approaches to the same problem. Finally, we study the linear response to external forces and obtain a generalization of the fluctuation-dissipation formula relating the response functions with the nonequilibrium correlation functions.     

Vapour pressure of non-polar fluids and their repulsive and attractive contributions

Using a simple equation of state, based on the Weeks-Chandler-Andersen separation of the intermolecular potential, we have obtained the contributions of repulsive and attractive intermolecular forces to the thermodynamic properties of coexisting vapour and liquid phases of a Lennard-Jones (LJ) fluid.

In order to obtain the vapour pressure of real non-polar fluids, we take the LJ fluid as a reference model, and propose a new perturbative contribution, which is dependant on the temperature and on the acentric factor of the substance. Using the complete perturbed equation, we determine the corresponding repulsive and attractive contributions to the vapour pressure of non-polar fluids. The results show that the attractive vapour pressure of non-polar fluids increases with increasing acentric factor, i.e., larger deviation of the molecular shape from spherical symmetry.

This procedure could be extended to separate the repulsive and attractive contributions of the intermolecular forces to other thermodynamic properties of non-polar fluids as well as of polar fluids and fluid mixtures.     

Thermal fluctuations in d-wave layered superconductors

We study the thermal fluctuations of anisotropic order parameters (OP) in layered superconductors. In particular, for copper oxides and a d-wave OP, we present some experimental consequences of fluctuations in the direction normal to the layers. It is shown that the c-axis penetration depth λ_c can have a “disorder-like” quadratic temperature dependence at low temperature. The fluctuations are analyzed in the framework of a Lawrence-Doniach model with an isotropic Fermi surface. Anisotropies pin the orientation of the OP to the crystallographic axes of the lattice. Then we study an extended t-J model that fits Fermi surface data of bilayers YBCO and BSCCO. This leads to a d-wave OP with two possible orientations and, including the thermal fluctuations, yields the announced temperature dependence of λ_c. Furthermore a reservoir layer is introduced. It implies a finite density of states at the Fermi energy which is successfully compared to conductance and specific heat measurements.     

Back-reaction equations for isotropic cosmologies when nonconformal particles are created

A model proposed some years ago by Hartle to study the back reaction in a cosmological model due to the creation of massless non-conformally coupled particles is reexamined. The model consists of a spatially flat FRW spacetime with a classical source made of two perfect fluids one a radiative fluid and the other a baryonic fluid with the equation of state of dust, and it is assumed that the ratio of baryons to photons is small. The back-reaction equations for the cosmological scale factor are derived using a CTP (closed time path) effective action method. Making use of the connection, in the semiclassical context, between the CTP effective action and the influence functional in quantum statistical mechanics, improved back-reaction equations are derived which take into account the fluctuations of the stress-energy tensor of the quantum field. These new dynamical equations are real and causal and predict stochastic fluctuations for the cosmological scale factor.     

Monte Carlo simulations in the vicinity of the critical point: Vapour-liquid coexistence curve

A generalized Gibbs-ensemble methodology with a fluctuating particle is used to determine the coexistence vapor-liquid densities for a square-well fluid. It is shown that the presence of the fluctuating particle in sub-states of the Markov chain of states supresses considerably density fluctuations which makes it possible to carry out simulations efficiently even for temperatures very close to the critical temperature.