Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave parameter”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Lagrangian particle simulation of hydrogen pellets and SPI into runaway electron beam in ITER

Numerical studies of the ablation of pellets and shattered pellet injection (SPI) fragments into a runaway electron beam in ITER have been performed using a time-dependent pellet ablation code [Samulyak et al., Nucl. Fusion, 61(4), 046007 (2021)]. The code resolves detailed ablation physics near pellet fragments and large-scale expansion of ablated clouds. The study of a single-fragment ablation quantifies the influence of various factors, in particular, the impact ionization by runaway electrons and cross-field transport models, on the dynamics of ablated plasma and its penetration into the runaway beam. Simulations of SPI performed using different numbers of pellet fragments study the formation and evolution of the ablation clouds and their large-scale dynamics in ITER. In conclusion, the penetration depth of the ablation clouds is found to be of the order of 50 cm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Hall interchange instability as a seed for helical magneto-Rayleigh–Taylor instabilities in magnetized liner inertial fusion Z-Pinches scaled from Z-Machine parameters to a next generation pulsed power facility

Magnetized liner inertial fusion (MagLIF) is a magneto-inertial-fusion concept that is studied on the 20-MA, 100-ns rise time Z Pulsed Power Facility at Sandia National Laboratories. Given the relative success of the platform, there is a wide interest in studying the scaled performance of this concept at a next-generation pulsed-power facility that may produce peak currents upward of 60 MA. An important aspect that requires more research is the instability dynamics of the imploding MagLIF liner, specifically how instabilities are initially seeded. It has been shown in magnetized 1-MA thin-foil liner Z-pinch implosion simulations that a Hall interchange instability (HII) effect can provide an independent seeding mechanism for helical magneto-Rayleigh–Taylor instabilities. Here in this paper, we explore this instability at higher peak currents for MagLIF using 2D discontinuous Galerkin PERSEUS simulations, an extended magneto-hydrodynamics code, which includes Hall physics. Our simulations of scaled MagLIF loads show that the growth rate of the HII is invariant to the peak current, suggesting that studies at 20-MA are directly relevant to 60-MA class machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Kinetic simulations comparing quasi-parallel and quasi-perpendicular piston-driven collisionless shock dynamics in magnetized laboratory plasmas

Magnetized collisionless shocks are common in astrophysical systems, and scaled versions can be created in laboratory experiments by utilizing laser-driven piston plasmas to create these shocks in a magnetized background plasma. A key parameter for these experiments is the angle θB between the shock propagation direction and the background magnetic field. We performed quasi-1D piston-driven shock simulations to explore shock formation, evolution, and key observables relevant to laboratory experiments for a range of shock angles between θB=90° to θB=30°. Our results show that the spatial and temporal scales of shock formation for all angles considered are similar when expressed in terms of the perpendicular component of the magnetic field. In a steady state, ion and electron temperatures become more isotropic, and the electron-to-ion temperature ratio is higher for smaller θB. At θB=30°, ion heating parallel to the magnetic field becomes dominant, associated with more ions being reflected at one discontinuity and subsequently trapped by the next discontinuity due to shock reformation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗
DOE OSTI2025

Spatially and temporally resolved plasma parameter estimations of laser heated MagLIF relevant gas pipes at NIF

The ability to control laser pre-heat is an integral part of the inertial confinement fusion concept known as Magnetized Liner Inertial Fusion. This process is studied at the National Ignition Facility (NIF) where 4 of the 192 laser beams are propagated through a 1-cm long gas cell where they deposit >20 kJ of energy into the gaseous fuel via inverse bremsstrahlung absorption. This process ionizes the gas, producing a plasma that follows behind the laser front and expands over the radius of the cell. Emission from this plasma, as viewed by a gated x-ray detector, can be used to build spatially and temporally resolved estimations of the pre-heat plasma's density and temperature profiles. This can then be used to estimate the plasma pressure, internal energy, and radiation losses. Estimations show the evolution of the plasma in magnetized and unmagnetized gas cells filled with ambient temperature neopentane (C5H12) +1% Ar, as well as unmagnetized cryogenically cooled (32 K) deuterium +1% Ne filled targets. This analysis shows the effects of initial gas-fill density, composition, and axial magnetization on the time-dependent plasma parameters. Previously, these parameters at the NIF had not been experimentally characterized, and these estimations provided a potential new means of testing radiation magneto-hydrodynamic predictive capability models. Results in unmagnetized targets have strong agreement with simulations. However, in targets with a 19 T applied axial magnetic field, this method yields electron temperatures up to 100% hotter than those predicted by HYDRA codes.

Bremsstrahlung↗
DOE OSTI2025

A theoretical model for sheath dynamics in the two-plasma mode of a cylindrical filament discharge

We present a refined understanding of cylindrical plasma-facing emitter cathodes by examining the limitations of the conventional space charge limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge limited regimes based on electron current constraints, our findings reveal that trapped ions can form a quasi-neutral “second plasma,” expanding upstream and interacting with the primary plasma. In cylindrical and other non-planar geometries, an expanding trapped-ions layer produces current enhancement by converting more of the emitted electron flux to passing flux. The current enhancement leads to a complex coupling between the trapped-ions plasma and upstream plasma through ionizing collisions in the upstream region and charge exchange collisions in the virtual cathode. The coupling results in mutual density enhancement and at the same time, antagonistic expansion dynamics between the two plasma layers. The process is understood using a new “aid-compete” model. Numerical solutions, validated against simulations, reveal the model's capacity to predict system evolution and complex discharge behaviors, albeit with qualitative approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Platform for 100 s Mbar equation of state measurements on the National Ignition Facility

Equation of state (EOS) measurements in the 100 s Mbar range are needed to underwrite models employed in the simulation of high energy density plasmas. To this end, a platform has been developed for fielding on the National Ignition Facility, capable of producing high-quality impedance match EOS data, wherein a planar, high-pressure, steady shock is driven into a sample package, and sample and reference standard shock velocities are measured. This platform, dubbed planar high pressure, or PHP, was fielded with an initial proof-of-concept shot in January 2023. The first PHP shot, aiming to study gold, demonstrated a pressure close to 400 Mbar, two orders of magnitude higher than previously reported gold EOS data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

QUBIC II: Spectral polarimetry with bolometric interferometry

Bolometric interferometry is a novel technique that has the ability to perform spectral imaging. A bolometric interferometer observes the sky in a wide frequency band and can reconstruct sky maps in several sub-bands within the physical band in post-processing of the data. This provides a powerful spectral method to discriminate between the cosmic microwave background (CMB) and astrophysical foregrounds. In this paper, the methodology is illustrated with examples based on the Q & U Bolometric Interferometer for Cosmology (QUBIC) which is a ground-based instrument designed to measure the B-mode polarization of the sky at millimeter wavelengths. We consider the specific cases of point source reconstruction and Galactic dust mapping and we characterize the point spread function as a function of frequency. We study the noise properties of spectral imaging, especially the correlations between sub-bands, using end-to-end simulations together with a fast noise simulator. Finally, we conclude showing that spectral imaging performance are nearly optimal up to five sub-bands in the case of QUBIC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Pairing tendencies in the doped Kitaev-Heisenberg model

Here, we study the impact of hole doping on the Kitaev-Heisenberg model on the honeycomb lattice. We investigate the pairing tendencies and correlation functions in the framework of a t - J - K model using density matrix renormalization group calculations on three-leg cylinders. In the case of the pure Kitaev model, which realizes a quantum spin-liquid phase at half-filling, we find that binding of two holes only occurs at low values of the hopping, where the holes are slow. We have theoretically verified that pair formation occurs in the limit of immobile holes, where the pure Kitaev model remains exactly solvable. When we instead fix the hopping at an intermediate, more realistic, value, and vary the Heisenberg and Kitaev interaction strengths, we find pairing tendencies only in the Néel phase. This is in contrast to prior mean-field calculations, highlighting the importance of accounting for the kinetic energy of dopants in generalized Kitaev models. Interestingly, we also find signatures of pair-density wave formation over the studied range of model parameters, namely, a periodic modulation of the charge density as well as the spin-spin and pair-pair correlations in real space. Moreover, we present a comparative study of the different correlations as a function of doping. We finally discuss the potential for experimentally observing the studied physics in quantum materials and heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Low-energy Scattering and Effective Interactions of Two Baryons at $m_{\pi}\sim$ 450 MeV from Lattice Quantum Chromodynamics

The interactions between two octet baryons are studied at low energies using lattice QCD (LQCD) with larger-than-physical quark masses corresponding to a pion mass of $m_{\pi}\sim 450$ MeV and a kaon mass of $m_{K}\sim 596$ MeV. The two-baryon systems that are analyzed range from strangeness $S=0$ to $S=-4$ and include the spin-singlet and triplet $NN$, $\Sigma N$ ($I=3/2$), and $\Xi\Xi$ states, the spin-singlet $\Sigma\Sigma$ ($I=2$) and $\Xi\Sigma$ ($I=3/2$) states, and the spin-triplet $\Xi N$ ($I=0$) state. The $s$-wave scattering phase shifts, low-energy scattering parameters, and binding energies when applicable, are extracted using L\"uscher's formalism. While the results are consistent with most of the systems being bound at this pion mass, the interactions in the spin-triplet $\Sigma N$ and $\Xi\Xi$ channels are found to be repulsive and do not support bound states. Using results from previous studies at a larger pion mass, an extrapolation of the binding energies to the physical point is performed and is compared with experimental values and phenomenological predictions. The low-energy coefficients in pionless EFT relevant for two-baryon interactions, including those responsible for $SU(3)$ flavor-symmetry breaking, are constrained. The $SU(3)$ symmetry is observed to hold approximately at the chosen values of the quark masses, as well as the $SU(6)$ spin-flavor symmetry, predicted at large $N_c$. A remnant of an accidental $SU(16)$ symmetry found previously at a larger pion mass is further observed. The $SU(6)$-symmetric EFT constrained by these LQCD calculations is used to make predictions for two-baryon systems for which the low-energy scattering parameters could not be determined with LQCD directly in this study, and to constrain the coefficients of all leading $SU(3)$ flavor-symmetric interactions, demonstrating the predictive power of two-baryon EFTs matched to LQCD.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Fine structure of pentaquark multiplets in the dynamical diquark model

In this work, we apply the dynamical diquark model to predict the spectrum of hidden-charm pentaquark states in both unflavored and open-strange sectors. Using only Hamiltonian parameters introduced in the tetraquark S - and P -wave multiplets, the model naturally produces the level spacing supported by the most recent LHCb results for P c structures. Furthermore, using model inputs obtained from data of hidden-charm, open-strange tetraquarks ( Z c s ), we predict the spectrum of P c s states, including the recently observed P c s ( 4459 ) . We find all pentaquark candidates observed to date to belong to the 1 P multiplet and hence have positive parity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Nonlinear ultrasonic technique for the characterization of microstructure in additive materials

This study employs nonlinear ultrasonic techniques to track microstructural changes in additively manufactured metals. The second harmonic generation technique based on the transmission of Rayleigh surface waves is used to measure the acoustic nonlinearity parameter, β. Stainless steel specimens are made through three procedures: traditional wrought manufacturing, laser-powder bed fusion, and laser engineered net shaping. The β parameter is measured through successive steps of an annealing heat treatment intended to decrease dislocation density. Dislocation density is known to be sensitive to manufacturing variables. In agreement with fundamental material models for the dislocation-acoustic nonlinearity relationship in the second harmonic generation, β drops in each specimen throughout the heat treatment before recrystallization. Geometrically necessary dislocations (GNDs) are measured from electron back-scatter diffraction as a quantitative indicator of dislocations; average GND density and β are found to have a statistical correlation coefficient of 0.852 showing the sensitivity of β to dislocations in additively manufactured metals. Moreover, β shows an excellent correlation with hardness, which is a measure of the macroscopic effect of dislocations.

42 ENGINEERING↗
DOE OSTI2021

A Shallow-Water Model for Convective Self-Aggregation

Randomly distributed convective storms can self-aggregate in the absence of large-scale forcings. Here we present a 1D shallow-water model to study the convective self-aggregation. This model simulates the dynamics of the planetary boundary layer and represents convection as a triggered process. Once triggered, convection lasts for finite time and occupies finite length. Here we show that the model can successfully simulate self-aggregation, and that the results are robust to a wide range of parameter values. In the simulations, convection excites gravity waves. The gravity waves then form a standing wave pattern, separating the domain into convectively active and inactive regions. We analyze the available potential energy (APE) budget and show that convection generates APE, providing energy for self-aggregation. By performing dimensional analysis, we develop a scaling theory for the size of convective aggregation, which is set by the gravity wave speed, damping time scale, and number density of convective storms. This paper provides a simple modeling framework to further study convective self-aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1972

Viscous stability theory for thermally stratified flow - Discontinuous jets and shear layers.

The method used by Drazin to investigate the stability of unbounded, viscous, homogeneous, parallel shear flow to small wavenumber disturbances is extended to study the effect of thermal stratification on the stability of unbounded jets and shear layers. By this method the stability characteristics of continuous profiles are inferred from the stability characteristics of discontinuous profiles. The characteristic value problem for discontinuous jet and shear layers is posed by the requirement that the solutions of the governing differential equation satisfy the matching conditions and boundedness conditions for layers that extend to infinity. The analysis leads to a characteristic determinant which is required to vanish for the characteristic values of the parameters: the Reynolds number, the wavenumber, and the wave speed. The stabilizing effect of the thermal stratification as parameterized by the Richardson number was found to be most stabilizing for small wavenumber (large-scale) disturbances.

Gage, K. S.↗