SEARCH · Engineering Papers

Results for “nonlinear wave propagation”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Dynamic Carrier Modulation via Nonlinear Acoustoelectric Transport in van der Waals Heterostructures

Dynamically manipulating carriers in van der Waals heterostructures could enable solid-state quantum simulators with tunable lattice parameters. A key requirement is the formation of deep potential wells to reliably trap excitations. Here, we report the observation of nonlinear acoustoelectric transport and dynamic carrier modulation in boron nitride-encapsulated graphene devices coupled to intense surface acoustic waves (SAWs) on LiNbO 3 substrates. SAWs generate strong acoustoelectric current densities ( J AE ), transitioning from linear to nonlinear regimes with increasing SAW intensity. In the nonlinear regime, periodic carrier (electrons, holes, or their mixtures) stripes emerge. Using counter-propagating SAWs, we create standing SAWs (SSAWs) to dynamically manipulate charge distributions without static gates. The saturation of J AE , attenuation transitions, and tunable resistance peaks confirms strong carrier localization. Finally, these results establish SAWs as a powerful tool for controlling carrier dynamics in two-dimensional (2D) materials, paving the way for the development of time-dependent quantum systems and acoustic lattices for quantum simulation.

2D materials↗

DOE OSTI2026

Magnonic spontaneous oscillation induced by parametric pumping

Spontaneous dynamic systems have attracted significant attention for their rich underlying physics such as phase-locking and synchronization. In this work, we report a new mechanism of generating magnetic spontaneous oscillation via parametric pumping. By applying a pump tone to excite propagating spin waves in a yttrium iron garnet delay line, four-wave mixing converts the pump mode into two phase-autonomous propagating magnon modes, i.e. a spontaneous mode with nearly twice the wavenumber of the pump mode and an idler mode with nearly zero wavenumber. This allows us to reliably generate ultrasharp spin wave dynamics with broad frequency tunability from the pump and magnetic field. We show that the spontaneous mode can be phase-locked to a probe tone, similar to an auto-oscillator. Furthermore, the spontaneous dynamics can be used to implement a high-gain magnonic parametric amplifier with a gain up to 40 dB. Our results open a new avenue for studying nonlinear magnonics and synchronization physics in propagating magnon geometry and for developing new magnonic devices.

Li, Yi↗

DOE OSTI2022

Asymptotics of transverse momentum broadening in dense QCD media

We study the asymptotic behaviour of the transverse momentum broadening distribution of an energetic quark or gluon propagating through dense QCD matter, in the large system size L L limit, taking into account radiative corrections in the double logarithmic approximation. Thanks to a connection between the evolution of the jet quenching parameter \hat{q} q ̂ and the formation of traveling wave fronts in nonlinear physics, we obtain a formula for the L L dependence of the characteristic transverse momentum scale Q_s Q s of the distribution valid up to terms of order 1/\ln(L) 1 / ln ( L ) . We briefly discuss the physical implications of this formula for jet quenching and small- x x phenomenology.

Caucal, Paul↗

DOE OSTI2022

An improved framework for the dynamic likelihood filtering approach to data assimilation

Here, we propose improvements to the Dynamic Likelihood Filter (DLF), a Bayesian data assimilation filtering approach, specifically tailored to wave problems. The DLF approach was developed to address the common challenge in the application of data assimilation to hyperbolic problems in the geosciences and in engineering, where observation systems are sparse in space and time. When these observations have low uncertainties, as compared to model uncertainties, the DLF exploits the inherent nature of information and uncertainties to propagate along characteristics to produce estimates that are phase aware as well as amplitude aware, as would be the case in the traditional data assimilation approach. Along characteristics, the stochastic partial differential equations underlying the linear or nonlinear stochastic dynamics are differential equations. This study focuses on developing the explicit challenges of relating dynamics and uncertainties in the Eulerian and Lagrangian frames via dynamic Gaussian processes. It also implements the approach using the ensemble Kalman filter (EnKF) and compares the DLF approach to the conventional one with respect to wave amplitude and phase estimates in linear and nonlinear wave problems. Numerical comparisons show that the DLF/EnKF outperforms the EnKF estimates, when applied to linear and nonlinear wave problems. This advantage is particularly noticeable when sparse, low uncertainty observations are used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Dynamics of the water-plasma interface in various discharge modes of atmospheric-pressure plasmas

This paper presents an experimental study of the dynamics of the water-plasma interface during the interaction high-voltage (HV) atmospheric pressure discharge with the water surface. Under the influence of the voltage applied to the pin type electrode, discharges are formed at the air layer between the sharp-tip HV electrode and deionized water. Three discharge regimes were identified from synchronized electrical waveforms and imaging: (i) a weak linear regime at 3–10.6 kV, (ii) a branching streamer regime at the maximum applied voltage of 12.6 kV, and (iii) a continuous arc-like regime occurring during the transition to lower voltage and higher current (e.g., U ≈ 3.2 kV, I ≈ 4.07 mA). High-speed shadowgraph images showed a symmetric interfacial cavity whose depth increased nonlinearly with voltage from h 0 ≈ 0.1 mm at 3 kV to h 0 ≈ 2.7 mm at 10.6 kV, reaching a maximum depth of h 0 ≈ 5.9 mm at 12.6 kV, while the cavity disappeared in the continuous-channel regime and was replaced by outward-propagating wave-like motion. The deformation of the water surface during the discharge-water interaction is governed by the balance between electric field forces, surface tension, gravitational forces, and electrohydrodynamic forces, whose relative contributions vary with the applied voltage. This effect was explained by quantitatively estimating the magnitudes of the acting forces and establishing the force balance. Optical emission spectroscopy in the continuous-channel regime indicated air plasma signatures and yielded a gas temperature of approximately 400 K, while prolonged operation (≈ 20 min) increased the water temperature to ~ 70 °C, reduced the water-layer thickness from 6 mm to 3 mm, and decreased pH from 7 to 4. These regime-resolved, quantitative results clarify how the dominant interfacial forcing shifts with discharge mode and provide a mechanistic basis for controlling plasma-water interactions in atmospheric-pressure applications.

atmospheric pressure plasma↗

DOE OSTI2022

Differences in Tropical Rainfall in Aquaplanet Simulations With Resolved or Parameterized Deep Convection

Abstract This study investigates the effects of resolved deep convection on tropical rainfall and its multi‐scale variability. A series of aquaplanet simulations are analyzed using the Model for Prediction Across Scales‐Atmosphere with horizontal cell spacings from 120 to 3 km. The 3‐km experiment uses a novel configuration with 3‐km cell spacing between 20°S and 20°N and 15‐km cell spacing poleward of 30°N/S. A comparison of those experiments shows that resolved deep convection yields a narrower, stronger, and more equatorward intertropical convergence zone, which is supported by stronger nonlinear horizontal momentum advection in the boundary layer. There is also twice as much tropical rainfall variance in the experiment with resolved deep convection than in the experiments with parameterized convection. All experiments show comparable precipitation variance associated with Kelvin waves; however, the experiment with resolved deep convection shows higher precipitation variance associated with westward propagating systems. Resolved deep convection also yields at least two orders of magnitude more frequent heavy rainfall rates (>2 mm hr −1 ) than the experiments with parameterized convection. A comparison of organized precipitation systems demonstrates that tropical convection organizes into linear systems that are associated with stronger and deeper cold pools and upgradient convective momentum fluxes when convection is resolved. In contrast, parameterized convection results in more circular systems, weaker cold pools, and downgradient convective momentum fluxes. These results suggest that simulations with parameterized convection are missing an important feedback loop between the mean state, convective organization, and meridional gradients of moisture and momentum.

Rios‐Berrios, Rosimar↗

DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Fine Structure of Magnetospheric Magnetosonic Waves: 1. Elementary Rising‐Tone Emissions Within Individual Harmonic

Abstract The present study uncovers the fine structures of magnetosonic waves by investigating the EFW waveforms measured by Van Allen Probes. We show that each harmonic of the magnetosonic wave may consist of a series of elementary rising‐tone emissions, implying a nonlinear mechanism for the wave generation. By investigating an elementary rising‐tone magnetosonic wave that spans a wide frequency range, we show that the frequency sweep rate is likely proportional to the wave frequency. We studied compound rising‐tone magnetosonic waves, and found that they typically consist of multiple harmonics in the source region, and may gradually become continuous in frequency as they propagate away from source. Both elementary and compound rising‐tone magnetosonic waves last for ∼1 min which is close to the bounce period of the ring proton distribution, but their relation is not fully understood.

Li, Jinxing↗

DOE OSTI2026

Stability and Characteristics of Lower-hybrid Drift Waves: Dependence on Electron Beta and Cross-field Relative Drift

Lower-hybrid drift waves (LHDWs) are frequently observed microinstabilities in both space and laboratory plasmas. Despite decades of study, the relationship between electrostatic (ES-LHDW) and electromagnetic (EM-LHDW) variants and the plasma parameters controlling their stability remains unclear. Here, we systematically examine LHDW behavior by solving the local linear dispersion relation over a wide range of plasma and field conditions. Our results demonstrate that ES-LHDWs and EM-LHDWs are not distinct modes but rather two different regimes of the same drift wave whose character evolves smoothly with electron beta (β e ) and the cross-field electron drift velocity relative to ions, normalized to the ion sound speed (u 0x /C s ). The nature of the waves changes from electrostatic to electromagnetic when β e increases. Growth rates increase with u 0x /C s but decrease with β e , while the most unstable wavelength remains nearly universal, with kρ e ∼ 0.8 (k is the magnitude of the wave vector and ρ e is the electron gyroradius). We further present quasi-linear estimates of nonlinear saturation properties, including energy partition among electric fields, magnetic fields, and particle kinetic responses. We show that ES-LHDWs reach higher electric-field saturation amplitudes, whereas EM-LHDWs generate strong magnetic perturbations and parallel electric fields that may enable efficient particle heating. Comparisons with the classical model reveal that retaining electromagnetic effects is essential for accurate predictions of frequency, growth rate, and the propagation angle. These findings provide a unified framework for understanding LHDWs across diverse collisionless plasma environments, including current sheets of magnetic reconnection, shear layers, collisionless shocks, and boundary regions.

Solar coronal waves↗

DOE OSTI2022

Linear and Nonlinear Kelvin–Helmholtz Instability and Magnetohydrodynamic Wave Emission in Sheared Astrophysical Plasma Flows

The evolution of the Kelvin–Helmholtz instability (KHI) and magnetohydrodynamic (MHD) wave emission is investigated at shear-flow boundaries of magnetized plasmas. While MHD wave emission has been suggested to be only possible during the nonlinear stages, we find that there is also significant wave emission during the KHI's linear stages. These emitted MHD waves may have stronger impacts than KHI surface waves since they can act to transport energy away from the local region of the shear flow. The removal of energy from the shear-flow region, instead of just the local redistribution considered in previous studies, and its propagation away from the interface could have major implications for the evolution of astrophysical objects characterized by fast plasma flow shears.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory

Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.

Drops and Bubbles: Bubble dynamics↗

DOE OSTI2024

High-Fidelity Numerical Wave Tank Verification & Validation Study: Wave Generation Through Paddle Motion: Preprint

This paper presents a numerical benchmark study of wave propagation due to a paddle motion using different high-fidelity numerical models, which are capable of replicating the nearly actual physical wave tank testing. A full time series of the measured wave generation paddle motion which was used to generate wave propagation in the physical wave tank will be utilized in each of the models contributed by IEA OES Task 10's participants, which includes both computational fluid dynamics (CFD) and smooth hydrodynamic particle (SPH). The high-fidelity simulations of the physical wave testcase will allow for the evaluation of the initial transient effects from wave ramp-up and its evolution in the wave tank over time for two representative regular waves with varying levels of nonlinearity. A couple of interesting metrics like the predicted wave surface elevation at select wave probes, wave period, and phase-shift in time will be assessed to evaluate the relative accuracy of numerical models versus experimental data within specified time intervals. These models will serve as a guide for modelers in the wave energy community and provide a base case to allow further and more detailed numerical modeling of the fixed Kramer Sphere Cases under wave excitation force wave tank testing.

HYDRO ENERGY,TIDAL AND WAVE POWER↗

DOE OSTI2022

Visualization and quantum control of light-accelerated condensates by terahertz multi-dimensional coherent spectroscopy

Abstract Characterizing and controlling high-order correlation of quantum systems is key for developing quantum devices and switching technologies. Although conventional static and ultrafast spectroscopy gives access to collective excitations characterizing quantum states, more exotic correlations cannot be easily separated from other contributions. Here we develop density matrix simulations to show that seventh-order-wave-mixing peaks with distinct temperature and field dependences in two-dimensional terahertz nonlinear spectra reveal light-induced correlations in non-equilibrium superconducting states. Above critical terahertz driving, these emerging peaks split from conventional peaks along the second axis introduced by pump-probe relative phase in two-dimensional frequency space. They are photo-generated by correlations between two-photon fluctuations and interacting quasi-particle and quasi-particle/Higgs superconductor excitations. By photo-inducing persistent symmetry breaking via light-wave propagation, we also demonstrate seventh-order-wave-mixing sensing of Higgs collective modes. Our theory suggests to use multi-dimensional spectroscopy for quantum sensing of light-driven superconductivity and paves a path for quantum operations by few-cycle-THz-periodic photocurrent modulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Propagation of partially spatially coherent laser beams in instantaneous Kerr media

The propagation of intense, partially spatially coherent laser beams in a medium with instantaneous third-order susceptibility is studied analytically and numerically. For sufficiently high power relative to that required for nonlinear self-focusing, the propagation initially proceeds in two stages. In the first stage, spatial coherence builds up, and in the second stage, the number of speckles reduces. Once the degree of coherence is sufficiently high, whole-beam self-focusing occurs. The beam power is mostly confined within the initial spot radius. Two analytical approaches for describing the evolution of the beam are presented. The method of moments leads to an analytical solution for the rms spot radius that is in excellent agreement with simulations. This method does not require any knowledge of the field statistics beyond the initial conditions and provides no information about the evolution of the individual speckles. The other approach employs a self-similar solution for the second-order coherence function of the field and assumes that the fourth-order coherence function is factorizable and obeys complex circular Gaussian random statistics. The latter method also leads to an analytical expression for the spot radius, but its predictions for the qualitative evolution of the speckles disagree with wave-optics simulations.

lasers↗

DOE OSTI2022

Numerical modeling of pedestal stability and broadband turbulence of wide-pedestal QH-mode plasmas on DIII-D

The wide-pedestal quiescent high confinement mode discovered on DIII-D in recent years is a stationary and quiescent H-mode (QH-mode) with the pedestal width exceeding EPED prediction by at least 25%. Its characteristics, such as low rotation, high energy confinement and edge localized mode-free operation, make it an attractive operation mode for future reactors. Linear and nonlinear simulations using BOUT++ reduced two fluid MHD models and awere carried out to investigate the bursty broadband turbulence often observed in the edge of wide-pedestal QH-mode plasmas. Two kinds of MHD-scale instabilities in different spatial locations within the pedestal were found in the simulations: one mild peeling–ballooning (PB) mode γ PB < 0.04ω A ) located near the minimum in E r</:sub> well propagating in ion diamagnetic drift direction; and one drift-Alfvén wave locates at smaller radius compared to E r</:sub> well propagating in the electron diamagnetic drift direction and unstable only when the parallel electron dynamics is included in the simulation. The coupling between drift wave and shear Alfvén wave provides a possible cause of the experimentally observed local profile flattening in the upper-pedestal. The rotation direction, mode location, as well as the wavenumber of these two modes from BOUT++ simulations agree reasonably well with the experimental measurements, while the lack of quantitative agreement is likely due to the lack of trapped electron physics in current fluid model. This work presents improved physics understanding of the pedestal stability and turbulence dynamics for wide-pedestal QH-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Multimodal nonlinear optical imaging via evanescent wave excitation

A system for multimodal nonlinear optical imaging is provided. Each mode uses a high NA objective to cause total internal reflection excitation at a sample-substrate interface. The system has a femtosecond oscillator to generate pulses used for two beams. The objective receives at least one beam, redirects the received at least one beam through a dielectric substrate to cause the TIR and produces corresponding evanescent waves in a portion of the sample adjacent to the sample-substrate interface, and collects a backward-propagating beam of pulses of responsive light. The portion of the sample illuminated by the evanescent waves emits responsive light. Different modes or combinations of the distinct modalities may be selected to access complementary chemical and structural information for various chemical species near the sample-substrate interface. Each mode may have mode-specific control such as selective beam blocking, power ratios and filtering.

Doughty, Benjamin L.↗

DOE OSTI2025

Time-refraction near the critical angle and angular streaking with attosecond resolution [Invited]

We study the time-refraction of optical waves as they propagate through a time-varying slab and find the frequency shift associated with the abrupt change in the refractive index. We discover a geometrical mechanism responsible for an ultrafast blueshift followed by a redshift, occurring while the refractive index is increased monotonically. Furthermore, we find high sensitivity to incident angles near the traditional critical angle and unravel a new critical angle that varies in time. Finally, we introduce an innovative concept of a streak camera design that utilizes ultrafast nonlinearities in time-varying media and offers attosecond resolution capabilities.

Fruhling, Colton (ORCID:0000000282968802)↗