Search NASA⌕ Search

SEARCH · Search NASA

Results for “multiphase modeling”

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 325 records · Page 18

Gas Evolution in Water Electrolysis

Gas bubbles generated by the hydrogen evolution reaction and oxygen evolution reaction during water electrolysis influence the energy conversion efficiency of hydrogen production. Here we survey what is known about the interaction of gas bubbles and electrode surfaces and the influence of gas evolution on practicable devices used for water electrolysis. Here, we survey the physical processes occurring during the life cycle of a bubble, summarize techniques used to characterize gas evolution phenomena in-situ and in practical device environments, and discuss ways that electrodes can be tailored to facilitate gas removal at high current densities. Lastly, we summarize efforts to model the behavior of individual gas bubbles and multiphase flows produced at gas evolving electrodes. We conclude our review with a short summary of outstanding questions which could be answered by future efforts to characterize gas evolution in electrochemical device environments or by improved simulations of multiphase flows.

Bubbles↗

A Relaxed PV Bus Model in Linear Power Flow

The emerging multimodal active energy resource units make it necessary to introduce constant voltage amplitude buses (PV bus) in the operation of modern distribution systems. To enable the analysis of PV bus using linear power flow models, this letter proposes a relaxed PV bus model for solving power flow related problems using an indirect modeling approach. The developed PV bus model can be applied to linear power flow models in both rectangular and polar coordinates. Besides, the approximation error of the relaxed PV bus model can also be controlled and quantified. Lastly, the proposed relaxed PV bus model is validated using different test systems. According to the numerical studies, the maximum error on the PV bus approximation in all case studies is only 0.61%, and the maximum error of the supporting power at the PV buses is 6.5%.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Multiphase representation of hotspot formation and growth

Multiphase reacting flow approach has been adopted for modeling of hotspot formation and growth on sub-zonal scale. Mechanical response of the high explosive (HE) is modeled by the Visco-Plastic SCRAM model which includes plastic work, crack growth, and frictional heating along the cracks. In this work a sub-zonal regions with high temperature, hotspots, are modeled as a new phase. Formation and growth of the hotspot phase from the cold phase is modeled as a phase change. Each phase has its own temperature and pressure. Chemical reactions also progress in both phases. The growth of hotspots are determined by the solution of the heat conduction equation with appropriate heat sources due to friction and plastic heating. Sample calculations using the BABBO code show successful detonation calculations with hotspot formation and growth.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Assessment of Radiation Scattering by Dust in Mars Entry

The impact of gas phase radiation scattering by dust particle presence in Mars entry scenarios is assessed using a suite of numerical tools. A hypersonic gas and dilute particle multiphase flow is simulated using a mixed continuum model and an efficient size distribution quadrature scheme. Radiative transport is simulated with a hybrid ray-tracing and P1 spherical harmonics model. The NASA HARA radiation code and a Lorentz-Mie scattering model are used to predict the gas-particle mixture radiative properties. Application of the numerical suite to a representative case based on the Mars 2020 entry shows that, even at conditions in excess of the strongest recorded dust storm, the augmentation in surface radiative heat flux due to scattering alone is lower than 1%, and therefore negligible. Predictions with even higher dust loading conditions confirm that even highly uncertain dust mass profiles are unlikely to be a concern for the scattering augmentation mechanism.

Andrew Hinkle↗

Using density functional theory to construct multiphase equations of state: tin as an example

We present a general methodology for using density functional theory (DFT) calculations as a basis for constructing multiphase equations of state (EOS). Focusing on tin as an example, we discuss the full process of generating a tabular EOS, starting from first-principles calculations and arriving at a full EOS in agreement with existing experimental data. We begin by describing DFT calculations for the five solid phases and liquid phase of tin, including cold curve, phonon, and DFT-based molecular dynamics calculations. We then discuss methods for incorporating DFT results into materials models used in OpenSesame, the program used to create a full tabular multiphase EOS. Next, we outline a general strategy for adjusting model parameters in OpenSesame to fit to a variety of experimental data, including isobaric data, isothermal data, shock data, solid-solid phase boundary measurements, and measurements of the melt curve. We end with a discussion of the advantages of using DFT data as a foundation for EOS generation and discuss potential strategies for automated EOS generation based on the fitting process highlighted in this work.

74 ATOMIC AND MOLECULAR PHYSICS↗

AMORE-Isoprene v1.0: a new reduced mechanism for gas-phase isoprene oxidation

Abstract. Gas-phase oxidation of isoprene by ozone (O3) and the hydroxyl (OH) and nitrate (NO3) radicals significantly impacts tropospheric oxidant levels and secondary organic aerosol formation. The most comprehensive and up-to-date chemical mechanism for isoprene oxidation consists of several hundred species and over 800 reactions. Therefore, the computational expense of including the entire mechanism in large-scale atmospheric chemical transport models is usually prohibitive, and most models employ reduced isoprene mechanisms ranging in size from ∼ 10 to ∼ 200 species. We have developed a new reduced isoprene oxidation mechanism using a directed-graph path-based automated model reduction approach, with minimal manual adjustment of the output mechanism. The approach takes as inputs a full isoprene oxidation mechanism, the environmental parameter space, and a list of priority species which are protected from elimination during the reduction process. Our reduced mechanism, AMORE-Isoprene (where AMORE stands for Automated Model Reduction), consists of 12 species which are unique to the isoprene mechanism as well as 22 reactions. We demonstrate its performance in a box model in comparison with experimental data from the literature and other current isoprene oxidation mechanisms. AMORE-Isoprene's performance with respect to predicting the time evolution of isoprene oxidation products, including isoprene epoxydiols (IEPOX) and formaldehyde, is favorable compared with other similarly sized mechanisms. When AMORE-Isoprene is included in the Community Regional Atmospheric Chemistry Multiphase Mechanism 1.0 (CRACMM1AMORE) in the Community Multiscale Air Quality Model (CMAQ, v5.3.3), O3 and formaldehyde agreement with Environmental Protection Agency (EPA) Air Quality System observations is improved. O3 bias is reduced by 3.4 ppb under daytime conditions for O3 concentrations over 50 ppb. Formaldehyde bias is reduced by 0.26 ppb on average for all formaldehyde measurements compared with the base CRACMM1. There was no significant change in computation time between CRACMM1AMORE and the base CRACMM. AMORE-Isoprene shows a 35 % improvement in agreement between simulated IEPOX concentrations and chamber data over the base CRACMM1 mechanism when compared in the Framework for 0-D Atmospheric Modeling (F0AM) box model framework. This work demonstrates a new highly reduced isoprene mechanism and shows the potential value of automated model reduction for complex reaction systems.

58 GEOSCIENCES↗

AMORE-Isoprene v1.0: A New Reduced Mechanism for Gas-Phase Isoprene Oxidation

Gas-phase oxidation of isoprene by ozone (O 3 ) and the hydroxyl (OH) and nitrate (NO 3 ) radicals significantly impacts tropospheric oxidant levels and secondary organic aerosol formation. The most comprehensive and up-to-date chemical mechanism for isoprene oxidation consists of several hundred species and over 800 reactions. Therefore, the computational expense of including the entire mechanism in large-scale atmospheric chemical transport models is usually prohibitive, and most models employ reduced isoprene mechanisms ranging in size from ∼ 10 to ∼ 200 species. We have developed a new reduced isoprene oxidation mechanism using a directed-graph path-based automated model reduction approach, with minimal manual adjustment of the output mechanism. The approach takes as inputs a full isoprene oxidation mechanism, the environmental parameter space, and a list of priority species which are protected from elimination during the reduction process. Our reduced mechanism, AMORE-Isoprene (where AMORE stands for Automated Model Reduction), consists of 12 species which are unique to the isoprene mechanism as well as 22 reactions. We demonstrate its performance in a box model in comparison with experimental data from the literature and other current isoprene oxidation mechanisms. AMORE-Isoprene's performance with respect to predicting the time evolution of isoprene oxidation products, including isoprene epoxydiols (IEPOX) and formaldehyde, is favorable compared with other similarly sized mechanisms. When AMORE-Isoprene is included in the Community Regional Atmospheric Chemistry Multiphase Mechanism 1.0 (CRACMM1AMORE) in the Community Multiscale Air Quality Model (CMAQ, v5.3.3), O 3 and formaldehyde agreement with Environmental Protection Agency (EPA) Air Quality System observations is improved. O 3 bias is reduced by 3.4 ppb under daytime conditions for O 3 concentrations over 50 ppb. Formaldehyde bias is reduced by 0.26 ppb on average for all formaldehyde measurements compared with the base CRACMM1. There was no significant change in computation time between CRACMM1AMORE and the base CRACMM. AMORE-Isoprene shows a 35 % improvement in agreement between simulated IEPOX concentrations and chamber data over the base CRACMM1 mechanism when compared in the Framework for 0-D Atmospheric Modeling (F0AM) box model framework. This work demonstrates a new highly reduced isoprene mechanism and shows the potential value of automated model reduction for complex reaction systems.

Automated Model Reduction↗

A model to predict the conditions for liquid drop breakup and the resultant mean fragment size

The potential significance of drop fragmentation in sprays and other propulsion-related multiphase flows has been noted in the literature. This has motivated recent experimental and theoretical works to: better understand the fundamentals of physics of drop breakup processes, and develop models of drop fragmentation suitable for use in multiphase flow codes. The works summarized below aim to contribute to both sides of this two-pronged attack.

Wert, K. L.↗

Experimentally Driven Automated Machine-Learned Interatomic Potential for a Refractory Oxide

Understanding the structure and properties of refractory oxides is critical for high temperature applications. In this work, a combined experimental and simulation approach uses an automated closed loop via an active learner, which is initialized by x-ray and neutron diffraction measurements, and sequentially improves a machine-learning model until the experimentally predetermined phase space is covered. Furthermore, a multiphase potential is generated for a canonical example of the archetypal refractory oxide, HfO 2 , by drawing a minimum number of training configurations from room temperature to the liquid state at similar to 2900 degrees C. The method significantly reduces model development time and human effort.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

TOUGH-RFPA: Coupled thermal-hydraulic-mechanical Rock Failure Process Analysis with application to deep geothermal wells

This paper presents a coupled thermal-hydraulic-mechanical (THM) model for simulating failure processes in geological media. The coupling between stress/strain and heat transfer describing how thermal expansion affects stress is based on the theory of thermo-elasticity. Additionally, the coupling between stress/strain and fluid flow describing how pore-pressure affects stress is governed by Biot's consolidation theory. Stress/strain-dependent porosity/permeability/capillary-pressure defines the influence of stress on fluid flow and thermal fields. The THM model is implemented in the solver module that links the multiphase, multicomponent fluid flow simulator TOUGH2 with the geomechanical simulator RFPA (Rock Failure Process Analysis), which is based on the concept that heterogeneity leads to non-linearity and progressive failure behavior. Capability of TOUGH-RFPA simulator was verified by two benchmark problems related to hydraulic-mechanical and thermal-mechanical response in a hollow cylinder. Then the applicability of TOUGH-RFPA was demonstrated for modeling cooling-induced damage around a very deep geothermal well with mesoscopic heterogeneous properties, which showed complex damage evolution and fracturing during well cooling and fluid pressurization.

42 ENGINEERING↗

Biomass Combustion in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing resulting in higher combustion efficiency and the ability to use lower quality fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions. SO 2 emissions are conveniently controlled through the addition of Ca-based sulphur sorbents within the combustor. This report summarizes the current modeling effort on a 50 kW th CFB combustor with a diameter of 10 cm and a height of 5 m; designed, built, and operated at CanmetENERGY in Ottawa, Canada employing 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 9 kg of 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 the physics and chemistry of a circulating fluidized bed combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

20 FOSSIL-FUELED POWER PLANTS↗

Evidence for a multiphase, inhomogeneous interstellar medium in damped Ly alpha systems

The large optical depth in absorption of the Ly-alpha line in damped Ly-alpha systems implies that the emission line produced by H2 regions is very broad, since the number of diffusions is very large. The emission lines observed in DLy-alpha systems are narrow, a fact that can only be explained if the emission comes from the outer layers alone, or if the interstellar medium (ISM) is clumpy and multiphase. It is shown how discrete, cold H1 clouds embedded in a hot, fully ionized phase, produce narrow lines. This might provide evidence for a supernova remnant SNR-dominated ISM at high redshift, in agreement with the observed metallicity and chemical evolution models. Other high redshift galaxies also offer evidence for such a multiphase, in homogeneous ISM.

Valls-Gabaud, David↗

Active motion of multiphase oil droplets: emergent dynamics of squirmers with evolving internal structure

Synthetic soft matter systems, when driven beyond equilibrium by active processes, offer the potential to achieve dynamical states and functions of a complexity found in living matter. Emulsions offer the basis of a simple yet versatile system for identification of the physicochemical principles underlying active soft matter, but how multiple internal phases within emulsion droplets (e.g., Janus morphologies) organize to impact emergent dynamics is not understood. Here, we create multiphase oil droplets with ultralow interfacial tensions but distinct viscosities, and drive them into motion in aqueous micellar solutions. Preferential solubilization of select components of the oil both drives the droplet motion and yields a progression of internal phase morphological states with distinct symmetries. In this work, we find the active droplets to exhibit five dynamical states during morphogenesis. By quantifying microscopic flow fields, we show that it is possible to map the diverse droplet behaviors to squirmer models of spherical microswimmers in Stokes flow, thus showing that multiphase droplets offer the basis of a versatile platform with which to study and engineer the hydrodynamics of microswimmers.

36 MATERIALS SCIENCE↗

Risk-based area of review estimation in overpressured reservoirs to support injection well storage facility permit requirements for CO 2 storage projects

This paper by the Energy & Environmental Research Center presents a workflow and modeling approach for delineating a risk-based area of review (AOR) to support a U.S. Environmental Protection Agency (EPA) Class VI permit for a carbon dioxide (CO 2 ) storage project. The approach combines semianalytical solutions for estimating formation fluid leakage through a hypothetical leaky wellbore with the results of numerical reservoir simulations to define the AOR. The modeling utilizes 1) semianalytical solutions from the peer-reviewed literature for formation fluid leakage through abandoned wellbores by Raven (1990) and Avci (1994), 2) a FORTRAN model compiled and described in Cihan et al. (2011, 2012) called ASLMA (Analytical Solution for Leakage in Multilayered Aquifers), and 3) a computational framework for estimating a risk-based AOR first proposed by Oldenburg et al. (2014, 2016). Therefore, the approach builds upon well-established research and underlying hydrogeological principles that have been upheld for nearly three decades. Moreover, the ASLMA model has been broadly applied to an array of storage projects. The work presented herein extends these earlier works using a custom wrapper written in the software environment, R (R Core Team, 2020), which was developed to perform multiple runs of the ASLMA model using given ranges for one or more input parameters. In addition, the current work simulates the pressure buildup within the storage reservoir in response to CO 2 injection using a compositional simulator to better accommodate the temporospatial evolution of pressure buildup within the storage reservoir that is more accurately modeled using a heterogeneous geologic model and a compositional simulator that accounts for the multiphase interactions. The workflow is demonstrated using a case study for a 180,000-metric-ton-per-year storage project located in the PCOR (Plains CO 2 Reduction) Partnership region. For the storage project evaluated here, under the scenario where the leaky wellbore is open to a saline aquifer (thief zone) between the overlying seal (cap rock) and the underground sources of drinking water (USDW), the risk-based AOR essentially collapses to the areal extent of the CO 2 plume in the storage reservoir because the pressure buildup in the storage reservoir beyond the CO 2 plume is insufficient to drive formation fluids up a hypothetical leaky wellbore into the USDW. However, even under the conservative assumption that the leaky wellbore is not open to a thief zone, beyond the areal extent of the CO 2 plume, the incremental leakage is less than 400 m 3 over 20 years, which represents ~0.0001% or less of the total volume of water contained within the USDW rock volume. As discussed in the text, the threshold criterion for defining the risk-based AOR is site-specific and should be informed by the results of the sensitivity analysis and available site characterization data. The approach outlined in this paper is designed to be protective of USDWs and, therefore, comply with the Safe Drinking Water Act requirements and provisions for the U.S. EPA Class VI Underground Injection Control (UIC) Program (Class VI Rule) and North Dakota Administrative Code Chapter 43-05-01.

54 ENVIRONMENTAL SCIENCES↗

Tight Coupling of Surface and In-Plant Biochemistry and Convection Governs Key Fine Particulate Components over the Amazon Rainforest

We report that throughout the troposphere over the Amazon rainforest, isoprene epoxydiol-secondary organic aerosols (IEPOX-SOA) are key components of submicrometer particle mass. IEPOX-SOA is thought to form by multiphase chemical pathways in the atmosphere, but using detailed model simulations, we show that these pathways are strongly inhibited by the solid particle phase state prevailing in the upper troposphere. We find that 2-methyltetrol gases formed by in-plant biochemical oxidation pathways, their partitioning to fine particles, and dissolution in cloud droplets can explain over 90% of IEPOX-SOA mass concentrations in the upper troposphere. In agreement with previous PM10 chiral measurements, our simulations indicate that plant biochemical processes are the dominant source of particulate 2-methyltetrols in the lower troposphere as well (>70%). Our results imply that in-plant biochemistry plays a central but previously unidentified role in aerosol processes and atmosphere-biosphere-climate interactions over the rainforest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density and Cavitating Flow Results from a Full-Scale Optical Multiphase Cryogenic Flowmeter

Liquid propulsion systems are hampered by poor flow measurements. The measurement of flow directly impacts safe motor operations, performance parameters as well as providing feedback from ground testing and developmental work. NASA Marshall Space Flight Center, in an effort to improve propulsion sensor technology, has developed an all optical flow meter that directly measures the density of the fluid. The full-scale sensor was tested in a transient, multiphase liquid nitrogen fluid environment. Comparison with traditional density models shows excellent agreement with fluid density with an error of approximately 0.8%. Further evaluation shows the sensor is able to detect cavitation or bubbles in the flow stream and separate out their resulting effects in fluid density.

Korman, Valentin↗

High-Fidelity Generalization Method of Cells for Inelastic Periodic Multiphase Materials

An extension of a recently-developed linear thermoelastic theory for multiphase periodic materials is presented which admits inelastic behavior of the constituent phases. The extended theory is capable of accurately estimating both the effective inelastic response of a periodic multiphase composite and the local stress and strain fields in the individual phases. The model is presently limited to materials characterized by constituent phases that are continuous in one direction, but arbitrarily distributed within the repeating unit cell which characterizes the material's periodic microstructure. The model's analytical framework is based on the homogenization technique for periodic media, but the method of solution for the local displacement and stress fields borrows concepts previously employed by the authors in constructing the higher-order theory for functionally graded materials, in contrast with the standard finite-element solution method typically used in conjunction with the homogenization technique. The present approach produces a closed-form macroscopic constitutive equation for a periodic multiphase material valid for both uniaxial and multiaxial loading. The model's predictive accuracy in generating both the effective inelastic stress-strain response and the local stress said inelastic strain fields is demonstrated by comparison with the results of an analytical inelastic solution for the axisymmetric and axial shear response of a unidirectional composite based on the concentric cylinder model, and with finite-element results for transverse loading.

Aboudi, Jacob↗

Coupled Interfacial and Bulk Kinetics Govern the Timescales of Multiphase Ozonolysis Reactions

Chemical transformations in aerosols impact the lifetime of particle phase species, the fate of atmospheric pollutants, and both climate- and health-relevant aerosol properties. Timescales for multiphase reactions of ozone in atmospheric aqueous phases are governed by coupled kinetic processes between the gas phase, the particle interface, and its bulk, which respond dynamically to reactive consumption of O 3 . However, models of atmospheric aerosol reactivity often do not account for the coupled nature of multiphase processes. Herein, to examine these dynamics, we use new and prior experimental observations of aqueous droplet reaction kinetics, including three systems with a range of surface affinities and ozonolysis rate coefficients ( trans -aconitic acid (C 6 H 6 O 6 ), maleic acid (C 4 H 4 O 4 ), and sodium nitrite (NaNO 2 )). Using literature rate coefficients and thermodynamic properties, we constrain a simple two-compartment stochastic kinetic model which resolves the interface from the particle bulk and represents O 3 partitioning, diffusion, and reaction as a coupled kinetic system. Our kinetic model accurately predicts decay kinetics across all three systems, demonstrating that both the thermodynamic properties of O 3 and the coupled kinetic and diffusion processes are key to making accurate predictions. An enhanced concentration of adsorbed O 3 , compared to gas and bulk phases is rapidly maintained and remains constant even as O 3 is consumed by reaction. Multiphase systems dynamically seek to achieve equilibrium in response to reactive O 3 loss, but this is hampered at solute concentrations relevant to aqueous aerosol by the rate of O 3 arrival in the bulk by diffusion. As a result, bulk-phase O 3 becomes depleted from its Henry's law solubility. This bulk-phase O 3 depletion limits reaction timescales for relatively slow-reacting organic solutes with low interfacial affinity (i.e., trans-aconitic and maleic acids, with k rxn ≈ 10 3 -10 4 M -1 s -1 ), which is in contrast to fast-reacting solutes with higher surface affinity (i.e., nitrite, with k rxn ≈ 10 5 M -1 s -1 ) where surface reactions strongly impact the observed decay kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗