DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Determination of minimum fluidization velocity for gas-solid beds by experimental data and numerical simulations

The purpose of this work is to predict the minimum fluidization velocity Umf in a gas-solid fluidized bed. The study was carried out with an experimental apparatus for sand particles with diameters between 310μm and 590μm, and density of 2,590kg/m3. The experimental results were compared with numerical simulations developed in MFIX (Multiphase Flow with Interphase eXchange) open source code [1], for three different sizes of particles: 310mum, 450μm and 590μm. A homogeneous mixture with the three kinds of particles was also studied. The influence of the particle diameter was presented and discussed. The Ergun equation was also used to describe the minimum fluidization velocity. The experimental data presented a good agreement with Ergun equation and numerical simulations. Copyright © 2011 by ASME.

Simulação numérica de escoamentos gás-sólido em leito fluidizado borbulhante utilizando a teoria cinética dos escoamentos granulares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Efeitos numéricos na simulação de escoamentos gás-sólido em leito fuidizado borbulhante utilizando a teoria cinética dos escoamentos granulares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Report on the multiphase flow of gas, oil, and water through vertical flow strings with application to the design of gas lift installations,

Cover title.

Multiphase flow in the vascular system of wood: From microscopic exploration to 3-D Lattice Boltzmann experiments

This paper provides insights into liquid free water dynamics in wood vessels based on Lattice Boltzmann experiments. The anatomy of real wood samples was reconstructed from systematic 3-D analyses of the vessel contours derived from successive microscopic images. This virtual vascular system was then used to supply fluid-solid boundary conditions to a two-phase Lattice Boltzmann scheme and investigate capillary invasion of this hydrophilic porous medium. Behavior of the liquid phase was strongly dependent on anatomical features, especially vessel bifurcations and reconnections. Various parameters were examined in numerical experiments with ideal vessel bifurcations, to clarify our interpretation of these features. (c) 2010 Elsevier Ltd. All rights reserved.

Accurate and efficient simulation of multiphase flow in a heterogeneous reservoir by using error estimate and control in the multiscale finite-volume framework

The multiscale finite-volume (MSFV) method is designed to reduce the computational cost of elliptic and parabolic problems with highly heterogeneous anisotropic coefficients. The reduction is achieved by splitting the original global problem into a set of local problems (with approximate local boundary conditions) coupled by a coarse global problem. It has been shown recently that the numerical errors in MSFV results can be reduced systematically with an iterative procedure that provides a conservative velocity field after any iteration step. The iterative MSFV (i-MSFV) method can be obtained with an improved (smoothed) multiscale solution to enhance the localization conditions, with a Krylov subspace method [e.g., the generalized-minimal-residual (GMRES) algorithm] preconditioned by the MSFV system, or with a combination of both. In a multiphase-flow system, a balance between accuracy and computational efficiency should be achieved by finding a minimum number of i-MSFV iterations (on pressure), which is necessary to achieve the desired accuracy in the saturation solution. In this work, we extend the i-MSFV method to sequential implicit simulation of time-dependent problems. To control the error of the coupled saturation/pressure system, we analyze the transport error caused by an approximate velocity field. We then propose an error-control strategy on the basis of the residual of the pressure equation. At the beginning of simulation, the pressure solution is iterated until a specified accuracy is achieved. To minimize the number of iterations in a multiphase-flow problem, the solution at the previous timestep is used to improve the localization assumption at the current timestep. Additional iterations are used only when the residual becomes larger than a specified threshold value. Numerical results show that only a few iterations on average are necessary to improve the MSFV results significantly, even for very challenging problems. Therefore, the proposed adaptive strategy yields efficient and accurate simulation of multiphase flow in heterogeneous porous media.

A high frequency electronic transducer for multiphase flow measurements

This paper describes an electronic transducer for multiphase flow measurement. Its high sensitivity, good signal to noise ratio and accuracy are achieved through an electrical impedance sensor with a special guard technique. The transducer consists of a wide bandwidth and high slew rate differentiator where the lead inductance and stray capacitance effects are compensated. The sensor edge effect is eliminated by using a guard electrode based on the virtual ground potential of the operational amplifier. A theoretical modeling and a calibration method are also presented. The results obtained seem to confirm the validity of the proposed technique.

Influence of impact angles on penetration rates of CRAs exposed to a high velocity multiphase flow

A combined flow loop - jet impingement pilot plant has been used to determine mass loss rates in a mixed gas - saltwater - sand multiphase flow at impact velocities up to 70 m/s. Artificial brine with a salt content of 27 g/1 was used as liquid phase. Sand content, with grain size below 150 µ, was 2.7 g/l brine. CO at a pressure of 15 bar was used as gas phase. The impact angle between jet stream (nozzle) and sample surface was varied between 30 and 90°. Rectangular stainless steel disc samples with a size of 20 × 15 × 5 mm were used. They were mechanically ground and polished prior to testing. Damaged surfaces of specimens exposed to the high velocity multiphase flow were investigated by stereo microscopy, scanning electron microscopy (SEM) and an optical device for 3D surface measurements. Furthermore, samples were investigated by applying atomic force microscopy (AFM), magnetic force microscopy (MFM) and nanoindentation. Influence of impact velocity and impact angle on penetration rates (mass loss rates) of two CRAs (UNS S30400 and N08028) are presented. Moreover effects of chemical composition and mechanical properties are critically discussed. © 2008 by NACE International.

Volume averaging of multiphase flows with hydrate formation in subsea pipelines

In oil and gas pipeline operations, the gas, oil, and water phases simultaneously move through pipe systems. The mixture cools as it flows through subsea pipelines, and forms a hydrate formation region, where the hydrate crystals start to grow and may eventually block the pipeline. The potential of pipe blockage due to hydrate formation is one of the most significant flow-assurance problems in deep-water subsea operations. Due to the catastrophic safety and economic implications of hydrate blockage, it is important to accurately predict the simultaneous flow of gas, water, and hydrate particles in flowlines. Currently, there are few or no studies that account for the simultaneous effects of hydrate growth and heat transfer on flow characteristics within pipelines. This thesis presents new and more accurate predictive models of multiphase flows in undersea pipelines to describe the simultaneous flow of gas, water, and hydrate particles through a pipeline. A growth rate model for the hydrate phase is presented and then used in the development of a new three-phase model. The conservation equations of mass, momentum, and energy are formulated to describe the physical phenomena of momentum and heat transfer between the fluid and the wall. The governing equations are solved based on an analytical-numerical approach using a Newton-Raphson method for the nonlinear equations. An algorithm was developed in Matlab software to solve the equations from the inlet to the outlet of the pipeline. The developed models are validated against a single-phase model with mixture properties, and the results of comparative studies show close agreement. The new model predicts the volume fraction and velocity of each phase, as well as the mixture pressure and temperature profiles along the length of the pipeline. The results from the hydrate growth model reveal the growth rate and location where the initial hydrates start to form. Finally, to assess the impact of certain parameters on the flow characteristics, parametric studies have been conducted. The results show the effect of a variation in the pipe diameter, mass flow rate, inlet pressure, and inlet temperature on the flow characteristics and hydrate growth rates.

Multiphase flow in porous media: meta-modeling techniques for sensitivity analysis and risk assessment

One of the biggest challenges that contaminant hydrogeology is facing, is how to adequately address the uncertainty associated with model predictions. Uncertainty arise from multiple sources, such as: interpretative error, calibration accuracy, parameter sensitivity and variability. This critical issue needs to be properly addressed in order to support environmental decision-making processes. In this study, we perform Global Sensitivity Analysis (GSA) on a contaminant transport model for the assessment of hydrocarbon concentration in groundwater. We provide a quantification of the environmental impact and, given the incomplete knowledge of hydrogeological parameters, we evaluate which are the most influential, requiring greater accuracy in the calibration process. Parameters are treated as random variables and a variance-based GSA is performed in a optimized numerical Monte Carlo framework. The Sobol indices are adopted as sensitivity measures and they are computed by employing meta-models to characterize the migration process, while reducing the computational cost of the analysis. The proposed methodology allows us to: extend the number of Monte Carlo iterations, identify the influence of uncertain parameters and lead to considerable saving computational time obtaining an acceptable accuracy.

DEVELOPMENT OF STATIC MIXERS FOR MISCIBLE FLUIDS IN LAMINAR FLOW WITH THE USE OF COMPUTATIONAL FLUID DYNAMICS (CFD)

Static mixers with improved performance were developed from CFD simulations in a stepwise approach. The relevant geometric features of simple mixer designs and the corresponding mixing mechanisms-laminar shear, elongational flow, and distributive mixing-were identified first. This information was used to formulate guidelines for the development of new geometries. The solid elements of the static mixer should: (a) provide restrictions to the flow; (b) deflect the flow; (c) be sequentially rotated around the flow direction to provide symmetry; (d) extend from the center of the pipe to the vicinity of the walls to avoid short-circuiting; and (e) distribute and remix the flow. Based on these guidelines, two improved mixer designs were developed: the DS A-I mixer has a good mixing efficiency and an acceptable pressure drop; the Fins 35 degrees mixer is more efficient and compact, but requires a larger pressure drop. Their performance indicates that their use is possible on industrial applications.

Measurement of colour flow with the jet pull angle in t(t)over-bar events using the ATLAS detector at root s=8 TeV

The distribution and orientation of energy inside jets is predicted to be an experimental handle on colour connections between the hard--scatter quarks and gluons initiating the jets. This Letter presents a measurement of the distribution of one such variable, the jet pull angle. The pull angle is measured for jets produced in tt¯ events with one W boson decaying leptonically and the other decaying to jets using 20.3 fb−1 of data recorded with the ATLAS detector at a centre--of--mass energy of s√=8 TeV at the LHC. The jet pull angle distribution is corrected for detector resolution and acceptance effects and is compared to various models.

Error estimate and control in the MSFV method for multiphase flow in porous media

n this paper the iterative MSFV method is extended to include the sequential implicit simulation of time dependent problems involving the solution of a system of pressure-saturation equations. To control numerical errors in simulation results, an error estimate, based on the residual of the MSFV approximate pressure field, is introduced. In the initial time steps in simulation iterations are employed until a specified accuracy in pressure is achieved. This initial solution is then used to improve the localization assumption at later time steps. Additional iterations in pressure solution are employed only when the pressure residual becomes larger than a specified threshold value. Efficiency of the strategy and the error control criteria are numerically investigated. This paper also shows that it is possible to derive an a-priori estimate and control based on the allowed pressure-equation residual to guarantee the desired accuracy in saturation calculation.

Quantification of myocardial blood flow with 82Rb positron emission tomography: clinical validation with 15O-water.

PURPOSE: Quantification of myocardial blood flow (MBF) with generator-produced (82)Rb is an attractive alternative for centres without an on-site cyclotron. Our aim was to validate (82)Rb-measured MBF in relation to that measured using (15)O-water, as a tracer 100% of which can be extracted from the circulation even at high flow rates, in healthy control subject and patients with mild coronary artery disease (CAD). METHODS: MBF was measured at rest and during adenosine-induced hyperaemia with (82)Rb and (15)O-water PET in 33 participants (22 control subjects, aged 30 ± 13 years; 11 CAD patients without transmural infarction, aged 60 ± 13 years). A one-tissue compartment (82)Rb model with ventricular spillover correction was used. The (82)Rb flow-dependent extraction rate was derived from (15)O-water measurements in a subset of 11 control subjects. Myocardial flow reserve (MFR) was defined as the hyperaemic/rest MBF. Pearson's correlation r, Bland-Altman 95% limits of agreement (LoA), and Lin's concordance correlation ρ (c) (measuring both precision and accuracy) were used. RESULTS: Over the entire MBF range (0.66-4.7 ml/min/g), concordance was excellent for MBF (r = 0.90, [(82)Rb-(15)O-water] mean difference ± SD = 0.04 ± 0.66 ml/min/g, LoA = -1.26 to 1.33 ml/min/g, ρ(c) = 0.88) and MFR (range 1.79-5.81, r = 0.83, mean difference = 0.14 ± 0.58, LoA = -0.99 to 1.28, ρ(c) = 0.82). Hyperaemic MBF was reduced in CAD patients compared with the subset of 11 control subjects (2.53 ± 0.74 vs. 3.62 ± 0.68 ml/min/g, p = 0.002, for (15)O-water; 2.53 ± 1.01 vs. 3.82 ± 1.21 ml/min/g, p = 0.013, for (82)Rb) and this was paralleled by a lower MFR (2.65 ± 0.62 vs. 3.79 ± 0.98, p = 0.004, for (15)O-water; 2.85 ± 0.91 vs. 3.88 ± 0.91, p = 0.012, for (82)Rb). Myocardial perfusion was homogeneous in 1,114 of 1,122 segments (99.3%) and there were no differences in MBF among the coronary artery territories (p > 0.31). CONCLUSION: Quantification of MBF with (82)Rb with a newly derived correction for the nonlinear extraction function was validated against MBF measured using (15)O-water in control subjects and patients with mild CAD, where it was found to be accurate at high flow rates. (82)Rb-derived MBF estimates seem robust for clinical research, advancing a step further towards its implementation in clinical routine.