Physical properties of sodium carboxymethyl cellulose molecules adsorbed on a polyacrylonitrile ultrafiltration membrane

Ultrafiltration experiments showed that the graphical relationship between flux and pressure was a straight line through the origin, provided that the wall shear rate of the bulk fluid was higher than a certain critical value or the pressure was below a critical value. A higher critical shear rate corresponded to a higher critical pressure. For these conditions the total hydrodynamic resistance was only slightly greater than the resistance of a clean membrane for pure water. This additional resistance is attributed to a (mono-) molecular layer of macromolecules which is adsorbed on the membrane in the absence of both a concentration polarization layer and a conventional gel layer.At steady state ultrafiltration conditions, an increase of the flux was obtained after replacing the bulk solution by distilled water at constant experimental conditions, which is attributed to the removal of the concentration polarization layer whereas a mono-molecular layer of macromolecules remained adsorbed on the membrane. For these conditions the flux vs. pressure relationship showed a qualitatively similar behaviour as for ultrafiltration conditions.At a constant shear rate the flux vs. pressure relationship was a straight line through the origin for pressures below the critical pressure, the value of which increased with the shear rate. This linear relationship was reversible, showing no hysteresis. However, if the pressure was higher than its critical value, the flux vs. pressure relationship was no longer a straight line as a consequence of the occurrence of an additional hydrodynamic resistance which did not disappear entirely upon lowering the pressure below its critical value. For the explanation of these phenomena it is assumed that freely moveable parts of the adsorbed macromolecules can block the entrance region of the pores in the membrane if the pressure is beyond its critical value.On the other hand, for pressures below the critical pressure or shear rates beyond the critical shear rate, the pores of the membrane are deblocked. This blocking and deblocking of pores by parts of adsorbed macromolecules apparently takes place in a partly reversible way.     

Friction of rodlike particles adsorbed to a planar surface in shear flow

A planar hard surface covered with elongated stiff rodlike particles in shear flow is considered in the low-Reynolds-number regime assuming low particle surface coverage. The particles are modeled as straight chains of spherical beads. Multipole expansion of the Stokes equations (the accurate HYDROMULTIPOLE algorithm) is applied to evaluate the hydrodynamic force exerted by the fluid on the rodlike particles, depending on their shape, i.e., on the number of beads and their orientation with respect to the wall and to the ambient shear flow.     

The interfacial tension and phase diagram of the Widom-Rowlinson mixture via Monte Carlo simulations

Results of Monte Carlo simulations are reported for the interfacial tension between two fluid phases in a binary mixture of penetrable spheres in which molecular pairs of like species do not interact, while those of unlike species interact as hard spheres. Semigrand canonical ensemble Monte Carlo simulations in a cubic cell with periodic boundary conditions are used to obtain histograms for various system sizes at various densities. At a given density, the interfacial tension and compositions of coexisting phases for an infinite system are evaluated via histogram analysis combined with finite-size scaling. The density dependence of the interfacial tension and phase diagram for an infinite system are thus obtained. The simulated behavior of the interfacial tension close to the critical density corroborates previous suggestions that the model belongs to the three-dimensional Ising universality class.     

Effects of pre-oxidation temperature on coal secondary spontaneous combustion

The hydrodynamic force (drag) on spherical and irregularly shaped particles significantly increases when the particles move close to solid and permeable boundaries. The overall effect of the increased hydrodynamic drag is to hinder the particle movement in the vicinity of boundaries and this includes the Brownian movement and electrophoresis. The Monte Carlo simulation method is used to model the Brownian movement, the resulting diffusion, and the electrophoresis of spherical particles in narrow, cylindrical pores, filled with Newtonian fluids. It is observed that the effect of the pore walls is a significant reduction of the space-averaged electrophoretic velocity of the particles, which implies reduced particle flux through the pores. The hindered electrophoresis is primarily a geometric phenomenon, caused by the increased drag and depends on the size of the particles and the pore-to-particle diameter ratio. The temperature of the fluid slightly affects the hindered electrophoresis through its effect on the viscosity, which is a determinant of the Brownian force, the diffusivity and the electrophoretic velocity. The hindered electrophoresis is almost independent of the other fluid and particle properties, such as density. Based on the simulation results a non-linear correlation for the flux of particles is derived, valid in the ranges 5?<?R/α?<?120, 5 nm?<?α?<?100 nm and 273 K?<?T?<?355 K.

     

Polymorphism in simple liquids: a Gibbs ensemble Monte Carlo study

We perform Gibbs ensemble Monte Carlo (GEMC) simulations of a one-component system of hard spheres with a repulsive shoulder and an attractive well. We show the existence of two distinct liquid-gas and liquid-liquid phase equilibria. The GEMC estimate of the critical parameters, as following from an interpolation of the binodal points, is only slightly influenced by finite size effects. The liquid-gas critical temperature and pressure are lower than those of the liquid-liquid phase separation. A discussion of our findings in comparison with those of previous numerical studies is also presented.     

Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension

We study the phase behavior and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive model (YRPM) of oppositely charged hard spheres with diameter sigma using Monte Carlo simulations. We determine the gas-liquid and gas-solid phase transitions using free energy calculations and grand-canonical Monte Carlo simulations for varying inverse Debye screening length kappa. We find that the gas-liquid phase separation is stable for kappasigma相似文献   

The ferroelectric transition of dipolar hard spheres

We investigate by Monte Carlo simulation the size dependence of the variation of the polarization and the dielectric constant with temperature for dipolar hard spheres at the two densities rho sigma3=0.80 and 0.88. From the crossing of the fourth-order cumulant for different system sizes first more precise estimates of the ferroelectric transition temperatures are obtained. Theoretical approaches, when predicting an ordering transition, are shown to generally overestimate the critical temperature.     

Nonequilibrium Monte Carlo simulation of lattice block copolymer chains subject to oscillatory shear flow

This paper has extended nonequilibrium Monte Carlo (MC) approach to simulate oscillatory shear flow in a lattice block copolymer system. Phase transition and associated rheological behaviors of multiple self-avoiding chains have been investigated. Stress tensor has been obtained based upon sampled configuration distribution functions. At low temperatures, micellar structures have been observed and the underlying frequency-dependent rheological properties exhibit different initial slopes. The simulation outputs are consistent with the experimental observations in literature. Chain deformation during oscillatory shear flow has also been revealed. Although MC simulation cannot account for hydrodynamic interaction, the highlight of our simulation approach is that it can, at small computing cost, investigate polymer chains simultaneously at different spatial scales, i.e., macroscopic rheological behaviors, mesoscopic self-assembled structures, and microscopic chain configurations.     

Study of the ST2 model of water close to the liquid-liquid critical point

We perform successive umbrella sampling grand canonical Monte Carlo computer simulations of the original ST2 model of water in the vicinity of the proposed liquid-liquid critical point, at temperatures above and below the critical temperature. Our results support the previous work of Y. Liu, A. Z. Panagiotopoulos and P. G. Debenedetti [J. Chem. Phys., 2009, 131, 104508], who provided evidence for the existence and location of the critical point for ST2 using the Ewald method to evaluate the long-range forces. Our results therefore demonstrate the robustness of the evidence for critical behavior with respect to the treatment of the electrostatic interactions. In addition, we verify that the liquid is equilibrated at all densities on the Monte Carlo time scale of our simulations, and also that there is no indication of crystal formation during our runs. These findings demonstrate that the processes of liquid-state relaxation and crystal nucleation are well separated in time. Therefore, the bimodal shape of the density of states, and hence the critical point itself, is a purely liquid-state phenomenon that is distinct from the crystal-liquid transition.     

Elasticity of two-dimensional crystals of polydisperse hard disks near close packing: surprising behavior of the Poisson's ratio

The equation of state, elastic constants, and Poisson's ratio of a crystalline two-dimensional polydisperse hard disk system were determined in the close packing limit. Monte Carlo simulations in the NpT ensemble with variable shape of the periodic box reveal that the pressure and elastic constants grow with increasing polydispersity. The equation of state and the bulk modulus are well described by the free volume approximation. The latter approximation fails, however, for the shear modulus. The simulations also show that the introduction of any amount of size polydispersity in the hard disk systems causes a discontinuous "jump" of the Poisson's ratio in the close packing limit from the value ν(δ=0) = 0.1308(22), obtained for equidiameter hard disks, to ν(δ>0) ≈ 1, estimated for the polydisperse disks.     

Modified free volume theory of self-diffusion and molecular theory of shear viscosity of liquid carbon dioxide

In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site-site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.     

On the behavior of the electrostatic colloidal interaction in the membrane filtration of latex suspensions

In order to investigate effects of the colloidal interaction in the membrane filtrations, the dead-end ultrafiltration of latex colloids was conducted with fully retentive membranes. Experimental results concerning the permeate flux during the filtration indicate that the void fraction of cake layer increased with the decrease of the ionic strength, due to the expanded Debye double layer thickness around the particles. The concentration dependence of the gradient diffusion coefficient of colloidal particles has been examined as a function of solution ionic strength. The NVT Monte Carlo simulation was applied on the bulk suspension so as to determine the thermodynamic coefficient, and the hydrodynamic coefficient was evaluated from the previously developed relation for an ordered system. The long-range electrostatic interactions between the particles are determined by using a singularity method, which provides accurate solutions to the linearized electrostatic field. The predictions on the variation of concentration polarization layer have been presented, from which we found that both the permeate flux and the particle diffusion are related to determine the concentration distribution above the cake layer.     

Chemical and physical aspects of organic fouling of forward osmosis membranes

The growing attention to forward osmosis (FO) membrane processes from various disciplines raises the demand for systematic research on FO membrane fouling. This study investigates the role of various physical and chemical interactions, such as intermolecular adhesion forces, calcium binding, initial permeate flux, and membrane orientation, in organic fouling of forward osmosis membranes. Alginate, bovine serum albumin (BSA), and Aldrich humic acid (AHA) were chosen as model organic foulants. Atomic force microscopy (AFM) was used to quantify the intermolecular adhesion forces between the foulant and the clean or fouled membrane in order to better understand the fouling mechanisms. A strong correlation between organic fouling and intermolecular adhesion was observed, indicating that foulant–foulant interaction plays an important role in determining the rate and extent of organic fouling. The fouling data showed that FO fouling is governed by the coupled influence of chemical and hydrodynamic interactions. Calcium binding, permeation drag, and hydrodynamic shear force are the major factors governing the development of a fouling layer on the membrane surface. However, the dominating factors controlling membrane fouling vary from foulant to foulant. With stronger intermolecular adhesion forces, hydrodynamic conditions for favorable foulant deposition leading to cake formation are more readily attained. Before a compact cake layer is formed, the fouling rate is affected by both the intermolecular adhesion forces and hydrodynamic conditions. However, once the cake layer forms, all three foulants have very similar flux decline rates, and further changes in hydrodynamic conditions do not influence fouling behavior.     

Orientational ordering in hard rectangles: The role of three-body correlations

We investigate the effect of three-body correlations on the phase behavior of hard rectangle two-dimensional fluids. The third virial coefficient B3 is incorporated via an equation of state that recovers scaled particle theory for parallel hard rectangles. This coefficient, a functional of the orientational distribution function, is calculated by Monte Carlo integration, using an accurate parametrized distribution function, for various particle aspect ratios in the range of 1-25. A bifurcation analysis of the free energy calculated from the obtained equation of state is applied to find the isotropic (I)-uniaxial nematic (N(u)) and isotropic-tetratic nematic (N(t)) spinodals and to study the order of these phase transitions. We find that the relative stability of the N(t) phase with respect to the isotropic phase is enhanced by the introduction of B3. Finally, we have calculated the complete phase diagram using a variational procedure and compared the results with those obtained from scaled particle theory and with Monte Carlo simulations carried out for hard rectangles with various aspect ratios. The predictions of our proposed equation of state as regards the transition densities between the isotropic and orientationally ordered phases for small aspect ratios are in fair agreement with simulations. Also, the critical aspect ratio below which the N(t) phase becomes stable is predicted to increase due to three-body correlations, although the corresponding value is underestimated with respect to simulation.     

GEMC和GDI方法计算流体气液相平衡的比较

考察采用TraPPE联合原子和OPLS全原子力场两种分子力场, Gibbs系综蒙特卡罗(GEMC)方法和Gibbs-Duhem积分(GDI)方法计算流体气液相平衡的适用性、计算速度、计算精度等问题. 结果表明, 在采用全原子力场情况下, GDI方法比GEMC方法极大地节省了计算时间. 从计算结果来看, 两种方法各有适用范围, 在使用时可互为补充. 在给定力场的前提下, 两种方法所得到的液相密度、蒸发焓、临界温度和临界密度相差不大, 而当力场中的缺陷导致蒸发焓的计算不够准确时, 两种计算方法得到的气体的压力和密度明显不同,进而导致预测的临界压力也明显不同.     

Orientational order in high density dipolar hard sphere fluids

Taking advantage of recent estimates, by one of us, of the critical temperature of the isotropic-ferroelectric transition of high density dipolar hard spheres, we performed new Monte Carlo simulations in the close vicinity of these estimates and applied histogram reweighting methods to obtain refined values of the critical temperatures from the crossing of the fourth-order cumulant for different system sizes. The ferroelectric line is determined in the density range rho*=0.80-0.95, and the onset of columnar ordering is located.     

Effect of Hydrodynamic Interaction on Flow-induced Polymer Translocation Through a Nanotube

We investigated the effect of hydrodynamic interaction(HI) on flow-induced polymer translocation through a nanotube by Brownian dynamics simulations. Whether there is HI in the simulation system is separately controlled by using different diffusion tensors. It is found that HI has no effect on critical velocity flux for long polymer chains due to the competition between more drag force and the hindrance of chain stretching from HI, however, HI broadens the transition interval. In addition, for flow-induced polymer translocation with HI, the critical velocity flux firstly slowly decreases with the increase of chain length and then becomes identical to that of it without HI, that is, the critical velocity flux is independent of chain length. At the same time, HI also accelerates the translocation process and makes the relative variation amplitude of single bead translocation time smaller. In fact, HI can enhance the intrachain cooperativity to make the whole chain obtain more drag force from fluid field and hinder chain stretching, both of which play an important role in translocation process.     

Shear-induced aggregation of carboxylated polymer latices

Orthokinetic flocculation of latex particles containing high levels of surface carboxylation was studied by shearing for short times in a rotational viscometer and measuring changes in particle size using photon correlation spectroscopy. Results can be interpreted in terms of a comparison between the DLVO repulsive force between a pair of particles and the hydrodynamic shear force opposing it. This enables a prediction of the critical shear rate, , required to initiate aggregation. Experimental and calculated values of showed good agreement.The stability of the carboxylated latices under shear was much reduced at pH values >7 when the surface groups are ionised. The increase in suspension viscosity with pH was shown to be critical in determining the onset of aggregation via hydrodynamic forces.     

A simple approach to treat anisotropic elastic collisions in Monte Carlo calculations of the electron energy distribution function in cold plasmas

A simplified method allowing one to treat anisotropic electron heavy species elastic scattering in MonteCarlo models of gas discharges with the proper value for collision frequency is proposed The method is applied to an electric discharge in a Ne · Xe/HCl mixture, and the results are compared with the solution of the two-term expansion of the Boltzmann equation under the same conditions. Methods for reduction of computational time in Monte Carlo codes and the use of the Monte Carlo flux method are also discussed.     

Entropic forces and directed alignment of hard squares in suspensions of rods and disks

We use Monte Carlo simulations in two dimensions to study the depletion forces between two hard squares in a suspension of hard rods or disks. We determine the effects of size and concentration of rods and disks on the potential of mean force between the squares. Both rods and disks produce a short-range depletion attraction between the two squares. The depletion interaction can be strong enough to outweigh the (rotational) entropic repulsion between the squares at certain sizes and concentrations of the rods and disks. We also probe the relative orientation that two squares adopt as they approach each other and we observe rich behavior, in which the relative orientation depends on the size, concentration, and shape of the depletion agent. Simple models based on the ideas of Asakura and Oosawa [J. Chem. Phys. 22, 1255 (1954)] can explain trends in the potentials of mean force obtained from the simulations.