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 955 records · Page 53

Effect of subsurface soil moisture variability and atmospheric conditions on methane gas migration in shallow subsurface

We report a major concern resulting from the increased use and production of natural gas has been how to mitigate fugitive greenhouse gas emissions (predominantly methane) from natural gas infrastructure (e.g., leaky shallow pipelines). Subsurface migration and atmospheric loading of methane from pipeline leakage is controlled by source configurations and subsurface soil conditions (e.g., soil heterogeneity and soil moisture) and are further affected by atmospheric conditions (e.g., wind and temperature). However, the transport and attenuation of methane under varying subsurface and atmospheric conditions are poorly understood, making it difficult to estimate leakage fluxes from methane concentration measurements at and above the soil surface. Based on a series of controlled bench-scale experiments using a large porous media tank interfaced with an open-return wind tunnel, this study investigated multiphase processes controlling migration of methane from a point source representing a buried pipeline leaking at fixed flow rate (kg/s) under various saturation and soil-texture conditions. In addition, potential effects of atmospheric boundary controls, wind (0.5 and 2.0 m s -1 ) and temperature (22 and 35 °C), were also examined. Results showed the distinct effects of soil heterogeneity and, to a varying degree, of soil moisture on surface methane concentrations. In addition, results also showed the pronounced effects of wind and, to a lesser degree, of temperature on surface methane concentrations in the presence of varying soil and moisture conditions. The observed subsurface methane profiles were simulated using the multiphase transport simulator TOUGH2-EOS7CA. Observed agreement between measured and simulated data demonstrates that for the conditions studied, multiphase migration of a multicomponent gas mixture (including methane) under density-dependent flow can be adequately represented with a Fickian advection-diffusion (or dispersion) model (ADM) framework. The dominant effect of saturation over the soil texture, could also be inferred from numerical characterization.

09 BIOMASS FUELS↗

Dynamic evolution of liquid-liquid phase separation during continuous cooling

Solidification from a multiphase fluid involves many unknown quantities due to the difficulty of predicting the impact of fluid flow on chemical partitioning. Real-time x-ray radiography has been used to observe liquideliquid phase separation in Al90In10 prior to solidification. Quantitative image analysis has been used to measure the motion and population characteristics of the dispersed indium-rich liquid phase during cooling. Here we determine that the droplet growth characteristics resemble well known steady-state coarsening laws with likely enhancement by concurrent growth due to supersaturation. Simplistic views of droplet motion are found to be insufficient until late in the reaction due to a hydrodynamic instability caused by the large density difference between the dispersed and matrix liquid phases.

Alloys↗

Long-term evolution of decaying magnetohydrodynamic turbulence in the multiphase interstellar medium

Supersonic turbulence in the interstellar medium (ISM) is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this general consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to take into account the roles of the global structures of magnetic fields and the ISM. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in a warm neutral medium. The early evolution of turbulent kinetic energy is consistent with previous results for decaying compressible MHD turbulence characterized by rapid energy decay with a power-law form of E∝t {sup –1} and by a short decay time compared with the flow crossing time. If initial magnetic fields are strong and perpendicular to the sheet, however, long-term evolution of the kinetic energy shows that a significant amount of turbulent energy (∼0.2E {sub 0}) still remains even after 10 flow crossing times for models with periodic boundary conditions. The decay rate is also greatly reduced as the field strength increases for such initial and boundary conditions, but not if the boundary conditions are those for a completely isolated sheet. We analyze velocity power spectra of the remaining turbulence to show that in-plane, incompressible motions parallel to the sheet dominate at later times.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Overloading Risk Mitigation With a Semi-Analytical Probabilistic Model on Branch Current

A semi-analytical formulation is presented in this paper for the probability computation of branch current in multiphase systems. The developed formula is derived based on the linear power flow model in rectangular coordinates. The system uncertainty injections can be renewable energy resources or loads and are modeled using a Gaussian mixture model (GMM). The developed formula can be used to compute the line current violation probability as well as integrate into optimal power flow problem as chance-constraint relaxation. Here, the proposed formula is first compared with the Matlab embedded numerical integration function to show its performance. Besides, the semi-analytical formula is validated and compared with the Monte Carlo simulation method using the IEEE 123-bus system, EPRI Ckt5, and Ckt7 systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of an Atomization Methodology for Spray Combustion

In liquid rocket propulsion, the knowledge and the understanding of liquid-gas interfacial phenomena are very important. This is important for predicting the onset of cavitation occurring in swirl injection elements used in STME, as well as atomization processes in shear-induced injectors (co-axial) and impinging injector elements. From the fact that all the physical processes including droplet size distribution, droplet dispersion, mixing and combustion are controlled by atomization processes, it is expected that the successful incorporation of the volume of fraction (VOF) will greatly enhance the analytical capability of predicting spray combustion processes in liquid-fueled engines. In this paper, a methodology is developed to define and track interfaces between two fluids in non-orthogonal, body-fitted grids using a single fractional volume of fluid (VOF) variable to describe the distribution of the liquid phase in a gas-liquid flow field. This method was implemented in a mature CFD code MAST (Multiphase All-Speed Transient) utilizing the general PISO-C algorithm. For the preliminary study on the analysis of spray combustion and tracking of the interface between two phases, we will report on the progress of the simulation of the instability on the liquid column; the surface wave instability and the droplet breakup from the liquid surface.

Seung, S. P.↗

Modeling of Thermal Performance of Multiphase Nuclear Fuel Cell Under Variable Gravity Conditions

A unique numerical method has been developed to model the dynamic processes of bulk evaporation and condensation processes, associated with internal heat generation and natural convection under different gravity levels. The internal energy formulation, for the bulk liquid-vapor phase change problems in an encapsulated container, was employed. The equations, governing the conservation of mass, momentum and energy for both phases involved in phase change, were solved. The thermal performance of a multiphase uranium tetra-fluoride fuel element under zero gravity, micro-gravity and normal gravity conditions has been investigated. The modeling yielded results including the evolution of the bulk liquid-vapor phase change process, the evolution of the liquid-vapor interface, the formation and development of the liquid film covering the side wall surface, the temperature distribution and the convection flow field in the fuel element. The strong dependence of the thermal performance of such multiphase nuclear fuel cell on the gravity condition has been revealed. Under all three gravity conditions, 0-g, 10(exp -3)-g, and 1-g, the liquid film is formed and covers the entire side wall. The liquid film covering the side wall is more isothermalized at the wall surface, which can prevent the side wall from being over-heated. As the gravity increases, the liquid film is thinner, the temperature gradient is larger across the liquid film and smaller across the vapor phase. This investigation provides valuable information about the thermal performance of multi-phase nuclear fuel element for the potential space and ground applications.

Ding, Z.↗

Individual phase constitutive properties of a TRIP-assisted QP980 steel from a combined synchrotron X-ray diffraction and crystal plasticity approach

Microstructure-based constitutive models for multiphase steels require accurate constitutive properties of the individual phases for component forming and performance simulations. We address this requirement with a combined experimental/theoretical methodology which determines the critical resolved shear stresses and hardening parameters of the constituent phases in QP980, a TRIP assisted steel subject to a two-step quenching and partitioning heat treatment. High energy X-Ray diffraction (HEXRD) from a synchrotron source provided the average lattice strains of the ferrite, martensite, and austenite phases from the measured volume during in situ tensile deformation. The HEXRD data was then input to a computationally efficient, elastic-plastic self-consistent (EPSC) crystal plasticity model which estimated the constitutive parameters of different slip systems for the three phases via a trial-and-error approach. The EPSC-estimated parameters are then input to a finite element crystal plasticity (CPFE) model representing the QP980 tensile sample. The predicted lattice strains and global stress versus strain curves are found to be 8% lower that the EPSC model predicted values and from the HEXRD measurements, respectively. This discrepancy, which is attributed to the stiff secant assumption in the EPSC formulation, is resolved with a second step in which CPFE is used to iteratively refine the EPSC-estimated parameters. Remarkably close agreement is obtained between the theoretically-predicted and experimentally derived flow curve for the QP980 material.

Hu, X. H. (ORCID:0000000277355091)↗

Simulation of Transport Phenomena in Bubble Column Reactors

The production of carbon and hydrogen from methane is a viable pathway for monetizing the vast reserves of US natural gas. Liquid (metals and/or salts) bubble column reactors operating at temperatures greater than 1000° C allow the conversion of methane into hydrogen without CO 2 emission. The low-density carbon floats on top of the liquid surface, allowing an easier removal, compared to other reactor concepts for methane pyrolysis. In this current report, we present results from simulations of transport phenomena in a liquid bubble column reactor. The goal is to investigate and understand parameters that are crucial for bubble column reactor design and scale up. The report presents the formulation and numerical methods used in the multiphase direct numerical simulation software, Quilt. This software was validated for the bubbly flow using numerical and physics based test problems. Quilt was then used to perform simulations of the bubble column reactor as a parametric study, varying the bubble injection parameters. The simulations were visualized for qualitative analysis and followed up with a quantitative analysis. The effect of injection parameters on various quantitative measures of the bubble residence time, interfacial area and motion through the column are presented.

03 NATURAL GAS↗

Comprehensive 3D Multiphysics Model on Electrochemical Recovery of O 2 from Metabolic CO 2 at the International Space Station (ISS)

The International Space Station (ISS) is presently equipped with an elaborate, heavy, and high-power consuming system that recovers approximately 50% of O 2 from metabolic CO 2 as part of the atmospheric revitalization (AR) at the ISS habitat. Future long-duration missions will require a sustainable and efficient system capable of yielding a minimum of 75% O 2 recovery to reach the self-sufficiency required for long space missions beyond earth’s low orbit. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to not only increase significantly current O 2 recovery efficiency, improving self-sufficiency on AR at the ISS habitat and future long missions, but also reduce the complexity of the system. The authors have developed and deployed a comprehensive 3D multiphysics model that thoroughly replicates the actual configuration and fluid/material domains of the MFECR. The coupled physics in this multiphysics model include multicomponent-multiphase electrochemical-driven reactions, non-ideal mass transport mechanism, free and porous flow, heat transfer, CO 2 solubility on alkaline electrolyte, water condensation on porous medium, and DC electrical current generation along with Joule heating effect. This model is aimed to conduct quantitive benchmark on three different MFECR’s layouts, one without serpentine paths (plain) and two with serpentines leading to four and twelve paths respectively. Once experimental data is generated via a test matrix of 200 tests, the model will be validated to conduct MFECR’s process optimization and revalidate the quantitive benchmark on three different MFECR’s layouts.

Jesus A. Dominguez↗

Preconditioners for multiphase poromechanics with strong capillarity

This paper aims to enhance the performance of Newton–Krylov solvers for coupled poromechanical problems with two-phase flow. In particular, we investigate the impact of capillary pressure on preconditioning strategies. Capillarity complicates the coupling between the solid deformation and fluid pressure degrees of freedom, as well as increases the nonlinearity of the system. Depending on the capillary pressure relation used in the constitutive formulation, the flow equations may exhibit a spectrum of advection-dominated to diffusion-dominated behavior. We propose preconditioning approaches that account for this behavior and lead to robust numerical performance within a broad range of regimes.

42 ENGINEERING↗

Elliptic flow coalescence to identify the $f_0$(980) content

Here we use a simple coalescence model to generate $f_0$(980) particles for three configurations: a $s\bar{s}$ meson, a $u\bar{u}s\bar{s}$ tetraquark, and a $K^+K^–$ molecule. The phase-space information of the coalescing constituents is taken from a multiphase transport (AMPT) simulation of heavy-ion collisions. It is shown that the number of constituent quarks scaling of the elliptic flow anisotropy can be used to discern $s\bar{s}$ from $u\bar{u}s\bar{s}$ and $K^+K^–$ configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A filter-dependent granular temperature model from large-scale CFD-DEM data

The computational study of strongly-coupled, gas–solid flows at scales relevant to most environmental and engineering applications requires the use of ‘coarse-grained’ methodologies such as the two-fluid model, particle-in-cell approach or the multiphase Reynolds Averaged Navier–Stokes equations. While these strategies enable computations at desirable length- and time-scales, they rely heavily on models to capture important flow physics that occur at scales smaller than the mesh. To date, the models that do exist are based on a limited set of flow conditions, such as very dilute particle phase. To this end, we leverage a large-scale repository of CFD-DEM data to develop filter-size dependent models for the mean variance in particle volume fraction, a quantity commonly used to assess the degree of clustering, and the granular temperature, a key quantity for accurately predicting gas–solid flows. In conclusion, because of its filter-size dependence, the granular temperature model can be directly translated to coarse-grained approaches and tied directly to grid size.

AMReX↗

Reanalysis of x ray emission from M87 II: Multiphase models

In a previous paper by the authors (hereafter TB), it was shown that the hot intracluster gas around M87 could not be adequately explained in terms of a spherically symmetric, single phase model. It was found that although data from the Einstein satellite High Resolution Imager (HRI), Imaging Proportional Counter (IPC), and Focal Point Crystal Spectrometer (FPCS) could indeed be simultaneously explained by a single phase model, data from the Solid State Spectrometer (SSS) and optically determined mass estimates could not be similarly explained. A qualitative discussion was then given to indicate how the adoption of a specific multiphase model could plausibly improve upon the single phase model, but a detailed discussion was delayed until now. The x ray data is reanalyzed for the gas around M87 in the same spirit as the work in TB. That is, by assuming several cooling flow models, the shortcomings are improved upon. The procedure of TB was adopted and the data reanalyzed in the following manner: a multiphase model is assumed which can compute the surface brightness as seen by HRI and IPC. Also, the line fluxes can be computed of the lines seen by the FPCS and the mass profile after assuming hydrostatic equilibrium. The parameters of the model are then adjusted to fit the data. A check is then made to see whether the resulting model is consistent with the SSS spectrum and the equivalent width of the 7 keV complex of Fe lines as seen by large field of view instruments.

Tsai, John C.↗

THE LAUNCHING OF COLD CLOUDS BY GALAXY OUTFLOWS. I. HYDRODYNAMIC INTERACTIONS WITH RADIATIVE COOLING

To better understand the nature of the multiphase material found in outflowing galaxies, we study the evolution of cold clouds embedded in flows of hot and fast material. Using a suite of adaptive mesh refinement simulations that include radiative cooling, we investigate both cloud mass loss and cloud acceleration under the full range of conditions observed in galaxy outflows. The simulations are designed to track the cloud center of mass, enabling us to study the cloud evolution at long disruption times. For supersonic flows, a Mach cone forms around the cloud, which damps the Kelvin–Helmholtz instability but also establishes a streamwise pressure gradient that stretches the cloud apart. If time is expressed in units of the cloud crushing time, both the cloud lifetime and the cloud acceleration rate are independent of cloud radius, and we find simple scalings for these quantities as a function of the Mach number of the external medium. A resolution study suggests that our simulations accurately describe the evolution of cold clouds in the absence of thermal conduction and magnetic fields, physical processes whose roles will be studied in forthcoming papers.

79 ASTRONOMY AND ASTROPHYSICS↗

Production of Cool Gas in Thermally Driven Outflows

Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the Cholla Galactic OutfLow Simulations suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ( v ~ 1000 km s -1 ), cool (10 4 K) gas observed in absorption-line studies of outflows around some star-forming galaxies. By varying the mass loading and geometry of the simulated winds, we identify a region of parameter space that leads to cool gas in outflows. In particular, when using an analytically motivated central feedback model, we find that cooling flows can be produced with reasonable mass-loading rates (${\dot{M}}_{\mathrm{wind}}/{\dot{M}}_{\mathrm{SFR}}\sim 0.5$), provided that the star formation rate surface density is high. When a more realistic clustered feedback model is applied, destruction of high-density clouds near the disk and interactions between different outflow regions indicate that lower mass-loading rates of the hot gas within the feedback region may still produce multiphase outflows. These results suggest an origin for fast-moving cool gas in outflows that does not rely on directly accelerating cool gas from the interstellar medium. Finally, these cooling flows may additionally provide an explanation for the multiphase gas ubiquitously observed in the halos of star-forming galaxies at low redshift.

79 ASTRONOMY AND ASTROPHYSICS↗