Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical mixing”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Simulating many-engine spacecraft: Exceeding 1 quadrillion degrees of freedom via information geometric regularization

We present an optimized implementation of the recently proposed information geometric regularization (IGR) for unprecedented scale simulation of compressible fluid flows applied to multi-engine spacecraft boosters. We improve upon state-of-the-art computational fluid dynamics (CFD) techniques in terms of computational cost, memory footprint, and energy-to-solution metrics. Unified memory on coupled CPU–GPU or APU platforms increases problem size with negligible overhead. Mixed half/single-precision storage and computation are used on well-conditioned numerics. We simulate flow at 200 trillion grid points and 1 quadrillion degrees of freedom, exceeding the current record by a factor of 20. A factor of 4 wall-time speedup is achieved over optimized baselines. Ideal weak scaling is observed on OLCF Frontier, LLNL El Capitan, and CSCS Alps using the full systems. Strong scaling is near ideal at extreme conditions, including 80% efficiency on CSCS Alps with an 8 node baseline and stretching to the full system.

Wilfong, Benjamin [Georgia Institute of Technology↗
DOE OSTI2023

Tidal energy extraction modifies tidal asymmetry and transport in a shallow, well-mixed estuary

Tidal energy extraction is increasingly being studied as a potential renewable energy resource in estuaries worldwide. Although it is understood that energy extraction via tidal stream turbines can modify currents and transport within estuaries, it is not clear how the underlying nonlinear physical mechanisms dictating tidal hydrodynamics are modulated. This research investigates the influence of a hypothetical tidal stream turbine array on barotropic tidal processes in a shallow, well-mixed system: the Piscataqua River – Great Bay (PRGB) estuary, using a numerical model. The modeled turbine farm includes 180 turbines which would extract an estimated 44.7 GWh of energy, annually. The tidal hydrodynamic model for the existing condition is validated with in-situ observations of currents and water level before analyzing tidal asymmetry and transport with and without tidal turbines. Results indicate that the tidal turbine array will decrease tidal elevation and current magnitudes system-wide, but generally reduce ebb currents and transport more than flood over most of the estuary footprint, thereby diminishing tidal asymmetry. The smaller asymmetric distortion compared to the no-turbine case is attributed to reductions in the storage volume of water over the estuary’s extensive tidal flat regions between low and high waters which decreases the associated nonlinear intertidal storage mechanism up to 25%. This leads to weakened ebb dominance over estuary sections from the mouth to mid-reaches, where depths are deep enough to keep the combined nonlinear shallow water and frictional effects from asserting control over the storage mechanism. Even in upstream shallow regions where depth-dependent friction controls asymmetry in both cases, the frictional mechanism is reduced only by 10% with turbines. Some environmental considerations of this work are discussed, with focus on sediment transport, water quality, and ecology.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Self-similarity of scalar isosurface area density in a temporal mixing layer

Understanding the behaviour of scalar isosurfaces in a turbulent flow is of particular interest for many problems in turbulent mixing that contain sharp interfaces between regions of the flow. Common examples include combustion, where the chemical reactions occur in thin regions within the flow, and the turbulent/non-turbulent interface in shear flows, where a thin region separates the rotational, turbulent motions from the irrotational, non-turbulent background. Recent advances in computing technology allow for in-depth analysis of these interface problems that are difficult to quantify in a laboratory setting. In this paper, the results of a direct numerical simulation of a passive scalar $\varPhi$ evolving on a turbulent, temporally developing mixing layer are described. A novel approach has been taken to calculate the surface area of individual scalar isosurfaces, $A_{iso}$ , throughout the simulation, as well as the mean isosurface area density, $\varSigma$ , as a function of the cross-stream direction and time. A notable finding is that the profiles of $\varSigma$ develop in a self-similar manner when scaled by the Taylor scale of the scalar field, $\lambda _\phi$ . Remarkably, the scaling appears to hold for a wide range of isovalues. A rough scaling argument based on the formal definition of $\varSigma$ and properties of a temporal mixing layer is presented which also exposes a dependence on $\lambda _\phi$ . Based on these results, a possible scaling for the isosurface area is presented as $A_{iso}/A_0 \sim (Re \, Sc)^{1/2}$ , where $Re$ and $Sc$ are local Reynolds and Schmidt numbers, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Review and Analysis of Heat Transfer Correlations For Horizontal Pseudocritical CO 2 Heat Exchanger Applications

There is significant interest in the development of supercritical carbon dioxide (sCO 2 ) power cycles because of the potential for smaller and more energy efficient systems than a steam Rankine cycle. Heat exchanger designs typically use empirical correlations, but the applicability of these correlations near the CO 2 critical point is a potential issue. Though numerous correlations have been proposed in the literature, there are some disagreements when it comes to the accuracy. The current work recognizes the role of thermophysical properties, and its impacts on the heat transfer correlations and cycle efficiency. Heat transfer correlations proposed for horizontal flow inside circular pipes were analyzed with the help of numerical simulations. Steady state RANS simulations were performed using SST k-w turbulence model to evaluate the Nusselt number empirical correlations. It was found that the most of correlations (except Yoon) produced a Nusselt number that differed significantly with the one predicted numerically. Some of the correlations were developed for pure forced convection regime and as mentioned in Lin et al. do not account for mixed convection or free convection effects. Based on the limited observation, it appears that Yoon et al. predictions match well with the numerically predicted Nusselt Numbers. However, further analysis is required understand the applicability of various correlations and is contingent on accurate measurements or predictions of wall temperature profiles in the axial and the circumferential directions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2025

Mixed-valence halide perovskites

Mixed-valence compounds—which feature an element in at least two different oxidation states—can display emergent optical and transport phenomena stemming from electron transfer between the different valences (intervalence charge-transfer; IVCT). As halide perovskites show promise as active materials in numerous optoelectronic devices, it is an opportune moment to incorporate and study the effects of mixed-valence in this versatile materials family, to access tunable electronic structures ranging from insulators to semiconductors to metals. Herein, we introduce the basic concepts of mixed-valence in molecules and discuss how these concepts may be extended to mixed-valence in extended solids. We then review the few studies of mixed valence in 3D and 2D halide perovskites and halide perovskites with mixed-valence impurities, ranging from studies in the early 1900s to the present day. In conclusion, through judicious choice of metal ion, its coordinating ligands and their geometry, and overall structural dimensionality, chemists can exert powerful synthetic control over electronic delocalization in mixed-valence perovskites, and we hope to see this intriguing materials class expand to encompass new compositions.

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

Impulse-driven release of gas-encapsulated drops

Gas-encapsulated drops, much like antibubbles, are drops enclosed in a bubble within a liquid. They show potential as payload carriers in fluid transport and mixing techniques where sound waves can be leveraged to induce the collapse of the gas core and the subsequent release of the drop. Here, the interaction of millimetre-sized gas-encapsulated drops with impulsive laser-induced shock waves is investigated to gain fundamental insights into the release process. Experimental synchrotron X-ray phase contrast imaging, which allows the drop dynamics to be visualised inside the encapsulating bubble, is complemented by numerical simulations to study the intricate physics at play. Three drop dynamical release regimes are discovered, namely the drop impact , partial deposition and jet impact regimes. The regime type is mainly dependent on the shape of the bubble interface impacting the drop and the associated Weber and Reynolds numbers. The drop dynamics of the drop impact and partial deposition regimes show similarities with the canonical configuration of drops impacting flat liquid surfaces, whereas the jet impact regime resembles binary drop collisions, which allows existing scaling laws to be applied to describe the underlying processes. The release of the drop is investigated numerically. The time evolution of the drop dissemination within the surrounding liquid discloses enhanced mixing for dynamics involving high Weber and Reynolds numbers such as the drop impact and jet impact regimes.

42 ENGINEERING↗
DOE OSTI2021

Top-level physics requirements and simulated performance of the MRSt on the National Ignition Facility

The time-resolving Magnetic Recoil Spectrometer (MRSt) for the National Ignition Facility (NIF) has been identified by the US National Diagnostic Working Group as one of the transformational diagnostics that will reshape the way inertial confinement fusion (ICF) implosions are diagnosed. The MRSt will measure the time-resolved neutron spectrum of an implosion, from which the time-resolved ion temperature, areal density, and yield will be inferred. Top-level physics requirements for the MRSt were determined based on simulations of numerous ICF implosions with varying degrees of alpha heating, P2 asymmetry, and mix. Synthetic MRSt data were subsequently generated for different configurations using Monte–Carlo methods to determine its performance in relation to the requirements. The system was found to meet most requirements at current neutron yields at the NIF. This work was supported by the DOE and LLNL.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modulation-slippage trade-off in resonant four-wave upconversion

Following up on a proposal to use four-wave mixing in an underdense plasma at mildly relativistic laser intensities to produce vastly more energetic x-ray pulses [V. M. Malkin and N. J. Fisch, Phys. Rev. E, 101, 023211 (2020)], we perform the first numerical simulations in one dimension to illustrate amplification of a short high frequency seed through four-wave mixing. Here, we illustrate how parasitic processes including phase modulation and spatial pulse slippage limit the amplification efficiency. Although the regimes studied were not where the optimal efficiencies were expected, these regimes do expose the basic physical processes at play, while still yielding not insignificant spectral power upshift.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simplifying the creation of thermal decomposition mechanisms for designing mixed fibre composites

Among the most challenging aspects of simulating thermal decomposition of fibre reinforced polymers is the determination of appropriate reaction parameters. Thermogravimetric analysis is typically used to generate decomposition data, to which the reaction parameters are then fit. When designing mixed fibre materials (e.g. combinations of glass and carbon fibres), the number of TGA experiments needed to explore the entire design space may be intractable. Here, we demonstrate the creation of two candidate proxy mixed fibre mechanisms and compare them to a mechanism created by fitting parameters from TGA on the mixed fibre composite. These mechanisms are then demonstrated in a 2D axisymmetric numerical decomposition, heat transfer, and porous flow model. We find a maximum 11% uncertainty in mass and 4% in temperature difference when using a proxy mechanism.

Scott, Sarah N. [Sandia National Lab. (SNL-CA), Li↗
DOE OSTI2023

Mixed-Precision S/DGEMM Using the TF32 and TF64 Frameworks on Low-Precision AI Tensor Cores

Using NVIDIA graphics processing units (GPUs) equipped with Tensor Cores has enabled the significant acceleration of general matrix multiplication (GEMM) for applications in machine learning (ML) and artificial intelligence (AI) and in high-performance computing (HPC) generally. The use of such power-efficient, specialized accelerators can provide a performance increase between 8 × and 20 ×, albeit with a loss in precision. However, a high level of precision is required in many large scientific and HPC applications, and computing in single or double precision is still necessary for many of these applications to maintain accuracy. Fortunately, mixed-precision methods can be employed to maintain a higher level of numerical precision while also taking advantage of the performance increases from computing with lower-precision AI cores. With this in mind, we extend the state of the art by using NVIDIA’s new TF32 framework. This new framework not only burdens some constraints of the previous frameworks, such as costly 32 16-bit castings but also provides an equivalent precision and performance by using a much simpler approach. We also propose a new framework called TF64 that attempts double-precision arithmetic with low-precision Tensor Cores. Although this framework does not exist yet, we validated the correctness of this idea and achieved an equivalent of 64-bit precision on 32-bit hardware.

Valero Lara, Pedro↗
DOE OSTI2022

On the Spectral Evolution of Hot White Dwarf Stars. II. Time-dependent Simulations of Element Transport in Evolving White Dwarfs with STELUM

White dwarf stars are subject to various element transport mechanisms that can cause their surface composition to change radically as they cool, a phenomenon known as spectral evolution. In this paper, we undertake a comprehensive theoretical investigation of the spectral evolution of white dwarfs. First, we introduce STELUM, a new implementation of the stellar evolutionary code developed at the Université de Montréal. We provide a thorough description of the physical content and numerical techniques of the code, covering the treatment of both stellar evolution and chemical transport. Then, we present two state-of-the-art numerical simulations of element transport in evolving white dwarfs. Atomic diffusion, convective mixing, and mass loss are considered simultaneously as time-dependent diffusive processes and are fully coupled to the cooling. We first model the PG 1159-DO-DB-DQ evolutionary channel: a helium-, carbon-, and oxygen-rich PG 1159 star transforms into a pure-helium DB white dwarf due to gravitational settling and then into a helium-dominated, carbon-polluted DQ white dwarf through convective dredge-up. We also compute for the first time the full DO-DA-DC evolutionary channel: a helium-rich DO white dwarf harboring residual hydrogen becomes a pure-hydrogen DA star through the float-up process and then a helium-dominated, hydrogen-bearing DC star due to convective mixing. We demonstrate that our results are in excellent agreement with available empirical constraints. In particular, our DO-DA-DC simulation perfectly reproduces the lower branch of the bifurcation observed in the Gaia color–magnitude diagram, which can therefore be interpreted as a signature of spectral evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Top-level physics requirements and simulated performance of the MRSt on the National Ignition Facility

The time-resolving Magnetic Recoil Spectrometer (MRSt) for the National Ignition Facility (NIF) has been identified by the US National Diagnostic Working Group as one of the transformational diagnostics that will reshape the way inertial confinement fusion (ICF) implosions are diagnosed. The MRSt will measure the time-resolved neutron spectrum of an implosion, from which the time-resolved ion temperature, areal density, and yield will be inferred. Top-level physics requirements for the MRSt were determined based on simulations of numerous ICF implosions with varying degrees of alpha heating, P2 asymmetry, and mix. Synthetic MRSt data were subsequently generated for different configurations using Monte–Carlo methods to determine its performance in relation to the requirements. The system was found to meet most requirements at current neutron yields at the NIF. This work was supported by the DOE and LLNL.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Low Precision for Lower Energy Consumption: Preprint

Low-precision numeric types offer significant efficiency and energy benefits for computing applications. Mixed-precision algorithms, combining low and high precision types, maintain accuracy while improving performance. Despite advantages, there exist challenges on adapting existing mixed-precision algorithms to new technologies, such as new hardware architectures and new low-precision data types. This paper presents current challenges and opportunities to advance science in this domain targetting more energy efficient solutions.

energy efficiency↗
DOE OSTI2022

Online Learning and Pricing for Network Revenue Management with Reusable Resources

We consider a price-based network revenue management problem with multiple products and multiple reusable resources. Each randomly arriving customer requests a product (service) that needs to occupy a sequence of reusable resources (servers). We adopt an incomplete information setting where the firm does not know the price-demand function for each product and the goal is to dynamically set prices of all products to maximize the total expected revenue of serving customers. We propose novel batched bandit learning algorithms for finding near-optimal pricing policies, and show that they admit a near-optimal cumulative regret bound of $O(J\sqrt{XT})$, where J, X, and T, are the numbers of products, candidate prices, and service periods, respectively. As part of our regret analysis, we develop the first finite-time mixing time analysis of an open network queueing system (i.e., the celebrated Jackson Network), which could be of independent interest. Our numerical studies show that the proposed approaches perform consistently well.

42 ENGINEERING↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2026

Review of Particle Physics

The Review summarizes much of particle physics and cosmology. Using data from previous editions, plus 3,200 new measurements from 903 papers, we list, evaluate, and average measured properties of gauge bosons and the recently discovered Higgs boson, leptons, quarks, mesons, and baryons. We summarize searches for hypothetical particles such as supersymmetric particles, heavy bosons, axions, dark photons, etc. Particle properties and search limits are listed in Summary Tables. We give numerous tables, figures, formulae, and reviews of topics such as Higgs Boson Physics, Supersymmetry, Grand Unified Theories, Neutrino Mixing, Dark Energy, Dark Matter, Cosmology, Particle Detectors, Colliders, Probability and Statistics. Most of the 118 reviews are updated, including many that are heavily revised.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Flame Characterization of a Hydrogen Micromix Burner using Experimental and Numerical Methods - Extended Abstract

Micromix combustion presents a promising solution for integrating hydrogen into gas turbine combustors, offering rapid and compact air-fuel mixing that effectively shortens the combustion zone and significantly reduces nitrogen oxides (NOX) emissions. A major challenge in hydrogen combustion systems is flashback, which is mitigated through short mixing lengths and high flow velocities. This study investigates three fuel injector designs – zero swirl, low swirl, and high swirl – for a 6-mm inner diameter nozzle, where hydrogen is injected through six 0.3 mm diameter holes in a jet-in-crossflow configuration relative to the airstream. Numerical simulations indicate stable flame formation, moderate peak flame temperatures, and subsequent low NOx potential for all designs. These injectors have been fabricated using additive manufacturing, and upcoming experiments will explore flame structure and NOx behavior across varying fuel injection angles, equivalence ratios and pre-mixer velocities. These efforts aim to advance the development of low-emission hydrogen micromix combustion.

hydrogen↗