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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Direct simulation of particle dispersion in a decaying isotropic turbulence

Results of a numerical investigation of the dispersion of solid particles in decaying isotropic turbulence are presented. The 3D time-dependent velocity field of a homogeneous nonstationary turbulence is computed using the method of direct numerical simulation (DNS). The dispersion characteristics of three different solid particles (corn, copper, and glass) injected in the flow are obtained by integrating the complete equation of particle motion along the instantaneous trajectories of 22-cubed particles for each particle type, and then by performing ensemble averaging. Good agreement was achieved between the present DNS results and the measured time development of the mean-square displacement of the particles. Questions of how and why the dispersion statistics of a solid particle differ from those of its corresponding fluid point and surrounding fluid and what influences inertia and gravity have on these statistics are also discussed.

Elghobashi, S.↗
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↗
NASA NTRS2011

Numerical Study of Ammonia Leak and Dispersion in the International Space Station

Management of off-nominal situations on-board the International Space Station (ISS) is important to its continuous operation. One situation of concern is an accidental release of a chemical into the ISS atmosphere. In particular, introduction of ammonia into the cabin atmosphere can occur via the interface heat exchangers (IFHX) between the external thermal control system containing ammonia and internal thermal control system that uses water as a coolant to remove heat from ISS subsystems. Breach of the water/ammonia barrier of the IFHX can lead to a catastrophic rupture. Once the liquid water/ammonia mixture exits the ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin atmosphere that results in rapid contamination of the cabin. The goal of the study is to assess the amount of ammonia in the Russian Segment by the time the crew is able to isolate the U.S. Segment. A Computational Fluid Dynamics (CFD) model for an accurate prediction of airflow and ammonia transport in the frozen flow field within the assembly complete ISS cabin was developed. The localized effects of ammonia dispersion are examined and discussed.

Son, Chang H.↗
NASA NTRS2012

Numerical Study of Ammonia Leak and Dispersion in the International Space Station

Release of ammonia into the International Space Station (ISS) cabin atmosphere can occur if the water/ammonia barrier breach of the active thermal control system (ATCS) interface heat exchanger (IFHX) happens. After IFHX breach liquid ammonia is introduced into the water-filled internal thermal control system (ITCS) and then to the cabin environment through a ruptured gas trap. Once the liquid water/ammonia mixture exits ITCS, it instantly vaporizes and mixes with the U.S. Laboratory cabin air that results in rapid deterioration of the cabin conditions. The goal of the study is to assess ammonia propagation in the Station after IFHX breach to plan the operation procedure. A Computational Fluid Dynamics (CFD) model for accurate prediction of airflow and ammonia transport within each of the modules in the ISS cabin was developed. CFD data on ammonia content in the cabin aisle way of the ISS and, in particular, in the Russian On- Orbit Segment during the period of 15 minutes after gas trap rupture are presented for four scenarios of rupture response. Localized effects of ammonia dispersion and risk mitigation are discussed.

Son, Chang H.↗
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↗
NASA NTRS1978

The effect of the London-van der Waals dispersion force on interline heat transfer

A theoretical procedure to determine the heat transfer characteristics of the interline region (junction of liquid-solid-vapor) from the macroscopic optical and thermophysical properties of the system is outlined. The analysis is based on the premise that the interline transport processes are controlled by the London-van der Waals dispersion force between condensed phases (solid and liquid). Numerical values of the dispersion constant are presented. The procedure is used to compare the relative size of the interline heat sink of various systems using a constant heat flux mode. This solution demonstrates the importance of the interline heat flow number, which is evaluated for various systems.

Wayner, P. C., Jr.↗
NASA NTRS1996

Excitation of Plasma Waves in Aurora by Electron Beams

In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

daSilva, C. E.↗
NASA NTRS2004

Computation of Gust-Cascade Interaction Using the CE/SE Method

The problem 2 in Category 3 of the 4th Computational Aeroacoustic(CAA) Workshop is solved using the space-time conservation element and solution element (CE/SE) method. This problem models rotor-stator interaction in a 2D cascade. It involves complex geometries and flow physics including vortex shedding and acoustic radiation. The parallel version of the 2D nonlinear Euler solver is used with an unstructured triangular mesh to solve this problem. The Giles approach is incorporated with the CE/SE method to handle non-equal pitches of the rotor and stator. Validation on the Giles approach is performed using Problem 3.1 in the 2nd CAA Workshop. The space-time CE/SE method is a finite volume method with second-order accuracy in both space and time. The flux conservation is enforced in both space and time instead of space only. It has low numerical dissipation and dispersion errors. It uses simple non-reflecting boundary conditions and is compatible with unstructured meshes. It is simple, flexible, and generate reasonably accurate solutions. The CE/SE method has been successfully applied to solve numerous practical problems, especially aeroacoustic problems. Some preliminary numerical results of the benchmark problem 3.2 of the 4th CAA Workshop are shown. The steady-state pressure contour is plotted. The mean pressure distribution on the blade surface is compared with Turbo solution showing a good agreement. The sound pressure level versus the rotor harmonic n at the six designated positions on the blade surface, three locations at inlet plane, and three locations at the outlet plane are plotted. It can be seen that the acoustic response exists only at the excitation frequencies (n = 1,2,3). On the blade surface, the acoustic wave at n = 1 is dominant, while at the inlet and outlet planes, the sound pressure level at n = 2 becomes the largest, which is similar to the results presented. The distribution of sound pressure level at different spatial modes along the z- direction is plotted for n = 1,2,3, respectively. It shows that the spatial modes m = -32 and 22 at n = 1 exponentially decay, and the spatial modes m = 10 at n = 2, m = -42 and 12 at n = 3 propagate both upstream and downstream, which agrees with the prediction based on the linearized theory. Some oscillations are observed, which needs to be investigated further. In the final paper, the numerical results will be compared with a frequency-domain solver LINFLUX solution if it is available.

Wang, X.-Y.↗
NASA NTRS1974

Scattering and absorption from poly-dispersed aerosols

Results from numerical computations of the diffuse sky radiances and degree of polarization are given. These results were obtained by solving the equation of radiative transfer in the formulation of Eschelbach. In these computations, multiple scattering and absorption by aerosol particles were considered. In addition, results from experimental measurements of the complex refractive index of aerosol samples at different relative humidities are presented. From the radiation flux divergences which were computed based on the determined properties of atmospheric aerosol particles, atmospheric heating rates were derived which were found to be comparable to the heating rates by water vapor.

Bullrich, K.↗
NASA NTRS2004

Carbon Dioxide Dispersion in the Combustion Integrated Rack Simulated Numerically

When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Wu, Ming-Shin↗
NASA NTRS1997

Numerical Boundary Conditions for Computational Aeroacoustics Benchmark Problems

Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.

Tam, Chritsopher K. W.↗
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)↗
NASA NTRS1988

Excitations of low-frequency hydromagnetic waves by freshly created ions in the solar wind

Low-frequency hydromagnetic waves excited by newborn ions in the solar wind plasma are studied. The freshly created ions appear in the solar wind frame with a ring beam distribution. Both Alfven and fast magnetosonic waves are made unstable by the presence of the newborn ions. The dependence of the growth rate of both waves on the newborn ion density, the angle between the interplanetary magnetic field (IMF) and solar wind flow, and the angle of wave propagation relative to the IMF is investigated. Analytic approximations for the growth rates are presented, and numerical solutions of the dispersion equation are shown. The approximations are quite close to the numerically determined growth rates. It is found that the waves grow preferentially in the direction parallel to the IMF, and that the growth rates increase with both newborn ion density and the angle between the IMF and the solar wind flow.

Price, C. P.↗
NASA NTRS1995

A compact solution to computational acoustics

This paper demonstrates that the linearized, dimensional Euler equations for acoustic computation can be accurately solved as a set of decoupled first-order wave equations, and that if ordered properly, this system of simple waves has unambiguous, easily implemented boundary conditions, allowing waves of same group speeds to pass through numerical boundaries or comply with wall conditions. Thus, the task of designing a complex multi-dimensional scheme with approximate far-field boundary conditions reduces to the design of higher order schemes for the one-dimensional simple wave equation. A compact finite-difference scheme and a characteristically exact but numerically n(th) order accurate boundary condition are introduced for solving the first order wave equation. Spanning a three-point two-level stencil, this low-dispersion implicit scheme has a third order spatial accuracy when used on nonuniform meshes, fourth order accurate on uniform meshes, and a temporal accuracy of second order due to the choice of trapezoidal integration for algorithmic simplicity. The robustness and accuracy of the scheme are demonstrated through a series of numerical experiments and comparisons with published results. When tested on the one-dimensional wave equation on a uniform grid, this scheme allows a Gaussian wave packet to pass through any finite domain with low numerical dispersion characteristic of a spatially fourth-order scheme and reflections at numerical boundaries maintained below truncation error. On highly stretched and irregular grids, only mild dispersions are found in the solution while solutions by other methods fail or are severely distorted. Yet, this scheme is no more sophisticated to solve or implement than the Crank-Nicolson scheme. This scheme has been tested on four categories of the ICASE/LaRC benchmark problems, which include propagation of acoustic and convective waves in Cartesian and cylindrical domains, reflection of acoustic wave at stationary/moving boundaries, and sound generation by gust-blade interaction.

Fung, K.-Y.↗
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↗