Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase flows”

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 415 records · Page 23

Numerical Analysis of Liquid Metal MHD Flow and Heat Transfer for Open-Surface Li Divertor in FNSF

Within the ongoing U.S.-based program on the development of liquid metal plasma-facing components, numerical simulations and analyses are performed to address the feasibility of the open-surface Li divertor. In the previous scoping studies (Smolentsev, 2021), heat-removal capabilities of the divertor were assessed using a simplified flow model for a slug-type velocity profile and constant flow thickness. Here, new analyses take into account forces acting on the flowing Li layer. Three reduced-order mathematical models are applied under the conditions of the U.S. Fusion Nuclear Science Facility (FNSF) to access magnetohydrodynamic (MHD) flow development effects, velocity distribution, and surface waves: 1) fully developed MHD flow; 2) quasi-2-D developing MHD flow; and 3) multiphase MHD flow. The obtained results for MHD flows and the surface heat flux computed with the plasma code scrape-off layer plasma simulation for ITER (SOLPS-ITER) are then used as input data to compute the temperature distribution in the divertor by solving the convection–diffusion energy equation.

Smolentsev, Sergey↗

Alquimia v1.0: a generic interface to biogeochemical codes – a tool for interoperable development, prototyping and benchmarking for multiphysics simulators

Alquimia v1.0 is a generic interface to geochemical solvers that facilitates development of multiphysics simulators by enabling code coupling, prototyping and benchmarking. The interface enforces the function arguments and their types for setting up, solving, serving up output data and carrying out other common auxiliary tasks while providing a set of structures for data transfer between the multiphysics code driving the simulation and the geochemical solver. Alquimia relies on a single-cell approach that permits operator splitting coupling and parallel computation. We describe the implementation in Alquimia of two widely used open-source codes that perform geochemical calculations: PFLOTRAN and CrunchFlow. We then exemplify its use for the implementation and simulation of reactive transport in porous media by two open-source flow and transport simulators: Amanzi and ParFlow. We also demonstrate its use for the simulation of coupled processes in novel multiphysics applications including the effect of multiphase flow on reaction rates at the pore scale with OpenFOAM, the role of complex biogeochemical processes in land surface models such as the E3SM Land Model (ELM) and the impact of surface–subsurface hydrological interactions on hydrogeochemical export from watersheds with the Advanced Terrestrial Simulator (ATS). These applications make it apparent that the availability of a well-defined yet flexible interface has the potential to improve the software development workflow, freeing up resources to focus on advances in process models and mechanistic understanding of coupled problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool for the analysis of flowfields near the injector face has been developed and used to analyze the Fastrac engine. The objective was to produce a computational design tool which was detailed enough to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows. To obtain a model which could be used to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe liquid and vapor sub- and super-critical flows, the model included thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. A homogeneous model was constructed such that the local state of the flow was determined directly, i.e. the quality of the flow was calculated. Such a model does not identify drops or their distribution, but it does allow the flow along the injector face and into the acoustic cavity to be predicted. It also allows the film coolant flow to be accurately described. The initial evaluation of the injector code was made by simulating cold flow from an unlike injector element and from a like-on-like overlapping fan (LOL) injector element. The predicted mass flux distributions of these injector elements compared well to cold flow test results. These are the same cold flow tests which serve as the data base for the JANNAF performance prediction codes. The flux distributions 1 inch downstream of the injector face are very similar; the differences were somewhat larger at further distances from the faceplate. Since the cold flow testing did not achieve good mass balances when integrations across the entire fan were made, the CFD simulation appears to be reasonable alternative to future cold flow testing. To simulate the Fastrac, an RP-1/LOX combustion model must be chosen. This submodel must be relatively simple to accomplish three-dimensional, multiphase flow simulations. Single RP-1 pyrolysis and partial oxidation steps were chosen and the combustion was completed with the wet CO mechanism. Soot was also formed with a single global reaction. To validate the combustion submodel, global data from gas generator tests and from subscale motor test were used to predict qualitatively correct mean molecular weights, temperature, and soot levels. Because such tests do not provide general kinetics rates, the methodology is not necessarily appropriate for other than rocket type flows conditions. Soot predictions were made so that radiation heating to the motor walls can be made. These initial studies of the Fastrac were for a small region close to the injector face and chamber wall which included a segment of the acoustic cavity. The region analyzed includes 11 individual orifice holes to represent the LOL elements and the H2 film coolant holes. Typical results of this simulation are shown in Figure 1. At this point the only available test data to verify the predictions are temperatures measured in the acoustic cavity. These temperatures are in reasonable agreement at about 2000R (1111 K). Future work is expected to include improving the computational efficiency or the CFD model and/or using more computer capacity than the single Pentium PC with which these simulations were made.

Farmer, Richard C.↗

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗

Second order upwind Lagrangian particle method for Euler equations

A new second order upwind Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwind second-order particle-based algorithm with limiter, providing accuracy and long term stability, and (c) accurate resolution of states at free interfaces. In conclusion, numerical verification tests demonstrating the convergence order for fixed domain and free surface problems are presented.

97 MATHEMATICS AND COMPUTING↗

Results of the Workshop on Two-Phase Flow, Fluid Stability and Dynamics: Issues in Power, Propulsion, and Advanced Life Support Systems

The Two-phase Flow, Fluid Stability and Dynamics Workshop was held on May 15, 2003 in Cleveland, Ohio to define a coherent scientific research plan and roadmap that addresses the multiphase fluid problems associated with NASA s technology development program. The workshop participants, from academia, industry and government, prioritized various multiphase issues and generated a research plan and roadmap to resolve them. This report presents a prioritization of the various multiphase flow and fluid stability phenomena related primarily to power, propulsion, fluid and thermal management and advanced life support; and a plan to address these issues in a logical and timely fashion using analysis, ground-based and space-flight experiments.

McQuillen, John↗

Solving multiphysics-based inverse problems with learned surrogates and constraints

Abstract Solving multiphysics-based inverse problems for geological carbon storage monitoring can be challenging when multimodal time-lapse data are expensive to collect and costly to simulate numerically. We overcome these challenges by combining computationally cheap learned surrogates with learned constraints. Not only does this combination lead to vastly improved inversions for the important fluid-flow property, permeability, it also provides a natural platform for inverting multimodal data including well measurements and active-source time-lapse seismic data. By adding a learned constraint, we arrive at a computationally feasible inversion approach that remains accurate. This is accomplished by including a trained deep neural network, known as a normalizing flow, which forces the model iterates to remain in-distribution, thereby safeguarding the accuracy of trained Fourier neural operators that act as surrogates for the computationally expensive multiphase flow simulations involving partial differential equation solves. By means of carefully selected experiments, centered around the problem of geological carbon storage, we demonstrate the efficacy of the proposed constrained optimization method on two different data modalities, namely time-lapse well and time-lapse seismic data. While permeability inversions from both these two modalities have their pluses and minuses, their joint inversion benefits from either, yielding valuable superior permeability inversions and CO 2 plume predictions near, and far away, from the monitoring wells.

Yin, Ziyi (ORCID:0000000250248771)↗

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗

A Design Tool for Liquid Rocket Engine Injectors

A practical design tool which emphasizes the analysis of flowfields near the injector face of liquid rocket engines has been developed and used to simulate preliminary configurations of NASA's Fastrac and vortex engines. This computational design tool is sufficiently detailed to predict the interactive effects of injector element impingement angles and points and the momenta of the individual orifice flows and the combusting flow which results. In order to simulate a significant number of individual orifices, a homogeneous computational fluid dynamics model was developed. To describe sub- and supercritical liquid and vapor flows, the model utilized thermal and caloric equations of state which were valid over a wide range of pressures and temperatures. The model was constructed such that the local quality of the flow was determined directly. Since both the Fastrac and vortex engines utilize RP-1/LOX propellants, a simplified hydrocarbon combustion model was devised in order to accomplish three-dimensional, multiphase flow simulations. Such a model does not identify drops or their distribution, but it does allow the recirculating flow along the injector face and into the acoustic cavity and the film coolant flow to be accurately predicted.

Farmer, R.↗

A Numerical Study of Two-Phase Flow and Interfacial Mass Transfer in a Wetted Wall Column for Counter-Current Gas Absorption

We present a numerical study to investigate the effects of surface waves of falling films on the mass transfer from flowing gas phase into the liquid film. The model employs the volume of fluid method to explicitly track the gas-liquid phase interface and a one-fluid formulation to model the multiphase flow coupled with mass transfer between the two phases. The surface waves are generated by the time-varying injection rate of the falling film. The effects of surface wave frequencies and amplitudes on the overall mass transfer also are investigated for a counter-current gas-liquid flow of oxygen and water. Because of the solubility difference in the two phases, the model can capture the concentration discontinuity across the liquid-gas interface. Good agreement is obtained between the numerical results and experimental data reported in the literature. By systematically varying the frequencies and amplitudes of flow rate for the falling film, the surface waves generated on the film surface enhance mass transfer in general. Here, the model prediction shows a proportional increase in mass transfer for small frequencies of surface waves, but saturation at higher frequencies.

42 ENGINEERING↗

CFD validation experiments for internal flows

Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis are described. The material presented summarized the research in 3 areas: Inlets, ducts and nozzles; Turbomachinery; and Chemically reacting flows. The specific validation activities are concerned with shock boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.

Povinelli, Louis A.↗

Direct Effect of Solvent Viscosity on the Physical Mass Transfer for Wavy Film Flow in a Packed Column

ABSTRACT Post-combustion carbon capture via solvent absorption in a packed column can be a promising technology for mitigating greenhouse gas emissions. The mass transfer between gas and liquid solvent phases plays a critical role in determining the packing height of the packed column. Recent experiments [1] involving physical mass transfer for aqueous amine in a packed column have distinguished the direct and indirect impacts of viscosity (µ_L) on the mass transfer coefficient ?( k?_L). The direct impact of viscosity on the mass transfer coefficient has been found to be significantly higher for a packed column ?(k?_L~µ_L^(-0.38)) than a wetted wall column (WWC) ?(k?_L~µ_L^(-0.17)). To explain the deviation from the well-controlled uniform film flow in a WWC, we propose a mechanism involving the unavoidable wavy film flow and eddy-enhanced mass transfer in a packed column. A mathematical analysis based on fundamental flow and mass transfer equations is presented to quantify this deviation. The correlation for mass transfer coefficient derived from the mathematical analysis is found to match the corresponding one derived from the experiments for a packed column. Multiphase flow simulations using the Volume of Fluid method also are conducted to verify the direct impact of viscosity on the mass transfer coefficient for both uniform and wavy film flows. The mass transfer coefficient calculated from simulations shows steeper variation with viscosity for wavy film compared to the uniform film, further confirming our proposed theory. A similar relation ?(k?_L~µ^(-0.38)) for wavy film flow in a packed column is observed for aqueous solvent from both theory, experiment, and computer simulations.

Mass transfer, viscosity, eddy diffusion↗

CFD validation experiments for internal flows

Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis Research Center are described. The material presented summarizes the research in three areas: Inlets, Ducts and Nozzles; Turbomachinery; and Chemically Reacting Flows. The specific validation activities are concerned with shock-boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.

Povinelli, Louis A.↗

The effect of CO 2 -induced dissolution on flow properties in Indiana Limestone: An in situ synchrotron X-ray micro-tomography study

The injection of CO 2 -rich fluids in carbonate rocks results in an evolution of the pore space, with consequent changes in the hydraulic properties of the reservoir; how these properties evolve, particularly for parameters relevant to multiphase flow e.g. P c (s), remains a topic of active research despite several decades of study. We have carried out an in situ synchrotron X-ray microtomography experiment to monitor pore structure evolution during dissolution of an Indiana Limestone core; the experiment involved flowing CO 2 -saturated water through the core for 36 h and resulted in 10 volumes corresponding to different temporal stages of the dissolution process. The injection parameters corresponded to the flow velocities expected near the well-bore region of a shallow aqueous CO 2 injection; fast flow rates with high reactant availability. Analysis of the tomographic data shows flow-enhanced dissolution i.e. channeling, and provides a time-resolved map of pore space alteration. Using the resulting 4D pore space volume, we modeled the evolution of capillary-pressure curves; this exercise demonstrates how pore structure evolution could impact the invasion and remobilization of non-wetting fluids, dramatically decreasing the entry pressure and the P C in some parts of the sample. The modeling of permeability, using a Stokes solver approach, quantified the relationship of porosity vs. permeability; we found that a modest increase in porosity, especially when the channeling system is more developed, greatly affects permeability. Furthermore, these results demonstrate how movement of CO 2 saturated brine near injected plumes might alter drainage dynamics near the plume boundary, thus leading to mobilization across subtle capillary barriers.

58 GEOSCIENCES↗

Cellular fluidics

The natural world provides many examples of multiphase transport and reaction processes that have been optimized by evolution. These phenomena take place at multiple length and time scales and typically include gas–liquid–solid interfaces and capillary phenomena in porous media. Many biological and living systems have evolved to optimize fluidic transport. However, living things are exceptionally complex and very difficult to replicate, and human-made microfluidic devices (which are typically planar and enclosed) are highly limited for multiphase process engineering. In this paper, we introduce the concept of cellular fluidics: a platform of unit-cell-based, three-dimensional structures—enabled by emerging 3D printing methods—for the deterministic control of multiphase flow, transport and reaction processes. We show that flow in these structures can be ‘programmed’ through architected design of cell type, size and relative density. We demonstrate gas–liquid transport processes such as transpiration and absorption, using evaporative cooling and CO 2 capture as examples. We design and demonstrate preferential liquid and gas transport pathways in three-dimensional cellular fluidic devices with capillary-driven and actively pumped liquid flow, and present examples of selective metallization of pre-programmed patterns. Our results show that the design and fabrication of architected cellular materials, coupled with analytical and numerical predictions of steady-state and dynamic behaviour of multiphase interfaces, provide deterministic control of fluidic transport in three dimensions. Cellular fluidics may transform the design space for spatial and temporal control of multiphase transport and reaction processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

How does the connectivity of open-framework conglomerates within multi-scale hierarchical fluvial architecture affect oil-sweep efficiency in waterflooding?

Understanding multi-phase fluid flow and transport processes within aquifers, candidate reservoirs for CO 2 sequestration, and petroleum reservoirs requires understanding a diverse set of geologic properties of the aquifer or reservoir, over a wide range of spatial and temporal scales. We focus on multiphase flow dynamics with wetting (e.g., water) and non-wetting (e.g., gas or oil) fluids, with one invading another. This problem is of general interest in a number of fields and is illustrated here by considering the sweep efficiency of oil during a waterflood. Using a relatively fine-resolution grid throughout a relatively large domain in these simulations and probing the results with advanced scientific visualization tools (Reservoir Visualization Analysis [RVA]/ ParaView software) promote a better understanding of how smaller-scale features affect the aggregate behavior at larger scales. We studied the effects on oil-sweep efficiency of the proportion, hierarchical organization, and connectivity of high-permeability open-framework conglomerate (OFC) cross-sets within the multi-scale stratal architecture found in fluvial deposits. We further analyzed oil production rate, water breakthrough time, and spatial and temporal distribution of residual oil saturation. As expected, the effective permeability of the reservoir exhibits large-scale anisotropy created by the organization of OFC cross-sets within unit bars, and the organization of unit bars within compound- bars. As a result, oil-sweep efficiency critically depends on the direction of the pressure gradient. However, contrary to expectations, the total amount of trapped oil due to the effect of capillary trapping does not depend on the magnitude of the pressure gradient within the examined range. Hence the pressure difference between production and injection wells does not affect sweep efficiency; although the spatial distribution of oil remaining in the reservoir depends on this value. Whether or not clusters of connected OFC span the domain affects only the absolute rate of oil production—not sweep efficiency.

58 GEOSCIENCES↗

SFAI and Related Technologies

Researchers at Los Alamos have developed Swept-Frequency Acoustic Interferometry (SFAI) a technology that enables businesses to have greater control over their industrial processes, allowing for significantly improved efficiency by giving operators vital composition and flow rate information in real-time. This technology fills the information gap by providing this data non-invasively and at a fraction of the cost of traditional composition meters. SFAI uses engineered ultrasound and advanced signal processing to characterize multiphase flow (e.g., liquids, gasses, multi-phase mixtures) and other flows inside sealed containers and flow pipes. By being both low-cost and non-invasive, this technology is set to disrupt the current oil field monitoring state-of-the-art by providing individual well-level composition data, compared to the multi-well macroscopic information that is currently available. We are seeking a commercialization partner to license the technology or fund us through a Cooperative Research and Development Agreement (CRADA) to a further develop a usable product for commercial purposes.

03 NATURAL GAS↗

Countercurrent flow characteristics of next generation solvent in novel 3D printed columns for carbon capture

Solvent based absorption for carbon dioxide capture in a packed column is being considered as an efficient technology for the decarbonization process of point source. Significant efforts are underway to improve the design of packings with the cutting edge 3D printing technology for efficient carbon capture. Accordingly, multiphase flow studies were conducted to assess the performance of novel 3D printed columns with various triply periodic minimal surface (TPMS) designs. The effects of solvent properties, liquid and gas loads on the performance of TPMS columns are extensively explored. Hydrodynamics of the potential water-lean solvents (EEMPA) as well as aqueous monoethanolamine (MEA) solvent for carbon capture are evaluated and compared. The interfacial area and liquid holdup increase with increasing liquid loads (q L ) for TPMS columns. Schwarz column consistently shows the highest liquid holdup value. The EEMPA exhibits higher values for the interfacial area (10–15 %) and liquid holdup (~4 %) in comparison to the MEA. TPMS columns exhibit the intermediate value of dry pressure drop between the random and the structure packed columns. Among selected TPMS columns, the gyroid packing shows the lowest pressure drop. The gas load has marginal impact on the interfacial area at lower value while a higher gas load leads to column’s flooding. Prior to the flooding, the interfacial area in TMPS packings rises with increased gas load at a fixed liquid load except Schwarz packing where interfacial area is incentive to the gas load. Additionally, the liquid holdup and wet pressure drop rise as gas load increases in TPMS columns. Overall, Fisher Koch packing is more susceptible to flood as compared to other TPMS packings. Flow regimes: loading and flooding are also delineated in the TPMS packings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗