Dewetting in associating lattice gas model confined by hydrophobic walls

The phase behavior of a two dimensional fluid confined within hydrophobic walls is obtained by Monte Carlo simulations. The fluid is described by the associating lattice gas model which reproduces the density and diffusion anomalous behavior of water.The confined fluid exhibits a liquid-liquid critical temperature which decreases with the decrease of the distance between the confining walls. In contact with the wall a dewetting is observed. The thickness of this interfacial layer is independent of the distance between the two walls. Even for very small distances between the two walls no total depletion is observed and consequently no drying transition is present.     

受限导致硬椭球液体的取向有序

Monte Carlo模拟研究了两平行硬墙受限下各向异相的硬椭球液体形成的有序结构. 墙壁对椭球粒子的过量吸附和椭球的长短径之比有关;对于高密度的椭球液体,它在一个临界的长短径之比 a/b=2.9处达到最大值,表明在墙表面附近发生了取向有序. 墙表面附近的密度分布和取向序参量进一步证实这结论. 相同密度条件下体相中的取向序参量计算表明体相结构均为各向同性,这是由于体系的密度依然低于各向同性相向相列向转变的临界密度. 该受限椭球液体中墙表面的取向有序是由墙面的几何受限所导致的.     

Effect of side walls on the motion of a viscous fluid induced by an infinite plate that applies an oscillating shear stress to the fluid

The velocity field corresponding to the unsteady motion of a viscous fluid between two side walls perpendicular to a plate is determined by means of the Fourier transforms. The motion of the fluid is produced by the plate which after the time t = 0, applies an oscillating shear stress to the fluid. The solutions that have been obtained, presented as a sum of the steady-state and transient solutions satisfy the governing equation and all imposed initial and boundary conditions. In the absence of the side walls they are reduced to the similar solutions corresponding to the motion over an infinite plate. Finally, the influence of the side walls on the fluid motion, the required time to reach the steady-state, as well as the distance between the walls for which the velocity of the fluid in the middle of the channel is unaffected by their presence, are established by means of graphical illustrations.     

Integral equation theory of a one dimensional hard-ellipse fluid between hard walls

Density functional theory is used to study the structure of a one dimensional fluid model of hard-ellipse molecules with their axes freely rotating in a plane, confined between hard walls. A simple Hypernetted chain (HNC) approximation is used for the density functional of the fluid and the integral equation for the density is obtained from the grand potential. The only required input is the direct correlation function of the one dimensional hard-ellipse fluid. For this model, the pressure, sum rule and the density at the walls are obtained. The Percus Yevick (PY), for lower density, and HNC, for higher density, integral equations are also solved to obtain the direct correlation function of hard-ellipse model introduced here. We obtain the average density at the wall as well as the radial density profile. We compare these with Monte Carlo simulations of the same model and find reasonable agreement.     

A novel fluid–wall heat transfer model for molecular dynamics simulations

The ‘fluid–wall thermal equilibrium model’, to numerically simulate heating/cooling of fluid atoms by wall atoms, is used to compare molecular dynamics simulation results to the analytical solution of 1-D heat equation. Liquid argon atoms are placed between two platinum walls and simultaneous heating and cooling is simulated at the walls. Temperature gradient in liquid argon is evaluated and the results are found to match well with the analytical solution showing the physical soundness of the proposed model. Additional simulations are done where liquid argon atoms are heated by both the walls for two different channel heights and it is shown that in such cases, heat transfer occurs at a faster rate than predicted by heat equation with decreasing channel heights.     

Studies of diffusion coefficients in disordered porous matrices confined in a slit-pore

We present a series of molecular dynamics simulations to study the structure of porous matrices confined in a slit-pore. The matrices were prepared by two different methods. In the first method we used direct simulations of a fluid at a fixed density and the matrix was taken from the last configuration of its particles. In the second method we simulated a binary mixture where one of the components served as a template material and the final porous matrix configuration was obtained by removing template particles from the mixture. In both methods the matrices were confined by two parallel walls (slit-pore) modeled by continuous solid surfaces. The results show that the matrix structure and porosity were affected by the method of preparation of the porous matrices. Moreover, we found smaller void cavities in these matrices than in matrices prepared without walls. Finally, diffusion of a fluid inside the matrices was investigated and it was found that the diffusion coefficient did not decrease with the fluid density, and presented a maximum at certain values of the fluid density.     

Phase behaviour of a Lennard-Jones fluid in a pore with permeable walls of a finite thickness: a density functional approach

A confinement of a Lennard-Jones fluid in a system of slitlike pores separated by semipermeable walls of a finite width is studied. The walls are modelled by square-well repulsive potential wells. The structure of the confined fluid is investigated by means of a density functional method. For high potential barriers separating the pores, the phase behaviour of the system is similar to that for a single slitlike pore with impenetrable walls. For intermediate and low potential barriers the system shows different phase behaviour. Within some temperature range the isotherms exhibit two hysteresis loops, which characterize the condensation of the fluid in different parts of the system, namely in the pore and inside the semipermeable walls. The systems characterized by low and intermediate potential barriers exhibit the triple point, such that at temperatures below that triple point the condensation instantaneously takes place in both the pore and inside the permeable wall.     

Magnetohydrodynamic plane Couette flow of an elasto-viscous fluid

Summary Using the Laplace transform, an analytical solution of the Navier-Stokes equation is obtained for a two-dimensional incompressible elasto-viscous fluid past between two infinite parallel walls. It is assumed that the lower wall is moving with velocity which is a function of any given free stream velocity. As an application of the solution, two cases for the stream velocity are studied.     

Second-order Percus-Yevick theory for a confined hard-sphere fluid

A fluid of hard spheres confined between two hard walls and in equilibrium with a bulk hard-sphere fluid is studied using a second-order Percus-Yevick approximation. We refer to this approximation as second-order because the correlations that are calculated depend upon the position of two hard spheres in the confined fluid. However, because the correlation functions depend upon the positions of four particles (two hard spheres and two walls treated as giant hard spheres), this is the most demanding application of the second-order theory that has been attempted. When the two walls are far apart, this calculation reduces to our earlier second-order approximation calculations of the properties of hard spheres near a single hard wall. Our earlier calculations showed this approach to be accurate for the single-wall case. In this work we calculate the density profiles and the pressure of the hard-sphere fluid on the walls. We find, by comparison with grand canonical Monte Carlo results, that the second-order approximation is very accurate, even when the two walls have a small separation. We compare with a singlet approximation (in the sense that correlation functions that depend on the position of only one hard sphere are considered). The singlet approach is fairly satisfactory when the two walls are far apart but becomes unsatisfactory when the two walls have a small separation. We also examine a simple theory of the pressure of the confined hard spheres, based on the usual Percus-Yevick theory of hard-sphere mixtures. Given the simplicity of the latter approach the results of this simple (and explicit) theory are surprisingly good.     

The influence of viscous heating and wall thermal conditions on the thermal ignition of a poiseuille/couette reactive flow

This paper investigates the stationary thermal stability of an exothermic reactive viscous non-Newtonian fluid between two parallel walls in the plane Poiseuille and generalized Couette flow configurations for two different thermal conditions. Assuming the system is adiabatic with negligible reactant consumption, the closed-form solution obtained from the momentum equation was inserted into the energy equation due to the dissipative effect of viscosity. The resulting energy equation was analyzed for criticality using a variational technique. The influence of the viscous heating parameter (Γ), wall dynamical (Λ) and the thermal conditions of the walls on the thermal ignition parameters were examined.     

X-Ray waveguiding studies of ordering phenomena in confined fluids

We have determined the structure of a colloidal fluid confined in a gap between two walls by making use of the waveguiding properties of the gap at x-ray wavelengths. The method is based on an analysis of the coupling of waveguide modes induced by the density variations in the confined fluid. Studies on suspensions confined within gaps of a few hundred nanometers showed strongly selective mode coupling effects, indicative of an ordering of the colloidal particles in layers parallel to the confining walls.     

Hydromagnetic elasto-viscous fluid flow between two parallel walls: Temperature field. — II

Summary A formula for the temperature field is obtained for elasto-viscous flow between two parallel walls, one is fixed and the other moves periodically in time. The effects of the elastic parameter on temperature and heat transfers at the walls are obtained.     

Classical Coulomb fluids in a confined geometry

It has already been argued that a classical (three-dimensional) Coulomb fluid confined between two parallel walls exhibits ideal gas features when the distance between the walls becomes small; this is confirmed in the present paper. Two-dimensional models of Coulomb fluids (with a logarithmic interaction), confined in a strip, are also studied. These models do not become ideal gases in the narrow strip limit. The correlation functions are also studied. There is a special temperature at which exact results are obtained. At that temperature, the two-dimensional, two-component plasma (two-dimensional Coulomb gas), which is a conductor when unconfined, becomes a dielectric as soon as it is confined in a strip of noninfinite width. This can be understood as a displacement of the Kosterlitz-Thouless transition by the confinement.     

Monte Carlo simulation of a confined hard ellipse fluid

The structural and thermodynamic properties of a confined hard ellipse fluid are studied using Monte Carlo simulation. The angular, average number densities and order parameters of hard ellipses confined between hard parallel walls are obtained for various bulk densities, aspect ratios and wall separations. The results show that the effect of the existence of the wall on the molecular fluid structure, either on their directions or their locations, with respect to the bulk, especially close to the walls, is significant. For this system the pressure is also obtained and it is shown that the average density at the wall is proportional to the pressure, βP=〈ρ_w〉. Our simulation results show that the order parameters depend on the number of the particles in the box unless it exceeds thousand.     

Dynamics of a polymer chain confined in a membrane

We present a Brownian dynamics theory with full hydrodynamics (Stokesian dynamics) for a Gaussian polymer chain embedded in a liquid membrane which is surrounded by bulk solvent and walls. The mobility tensors are derived in Fourier space for the two geometries, namely, a free membrane embedded in a bulk fluid, and a membrane sandwiched by the two walls. Within the preaveraging approximation, a new expression for the diffusion coefficient of the polymer is obtained for the free-membrane geometry. We also carry out a Rouse normal mode analysis to obtain the relaxation time and the dynamical structure factor. For large polymer size, both quantities show Zimm-like behavior in the free-membrane case, whereas they are Rouse-like for the sandwiched membrane geometry. We use the scaling argument to discuss the effect of excluded-volume interactions on the polymer relaxation time.     

Eulerian simulation of sedimentation flows in vertical and inclined vessels

Sedimentation of particles in inclined and vertical vessels is numerically simulated using a finite volume method where the Eulerian multiphase model is applied. The particulate phase as well as the fluid phase is regarded as a continuum while the viscosity and solid stress of the particulate phase are modelled by the kinetic theory of granular flows. The numerical results show an interesting phenomenon of the emergence of two circulation vortices of the sedimentation flow in a vertical vessel but only one in the inclined vessel. Several sensitivity tests are simulated to understand the factors that influence the dual-vortex flow structure in vertical sedimentation. Results show that a larger fluid viscosity makes the two vortex centres much closer to each other and the boundary layer effect at lateral walls is the key factor to induce this phenomenon. In the fluid boundary layer particles settle down more rapidly and drag the local carrier fluid to flow downward near the lateral walls and thus form the dual-vortex flow pattern.     

Analytical approximate solutions for two-dimensional viscous flow through expanding or contracting gaps with permeable walls

In this paper, the problem of laminar, isothermal, incompressible and viscous flow in a rectangular domain bounded by two moving porous walls, which enable the fluid to enter or exit during successive expansions or contractions is solved analytically by using the homotopy analysis method (HAM). Graphical results are presented to investigate the influence of the nondimensional wall dilation rate α and permeation Reynolds number Re on the velocity, normal pressure distribution and wall shear stress. The obtained solutions, in comparison with the numerical solutions, demonstrate remarkable accuracy. The present problem for slowly expanding or contracting walls with weak permeability is a simple model for the transport of biological fluids through contracting or expanding vessels.      

Effective interaction between two giant spheres suspended in a size polydisperse hard-sphere fluid

Analytic expressions for the Laplace transform of the interaction energy and force between two exceedingly large hard spheres at infinite dilution in a polydisperse hard-sphere suspending fluid are presented. The equations are based on the Percus–Yevick approximation for the many-component suspending fluid, supplemented by the hypernetted chain approximation for the correlation function of the suspended spheres. By applying the Derjaguin approximation, the energy and force results for two spheres are related to the energy per unit area and the disjoining pressure between two flat walls suspended in a polydisperse fluid. Numerical results for the representative Schultz distributions of the diameters of the species comprising the suspending fluid are presented and discussed.     

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.     

Molecular simulation study of triangle-well fluids confined in slit pores

Grand canonical Monte Carlo simulations are used to study the behaviour of triangle-well (TW) fluids with variable well widths confined inside slit pores. The effect of individual factors influencing the properties of confined fluids such as fluid–fluid interactions, pore size and pore wall–fluid interactions are obtained using simulations as it is difficult to experimentally determine the same. An interesting observation of this study is that inside the narrow pore of slit height h* = 5 at the high-pressure condition of P* = 0.8, for the TW fluid with long-range attraction or for the fluid at a low temperature for even a short-range attraction, the density profiles show layering such that there is a sticking tendency of the particles at centre, while there is a depletion of particles near the wall (as the layers at the centre have higher density peak heights than near the walls).