Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Device-scale CFD modeling of gas-liquid multiphase flow and amine absorption for CO2 capture

In this paper we aim to develop a validated device-scale CFD model that can predict quantitatively both hydrodynamics and CO2 capture efficiency for an amine-based solvent absorber column with random Pall ring packing. A Eulerian porous-media approach and a two-fluid model were employed, in which the momentum and mass transfer equations were closed by literature-based empirical closure models. We proposed a hierarchical approach for calibrating the parameters in the closure models to make them accurate for the packed column. Specifically, a parameter for momentum transfer in the closure was first calibrated based on data from a single experiment. With this calibrated parameter, a parameter in the closure for mass transfer was next calibrated under a single operating condition. Last, the closure of the wetting area was calibrated for each gas velocity at three different liquid flow rates. For each calibration, cross validations were pursued using the experimental data under operating conditions different from those used for calibrations. This hierarchical approach can be generally applied to develop validated device-scale CFD models for different absorption columns.

Pan, Wenxiao↗
DOE OSTI2016

Simulation of Subsurface Multiphase Contaminant Extraction Using a Bioslurping Well Model

Subsurface simulation of multiphase extraction from wells is notoriously difficult. Explicit representation of well geometry requires small grid resolution, potentially leading to large computational demands. To reduce the problem dimensionality, multiphase extraction is mostly modeled using vertically-averaged approaches. In this paper, a multiphase well model approach is presented as an alternative to simplify the application. The well model, a multiphase extension of the classic Peaceman model, has been implemented in the STOMP simulator. The numerical solution approach accounts for local conditions and gradients in the exchange of fluids between the well and the aquifer. Advantages of this well model implementation include the option to simulate the effects of well characteristics and operation. Simulations were conducted investigating the effects of extraction location, applied vacuum pressure, and a number of hydraulic properties. The obtained results were all consistent and logical. A major outcome of the test simulations is that, in contrast with common recommendations to extract from either the gas-NAPL or the NAPL-aqueous phase interface, the optimum extraction location should be in between these two levels. Further, the new model implementation was also used to simulate extraction at a field site in Brazil. The simulation shows a good match with the field data, suggesting that the new STOMP well module may correctly represent oil removal. The field simulations depend on the quality of the site conceptual model, including the porous media and contaminant properties and the boundary and extraction conditions adopted. The new module may potentially be used to design field applications and analyze extraction data.

42 ENGINEERING↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2019

A scalable Euler–Lagrange approach for multiphase flow simulation on spectral elements

Multiphase flow can be difficult to simulate with high accuracy due to the wide range of scales associated with various multiphase phenomena. These scales may range from the size of individual particles to the entire domain of interest. Traditionally, large scale systems can only be simulated using averaging approaches that filter out the locations of individual particles. In this work, the Euler–Lagrange method is used to simulate large-scale dense particle systems in which each individual particle is tracked. In order to accomplish this, the highly scalable spectral element code nek5000 has been extended to handle the multiple levels of multiphase coupling in these systems. These levels include what has been called one-, two-, and four-way coupling. Here, each level has been separated to detail the computational impact of each stage. A binned ghost particle algorithm has also been developed to efficiently handle the challenges of two- and four-way coupling in a parallel processing context. The algorithms and their implementations are then shown to scale to 65,536 Message Passing Interface (MPI) ranks in both the strong and weak limits. After this, validation is performed through simulation of a small-scale fluidized bed. Lastly, a large-scale fluidized bed is simulated with 65,536 MPI ranks and is able to capture the unique physics of the onset of fluidization.

Computer Science↗
DOE OSTI2020

Geochemical narrowing of cement fracture aperture during multiphase flow of supercritical CO2 and brine

For carbon capture and storage operations, wells are used as a necessary conduit for injecting CO2 at depth, but they can also act as leakage conduits if the cement seals are compromised. The specific objective of this research was to investigate the nature of self-healing of a fracture within cement under multiphase flow of CO2 and brine, and to compare this to a newly developed model. This was accomplished through a laboratory investigation, flowing reacting supercritical CO2 and brine through a cement fracture with multiphase flow. The influent end of the fracture was degraded as evidenced by the formation of cracks across the surface. At the effluent end of the fracture, the initial aperture of 137 µm was reduced to 50 µm. This reduction by 87 µm was similar to results from a recently developed model tested in this study, which predicted an aperture reduction of 80 µm. Our results demonstrated that the model accurately predicted the reactions which occurred in this multiphase flow system.

wellbore, cement, fracture, multiphase-flow, super↗
DOE OSTI2017

Modeling variability in porescale multiphase flow experiments

Microfluidic devices and porescale numerical models are commonly used to study multiphase flow in biological, geological, and engineered porous materials. In this work, we perform a set of drainage and imbibition experiments in six identical microfluidic cells to study the reproducibility of multiphase flow experiments. We observe significant variations in the experimental results, which are smaller during the drainage stage and larger during the imbibition stage. We demonstrate that these variations are due to sub-porescale geometry differences in microcells (because of manufacturing defects) and variations in the boundary condition (i.e.,fluctuations in the injection rate inherent to syringe pumps). Computational simulations are conducted using commercial software STAR-CCM+, both with constant and randomly varying injection rate. Stochastic simulations are able to capture variability in the experiments associated with the varying pump injection rate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Liquid-based gating mechanism with tunable multiphase selectivity and antifouling behaviour

Living organisms make extensive use of micro- and nanometre-sized pores as gatekeepers for controlling the movement of fluids, vapours and solids between complex environments. The ability of such pores to coordinate multiphase transport, in a highly selective and subtly triggered fashion and without clogging, has inspired interest in synthetic gated pores for applications ranging from fluid processing to 3D printing and lab-on-chip systems(1-10). But although specific gating and transport behaviours have been realized by precisely tailoring pore surface chemistries and pore geometries(6,11-17), a single system capable of controlling complex, selective multiphase transport has remained a distant prospect, and fouling is nearly inevitable(11,12). Here we introduce a gating mechanism that uses a capillary-stabilized liquid as a reversible, reconfigurable gate that fills and seals pores in the closed state, and creates a non-fouling, liquid-lined pore in the open state. Theoretical modelling and experiments demonstrate that for each transport substance, the gating threshold-the pressure needed to open the pores-can be rationally tuned over a wide pressure range. This enables us to realize in one system differential response profiles for a variety of liquids and gases, even letting liquids flow through the pore while preventing gas from escaping. These capabilities allow us to dynamically modulate gas-liquid sorting in a microfluidic flow and to separate a three-phase air-water-oil mixture, with the liquid lining ensuring sustained antifouling behaviour. Because the liquid gating strategy enables efficient long-term operation and can be applied to a variety of pore structures and membrane materials, and to micro- as well as macroscale fluid systems, we expect it to prove useful in a wide range of applications.

Hou, X↗
DOE OSTI2018

A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state: A localized artificial diffusivity method to simulate compressible multiphase flows using the stiffened gas equation of state

The development of numerical approaches to perform direct numerical simulations of compressible multiphase flows has been an active field of research for several years. Proper treatment of fluid interfaces is crucial as important physics occur in this infinitesimally small region. Furthermore, the compressibility of the fluid requires proper treatment of discontinuities. Artificial diffusivity is among a number of methods widely used for compressible flows. This study develops a general form of consistent artificial diffusion fluxes and extends the localized artificial diffusivity method for high-order central schemes to solve multiphase flows with an interface-capturing method. These fluxes ensure an oscillation-free interface for pressure, velocity, and temperature without employing a sharpening technique. Moreover, the high-order representation of all scales in the flow helps capture the wide range of instabilities inherent in these flows. The goal is to develop an approach capable of performing high-fidelity simulations supported by physics-driven validation. lastly, this is achieved by solving the five-equation model with the stiffened-gas equation of state using the proposed method for multicomponent and multiphase flows on a variety of 1D and 2D problems.

42 ENGINEERING↗
DOE OSTI2019

Comprehensive study of the interactions between the critical dimensionless numbers associated with multiphase flow in 3D porous media

Multiphase flow in porous media is of signifcant interest in many engineering applications, such as geologic carbon sequestration, enhanced oil recovery, and groundwater contamination and remediation. In order to advance the fundamental understanding of multiphase flow in complex three-dimensional (3D) porous media, the interactions between the critical dimensionless numbers, including the contact angle, viscosity ratio, and capillary (Ca) number, were explored using X-ray micro-computed tomography (micro-CT) scanning and lattice Boltzmann (LB) modeling. In this study, the 3D pore structure information was extracted from micro-CT images and then used as interior boundary conditions of flow modeling in a pore-scale LB simulator to simulate multiphase flow within the pore space. A Berea sandstone sample was scanned and then two-phase flow LB simulations were performed based on the micro-CT images. The LB-simulated water/CO 2 distributions agreed well with the micro-CT scanned images. Simulation findings showed that a decreasing contact angle causes a decrease in wetting-fluid relative permeability and an increase in non-wetting fluid relative permeability. A rising Ca number increases both wetting and non-wetting fluid relative permeabilities. An increasing viscosity ratio (the ratio of non-wetting fluid viscosity to wetting fluid viscosity) facilitates the increase of non-wetting fluid relative permeability and mitigates the reduction of wetting fluid relative permeability, when the contact angle decreases continuously. The primary novel finding of this study is that the viscosity ratio affects the rate of change of the relative permeability curves for both phases when the contact angle changes continuously. To date, we believe it is the first time that comprehensive interactions between these dimensionless numbers are demonstrated in a sandstone sample based on real 3D structures. We also investigated the role of the changes of flow direction and sample location on relative permeability curves. Simulation results showed that the change in non-wetting fluid relative permeability was larger when the flow direction was switched from vertical to horizontal, which indicated that there was stronger anisotropy in larger pore networks that were primarily occupied by the non-wetting fluid.

42 ENGINEERING↗
DOE OSTI2016

Smoothed dissipative particle dynamics model for mesoscopic multiphase flows in the presence of thermal fluctuations

Thermal fluctuations cause perturbations of fluid-fluid interfaces and highly nonlinear hydrodynamics in multiphase flows. In this work, we develop a novel multiphase smoothed dissipative particle dynamics model. This model accounts for both bulk hydrodynamics and interfacial fluctuations. Interfacial surface tension is modeled by imposing a pairwise force between SDPD particles. We show that the relationship between the model parameters and surface tension, previously derived under the assumption of zero thermal fluctuation, is accurate for fluid systems at low temperature but overestimates the surface tension for intermediate and large thermal fluctuations. To analyze the effect of thermal fluctuations on surface tension, we construct a coarse-grained Euler lattice model based on the mean field theory and derive a semi-analytical formula to directly relate the surface tension to model parameters for a wide range of temperatures and model resolutions. We demonstrate that the present method correctly models the dynamic processes, such as bubble coalescence and capillary spectra across the interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Enhancing Oxygen Reduction Activity and Cr Tolerance of Solid Oxide Fuel Cell Cathodes by a Multiphase Catalyst Coating

Intermediate temperature solid oxide fuel cells (IT-SOFCs) are cost-effective and efficient energy conversion systems. Here, the sluggish oxygen reduction reaction (ORR) and the degradation of cathodes are critical challenges to the commercialization of IT-SOFCs. Here, a highly efficient multiphase (MP) catalyst coating, consisting of Ba 1–x Co 0.7 Fe 0.2 Nb 0.1 O 3–δ (BCFN) and BaCO 3 , to enhance the ORR activity and durability of the state-of-the-art lanthanum strontium cobalt ferrite (La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ , LSCF) cathode is reported. The conformal MP catalyst-coated LSCF cathode shows a polarization resistance (R p ) of 0.048 Ω cm 2 at 650 °C, about one order of magnitude smaller than that of the bare LSCF. In an accelerated Cr-poisoning test, the degradation rate of the catalyst-coated LSCF electrode is 10 –3 Ω cm 2 h –1 (0.59% h –1 ) over 200 h, only one fifth of the degradation rate of the bare LSCF electrode at 750 °C. In addition, anode-supported single cells with the MP catalyst-coated LSCF cathode show a dramatically enhanced peak power density (1.4 W cm –2 vs 0.67 W cm –2 at 750 °C) and increased durability against Cr and H 2 O. Both experimental results and density functional theory-based calculations indicate that the BCFN phase improves the ORR activity while the BaCO 3 phase enhances the stability of the LSCF cathode.

solid oxide fuel cells↗
DOE OSTI2013

Multiphase Resistivity Model for Magnetic Nanocomposites Developed for High Frequency, High Power Transformation

New power conversion systems that offer promise to transform electricity grids into unified interactive supply networks require high-resistivity soft-magnetic materials to allow for switching of magnetic materials at frequencies approaching 100 kHz for power transformation in the megawatt range. Amorphous and nanocomposite soft-magnetic materials, which represent the state of the art in terms of high power densities and low losses at high frequencies, have resistivities that depend on the structures and spatial distributions of multiple phases in thin ribbons. We present a multiphase resistivity model applicable to nanocomposite materials by considering an equivalent circuit approach considering paths through an amorphous, crystalline, and growth inhibitor shell phase. We detail: (a) identification of amorphous, crystalline, and shell phases; (b) consideration of the role of the morphology of each phase in an equivalent circuit model for the resistance; (c) a two-band model for the Fe/Co composition dependence of the resistivity in crystalline and amorphous phases; (d) a virtual bound state model for resistivity to explain increased resistivity due to early transition-metal growth inhibitors in the shell surrounding the nanocrystalline phase; and (e) disorder effects on amorphous phase resistivity. Experimental design and results for systems of interest in high-frequency power transformation are discussed in the context of our model including: (a) techniques for measurements of cross-section and density, (b) four-point probe and surface resistivity measurements, and (c) measurements in Fe- and Co-rich systems comparing amorphous and nanocomposite materials.

DeGeorge, V↗
DOE OSTI2017

Cesium incorporation in hollandite-rich multiphasic ceramic waste forms

Hollandite-rich multiphase waste form compositions processed by melt-solidification and spark plasma sintering (SPS) were characterized, compared, and validated for nuclear waste incorporation. Phase identification by x-ray diffraction (XRD) and electron back-scattered diffraction (EBSD) confirmed hollandite as the major phase present in these samples along with perovskite, pyrochlore and zirconolite. Distribution of selected elements observed by wavelength dispersive spectroscopy (WDS) maps indicated that Cs formed a secondary phase during SPS processing, which was considered undesirable. On the other hand, Cs partitioned into the hollandite phase in melt-processed samples. Further analysis of hollandite structure in melt-processed composition by selected area electron diffraction (SAED) revealed ordered arrangement of tunnel ions (Ba/Cs) and vacancies, suggesting efficient Cs incorporation into the lattice.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Trace Uranium Partitioning in a Multiphase Nano-FeOOH System

The characterization of trace elements in nanomaterials using extended X-ray absorption fine structure (EXAFS) spectroscopy constitutes a first step toward understanding how impurities or dopants affect the properties of the host phase. However, limitations to EXAFS interpretation complicate the analysis of trace concentrations of impurities that are distributed across multiple phases in a heterogeneous system. Ab initio molecular dynamics (AIMD)-informed EXAFS analysis was employed to investigate the immobilization of trace uranium associated with nanophase iron (oxyhydr)oxides, a model system for the geochemical sequestration of radiotoxic contaminants. The reductive transformation of ferrihydrite (Fe(OH)3) to nano-particulate iron oxyhydroxide minerals in the presence of uranyl (UO2)2+(aq) resulted in the preferential incorporation of U into goethite (a-FeOOH) over lepidocrocite (g-FeOOH), even though reaction conditions favored the formation of excess lepidocrocite. This unexpected result is supported by atomically resolved transmission electron microscopy. Using this model system, we demonstrate how AIMD-informed EXAFS analysis lifts the strict statistical limitations of traditional shell-by-shell EXAFS modeling, enabling the detailed analysis of the local bonding environment, charge compensation mechanisms, and oxidation states of polyvalent impurities in complex multi-phase nano-systems.

ab initio molecular dynamics↗