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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Automated contact angle estimation for three-dimensional X-ray microtomography data

Multiphase flow in capillary regimes is a fundamental process in a number of geoscience applications. The ability to accurately define wetting characteristics of porous media can have a large impact on numerical models. In this paper, a newly developed automated three-dimensional contact angle algorithm is described and applied to high-resolution X-ray microtomography data from multiphase bead pack experiments with varying wettability characteristics. The algorithm calculates the contact angle by finding the angle between planes fit to each solid/fluid and fluid/fluid interface in the region surrounding each solid/fluid/fluid contact point. Results show that the algorithm is able to reliably compute contact angles using the experimental data. The in situ contact angles are typically larger than flat surface laboratory measurements using the same material. Furthermore, wetting characteristics in mixed-wet systems also change significantly after displacement cycles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

An extension of Darcy’s law incorporating dynamic length scales

We propose a physics-based, macroscale formulation of multiphase porous-media flows that both honors the validity of Darcy’s law in steady or near-steady flows and accommodates the effects of heterogeneities and nonlinearities in unsteady flows. The new formulation recognizes that parameters characterizing the system operate at different length-scales. In particular, the use of Darcy’s law, predicated on the assumption of near-steady flows, requires dynamic length scales, owing to the possibility of rapid fluctuations in rock properties and fluid saturations attributable to heterogeneity and nonlinearity. We accommodate these dynamic length-scales through dynamic spatial averaging. The length and position of the averaging window are characterized by the length of the mixing zone and the direction of propagation of information in the transport process. Here, we validate the proposed formulation by comparing highly accurate,two dimensional numerical solutions against core-scale displacement experiments. The proposed paradigm is consistent with the classical multiphase Darcy formulation, in the sense that the latter also represents the results of an averaging approach.

42 ENGINEERING↗
DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2022

Transient simulation of biomass combustion in a circulating fluidized bed riser

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

09 BIOMASS FUELS↗
DOE OSTI2018

Burial-driven methane recycling in marine gas hydrate systems

Natural gas hydrate may be buried with sediments until it is no longer stable at a given pressure and temperature, resulting in conversion of hydrate into free gas. This gas may migrate upward and recycle back into the hydrate stability zone to form hydrate. As of yet, however, no quantitative description of the methane recycling process has been developed using multiphase flow simulations to model burial-driven gas hydrate recycling. In this study, we present a series of 1D multiphase transport simulations to investigate the methane recycling process in detail. By invoking the effects of capillary phenomena on hydrate and gas formation in pores of varying size, we find that a free gas phase can migrate a significant distance above the bulk base of hydrate stability. Since the top of the free gas occurrence is often identified as the base of the hydrate stability zone from seismic data, our results demonstrate that not only could this assumption mischaracterize a hydrate system, but that under recycling conditions the highest hydrate saturations can occur beneath the top of the free gas occurrence. We show that the presence of pore size distributions requires a replacement zone through which hydrate saturations progressively decrease with depth and are replaced with free gas. This replacement zone works to buffer against significant gas buildup that could lead to fracturing of overlying sediments. In conclusion, this work provides a framework for simulating flow and transport of methane within the 3-phase stability zone from a mass conservation perspective.

58 GEOSCIENCES↗
DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows

A mass-conserving lattice Boltzmann method (LBM) for multiphase flows is presented in this paper. The proposed LBM improves a previous model (Lee and Liu, 2010 [21]) in terms of mass conservation, speed-up, and efficiency, and also extends its capabilities for implementation on non-uniform grids. The presented model consists of a phase-field lattice Boltzmann equation (LBE) for tracking the interface between different fluids and a pressure-evolution LBM for recovering the hydrodynamic properties. In addition to the mass conservation property and the simplicity of the algorithm, the advantages of the current phase-field LBE are that it is an order of magnitude faster than the previous interface tracking LBE proposed by Lee and Liu (2010) [21] and it requires less memory resources for data storage. Meanwhile, the pressure-evolution LBM is equipped with a multi-relaxation-time (MRT) collision operator to facilitate attainability of small relaxation rates thereby allowing simulation of multiphase flows at higher Reynolds numbers. Additionally, we reformulate the presented MRT-LBM on nonuniform grids within an adaptive mesh refinement (AMR) framework. Various benchmark studies such as a rising bubble and a falling drop under buoyancy, droplet splashing on a wet surface, and droplet coalescence onto a fluid interface are conducted to examine the accuracy and versatility of the proposed AMR-LBM. The proposed model is further validated by comparing the results with other LB models on uniform grids. A factor of about 20 in savings of computational resources is achieved by using the proposed AMR-LBM. As a more demanding application, the Kelvin–Helmholtz instability (KHI) of a shear-layer flow is investigated for both density-matched and density-stratified binary fluids. The KHI results of the density-matched fluids are shown to be in good agreement with the benchmark AMR results based on the sharp-interface approach. When a density contrast between the two fluids exists, a typical chaotic structure in the flow field is observed at a Reynolds number of 10000, which indicates that the proposed model is a promising tool for direct numerical simulation of two-phase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Wettability Effects on Primary Drainage Mechanisms and NAPL Distribution: A Pore-Scale Study

The pore-scale processes governing water drainage behavior in porous media have implications for geoscience multiphase scenarios including carbon capture and storage, contaminant site remediation, oil recovery, and vadose zone processes. However, few studies report directly observed pore-scale water drainage phenomena in 3-D soils. This knowledge gap limits our ability to verify assumptions underlying existing models and develop optimal solutions. As such, this paper utilizes synchrotron X-ray microtomography to present an experimental pore-scale examination of nonaqueous phase liquid (NAPL)/water distribution along a primary drainage front as dense NAPL was injected upward into water wetting (WW) and intermediate wetting (IW) sand-packed columns. Pore-network structures were extracted from imaged data sets and mapped onto segmented NAPL/water data sets which allowed quantitative examinations of wettability impacts on (a) the extent to which NAPL fills individual pore bodies and (b) relationships between pore size and the phase occupying the pore, with both considered as a function of distance (and capillary pressure) relative to the NAPL front. These results revealed that several hypotheses treating IW sand similarly to WW sands are simplistic. IW systems exhibited a sequence of pore filling that deviated from traditional capillary pressure-based model predictions: NAPL invades smaller pores, while larger, adjacent pores are bypassed leaving multipore residual water ganglia. NAPL pore saturations were close to 1 and did not change with capillary pressure in IW systems. Overall, the results illustrate how a relatively small change in operative contact angle alters NAPL distribution during water drainage, with important implications for geoscience multiphase flow scenarios.

58 GEOSCIENCES↗
DOE OSTI2022

Preferential Solute Transport in Low Permeability Zones During Spontaneous Imbibition in Heterogeneous Porous Media

Multiphase flows in porous media|and the associated solute transport processes|are controlled by a combination of gravity, capillary, and viscous forces. Geologic heterogeneity influences flow and transport processes, yet gaps exist in our understanding of how heterogeneity impacts capillary-driven transport, such as during spontaneous imbibition. Here we use positron emission tomography, combined with a newly developed method for conducting spontaneous imbibition experiments, to observe solute transport into low permeability regions during imbibition. Using an experimentally parameterized 2D numerical model, we demonstrate how capillary-driven ow near the imbibition front can carry solutes into low permeability regions. This process displaces solutes from high permeability zones while the cumulative amount of solute in low permeability zones increases- opposite of what is observed under fully saturated solute transport conditions. These results highlight the complex flow and solute transport behavior that arise during multiphase displacements in heterogeneous geologic porous media.

58 GEOSCIENCES↗
DOE OSTI2022

Potential Seismicity Along Basement Faults Induced by Geological Carbon Sequestration

Large-scale CO 2 sequestration into geological formations has been suggested to reduce CO 2 emissions from industrial activities. However, much like enhanced geothermal stimulation and wastewater injection, CO 2 sequestration has a potential to induce earthquake along weak faults, which can be considered a negative impact on safety and public opinion. This research shows the physical mechanisms of potential seismic hazards along basement faults driven by CO 2 sequestration under variation in geological and operational constraints. Specifically we compare the poroelastic behaviors between multiphase flow and single-phase flow cases, highlighting specific needs of evaluating induced seismicity associated with CO 2 sequestration. In contrast to single-phase injection scenario, slower migration of the CO 2 plume than pressure pulse may delay accumulation of pressure and stress along basement faults that may not be mitigated immediately by shut-in of injection. The impact of multiphase flow system, therefore, needs to be considered for proper monitoring and mitigation strategies.

58 GEOSCIENCES↗
DOE OSTI2014

A comparative study of lock-in-amplifiers and improved duffing chaotic oscillators for the detection and processing of weak signals

The measurement of multiphase flow parameters is of great importance in a wide range of industries. In the measurement of multiphase, the signals from the sensors are extremely weak and often buried in strong background noise. It is thus desirable to develop effective signal processing techniques that can detect the weak signal from the sensor outputs. In this paper, two methods, i.e., lock-in-amplifier (LIA) and improved Duffing chaotic oscillator are compared to detect and process the weak signal. For sinusoidal signal buried in noise, the correlation detection with sinusoidal reference signal is simulated by using LIA. The improved Duffing chaotic oscillator method, which based on the Wigner transformation, can restore the signal waveform and detect the frequency. Two methods are combined to detect and extract the weak signal. Simulation results show the effectiveness and accuracy of the proposed improved method. The comparative analysis shows that the improved Duffing chaotic oscillator method can restrain noise strongly since it is sensitive to initial conditions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Integrated rheology model: Explosive Composition B-3

Composition B-3 (Comp B-3) is a high explosive formulation composed of 60/40wt% RDX (1,3,5-trinitroperhydro-1,3,5-triazine) /TNT (2,4,6 trinitrotoluene). Above approximately 78°C this formulation partially melts to form a multiphase system with solid RDX particles in a molten TNT matrix. This multiphase system presents a number of phenomena that influence its apparent viscosity. In an earlier study explosive Composition B-3 (Comp B-3, 60/40wt% RDX/TNT) was examined for evidence of yield stress using a non-isothermal falling ball viscometer and a yield stress model was proposed in this paper. An integrated viscosity model suitable for use in computational fluid dynamics (CFD) simulations is developed to capture the transition from a heterogeneous solid to a Bingham viscoplastic fluid. This viscosity model is used to simulate the motion of imbedded spheres falling through molten Comp B-3. Finally, comparison of the simulations to physical tests show agreement between the positions predicted by the model and the measured locations of the spheres as a function of temperature between 90C and 165C.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Thermodynamic approach to the stability of multi-phase systems. Application to the Y 2 O 3 –Fe system

Oxide-metal systems (OMSs) are important in many practical applications, and therefore, are under extensive studies using a wide range of techniques. The most accurate theoretical approaches are based on density functional theory (DFT), which are limited to ~10 2 atoms. Multi-scale approaches, e.g., DFT+Monte Carlo, are often used to model OMSs at the atomic level. These approaches can describe qualitatively the kinetics of some processes but not the overall stability of OMSs. In this paper, we propose a thermodynamic approach to study equilibrium in multiphase systems, which can be sequentially enhanced by considering different defects and microstructures. We estimate the thermodynamic equilibrium by minimization the free energy of the whole multiphase system using a limited set of defects and microstructural objects for which the properties are calculated by DFT. As an example, we consider Y 2 O 3 +bcc Fe with vacancies in both the Y 2 O 3 and bcc Fe phases, Y substitutions and O interstitials in Fe, Fe impurities and antisite defects in Y 2 O 3 . The output of these calculations is the thermal equilibrium concentration of all the defects for a particular temperature and composition. The results obtained confirmed the high temperature stability of yttria in iron. As a result, model development towards more accurate calculations is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Data Driven Smart Proxy for CFD: Application of Big Data Analytics & Machine Learning in Computational Fluid Dynamics (Part 1)

Simulation technologies can reduce the time and cost of the development and deployment of advanced technologies and allow rapid scale-up of these technologies for fossil fuel based energy systems. However, to ensure their usefulness in practice, the credibility of the simulations needs to be established with Uncertainty Quantification (UQ) methods. National Energy Technology Laboratory (NETL) has been applying non-intrusive UQ methodologies to categorize and quantify uncertainties in CFD simulations of gas-solid multiphase flows. To reduce the computational cost associated with gas-solid flow simulations required for UQ analysis, techniques commonly used in the area of Artificial Intelligence (AI) and Data Mining (DM) are used to construct smart proxy models, which can reduce the computational cost of conducting large number of multiphase CFD simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗
NASA NTRS2018

Aerosols in Atmospheric Chemistry and Biogeochemical Cycles of Nutrients

Atmospheric aerosols have complex and variable compositions and properties.While scientific interest is centered on the health and climatic effects of atmospheric aerosols, insufficient attention is given to their involvement in multiphase chemistry that alters their contribution as carriers of nutrients in ecosystems. However, there is experimental proof that the nutrient equilibria of both land and marine ecosystems have been disturbed during the Anthropocene period. This review study first summarizes our current understanding of aerosol chemical processing in the atmosphere as relevant to biogeochemical cycles. Then it binds together results of recent modeling studies based on laboratory and field experiments, focusing on the organic and dust components of aerosols that account for multiphase chemistry, aerosol ageing in the atmosphere, nutrient (N, P, Fe) emissions, atmospheric transport, transformation and deposition. The human-driven contribution to atmospheric deposition of these nutrients, derived by global simulations using past and future anthropogenic emissions of pollutants, is put into perspective with regard to potential changes in nutrient limitations and biodiversity. Atmospheric deposition of nutrients has been suggested to result in human-induced ecosystem limitations with regard to specific nutrients. Such modifications favor the development of certain species against others and affect the overall functioning of ecosystems. Organic forms of nutrients are found to contribute to the atmospheric deposition of the nutrients N, P and Fe by 20%-40%, 35%-45% and 7%-18%, respectively. These have the potential to be key components of the biogeochemical cycles since there is initial proof of their bioavailability to ecosystems. Bioaerosols have been found to make a significant contribution to atmospheric sources of N and P, indicating potentially significant interactions between terrestrial and marine ecosystems. These results deserve further experimental and modeling studies to reduce uncertainties and understand the feedbacks induced by atmospheric deposition of nutrients to ecosystems.

dust↗