Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear waves”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Rayleigh–Taylor and Richtmyer–Meshkov instabilities: A journey through scales

Hydrodynamic instabilities such as Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities usually appear in conjunction with the Kelvin–Helmholtz (KH) instability and are found in many natural phenomena and engineering applications. They frequently result in turbulent mixing, which has a major impact on the overall flow development and other effective material properties. This can either be a desired outcome, an unwelcome side effect, or just an unavoidable consequence, but must in all cases be characterized in any model. The RT instability occurs at an interface between different fluids, when the light fluid is accelerated into the heavy. The RM instability may be considered a special case of the RT instability, when the acceleration provided is impulsive in nature such as that resulting from a shock wave. Here in this pedagogical review, we provide an extensive survey of the applications and examples where such instabilities play a central role. First, fundamental aspects of the instabilities are reviewed including the underlying flow physics at different stages of development, followed by an overview of analytical models describing the linear, nonlinear and fully turbulent stages. RT and RM instabilities pose special challenges to numerical modeling, due to the requirement that the sharp interface separating the fluids be captured with fidelity. These challenges are discussed at length here, followed by a summary of the significant progress in recent years in addressing them. Examples of the pivotal roles played by the instabilities in applications are given in the context of solar prominences, ionospheric flows in space, supernovae, inertial fusion and pulsed-power experiments, pulsed detonation engines and Scramjets. Progress in our understanding of special cases of RT/RM instabilities is reviewed, including the effects of material strength, chemical reactions, magnetic fields, as well as the roles the instabilities play in ejecta formation and transport, and explosively expanding flows. The article is addressed to a broad audience, but with particular attention to graduate students and researchers who are interested in the state-of-the-art in our understanding of the instabilities and the unique issues they present in the applications in which they are prominent.

42 ENGINEERING↗
DOE OSTI2026

Hierarchical Testing of a Hybrid Machine Learning‐Physics Global Atmosphere Model

Machine learning (ML)-based models have demonstrated high skill and computational efficiency, often outperforming conventional physics-based models in weather and subseasonal predictions. While prior studies have assessed their fidelity in capturing synoptic-scale atmospheric dynamics, their performance across timescales and under out-of-distribution forcing, such as +3K or +4K uniform-warming forcings, and the sources of biases remain elusive, to establish the model's reliability for Earth science. Here, we design three sets of experiments targeting synoptic-scale phenomena, interannual variability, and out-of-distribution uniform-warming forcings. We evaluate the Neural General Circulation Model (NeuralGCM), a hybrid model integrating a dynamical core with ML-based component, against observations and physics-based Earth system models (ESMs). At the synoptic scale, NeuralGCM captures the evolution and propagation of extratropical cyclones with performance comparable to ESMs. At the interannual scale, when forced by El Niño-Southern Oscillation sea surface temperature (SST) anomalies, NeuralGCM successfully reproduces associated teleconnection patterns but exhibits deficiencies in capturing nonlinear response. Under out-of-distribution uniform-warming forcings, NeuralGCM simulates similar responses in global-average temperature and precipitation and reproduces large-scale tropospheric circulation features similar to those in ESMs. Notable weaknesses include overestimating the tracks and spatial extent of extratropical cyclones, biases in the teleconnected wave train triggered by tropical SST anomalies, and differences in upper-level warming and stratospheric circulation responses to SST warming compared to physics-based ESMs. The causes of these weaknesses were explored. Despite the noted weaknesses, NeuralGCM reproduces responses across experiments reasonably and performs comparably to ESMs. By integrating a dynamical core with ML, NeuralGCM shows potential for developing ML-based ESMs.

global warming↗
DOE OSTI2025

Suppressing trapped-electron-mode-driven turbulence via optimization of three-dimensional shaping

Abstract Turbulent transport driven by trapped electron modes (TEMs) is believed to drive significant heat and particle transport in quasihelically symmetric stellarators. Two three-dimensionally-shaped magnetic configurations with suppressed TEM-driven turbulence were generated through optimization that targeted quasihelical symmetry and the available energy of trapped electrons. Initial equilibria have flux surface shapes with a helically rotating negative triangularity (NT) and positive triangularity (PT). In gyrokinetic simulations, TEMs are suppressed in the reduced-TEM NT and PT configurations, showing that NT does not have the same beneficial turbulence properties over PT as seen in tokamaks. Heat fluxes from TEMs are also suppressed. Without temperature gradients and with a strong density gradient, the most unstable modes at low k y were consistent with toroidal universal instabilities (UIs) in the NT case and slab UIs in the PT case. Nonlinear simulations show that UIs drive substantial heat flux in both the NT and PT configurations. A moderate increase in β halves the heat flux in the NT configuration, while suppressing the heat flux in the PT geometry. Based on the present work, future optimizations aimed at reducing electrostatic drift wave-driven turbulent transport will need to consider UIs if β is sufficiently small.

Duff, J. M. (ORCID:0000000153086694)↗
DOE OSTI2024

Rotational Millimeter-Wave Shoe Scanner Using the Discrete Fourier Transform for Backprojection-Based Image Reconstruction

An active 3D microwave / millimeter-wave shoe scanner was previously developed at the Pacific Northwest National Laboratory (PNNL) using two linear arrays scanned over a rectilinear aperture. The radar system chirps a frequency sweep from 10-40 GHz. These frequencies allow imaging through optically opaque material such as leather, rubber, plastics, and other dielectrics. The system was designed to detect concealed items in the soles of shoes while allowing people to leave their shoes on through a security checkpoint. To shrink the footprint of the system, a new iteration of the design has been developed that scans the two linear arrays over a circular aperture. This new footprint opens the possibility of it being installed in the floor of a cylindrical millimeter-wave body scanner. The backprojection-based multilayer dielectric image reconstruction developed at PNNL can easily handle arbitrary spatial sampling, accommodating the new rotational shoe scanner design. Commonly, the fast Fourier transform (FFT) is used to efficiently compute the range response from the data collected by the system as a preprocessing step to the backprojection algorithm. It was found that converting to range using the discrete Fourier transform (DFT) directly has some advantages over the FFT. For example, nonlinear and non-uniform frequency sweeps can easily be compensated for during the computation of the DFT and only the range bins of interest need to be computed and their spacing can be chosen arbitrarily. Because the range conversion step of the image reconstruction is the fastest part of the process there is very little speed penalty for using the DFT over the FFT and it can even increase the speed of image reconstruction when the ranges of interest are fewer than the total span that is calculated in the FFT.

Millimeter-wave imaging, microwave imaging, shoe s↗
DOE OSTI2020

OC6 Phase I: Investigating the underprediction of low-frequency hydrodynamic loads and responses of a floating wind turbine

Phase I of the OC6 project is focused on examining why offshore wind design tools underpredict the response (loads/motion) of the OC5-DeepCwind semisubmersible at its surge and pitch natural frequencies. Previous investigations showed that the underprediction was primarily related to nonlinear hydrodynamic loading, so two new validation campaigns were performed to separately examine the different hydrodynamic load components. In this paper, we validate a variety of tools against this new test data, focusing on the ability to accurately model the low-frequency loads on a semisubmersible floater when held fixed under wave excitation and when forced to oscillate in the surge direction. However, it is observed that models providing better load predictions in these two scenarios do not necessarily produce a more accurate motion response in a moored configuration.

17 WIND ENERGY↗
DOE OSTI2024

Shock tube experiments on the three-layer Richtmyer–Meshkov instability

Here, a vertical shock tube is used for experiments on the three-layer Richtmyer–Meshkov instability. Two closely spaced membrane-less interfaces are formed by the flow of two different sects of three gases: one with air above CO 2 above SF 6 and the other with helium above air above SF 6 . The lightest of the three gases enters the shock tube at the top of the driven section and flows downward. Conversely, the heaviest gas enters at the bottom of the shock tube and flows upward while the intermediate density gas enters at the middle through porous plates. All three gases are allowed to escape through holes at the layer location, leaving an approximately 30-mm layer of intermediate-density gas suspended between the lightest gas from above and the heaviest gas from below. A single-mode, two-dimensional initial perturbation is then imposed on the lower interface by oscillating the shock tube in the horizontal direction. The flow is visualized by seeding the intermediate gas with particles and illuminating it with a pulsed laser. Image sequences are then captured using high-speed video cameras. Perturbation amplitude measurements are made from the three-layer system and compared with measurements from 2, two-layer systems. It is observed that the presence of the upper, initially flat interface produces a decrease in growth of instability amplitude in the nonlinear phase over an equivalent single-interface configuration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

TRUST Nonlinear Dynamics Testbed Assessment

The following assessment evaluates the efficacy of the control script for carrying out a linear signal to generate a linear mechanical response of the system: the TRUST nonlinear dynamics (TRUSTND) testbed. The hardware has three main components: the controller (NI PXIe 8861) embedded in the National Instruments chassis (NI PXIe 1092), the signal amplifier (The Modal Shop Linear Power Amplifier 2050E09) to amplify the output of a custom LabVIEW script, and a shaker (The Modal Shop Electrodynamic Exciter 2075E) where wave spring specimens live and are tested within the attached aluminum testbed (center mass). Different kinds of signals (single-tone, swept frequency and white noise) were sent through this hardware in the E-1 lab space at TA-53. First, the linearity of the setup was tested by specifying a sine wave in the control script and checking the quality of the oscillations in a wave spring (McMaster-Carr 9714K19 [1]). When passing a sine wave through the amplifier was succeeded, the frequency response of the wave spring was tested with white noise for comparison to finite-element model predictions. The testbed outputs a repeatable, linear response to the input excitation when used with the wave spring. These measurements reduce the uncertainty associated with the testbed and test procedure to better characterize uncertainty associated with material behavior (SX358 foam samples) in future tests

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Convergence of sum-up rounding schemes for cloaking problems governed by the Helmholtz equation

In this work, we consider the problem of designing a cloak for waves described by the Helmholtz equation from an integer programming point of view. The problem can be modeled as a PDE-constrained optimization problem with integer-valued control inputs that are distributed in the computational domain. A first-discretize-then-optimize approach results in a large-scale mixed-integer nonlinear program that is in general intractable because of the large number of integer variables that arise from the discretization of the domain. Instead, we propose an efficient algorithm that is able to approximate the local infima of the underlying nonconvex infinite-dimensional problem arbitrarily close without the need to solve the discretized finite-dimensional integer programs to optimality. We optimize only the continuous relaxations of the approximations for local minima and then apply the sum-up rounding methodology to obtain integer-valued controls. If the solutions of the discretized continuous relaxations converge to a local minimizer of the continuous relaxation, then the resulting discrete-valued control sequence converges weakly \(^*\) in \(L^\infty\) to the same local minimizer. These approximation properties follow under suitable refinements of the involved discretization grids. Our results use familiar concepts arising from the analytical properties of the underlying PDE and complement previous results, derived from a topology optimization point of view.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2020

Investigating the use of field solvers for simulating classical systems

We explore the use of field solvers as approximations of classical Vlasov-Poisson systems. This correspondence is investigated in both electrostatic and gravitational contexts. We demonstrate the ability of field solvers to be excellent approximations of problems with cold initial condition into the nonlinear regime. We also investigate extensions of the Schrödinger-Poisson system that employ multiple stacked cold streams, and the von Neumann–Poisson equation as methods that can successfully reproduce the classical evolution of warm initial conditions. We then discuss how appropriate simulation parameters need to be chosen to avoid interference terms, aliasing, and wave behavior in the field solver solutions. Finally, we present a series of criteria clarifying how parameters need to be chosen in order to effectively approximate classical solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Extreme-value statistics in nonlinear optics

We show that, although nonlinear optics may give rise to a vast multitude of statistics, all these statistics converge, in their extreme-value limit, to one of a few universal extreme-value statistics. Specifically, in the class of polynomial nonlinearities, such as those found in the Kerr effect, weak-field harmonic generation, and multiphoton ionization, the statistics of the nonlinear-optical output converges, in the extreme-value limit, to the exponentially tailed, Gumbel distribution. Exponentially growing nonlinear signals, on the other hand, such as those induced by parametric instabilities and stimulated scattering, are shown to reach their extreme-value limits in the class of the Fréchet statistics, giving rise to extreme-value distributions (EVDs) with heavy, manifestly nonexponential tails, thus favoring extreme-event outcomes and rogue-wave buildup.

Zheltikov, Aleksei M. (ORCID:0000000291380576)↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The polarized-signal density matrix: A practical way to recover molecular frame information from isotropic samples

We present a novel approach to model ultrafast time-dependent nonlinear optical polarization sensitive signals emitted from randomly oriented molecules. By projecting the laboratory-frame analyzer polarization axis into the molecular-frame and linking that axis with the density matrix through a tensor product, we demonstrate an approach to find a specific molecular orientation that yields a good approximation to simulated four-wave mixing signals produced by the same model but with averaging over molecular orientation.

Thurston, Richard L↗
DOE OSTI2025

Single-shot electrical detection of short-wavelength magnon pulse transmission in a magnonic thin-film waveguide

Abstract The advance of magnon spintronics requires understanding of time-domain magnon pulse transmission in order to develop high-speed information processing protocols. In this work, we demonstrate single-shot electrical detection of narrow-band magnon pulse transmission in a yttrium iron garnet thin-film delay line. The high signal-to-background ratio of magnon transmission band allows us to directly probe the magnon transmission electrically using a fast oscilloscope and to study its spectral evolution using Fast Fourier Transform (FFT) of the time-domain transmitted signal. At elevated input power, we show a magnon transmission reduction and a spectral distortion, which can be understood by the nonlinear magnon excitation in the transmission band defined by the antenna geometry. In addition, we also find that the higher- (lower-) frequency magnon spectral component exhibits a lower (higher) magnon group velocity, showing a dispersion agreeing with the Damon-Eshbach dependence. Our results provide important guidance of magnon pulse engineering for their applications in spin wave computing and coherent magnon information processing.

Song, Moojune↗
DOE OSTI2024

A semi-analytical model of RF condensation that can handle localized power depositions

A nonlinear effect, RF (radio frequency) condensation, can be used to facilitate RF stabilization of magnetic islands. Previously studied semi-analytical models for RF condensation are suited mainly for broad deposition profiles and are unable to handle power depositions that are localized in the interior of a magnetic island. The power deposition profile of electron cyclotron waves, in particular, can be highly localized near a cyclotron resonance. Here, a model is developed that can handle both localized profiles and broad profiles. This allows a comparison of RF condensation for narrow vs broad deposition profiles, and it allows a study of the dependence of RF condensation of localized deposition profiles on key parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Negative Refraction in Time-Varying Strongly Coupled Plasmonic-Antenna–Epsilon-Near-Zero Systems

Time-varying metasurfaces are emerging as a powerful instrument for the dynamical control of the electromagnetic properties of a propagating wave. Here we demonstrate an efficient time-varying metasurface based on plasmonic nano-antennas strongly coupled to an epsilon-near-zero (ENZ) deeply subwavelength film. The plasmonic resonance of the metal resonators strongly interacts with the optical ENZ modes, providing a Rabi level spitting of ~30%. Optical pumping at frequency ω induces a nonlinear polarization oscillating at 2ω responsible for an efficient generation of a phase conjugate and a negative refracted beam with a conversion efficiency that is more than 4 orders of magnitude greater compared to the bare ENZ film. In this work, the introduction of a strongly coupled plasmonic system therefore provides a simple and effective route towards the implementation of ENZ physics at the nanoscale.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

EPR of Photoexcited Triplet-State Acceptor Porphyrins

The photoexcited triplet states of porphyrin architectures are of significant interest in a wide range of fields including molecular wires, nonlinear optics, and molecular spintronics. Electron paramagnetic resonance (EPR) is a key spectroscopic tool in the characterization of these transient paramagnetic states singularly well suited to quantify spin delocalization. Previous work proposed a means of extracting the absolute signs of the zero-field splitting (ZFS) parameters, D and E, and triplet sublevel populations by transient continuous wave, hyperfine measurements, and magnetophotoselection. Here, we present challenges of this methodology for a series of mesoperfluoroalkyl-substituted zinc porphyrin monomers with orthorhombic symmetries, where interpretation of experimental data must proceed with caution and the validity of the assumptions used in the analysis must be scrutinized. The EPR data are discussed alongside quantum chemical calculations, employing both DFT and CASSCF methodologies. Despite some success of the latter in quantifying the magnitude of the ZFS interaction, the results clearly provide motivation to develop improved methods for ZFS calculations of highly delocalized organic triplet states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Pump-depletion dynamics and saturation of stimulated Brillouin scattering in shock ignition relevant experiments

As an alternative inertial confinement fusion scheme, shock ignition requires a strong converging shock driven by a high-intensity laser pulse to ignite a pre-compressed fusion capsule. Understanding nonlinear laser-plasma instabilities is crucial to assess and improve the laser-shock energy coupling. Recent experiments conducted on the OMEGA EP laser facility have for the first time demonstrated that such instabilities can 100% deplete the first 0.5 ns of the high-intensity laser. Here, analyses of the observed laser-generated blast wave suggest that this pump-depletion starts at 0.02 critical density and progresses to 0.1-0.2 critical density, which is also confirmed by the time-resolved stimulated Raman backscattering spectra. The pump-depletion dynamics can be explained by the breaking of ion-acoustic waves in stimulated Brillouin scattering. Such pump-depletion would inhibit the collisional laser energy absorption but may benefit the generation of hot electrons with moderate temperatures for electron shock ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗