Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase flow”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Multiphase Modeling in a Parallel Plate Fluidized Bed Receiver for Concentrating Solar Power

A novel high temperature particle solar receiver is developed by using a light trapping planar cavity configuration. As particles fall through the cavity, the concentrated solar radiation warms the boundaries of the receiver and in turn heats the particles. Particles flow through the system, forming a packed bed at the lower end, leaving the system from the bottom at a constant flow rate. Air is introduced to the system as the fluidizing medium to improve particle heat transfer and mixing. A laboratory scale cavity receiver is built and a near IR quartz lamp is used to provide flux to the vertical wall of the heat exchanger. The system is modeled using a continuum two-fluid method. The computational model matches the experimental system size and the particle size distribution is assumed monodisperse. A conduction model that accounts for the effects of solid concentration is implemented, and the heat flux boundary condition matches the experimental setup. Radiative heat transfer is estimated using a widely used correlation during the post-processing step to determine an overall heat transfer coefficient. The model is validated against testing data and achieves less than 30% discrepancy and a heat transfer coefficient greater than 1000 W/m2K.

CSP↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1984

Purification and cultivation of human pituitary growth hormone secreting cells

A multiphase study was conducted to examine the properties of growth hormone cells. Topics investigated included: (1) to determine if growth hormone (GH) cells contained within the rat pituitary gland can be separated from the other hormone producing cell types by continuous flow electrophoresis (CFE); (2) to determine what role, if any, gravity plays in the electrophoretic separation of GH cells; (3) to compare in vitro GH release from rat pituitary cells previously exposed to microgravity conditions vs release from cells not exposed to microgravity; (4) to determine if the frequency of different hormone producing pituitary cell types contained in cell suspensions can be quantitated by flow cytometry; and (5) to determine if GH contained within the human post mortem pituitary gland can be purified by CFE. Specific experimental procedures and results are included.

Hymer, W. C.↗
DOE OSTI2026

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗
DOE OSTI2023

Brewster angle-cavity ringdown spectroscopy for low temperature plasma measurements in multiphases

Here we report on the development of a Brewster angle-cavity ringdown spectroscopy (BA-CRDS) system for low temperature plasma diagnostics. The system can measure gas species in solutions, with a detection limit (minimum detectable absorbance) of 9.1 × 10 -5 , which is equivalent to a detection limit of 0.04 parts per billion for measuring OH radicals in water at 308 nm. With higher reflectivity ringdown mirrors and improved design of a Brewster angle cell, the detection limit can potentially be up to 10 -6 or lower. In this exploratory study, the absorption cross sections of HgBr 2 and H 2 O 2 in the aqueous phase at 256 nm are measured to be (1.8 ± 0.1) × 10 -18 cm 2 and (5.2 ± 0.5) × 10 -20 cm 2 , respectively. Furthermore, temporal profiles of absorbance from distilled water, HgBr 2 , and H 2 O 2 solutions when interacting with a helium atmospheric plasma jet are individually characterized at different plasma powers, gas flow rates, and/or solute concentrations. The observed linear temporal profiles of absorbance from the plasma-interacted water suggest formation of H 2 O 2 from plasma-generated OH radicals, while the nonlinear temporal profiles from the plasma-treated HgBr 2 solutions reveal possible removal of HgBr 2 by OH radicals. Our results demonstrate that the new BA-CRDS system is a powerful tool for quantification of reactive plasma species in multiphases or other complex settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1979

The analysis and design of transonic two-element airfoil systems

The multiphase effort in the development of tools for the analysis and design of two-element airfoil systems, that is, airfoils with a slat or a flap at transonic speeds is described. The first phase involved the development of a method to compute the inviscid flow over such configurations. In the second phase the inviscid code was coupled to a boundary layer calculation program in order to compute the loss in performance due to viscous effects. An inverse code that constructs the airfoil system corresponding to a desired pressure distribution is described.

Volpe, G.↗
NASA NTRS1986

Numerical study of combustion processes in afterburners

Mathematical models and numerical methods are presented for computer modeling of aeroengine afterburners. A computer code GEMCHIP is described briefly. The algorithms SIMPLER, for gas flow predictions, and DROPLET, for droplet flow calculations, are incorporated in this code. The block correction technique is adopted to facilitate convergence. The method of handling irregular shapes of combustors and flameholders is described. The predicted results for a low-bypass-ratio turbofan afterburner in the cases of gaseous combustion and multiphase spray combustion are provided and analyzed, and engineering guides for afterburner optimization are presented.

Zhou, Xiaoqing↗
DOE OSTI2020

Progress Toward Simulating Departure from Nucleate Boiling at High-Pressure Applications with Selected Wall Boiling Closures

Recently, a Eulerian-based two-fluid computational fluid dynamics (CFD) framework with a wall heat flux partitioning approach has been intensively investigated for departure from nucleate boiling (DNB) simulation under the U.S. Department of Energy–funded Consortium for Advanced Simulation of Light Water Reactors (CASL) program. Understanding of the DNB characteristics over a range of pressurized water reactor–like operating conditions and accurate prediction of boiling crisis in the nuclear power system have been grand challenges because of the large impact of DNB on reactor safety and operational economics. The ultimate goal of this task in the CASL program is to introduce a robust multiphase CFD–based DNB modeling framework that is capable of characterizing an entire boiling history in which the wall boiling mode experiences the following through multiple stages of heat transfer mode: (1) single-phase convective heat transfer, (2) nucleate boiling heat transfer, and (3) identification of the departure of nucleate boiling. To validate the CASL boiling model, we have benchmarked simulated DNB over three different flow channel configurations (pipe flow, 5 × 5 fuel bundle with mixing vane tests, and 5 × 5 fuel bundle without mixing vane tests) against experimental measurements, and the validation result with open literature is reported. The DNB detection criteria in the simulation are checked by monitoring the peak wall temperature, wall dryout factor, and net energy balance. In addition to the DNB performance test, some preliminary sensitivity results on closure model selection are reported to address the prediction capability of local void profile against measurements. The boiling simulation tested in this study exhibits a maximum deviation of 24% from the measured DNB value in a high-pressure (i.e., 138 bars) subcooled pipe flow test. The ranges of operating conditions are as follows: 1650 to 2650 kg/m 2 ·s for mass flux and 8.5 to 96 K for subcooled inlet temperature. The deviation is even reduced to 7% when the subcooled temperature is less than 40 K. Besides accuracy, base practice guidelines for DNB detection criteria are tested by monitoring three simulation variables: (1) maximum wall temperature, (2) wall dryout factor (i.e., K-value), and (3) energy balance. Numerical robustness of DNB simulation is largely achieved in most of the validation test except for a few high subcooled test cases.

42 ENGINEERING↗
DOE OSTI2021

Diffusion-to-Imbibition Transition in Water Sorption in Nanoporous Media: Theoretical Studies

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. Here, the model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rh < rh crit and the imbibition-like flow regimes observed at rh > rh crit . The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

58 GEOSCIENCES↗
NASA NTRS2011

Systems and methods for separating a multiphase fluid

Apparatus and methods for separating a fluid are provided. The apparatus can include a separator and a collector having an internal volume defined at least in part by one or more surfaces narrowing toward a bottom portion of the volume. The separator can include an exit port oriented toward the bottom portion of the volume. The internal volume can receive a fluid expelled from the separator into a flow path in the collector and the flow path can include at least two directional transitions within the collector.

Weislogel, Mark M.↗
DOE OSTI2026

Model‐Based Interpretation of Solute Exports and Carbon Partitioning During Shale Weathering in a Mountainous Hillslope

The weathering of sedimentary rocks in high-elevation catchments influences freshwater quality and the global carbon cycle. While individual biogeochemical mechanisms involved in this process are relatively well understood, quantifying their contributions to solute export and carbon fluxes under natural, transient conditions remains challenging. Here, we implement a numerical multidimensional and multiphase model to simulate coupled hydrological and biogeochemical processes in a shale-underlain, snow-dominated hillslope in the Rocky Mountains, Colorado. The model captures the dynamic interplay between soil respiration, mineral weathering, and climate-driven hydrological forcing, reproducing observed soil CO 2 dynamics, groundwater chemistry, and subsurface flow. Our results reveal that seasonal snowmelt enhances carbonate weathering by promoting the infiltration of CO 2 -rich water to depth, while pyrite oxidation is primarily sensitive to low water saturation that facilitates O 2 diffusion through the regolith. Topography modulates the spatial distribution of shale weathering, as steeper slopes enhance lateral drainage, favoring the delivery of reactants to greater depths. While shale weathering at our site acts as a transient carbon sink, with silicates and carbonates buffering acidity and promoting atmospheric CO 2 consumption (1% of soil-derived CO 2 ), the exported dissolved inorganic carbon is predominantly geogenic (∼73%). Consequently, when accounting for long-term marine carbonate precipitation. The current weathering regime represents a net source of carbon to the atmosphere. The oxidation of pyrite and petrogenic organic carbon together release approximately 0.9 mol·m −2 ·yr −1 of CO 2 . Our findings highlight the role of topography, hydroclimate, and the coupling between acid-base reactions in shaping the carbon balance and the solute exports in mountainous critical zones.

carbon cycling↗
DOE OSTI2020

Heterogeneous oxidation of amorphous organic aerosol surrogates by O 3 , NO 3 , and OH at typical tropospheric temperatures

Typical tropospheric temperatures render possible phase states of amorphous organic aerosol (OA) particles of solid, semisolid, and liquid. This will affect the multiphase oxidation kinetics involving the organic condensed-phase and gaseous oxidants and radicals. To quantify this effect, we determined the reactive uptake coefficients (γ) of O 3 , NO 3 , and OH by substrate films composed of single and binary OA surrogate species under dry conditions for temperatures from 213 to 313 K. A temperature-controlled coated-wall flow reactor coupled to a chemical ionization mass spectrometer was applied to determine γ with consideration of gas diffusion transport limitation and gas flow entrance effects, which can impact heterogeneous reaction kinetics. The phase state of the organic substrates was probed via the poke-flow technique, allowing the estimation of the substrates' glass transition temperatures. γ values for O 3 and OH uptake to a canola oil substrate, NO 3 uptake to a levoglucosan and a levoglucosan / xylitol substrate, and OH uptake to a glucose and glucose / 1,2,6-hexanetriol substrate have been determined as a function of temperature. We observed the greatest changes in γ with temperature for substrates that experienced the largest changes in viscosity as a result of a solid-to-liquid phase transition. Organic substrates that maintain a semisolid or solid phase state and as such a relatively higher viscosity do not display large variations in heterogeneous reactivity. From 213 to 293 K, γ values of O 3 with canola oil, of NO 3 with a levoglucosan / xylitol mixture, and of OH with a glucose / 1,2,6-hexanetriol mixture and canola oil, increase by about a factor of 34, 3, 2, and 5, respectively, due to a solid-to-liquid phase transition of the substrate. These results demonstrate that the surface and bulk lifetime of the OA surrogate species can significantly increase due to the slowed heterogeneous kinetics when OA species are solid or highly viscous in the middle and upper troposphere. This experimental study will further our understanding of the chemical evolution of OA particles with subsequent important consequences for source apportionment, air quality, and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Enhanced energy conversion efficiency promoted by cavitation in gasoline direct injection

High-pressure direct fuel injection plays the most crucial role in energy conversion and improving engine combustion efficiency and emission. The optimization of turbulent and multiphase fuel injection has focused on controlling hydrodynamic parameters such as injection pressure. While the thermodynamic influence is often considered in the flash boiling situation, we inquire into how gasoline-type fuel’s hydro- and thermodynamic properties impact the injection dynamics by fuel-temperature induced cavitation. The turbulent and cavitating flows emanating from the direct-injection nozzle is visualized by ultrafast x-ray imaging with an unprecedented spatiotemporal resolution. The ultrafast liquid-fuel dynamics are dominated by injection pressure as well as fuel temperature through cavitation, an important thermodynamic parameter but often difficult to control in engine combustion. With the most direct and quantitative measurement, we discovered that the near-nozzle fuel-jet dynamics can be perfectly scaled by a single dimensionless parameter, cavitation number, particularly sensitive to the fuel temperature, in a wide operation range. This universal scaling shows that cavitation can be harnessed to elevate the pneumatic-hydraulic to kinetic energy conversion efficiency, critical for promoting fuel atomization and engine combustion performance. This enhancement effect will have even more impact on engine combustion using alternative low-emission fuels with higher saturated vapor pressure.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Mathematical modeling of novel porous transport layer architectures for proton exchange membrane electrolysis cells

Thin foil based porous transport layers (PTLs) that contain highly structured pore arrays have shown promise as anode PTLs in proton exchange membrane electrolysis cells. These novel PTLs, fabricated with advanced manufacturing techniques, produce thin, tunable, multifunctional layers with reduced flow and interfacial resistances and high thermal and electric conductivities. To further optimize their design, it is important to understand their fundamental impact on the transport of protons, electrons, and liquid/vapor mixtures in the electrode. In this work, we develop a two-dimensional multiphysics model to simulate the coupled electrochemistry and multiphase transport in an electrolysis cell operated with the novel PTL architecture. The results show that larger pores improve access of water to the anode catalyst layer, which is beneficial for both the oxygen evolution reaction and membrane hydration. Larger pore sizes also improve oxygen gas transport from the catalyst layer, because generated oxygen gas is forced to travel in-plane through the anode catalyst layer until it reaches a pore opening that is connected to a channel. The discussed results confirm that the proposed thin foil based PTLs are fundamentally different from conventional PTLs, such as felts or layered meshes. The model developed in this work also provides generalizable insight into fundamental PEMEC phenomena, such as the competition between liquid and gas phase transport, membrane hydration and water management, and nonuniform electrochemical reactions, which are processes relevant to all PEMEC designs.

25 ENERGY STORAGE↗
NASA NTRS2022

An Investigation Into Some Important Aspects of Droplet Breakup/Vaporization Behavior Caused By a Gas Flow (Both Nonreacting & Detonative Combustion) in a Shock Tube

As a part of the rotating detonation engine (RDE) technology enablement project at NASA Glenn Research Center (GRC), an effort was undertaken to extend our current computational capabilities of OpenNCC in some important ways with the implementation of a modeling approach to account for the droplet breakup caused by a shock-induced gas motion, & a vaporization model valid over a wide range of pressure conditions encountered in multiphase detonation. With the modified code, a study was undertaken to investigate the individual droplet behavior followed by the passage of a shock front. The study is carried out by tracking a sparse group of droplets to gain some understanding of shock induced droplet behavior under various shock strengths & fuel injector conditions. The study also looks into the effect of randomization involved in determining the droplet breakup outcomes. Over a wide range of sub-critical conditions examined, larger droplets are observed to undergo significant changes in droplet behavior following their breakup. However, smaller droplets (10 µm or less ) remain unaffected by any shock induced breakup. In a follow-on work, we investigated the impact of shock and droplet interaction in a detonation study involving both gaseous as well as gas/liquid (droplet clouds) fuel/air stoichiometric mixtures in a simple 3D shock-tube configuration. The droplet clouds are made up of different initial droplet sizes of either 6, 10, or 30 µm. We also investigated the individual droplet behavior followed by the passage of a detonation front. The results represent conditions that lead to both overdriven and C-J (ChapmanJouguet) detonations. Under both test conditions, most of the droplet vaporization is completed within a short distance (duration) behind the detonation front & well within the region of complete combustion observed in a corresponding equivalent gas-phase fuel/air mixture. The impact of the shock-induced droplet breakup is found to be significant in the calculations involving the 30-µm droplets. Subsequent to the breakup, the drop sizes vary from 1 to 10 µm. Another factor that contributed to the observed rapid vaporization is the result of vaporization taking place under supercritical conditions. The overall detonation properties of various droplet clouds (made up of different initial sizes) are similar to those observed in a corresponding gaseous fuel/air mixture. In the calculation involving a gaseous fuel, the calculated C-J detonation velocity is 1822 m/s involving Jet-A/air and φ = 1. In the overdriven detonation, it is 2044 m/s. In the calculations involving droplet clouds, the corresponding detonation velocities are lower. The impact of increased droplet size is primarily seen in a higher reduction in the detonation velocity.

shock/droplet interaction↗
NASA NTRSDate not supplied

CFD Model Development of a Cryogenic Storage Tank Self-Pressurization in Normal Gravity and Validation against SHIIVER Experiment

Two-Phase flow and heat transfer simulations with interfacial phase change of the Structural Heat Intercept, Insulation, and Vibration Evaluation Rig (SHIIVER) self-pressurization experiment were conducted using storage tank CFD model in the framework of the ANSYS Fluent CFD code. The simulations were performed for the 70% fill level case with MLI on domes and no vapor cooling. All the phase change calculations in these simulations were generated by in-house Schrage-based evaporation-condensation model. The calculations were performed using both the VOF and Sharp Interface multiphase 2Daxisymmetric models. A number of parametric and sensitivity studies were performed to check the various aspects of the CFD model. These studies helped to understand the effects of varying several parameters on the tank pressure and temperature during self-pressurization. These parameters include: the effect of turbulence modeling and turbulence damping at the interface; the effect of constant vs. temperature dependent fluid properties; the effect of accommodation coefficient; the effect of modeling phase change at the interface and the effect of the fill level.

Computational Fluid Dynamics↗
DOE OSTI2024

Surface Structure of Reacting Solid Ramjet Fuels

Solid fuel ramjets are a developing technology that will enable extended range for supersonic vehicles. The simplicity of the system in terms of moving parts results in the performance being dictated by the solid fuel combustion processes. The physics near the surface, particularly in the fuel-rich recirculation zone, are challenging to resolve due to flow opacity. Here, regression rate measurements of hydroxyl-terminated polybutadiene from a miniature, portable slab burner are reported. Mass-averaged regression rates agree with correlations from the literature, suggesting that the proper physics in the burner are captured. Synchrotron-based phase-contrast imaging is used to quantify local surface features, including a multiphase froth layer, molten polymer transport, and local regression rate measurements in the recirculation zone. Regression rates near the inlet are less than the mass-averaged values, whereas near the reattachment point the local value exceeds the mass-averaged values. Further, the observed froth layer is as large as 2 mm near the inlet and decreases to sub-millimeter thicknesses near the reattachment point. An analytical heat transfer model is used to quantify the influence of the surface froth layer on the regression rates. Results indicate that the froth layer has an insulating effect and will reduce the local regression rate.

33 ADVANCED PROPULSION SYSTEMS↗