SEARCH · Search NASA

Results for “dispersive equations”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1991

I-BIEM calculations of the frequency dispersion and AC current distribution at disk and ring-disk electrodes

The Iterative Boundary Integral Equation Method (I-BIEM) has been applied to the problem of frequency dispersion at a disk electrode in a finite geometry. The I-BIEM permits the direct evaluation of the AC potential (a complex variable) using complex boundary conditions. The point spacing was made highly nonuniform, to give extremely high resolution in those regions where the variables change most rapidly, i.e., in the vicinity of the edge of the disk. Results are analyzed with respect to IR correction, equipotential surfaces, and reference electrode placement. The current distribution is also examined for a ring-disk configuration, with the ring and the disk at the same AC potential. It is shown that the apparent impedance of the disk is inductive at higher frequencies. The results are compared to analytic calculations from the literature, and usually agree to better than 0.001%.

Boris D Cahan↗

NASA NTRS1991

I-BIEM calculations of the frequency dispersion and ac current distribution at disk and ring-disk electrodes

Cahan, Boris D.↗

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↗

NASA NTRS1994

Direct calculations of waves in fluid flows using high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.

Yu, Sheng-Tao↗

DOE OSTI2024

Dispersive nonreciprocity between a qubit and a cavity

The dispersive interaction between a qubit and a cavity is ubiquitous in circuit and cavity quantum electrodynamics. It describes the frequency shift of one quantum mode in response to excitations in the other and, in closed systems, is necessarily bidirectional, i.e., reciprocal. Here, we present an experimental study of a nonreciprocal dispersive-type interaction between a transmon qubit and a superconducting cavity, arising from a common coupling to dissipative intermediary modes with broken time reversal symmetry. We characterize the qubit-cavity dynamics, including asymmetric frequency pulls and photon shot noise dephasing, under varying degrees of nonreciprocity by tuning the magnetic field bias of a ferrite component in situ. We introduce a general master equation model for nonreciprocal interactions in the dispersive regime, providing a compact description of the observed qubit-cavity dynamics agnostic to the intermediary system. Our result provides an example of quantum nonreciprocal phenomena beyond the typical paradigms of non-Hermitian Hamiltonians and cascaded systems.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2020

Kinetic-energy dissipation and fluctuations in strongly damped heavy-ion collisions within the stochastic mean-field approach

Microscopic mean-field approaches have been successful in describing the most probable reaction outcomes in low-energy heavy-ion reactions. However, those approaches are known to severely underestimate dispersions of observables around the average values that has limited their applicability. Recently it has been shown that a quantal transport approach based on the stochastic mean-field (SMF) theory significantly improves the description, while its application has been limited so far to fragment mass and charge dispersions. In this work, we extend the quantal transport approach based on the SMF theory for relative kinetic energy dissipation and angular momentum transfer in low-energy heavy-ion reactions. Based on the SMF concept, analytical expressions are derived for the radial and tangential friction and associated diffusion coefficients. Those quantal transport coefficients are calculated microscopically in terms of single-particle orbitals within the time-dependent Hartree-Fock (TDHF) approach. As the first application of the proposed formalism, we consider the radial linear momentum dispersion, neglecting the coupling between radial and angular momenta. We analyze the total kinetic energy (TKE) distribution of binary reaction products in the 136 Xe + 208 Pb reaction at E c.m. =526 MeV and compare with experimental data. From time evolution of single-particle orbitals in TDHF, the radial diffusion coefficient is computed on a microscopic basis, while a phenomenological treatment is introduced for the radial friction coefficient. By solving the quantal diffusion equation for the radial linear momentum, the dispersion of the radial linear momentum is obtained, from which one can construct the TKE distribution. We find that the calculations provide a good description of the TKE distribution for strongly damped events with large energy losses, TKEL ≳ 150 MeV. However, the calculations underestimate the TKE distribution for smaller energy losses. Further studies are needed to improve the technical details of calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Quantum space, quantum time, and relativistic quantum mechanics

We treat space and time as bona fide quantum degrees of freedom on an equal footing in Hilbert space. Motivated by considerations in quantum gravity, we focus on a paradigm dealing with linear, first-order Hamiltonian and momentum constraints that lead to emergent features of temporal and spatial translations. Unlike the conventional treatment, we show that Klein-Gordon and Dirac equations in relativistic quantum mechanics can be unified in our paradigm by applying relativistic dispersion relations to eigenvalues rather than treating them as operator-valued equations. With time and space being treated on an equal footing in Hilbert space, we show symmetry transformations to be implemented by unitary basis changes in Hilbert space, giving them a stronger quantum mechanical footing. Global symmetries, such as Lorentz transformations, modify the decomposition of Hilbert space; and local symmetries, such as U(1) gauge symmetry are diagonal in coordinate basis and do not alter the decomposition of Hilbert space. Here, we briefly discuss extensions of this paradigm to quantum field theory and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Spatiotemporal predictions of toxic urban plumes using deep learning

Industrial accidents, chemical spills, and structural fires can release large amounts of harmful materials that disperse into urban atmospheres and impact populated areas. Computer models are typically used to predict the transport of toxic plumes by solving fluid dynamical equations. However, these models can be computationally expensive due to the need for many grid cells to simulate turbulent flow and resolve individual buildings and streets. In emergency response situations, alternative methods are needed that can run quickly and adequately capture important spatiotemporal features. Here, we present a novel deep learning model called ST-GasNet inspired by the mathematical equations that govern the behavior of plumes as they disperse through the atmosphere. ST-GasNet learns the spatiotemporal dependencies from a limited set of temporal sequences of ground-level toxic urban plumes generated by a high-resolution large eddy simulation model. On independent sequences, ST-GasNet accurately predicts the late-time spatiotemporal evolution, given the early-time behavior as an input, even when a building splits a large plume into smaller plumes. By incorporating large-scale wind boundary condition information, ST-GasNet achieves a prediction accuracy of at least 90% on test data for the entire prediction period.

Civil and Environmental Engineering↗

NASA NTRS1978

Segregation of acid plume pixels from background water pixels, signatures of background water and dispersed acid plumes, and implications for calculation of iron concentration in dense plumes

Two files of data, obtained with a modular multiband scanner, for an acid waste dump into ocean water, were analyzed intensively. Signatures were derived for background water at different levels of effective sunlight intensity, and for different iron concentrations in the dispersed plume from the dump. The effect of increased sunlight intensity on the calculated iron concentration was found to be relatively important at low iron concentrations and relatively unimportant at high values of iron concentration in dispersed plumes. It was concluded that the basic equation for iron concentration is not applicable to dense plumes, particularly because lower values are indicated at the very core of the plume, than in the surrounding sheath, whereas radiances increase consistently from background water to dispersed plume to inner sheath to innermost core. It was likewise concluded that in the dense plume the iron concentration would probably best be measured by the higher wave length radiances, although the suitable relationship remains unknown.

Bahn, G. S.↗

DOE OSTI2022

Spontaneous Exciton Collapse in a Strongly Flattened Ellipsoidal InSb Quantum Dot

Abstract Electronic and excitonic states in an InSb strongly flattened ellipsoidal quantum dot (QD) with complicated dispersion law are theoretically investigated within the framework of the geometric adiabatic approximation in the strong, intermediate, and weak quantum confinement regimes. For the lower levels of the spectrum, the square root dependence of energy on QD sizes is revealed in the case of Kane’s dispersion law. The obtained results are compared to the case of a parabolic (standard) dispersion law of charge carriers. The possibility of the accidental exciton instability is revealed for the intermediate quantum confinement regime. For the weak quantum confinement regime, the motion of the exciton's center-of-gravity is quantized, which leads to the appearance of additional Coulomb-like sub-levels. It is revealed that in the case of the Kane dispersion law, the Coulomb levels shift into the depth of the forbidden band gap, moving away from the quantum confined level, whereas in the case of the parabolic dispersion law, the opposite picture is observed. The corresponding selection rules of quantum transitions for the interband absorption of light are obtained. New selection rules of quantum transitions between levels conditioned by 2D exciton center of mass vertical motion quantization in a QD are revealed. The absorption threshold behavior characteristics depending on the QDs geometrical sizes are also revealed.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗

NASA NTRS1993

Direct calculations of waves in fluid flows using a high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expressions for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations, and the numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the CFL limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about 6 to 8 grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 to 300 time steps. Three numerical examples of waves in compressible flow are included.

Yu, Sheng-Tao↗

DOE OSTI2023

A numerical extension of the spatially-filtered Euler equations for contact discontinuities

Solving the spatially-filtered Euler equations with a kinetic energy-preserving centered discretization is an efficient framework for simulating shock-dominated flows. Here, the present work describes a numerical treatment of the spatially-filtered Euler equations to minimize oscillations when simulating contact discontinuities. To counteract dispersive errors inherent to centered schemes, a WENO-like correction term is applied to the enthalpy transport in the energy equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

NASA NTRS1996

[Kelvin-Helmholtz instability in protostellar jets]

NASA grant NAG 5 2866, funded by the Astrophysics Theory Program, enabled the study the Kelvin-Helmholtz instability in protostellar jets. In collaboration with co-investigator Philip Hardee, the PI derived the analytic dispersion relation for the instability in including a cooling term in the energy equation which was modeled as one of two different power laws. Numerical solutions to this dispersion relation over a wide range of perturbation frequencies, and for a variety of parameter values characterizing the jet (such as Mach number, and density ratio) were found It was found that the growth rates and wavelengths associated with unstable roots of the dispersion relation in cooling jets are significantly different than those associated with adiabatic jets, which have been studied previously. In collaboration with graduate student Jianjun Xu (funded as a research associate under this grant), hydrodynamical simulations were used to follow the growth of the instability into the nonlinear regime. It was found that asymmetric surface waves lead to large amplitude, sinusoidal distortions of the jet, and ultimately to disruption Asymmetric body waves, on the other hand, result in the formation of shocks in the jet beam in the nonlinear regime. In cooling jets, these shocks lead to the formation of dense knots and filaments of gas within the jet. For sufficiently high perturbation frequencies, however, the jet cannot respond and it remains symmetric. Applying these results to observed systems, such as the Herbig-Haro jets HH34, HH111 and HH47 which have been observed with the Hubble Space Telescope, we predicted that some of the asymmetric structures observed in these systems could be attributed to the K-H modes, but that perturbations on timescales associated with the inner disk (about 1 year) would be too rapid to cause disruption. Moreover, it was discovered that weak shock 'spurs' in the ambient gas produced by ripples in the jet surface due to nonlinear, modes of surface and/or body waves could accelerate the ambient gas to low velocity. This latter effect represents a new mechanism by which supersonic jets can accelerate low velocity outflows.

Stone, James↗

DOE OSTI2020

A class of solutions of the two-dimensional Toda lattice equation

In this study, a method is proposed to systematically generate solutions of the two-dimensional Toda lattice equation in terms of previously known solutions &straightphi;(x,y) of the two-dimensional Laplace's equation. The solitonic solution of Nakamura's [J. Phys. Soc. Jpn. 52 (1983) 380] is shown to correspond to one particular choice of &straightphi;(x,y).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Development and assessment of hierarchical multi-reward reinforcement learning based potential for silicene with state-of-the-art models

We develop a new interatomic force field for Silicene, a 2D material with a buckled hexagonal lattice structure with high polymorphism. We introduce new parameterizations of a Tersoff model using a hierarchical multi-reward reinforcement learning (RL) methodology coupled with a continuous Monte Carlo Tree Search optimization. Our model significantly outperforms existing methods by enhancing the accuracy of predictions for the structural and thermodynamic properties of seven silicene polymorphs-including structure, energy, equation of state, elasticity, and phonon dispersion-when compared to established models. We further make a comprehensive comparison of the various models in predicting the mechanical and thermal properties of silicene. We trace the origin of the improved performance to the description of the angular dependence in the bond-order term, suggesting that modifying the angular terms in short-range models is essential to capture the structural diversity in low dimensional systems.

2D materials↗

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)↗