SEARCH · Engineering Papers

Results for “Damping”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Anomalous resistivity and electron heating by lower hybrid drift waves inside reconnecting current sheets

Inside an electron diffusion region of laboratory reconnection experiments, the quasi-electrostatic lower hybrid drift wave (ES-LHDW) is observed when a significant guide field component is present. Through direct measurement of the anomalous drag term and quasilinear analysis, it is shown that ES-LHDW can account for approximately 20% of the mean reconnection electric field in a case with moderate guide field. This value exceeds the contribution from classical resistivity, which is around 10%. The effects of the Lorentz force term, often neglected for electrostatic waves, are crucial for the observed correlation between electric field and density fluctuations. Anomalous electron heating by the perturbed current and resistivity (2.6 MW/m 3 ) also surpasses the classical Ohmic heating, which is about 2.0 MW/m 3 . For the case with a high guide field, significantly higher local electron temperatures were observed during periods of strong ES-LHDW activity. A statistical analysis further supports electron heating by LHDW, showing a larger increase in electron temperature with a high guide field. Finally, data from the Magnetospheric Multiscale mission provide evidence of Landau damping of ES-LHDW, suggesting that ES-LHDW may contribute to the generation of nonthermal electrons along the direction parallel to the magnetic field.

Magnetic reconnection↗

DOE OSTI2025

Conservative projection-based data-driven model order reduction of a fluid-kinetic spectral solver

Kinetic simulations are computationally intensive due to six-dimensional phase space discretization. Many kinetic spectral solvers use the asymmetrically weighted Hermite expansion due to its conservation and fluid-kinetic coupling properties, i.e., the lower-order Hermite moments capture and describe the macroscopic fluid dynamics, and higher-order Hermite moments describe the microscopic kinetic dynamics. We leverage this structure by developing a parametric data-driven reduced-order model based on the proper orthogonal decomposition, which projects the higher-order kinetic moments while retaining the fluid moments intact. We demonstrate analytically and numerically that the method ensures local and global mass, momentum, and energy conservation. The numerical results show that the proposed method effectively replicates the high-dimensional spectral simulations at a fraction of the computational cost and memory, as validated on the weak Landau damping and two-stream instability benchmark problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Diffusion-mediated passing of molecular species in linear nanopores constrained by orientational alignment

For diffusion-mediated catalytic conversion reactions in materials with narrow linear nanopores, e.g., mesoporous silica MCM-41, a key parameter is the propensity for product species to be able to pass reactant species and thus to efficiently exit the pore. For elongated species, this can require orientational alignment with the pore axis. In this work, we perform benchmark analyses for such solution-phase systems where one of these species is elongated in order to quantify the dependence of this passing propensity, P, on pore diameter and on the rotational diffusion coefficient, Dr, of the elongated species. In particular, we consider the passing of a spherical and an elongated spherocylindrical shaped species in a cylindrical pore in an implicit solvent, where these species cannot overlap. Passing is mediated by diffusive Brownian motion of these species as described by strongly damped Langevin dynamics. We quantify scaling of P for pore width just above the threshold where passing is sterically blocked, and also reveal a significant decrease in P for lower Dr. We also consider the dependence of P on the aspect ratio of the elongated species and obtain an exact result in the limiting regime of large aspect ratio.

Brownian motion↗

DOE OSTI2025

Saturation and decay of an induced scattering instability driven by static field asymmetries in a non-neutral plasma

Experimental data on the saturation and decay of a previously described induced scattering instability are presented. The instability is driven by static background field asymmetries and initially produces a single growing mode (ω 1 , k 1 ). As the instability develops, several harmonically related (ω n =nω 1 , k n =nk 1 ) eigenmodes of the plasma column grow up, resulting in a spatially localized moving potential perturbation with δϕ/T ≈ 20%. After a period of rapid growth, the mode amplitudes saturate and vary slowly for roughly 10 5 ω 1 −1 . This quasi-equilibrium ends with an unusual change in the eigenmode phases, wherein δϕ changes from positive to negative. Plasma heating during this period increases the damping of all modes and rapidly quenches the instability. This instability is of particular interest because the coupling with the field asymmetry produces an enhanced level of radial transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Peeling-ballooning modes in spherical tokamaks: Multi-branch instabilities and effects beyond ideal MHD

A number of important physics effects on the stability of relatively high-n (n is the toroidal mode number) peeling-ballooning modes (PBMs) are investigated based on an equilibrium reconstructed from a NSTX discharge, utilizing extended magnetohydrodynamic (MHD) eigenvalue solvers. For a given toroidal mode number n, multiple branches of instabilities are computed, with the total number of unstable branches roughly linearly scaling with n. Most of the unstable branches are located in the plasma core region, but edge-localized branches, i.e., PBMs, are also identified at higher n-numbers. For the single-fluid-wise most unstable PBM with n = 19⁠, stabilizing/destabilizing effects due to various physics beyond ideal MHD are systematically investigated. Plasma toroidal flow is found to be weakly stabilizing. Local flow shear is generally stabilizing as well, with the degree of stabilization depending on the initial growth rate (without flow shear) of the mode. The plasma resistivity can strongly destabilize the PBM within the single-fluid framework. Anisotropic thermal transport, strong parallel sound wave damping, as well as two-fluid effects are all stabilizing to the mode. In particular, diamagnetic stabilization (within the two-fluid model) is found to be very strong for this mode.

Linear stability analysis↗

DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Magnetic anisotropy and spin scattering in (La 2/3 Sr 1/3 )MnO 3 /CaRuO 3 bilayers

Complex ferromagnetic oxides have been identified as possible candidate materials for sources of spin currents. Here we study bilayers of ferromagnetic (La 2/3 Sr 1/3 )MnO 3 (LSMO) and metallic CaRuO 3 (CRO) on LSAT substrates as a model system for spin pumping. Ferromagnetic resonance (FMR) measurements of these bilayers show evidence of spin pumping across the interface in the form of an increase in Gilbert damping with the addition of CRO. FMR indicates that the presence of CRO modifies the magnetic anisotropy of the LSMO. By increasing CRO thickness, we find a reduction of the out-of-plane anisotropy and simultaneous rotation of the easy axis within the plane, from the $\langle$110$\rangle$ to $\langle$100$\rangle$ axis. The evolution of magnetic anisotropy determined by FMR disagrees with that measured by bulk SQUID magnetometry and is accompanied by structural distortions in the LSMO layer as measured by x-ray diffraction, thus suggesting a change in magnetic anisotropy attributed to structural distortions imposed on LSMO by CRO. These results suggest that while LSMO and CRO remain promising candidates for efficient pure spin current generation and detection, respectively, epitaxial integration of perovskites will cause additional changes which must be accounted for in spintronics applications.

36 MATERIALS SCIENCE↗

DOE OSTI2020

The influence of packing structure and interparticle forces on ultrasound transmission in granular media

Ultrasound propagation through externally stressed, disordered granular materials was experimentally and numerically investigated. Experiments employed piezoelectric transducers to excite and detect longitudinal ultrasound waves of various frequencies traveling through randomly packed sapphire spheres subjected to uniaxial compression. The experiments featured in situ X-ray tomography and diffraction measurements of contact fabric, particle kinematics, average per-particle stress tensors, and interparticle forces. The experimentally measured packing configuration and inferred interparticle forces at different sample stresses were used to construct spring networks characterized by Hessian and damping matrices. The ultrasound responses of these network were simulated to investigate the origins of wave velocity, acoustic paths, dispersion, and attenuation. Results revealed that both packing structure and interparticle force heterogeneity played an important role in controlling wave velocity and dispersion, while packing structure alone quantitatively explained most of the observed wave attenuation. This research provides insight into time- and frequency-domain features of wave propagation in randomly packed granular materials, shedding light on the fundamental mechanisms controlling wave velocities, dispersion, and attenuation in such systems.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Expanded Accident Scenarios and Frequency Characteristics for a Fast Spectrum Molten Salt Reactor in an Integrated Energy System

Integration of advanced nuclear reactors to industrial processes in an integrated energy system (IES) has many potential advantages. The coupling of nuclear power to varied industrial processes introduces additional accident scenarios unique to the specific systems included in the IES. Similarly, noise in one system of the IES could potentially create negative effects in the other systems. Because of the limited amount of comparable nuclear IES operating experience, there is added value in scoping the inherent behavior of these systems in response to these transients. Here, this paper investigates the safety-related behavior of an open-source IES dynamic model of a fast spectrum molten salt reactor (MSR) powering a regenerative Rankine cycle for electricity production and a hybrid sulfur (HyS) cycle for hydrogen production. The simulations include a loss of cooling water to the Rankine cycle, a loss of sulfur in the HyS cycle, and the frequency characteristics of the entire IES across six orders of magnitude of frequency. The results show safe behavior in response to the two accident scenarios, wherein a disruption in heat transfer leads to increased salt temperatures and a decrease in reactor power. The frequency analysis shows that at higher frequencies, temperature changes propagating through an IES are increasingly damped to the point of negligibility.

Dunkle, Nicholas [Univ. of Tennessee, Knoxville, T↗

DOE OSTI2024

Explicit simulation of the Brownian rotation of arbitrary shaped aerosol particles using quaternions

The shape of an aerosol particle strongly influences its mass and momentum transfer cross-sections, charging properties, and other physical properties. Here, we present an explicit time-stepping procedure to simulate the rotational Brownian motion of arbitrary shaped aerosol particles by solving Euler’s equation of rotation. A Langevin formulation of the rotation equations is used, wherein Brownian motion due to thermal collisions between a particle and background gas molecules is represented using a stochastic fluctuating torque and fluid resistance is included as a drag torque. To avoid singularities associated with describing the orientation of a shape with Euler angles, we employ a quaternion formulation that leads to first-order stochastic differential equations to describe the evolution of the angular position and angular velocity of a rigid body. We perform all the rotational dynamics calculations in the body-fixed frame of reference attached to the rotating shape whose basis vectors are the normalized eigenvectors of the inertia tensor of the particle. Numerical solutions to rotation under torque-free conditions, damped rotation without Brownian motion, and stochastic rotation for arbitrary shapes are presented and discussed. The presented method enables time-resolved simulation of Brownian rotation for direct comparison with experimentally measured trajectories or statistical measures. The second order accuracy of the used time-stepping procedure places a severe restriction on the timestep that can be used for obtaining accurate results. Animations of presented simulations are included for visualizing rotational motion at various gas pressures. To aid implementation, MATLAB ® codes are also provided. Extension to include translation Brownian motion is straightforward.

Roy, Mrittika↗

DOE OSTI2022

Temperature-dependent spectral emittance of bauxite and silica particle beds

Bauxite and silica particles are candidate materials for solar thermal energy storage at high temperatures. The temperature-dependent emittance of packed beds with bauxite and silica particles was measured using a newly upgraded emissometer at wavelengths 2 μm ≤ λ ≤ 16 μm and temperatures up to ~730 K. The room-temperature emittance was obtained from the measured directional-hemispherical reflectance. A fused silica disc was used to test the emissometer by comparing the measured spectral emittance with the calculated emittance from a fitted Lorentz oscillator model. For the polycrystalline silica particles and the fused silica disc, the measured emittance increases with temperature in the mid-infrared region. The underlying mechanism is interpreted as the temperature-dependent damping coefficient in the Lorentz oscillator model. Two types of bauxite particles with different compositions and sizes were investigated. For λ > 10 μm, the measured emittance at elevated temperatures is higher than that at room temperature. In the region 2 μm < λ < 6 μm, the temperature dependence varies for different types of particles. The total emittance of bauxite particle beds was calculated by spectral integration using Planck’s distribution at the prescribed temperature. As a result, the calculated total emittance is between 0.89 and 0.96, but it does not change monotonically with temperature.

42 ENGINEERING↗

DOE OSTI2024

Do black holes remember what they are made of?

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ($l$ = 2, $m$ = 2) and ($l$ = 2, $m$ = 1) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio q and reduced tidal parameter $\overline{Λ}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Role of gas composition in weakened nonlinear standing wave excitation and improved plasma radial uniformity in very-high-frequency asymmetric capacitive Ar/CF 4 discharges

Abstract The higher harmonics generated by nonlinear sheath motion would enhance the standing wave effect, and thus lead to center-peaked plasma density profile in very-high-frequency (VHF) capacitive discharges. In this work, a nonlinear transmission line (NTL) model introduced in Zhouet al(2021Plasma Sources Sci. Technol.30125017) has been extended, with radial transport of various particles and nonlinear sheath motion into account, to investigate the effects of CF 4 fraction α on the nonlinear standing wave excitation and plasma radial uniformity in VHF (60 MHz) capacitively coupled Ar/CF 4 plasmas at 3 Pa. The results indicate that for pure Ar discharges (i.e. α = 0%), the nonlinearly excited harmonics with short wavelength significantly enhance the electron power absorption at the radial center, resulting in a pronounced central-high plasma density profile. As α increases, the high-order harmonics are gradually damped due to the increase of resistance, as well as the longer wavelength caused by thicker sheath thickness. Thus, the radial profile of the electron absorbed power density shifts from center-peak to edge-high. Besides, at the radial center, the electron density and Ar + ion density decrease with α , CF 3 + ion density shows an increasing trend, while F − ion density initially rises and then decreases. Moreover, the density profiles of all the species become more uniform at higher CF 4 fraction, due to the suppressed nonlinear standing wave excitation and the longer wavelength of the nonlinear harmonics.

Physics↗

DOE OSTI2025

Particle acceleration during classical phase transitions on a spherical lattice

Abstract When compressed, certain lattices undergo phase transitions that may allow nuclei to gain significant kinetic energy. To explore the dynamics of this phenomenon, we develop a methodology to study Coulomb coupled N -body systems constrained to a sphere, as in the Thomson problem. We initialize N total Boron nuclei as point particles on the surface of the sphere, allowing them to equilibrate via Coulomb scattering with a viscous damping term. To simulate a phase transition, we remove N rm particles, forcing the system to rearrange into a new equilibrium. With this model, we consider the Thomson problem as a dynamical system, providing a framework to explore how non-zero temperature affects structural imperfections in Thomson minima. We develop a scaling relation for the average peak kinetic energy attained by a single particle as a function of N and N rm . For certain values of N , we find an order of magnitude energy gain when increasing N rm from 1 to 6. The model may help to design a lattice that maximizes the energy output.

Bachmann, Aidan M. (ORCID:0000000294706404)↗