Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical dispersion”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Investigation of finite-volume methods to capture shocks and turbulence spectra in compressible flows

The aim of the present paper is to provide a comparison between several finite-volume methods of different numerical accuracy: the second-order Godunov method with PPM interpolation and the high-order finite-volume WENO method. In this work, the results show that while on a smooth problem the high-order method performs better than the second-order one, when the solution contains a shock all the methods collapse to first-order accuracy. In the context of the decay of compressible homogeneous isotropic turbulence with shocklets, the actual overall order of accuracy of the methods reduces to second-order, despite the use of fifth-order reconstruction schemes at cell interfaces. Most important, results in terms of turbulent spectra are similar regardless of the numerical methods employed, except that the PPM method fails to provide an accurate representation in the high-frequency range of the spectra. It is found that this specific issue comes from the slope-limiting procedure and a novel hybrid PPM/WENO method is developed that has the ability to capture the turbulent spectra with the accuracy of a high-order method, but at the cost of the second-order Godunov method. Overall, it is shown that virtually the same physical solution can be obtained much faster by refining a simulation with the second-order method and carefully chosen numerical procedures, rather than running a coarse high-order simulation. Our results demonstrate the importance of evaluating the accuracy of a numerical method in terms of its actual spectral dissipation and dispersion properties on mixed smooth/shock cases, rather than by the theoretical formal order of convergence rate.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Vibrational Fermi Resonance in Atomically Thin Black Phosphorus

Fermi resonance is a phenomenon involving the hybridization of two coincidentally quasi-degenerate states that is observed in the vibrational or electronic spectra of molecules. Despite numerous examples in molecular systems, vibrational Fermi resonances in dispersive semiconducting systems remain largely unexplored due to the rarity of occurrence. Here we report a vibrational Fermi resonance in atomically thin black phosphorus. The Fermi resonance arises via anharmonic mixing of a fundamental Raman mode and a Davydov component of an infrared mode, leading to a doublet with mixed character. The extent of Fermi coupling can be modulated by the application of external biaxial strain. The consequences of Fermi hybridization are revealed by electronic resonance effects in the thickness-dependent and excitation-wavelength-dependent Raman spectrum, which is predicted by ab initio hybrid functional simulations including excitonic interactions. This work reveals new insight into electron–phonon coupling in black phosphorus and demonstrates a novel method for modulating Fermi resonances in 2D semiconductors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Shock formation in flowing plasmas by temporally and spatially smoothed laser beams

The cumulative impact of multiple laser speckles on a supersonic plasma flow across optically smoothed laser beams is investigated. The bending of laser beams caused by ponderomotive laser–plasma coupling, together with flow, leads to plasma a momentum-conserving response that results in a deceleration of the flow. Once the flow velocity decreases to a subsonic level, the action of the laser beams can generate a shock within the plasma. This scenario has been predicted theoretically and confirmed by hydrodynamic simulations. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size, and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging only in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser–plasma coupling in the inertial confinement fusion experiments are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗
DOE OSTI2025

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗
DOE OSTI2021

X-ray free electron laser heating of water and gold at high static pressure

Probing of reactive materials such as H 2 O ices and fluids at the high pressures and temperatures of planetary interiors is limited by unwanted chemical reactions and confinement failure. Faster experiments can mitigate such issues, but the common approach of adiabatic compression limits the conditions achieved. This study demonstrates a fast experimental strategy for the creation and probing of selected extreme states using static compression coupled with ultrafast X-ray laser heating. Indirect X-ray heating of H 2 O through the use of a gold absorber is evidenced by sample melting inferred from textural changes in the H 2 O diffraction lines and inter-dispersion of gold and H 2 O melts. Coupled with numerical analysis of femtosecond energy absorption, thermal equilibration, and heat transfer, all evidence indicates that temperatures in excess of an electron volt have been reached in the H 2 O at high pressure. Even after repeated heating, samples stayed chemically unchanged from the starting material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Latest Updates on JSPEC - An IBS and Electron Cooling Simulation Program

JSPEC (JLab Simulation Package on Electron Cooling) is an open-source C++ program developed at Jefferson Lab, which simulates the evolution of the ion beam under the influence of both IBS and electron cooling effects. In this paper, we will report the latest updates to JSPEC. Firstly, we have added theoretical and numerical models that simulate the effect of the electron beam dispersion on non-magnetized cooling. Secondly, the cooler can now be treated as an element with length, rather than a thin lens. This change will impact the modelling of the ions and the electrons in the cooling rate calculation for both magnetized and non-magnetized cooling. Numeri-cal results will be provided to demonstrate the perfor-mance of the new models.

Zhang, H.↗
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗
DOE OSTI2019

Effects of photoevaporation on protoplanetary disc ‘isochrones’

ABSTRACT Protoplanetary discs are the site of star and planet formation, and their evolution and consequent dispersal deeply affect the formation of planetary systems. In the standard scenario they evolve on time-scales ∼Myr due to the viscous transport of angular momentum. The analytical self-similar solution for their evolution predicts specific disc isochrones in the accretion rate–disc mass plane. However, photoevaporation by radiation emitted by the central star is likely to dominate the gas disc dispersal of the innermost region, introducing another (shorter) time-scale for this process. In this paper, we include the effect of internal (X and EUV) photoevaporation on the disc evolution, finding numerical solutions for a population of protoplanetary discs. Our models naturally reproduce the expected quick dispersal of the inner region of discs when their accretion rates match the rate of photoevaporative mass loss, in line with previous studies. We find that photoevaporation preferentially removes the lightest discs in the sample. The net result is that, counter-intuitively, photoevaporation increases the average disc mass in the sample, by dispersing the lightest discs. At the same time, photoevaporation also reduces the mass accretion rate by cutting the supply of material from the outer to the inner disc. In a purely viscous framework, this would be interpreted as the result of a longer viscous evolution, leading to an overestimate of the disc age. Our results thus show that photoevaporation is a necessary ingredient to include when interpreting observations of large disc samples with measured mass accretion rates and disc masses. Photoevaporation leaves a characteristic imprint on the shape of the isochrone. Accurate data of the accretion rate–disc mass plane in the low disc mass region therefore give clues on the typical photoevaporation rate.

Somigliana, Alice↗
DOE OSTI2022

Dispersion interactions in proposed covalent superhydride superconductors

Recent developments in high temperature superconducting materials under high pressure have made numerical evaluation of the superconducting transition temperature (T c ) of predicted materials critically important as a means to identification. Existing methods of calculating T c often do not agree with each other or experiment, often due to the large number of complex factors that contribute to this property; among them is the neglect of dispersion interactions in commonly used density functionals. We evaluate the effect of including dispersion interactions on the predicted superconducting properties of two examples of the covalent superhydride class of very high-T c superconducting materials. In both cases, dispersion is found to have sizable effects, increasing the electron-phonon coupling as compared to the reference case of elemental niobium. A detailed investigation traces the origin of this effect in a 270 GPa $\mathrm {R3m CSH_7}$ structure to structural distortions driven by long-range electron-phonon interactions rather than novel bonding networks.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gallium-doped zinc oxide: nonlinear reflection and transmission measurements and modeling in the ENZ region

Abstract Strong nonlinear materials have been sought after for decades for applications in telecommunications, sensing, and quantum optics. Gallium-doped zinc oxide is a II–VI transparent conducting oxide that shows promising nonlinearities similar to indium tin oxide and aluminum-doped zinc oxide for the telecommunications band. Here we explore its nonlinearities in the epsilon near zero (ENZ) region and show n 2,eff values on the order of 4.5 × 10 −3 cm 2 GW −1 for IR pumping on 200–300 nm thin films. Measuring nonlinear changes in transmission and reflection with a white light source probe in the near-IR while exciting in the near-IR provides data in both time and wavelength. Three films varying in thickness, optical loss, and ENZ crossover wavelength are numerically modeled and compared to experimental data showing agreement for both dispersion and temporal relaxation. In addition, we discuss optimal excitation and probing wavelengths occur around ENZ for thick films but are red-shifted for thin films where our model provides an additional degree of freedom to explore. Obtaining accurate nonlinear measurements is a difficult and time-consuming task where our method in this paper provides experimental and modeled data to the community for an ENZ material of interest.

Ball, Adam (ORCID:0000000318167813)↗
DOE OSTI2024

Toward a unified description of hadron scattering at all energies

The construction of general amplitudes satisfying symmetries and S -matrix constraints has been the primary tool in studying the spectrum of hadrons for over half a century. In this work, we present a new parametrization, which can fulfill many expectations of S -matrix and Regge theory and connects the essential physics of hadron scattering in the resonance region and in asymptotic limits. In this construction, dynamical information is entirely contained in Regge trajectories that generalize resonance poles in the complex energy plane to moving poles in the angular momentum plane. We highlight the salient features of the model, compare with existing literature on dispersive and dual amplitudes, and benchmark the formalism with an initial numerical application to the ρ and σ / f 0 ( 500 ) mesons in π π scattering. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Spatiotemporal spread of perturbations in power-law models at low temperatures: Exact results for classical out-of-time-order correlators

Here, we present exact results for the classical version of the out-of-time-order commutator (OTOC) for a family of power-law models consisting of N particles in one dimension and confined by an external harmonic potential. These particles are interacting via power-law interaction of the form ∝ $Σ_{i,j=(i≠j)}^{N}$ |x i -x j | -k for every k>1 where x i is the position of the i th particle. We present numerical results for the OTOC for finite N at low temperatures and short enough times so that the system is well approximated by the linearized dynamics around the many-body ground state. In the large-N limit, we compute the ground-state dispersion relation in the absence of external harmonic potential exactly and use it to arrive at analytical results for OTOC. We find excellent agreement between our analytical results and the numerics. We further obtain analytical results in the limit where only linear and leading nonlinear (in momentum) terms in the dispersion relation are included. The resulting OTOC is in agreement with numerics in the vicinity of the edge of the “light cone.” We find remarkably distinct features in OTOC below and above k=3 in terms of going from non-Airy behavior (1 < k < 3) to an Airy universality class (k>3). We present certain additional rich features for the case k=2 that stem from the underlying integrability of the Calogero-Moser model. We present a field theory approach that also assists in understanding certain aspects of OTOC such as the sound speed. Our findings are a step forward towards a more general understanding of the spatiotemporal spread of perturbations in long-range interacting systems.

1-dimensional systems↗
DOE OSTI2022

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗
DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗