Verifying scalings for bending rigidity of bilayer membranes using mesoscale models

The bending rigidity kappa of bilayer membranes was studied with coarse grained soft repulsive potentials using dissipative particle dynamics (DPD) simulations. Using a modified Andersen barostat to maintain the bilayers in a tensionless state, the bending rigidity was obtained from a Fourier analysis of the height fluctuations. From simulations carried out over a wide range of membrane thickness, the continuum scaling relation kappa proportional to d(2) was captured for both the L-alpha and L-beta phases. For membranes with 4 to 6 tail beads, the bending rigidity in the L-beta phase was found to be 10-15 times higher than that observed for the L-alpha phase. From the quadratic scalings obtained, a six fold increase in the area stretch modulus, k(A) was observed across the transition. The magnitude of increase in both kappa and k(A) from the L-alpha to the L-beta phase is consistent with current experimental observations in lipid bilayers and to our knowledge provides for the first time a direct evaluation of the mechanical properties in the L-beta phase.

Methanation of CO2 over Ru?SiO2 catalyst

Using Ru - SiO2 catalyst, the kinetics of methanation of carbon dioxide has been studied. In the temperature range of 320-460-degrees-C a simple power law model is found to predict experimental results with a good agreement over the range of variables studied.

Modeling Of Hydrolytic Polymerization In A Semibatch Nylon 6 Reactor

Hydrolytic polymerization of caprolactam to Nylon 6 in a semibatch reactor is carried out by heating a mixture of water and caprolactam. Evaporation of volatiles caused by heating results in a pressure build-up. After the pressure reaches a predetermined value, vapors are vented to keep the pressure constant for some time, and thereafter, to lower the pressure to a value slightly above atmospheric in a preprogrammed manner. The characteristics of the polymer are determined by the chemical reactions and the vaporization of water and caprolactam. The semibatch operation has been simulated and the predictions have been compared with industria data. The observed temperature and pressure histories were predicted with a fair degree of accuracy. It was found that the predictions of the degree of polymerization however are sensitive to the vapor-liquid equilibrium relations. A comparison with an earlier model, which neglected mass transfer resistance, indicates that simulation using the VLE data of Giori and Hayes and accounting for mass transfer resistance is more reliable.

Rapid Synthesis of Ultrahigh Adsorption Capacity Zirconia by a Solution Combustion Technique

Tetragonal ZrO2 was synthesized by the solution combustion technique using glycine as the fuel. The compound was characterized by X-ray diffraction, scanning electron microscopy, X-ray photoelectron spectroscopy, infrared spectroscopy, and BET surface area analysis. The ability of this compound to adsorb dyes was investigated, and the compound had a higher adsorption capacity than commercially activated carbon. Infrared spectroscopic observations were used to determine the various interactions and the groups responsible for the adsorption activity of the compound. The effects of the initial concentration of the dye, temperature, adsorbent concentration, and pH of the solution were studied. The kinetics of adsorption was described as a first-order process, and the relative magnitudes of internal and external mass transfer processes were determined. The equilibrium adsorption was also determined and modeled by a composite Langmuir-Freundlich isotherm.

Solubilities of Benzene Derivatives in Supercritical Carbon Dioxide

The solubilities of benzene derivatives in supercritical carbon dioxide was determined by the saturation method over the pressure range (9.5 to 14.5) MPa. The solubilities were determined at (308 and 313) K for 1-chloro-2,4-dinitrobenzene and (308, 318, and 328) K for m-dinitrobenzene. At 308K, the solubility (in mole fraction) of 1-chloro-2,4-dinitrobenzene varied from (2.83 to 5.88).10(-3) while the solubility of m-dinitrobenzene increased from (2.05 to 5.54).10(-3) as the pressure increased from (9.5 to 14.5) MPa. However, the solubilities of both compounds decreased with increasing temperature. Models based on the solubility parameter and semiempirical models such as the Mendez-Santiago-Teja model, the Gordillo model, and the association model, were used to correlate the experimental solubility data for the benzene derivatives.

Monte Carlo simulation of network formation based on structural fragments in epoxy-anhydride systems

A method combining the Monte Carlo technique and the simple fragment approach has been developed for simulating network formation in amine-catalysed epoxy-anhydride systems. The method affords a detailed insight into the nature and composition of the network, showing the distribution of various fragments. It has been used to characterize the network formation in the reaction of the diglycidyl ester of isophthalic acid with hexahydrophthalic anhydride, catalysed by benzyldimethylamine. Pre-gel properties like number and weight distributions and average molecular weights have been calculated as a function of epoxy conversion, leading to a prediction of the gel-point conversion. Analysis of the simulated network further yields other characteristic properties such as concentration of crosslink points, distribution and concentration of elastically active chains, average molecular weight between crosslinks, sol content and mass fraction of pendent chains. A comparison has been made of the properties obtained through simulation with those predicted by the fragment approach alone, which, however, gives only average properties. The Monte Carlo simulation results clearly show that loops and other cyclic structures occur in the gel. This may account for the differences observed between the results of the simulation and the fragment model in the post-gel phase. Copyright (C) 1996 Elsevier Science Ltd.

Structural dependence of density in high molecular weight esters

The effect of molecular structure on density has been examined in high molecular weight esters (molecular weight 300-900), having varying degrees of branching. Densities were calculated from an empirical equation, which agrees well with the experimental values (error +/-1.5%), irrespective of branching. Since density is related to molecular packing and hence to the molecular rotation, in n-alkanes, the glass transition temperature (T(g)) and density both increase with molecular weight, and hence T(g) is directly related to the density. The esters exhibit a complex behavior. In linear esters the T(g) decreases with molecular weight which is explained from group contribution and molecular interactions. In the +-branched esters, however, T(g) decreases with molecular weight until the molecular weight reaches 600 and increases sharply thereafter. The Y-branched esters show an intermediate behavior. The dependence of T(g) on molecular weight has been explained from the segmental motion.

Evaluation of relative growth limitation due to depletion of glucose and oxygen during fungal growth on a spherical solid particle

A mathematical model for glucose and oxygen consumption, and cell growth during fungal growth on a single solid particle is developed. A moving biofilm is assumed to be present on the surface of the solid particle. Initially only glucose is assumed to be growth limiting and later oxygen transferred from the gas phase on to the biofilm is also assumed to be growth limiting. Glucose is found to be severely growth limiting when assumed to be the only growth limiting factor and its limiting levels far less severe when oxygen limitation is also included. The objective of the model is to gain a better understanding of the mass transfer and relative growth limiting characteristics of glucose and oxygen in fungal growth systems. The results obtained from the model proposed here will be the subject of future work.

Fluctuation-dissipation relation for nonlinear Langevin equations

It is shown that the fluctuation-dissipation theorem is satisfied by the solutions of a general set of nonlinear Langevin equations with a quadratic free-energy functional (constant susceptibility) and field-dependent kinetic coefficients, provided the kinetic coefficients satisfy the Onsager reciprocal relations for the irreversible terms and the antisymmetry relations for the reversible terms. The analysis employs a perturbation expansion of the nonlinear terms, and a functional integral calculation of the correlation and response functions, and it is shown that the fluctuation-dissipation relation is satisfied at each order in the expansion.

Oxidation of Sulfur Dioxide in Aqueous Suspensions of Activated Carbon

Sulfur dioxide in aqueous solutions at low pH levels exists both in molecular SO2(aq) and in hydrolyzed ionic form HSO3-. Experiments indicate that only HSO3- is the reacting species in the oxidation catalyzed by activated carbon, while SO2(aq) deactivates by competing with HSO3 for the active sites of the catalyst particles. A mechanism is proposed and a rate model is developed that also accounts for the effect of sulfuric acid (product of the oxidation) on the solubility of sulfur dioxide. It predicts first order in HSO3-, half order in dissolved oxygen, and a linear deactivation effect of SO2(aq), which are confirmed by experimental data. The deactivation reaches a constant level corresponding to saturation of the active sites by SO2(aq). Activation energy for the oxidation is 93.55 kJ mol(-1) and for deactivation is 21.4 kJ mol(-1).

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Treatment of electroplating effluent through emulsion-free liquid embrane

Separation of dissolved heavy metals such-as Cr(VI) and Cu(II) from electroplating effluents using a new technique of emulsion-free liquid membrane (EFLM) has been studied. Experimental results show that nearly 95% extraction is obtained resulting in stripping phase enrichment up to 50 times relative to feed. It is also found that emulsion-free liquid membranes are highly efficient and superior to other types of liquid membranes.

Gas-phase controlled mass transfer in a foam-bed reactor

Gas-phase controlled absorption of ammonia in foams made of solutions of sulphuric acid has been studied experimentally. Effects of gas-phase concentration of ammonia and type of surfactant on the performance of the foam-bed reactor are investigated. Gas-phase controlled absorption from a spherical bubble is anaylzed using the asymptotic value of Sherwood number (Sh = 6.58), for both negligible as well as significant changes in the volume of the bubble. The experimental data are shown to be in good agreement with the single-stage model of the foam-bed reactor using these asymptotic sub-models, as well as the diffusion-in-sphere analysis available in literature. Influence of effective diffusivity on the time dependence of fractional gas absorption has been found to be unimportant for foam columns with large times of contact. The asymptotic sub-models have been compared and use of the rigid-sphere asymptotic sub-model is recommended for foam columns of practical relevence.