Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow regime”

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 271 records · Page 15

Dynamic control of molecular transport MXene transistor membranes

Controlled spatial confinement and surface properties of lamellar 2D nanomaterial membranes could enhance many precision separation processes. Traditionally, researchers view channel dimensions, surface properties, and permeation rates of these membranes as intrinsic properties that cannot be modulated in operando. We report that gate voltage applied to the conducting laminar MXene membrane can modulate the permeation rate of ions and neutral solutes, as well as its effective size rejection. In operando wide-angle x-ray scattering measurements reveal that these changes are not driven by electrically induced variations in the d spacing of the MXene layers. Instead, experimental data and continuum electrokinetic modeling reveal that ion transport through the MXene channels is primarily affected by Donnan equilibrium at the membrane-solution interface. We also report a strong increase in the permeation rates through the membrane under a low-frequency ac voltage gating regime that we attribute to diffusioosmotic flow oscillations induced in the membrane. Overall, MXene “transistor” membranes provide a previously unidentified approach to dynamic control of molecular separations.

Pendse, Aaditya [Lawrence Livermore National Labor↗

Power Generation Cycle with RELAP-7

Balance of plant of a High-Temperature gas-cooled Reactor – Pebble bed Module (HTR-PM) is modeled using RELAP-7. The model includes the primary helium loop where the heat produced by the core is transferred to a helical-coil steam generator and the secondary loop where super-heated steam is fed to a turbine to extract power. The steam is then condensed to subcooled water and pumped back into the steam generator. Results demonstrate that RELAP-7 has the capability to model different regimes encountered in two-phase flow with wall boiling, superheated steam, subcooled liquid and condensation. RELAP-7 also has the capability of modeling different fluids with their own specific models in a very flexible way, without any code change. This also demonstrates the capability of modeling the transfer of fission power into electrical power.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimizing the Heisenberg Vortex Tube for Hydrogen Cooling (Final Technical Report)

Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.

08 HYDROGEN↗

FY-24 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the computational modeling of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY24. This year’s modeling campaign focuses on the continual testing of the RELAP5 model against experimental data. Several fault cases were simulated with the RELAP5 model in FY23 to study the accuracy of the model under complex flow conditions. Similarly, in FY24, a case is simulated where the chimney outlet is throttled with a valve while the flow is undergoing two-phase oscillations. The purpose of this case is to study the effects of increasing pressure drop in the two-phase region on the system behavior. Even though the model is able to capture the experimental data qualitatively, it overpredicts the pressure drops experienced by the flow in the two-phase regime. Nevertheless, the model and the experiment show that oscillations are stabilized with increased pressure drop in the two-phase region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Transport Simulations on Scanning Transmission Electron Microscope Images of Nanoporous Shale

Digital rock physics is an often-mentioned approach to better understand and model transport processes occurring in tight nanoporous media including the organic and inorganic matrix of shale. Workflows integrating nanometer-scale image data and pore-scale simulations are relatively undeveloped, however. In this paper, a workflow is demonstrated progressing from sample acquisition and preparation, to image acquisition by Scanning Transmission Electron Microscopy (STEM) tomography, to volumetric reconstruction to pore-space discretization to numerical simulation of pore-scale transport. Key aspects of the workflow include (i) STEM tomography in high angle annular dark field (HAADF) mode to image three-dimensional pore networks in µm-sized samples with nanometer resolution and (ii) lattice Boltzmann method (LBM) simulations to describe gas flow in slip, transitional, and Knudsen diffusion regimes. It is shown that STEM tomography with nanoscale resolution yields excellent representation of the size and connectivity of organic nanopore networks. In turn, pore-scale simulation on such networks contributes to understanding of transport and storage properties of nanoporous shale. Interestingly, flow occurs primarily along pore networks with pore dimensions on the order of tens of nanometers. Smaller pores do not form percolating pathways in the sample volume imaged. Apparent gas permeability in the range of 10−19 to 10−16 m2 is computed.

42 ENGINEERING↗

Infrared fixed point of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors

We use lattice simulations and the continuous renormalization group method, based on the gradient flow, to calculate the 𝛽 function and anomalous dimensions of the SU(3) gauge theory with 𝑁 𝑓 = 10 flavors of fermions in the fundamental representation. We employ several improvements to extend the range of available renormalized couplings, including the addition of heavy Pauli-Villars bosons to reduce cutoff effects and the combination of a range of gradient flow transformations. While in the weak coupling regime our result is consistent with those of earlier studies, our techniques allow us to study the system at much stronger couplings than previously possible. We find that the renormalization group 𝛽 function develops a zero, corresponding to an infrared-stable fixed point, at gradient-flow coupling 𝑔 2 = 15.0⁢(5). We also determine the mass and tensor anomalous dimensions: At the fixed point we find 𝛾 𝑚 ≃0.6, suggesting that this system might be deep inside the conformal window.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Detecting depinning and nonequilibrium transitions with unsupervised machine learning

Using numerical simulations of a model disk system, we demonstrate that a machine learning generated order-parameter-like measure can detect depinning transitions and different dynamic flow phases in systems driven far from equilibrium. Furthermore, we specifically consider monodisperse passive disks with short range interactions undergoing a depinning phase transition when driven over quenched disorder. The machine learning derived order-parameter-like measure identifies the depinning transition as well as different dynamical regimes, such as the transition from a flowing liquid to a phase separated liquid-solid state that is not readily distinguished with traditional measures such as velocity-force curves or Voronoi tessellation. The order-parameter-like measure also shows markedly distinct behavior in the limit of high density where jamming effects occur. Our results should be general to the broad class of particle-based systems that exhibit depinning transitions and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗

Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel

In this LDRD “Materials response under hypersonic flow conditions probed by multi-modal diagnostics in a benchtop wind tunnel” a peer-reviewed lab-scale hypersonic wind tunnel (WT) was planned to be assembled in-house, along with time and space resolved and related flow and materials diagnostics to derive high quality measurements and the science that is able to inform the development of our modeling and simulation codes for the hypersonic regime. The baseline capability to achieve Mach 3 flow was scheduled to be completed in Q4 2020, in the first year of the LDRD. However due to the COVID-19 delays and travel restrictions starting in March, 2020 this task could not be completed to carry out the materials response studies in-house, or those planned with external partners that were engaged during this LDRD (NASA Langley, NASA Ames, Texas A&M university, and Virginia Tech). Still, using our hypersonic computational fluid simulation codes, many of the hypersonic components have been designed (for Mach 3, 5, and 7) and procured (or in procurement phase), and a new dedicated hypersonic laboratory was facilitized, 2 WT test sections have been manufactured, and a Mach 3 nozzle was 3D printed for prototyping. In parallel, unique materials diagnostics have been tested under static conditions and fully specified and awaiting validation in our partner facilities. Due to the new hypersonic initiative announced in January, 2020 we proposed a related scope increase in hypersonics research requiring energy interactions and new diagnostics well beyond what could be supported by an LDRD ER. A decision was made to end the current LDRD ER and formally propose a new LDRD SI that leverages the current effort, which will be starting in FY21. As a result of the uniqueness of this LDRD work, our early efforts have been encouraged and well received in the hypersonic research community, and is now part of multiple unsolicited internal and external project proposals expected to be funded in the coming years, which is encouraging to make LLNL a leader in experimentally probing extreme physics of local energetics critical to hypersonics.

36 MATERIALS SCIENCE↗

Controlling solute channel formation using magnetic fields

Solute channel formation introduces compositional and microstructural variations in a range of processes, from metallic alloy solidification, to salt fingers in ocean and water reservoir flows. Applying an external magnetic field interacts with thermoelectric currents at solid/liquid interfaces generating additional flow fields. This thermoelectric (TE) magnetohydrodynamic (TEMHD) effect can impact on solute channel formation, via a mechanism recently drawing increasing attention. To investigate this phenomenon, we combined in situ synchrotron X-ray imaging and Parallel-Cellular-Automata-Lattice-Boltzmann based numerical simulations to study the characteristics of flow and solute transport under TEMHD. Observations suggest the macroscopic TEMHD flow appearing ahead of the solidification front, coupled with the microscopic TEMHD flow arising within the mushy zone are the primary mechanisms controlling plume migration and channel bias. Two TE regimes were revealed, each with distinctive mechanisms that dominate the flow. Further, we show that grain orientation modifies solute flow through anisotropic permeability. These insights led to a proposed strategy for producing solute channel-free solidification using a time-modulated magnetic field.

36 MATERIALS SCIENCE↗

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The reactive flow evolution of the polymer-bonded explosive PBX 9502: Experiments and model validation in extreme pressure regimes

The shock-to-detonation transition properties of the triaminotrinitrobenzene based PBX 9502 high explosive (HE) are experimentally and computationally explored in extremely high input pressure conditions. These include both slightly sub-Chapman–Jouguet and overdriven input pressure conditions, namely, ∼25 and ∼31 GPa, respectively. Our experiments capture the transient buildup of a shock-induced reaction via measurement of HE and polymethyl methacrylate window interface particle velocity profiles for a variety of sample thicknesses for this insensitive HE. These observations necessitate extremely thin explosive samples, and the high rates of reaction provide a considerable challenge to optical diagnostics. Samples at these thicknesses also provide an opportunity for evaluation of potential micro-structure effects on the resulting shock-to-detonation-transition measurements. To address this, the thin samples are also characterized via x-ray micro-computed tomography. Finally, a pair of previously established continuum-level detonation performance modeling approaches for PBX 9502 were used to analyze the experiments. The employed model variants crucially differ in their definition of each model’s empirical reaction rate functional form, utilization of shock state quantities, and local flow variable dependencies. As a result, the present experiments provide a novel platform to evaluate the quantitative and qualitative consequences stemming from these modeling choices in a challenging initiation scenario, largely beyond the chosen calibration range of either model. This new experimental information will provide a platform for both improved physics and model parameterizations for this well-studied explosive.

36 MATERIALS SCIENCE↗

Advanced manufacturing for electrosynthesis of fuels and chemicals from CO 2

Advanced manufacturing (AM) represents an appealing approach for creating novel electrochemical systems for chemicals synthesis. In this work, we demonstrate AM for rapid development and testing for improved performance for the carbon dioxide reduction reaction across an evolution of vapor-fed reactor designs. In our final design, we observe activation- and mixed-control regimes for a variety of operating conditions via inlet CO 2 flow rate and electrochemical potential. Furthermore, we define a dimensionless number (Da) to identify mass transport regimes by exploring the impact of hypothesized underlying mass transport mechanisms, including consumption of CO 2 via OH - , increased local temperatures, and partial penetration of electrolyte into gas diffusion layer. The accelerated pace of reactor design and development led to high geometric current densities (500 mA cm -2 ), heightened selectivity (85.5% FE C 2+ products), and increased carbon dioxide conversion (16.6%) and cathodic energy efficiency (49.6% CO 2 R). Using AM vapor-fed reactors, we attain high ethylene (3.67%) and record ethanol (3.66%) yields compared to the literature. Finally, this work underscores the promise of AM for accelerating reactor design, understanding of governing phenomena, and improving the performance of catalytic systems.

36 MATERIALS SCIENCE↗

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

Dimensionless Regime Mapping and Chaos Analysis of Solids Volume Fraction Data in a Circulating Fluidized Bed Riser

Presentation of experimental results of solids volume fraction data measurements in a 12-inch diameter cold flow circulating fluidized bed in relation to a dimensionless fluidization regime map. Statistical and chaotic system analysis results are shown to illustrate differences in hydrodynamic behavior at different operating conditions. Presentation was made at the ASTFE 10th Thermal and Fluids Engineering Conference, March 09-12, 2025.

chaos↗

An attention-based neural ordinary differential equation framework for modeling inelastic processes

To preserve strictly conservative behavior as well as model the variety of dissipative behavior displayed by solid materials, we propose a significant enhancement to the internal state variable-neural ordinary differential equation (ISV-NODE) framework. In this data-driven, physics-constrained modeling framework internal states are inferred rather than prescribed. The ISV-NODE consists of: (a) a stress model dependent on observable deformation and inferred internal state, and (b) a model of the evolution of the internal states. The enhancements to ISV-NODE proposed in this work are multifold: (a) a partially input convex neural network stress potential provides polyconvexity in terms of observed strain while leaving the inferred state unconstrained, and (b) an internal state flow model uses common latent features to inform novel attention-based gating and drives the flow of internal state only in dissipative regimes. We demonstrated that this architecture can accurately model dissipative and conservative behavior across an isotropic, isothermal elastic-viscoelastic-elastoplastic spectrum with three exemplars, while maintaining fundamental principles by design.

97 MATHEMATICS AND COMPUTING↗