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

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A Comparison of Adverse Meteorology Inputs for Emergency Management and Safety Basis Dispersion Modeling Applications at SRS

The MELCOR Accident Consequence Calculation System (MACCS), used for safety basis accident analyses at SRS, implements adjustment factors to consider the increased vertical turbulence associated with local surface roughness and the increased lateral turbulence associated with time-based plume meander. The two models used for SRS emergency planning radiological consequence assessment, i.e., Hotspot and Puff/Plume, currently do not have an option to input or apply diffusion adjustments for meander or surface roughness. Adverse meteorology conditions at SRS have been determined as E stability class and 1.3 meters per second wind speed for ground level releases. This is the input used for emergency planning models. Contrastingly safety basis analyses use a different input, representing adverse meteorology values as F stability class and wind speeds of 1.3 meters per second (for a ground level release). This results in an apparent discrepancy between the two approaches. In this report we document the differing modeling assumptions, discuss the logic behind each modeling approach used to represent local roughness and plume meander conditions, and numerically compare dispersion coefficients and concentration estimates obtained from each set of model inputs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A cancellation problem in hybrid particle-in-cell schemes due to finite particle size

The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the different spatial discretizations of ions as finite-spatial-size particles and electrons as a grid-based fluid can lead to significant numerical wave dispersion errors in the long wavelength limit (kd i «1, where k is the wavenumber and di is the ion skin-depth). The problem occurs when high-order particle-grid interpolations, or grid-based smoothing, spreads the electric field experienced by the ions across multiple spatial cells and leads to inexact cancellation of electric field terms in the total (ion + electron) momentum equation. Practical requirements on the mesh spacing Δx/d i are suggested to bound these errors from above. The accuracy impact of not respecting these resolution constraints is shown for a non-linear shock problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Infrared Near-Field Spectroscopy of Gold Nanotriangle Fabry-Pérot Resonances

Here, in this work, we characterize the near-field response of individual gold nanotriangles over a broad, visible-to-infrared, spectral region (200–1500 meV) using high-resolution electron energy-loss spectroscopy (EELS) performed inside of a scanning transmission electron microscope (STEM). We begin by experimentally imaging the spatial and spectral extent of each nanotriangle’s plasmonic Fabry-Pérot modes and measuring the evolution of their resonance energies with increasing edge length; thereby providing detailed information on infrared plasmon dephasing times and dispersion relations. Numerical electromagnetic simulations of the electron probe are used to interpret these experimental results and to compare the near-field electromagnetic enhancement factors of gold nanotriangles and nanorods of equal resonant energy. Taken together, this combined experimental and theoretical study provides unique insights relevant to designing noble metal plasmonic nanoparticle systems for solar energy harvesting and sensing applications in the near- and mid-infrared.

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

Magnetorotational instability in a swirling partially ionized gas

ABSTRACT The magnetorotational instability (MRI) has been proposed as the method of angular momentum transport that enables accretion in astrophysical discs. However, for weakly ionized discs, such as protoplanetary discs, it remains unclear whether the combined non-ideal magnetohydrodynamic (MHD) effects of Ohmic resistivity, ambipolar diffusion, and the Hall effect make these discs MRI stable. While much effort has been made to simulate non-ideal MHD MRI, these simulations make simplifying assumptions and are not always in agreement with each other. Furthermore, it is difficult to directly observe the MRI astrophysically because it occurs on small scales. Here, we propose the concept of a swirling gas experiment of weakly ionized argon gas between two concentric cylinders threaded with an axial magnetic field that can be used to study non-ideal MHD MRI. For our proposed experiment, we derive the hydrodynamic equilibrium flow and a dispersion relation for MRI that includes the three non-ideal effects. We solve this dispersion relation numerically for the parameters of our proposed experiment. We find it should be possible to produce a non-ideal MRI in such an experiment because of the Hall effect, which increases the MRI growth rate when the vertical magnetic field is anti-aligned with the rotation axis. As a proof of concept, we also present experimental results for a hydrodynamic flow in an unmagnetized prototype. We find that our prototype has a small, but non-negligible, α-parameter that could serve as a baseline for comparison to our proposed magnetized experiment, which could be subject to additional turbulence from the MRI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Experimental observation of a Raman-induced temporal waveguide

The formation of a Raman-induced temporal waveguide is demonstrated by launching short pump and probe pulses inside a photonic crystal fiber. The pump pulse creates a fundamental soliton whose speed changes continuously owing to its deceleration through the Raman-induced red shift of its spectrum. The spectrum of the probe pulse is blue-shifted to ensure that the two pulses move at the same speed and follow the same trajectory over the entire length of the fiber. Furthermore, the output wavelengths of the pump and probe pulses depend on the peak power of input pump pulses and their measured values agree with the predictions based on the dispersion data. Numerical modeling also shows good agreement with the experimental results.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Colloidal Stability of PFSA-Ionomer Dispersions Part II: Determination of Suspension pH Using Single-Ion Potential Energies

Perfluorosulfonic acid (PFSA) ionomers serve a vital role in the performance and stability of fuel-cell catalyst layers. These properties, in turn, depend on the colloidal processing of precursor inks. Here, to understand the colloidal structure of fuel-cell catalyst layers, we explore the aggregation of PFSA ionomers dissolved in water/alcohol solutions and relate the predicted aggregation to experimental measurements of solution pH. Not all side chains contribute to measured pH because of burying inside particle aggregates. To account for the measured degree of dissociation, a new description is developed for how PFSA aggregates interact with each other. The developed single-counterion electrostatic repulsive pair potential from Part I is incorporated into the Smoluchowski collision-based kinetics of interacting aggregates with buried side chains. We demonstrate that the surrounding solvent mixture affects the degree of aggregation as well as the pH of the system primarily through the solution dielectric permittivity, which drives the strength of the interparticle repulsive energies. Successful pH prediction of Nafion ionomer dispersions in water/n-propanol solutions validates the numerical calculations. Nafion-dispersion pH measurements serve as a surrogate for Nafion particle-size distributions. The model and framework can be leveraged to explore different ink formulations.

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

Implicit highly-coupled single-ion Hall-MHD formulation for hybrid particle-in-cell codes

The rudiments of a particle-based single-fluid two-temperature magnetohydrodynamic (MHD) algorithm have been outlined in Thoma et al. (2013). The extension of this algorithm to include the effect of Hall physics is described in this paper. An implicit leapfrog version of the algorithm, which allows timesteps large compared to the resistive decay time and other relevant timescales, has recently been added to a hybrid particle-in-cell code. In standard MHD the Hall term in the generalized Ohm’s law can often be neglected when the Hall parameter is small. This term must, however, be retained in regimes where it is non-negligible. The retention of displacement current in Maxwell’s equations avoids the numerical difficulties associated with the whistler mode, which are encountered in standard explicit Hall-MHD codes, and allows the algorithm to be incorporated into hybrid particle-in-cell codes, for which particles may migrate from a kinetic to fluid to MHD description based upon local ambient plasma conditions. A highly-coupled implicit Hall-MHD formalism is presented, in which displacement current can either be retained or neglected. Even when displacement current is neglected, the highly-coupled implicit formalism avoids the restrictive timesteps for the whistler mode in explicit Hall-MHD codes. A comparison of numerical and analytic dispersion analysis demonstrates the feasibility of this approach and establishes relevant constraints to assure numerical stability. The implementation of the algorithm is described, and test simulation results in 1D and 2D in both linear and nonlinear regimes are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Numerical investigation of gas-filled multipass cells in the enhanced dispersion regime for clean spectral broadening and pulse compression

We show via numerical simulations that the regime of enhanced frequency chirp can be achieved in gas-filled multipass cells. Our results demonstrate that there exists a region of pulse and cell parameters for which a broad and flat spectrum with a smooth parabolic-like phase can be generated. This spectrum is compatible with clean ultrashort pulses, whose secondary structures are always below the 0.5% of its peak intensity such that the energy ratio (the energy contained within the main peak of the pulse) is above 98%. This regime makes multipass cell post-compression one of the most versatile schemes to sculpt a clean intense ultrashort optical pulse.

Segundo Staels, V. W. (ORCID:0000000200056157)↗
DOE OSTI2021

Modeling Explosive Dispersal of Radiological Material: A Nonparametric Analysis

The National Atmospheric Release Advisory Center (NARAC) is a Department of Energy (DOE) asset that provides guidance on the potential effects of atmospherically released hazards based on dispersion simulation results. Numerical simulation results are compared to well-characterized experimental measurements whenever possible to validate and refine the model predictions. In 2012, such measurements were obtained by Defence Research and Development Canada (DRDC) and collaborating organizations for the case of explosively released 140 La – a radioactive isotope with a half-life of 1.6781 days.

42 ENGINEERING↗
DOE OSTI2026

Artificial Intelligence in Nuclear Safeguards; Evaluating Safeguards and Security Risks and Benefits for Advanced and Small Modular Reactor Deployments

Rapidly growing interest in advanced and small modular reactor (A/SMR) technologies presents challenges as well as opportunities for implementing international safeguards and security. A/SMR deployments are expected to be more numerous, more geographically dispersed, and more varied in their designs, placing new demands on the data systems and analytical tools used to support oversight (Alberti et al., 2023; Canadian Nuclear Safety Commission et al., 2024). Because of this variability, the importance and reliance on data systems for A/SMR deployments is expected to be higher than for previous reactor generations. Artificial Intelligence and Machine Learning (AI/ML) offer potential capabilities to address the high variability inherent in A/SMR technology. The beneficiaries of AI-assisted tools include facility operators, government regulators, IAEA inspectors, and A/SMR vendors. This report analyzes how AI/ML-assisted technologies can strengthen the implementation of IAEA safeguards and security measures. It also identifies AI-assisted tools to strengthen operator, facility, and regulator knowledge management practices and examines the potential risks AI/ML-based tools may introduce to IAEA safeguards and security efforts. It concludes with a set of hypothetical, standards-style requirements for AI/ML systems used in safeguards contexts, grounded in an inspector-centric view of system verification. Despite the potential benefits of AI/ML systems, understanding potential intentional and unintentional failure modes is critical for ensuring adequate protection of nuclear materials and facilities. Unique features of A/SMRs including sealed cores, remote and novel paradigms of operation, off-site reactor fabrication, novel fuel forms, and varied refueling requirements, introduce challenges for traditional safeguards technological approaches (Pensado et al., 2024; Federation of American Scientists, 2025). AI/ML systems deployed to address these challenges may introduce new risks requiring systematic evaluation rooted in both AI-specific risk frameworks, such as the NIST AI Risk Management Framework (NIST AI RMF), and established cyber risk management standards such as NIST SP 800-30 (National Institute of Standards and Technology [NIST], 2023; NIST, 2012).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Exascale Multiphysics Nuclear Reactor Simulations for Advanced Designs

ENRICO is a coupled application developed under the U.S. Department of Energy's Exascale Computing Project (ECP) targeting the modeling of advanced nuclear reactors. It couples radiation transport with heat and fluid simulation, including the high-fidelity, highresolution Monte-Carlo code Shift and the Computational fluid dynamics code NekRS. NekRS is a highly-performant open-source code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on rapidly convergent high-order spectral element discretizations that feature minimal numerical dissipation and dispersion. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Frontier, Nek5000/RS has recently achieved an unprecedented milestone in breaching over 1 billion spectral elements and 350 billion degrees of freedom. Shift has demonstrated the capability to transport upwards of 1 billion particles per second in full core nuclear reactor simulations featuring complete temperature-dependent, continuous-energy physics on Frontier. Shift achieved a weak-scaling efficiency of 97.8% on 8192 nodes of Frontier and calculated 6 reactions in 214,896 fuel pin regions below 1% statistical error yielding first-of-a-kind resolution for a Monte Carlo transport application.

Hamilton, Steven P.↗
DOE OSTI2023

Explicit dispersion relations for warm fluid waves in a uniform plasma (invited)

Classical dispersion relations for waves in a fluid plasma are expressed as implicit functions of wave frequency and wave- number. The implicit dispersion relation must in general be solved numerically. This work introduces an astute method for obtaining an explicit dispersion relation for waves in a fluid plasma. The explicit expression has an advantage of providing eigenmodes without numerical computations and giving a dispersion relation more detailed than the implicit expression. As in the case of cold waves, the wavenumber for an arbitrary frequency can be directly obtained from the explicit formula. The explicit formula also enables an accurate evaluation of finite-temperature effects on a dispersion relation even near the cyclotron resonance, where a numerical analysis of the implicit relation is impractical. As a result, the analytic formula can be used to investigate temperature effects on electromagnetic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

BASS. XXVI. DR2 Host Galaxy Stellar Velocity Dispersions

We present new central stellar velocity dispersions for 484 Sy 1.9 and Sy 2 from the second data release of the Swift/BAT AGN Spectroscopic Survey (BASS DR2). This constitutes the largest study of velocity dispersion measurements in X-ray-selected obscured active galactic nuclei (AGN) with 956 independent measurements of the Ca II H and K λ3969, 3934 and Mg I λ5175 region (3880–5550 Å) and the calcium triplet region (8350–8730 Å) from 642 spectra mainly from VLT/X-Shooter or Palomar/DoubleSpec. Our sample spans velocity dispersions of 40–360 km s 1 , corresponding to 4–5 orders of magnitude in black hole mass (M BH = 10 5.5-9.6 M ⊙ ), bolometric luminosity (L bol ~ 10 42–46 erg s -1 ), and Eddington ratio (L/L Edd ~ 10 -5 to 2). For 281 AGN, our data and analysis provide the first published central velocity dispersions, including six AGN with low-mass black holes (M BH = 10 5.5-6.5 M ⊙ ), discovered thanks to high spectral resolution observations (σ inst ~ 25 km s -1 ). The survey represents a significant advance with a nearly complete census of velocity dispersions of hard X-ray–selected obscured AGN with measurements for 99% of nearby AGN (z < 0.1) outside the Galactic plane (|b| > 10°). The BASS AGN have much higher velocity dispersions than the more numerous optically selected narrow-line AGN (i.e., ~150 versus ~100 km s -1 ) but are not biased toward the highest velocity dispersions of massive ellipticals (i.e., >250 km s -1 ). Despite sufficient spectral resolution to resolve the velocity dispersions associated with the bulges of small black holes (~10 4–5 M ⊙ ), we do not find a significant population of super-Eddington AGN. Using estimates of the black hole sphere of influence from velocity dispersion, direct stellar and gas black hole mass measurements could be obtained with existing facilities for more than ~100 BASS AGN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Nonlinear noise suppression in a high-rate optical channel with phase modulation and dispersion compensation

We have performed a numerical simulation of transmission of optical DP QPSK signals at a rate of 40 Gbaud over a distance of 2000 km. It is shown that a channel with large negative dispersion and preliminary large positive pulse chirping significantly reduces nonlinear distortion and improves the signal quality. To further reduce the probability of recognition error, it is proposed, instead of a negative dispersion fibre, to use compensators with low attenuation. (fibre-optic communication lines)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Machined phononic crystals to block high-order Lamb waves and crosstalk in through-metal ultrasonic communication systems

For systems that require complete metallic enclosures (e.g., containment buildings for nuclear reactors), it is impossible to access interior sensors and equipment using standard electromagnetic techniques. A viable way to communicate and supply power through metallic barriers is the use of elastic waves and ultrasonic transducers, introducing several design challenges that must be addressed. Specifically, the use of multiple communication channels on the same enclosure introduces an additional mechanism for signal crosstalk between channels: guided waves propagating in the barrier between channels. This work numerically and experimentally investigates a machined phononic crystal to block MHz Lamb wave propagation between ultrasonic communication channels, greatly reducing wave propagation and the resulting crosstalk voltage. Blind grooves are machined into one or both sides of a metallic barrier to introduce a periodic unit cell, greatly altering the guided wave dispersion in the barrier. Numerical simulations are used to determine a set of groove geometries for testing, and experiments were performed to characterize the wave-blocking performance of each design. The best-performing design was tested using piezoelectric transducers bonded to the barrier, showing a 14.4 dB reduction in crosstalk voltage. Overall, the proposed periodic grooving method is a promising technique for completely isolating ultrasonic power/data transfer systems operating in a narrow frequency range.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗