Arbitrary Lagrangian-Eulerian finite-element method for computation of two-phase flows with soluble surfactants

A finite-element scheme based on a coupled arbitrary Lagrangian-Eulerian and Lagrangian approach is developed for the computation of interface flows with soluble surfactants. The numerical scheme is designed to solve the time-dependent Navier-Stokes equations and an evolution equation for the surfactant concentration in the bulk phase, and simultaneously, an evolution equation for the surfactant concentration on the interface. Second-order isoparametric finite elements on moving meshes and second-order isoparametric surface finite elements are used to solve these equations. The interface-resolved moving meshes allow the accurate incorporation of surface forces, Marangoni forces and jumps in the material parameters. The lower-dimensional finite-element meshes for solving the surface evolution equation are part of the interface-resolved moving meshes. The numerical scheme is validated for problems with known analytical solutions. A number of computations to study the influence of the surfactants in 3D-axisymmetric rising bubbles have been performed. The proposed scheme shows excellent conservation of fluid mass and of the total mass of the surfactant. (C) 2012 Elsevier Inc. All rights reserved.

Two-Phase Flow With Phase Change in Porous Channels

Phase change heat transfer in porous media finds applications in various geological flows and modern heat pipes. We present a study to show the effect of phase change on heat transfer in a porous channel. We show that the ratio of Jakob numbers based on wall superheat and inlet fluid subcooling governs the liquid-vapor interface location in the porous channel and below a critical value of the ratio, the liquid penetrates all the way to the extent of the channel in the flow direction. In such cases, the Nusselt number is higher due to the proximity of the liquid-vapor interface to the heat loads. For higher heat loads or lower subcooling of the liquid, the liquid-vapor interface is pushed toward the inlet, and heat transfer occurs through a wider vapor region thus resulting in a lower Nusselt number. This study is relevant in the designing of efficient two-phase heat exchangers such as capillary suction based heat pipes where a prior estimation of the interface location for the maximum heat load is required to ensure that the liquid-vapor interface is always inside the porous block for its operation.

Evolution of crystallographic texture and microstructure in the orthorhombic phase of a two-phase alloy Ti-22Al-25Nb

Evolution of crystallographic texture in the orthorhombic phase of a two-phase alloy Ti–22Al–25Nb (at%), consisting of orthorhombic (O) and bcc (β/B2) phases, was studied. The material was subjected to deformation in two-phase field as well as in the single β phase field. The resulting evolution of microstructure and crystallographic texture were recorded using scanning electron microscopy and X-ray diffraction. The orthorhombic phase underwent change in morphology (from platelets to equiaxed) on rolling in the two-phase field with the texture getting sharper with the amount of deformation. Rolling above β transus temperature led to hot deformation of single β phase microstructure and its subsequent cooling produced transformed coarse platelets of orthorhombic phase with texture in orientation relation with the high temperature deformed β phase.

Theoretical study of the optimum performance of a two-phase flow CO2 gasdynamic laser

Theoretical optimization studies of the performance of a combustion driven premixed two-phase flow gasdynamic laser are presented. The steady inviscid nonreacting quasi-one-dimensional two-phase flow model including appropriate finite rate vibrational kinetic rates has been used in the analysis. The analysis shows that the effect of the particles on the optimum performance of the two-phase laser is very small. The results are presented in graphical form. Applied Physics Letters is copyrighted by The American Institute of Physics.

Analysis of EXCOBULLE two-phase expansion tests

analysis of a complex physical problem and the close agreement they achieved with observations. However, the following points need to be clarified. First of all the authors assume that during the initial phases of expansion, the Tayior's instability sets in due to the acceleraacceleration of lighter fluid against the more dense cold water.

A new kinetic model for bulk polymerization of vinyl chloride based on two-phase hypothesis

A kinetic model has been developed for the bulk polymerization of vinyl chloride using Talamini's hypothesis of two-phase polymerization and a new concept of kinetic solubility which assumes that rapidly growing polymer chains have considerably greater solubility than the thermodynamic solubility of preformed polymer molecules of the same size and so can remain in solution even under thermodynamically unfavourable conditions. It is further assumed that this kinetic solubility is a function of chain length. The model yields a rate expression consistent with the experimental data for vinyl chloride bulk polymerization and moreover is able to explain several characteristic kinetic features of this system. Application of the model rate expression to the available rate data has yielded 2.36 × 108l mol−1 sec−1 for the termination rate constant in the polymer-rich phase; as expected, this value is smaller than that reported for homogenous polymerization by a factor of 10–30.

Two-Phase Thermodynamic Model for Efficient and Accurate Absolute Entropy of Water from Molecular Dynamics Simulations

Presented here is the two-phase thermodynamic (2PT) model for the calculation of energy and entropy of molecular fluids from the trajectory of molecular dynamics (MD) simulations. In this method, the density of state (DoS) functions (including the normal modes of translation, rotation, and intramolecular vibration motions) are determined from the Fourier transform of the corresponding velocity autocorrelation functions. A fluidicity parameter (f), extracted from the thermodynamic state of the system derived from the same MD, is used to partition the translation and rotation modes into a diffusive, gas-like component (with 3Nf degrees of freedom) and a nondiffusive, solid-like component. The thermodynamic properties, including the absolute value of entropy, are then obtained by applying quantum statistics to the solid component and applying hard sphere/rigid rotor thermodynamics to the gas component. The 2PT method produces exact thermodynamic properties of the system in two limiting states: the nondiffusive solid state (where the fluidicity is zero) and the ideal gas state (where the fluidicity becomes unity). We examine the 2PT entropy for various water models (F3C, SPC, SPC/E, TIP3P, and TIP4P-Ew) at ambient conditions and find good agreement with literature results obtained based on other simulation techniques. We also validate the entropy of water in the liquid and vapor phases along the vapor-liquid equilibrium curve from the triple point to the critical point. We show that this method produces converged liquid phase entropy in tens of picoseconds, making it an efficient means for extracting thermodynamic properties from MD simulations.

Effect of mode of rolling on development of texture and microstructure in two-phase (alpha plus beta) brass

The evolution of microstructure and texture during deformation of two-phase (alpha + beta) brass was studied for different initial microstructure and texture. The deformation processing involved unidirectional and multi-step cross-rolling. The bulk textures were determined by measuring the pole figures and calculating the orientation distribution function ODF for both alpha (fcc) and beta (bcc) phases, while the microstructure and other microstructural parameters were measured through optical microscopy and scanning electron microscopy with electron back scatter diffraction (SEM-EBSD). Results indicate that textures developed after unidirectional rolling and multi-step cross-rolling are significantly different. The variation in initial texture had a pronounced effect on the development of texture in the alpha phase during subsequent deformation. (C) 2010 Elsevier B.V. All rights reserved.

Potentiometric determination of activities in the two-phase fields of the system Na2O---(α)Al2O3

Three compounds have been found to be stable in the pseudobinary system Na2O---(α)Al2O3 between 825 and 1400 K; two nonstoichiometric phases, β-alumina and β″-alumina, and NaAlO2. The homogeneity of β-alumina ranges from 9.5 to 11 mol% Na2O, while that of β″-alumina from 13.3 to 15.9 mol% Na2O at 1173 K. The activity of Na2O in the two-phase fields has been determined by a solid-state potentiometric technique. Since both β- and β″-alumina are fast sodium ion conductors, biphasic solid electrolyte tubes were used in these electrochemical measurements. The open circuit emf of the following cells were measured from 790 to 980 K: [GRAPHICS] The partial molar Gibbs' energy of Na2O relative to gamma-Na2O in the two-phase regions can be represented as: DELTA-GBAR(Na2O)(alpha- + beta-alumina) = -270,900 + 24.03 T, DELTA-GBAR(Na2O)(beta- + beta"-alumina) = -232,700 + 56.19 T, and DELTA-GBAR(Na2O)(beta"-alumina + NaAlO2) = -13,100 - 4.51 T J mol-1. Similar galvanic cells using a Au-Na alloy and a mixture of Co + CoAl(2+2x)O4+3x + (alpha)Al2O3 as electrodes were used at 1400 K. Thermodynamic data obtained in these studies are used to evaluate phase relations and partial pressure of sodium in the Na2O-(alpha) Al2O3 system as a function of oxygen partial pressure, composition and temperature.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Phase nonequilibrium effects on the gain of a two-phase flow gasdynamic laser

Analysis of gas-particle nozzle flow is carried out with attention to the effect of dust particles on the vibrational relaxation phenomena and consequent effects on the gain of a gasdynamic laser. The phase nonequilibrium between the gas mixture and the particles during the nozzle expansion process is taken into account simultaneously. The governing equations of the two-phase nozzle flow have been transformed into similar form, and general correlating parameters have been obtained. It is shown from the present analysis that the particles present in the mixture affect the optimum gain obtainable from a gasdynamic laser adversely, and the effect depends on the size and loading of the particles in the mixture.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.

The Behaviour of Two-Phase Flow of DNAPL and Water through a Fractured Rock under Confining Pressure

This study presents the characterization of DNAPL and water flow in a fracture under confining pressure. A comprehensive mathematical model and the conditions under which DNAPL will enter an initially water-saturated deforming rock fracture are discussed. A numerical model with which to predict the quantity of each phase in terms of their saturations in deforming rock joint is developed. The effect of varying confining stresses on the traverse time of DNAPL across a fractured aquitard is studied. The sensitivity analysis for physical and hydraulic properties like initial fracture apertures, fracture dips, equivalent fracture aperture and confining pressures are performed and discussed.

Texture and grain boundary character distribution during Equal Channel Angular Extrusion of some two-phase copper alloys

In the present work, a thorough investigation of evolution of microstructure and texture has been carried out to elucidate the evolution of texture and grain boundary character distribution (GBCD) during Equal Channel Angular Extrusion (ECAE) of some model two-phase materials, namely Cu-0.3Cr and Cu-40Zn. Texture of Cu-0.3Cr alloy is similar to that reported for pure copper. On the other hand, in Cu-40Zn alloy, texture evolution in α and β (B2) phases are interdependent. In Cu-0.3Cr alloy, there is a considerable decreases in volume fraction of low angle boundaries (LAGBs), only a slight increase in CSL boundaries, but increase in high angle grain boundaries (HAGBs) from 1 pass to 4 passes for both the routes. In the case of Cu-40Zn alloy, there is an appreciable increase in CSL volume fraction.

Sulphur potential measurements with a two phase sulphide oxide electrolyte

The open circuit potentials of the galvanic cell,Pt (or Au)¦(Ar + H2S + H2)primeparCaS + ZrO2(CaO)par (Ar + H2S+ H2)Prime£t (or Au) has been measured in the temperature range 1000 to 1660 K and PH2S:PH 2 ratios from 1.73×10–5 to 2.65×10–1. The solid electrolyte consists of a dispersion of calcium sulphide in a matrix of calcia-stabilized zirconia. The surface of the electrolyte is coated with a thin layer of calcium sulphide to prevent the formation of water vapour by reaction of hydrogen sulphide with calcium oxide or zirconia present in the electrolyte. The use of a lsquopoint electrodersquo with a catalytically active tip was necessary to obtain steady emfs. At low temperatures and high sulphur potentials the emfs agreed with the Nernst equation. Deviations were observed at high temperatures and low sulphur potentials, probably due to the onset of significant electronic conduction in the oxide matrix of the electrolyte. The values of oxygen and sulphur potentials at which the electronic conductivity is equal to ionic conductivity in the two-phase electrolyte have been evaluated from the emf response of the cell. The sulphide-oxide electrolyte is unsuitable for sulphur potential measurements in atmospheres with high oxygen potentials, where oxidation of calcium sulphide may be expected.