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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Data-driven Mori–Zwanzig modeling of Lagrangian particle dynamics in turbulent flows

The dynamics of Lagrangian particles in turbulence play a crucial role in mixing, transport, and dispersion in complex flows. Their trajectories exhibit highly nontrivial statistical behavior, motivating the development of surrogate models that can reproduce these trajectories without incurring the high computational cost of direct numerical simulations of the full Eulerian field. This task is particularly challenging because reduced-order models typically lack access to the full set of interactions with the underlying turbulent field. Novel data-driven machine learning techniques can be powerful in capturing and reproducing complex statistics of the reduced-order/surrogate dynamics. In this work, we show how one can learn a surrogate dynamical system that is able to evolve a turbulent Lagrangian trajectory in a way that is point-wise accurate for short-time predictions (with respect to Kolmogorov time) and stable and statistically accurate at long times. This approach is based on the Mori–Zwanzig formalism, which prescribes a mathematical decomposition of the full dynamical system into resolved dynamics that depend on the current state and the past history of a reduced set of observables, and the unresolved orthogonal dynamics due to unresolved degrees of freedom of the initial state. We show how by training this reduced order model on a point-wise error metric on short time-prediction, we are able to correctly learn the dynamics of Lagrangian turbulence, such that also the long-time statistical behavior is stably recovered at test time. This opens up a range of applications, for example, for the control of active Lagrangian agents in turbulence.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Rayleigh–Taylor Instability With Varying Periods of Zero Acceleration

We present our findings from a numerical investigation of the acceleration-driven Rayleigh–Taylor Instability, modulated by varying periods without an applied acceleration field. It is well known from studies on shock-driven Richtmyer–Meshkov instability that mixing without external forcing grows with a scaling exponent as ≈ t 0.20-0.28 When the Rayleigh–Taylor Instability is subjected to varying periods of “zero” acceleration, the structural changes to the mixing layer remain remarkably small. After the acceleration is re-applied, the mixing layer quickly resumes the profile of development it would have had if there had been no intermission. As a result, this behavior contrasts in particular with the strong sensitivity that is found to other variable acceleration profiles examined previously in the literature.

42 ENGINEERING↗

DOE OSTI2021

Research Priorities and Opportunities in United States Wholesale Electricity Markets

The power system is currently undergoing a number of changes, including a rapidly evolving resource mix, growth of distributed energy resources (DERs), more active consumer participation, increased deployment of energy storage and hybrid resources, and more advanced communication and control requirements. These changes in the power system present numerous technical, economic, implementation, and policy challenges and research opportunities for power system operators. To help address these challenges, a collaboration among five research institutions—Argonne National Laboratory, the National Renewable Energy Laboratory, Lawrence Berkeley National Laboratory, the Electric Power Research Institute, and Johns Hopkins University—has been established to provide technical assistance to the seven U.S. Independent System Operators (ISOs) and Regional Transmission Organizations (RTOs). The 3-year project aims to leverage the advanced methods, tools, datasets, and resources of the collaborators to provide robust analytical support to address the high-priority market challenges that will be faced in the 2- to 10-year time horizon.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2022

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

A Predictive Deep-Reinforcement-Learning-Based Connected Automated Vehicle Anticipatory Longitudinal Control in a Mixed Traffic Lane Change Condition

Maintaining safety and efficiency for mixed traffic consisting of connected automated vehicles (CAVs) and human-driven vehicles (HDVs) is an arduous task due to the inherent HDVs’ stochasticity. Especially for longitudinal control, which is the basic function of vehicle automation, prevailing research primarily considers CAV’s car-following control merely the acceleration and deceleration of leading vehicles. However, this approach overlooks the potential disruptions caused by surrounding vehicles executing lane changes, which can significantly impact the control vehicle’s stability and overall safety. Hence, our study introduces a predictive deep reinforcement learning (DRL) longitudinal CAV controller. This innovative approach leverages prediction from a physics-informed neural network as well as the control capability of DRL to better anticipate and mitigate issues arising from lane-changing, enhancing the safety and efficiency of CAVs in such scenarios. Finally, validated by the numerical simulations embedded with the real-world data, the results indicate that the proposed controller significantly enhances the safety and efficiency of CAVs in situations involving lane changes by other vehicles, showcasing its potential as a valuable tool in advancing CAV technology in mixed traffic.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2021

Feasibility of full-core pin resolved CFD simulations of small modular reactor with momentum sources

Complex flow structure interactions and heat transfer processes take place in nuclear reactor cores. Given the extreme pressure/temperature and radioactive conditions inside the core, numerical simulations offer an attractive and sometimes more feasible approach to study the related flow and heat transfer phenomena in addition to the experiments. Under the Exascale Computing Project, the full-core simulation of a small modular reactor (SMR) has been pursued coupling Computational Fluid Dynamics (CFD) and neutronics. A key aspect of the modeling of SMR fuel assemblies is the presence of spacer grids and the mixing promoted by mixing vanes or the equivalent. A reduced order methodology is adopted based on momentum sources to mimic the mixing of the vanes. The momentum sources have been carefully calibrated with detailed Large Eddy Simulations (LES) of spacer grids performed with Nek5000. Modeling the spacer grid and mixing vanes (SGMV) effect without body-fitted computational grid avoids the excessive costs in resolving the local geometric details, and thus supports the simulation to be scaled up to the full core. Besides the progress on momentum source modeling, this paper also features the first full-core pin resolved CFD simulation ever performed to the authors' knowledge. This represents a significant advancement in capability for the CFD of nuclear reactors, which will hopefully serve as an inspiration for further integrating high-fidelity numerical simulations in actual engineering designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Particle-Hole Asymmetric Ferromagnetism and Spin Textures in the Triangular Hubbard-Hofstadter Model

In a lattice model subject to a perpendicular magnetic field, when the lattice constant is comparable to the magnetic length, one enters the “Hofstadter regime,” where continuum Landau levels become fractal magnetic Bloch bands. Strong mixing between bands alters the nature of the resulting quantum phases compared to the continuum limit; lattice potential, magnetic field, and Coulomb interaction must be treated on equal footing. Using determinant quantum Monte Carlo and density matrix renormalization group techniques, we study this regime numerically in the context of the Hubbard-Hofstadter model on a triangular lattice. In the field-filling phase diagram, we find a broad wedge-shaped region of ferromagnetic ground states for filling factor ν ≤ 1 , bounded below by filling factor ν = 1 and bounded above by half filling the lowest Hofstadter subband. We observe signatures of SU(2) quantum Hall ferromagnetism at filling factors ν = 1 and ν = 3 . The phases near ν = 1 are particle-hole asymmetric, and we observe a rapid decrease in ground-state spin polarization consistent with the formation of skyrmions only on the electron doped side. At large fields, above the ferromagnetic wedge, we observe a low-spin metallic region with spin correlations peaked at small momenta. We argue that the phenomenology of this region likely results from exchange interaction mixing fractal Hofstadter subbands. The phase diagram derived beyond the continuum limit points to a rich landscape to explore interaction effects in magnetic Bloch bands. Published by the American Physical Society 2024

Ding, Jixun K.↗

DOE OSTI2022

Microgrid Assisted Design for Remote Areas

In this work, we present a three-stage multiobjective mixed-integer linear programming (MILP) for the optimal expansion planning and operation of isolated multienergy microgrids in remote areas. By selecting the optimal distributed generators (DGs) and energy storage systems (ESSs) mix selection, siting, sizing, and scheduling in the remote microgrid, the proposed model is targeted to minimize the annualized total cost of microgrids while enhancing the performance of the system, i.e., minimizing the voltage deviations and line power loss. To represent the electricity and heat flow between generation resources and various electrical, heating, and cooling loads in the isolated microgrid, linearized power flow, and heat flow constraints are employed in the proposed optimization model. The available capacity of DGs and ESSs are modeled as discrete constants instead of continuous variables for practical purpose. Numerical simulation results on a remote microgrid consisting of DGs, ESSs, and various loads validate the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Less can be more: Insights on the role of electrode microstructure in redox flow batteries from two-dimensional direct numerical simulations

Understanding how to structure a porous electrode to facilitate fluid, mass, and charge transport is key to enhancing the performance of electrochemical devices, such as fuel cells, electrolyzers, and redox flow batteries (RFBs). Here, using a parallel computational framework, direct numerical simulations are carried out on idealized porous electrode microstructures for RFBs. Strategies to improve an electrode design starting from a regular lattice are explored. By introducing vacancies in the ordered arrangement, it is possible to achieve higher voltage efficiency at a given current density, thanks to improved mixing of reactive species, despite reducing the total reactive surface. Careful engineering of the location of vacancies, resulting in a density gradient, outperforms disordered configurations. Our simulation framework is a new tool to explore transport phenomena in RFBs, and our findings suggest new ways to design performant electrodes.

25 ENERGY STORAGE↗

DOE OSTI2022

Nonlinear post-compression in multi-pass cells in the mid-IR region using bulk materials

We numerically investigate the regime of nonlinear pulse compression at mid-IR wavelengths in a multi-pass cell (MPC) containing a dielectric plate. This post-compression setup allows for ionization-free spectral broadening and self-compression while mitigating self-focusing effects. We find that self-compression occurs for a wide range of MPC and pulse parameters and derive scaling rules that enable its optimization. We also reveal the solitonic dynamics of the pulse propagation in the MPC and its limitations and show that spatiotemporal/spectral couplings can be mitigated for appropriately chosen parameters. In addition, we reveal the formation of spectral features akin to quasi-phase matched degenerate four-wave mixing. Finally, we present two case studies of self-compression at 3-μm and 6-μm wavelengths using pulse parameters compatible with driving high-field physics experiments. The simulations presented in this paper set a framework for future experimental work using few-cycle pulses at mid-IR wavelengths.

Carlson, D. (ORCID:0000000346836463)↗

DOE OSTI2022

Status on lattice calculations of the proton spin decomposition

Abstract Lattice calculations of the proton spin components is reviewed. The lattice results of the quark spin from the axial-vector current matrix element at ∼ 0.3−0.4 is smaller than those from the constituent quark models. This is largely due to the fact that the vacuum polarization contribution from the disconnected insertion is negative. Its connection with the anomalous Ward identity is clarified and verified numerically. This resolves the contentious issue in the “proton spin crisis.” The glue spin and angular momentum are found to be large and there is notable contribution from the quark orbital angular momentum. Renormalization, mixing, and normalization of the quark and glue angular momenta are discussed. With sufficient precision, they can be compared with more precise experimental measurements when the electron-ion collider facility is available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Unravelling the orbits of cluster galaxy populations according to their dominant gas ionization source

ABSTRACT We investigate the kinematical and dynamical properties of cluster galaxy populations classified according to their dominant source of gas ionization, namely: star-forming (SF) galaxies, optical active galactic nuclei (AGNs), mixed SF plus AGN ionization (transition objects, T), and quiescent (Q) galaxies. We stack 8892 member galaxies from 336 relaxed galaxy clusters to build an ensemble cluster and estimate the observed projected profiles of numerical density and velocity dispersion, $\sigma _P(R)$, of each galaxy population. The MAMPOSSt code and the Jeans equations inversion technique are used to constrain the velocity anisotropy profiles of the galaxy populations in both parametric and non-parametric ways. We find that Q (SF) galaxies display the lowest (highest) typical cluster-centric distances and velocity dispersion values. Transition galaxies are more concentrated and tend to exhibit lower velocity dispersion values than SF galaxies. Galaxies that host an optical AGN are as concentrated as Q galaxies but display velocity dispersion values similar to those of the SF population. MAMPOSSt is able to find equilibrium solutions that successfully recover the observed $\sigma _P(R)$ profile only for the Q, T, and AGN populations. We find that the orbits of all populations are consistent with isotropy in the inner regions, becoming increasingly radial with the distance from the cluster centre. These results suggest that Q galaxies are in equilibrium within their clusters, while SF galaxies have more recently arrived in the cluster environment. Finally, the T and AGN populations appear to be in an intermediate dynamical state between those of the SF and Q populations.

Valk, Greique A. (ORCID:0009000827731299)↗

DOE OSTI2026

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗

DOE OSTI2024

Implosion simulations and optimum times of radiographic Images in OMEGA high-density-carbon cylinder experiments

The Los Alamos Ignition Threshold team is planning direct-drive cylinder implosion experiments (CylDRT24B) at the OMEGA laser, scheduled for Thursday 23 May 2024, to investigate the effect of granular microstructure in high density carbon (HDC) shells. Granular microstructure is believed to play a crucial role in seeding perturbations, reducing compression, and enhancing mix in the implosion of HDC shells at the National Ignition Facility (NIF). But the ICF community lacks a detailed understanding of how grains induce these effects, owing partly to the practical difficulty of numerically simulating the behavior of tiny nanometer-scale granular structures in laboratory-scale experiments. Our planned experiments are aimed at acquiring data to help constrain simulations and calibrate reduced models.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

DOE OSTI2023

Record-breaking 2020 summer marine heatwaves in the western North Pacific

Record-breaking marine heatwaves (MHWs) occurred in the western North Pacific during the summer of 2020. These unprecedented MHWs were consistent with favorable large-scale conditions that are linked to an anomalous western North Pacific subtropical high (WNPSH), resulting mainly from sea surface temperature (SST) anomalies across the tropical oceans. In addition, a moderate La Niña-like pattern was also conducive to transporting warm seawater to the western North Pacific. Mixed-layer heat budgets suggest that surface heat flux contributed to the SST anomaly in the western subtropical Pacific. In contrast, oceanic heat advection dominanted in the South China Sea and the western equatorial Pacific. Numerical model experiments indicated that the tropical Indian Ocean SST anomalies were responsible for the enhanced WNPSH. The increased zonal SST gradient across the tropical Pacific also played an important role. Inter-ocean interactions can modulate climate variability through ocean-atmospheric coupling and deserve more attention when predicting MHWs within the context of global warming. Additionally, it is critical to consider MHWs as a powerful tool in detecting acute, intense thermal stress events in the coral bleaching pre-warning system.

Coral bleaching↗

DOE OSTI2026

A robust framework for frictional fault contact in geological formations using a stabilized augmented Lagrangian approach

Numerical simulations are essential to evaluate the performance and safety of engineered subsurface systems such as geological carbon storage sites, enhanced geothermal fields, and oil and gas reservoirs. A key challenge lies in accurately modeling the frictional contact behavior along fault surfaces. This problem involves inequality constraints that arise from the physics of frictional slip, requiring specialized numerical methods to handle the resulting highly nonlinear and path-dependent behavior. Here, in this work, we address this challenge using an Augmented Lagrangian Method (ALM) implemented via the Uzawa algorithm. The formulation employs mixed finite element spaces, combining low-order piecewise linear displacements within the 3D domain cells with piecewise constant tractions defined on the fault surfaces. Furthermore, to ensure stability and satisfy the inf-sup condition, the discrete displacement space is enriched with face bubble functions on both sides of the contact interfaces. This approach offers several advantages over other stabilization techniques that rely on additional terms, and it integrates naturally in the Uzawa framework.

58 GEOSCIENCES↗