A framework for hybrid simulations of two-phase flow in porous media

Rough a global coarse problem. Although these techniques are usually employed for problems in which the fine-scale processes are described by Darcy's law, they can also be applied to pore-scale simulations and used as a mathematical framework for hybrid methods that couples a Darcy and pore scales. In this work, we consider a pore-scale description of fine-scale processes. The Navier-Stokes equations are numerically solved in the pore geometry to compute the velocity field and obtain generalized permeabilities. In the case of two-phase flow, the dynamics of the phase interface is described by the volume of fluid method with the continuum surface force model. The MsFV method is employed to construct an algorithm that couples a Darcy macro-scale description with a pore-scale description at the fine scale. The hybrid simulations results presented are in good agreement with the fine-scale reference solutions. As the reconstruction of the fine-scale details can be done adaptively, the presented method offers a flexible framework for hybrid modeling.

Effect of NaCl additive on solute-solvent interactions in aqueous polyethylene glycol-sodium sulfate two-phase systems

Abstract Partition behavior of eight small organic compounds and six proteins was examined in poly(ethylene glycol)-8000-sodium sulfate aqueous two-phase systems containing 0.215 M NaCl and 0.5 M osmolyte (sorbitol, sucrose, TMAO) and poly(ethylene glycol)-10000-sodium sulfate-0.215 M NaCl system, all in 0.01 M sodium phosphate buffer, pH 6.8. The differences between the solvent properties of the coexisting phases (solvent dipolarity/polarizability, hydrogen bond donor acidity, and hydrogen bond acceptor basicity) were characterized with solvatochromic dyes using the solvatochromic comparison method. Differences between the electrostatic properties of the phases were determined by analysis of partitioning of sodium salts of dinitrophenylated (DNP-) amino acids with aliphatic alkyl side-chain. The partition coefficients of all compounds examined (including proteins) were described in terms of solute-solvent interactions. The results obtained in the study show that solute-solvent interactions of nonionic organic compounds and proteins in polyethylene glycol-sodium sulfate aqueous two-phase system change in the presence of NaCl additive.

Effect of sodium chloride on solutesolvent interactions in aqueous polyethylene glycolsodium sulfate two-phase systems

Partition behavior of eight small organic compounds and six proteins was examined in poly(ethylene glycol)-8000sodium sulfate aqueous two-phase systems containing 0.215 M NaCl and 0.5 M osmolyte (sorbitol, sucrose, TMAO) and poly(ethylene glycol)-10000sodium sulfate0.215 M NaCl system, all in 0.01 M sodium phosphate buffer, pH 6.8. The differences between the solvent properties of the coexisting phases (solvent dipolarity/polarizability, hydrogen bond donor acidity, and hydrogen bond acceptor basicity) were characterized with solvatochromic dyes using the solvatochromic comparison method. Differences between the electrostatic properties of the phases were determined by analysis of partitioning of sodium salts of dinitrophenylated (DNP-) amino acids with aliphatic alkyl side-chain. The partition coefficients of all compounds examined (including proteins) were described in terms of solutesolvent interactions. The results obtained in the study show that solutesolvent interactions of nonionic organic compounds and proteins in polyethylene glycolsodium sulfate aqueous two-phase system change in the presence of NaCl additive.

Interrelationship between partition behavior of organic compounds and proteins in aqueous dextran-Polyethylene glycol and Polyethylene glycol-sodium sulfate two-phase systems

Partition behavior of adenosine and guanine mononucleotides was examined in aqueous dextran-polyethylene glycol (PEG) and PEG-sodium sulfate two-phase systems. The partition coefficients for each series of mononucleotides were analyzed as a functions of the number of phosphate groups and found to be dependent on the nature of nucleic base and on the type of \ATPS\ utilized. It was concluded that an average contribution of a phosphate group into logarithm of partition coefficient of a mononucleotide cannot be used to estimate the difference between the electrostatic properties of the coexisting phases of ATPS. The data obtained in this study were considered together with those for other organic compounds and proteins reported previously, and the linear interrelationship between logarithms of partition coefficients in dextran-PEG, PEG-Na2SO4 and PEG-Na2SO4-0.215 M NaCl (all in 0.01 M Na- or K/Na-phosphate buffer, pH 7.4 or 6.8) was established. Similar relationship was found for the previously reported data for proteins in Dex-PEG, PEG-600-Na2SO4, and PEG-8000-Na2SO4 ATPS. It is suggested that the linear relationships of the kind established in \ATPS\ may be observed for biological properties of compounds as well.

On a size-structured two-phase population model with infinite states-at-birth

In this work we introduce and analyze a linear size-structured population model with infinite states-at-birth. We model the dynamics of a population in which individuals have two distinct life-stages: an “active” phase when individuals grow, reproduce and die and a second “resting” phase when individuals only grow. Transition between these two phases depends on individuals’ size. First we show that the problem is governed by a positive quasicontractive semigroup on the biologically relevant state space. Then we investigate, in the framework of the spectral theory of linear operators, the asymptotic behavior of solutions of the model. We prove that the associated semigroup has, under biologically plausible assumptions, the property of asynchronous exponential growth.

Spatiotemporal adaptive multiscale multiphysics simulations of two-phase flow

We present a spatiotemporal adaptive multiscale algorithm, which is based on the Multiscale Finite Volume method. The algorithm offers a very efficient framework to deal with multiphysics problems and to couple regions with different spatial resolution. We employ the method to simulate two-phase flow through porous media. At the fine scale, we consider a pore-scale description of the flow based on the Volume Of Fluid method. In order to construct a global problem that describes the coarse-scale behavior, the equations are averaged numerically with respect to auxiliary control volumes, and a Darcy-like coarse-scale model is obtained. The space adaptivity is based on the idea that a fine-scale description is only required in the front region, whereas the resolution can be coarsened elsewhere. Temporal adaptivity relies on the fact that the fine-scale and the coarse-scale problems can be solved with different temporal resolution (longer time steps can be used at the coarse scale). By simulating drainage under unstable flow conditions, we show that the method is able to capture the coarse-scale behavior outside the front region and to reproduce complex fluid patterns in the front region.

Partitioning the relative contribution of one-phase and two-phase seed dispersal when evaluating seed dispersal effectiveness

Seed dispersal effectiveness (SDE) is a conceptual framework that aims at quantifying the contribution of seed dispersal vectors to plant fitness. While it is well recognized that diplochorous dispersal systems, characterized by two successive dispersal steps performed by two different vectors (Phase I=primary seed dispersal and Phase II=secondary seed dispersal) which are common in temperate and tropical regions, little attention has been given to distinguishing the relative contribution of one-phase and two-phase dispersal to overall SDE. This conceptual gap probably results from the lack of a clear methodology to include Phase II dispersal into the calculation of SDE and to quantify its relative contribution. We propose a method to evaluate the relative contribution of one-phase and two-phase dispersal to SDE and determine whether two seed dispersers are better than one. To do so, we used the SDE landscape and an extension of the SDE landscape, the Phase II effect landscape, which measures the direction and magnitude of the Phase II dispersal effect on overall SDE. We used simulated and empirical data from a diplochorous dispersal system in the Peruvian Amazon to illustrate this new approach. Our approach provides the relative contribution of one-phase SDE (SDE1) and two-phase SDE (SDE2) to overall SDE and quantifies how much SDE changes with the addition of Phase II dispersal. Considering that the seed dispersal process is context dependent so that Phase II depends on Phase I, we predict the possible range of variation of SDE according to the variation of the probability of Phase II dispersal. In our specific study system composed of two primate species as primary dispersal vectors and different species of dung beetles as secondary dispersal vectors, the relative contribution of SDE1 and SDE2 to overall SDE varied between plant species. We discuss the context dependency of the Phase II dispersal and the potential applications of our approach. This extension to the conceptual framework of SDE enables quantitative evaluation of the effect of Phase II dispersal on plant fitness and can be easily adapted to other biotic and/or abiotic diplochorous dispersal systems.

CFD modeling of two-phase boiling flows in the slug flow regime with an interface capturing technique

The objective of this thesis was to improve the commercial CFD software Ansys Fluent to obtain a tool able to perform accurate simulations of flow boiling in the slug flow regime. The achievement of a reliable numerical framework allows a better understanding of the bubble and flow dynamics induced by the evaporation and makes possible the prediction of the wall heat transfer trends. In order to save computational time, the flow is modeled with an axisymmetrical formulation. Vapor and liquid phases are treated as incompressible and in laminar flow. By means of a single fluid approach, the flow equations are written as for a single phase flow, but discontinuities at the interface and interfacial effects need to be accounted for and discretized properly. Ansys Fluent provides a Volume Of Fluid technique to advect the interface and to map the discontinuous fluid properties throughout the flow domain. The interfacial effects are dominant in the boiling slug flow and the accuracy of their estimation is fundamental for the reliability of the solver. Self-implemented functions, developed ad-hoc, are introduced within the numerical code to compute the surface tension force and the rates of mass and energy exchange at the interface related to the evaporation. Several validation benchmarks assess the better performances of the improved software. Various adiabatic configurations are simulated in order to test the capability of the numerical framework in modeling actual flows and the comparison with experimental results is very positive. The simulation of a single evaporating bubble underlines the dominant effect on the global heat transfer rate of the local transient heat convection in the liquid after the bubble transit. The simulation of multiple evaporating bubbles flowing in sequence shows that their mutual influence can strongly enhance the heat transfer coefficient, up to twice the single phase flow value.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

Aqueous two-phase extraction using thermoseparating copolymer: a new system for phenolic compounds removal from hemicelullosic hydrolysate

BACKGROUND: The hydrolysis of hemicellulosic material can provide liquor with high xylose concentration (which can be used as a fermentation medium) and phenolic compounds (Phs), potentially immunostimulating compounds. However, these hydrolysates must be detoxified in order to remove the Phs that can act as inhibitors in bioconversions. RESULTS: Aqueous two-phase systems composed of thermoseparating copolymers were used for rice straw hydrolysate detoxification. The hydrolysis process was able to promote chemical breakdown of 85% of the total hemicellulose content, 14% of the cellulose, and 2% of the lignin. The hydrolysate obtained contained 19.7 g L-1 of xylose and several phenolic compounds, such as vanillin, vanillic acid, ferullic acid, etc. The phenolics extraction was studied as a function of copolymer molar mass (1100 g mol(-1), 2000 g mol(-1) and 2800 g mol(-1)), their percentages (from 5% to 50%) and Phs initial concentration. Phenolic compounds extraction of around 80% was obtained under the following conditions: 20% (w/w) and 35% (w/w) copolymer 1100 g mol-1, 35% (w/w) copolymer 2000 g mol(-1) and 35% (w/w) copolymer 2800 g mol(-1) at 25 degrees C. CONCLUSIONS: The results demonstrated the viability of this method for the removal of Phs from rice straw hydrolysate, which has potential uses in bioconversion processes. (c) 2007 Society of Chemical Industry.

Film thickness measurement techniques applied to micro-scale two-phase flow systems

Recently semi-empirical models to estimate flow boiling heat transfer coefficient, saturated CHF and pressure drop in micro-scale channels have been proposed. Most of the models were developed based on elongated bubbles and annular flows in the view of the fact that these flow patterns are predominant in smaller channels. In these models, the liquid film thickness plays an important role and such a fact emphasizes that the accurate measurement of the liquid film thickness is a key point to validate them. On the other hand, several techniques have been successfully applied to measure liquid film thicknesses during condensation and evaporation under macro-scale conditions. However, although this subject has been targeted by several leading laboratories around the world, it seems that there is no conclusive result describing a successful technique capable of measuring dynamic liquid film thickness during evaporation inside micro-scale round channels. This work presents a comprehensive literature review of the methods used to measure liquid film thickness in macro- and micro-scale systems. The methods are described and the main difficulties related to their use in micro-scale systems are identified. Based on this discussion, the most promising methods to measure dynamic liquid film thickness in micro-scale channels are identified. (C) 2009 Elsevier Inc. All rights reserved.

New prediction methods for CO2 evaporation inside tubes: Part I - A two-phase flow pattern map and a flow pattern based phenomenological model for two-phase flow frictional pressure drops

An updated flow pattern map was developed for CO2 on the basis of the previous Cheng-Ribatski-Wojtan-Thome CO2 flow pattern map [1,2] to extend the flow pattern map to a wider range of conditions. A new annular flow to dryout transition (A-D) and a new dryout to mist flow transition (D-M) were proposed here. In addition, a bubbly flow region which generally occurs at high mass velocities and low vapor qualities was added to the updated flow pattern map. The updated flow pattern map is applicable to a much wider range of conditions: tube diameters from 0.6 to 10 mm, mass velocities from 50 to 1500 kg/m(2) s, heat fluxes from 1.8 to 46 kW/m(2) and saturation temperatures from -28 to +25 degrees C (reduced pressures from 0.21 to 0.87). The updated flow pattern map was compared to independent experimental data of flow patterns for CO2 in the literature and it predicts the flow patterns well. Then, a database of CO2 two-phase flow pressure drop results from the literature was set up and the database was compared to the leading empirical pressure drop models: the correlations by Chisholm [3], Friedel [4], Gronnerud [5] and Muller-Steinhagen and Heck [6], a modified Chisholm correlation by Yoon et al. [7] and the flow pattern based model of Moreno Quiben and Thome [8-10]. None of these models was able to predict the CO2 pressure drop data well. Therefore, a new flow pattern based phenomenological model of two-phase flow frictional pressure drop for CO2 was developed by modifying the model of Moreno Quiben and Thome using the updated flow pattern map in this study and it predicts the CO2 pressure drop database quite well overall. (C) 2007 Elsevier Ltd. All rights reserved.

Numerical solution of the two-phase expansion of a metastable flashing liquid jet using the dispersion-controlled dissipative scheme

This paper presents a study of the stationary phenomenon of superheated or metastable liquid jets, flashing into a two-dimensional axisymmetric domain, while in the two-phase region. In general, the phenomenon starts off when a high-pressure, high-temperature liquid jet emerges from a small nozzle or orifice expanding into a low-pressure chamber, below its saturation pressure taken at the injection temperature. As the process evolves, crossing the saturation curve, one observes that the fluid remains in the liquid phase reaching a superheated condition. Then, the liquid undergoes an abrupt phase change by means of an oblique evaporation wave. Across this phase change the superheated liquid becomes a two-phase high-speed mixture in various directions, expanding to supersonic velocities. In order to reach the downstream pressure, the supersonic fluid continues to expand, crossing a complex bow shock wave. The balance equations that govern the phenomenon are mass conservation, momentum conservation, and energy conservation, plus an equation-of-state for the substance. A false-transient model is implemented using the shock capturing scheme: dispersion-controlled dissipative (DCD), which was used to calculate the flow conditions as the steady-state condition is reached. Numerical results with computational code DCD-2D vI have been analyzed. Copyright (C) 2009 John Wiley & Sons, Ltd.

Liquid-Liquid Equilibrium of Aqueous Two-Phase Systems Containing Some Synthetic Polyelectrolytes and Polyethylene Glycol

Experimental results are presented for the liquid-liquid equilibrium of aqueous two-phase systems containing a synthetic polyelectrolyte (polysodium acrylate, polysodium methacrylate, and polysodium ethylene sulfonate) and polyethylene glycol at (298.2 and 323.2) K. A total of 40 phase diagrams were obtained, comprising data both of the binodal curve (obtained through cloud-point measurements) and of equilibrium compositions. The influences of temperature, the nature of the polyelectrolyte monomer unit, and the chain length of both types of polymers are analyzed and discussed.

Extractive Fermentation of Clavulanic Acid by Streptomyces DAUFPE 3060 Using Aqueous Two-Phase System

The influence of four variables, specifically PEG molar mass (400, 1,000, and 8,000 g/mol), concentrations of PEG and phosphate salts (15, 20, and 25% for both), and agitation intensity (110, 150, and 200 rpm), on clavulanic acid (CA) extraction by extractive fermentation with PEG/phosphate salts aqueous two-phase system was investigated in shaken flasks using a 2(4-1)-fractional factorial design. After selection of the two most significant variables (agitation intensity and PEG molar mass), an optimization study conducted according to a 2(2)-central composite design revealed that 25% PEG 8,000 g/mol and phosphate salts at 240 rpm (run 6) were the best conditions for the extractive fermentation, leading to the best results in terms of partition coefficient (k = 8.2), yield of CA in the PEG-rich phase (eta(T) = 93%) and productivity (P = 5.3 mg/Lh). As a first attempt to make a scale-up of these results, the effectiveness of the extractive fermentation was then checked in a bench-scale bioreactor under conditions as close as possible to the optimum ones determined in flasks. The highest CA concentration obtained in the PEG-rich phase (691 mg/L) was 30% higher than in flasks, thus demonstrating the potential of such a new process, integrating the production and extraction steps, as a promising, low-cost tool to obtain high yields of this and similar products. (C) 2010 American Institute of Chemical Engineers Biotechnol. Prog., 27: 95-103, 2011