Search NASA⌕ Search

SEARCH · Search NASA

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Review—Drag Coefficients of Non-Spherical and Irregularly Shaped Particles

Abstract The knowledge of simple and relatively accurate closure equations for the drag coefficients of nonspherical particles is very important for Eulerian multiphase numerical codes that require such inputs for the solution of the momentum equation. However, the existing correlations for the drag coefficients are laden with high uncertainty, which propagates in the numerical results. This review paper examines critically the state of our knowledge of the drag coefficients of nonspherical particles starting with our understanding of the “size” and “shape” of particles. Lengthscales and other parameters that have been used to characterize sizes of nonspherical particles are presented. Shape factors and shape descriptors that have been used in past correlations are also presented together with an appraisal of their applicability. Several correlations that have been developed for the determination of the drag coefficients of nonspherical particles with regular and irregular shapes are also presented together with available information on their accuracy and applicability. Based on the review of the available correlations general recommendations are offered on the use of the correlations as closure equations and inputs in computatioal fluid dynamics (CFD) codes.

Engineering↗
DOE OSTI2024

PeleMP: The Multiphysics Solver for the Combustion Pele Adaptive Mesh Refinement Code Suite

Combustion encompasses multiscale, multiphase reacting flow physics spanning a wide range of scales from the molecular scales, where chemical reactions occur, to the device scales, where the turbulent flow is affected by the geometry of the combustor. This scale disparity and the limited measurement capabilities from experiments make modeling combustion a significant challenge. Recent advancements in high-performance computing (HPC), particularly with the Department of Energy's Exascale Computing Project (ECP), have enabled high-fidelity simulations of practical applications to be performed. The major physics submodels, including chemical reactions, turbulence, sprays, soot, and thermal radiation, exhibit distinctive computational characteristics that need to be examined separately to ensure efficient utilization of computational resources. This paper presents the multiphysics solver for the Pele code suite, called PeleMP, which consists of models for spray, soot, and thermal radiation. Here, the mathematical and algorithmic aspects of the model implementations are described in detail as well as the verification process. The computational performance of these models is benchmarked on multiple supercomputers, including Frontier, an exascale machine. Results are presented from production simulations of a turbulent sooting ethylene flame and a bluff-body swirl stabilized spray flame with sustainable aviation fuels to demonstrate the capability of the Pele codes for modeling practical combustion problems with multiphysics. This work is an important step toward the exascale computing era for high-fidelity combustion simulations providing physical insights and data for predictive modeling of real-world devices.

42 ENGINEERING↗
DOE OSTI2021

An Investigation of the Effects of Volume Fraction on Drag Coefficient of Non-Spherical Particles Using PR-DNS

Prediction of the drag coefficient is required in gas-particle multiphase flow modeling and simulation. Experimental data and correlations on the fixed-bed system of spherical particles with high volume fractions for various possible arrangements are available in the literature. However, the effect of volume fraction on the drag coefficient of non-spherical particles is not well studied. In solving the momentum equation, the volume fraction plays a vital role in determining the flow resistances. In this paper, we study the impact of volume fraction in the range of 0.069 to 0.65 on the drag coefficient using the computational fluid dynamics (CFD) simulation of air for Reynold number in the range of 10 to 10000 using particle resolved direct numerical solution (PR-DNS). Regular non-spherical particles such as a cube, tetrahedron, and spheroids are used in this study since their single particle’s drag coefficient data are available in the literature for comparison. For this work, the simulations are carried out in the Ansys Fluent using polyhedral mesh, which consumes significantly less computational time and power. The study showed the sphericity and volume fraction have significant impact on the bed pressure drop and average drag coefficient of the particles in the bed especially in high Reynolds number regime. The bed of the spheroid experiences the lowest drag being the most streamlined particle, and the particles with the edges result in a large drag coefficient due to flow separation at the discontinuity. The vector plots verify this behavior where large wake regions are observed behind the tetrahedron particle.

Mahyawansi, Pratik↗
DOE OSTI2023

Numerical Study of Liquid Piston Compression Using Large-Eddy Simulation and Volume-of-Fluid Approach

Efforts to increase the efficiency of residential and commercial air conditioners and heat pumps have demonstrated that the compressor accounts for most of the system’s electrical energy usage. Therefore, the efficiency of this component should be improved to reduce its energy usage. The US Department of Energy’s Oak Ridge National Laboratory developed a near-isothermal liquid piston compressor (LPC) that uses propylene glycol (PG) to compress CO2. This report presents numerical studies of the LPC in which the compression chamber fills with injected PG from the bottom inlet. Numerical simulations were performed using the large-eddy simulation (LES) with the wall-adapting local eddy-viscosity (WALE) subgrid-scale model coupled with the multiphase volume-of-fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfer within the compression chamber. In this effort, effects of boundary conditions applied to the LES-VOF calculations (e.g., no wall, an adiabatic wall, and a wall with a heat flux subscribed) to the overall pressure and temperature of CO2 gas as well as the transient evolution of flow and heat transfer evolution within the compression chamber are investigated and discussed. It was found that the LES calculation with no wall have shown no dynamical flow patterns and the volume-averaged temperature of CO2 increased from 305 to 392.7 K, while LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO2 temperatures. Results of LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO2 gas layers and circulation flow patterns within the chamber that contribute to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor.

Nguyen, Thien D.↗
DOE OSTI2024

Mesofluidic Oil-Water Separation

Separating flowing multiphase solutions in process at industrial scale is challenging. A rapid continuous approach is desired. These mesofluidic oil-water separation tests were motivated by reports by Dijkshoorn et al [1] who reported separating droplets and particles using displacement with regularly spaced channels using a cross-flow microsieve.. Droplets, including deformable particles, larger than the channel migrate in one direction and thereby deplete in the other direction. Test flowrates were motivated by Burns et al. [4] who showed that separation could occur under laminar conditions, and Pease et al [8, 9] who extended the range to turbulent flow. Testing was conducted using two two-component oil-water mixtures at room temperature. Analysis showed separation of both mineral oils from their oil water mixtures at flowrates from 1-2 L/min and concentrations from 2-10 wt% oil in water.

cloud height, effective cleaning radius, effective↗
DOE OSTI2018

Review Article: Stress in thin films and coatings: Current status, challenges, and prospects

The issue of stress in thin films and functional coatings is a persistent problem in materials science and technology that has congregated many efforts, both from experimental and fundamental points of view, to get a better understanding on how to deal with, how to tailor, and how to manage stress in many areas of applications. With the miniaturization of device components, the quest for increasingly complex film architectures and multiphase systems and the continuous demands for enhanced performance, there is a need toward the reliable assessment of stress on a submicron scale from spatially resolved techniques. Also, the stress evolution during film and coating synthesis using physical vapor deposition (PVD), chemical vapor deposition, plasma enhanced chemical vapor deposition (PECVD), and related processes is the result of many interrelated factors and competing stress sources so that the task to provide a unified picture and a comprehensive model from the vast amount of stress data remains very challenging. This article summarizes the recent advances, challenges, and prospects of both fundamental and applied aspects of stress in thin films and engineering coatings and systems, based on recent achievements presented during the 2016 Stress Workshop entitled “Stress Evolution in Thin Films and Coatings: from Fundamental Understanding to Control.” Evaluation methods, implying wafer curvature, x-ray diffraction, or focused ion beam removal techniques, are reviewed. Selected examples of stress evolution in elemental and alloyed systems, graded layers, and multilayer-stacks as well as amorphous films deposited using a variety of PVD and PECVD techniques are highlighted. Based on mechanisms uncovered by in situ and real-time diagnostics, a kinetic model is outlined that is capable of reproducing the dependence of intrinsic (growth) stress on the grain size, growth rate, and deposited energy. The problems and solutions related to stress in the context of optical coatings, inorganic coatings on plastic substrates, and tribological coatings for aerospace applications are critically examined. This review also suggests strategies to mitigate excessive stress levels from novel coating synthesis perspectives to microstructural design approaches, including the ability to empower crack-based fabrication processes, pathways leading to stress relaxation and compensation, as well as management of the film and coating growth conditions with respect to energetic ion bombardment. Future opportunities and challenges for stress engineering and stress modeling are considered and outlined.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase composition of polycrystalline HfNx (0.45 ≤ x ≤ 1.60) and effects of low-energy ion irradiation on microstructure, texture, and physical properties

We have investigated the phase composition of HfNx as a function of x and the effects of low-energy ion irradiation on the microstructure and physical properties of polycrystalline layers grown on SiO2 at 350 °C by ultrahigh vacuum reactive dc magnetron sputtering of Hf in mixed N2/Ar discharges. X-ray diffraction and Rutherford backscattering spectrometry results show that the phases obtained in polycrystalline HfNx layers with increasing x are hcp-structure α-Hf:N (x ≲ 0.6); multiphase mixtures consisting of α-Hf, NaCl-structure δ-HfN, rhombohedral ɛ-Hf3N2, and/or ζ-Hf4N3 (0.6 ≲ x ≲ 0.9); δ-HfN single phase (0.9 ≲ x ≲ 1.3); and mixtures of δ-HfN and higher nitrides (x ≳ 1.3). HfNx layers with 0.9 ≲ x ≲ 1.2 grown under mild ion irradiation (incident ion energy Ei ≃ 7 eV and ion-to-Hf flux ratios Ji/JHf = 1−3) are underdense with mixed orientation, low in-plane stress, and rough surface morphology due to limited adatom mobilities resulting in kinetic roughening and atomic shadowing during film growth. However, the use of intense ion irradiation (Ei = 25 eV and Ji/JHf = 4−20) results in HfNx layers, which are fully dense with strongly 111-oriented texture, compressive in-plane stress, and smooth surfaces due to ion irradiation enhanced adatom surface mobilities. In addition, the latter films have lower resistivity and higher hardness. For stoichiometric δ-HfN layers, ρ decreases from 69.7 to 35.2 μΩ cm and H increases from 22.1 to 27.4 GPa, with increasing ion-irradiation intensity. However, for HfNx layers with 1.2 ≲ x ≲ 1.6, the correspondingly higher steady state atomic N surface coverages during deposition alter growth kinetics in favor of 001 texture with a fully dense structure and compressive in-plane stress.

Materials Science↗
DOE OSTI2021

Design of tough adhesive from commodity thermoplastics through dynamic crosslinking

Tough adhesives provide resistance against high debonding forces, and these adhesives are difficult to design because of the simultaneous requirement of strength and ductility. Here, we report a design of tough reversible/recyclable adhesive materials enabled by incorporating dynamic covalent bonds of boronic ester into commodity triblock thermoplastic elastomers that reversibly bind with various fillers and substrates. The spectroscopic measurements and density functional theory calculations unveil versatile dynamic covalent binding of boronic ester with various hydroxy-terminated surfaces such as silica nanoparticles, aluminum, steel, and glass. The designed multiphase material exhibits exceptionally high adhesion strength and work of debonding with a rebonding capability, as well as outstanding mechanical, thermal, and chemical resistance properties. Bonding and debonding at the interfaces dictate hybrid material properties, and this revelation of tailored dynamic interactions with multiple interfaces will open up a new design of adhesives and hybrid materials.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Giant piezoelectricity in oxide thin films with nanopillar structure

High-performance piezoelectric materials are critical components for electromechanical sensors and actuators. For more than 60 years, the main strategy for obtaining large piezoelectric response has been to construct multiphase boundaries, where nanoscale domains with local structural and polar heterogeneity are formed, by tuning complex chemical compositions. We used a different strategy to emulate such local heterogeneity by forming nanopillar regions in perovskite oxide thin films. We obtained a giant effective piezoelectric coefficient d 33 , f * of ~1098 picometers per volt with a high Curie temperature of ~450°C. Our lead-free composition of sodium-deficient sodium niobate contains only three elements (Na, Nb, and O). The formation of local heterogeneity with nanopillars in the perovskite structure could be the basis for a general approach to designing and optimizing various functional materials.

Science & Technology - Other Topics↗
DOE OSTI2019

Synthetic Methane-Consuming Communities from a Natural Lake Sediment

The factors and processes that influence the behavior and functionality of ecosystems inhabited by complex microbiomes are still far from being clearly understood. Synthetic microbial communities provide reduced-complexity models that allow an examination of ecological theories under defined and controlled conditions. In this study, we applied a multiphasic approach to study synthetic methane-oxidizing communities and species interactions as proxies to the natural communities. Our results confirm that, under selective pressures, natural-sediment communities of high complexity simplify rapidly, selecting for several major functional guilds, the major partners in methane oxidation being the Methylococcaceae methanotrophs and the Methylophilaceae methylotrophs, along with minor but persistent partners, members of Burkholderiales and Flavobacteriales. As a proof of concept, we established minimalist synthetic communities that were representative of the four functional guilds to demonstrate the dependency of the non-methane-utilizing species on the methanotrophs as the primary carbon-providing species. We observed that in communities consisting of multiple representatives of the key guilds, members of the same guild appeared to compete for resources. For example, when two methanotrophs of the same family were present, the two expressed similar key methanotrophy pathways and responded similarly to changing environmental conditions, suggesting that they perform a similar keystone function in situ. Similar observations were made for the Methylophilaceae. However, differences were noted in the expression of auxiliary and unique genes among strains of the same functional guild, reflecting differential adaptation and suggesting mechanisms for competition. At the same time, differences were also noted in the performances of partners with specific metabolic schemes. For example, a mutant of Methylotenera mobilis impaired in nitrate utilization behaved as a more efficient cooperator in methane consumption, suggesting that the loss of function may lead to changes in communal behavior. Overall, we demonstrate the robust nature of synthetic communities built of native lake sediment strains and their utility in addressing important ecological questions while using a simplified model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Evenly spaced columns in the Bishop Tuff (California, USA) as relicts of hydrothermal cooling

A few square kilometers of the Bishop Tuff in eastern California (USA) have evenly spaced columns that are more resistant to erosion than the surrounding tuff owing to the precipitation of mordenite, a low-temperature (100-130 °C) zeolite. We hypothesize that the columns are a result of instabilities at the liquid water and steam interface as cold water seeped into the still-cooling Bishop Tuff. We use two methods to quantitatively assess this hypothesis. First, scaling shows which hydrodynamic instabilities exist in the system. Second, to account for the effects of multiphase flow, latent heat, and the finite amplitude and temporal evolution of these instabilities we use two-dimensional numerical models of liquid water infiltrating hot tuff. These tests highlight several features of boiling hydrothermal systems. (1) The geometry of at least some convection appears to be broadly captured by linear stability theory that neglects reactive transport, heterogeneity of the host rock, and the finite amplitude of instabilities. (2) Slopes > 10% set the wavelength of convection, meaning that these columns formed somewhere with relatively gentle topography. (3) For permeabilities of > 10 -13 m 2 , the wavelength of the instability changes through time, slowing infiltration, while for permeabilities < 10 -15 m 2 , cooling is dominated by conduction. The spacing and stability of columns increase with higher vertical permeability and decrease with higher horizontal permeability. These columns are a rare window into hydrothermal processes that may be widespread.

58 GEOSCIENCES↗
DOE OSTI2019

Effect of Ion (Ar{sup +}) Irradiation on Cluster Magnetism and Magnetic Interactions in Fe{sub 67}Cr{sub 18}B{sub 15} Amorphous Alloy

The comparison of the resonance magnetic absorption spectra with the data on structure and magnetic properties of Fe{sub 67}Cr{sub 18}B{sub 15} amorphous samples has shown that a randomly inhomogeneous medium of amorphous alloy is multiphase. Two types of samples with the same chemical composition, but with a different clusterized structure, and the use of (Ar{sup +}) ion irradiation in the experiment made it possible to separate four basic phases with different magnetic properties. The basic phase α-(Fe,Cr) contains large clusters of D < 400 Å in size with a tendency to phase separation of iron (α-Fe) and chromium (α-Cr), which under irradiation, enhances the competing FM and AFM magnetic interactions and increases the intensity of the resonance magnetic absorption lines associated with α-Fe and α-Cr. The phase Fe{sub 3}B (D < 100 Å) is detected by independent line in the resonance absorption spectra only when its groupings are isolated from the α-(Fe,Cr) metal phase clusters. The fourth phase contains precluster groupings and disappears when the sample is irradiated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Efficient Operator-Coarsening Multigrid Schemes for Local Discontinuous Galerkin Methods

An efficient hp-multigrid scheme is presented for local discontinuous Galerkin (LDG) discretizations of elliptic problems, formulated around the idea of separately coarsening the underlying discrete gradient and divergence operators. We show that traditional multigrid coarsening of the primal formulation leads to poor and suboptimal multigrid performance, whereas coarsening of the flux formulation leads to essentially optimal convergence and is equivalent to a purely geometric multigrid method. The resulting operator-coarsening schemes do not require the entire mesh hierarchy to be explicitly built, thereby obviating the need to compute quadrature rules, lifting operators, and other mesh-related quantities on coarse meshes. We show that good multigrid convergence rates are achieved in a variety of numerical tests on two-dimensional and three-dimensional uniform and adaptive Cartesian grids, as well as for curved domains using implicitly defined meshes and for multiphase elliptic interface problems with complex geometry. Extensions, e.g., to non-LDG discretizations and fully unstructured meshes, are briefly discussed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mechanistic and Mitigation-Strategy Insights into NaCl and CaCl 2 Contamination of Proton-Exchange-Membrane Water Electrolysis Using Continuum Modeling

Cationic contaminants are detrimental to proton-exchange-membrane water electrolyzers (PEMWEs). To obtain insight, a 1-D, nonisothermal, multiphase continuum cell model including cationic contamination is developed. Simulations of steady-state cell performance predict decreased performance due to an increase in kinetic overpotential associated with the hydrogen-evolution reaction, which was attributed to decreased protonic-activity within the cathode catalyst layer from proton supplantation with contaminant cations. The accumulation and extent of cation exchange in the cathode catalyst layer depends on the operating current density due to migration. Simulations of cell recovery of potential suggest that a contaminated cell can recover approximately 78% (450 mV) with 24 h of constant current density operation at 2 A cm –2 , with higher current densities accelerating reduced recovery times. Parametric studies show that anode-side acidification at lower current densities inhibit cation contaminant adsorption, and cathode-side acidification at larger current densities facilitate the expulsion of adsorbed cations; for a cathode-side pH of 6 and 5, the cell can recover an additional 10% and 100% performance, respectively. Overall, the model serves as a framework for modeling other aspects of PEMWE systems to address durability and performance aspects, which can assist in improving the viability of the technology.

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

Elucidating Performance Limitations in Alkaline-Exchange- Membrane Fuel Cells

Water management is a serious concern for alkaline-exchange-membrane fuel cells (AEMFCs) because water is a reactant in the alkaline oxygen-reduction reaction and hydroxide conduction in alkaline-exchange membranes is highly hydration dependent. Here in this article, we develop and use a multiphysics, multiphase model to explore water management in AEMFCs. We demonstrate that the low performance is mostly caused by extremely non-uniform distribution of water in the ionomer phase. A sensitivity analysis of design parameters including humidification strategies, membrane properties, and water transport resistance was undertaken to explore possible optimization strategies. Furthermore, the strategy and issues of reducing bicarbonate/carbonate buildup in the membrane-electrode assembly with CO 2 from air is demonstrated based on the model prediction. Overall, mathematical modeling is used to explore trends and strategies to overcome performance bottlenecks and help enable AEMFC commercialization.

25 ENERGY STORAGE↗
DOE OSTI2020

(Invited) Virtual Analysis of Gas-Diffusion-Electrode CO 2 Electrolyzers

The electrochemical reduction of CO 2 (CO 2 R) to value-added products is an attractive technology for tackling the rising atmospheric CO 2 levels and storing intermittent renewable energy into chemical bonds. Fundamental understanding of CO 2 R has progressed significantly in recent years and is critical in the development of CO 2 R to liquid-fuel electrolyzers, where gas-diffusion electrodes (GDEs) have been shown to be key enabling architectures. Various designs have been proposed and studied in the literature to enhance overall selectivity, rates, and maximize the conversion of CO 2 , the latter of which is only now being recognized as a critical issue. In this respect, there is a need to explore the governing phenomena inherent in these architectures to enable optimization. Mathematical modeling is ideally suited to tackle and explore these multiphysics interactions and provide virtual design analysis. In this talk, we discuss modeling methodologies and physics inherent in these devices and present our recent modeling of GDEs for CO 2 reduction. We specifically examine the impacts of multiphase flow and related phenomena on overall cell performance. We then explore the performance and limitations of various cell designs guided by simulation results and examine potential methods for improving water management and tuning catalyst selectivity including the use of different anion-exchange and bipolar membranes. Finally, we discuss the disparities in local environments between aqueous and GDE devices and propose strategies to reduce the gap in knowledge between the two systems.

42 ENGINEERING↗
DOE OSTI2020

Virtual Analysis of Gas-Diffusion-Electrode CO 2 Electrolyzers

The electrochemical reduction of CO 2 (CO2R) to value-added products is an attractive technology for tackling the rising atmospheric CO 2 levels and storing intermittent renewable energy into chemical bonds. Fundamental understanding of CO 2 R has progressed significantly in recent years and is critical in the development of CO 2 R to liquid-fuel electrolyzers, where gas-diffusion electrodes (GDEs) have been shown to be key enabling architectures. Various designs have been proposed and studied in the literature to enhance overall selectivity, rates, and maximize the conversion of CO 2 , the latter of which is only now being recognized as a critical issue. In this respect, there is a need to explore the governing phenomena inherent in these architectures to enable optimization. Mathematical modeling is ideally suited to tackle and explore these multiphysics interactions and provide virtual design analysis. In this talk, we discuss modeling methodologies and physics inherent in these devices and present our recent modeling of GDEs for CO 2 reduction. Here, we specifically examine the impacts of multiphase flow and related phenomena on overall cell performance. We then explore the performance and limitations of various cell designs guided by simulation results and examine potential methods for improving water management and tuning catalyst selectivity including the use of different anion-exchange and bipolar membranes. Finally, we discuss the disparities in local environments between aqueous and GDE devices and propose strategies to reduce the gap in knowledge between the two systems.

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

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗