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 901 records · Page 50

Nonlinear Elastic Effects on the Energy Flux Deviation of Ultrasonic Waves in GR/EP Composites

In isotropic materials, the direction of the energy flux (energy per unit time per unit area) of an ultrasonic plane wave is always along the same direction as the normal to the wave front. In anisotropic materials, however, this is true only along symmetry directions. Along other directions, the energy flux of the wave deviates from the intended direction of propagation. This phenomenon is known as energy flux deviation and is illustrated. The direction of the energy flux is dependent on the elastic coefficients of the material. This effect has been demonstrated in many anisotropic crystalline materials. In transparent quartz crystals, Schlieren photographs have been obtained which allow visualization of the ultrasonic waves and the energy flux deviation. The energy flux deviation in graphite/epoxy (gr/ep) composite materials can be quite large because of their high anisotropy. The flux deviation angle has been calculated for unidirectional gr/ep composites as a function of both fiber orientation and fiber volume content. Experimental measurements have also been made in unidirectional composites. It has been further demonstrated that changes in composite materials which alter the elastic properties such as moisture absorption by the matrix or fiber degradation, can be detected nondestructively by measurements of the energy flux shift. In this research, the effects of nonlinear elasticity on energy flux deviation in unidirectional gr/ep composites were studied. Because of elastic nonlinearity, the angle of the energy flux deviation was shown to be a function of applied stress. This shift in flux deviation was modeled using acoustoelastic theory and the previously measured second and third order elastic stiffness coefficients for T300/5208 gr/ep. Two conditions of applied uniaxial stress were considered. In the first case, the direction of applied uniaxial stress was along the fiber axis (x3) while in the second case it was perpendicular to the fiber axis along the laminate stacking direction (x1).

Prosser, William H.↗

Numerical modeling of the atmosphere with an isentropic vertical coordinate

A theta-coordinate model simulating the nonlinear evolution of a baroclinic wave is presented. In the model, vertical discretization maintains important integral constraints such as conservation of the angular momentum and total energy. A massless-layer approach is used in the treatment of the intersections of coordinate surfaces with the lower boundary. This formally eliminates the intersection problem, but raises other computational problems. Horizontal discretization of the continuity and momentum equations in the model are designed to overcome these problems. Selected results from a 10-day integration with the 25-layer, beta-plane version of the model are presented. It is concluded that the model can simulate the nonlinear evolution of a baroclinic wave and associated dynamical processes without major computational difficulties.

Hsu, Yueh-Jiuan G.↗

Two-layer viscous instability in a rotating couette device

A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Gallagher, Christopher↗

Dependence of the kinetic energy absorption capacity of bistable mechanical metamaterials on impactor mass and velocity

Using an alternative mechanism to dissipation or scattering, bistable structures and mechanical metamaterials have shown promise for mitigating the detrimental effects of impact by reversibly locking energy into strained material. Herein, we extend prior works on impact absorption via bistable metamaterials to computationally explore the dependence of kinetic energy transmission on the velocity and mass of the impactor, with strain rates exceeding 10 2 s –1 . We observe a large dependence on both impactor parameters, ranging from significantly better to worse performance than a comparative linear material. We then correlate the variability in performance to solitary wave formation in the system and give analytical estimates of idealized energy absorption capacity under dynamic loading. In addition, we find a significant dependence on damping accompanied by a qualitative difference in solitary wave propagation within the system. The complex dynamics revealed in this study offer potential future guidance for the application of bistable metamaterials to applications including human and engineered system shock and impact protection devices.

42 ENGINEERING↗

Self-compression of laser pulses in a discrete medium

Using the discrete nonlinear Schrödinger equation, we study the laser radiation self-action dynamics in a nanostructured waveguide system. It is shown that for a laser pulse energy exceeding the critical value, the nonlinear evolution of the wave field markedly differs from the corresponding process in a continuous medium. While the pulse propagates in a discrete medium, its length decreases to values comparable to the scale of the structure and the velocity decreases to zero. The system parameters are determined at which the wave field after the laser pulse termination begins to propagate in the opposite direction in a compressed form. The process of slowing down and stopping the pulse is accompanied by strong radiation losses, so that about a third of the energy of the initial field distribution remains in the final compressed state. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Quantum computing for fusion energy science applications

This is a review of recent research exploring and extending present-day quantum computing capabilities for fusion energy science applications. We begin with a brief tutorial on both ideal and open quantum dynamics, universal quantum computation, and quantum algorithms. Then, we explore the topic of using quantum computers to simulate both linear and nonlinear dynamics in greater detail. Because quantum computers can only efficiently perform linear operations on the quantum state, it is challenging to perform nonlinear operations that are generically required to describe the nonlinear differential equations of interest. In this work, we extend previous results on embedding nonlinear systems within linear systems by explicitly deriving the connection between the Koopman evolution operator, the Perron–Frobenius evolution operator, and the Koopman–von Neumann evolution (KvN) operator. We also explicitly derive the connection between the Koopman and Carleman approaches to embedding. Extension of the KvN framework to the complex-analytic setting relevant to Carleman embedding, and the proof that different choices of complex analytic reproducing kernel Hilbert spaces depend on the choice of Hilbert space metric are covered in the appendixes. Finally, we conclude with a review of recent quantum hardware implementations of algorithms on present-day quantum hardware platforms that may one day be accelerated through Hamiltonian simulation. We discuss the simulation of toy models of wave–particle interactions through the simulation of quantum maps and of wave–wave interactions important in nonlinear plasma dynamics.

Joseph, I. (ORCID:0000000255400840)↗

Oblique instability of quasi-parallel whistler waves in the presence of cold and warm electron populations

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

79 ASTRONOMY AND ASTROPHYSICS↗

Auroral arc formation - Kinetic and MHD effects

Recent observations by the Dynamics Explorer satellites indicate that auroral field lines are not equipotentials. The field-aligned 'conductivity' is the same for up and downgoing currents and of the order of several 10 to the -8th to several 10 to the -9th/ohms sq m. This is not in agreement with models based on the resistance provided by the magnetic mirror force alone. It is suggested that the resistivity is due to damping of kinetic shear-Alfven wave packets. This damping may be nonlinear in nature or Landau damping. The properties of kinetic Alfven waves are reviewed.

Goertz, C. K.↗

Boundary Layer Instabilities Generated by Freestream Laser Perturbations

A controlled, laser-generated, freestream perturbation was created in the freestream of the Boeing/AFOSR Mach-6 Quiet Tunnel (BAM6QT). The freestream perturbation convected downstream in the Mach-6 wind tunnel to interact with a flared cone model. The geometry of the flared cone is a body of revolution bounded by a circular arc with a 3-meter radius. Fourteen PCB 132A31 pressure transducers were used to measure a wave packet generated in the cone boundary layer by the freestream perturbation. This wave packet grew large and became nonlinear before experiencing natural transition in quiet flow. Breakdown of this wave packet occurred when the amplitude of the pressure fluctuations was approximately 10% of the surface pressure for a nominally sharp nosetip. The initial amplitude of the second mode instability on the blunt flared cone is estimated to be on the order of 10 −6 times the freestream static pressure. The freestream laser-generated perturbation was positioned upstream of the model in three different configurations: on the centerline, offset from the centerline by 1.5 mm, and offset from the centerline by 3.0 mm. When the perturbation was offset from the centerline of a blunt flared cone, a larger wave packet was generated on the side toward which the perturbation was offset. The offset perturbation did not show as much of an effect on the wave packet on a sharp flared cone as it did on a blunt flared cone.

Chou, Amanda↗

Oblique Instability of Quasi-Parallel Whistler Waves in the Presence of Cold and Warm Electron Populations

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

Vadim Roytershteyn↗

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗

Understanding Dynamics and Thermodynamics of ENSO and Its Complexity Simulated by E3SM and Other Climate Models

Despite the seeming success of most state-of-the-art climate models in simulating the El Niño-Southern Oscillation (ENSO), there is strong evidence that models achieve realistic levels of ENSO activity do so often for wrong reasons. This is owing to an often occurred near cancelation of large errors in terms of contributions to ENSO growth rate from coupled dynamic and thermodynamic feedback processes. Climate models remain deficient in simulating the observed ENSO’s spatial and temporal complexity that involves interplays of coupled dynamic and thermodynamic feedbacks, interactions across multiple scales, nonlinear processes in the tropical atmosphere and ocean system, biases in mean sate and physical processes, and influences external to equatorial Pacific coupled ENSO dynamics. Our proposed research aims at advancing predictive and process-level understandings of ENSO simulated in E3SM and other climate models under current and future climate conditions with two main objectives: (i) better understanding the aforementioned broad range interactive processes and sources that control fundamental properties of ENSO in E3SM and CIMP6 outputs as well as in observational (reanalysis) data sets, using a hierarchical of coupled dynamical frameworks consisting of theoretical analysis, intermediate complexity modeling, and coupled dynamic diagnostics; (ii) to use this understanding to explore pathways towards improving E3SM’s capability of simulating ENSO and its complexity. More specifically, we will focus on four main thrusts of research: (1) ENSO’s dynamic and thermodynamic feedbacks; (2) the across-scale interactions of ENSO with annual cycle and MJO/WWB/TIW (Madden Julian Oscillation/Westerly Wind Burst/Tropical Instability Wave) activity; (3) key nonlinear processes of ENSO involving atmospheric convective thresholds, nonlinear ocean dynamic heating, and thermocline outcropping; and (4) the impacts of climate mean-state biases/changes and perturbed physical processes on simulated ENSO and its complexity.

54 ENVIRONMENTAL SCIENCES↗

Ion-acoustic solitary waves in a magnetized plasma with arbitrary electron equation of state

The oblique propagation of fully nonlinear, slow ion-acoustic solitary waves in a collisionless, low-beta, magnetized plasma is examined. The analysis includes the effects of a finite ion pressure, electron trapping, and multicomponent particle populations. The existence of both compressive and rarefactive modes propagating obliquely to the magnetic field in a plasma with two distinct Boltzmann electron populations and cold ions is demonstrated. It is shown that paired electrostatic shocks observed over the earth's auroral zone may be closely related to the rarefactive modes. As a measure of the collisionless dissipation encountered by the solitary waves, the linear response of the plasma to slow ion-acoustic waves is also examined.

Witt, E.↗

Extreme terahertz magnon multiplication induced by resonant magnetic pulse pairs

Nonlinear interactions of spin-waves and their quanta, magnons, have emerged as prominent candidates for interference-based technology, ranging from quantum transduction to antiferromagnetic spintronics. Yet magnon multiplication in the terahertz (THz) spectral region represents a major challenge. Intense, resonant magnetic fields from THz pulse-pairs with controllable phases and amplitudes enable high order THz magnon multiplication, distinct from non-resonant nonlinearities such as the high harmonic generation by below-band gap electric fields. Here, we demonstrate exceptionally high-order THz nonlinear magnonics. It manifests as 7th-order spin-wave-mixing and 6th harmonic magnon generation in an antiferromagnetic orthoferrite. We use THz two-dimensional coherent spectroscopy to achieve high-sensitivity detection of nonlinear magnon interactions up to six-magnon quanta in strongly-driven many-magnon correlated states. The high-order magnon multiplication, supported by classical and quantum spin simulations, elucidates the significance of four-fold magnetic anisotropy and Dzyaloshinskii-Moriya symmetry breaking. Moreover, our results shed light on the potential quantum fluctuation properties inherent in nonlinear magnons.

36 MATERIALS SCIENCE↗

Computational procedures for mixed equations with shock waves

This paper discusses the procedures we have developed to treat a canonical problem involving a mixed nonlinear equation with boundary data that imply a discontinuous solution. This equation arises in various physical contexts and is basic to the description of the nonlinear acoustic behavior of a shock wave near a caustic. The numerical scheme developed is of second order, treats discontinuities as such by applying the appropriate jump conditions across them, and eliminates the numerical dissipation and dispersion associated with large gradients. Our results are compared with the results of a first-order scheme and with those of a second-order scheme we have developed. The algorithm used here can easily be generalized to more complicated problems, including transonic flows with imbedded shocks.

Yu, N. J.↗

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Generating optical supercontinuum and frequency comb in tenuous plasmas

There are several mechanisms by which the frequency spectrum of a laser broadens when it propagates at near-relativistic intensity in a tenuous plasma. Focusing on one-dimensional effects, we identify two strong optical nonlinearities, namely, four-wave mixing (FWM) and forward Raman scattering (FRS), for creating octave-wide spectra. FWM dominates the interaction when the laser pulse is short and intense, and its combination with phase modulation produces a symmetrically broadened supercontinuum. FRS dominates when the laser pulse is long and relatively weak, and it broadens the laser spectrum mainly toward lower frequencies and produces a frequency comb. The frequency chirping combined with group velocity dispersion compresses the laser pulse, causing higher peak intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗