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.

Evolution of texture and grain boundary microstructure in two-phase (α + β) brass during recrystallization

The evolution of texture and microstructure during recrystallization is studied for two-phase copper alloy (Cu–40Zn) with a variation of the initial texture and microstructure (hot rolled and solution treated) as well as the mode of rolling (deformation path: uni-directional rolling and cross rolling). The results of bulk texture have been supported by micro-texture and microstructure studies carried out using electron back scatter diffraction (EBSD). The initial microstructural condition as well as the mode of rolling has been found to alter the recrystallization texture and microstructure. The uni-directionally rolled samples showed a strong Goss and BR {236}385 component while a weaker texture similar to that of rolling evolved for the cross-rolled samples in the α phase on recrystallization. The recrystallization texture of the β phase was similar to that of the rolling texture with discontinuous 101 α and {111} γ fiber with high intensity at {111}101. For a given microstructure, the cross-rolled samples showed a higher fraction of coincident site lattice Σ3 twin boundaries in the α phase. The higher fraction of Σ3 boundaries is explained on the basis of the higher propensity of growth accidents during annealing of the cross-rolled samples. The present investigation demonstrates that change in strain path, as introduced during cross-rolling, could be a viable tool for grain boundary engineering of low SFE fcc materials.

Eulerian two-phase flow simulation and experimental validation of semisolid slurry generation process using cooling slope

Experimental and numerical studies of slurry generation using a cooling slope are presented in the paper. The slope having stainless steel body has been designed and constructed to produce semisolid A356 Al alloy slurry. The pouring temperature of molten metal, slope angle of the cooling slope and slope wall temperature were varied during the experiment. A multiphase numerical model, considering liquid metal and air, has been developed to simulate the liquid metal flow along the cooling channel using an Eulerian two-phase flow approach. Solid fraction evolution of the solidifying melt is tracked at different locations of the cooling channel following Schiel's equation. The continuity, momentum and energy equations are solved considering thin wall boundary condition approach. During solidification of the melt, based on the liquid fraction and latent heat of the alloy, temperature of the alloy is modified continuously by introducing a modified temperature recovery method. Numerical simulations has been carried out for semisolid slurry formation by varying the process parameters such as angle of the cooling slope, cooling slope wall temperature and melt superheat temperature, to understand the effect of process variables on cooling slope semisolid slurry generation process such as temperature distribution, velocity distribution and solid fraction of the solidifying melt. Experimental validation performed for some chosen cases reveals good agreement with the numerical simulations.

Two phase microstructure for Ag-Ni nanowires

In the present study, electrodeposition technique was used to produce Ag-Ni nanowires. Ag-Ni system shows extremely high bulk immiscibility. Nanowire morphology was achieved by employing an anodic alumina membrane having pores of similar to 200 nm diameter. Microstructure of as-deposited wire was composed of nano-sized solid solution structured Ag-Ni nanoparticles embedded in a matrix of pure Ag phase. It is proposed that the two phase microstructure resulted from an initial formation of solid solution structured nanoparticles in the alumina template pore followed by nucleation of pure Ag phase over the particles which eventually grew to form the matrix phase. (C) 2013 Elsevier B. V. All rights reserved.

CFD Analysis of two phase flow with a prediction of torque demand of a Bi-lobe rotor of a mixing chamber

An attempt to study the fluid dynamic behavior of two phase flow comprising of solid and liquid with nearly equal density in a geometrical case that has an industrial significance in theareas like processing of polymers, food, pharma ceutical, paints. In this work,crystalline silica is considered as the dispersed medium in glycerin. In the CFD analysis carried out,the two phase components are considered to be premixed homogeneously at the initial state. The flow in a cylinder that has an axially driven bi-lobe rotor, a typical blender used in polymer industry for mixing or kneading to render the multi-component mixture to homogeneous condition is considered. A viscous, incompressible, isothermal flow is considered with an assumption that the components do not undergo any physical change and the solids are rigid and mix in fully wetting conditions. Silica with a particle diameter of 0.4 mm is considered and flow is analyzed for different mixing fractions. An industry standard CFD code is used for solving 3D-RANS equations. As the outcome of the study the torque demand by the bi-lobe rotor for different mixture fractions which are estimated show a behavioral consistency to the expected physical phenomena occurring in the domain considered.

Thermoelectric Properties of Two-Phase PbTe with Indium Inclusions

The thermoelectric figure of merit (zT) can be increased by introduction of additional interfaces in the bulk to reduce the thermal conductivity. In this work, PbTe with a dispersed indium (In) phase was synthesized by a matrix encapsulation technique for different In concentrations. x-Ray diffraction analysis showed single-phase PbTe with In secondary phase. Rietveld analysis did not show In substitution at either the Pb or Te site, and this was further confirmed by room-temperature Raman data. Low-magnification (similar to 1500x) scanning electron microscopy images showed micrometer-sized In dispersed throughout the PbTe matrix, while at high magnification (150,000x) an agglomeration of PbTe particles in the hot-pressed samples could be seen. The electrical resistivity (rho) and Seebeck coefficient (S) were measured from 300 K to 723 K. Negative Seebeck values showed all the samples to be n-type. A systematic increase in resistivity and higher Seebeck coefficient values with increasing In content indicated the role of PbTe-In interfaces in the scattering of electrons. This was further confirmed by the thermal conductivity (kappa), measured from 423 K to 723 K, where a greater reduction in the electronic as compared with the lattice contribution was found for In-added samples. It was found that, despite the high lattice mismatch at the PbTe-In interface, phonons were not scattered as effectively as electrons. The highest zT obtained was 0.78 at 723 K for the sample with the lowest In content.

The X-ray photoelectron spectroscopy and high-temperature structural studies of Zn1-xNixO/NiO two-phase composites

Detailed investigation of the chemical states and local atomic environment of Ni and Zn in the two-phase composites of Zn1-xNixO/NiO was reported. The X-ray photoelectron spectra of both Ni-2p and Zn-2p revealed the existence of a doublet with spin-orbit splitting approximate to 17.9 and 23.2eV, respectively confirming the divalent oxidation state of both Ni and Zn. However, the samples fabricated under oxygen-rich conditions exhibit significant difference in the binding energy approximate to 18.75eV between the 2p3/2 and 2p1/2 states of Ni. The shift in the satellite peaks of Ni-2p with increasing the Ni composition x within the Zn1-xNixO/NiO matrix signifies the attenuation of nonlocal screening because of reduced site occupancy of two adjacent Zn ions. The temperature dependence of X-ray diffraction analysis reveals a large distortion in the axial-rhombohedral angle for oxygen-rich NiO. Conversely, no significant distortion was noticed in the NiO system present as a secondary phase within Zn1-xNixO. Nevertheless, the unit-cell volume of both wurtzite h.c.p. Zn1-xNixO and f.c.c. NiO exhibits an anomalous behavior between 150 and 300 degrees C. The origin of such unusual change in the unit-cell volume was discussed in terms of oxygen stoichiometry.