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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Identification and Characterization of Persistent Cold Pool Events from Temperature and Wind Profilers in the Columbia River Basin

Cold pool events occur when deep layers of stable, cold air remain trapped in a valley or basin for multiple days, without mixing out from daytime heating. With large impacts on air quality, freezing events, and es-pecially on wind energy production, they are often poorly forecast by modern mesoscale numerical weather prediction (NWP) models. Understanding the characteristics of cold pools is, therefore, important to provide more accurate forecasts. This study analyzes cold pool characteristics with data collected during the Second Wind Forecast Improvement Project (WFIP2), which took place in the Columbia River basin and Gorge of Oregon and Washington from fall 2015 until spring 2017. A subset of the instrumentation included three microwave radiometer pro?lers, six radar wind pro?lers with radio acoustic sounding systems, and seven sodars, which together provided seven sites with collocated vertical pro?les of temperature, humidity, wind speed, and wind direction. Using these collocated observations, we developed a set of criteria to determine if a cold pool was present based on stability, wind speed, direction, and temporal continuity, and then developed an automated algorithm based on these criteria to identify all cold pool events over the 18 months of the ?eld project. Characteristics of these events are described, including statistics of the wind speed distributions and pro?les, stability conditions, cold pool depths, and descent rates of the cold pool top. The goal of this study is a better understanding of these characteristics and their processes to ultimately lead to improved physical parameterizations in NWP models, and consequently improve forecasts of cold pool events in the study region as well at other locations that experiences similar events.

McCaffrey, Katherine L.↗
DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗
DOE OSTI2021

Observed and Simulated Variability of Droplet Spectral Dispersion in Convective Clouds Over the Amazon

Abstract In this study, the variability of the spectral dispersion of droplet size distributions (DSDs) in convective clouds is investigated. Analyses are based on aircraft measurements of growing cumuli near the Amazon basin, and on numerical simulations of an idealized ice‐free cumulus. In cleaner clouds, the relative dispersion , defined as the ratio of the standard deviation to the mean value of the droplet diameter, is negatively correlated with the ratio of the cloud water content ( ) to the adiabatic liquid water content ( ), while no strong correlation between and is seen in polluted clouds. Bin microphysics numerical simulations suggest that these contrasting behaviors are associated with the effect of collision‐coalescence in cleaner clouds, and secondary droplet activation in polluted clouds, in addition to the turbulent mixing of parcels that experienced different paths within the cloud. Collision‐coalescence simultaneously broadens the DSDs and decreases , explaining the inverse relationship between and in cleaner clouds. Secondary droplet activation broadens the DSDs but has little direct impact on . The combination of a rather modest DSD broadening due to weak collision‐coalescence with enhanced droplet activation in both diluted and highly undiluted cloud regions may contribute to maintain a relatively uniform within polluted clouds. These findings can be useful for parameterizing the shape parameter ( ) of gamma DSDs in bulk microphysics cloud‐resolving models. It is shown that emulating the observed relationship improves the estimation of the collision‐coalescence rate in bulk microphysics simulations compared to the bin simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multimode description of self-mode locking in a single-section quantum-dot laser

This paper describes a theory for mode locking and frequency comb generation by four-wave mixing in a semiconductor quantum-dot active medium. The derivation uses a multimode semiclassical laser theory that accounts for fast carrier collisions within an inhomogeneous distribution of quantum dots. Numerical simulations are presented to illustrate the role of active medium nonlinearities in mode competition, gain saturation, carrier-induced refractive index and creation of combination tones that lead to locking of beat frequencies among lasing modes in the presence of cavity material dispersion.

Chow, Weng W.↗
DOE OSTI2021

A Hybrid Energy System Workflow for Energy Portfolio Optimization

This manuscript develops a workflow, driven by data analytics algorithms, to support the optimization of the economic performance of an Integrated Energy System. The goal is to determine the optimum mix of capacities from a set of different energy producers (e.g., nuclear, gas, wind and solar). A stochastic-based optimizer is employed, based on Gaussian Process Modeling, which requires numerous samples for its training. Each sample represents a time series describing the demand, load, or other operational and economic profiles for various types of energy producers. These samples are synthetically generated using a reduced order modeling algorithm that reads a limited set of historical data, such as demand and load data from past years. Numerous data analysis methods are employed to construct the reduced order models, including, for example, the Auto Regressive Moving Average, Fourier series decomposition, and the peak detection algorithm. All these algorithms are designed to detrend the data and extract features that can be employed to generate synthetic time histories that preserve the statistical properties of the original limited historical data. The optimization cost function is based on an economic model that assesses the effective cost of energy based on two figures of merit: the specific cash flow stream for each energy producer and the total Net Present Value. An initial guess for the optimal capacities is obtained using the screening curve method. The results of the Gaussian Process model-based optimization are assessed using an exhaustive Monte Carlo search, with the results indicating reasonable optimization results. The workflow has been implemented inside the Idaho National Laboratory’s Risk Analysis and Virtual Environment (RAVEN) framework. The main contribution of this study addresses several challenges in the current optimization methods of the energy portfolios in IES: First, the feasibility of generating the synthetic time series of the periodic peak data; Second, the computational burden of the conventional stochastic optimization of the energy portfolio, associated with the need for repeated executions of system models; Third, the inadequacies of previous studies in terms of the comparisons of the impact of the economic parameters. The proposed workflow can provide a scientifically defendable strategy to support decision-making in the electricity market and to help energy distributors develop a better understanding of the performance of integrated energy systems.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

Mesoscopic structure of mixed type domain walls in multiaxial ferroelectrics

The structure of a 180° uncharged rotational domain wall in a multiaxial ferroelectric film is studied in the framework of an analytical Landau-Ginzburg-Devonshire (LGD) approach. Finite element modeling (FEM) is used to solve numerically the system of the coupled nonlinear Euler-Lagrange (EL) second-order differential equations for two components of polarization. We show that the structure of the domain wall and corresponding metastable or stable phase of the film are controlled by a single parameter—the dimensionless ferroelectric anisotropy factor μ. We fitted the static profile of a solitary domain wall, calculated by FEM, with kinklike functions for polarization components, and extracted the five μ-dependent parameters from the fitting to FEM curves. The surprisingly high accuracy of the fitting results for two polarization components in the entire μ range allows us to conclude that the analytical functions, which are trial functions in the direct variational method, can be treated as a high-accuracy variational solution of the static EL equations. We further derive exact two-component analytical solutions of the static EL equations for a polydomain 180° domain structure in a multiaxial ferroelectric film. Using these, we derive analytical expressions for the system free energy and analyze its dependence on the film thickness and boundary conditions at the film surfaces. The single-domain state is ground for zero polarization derivative at the surfaces, while the polydomain states minimize the system energy for zero polarization at the surfaces. Counterintuitively, the energy of the polydomain states split into two levels, “0” and “1,” for zero polarization at the surfaces, and each of the levels contains a large number of close-energy sublevels, whose structure is characterized by a different number and type of domain walls. The analytical solutions can become a useful tool for Bayesian analysis of high-resolution scanning transmission electron microscopy images in ferroelectric films.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ice rule breakdown and frustrated antiferrotoroidicity in an artificial colloidal Cairo ice

Abstract We combine experiments and numerical simulations to investigate the low energy states and the emergence of topological defects in an artificial colloidal ice in the Cairo geometry. This type of geometry is characterized by a mixed coordination ( z ), with coexistence of both z = 3 and z = 4 vertices. We realize this particle ice by confining field tunable paramagnetic colloidal particles within a lattice of topographic double wells at a one to one filling using optical tweezers. By raising the interaction strength via an applied magnetic field, we find that the ice rule breaks down, and positive monopoles with charge q = + 2 accumulate in the z = 4 vertices and are screened by negative ones ( q = − 1 ) in the z = 3. The resulting, strongly coupled state remains disordered. Further, via analysis of the mean chirality associated to each pentagonal plaquette, we find that the disordered ensemble for this geometry is massively degenerate and it corresponds to a frustrated antiferrotoroid.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Nonlinear quantum behavior of ultrashort-pulse optical parametric oscillators

The quantum features of ultrashort-pulse optical parametric oscillators (OPOs) are investigated theoretically in the nonlinear regime near and above threshold. Viewing the pulsed OPO as a multimode open quantum system, we rigorously derive a general input-output model that features nonlinear coupling among many cavity (i.e., system) signal modes and a broadband single-pass (i.e., reservoir) pump field. Under appropriate assumptions, our model produces a Lindblad master equation with multimode nonlinear Lindblad operators describing two-photon dissipation and a multimode four-wave-mixing Hamiltonian describing a broadband, dispersive optical cascade, which we show is required to preserve causality. To simplify the multimode complexity of the model, we employ a supermode decomposition to perform numerical simulations in the regime where the pulsed supermodes experience strong single-photon nonlinearity. We find that the quantum nonlinear dynamics induces pump depletion as well as corrections to the below-threshold squeezing spectrum predicted by linearized models. We also observe the formation of non-Gaussian states with Wigner-function negativity and show that the multimode interactions with the pump, both dissipative and dispersive, can act as effective decoherence channels. Finally, we briefly discuss some experimental considerations for potentially observing such quantum nonlinear phenomena with ultrashort-pulse OPOs on nonlinear nanophotonic platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Entanglement and seniority

In quantum mechanics entanglement is the most striking phenomenon which has no counterpart in classical systems. Although, different approaches have already been developed to study correlations in the case of indistinguishable particles, the exploration of the so-called mode entanglement is still in its initial stage in nuclear physics. Study of mode-entanglement in the seniority model, derivation of analytic formulas for the one-body reduced density matrix of states with seniority $ν$ = 0, 1, 2, and $ν$ = 3, and also determination of the particle number dependence of the one-body reduced density matrix for arbitrary seniority. In addition, comparison of the predictions of the seniority model with the results of the shell model to gain insight into the structure and correlations of the ground and lowest yrast states. In the seniority model the analytic results are calculated using the quasi-spin formalism. The numerical shell model calculations are carried out in the standard shell model framework using an inert core. The applied realistic effective interactions are derived earlier using the $G$-matrix formalism. The density matrix renormalization group method is also applied in order to directly calculate the mode entropies. In the seniority model simple analytical expressions are given for the mode entropies. The peculiar behavior of the half-filled shells and the seniority-zero states are revealed. Numerical results are presented for the lightest stable calcium isotopes and for 94 Ru nucleus. For 94 Ru, the seniority model accounts for the 0$g$ 9/2 mode entropies, but seniority mixing is important for certain yrast states. Furthermore, interaction induced quantum fluctuations decrease the occupation of the 0$f$ 5/2 , 1$p$ 3/2 and 1$p$ 1/2 shells, and amount in non-negligible mode entropies on these shells, too, clearly outside the scope of the simple (0$g$ 9/2 ) 4 seniority model. The 0$f$ 7/2 shell-based seniority model is more accurate for the Ca isotopes, but seniority mixing is substantial for some 44 Ca yrast states, too. Mode and one-body entanglement entropies are useful tools to investigate the structure of quantum correlations in nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

The role of equatorial waves in triggering precipitation extremes in the Maritime Continent

This review offers a comprehensive analysis of convectively coupled equatorial waves (CCEWs) and their pivotal role in driving precipitation extremes across the Maritime Continent. It examines the current understanding of CCEWs, evaluates the performance of numerical models and forecasting techniques in predicting these phenomena, and pinpoints critical areas for improvement. The discussion centers on three key types of equatorial waves: equatorial Rossby waves, Kelvin waves, and mixed Rossby–gravity waves. By connecting scientific insights with practical forecasting applications, the review sheds light on the challenges of predicting these waves while identifying opportunities to advance both fundamental knowledge and forecasting accuracy. Designed as an educational resource, it targets operational forecasting centers, meteorologists, and researchers, aiming to enhance the prediction of extreme weather events in the region.

Maritime Continent↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
DOE OSTI2023

Convectively coupled Rossby–Gravity waves in a field campaign: How they are captured in reanalysis products

Abstract Convectively coupled equatorial waves are a significant source of atmospheric variability in the tropics. Current numerical models continue to struggle in simulating the coupled diabatic heating fields that are responsible for the development and maintenance of these waves. This study investigates how the diabatic fields associated with Mixed Rossby–Gravity waves (MRGs) are represented in four reanalysis products by using a unique observational dataset from the TRMM‐KWAJEX (Tropical Rainfall Measuring Mission‐Kwajalein Experiment) field campaign. These reanalyses include ERA5, Japanese 55‐year Reanalysis (JRA‐55), Climate Forecast System Reanalysis (CFSR), and Modern‐Era Retrospective Analysis for Research and Applications (MERRA). We found that all four reanalyses captured the MRG structures in winds and temperature, and to a lesser degree in the humidity field except in the boundary layer. However, only the ERA5 and MERRA reanalyses captured the gradual rise and succession of the diabatic heating from boundary layer turbulence, shallow convection, cumulus congestus, and deep convection within the waves. ERA5 is the only product that also captured the gradual rise of the subgrid‐scale vertical transport of moist static energy. All reanalysis products underestimated the diabatic heating from cumulus congestus. Results provide observational basis on what aspects of MRG can be trusted and what cannot in the reanalysis products.

Wang, Xiaocong↗
DOE OSTI2025

Large eddy simulation of wind farm performance in horizontally and vertically staggered layouts

This numerical investigation employs Large Eddy Simulation (LES) coupled with Actuator Disk Model (ADM) to evaluate wind farm layout optimization strategies. The study presents a systematic analysis of aligned, horizontal staggering, vertical staggering, and mixed (combination of horizontal and vertical) staggering configurations, aiming to establish optimal design parameters for enhanced power production. The investigation examines key performance metrics including mean velocity distributions, turbulence intensity characteristics, and power generation efficiency. Results demonstrate better performance of both horizontal and vertical staggering patterns compared to conventional aligned configurations, with horizontal staggering exhibiting notably higher power output than vertical arrangements. Our findings also suggest that mixed configurations, incorporating both horizontal and vertical staggering, can offer optimal performance characteristics. As a result, this research advances the understanding of wake interactions in complex wind farm layouts and provides design guidelines for maximizing wind farm power generation efficiency through strategic turbine positioning.

17 WIND ENERGY↗
DOE OSTI2021

Ethane Dehydrogenation on Single and Dual Centers of Ga-modified γ-Al 2 O 3

Density-functional theory calculations and microkinetic analysis are used to investigate the efficacy of Ga-modified γ-Al 2 O 3 (110) surfaces for the catalytic dehydrogenation of ethane and elucidate the synergy between Ga and Al sites. The model surfaces are modified by either Ga grafting or doping. We consider and analyze numerous active sites and rank them using microkinetic analysis. The kinetic parameters obtained from microkinetic modelling are compared with experimental values for ethane dehydrogenation over Ga 2 O 3 -Al 2 O 3 mixed oxides prepared by co-precipitation. The dominant reaction pathway proceeds via heterolytic C-H bond dissociation to a surface proton and a metal-carbanion intermediate that undergoes β-hydride elimination. We find that grafted Ga sites are catalytically inactive. In contrast, Ga-doped sites exhibit fivefold enhancement in catalytic activity when compared to the sites on pristine Al 2 O 3 , owed to the synergy between neighboring AlIII and GaIV sites. Furthermore, we model and investigate the effect of surface hydroxylation, demonstrate how surface water interferes with the aforementioned synergy between AlIII and GaIV sites and discuss the implications for the catalytic activity of the modified surfaces. Here, an increase in the partial pressure of H 2 O significantly increases the apparent activation energies of dehydrogenation and interestingly changes the most active site.

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

Scaling of the mean transverse flow and Reynolds shear stress in turbulent plane jet

Here, proper scaling for the mean transverse flow and Reynolds shear stress in a turbulent plane jet is determined using a scaling patch approach. By seeking an admissible scaling, a key concept in the scaling patch approach, for the mean continuity equation, a proper scale for the mean transverse flow in a turbulent plane jet is found as V ref = –δdU ctr /dx, where δ is the jet half width and dU ctr /dx is the decay rate of the mean axial velocity at the jet centerline. By seeking an admissible scaling for the mean axial momentum equation, a proper scale for the kinematic Reynolds shear stress is found as R uv,ref =U ctr V ref , which is a mix of the velocity scales in the axial and transverse directions. Approximation functions for the scaled mean transverse flow and Reynolds shear stress are developed and found to agree well with experimental and numerical data. Similarities and differences between the scales of the mean transverse flow and Reynolds shear stress in turbulent plane jets and zero-pressure-gradient turbulent boundary layer flows are clarified.

42 ENGINEERING↗