Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear normal modes”

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.

48 records · Page 3
All sourcesNASA NTRSDOE OSTI
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 OSTI2023

ETG turbulence in a tokamak pedestal

This paper explores the fundamental characteristics of electron-temperature-gradient (ETG)-driven turbulence in the tokamak pedestal. The extreme gradients in the pedestal produce linear instabilities and nonlinear turbulence that are distinct from the corresponding ETG phenomenology in the core plasma. The linear system exhibits multiple (greater than ten) unstable eigenmodes at each perpendicular wave vector, representing different toroidal and slab branches of the ETG instability. Proper orthogonal decomposition of the nonlinear fluctuations reveals no clear one-to-one correspondence between the linear and nonlinear modes for most wave vectors. Moreover, nonlinear frequencies deviate strongly from those of the linear instabilities, with spectra peaking at positive frequencies, which is opposite the sign of the ETG instability. The picture that emerges is one in which the linear properties are preserved only in a narrow range of k-space. Outside this range, nonlinear processes produce strong deviations from both the linear frequencies and eigenmode structures. This is interpreted in the context of critical balance, which enforces alignment between the parallel scales and fluctuation frequencies. We also investigate the nonlinear saturation processes. We observe a direct energy cascade from the injection scale to smaller scales in both perpendicular directions. However, in the bi-normal direction, there is also nonlocal inverse energy transfer to larger scales. Neither streamers nor zonal flows dominate the saturation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum Fluctuations and Multifractally enhanced Superconductivity in Disordered Thin Films

The interplay between electron–electron interactions and weak localization (or anti-localization) phenomena in two-dimensional systems can significantly enhance the superconducting transition temperature. For this work, we develop the theory of quantum fluctuations within such multifractally enhanced superconducting states in thin films. In conditions of weak disorder, we employ the Finkel’stein nonlinear sigma model to derive an effective action for the superconducting order parameter and the quasiclassical Green’s function, meticulously accounting for the influence of quantum fluctuations. This effective action, applicable for interactions of any strength, reveals the critical role of well-known collective modes in a dirty superconductor, and its saddle-point analysis leads to modified Usadel and gap equations. These equations comprehensively incorporate the renormalizations stemming from the interplay between interactions and disorder, resulting in the non-trivial energy dependence of the gap function. Notably, our analysis establishes a direct relation between the self-consistent gap equation at the superconducting transition temperature and the known renormalization group equations for interaction parameters in the normal state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2022

Control of second-harmonic generation in all-dielectric intersubband metasurfaces by controlling the polarity of χ (2)

All-dielectric metasurfaces have recently led to a paradigm shift in nonlinear optics as they allow for circumventing the phase matching constraints of bulk crystals and offer high nonlinear conversion efficiencies when normalized by the light-matter interaction volume. Unlike bulk crystals, in all-dielectric metasurfaces nonlinear conversion efficiencies primarily rely on the material nonlinearity, field enhancements, and the modal overlaps, therefore most efforts to date have only focused on utilizing these degrees of freedom. In this work, we demonstrate that for second-harmonic generation in all-dielectric metasurfaces, an additional degree of freedom is the control of the polarity of the nonlinear susceptibility. We demonstrate that semiconductor heterostructures that support resonant nonlinearities based on quantum-engineered intersubband transitions provide this new degree of freedom. We can flip and control the polarity of the nonlinear susceptibility of the dielectric medium along the growth direction and couple it to the Mie-type photonic modes. Here we demonstrate that engineering the χ (2) polarity in the meta-atom enables the control of the second-harmonic radiation pattern and conversion efficiency. Our results therefore open a new direction for engineering and optimizing second-harmonic generation using all-dielectric intersubband nonlinear metasurfaces.

42 ENGINEERING↗
DOE OSTI2024

Nonperturbative Zou-Wang-Mandel effect

The Zou-Wang-Mandel (ZWM) effect is a remarkable consequence of photon indistinguishability and continuous-variable entanglement in which an optical phase shift is imprinted on photonic modes associated with optical paths that that do not pass through the phase shift source. By bringing the canonical formalism of continuous-variable Gaussian states to bear on the mode-structure of the ZWM experiment, we show that the physical consequence of implementing optical path identity is a renormalization of quadrature squeezing which governs the entanglement of four effective optical modes. Nonperturbative expressions for the ZWM interference patterns and normalized first-order coherence function are derived. Generalizations to $\mathscr{H}$-graph states with more than four modes directly follow from the general method used to analyze the minimal example. Here, we show that a ZWM interferometer with a laser-seeded signal mode, which estimates an idler phase shift by detecting photons that did not propagate through the phase shift, exhibits an optimal sensitivity comparable to that of a laser-seeded SU(1,1) interferometer if path identity is implemented with high fidelity.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Suhl instabilities in nanoscopic spheroids

We simulate the magnetization dynamics of a permalloy spheroid of nanoscopic size in zero external field, such that both dipolar and exchange interactions are important. Low excitation power is used to obtain the frequencies and mode patterns of many normal modes. At higher power, non-linear three and four mode couplings between magnons carrying orbital angular momentum are observed to give rise to Suhl instabilities. Suhl’s analysis of the selection rules governing the allowed processes is extended to initial states other than uniform precession. These rules are studied and confirmed by the simulations. Both down- and up-conversion are seen as well as three and four-mode processes. As a result, general trends are inferred for preferred instabilities among those that are allowed, although the thresholds for some instabilities appear to be very high.

3D nanomagnetism↗
DOE OSTI2021

Machine learning methods for probabilistic locked-mode predictors in tokamak plasmas

A rotating tokamak plasma can interact resonantly with the external helical magnetic perturbations, also known as error fields. This can lead to locking and then to disruptions. We leverage machine learning (ML) methods to predict the locking events. We use a coupled third-order nonlinear ordinary differential equation model to represent the interaction of the magnetic perturbation and the plasma rotation with the error field. This model is sufficient to describe qualitatively the locking and unlocking bifurcations. Here, we explore using ML algorithms with the simulation data and experimental data, focusing on the methods that can be used with sparse datasets. These methods lead to the possibility of the avoidance of locking in real-time operations. We describe the operational space in terms of two control parameters: the magnitude of the error field and the rotation frequency associated with the momentum source that maintains the plasma rotation. The outcomes are quan- tified by order parameters that completely characterize the state, whether locked or unlocked. We use unsupervised ML methods to classify locked/unlocked states and note the usefulness of a certain normalization of the order parameters. Three supervised ML classifiers are used in suite to estimate the probability of locking in the region of control parameter space with hysteresis, i.e., the set of control parameters for which both locked and unlocked states can exist. The results show that a neural network gives the best estimate of the locking probability. An analogy of the present locking model with the van der Waals equation of state is also provided.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Mean-field equation for phase-modulated optical parametric oscillator

The widely established techniques for the generation of ultrashort optical pulses rely on passive mode locking of lasers, with the output pulse duration and emission spectrum determined by the intrinsic lifetime of laser transition in the gain medium. Due to the instantaneous nature of nonlinear gain, optical parametric oscillators (OPOs) are capable of generating optical radiation in all timescales from continuous-wave (cw) to ultrashort femtosecond regime, if driven by laser pump sources in the corresponding time domain. In the ultrashort timescale, operation of OPOs conventionally relies on mode-locked pump lasers, with the concomitant disadvantages of large footprint and high cost. At the same time, the lack of gain storage mandates the use of synchronous pumping, resulting in increased complexity. In this paper, we present the concept of phase-modulated OPO driven by cw pump laser. The approach overcomes the traditional drawbacks of ultrafast OPOs, enabling femtosecond pulse generation without the need for synchronous pumping, resulting in a simplified, compact, and cost-effective architecture using cw input pump lasers. We derive a mean-field equation for a degenerate χ ( 2 ) OPO driven by a cw laser with intracavity electro-optic modulator (EOM), and also including dispersion compensation. The equation predicts the formation of stable femtosecond pulses ( < 200 fs ) , in both normal and anomalous dispersion regimes, with a controllable repetition rate determined by the frequency of the EOM. The remarkable functionality of the proposed scheme paves the way for the development of an alternative class of widely tunable coherent femtosecond light sources in both bulk and integrated format based on χ ( 2 ) OPOs using cw pump lasers. Published by the American Physical Society 2024

Sanchez, A. D. (ORCID:0000000327630339)↗
DOE OSTI2023

DFT-Based Calculation of Molecular Hyperpolarizability and SFG Intensity of Symmetric and Asymmetric Stretch Modes of Alkyl Groups

Vibrational sum frequency generation (SFG) spectroscopy has been extensively used for obtaining structural information of molecular functional groups at two-dimensional (2D) interfaces buried in the gas or liquid medium. Although the SFG experiment can be done elegantly, interpreting the measured intensity in terms of molecular orientation with respect to the lab coordinate is quite complicated. One of the main reasons is the difficulty of determining the hyperpolarizability tensors of even simple molecules that govern their SFG responses. The single-bond polarizability derivative (SBPD) model has been proposed to estimate the relative magnitude of SFG-active hyperpolarizability by assuming that the perturbation associated to each vibration is negligible. In this study, density functional theory (DFT) was used to calculate the polarizability and dipole derivative tensors of the CH 3 stretch mode of CH 3 I, CH 3 CH 2 I, CH 3 OH, CH 3 CH 2 OH. Then, the hyperpolarizability tensors of symmetric and asymmetric vibration modes were calculated considering the Boltzmann distribution of representative conformers, which allowed theoretical calculation of their SFG intensities at all polarization combinations as a function of the tilt angle of the CH 3 group with respect to the surface normal direction. Then the ratios of the calculated SFG intensities for the CH 3 peaks used in experimental studies for the CH 3 tilt angle determination were compared. This comparison clearly showed the effect of vibrational coupling among neighboring functional groups. As a result, this study presents new parameters that can be used in determining the average tilt angle of the CH 3 group at the 2D interface with SFG measurements as well as limitations of the method.

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

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗
DOE OSTI2021

Kinetic-ballooning-mode turbulence in low-average-magnetic-shear equilibria

In this work, kinetic-ballooning-mode (KBM) turbulence is studied via gyrokinetic flux-tube simulations in three magnetic equilibria that exhibit small average magnetic shear: the Helically Symmetric eXperiment (HSX), the helical-axis Heliotron-J and a circular tokamak geometry. For HSX, the onset of KBM being the dominant instability at low wavenumber occurs at a critical value of normalized plasma pressure $\beta ^{\rm KBM}_{\rm crit}$ that is an order of magnitude smaller than the magnetohydrodynamic (MHD) ballooning limit $\beta ^{\rm MHD}_{\rm crit}$ when a strong ion temperature gradient (ITG) is present. However, $\beta ^{\rm KBM}_{\rm crit}$ increases and approaches the MHD ballooning limit as the ITG tends to zero. For these configurations, $\beta ^{\rm KBM}_{\rm crit}$ also increases as the magnitude of the average magnetic shear increases, regardless of the sign of the normalized magnetic shear. Simulations of Heliotron-J and a circular axisymmetric geometry display behaviour similar to HSX with respect to $\beta ^{\rm KBM}_{\rm crit}$ . Despite large KBM growth rates at long wavelengths in HSX, saturation of KBM turbulence with $\beta > \beta _{\rm crit}^{\rm KBM}$ is achievable in HSX and results in lower heat transport relative to the electrostatic limit by a factor of roughly five. Nonlinear simulations also show that KBM transport dominates the dynamics when KBMs are destabilized linearly, even if KBM growth rates are subdominant to ITG growth rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗