Engineering Papers⌕ Search

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.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Plasmon-enhanced optical nonlinearity in graphene nanomeshes

Using the density-matrix formalism, here we show that graphene nanomeshes (GNMs)—graphene sheets patterned with antidots—have large plasmon-enhanced nonlinear optical response. GNMs can be designed to behave as quasi-one-dimensional plasmonic crystals in which plasmons with large propagation lengths are efficiently excited. The associated third-order Kerr and third-harmonic-generation susceptibility can be as high as 10 –7 and 10 –9 m 2 V –2 , respectively, over the mid-to-near-infrared frequency range. Furthermore, carrier-density tuning in GNMs can flip the propagation direction of plasmonic waves and enables bidirectional switching of optical signals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Ultrafast-induced coherent acoustic phonons in the two-dimensional magnet CrSBr

Magnetism in two-dimensional (2D) van der Waals (vdW) crystals offers promising new directions for low-dimensional physics and devices. In this work, mega-electron volt (MeV) ultrafast electron diffraction was employed to investigate the ultrafast atomic dynamics of a novel, 2D vdW magnetic single-crystal CrSBr. Femtosecond (fs) optical pump pulses excited non-equilibrium atomic displacements shown to be coherent acoustic phonons (CAPs). Phonon frequencies were extracted by analyzing oscillations of different Bragg peak (BP) intensities and were determined to be GHz acoustic disturbances that propagated as strain waves. Phonon modes exhibit anisotropy with respect to the a and b crystal axes. Subharmonic phonon frequencies were also observed, and this provided a signature of nonlinear oscillatory coupling between the laser-induced pumping phonon frequency and secondary phonon frequencies. Thus, CrSBr was found to serve as a nonlinear phononic frequency converter. The ultrafast time dependence of the Bragg intensity was simulated by incorporating an oscillating deviation parameter ansatz into expressions for the dynamical scattering intensity yielded excellent modeling of the ultrafast structural dynamics of the photo-excited 2D crystal. Our work provides a foundation for exploring how fs light pulses can influence phonon dynamics in materials with strong spin-lattice coupling. These results suggest that CAPs can match the magnon frequencies and show the promise of CrSBr for use in optical-to-microwave transducers and phononic devices.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Enhanced stimulated Raman scattering during intense laser propagation

Stimulated Raman scattering is ubiquitous in many high-intensity laser environments. Parametric four-wave mixing between the pump and Raman sidebands can affect the Raman gain, but stringent phase matching requirements and strongly nonlinear dynamics obscure clear understanding of its effects at high laser powers. Here we investigate four-wave mixing in the presence of strong self-focusing and weak ionization at laser powers above the Kerr critical power. Theoretical analysis shows that the plasma generated at focus naturally leads to phase matching conditions suitable for enhanced Raman gain, almost without regard to the initial phase mismatch. Multidimensional nonlinear optical simulations with multiphoton and collisional ionization confirm the enhancement and suggest that it may lead to significantly higher Raman losses in some high-intensity laser environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Simulations of non-integer upconversion in resonant six-wave scattering

Resonant upconversion through a sixth order relativistic nonlinearity resulting in a unique resonance was recently proposed [Malkin and Fisch, Phys. Rev. E 108, 045208 (2023)]. The high order resonance is a unique non-integer multiple of a driving pump frequency resulting in a frequency upshift by a factor of ≈3.73. We demonstrate the presence, unique requirements, and growth of this mode numerically. Through tuning waves to high amplitude, in a mildly underdense plasma, the six-photon process may grow more than other non-resonant but lower order processes. The growth of the high frequency mode remains below the nonlinear growth regime. However, extending current numerical results to more strongly coupled resonances with longer pulse propagation distances suggests a pathway to significant upconversion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Detecting Large Explosions With Machine Learning Models Trained on Synthetic Infrasound Data

Explosions produce low-frequency acoustic (infrasound) waves capable of propagating globally, but the spatio-temporal variability of the atmosphere makes detecting events difficult. Machine learning (ML) is well-suited to identify the subtle and nonlinear patterns in explosion infrasound signals, but a previous lack of ground-truth data inhibited training of generalized models. We introduce a physics-based method that propagates infrasound sources through realistic atmospheres to create 28,000 synthetic events, which are used to train ML classifiers. A simple artificial neural network and modern temporal convolutional network discriminate synthetic events from background noise with >90% accuracy and, more importantly, successfully identify the majority of real-world explosion signals recorded during the Humming Road Runner experiment. ML models trained entirely on physics-based synthetics advance explosion detection capabilities and make ML more viable to related fields lacking training data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2022

Integrating a ponderomotive guiding center algorithm into a quasi-static particle-in-cell code based on azimuthal mode decomposition

High fidelity modeling of plasma based acceleration (PBA) requires the use of three dimensional, fully nonlinear, and kinetic descriptions based on the particle-in-cell (PIC) method. In PBA an intense particle beam or laser (driver) propagates through a tenuous plasma whereby it excites a plasma wave wake. Three-dimensional PIC algorithms based on the quasi-static approximation (QSA) have been successfully applied to efficiently model the interaction between relativistic charged particle beams and plasma. In a QSA PIC algorithm, the plasma response to a charged particle beam or laser driver is calculated based on forces from the driver and self-consistent forces from the QSA form of Maxwell's equations. These fields are then used to advance the charged particle beam or laser forward by a large time step. Since the time step is not limited by the regular Courant-Friedrichs-Lewy (CFL) condition that constrains a standard 3D fully electromagnetic PIC code, a 3D QSA PIC code can achieve orders of magnitude speedup in performance. Recently, a new hybrid QSA PIC algorithm that combines another speedup technique known as an azimuthal Fourier decomposition has been proposed and implemented. This hybrid algorithm decomposes the electromagnetic fields, charge and current density into azimuthal harmonics and only the Fourier coefficients need to be updated, which can reduce the algorithmic complexity of a 3D code to that of a 2D code. Modeling the laser-plasma interaction in a full 3D electromagnetic PIC algorithm is very computationally expensive due the enormous disparity of physical scales to be resolved. In the QSA the laser is modeled using the ponderomotive guiding center (PGC) approach. We describe how to implement a PGC algorithm compatible for the QSA PIC algorithms based on the azimuthal mode expansion. Here this algorithm permits time steps orders of magnitude larger than the cell size and it can be asynchronously parallelized. Details on how this is implemented into the QSA PIC code that utilizes an azimuthal mode expansion, QPAD, are also described. Benchmarks and comparisons between a fully 3D explicit PIC code (OSIRIS), as well as a few examples related to laser wakefield acceleration, are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Turbulent Regimes in Collisions of 3D Alfvén-wave Packets

Using three-dimensional gyrofluid simulations, we revisit the problem of Alfvén-wave (AW) collisions as building blocks of the Alfvénic turbulent cascade and their interplay with magnetic reconnection at magnetohydrodynamic (MHD) scales. Depending on the large-scale value of the nonlinearity parameter χ 0 (the ratio between the AW linear propagation time and nonlinear turnover time), different regimes are observed. For strong nonlinearities (χ 0 ~ 1), turbulence is consistent with a dynamically aligned, critically balanced cascade—fluctuations exhibit a scale-dependent alignment $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/4}$, resulting in a ${k}_{\perp }^{-3/2}$ spectrum and ${k}_{\parallel }\propto {k}_{\perp }^{1/2}$ spectral anisotropy. At weaker nonlinearities (small χ 0 ), a spectral break marking the transition between a large-scale weak regime and a small-scale ${k}_{\perp }^{-11/5}$ tearing-mediated range emerges, implying that dynamic alignment occurs also for weak nonlinearities. At χ 0 < 1 the alignment angle ${\theta }_{{k}_{\perp }}$ shows a stronger scale dependence than in the χ 0 ~ 1 regime, namely $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1/2}$ at χ 0 ~ 0.5, and $\sin {\theta }_{{k}_{\perp }}\propto {k}_{\perp }^{-1}$ at χ 0 ~ 0.1. Dynamic alignment in the weak regime also modifies the large-scale spectrum, scaling approximately as ${k}_{\perp }^{-3/2}$ for χ 0 ~ 0.5 and as ${k}_{\perp }^{-1}$ for χ 0 ~ 0.1. A phenomenological theory of dynamically aligned turbulence at weak nonlinearities that can explain these spectra and the transition to the tearing-mediated regime is provided; at small χ 0 , the strong scale dependence of the alignment angle combines with the increased lifetime of turbulent eddies to allow tearing to onset and mediate the cascade at scales that can be larger than those predicted for a critically balanced cascade by several orders of magnitude. Such a transition to tearing-mediated turbulence may even supplant the usual weak-to-strong transition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Nonlinear wavepacket dynamics in proximity to a stationary inflection point

A stationary inflection point (SIP) in the Bloch dispersion relation of a periodic waveguide is an exceptional point degeneracy where three Bloch eigenmodes coalesce, forming the so-called frozen mode with a divergent amplitude and vanishing group velocity of its propagating component. For this work, we have developed a theoretical framework to study the time evolution of wavepackets centered at an SIP. Analysis of the evolution of the statistical moment distribution of linear pulses shows a strong deviation from the conventional ballistic wavepacket dynamics in dispersive media. The presence of nonlinear interactions dramatically changes the situation, resulting in a mostly ballistic propagation of nonlinear wavepackets with the speed and even the direction of propagation essentially dependent on the wavepacket amplitude. Such a behavior is unique to nonlinear wavepackets centered at an SIP and can be used for the realization of a novel family of beam power routers for classical waves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Unidirectional amplification in the frozen mode regime enabled by a nonlinear defect

A stationary inflection point (SIP) is a spectral singularity of the Bloch dispersion relation ω(k) of a periodic structure where the first and the second derivatives of ω with respect to k vanish. An SIP is associated with a third-order exceptional point degeneracy in the spectrum of the unit-cell transfer matrix, where there is a collapse of one propagating and two evanescent Bloch modes. At the SIP frequency, the incident wave can be efficiently converted into the frozen mode with greatly enhanced amplitude and vanishing group velocity. This can be very attractive for applications, including light amplification. Due to its non-resonant nature, the frozen mode regime (FMR) has fundamental advantages over common cavity resonances. Furthermore, we propose, a novel, to the best of our knowledge, scheme for FMR-based unidirectional amplifiers by leveraging a tailored amplification/attenuation mechanism and a single nonlinear defect. The defect breaks the directional symmetry of the periodic structure and enables nonlinearity-related unidirectional amplification/attenuation in the vicinity of the SIP frequency. We demonstrate the robustness of the amplification mechanism to local impurities and parasitic nonlinearity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 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 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 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↗