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.

At least 163 records · Page 9

Dynamic Instability of Solidification Fronts and Formation of Transverse Layered Structures

A plane front of crystallization in a unidirectional process of solidification of a binary melt or a solution is shown to be unstable with respect to virtual perturbations of the front velocity within some region of physical and operating parameters. This dynamic instability has nothing in common with morphological instability of the front also possible under certain conditions. It gives rise to the origination of a self-oscillating regime of solidification. The oscillations of the front velocity are accompanied by oscillations of other relevant variables, such as the melt concentration and temperature at the front. What is most important is that this provides for a natural explanation of phenomena of the formation of ordered striations commonly observed when producing alloys and withdrawing crystals. A theory of nonlinear instability of the solidification front is developed which also allows for properties of both self-oscillating regimes of solidification and transverse layered structures it allows to be described in an explicit form, without introducing any additional hypotheses and empirical parameters. Apart from other things, the properties are strongly dependent on convective transport originated in consequence of a difference between the densities of the liquid melt and solid alloy. The theory results in a remarkably good agreement with experimental evidence on crystal withdrawal obtained in normal gravity and under space flight conditions. The problem of interplay between modes of the dynamic and morphological types of instability is also discussed. There are parametric regions in which either only one of these types occurs or both of them arise simultaneously.

Buyevich, Yu A.↗

Co‐Evolution of Specific Stiffness and Permeability of Rock Fractures Offset in Shear

Abstract Fractures and faults represent planes of weakness and compliance in rock masses that serve as focal points for both microearthquakes and fluid transport, with seismicity and permeability evolution closely linked. Contact stiffness is highly stress‐sensitive and directly influences permeability. We explore the co‐evolution of specific stiffness and permeability of rough fractures under normal stress and shear offset using numerical simulations. Individual rough fractures are represented by variable amplitude (Root mean square) and wavelength ( λ ) using a granular mechanics model. Contacting rough surfaces are mated, offset in shear, and then compacted in displacement mode. The compacting fractures generate stress‐dependent changes in contact porosity, which govern both permeability and stiffness evolution. We establish a universal dimensionless relationship linking specific stiffness and permeability that inherently incorporates the effects of surface roughness, shear offset, and microcracking. The observed cracking effect—where local stress redistribution and pressure‐driven microcrack propagation dynamically alter the aperture field—introduces a nonlinear permeability response at high stress. Increased roughness amplitude and larger shear offsets reduce stiffness while dampening permeability sensitivity to stress, demonstrating a strong interplay between surface texture and hydro‐mechanical behavior. While the model captures this behavior effectively, deviations emerge at very low porosities due to extreme aperture sensitivity in this limit.

He, Xinxin [Energy and Mineral Engineering Pennsyl↗

Nonlinear dynamic analysis of quasi-symmetric anisotropic structures

An efficient computational method for the nonlinear dynamic analysis of quasi-symmetric anisotropic structures is proposed. The application of mixed models simplifies the analytical development and improves the accuracy of the response predictions, and operator splitting allows the reduction of the analysis model of the quasi-symmetric structure to that of the corresponding symmetric structure. The preconditoned conjugate gradient provides a stable and effective technique for generating the unsymmetric response of the structure as the sum of a symmetrized response plus correction modes. The effectiveness of the strategy is demonstrated with the example of a laminated anisotropic shallow shell of quadrilateral planform subjected to uniform normal loading.

Noor, Ahmed K.↗

Influence of Stationary Crossflow Modulation on Secondary Instability

A likely scenario for swept wing transition on subsonic aircraft with natural laminar flow involves the breakdown of stationary crossflow vortices via high frequency secondary instability. A majority of the prior research on this secondary instability has focused on crossflow vortices with a single dominant spanwise wavelength. This paper investigates the effects of the spanwise modulation of stationary crossflow vortices at a specified wavelength by a subharmonic stationary mode. Secondary instability of the modulated crossflow pattern is studied using planar, partial-differential-equation based eigenvalue analysis. Computations reveal that weak modulation by the first subharmonic of the input stationary mode leads to mode splitting that is particularly obvious for Y-type secondary modes that are driven by the wall-normal shear of the basic state. Thus, for each Y mode corresponding to the fundamental wavelength of results in unmodulated train of crossflow vortices, the modulated flow supports a pair of secondary modes with somewhat different amplification rates. The mode splitting phenomenon suggests that a more complex stationary modulation such as that induced by natural surface roughness would yield a considerably richer spectrum of secondary instability modes. Even modest levels of subharmonic modulation are shown to have a strong effect on the overall amplification of secondary disturbances, particularly the Z-modes driven by the spanwise shear of the basic state. Preliminary computations related to the nonlinear breakdown of these secondary disturbances provide interesting insights into the process of crossflow transition in the presence of the first subharmonic of the dominant stationary vortex.

Vortex structure↗

Direct Multipoint Observations Capturing the Reformation of a Supercritical Fast Magnetosonic Shock

Using multipoint Magnetospheric Multiscale (MMS) observations in an unusual string-of-pearls configuration, we examine in detail observations of the reformation of a fast magnetosonic shock observed on the upstream edge of a foreshock transient structure upstream of Earth’s bow shock. The four MMS spacecraft were separated by several hundred km, comparable to suprathermal ion gyro-radius scales or several ion inertial lengths. At least half of the shock reformation cycle was observed, with a new shock ramp rising up out of the “foot” region of the original shock ramp. Using the multipoint observations, we convert the observed time-series data into distance along the shock normal in the shock’s rest frame. That conversion allows for a unique study of the relative spatial scales of the shock’s various features, including the shock’s growth rate, and how they evolve during the reformation cycle. Analysis indicates that: the growth rate increases during reformation, electron-scale physics play an important role in the shock reformation, and energy conversion processes also undergo the same cyclical periodicity as reformation. Strong, thin electron-kinetic-scale current sheets and large-amplitude electrostatic and electromagnetic waves are reported. Evidence is also presented of nonlinear wave decay from electromagnetic whistler-mode “lion roars” to electrostatic solitary waves in the downstream plasma regime. Results highlight the critical cross-scale coupling between electron-kinetic-and ion-kinetic-scale processes and details of the nature of nonstationarity, shock-front reformation at collisionless, fast magnetosonic shocks.

D L Turner↗

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↗

Current Sheets in the Heliosheath: Voyager 1, 2009

We identified all of the current sheets for which we have relatively complete and accurate magnetic field (B) data from Voyager 1 (V1) from days of year (DOYs) 1 to 331, 2009, which were obtained deep in the heliosheath between 108.5 and 111.8 AU. Three types of current sheets were found: (1) 15 proton boundary layers (PBLs), (2) 10 and 3 magnetic holes and magnetic humps, respectively, and (3) 3 sector boundaries. The magnetic field strength changes across PBL, and the profile B(t) is linearly related to the hyperbolic tangent function, but the direction of B does not change. For each of the three sector boundaries, B rotated in a plane normal to the minimum variance direction, and the component of B along the minimum variance direction was zero within the uncertainties, indicating that the sector boundaries were tangential discontinuities. The structure of the sector boundaries was not as simple as that for PBLs. The average thickness of magnetic holes and humps (approx.30 RL) was twice that of the PBLs (approx.15 RL). The average thickness of the current sheets associated with sector boundaries was close to the thickness of the PBLs. Our observations are consistent with the hypothesis that magnetic holes and humps are solitons, which are initiated by the mirror mode instability, and evolve by nonlinear kinetic plasma processes to pressure balanced structures maintained by magnetization currents and proton drift currents in the gradients of B.

Burlaga, L. F.↗

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↗

Electrostatic Kelvin-Helmholtz instability in a radially injected plasma cloud

An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.

Sydora, R. D.↗

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↗

Drift wave instability in the Io plasma torus

A linear normal mode analysis of the drift wave instability in the Io plasma torus was carried out on the basis of the Richmond (1973) and Huang et al. (1990) analyses of drift waves in the vicinity of the earth's plasmapause. Results indicate that the outer torus boundary is linearly unstable to the growth of electrostatic drift waves. It is shown that the linear growth rate is proportional to the ion drift frequency and to the ratio of the flux tube charge content to the Jovian ionospheric Pedersen conductance. It is also shown that various theoretical models of global radial transport in Jupiter's atmosphere (including corotating convection, interchange diffusion, and transient flux tube convection) can be understood as plausible nonlinear evolutions of electrostatic drift waves.

Huang, T. S.↗

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↗

Multiple switch-off slow shock solutions

A new electromagnetic quasi-neutral time-dependent implicit hybrid model for simulating the slow-mode switch-off shock is developed, in which the electrons are treated as a finite-mass adiabatic fluid (gamma sub e = 5/3), whose motion along the shock normal provides local charge neutrality. The model is used to examine the effects of various boundary conditions applied during shock formation. It is shown that, for each upstream condition, two distinctly different slow-shock structures can be obtained: a low-shock structure with a trailing magnetic wave (TMW) and a slow-shock structure without a TMW, depending strongly on the boundary conditions for the electric and magnetic fields at the physical piston which generates the shock. The results indicate that the slow-shock structure with a TMW is nonlinearly unstable with respect to changes in the upstream magnetic field, suggesting that the downstream structure of slow shocks need not be associated with a TMW.

Vu, H. X.↗

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↗

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↗

Coupling of radiation into thin film modes by means of localized plasma resonances

The interaction between the surface plasmon mode that propagates at a metal dielectric interface and the localized plasma resonances (LPR) is investigated experimentally in Ag-island films. A stair-stepped sample geometry comprising a glass substrate, a continuous 50-nm Ag film, an LiF spacer film of thickness d = 5-60 nm, and an Ag-island film of mass thickness 3 nm is used in near-normal-reflectivity and plasmon-propagation-constant (k) determinations. The results are presented graphically and discussed. The overall shape of the reflectivity curves is found to be characteristic of Ag films, but with a dip at about 400 nm (corresponding to the absorption resonance of the island film) which is most pronounced with d = 25 nm. It is inferred that the island resonances are strongly coupled to a continuous-film dissipative mechanism at this d value. This inference is supported by the fact that the variation in k, correctd for LiF effects and plotted as a function of d, is greatest at around d = 25 nm. The implications of this finding for broad-band coupling into a thin-film mode, LPR enhancement of waveguide nonlinear effects, and new surface-enhanced-Raman-scattering geometries are indicated.

Holland, W. R.↗

A Thermodynamic Theory Of Solid Viscoelasticity. Part 1: Linear Viscoelasticity.

The present series of three consecutive papers develops a general theory for linear and finite solid viscoelasticity. Because the most important object for nonlinear studies are rubber-like materials, the general approach is specified in a form convenient for solving problems important for many industries that involve rubber-like materials. General linear and nonlinear theories for non-isothermal deformations of viscoelastic solids are developed based on the quasi-linear approach of non-equilibrium thermodynamics. In this, the first paper of the series, we analyze non-isothermal linear viscoelasticity, which is applicable in a range of small strains not only to all synthetic polymers and bio-polymers but also to some non-polymeric materials. Although the linear case seems to be well developed, there still are some reasons to implement a thermodynamic derivation of constitutive equations for solid-like, non-isothermal, linear viscoelasticity. The most important is the thermodynamic modeling of thermo-rheological complexity , i.e. different temperature dependences of relaxation parameters in various parts of relaxation spectrum. A special structure of interaction matrices is established for different physical mechanisms contributed to the normal relaxation modes. This structure seems to be in accord with observations, and creates a simple mathematical framework for both continuum and molecular theories of the thermo-rheological complex relaxation phenomena. Finally, a unified approach is briefly discussed that, in principle, allows combining both the long time (discrete) and short time (continuous) descriptions of relaxation behaviors for polymers in the rubbery and glassy regions.

Freed, Alan D.↗

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↗