Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear wave propagation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗
DOE OSTI2023

Characteristics and selection of near-fault simulated earthquake ground motions for nonlinear analysis of buildings

Earthquake-induced ground shaking near rupturing faults is highly sensitive to the rupture characteristics, seismic wave propagation patterns and site conditions, and field recordings of near-fault shaking are relatively sparse. These challenges complicate the assessment of the seismic performance of near-fault structures. A common approach to representing near-fault ground motion in engineering analysis is to explicitly consider and select records with strong directivity pulses (pulse records). We use three-dimensional high-resolution physics-based earthquake simulations to test this approach in the context of scenario-based ground motion record selection, and to study the important characteristics of near-fault ground shaking. We highlight the deficiencies associated with classifying near-fault simulated records as “pulse” or “non-pulse,” based on the presence of a single dominating pulse in the velocity time history. We show that this approach is inadequate for characterizing near-fault shaking on soft soils which can be dominated by both forward rupture directivity and basin amplification effects. We conduct ground motion selection experiments for the analysis of near-fault structures with and without explicit classification of the pulse features in the records, and evaluate the bias in the predicted structural demands. We find that the maximum interstory drift demands on building structures imposed by unscaled site-specific simulated ground motion records selected based on relevant spectral shape features are not sensitive to the classification of records as pulse/non-pulse. Therefore, with regard to predicting the maximum interstory drifts in near-fault buildings, we do not find justification for the binary pulse classification of near-fault records.

42 ENGINEERING↗
NASA NTRS1985

Local effects of gravity wave propagation and saturation

In recent years, gravity waves were recognized to play a major role in the dynamics of the middle atmosphere. Perhaps the major effect of such motions are the reversal of the vertical shear of the mean zonal wind and the occurrence of a large turbulent diffusivity in the mesosphere due to gravity wave saturation. Yet, despite the importance of these gravity wave effects, the processes and the consequences of gravity wave propagation and saturation are only beginning to be understood in detail. The linear saturation theory predicts drag and turbulent diffusion due to saturating wave motion. This theory, however, fails to address a number of issues that are certain to be important for gravity wave propagation and saturation in the middle atmosphere. These issues, including wave transience, wave superposition, local convective adjustment, and nonlinearity, are discussed.

Fritts, D. C.↗
NASA NTRS1986

The wave magnetic field amplitude threshold for nonlinear trapping of energetic gyroresonant and Landau resonant electrons by nonducted VLF waves in the magnetosphere

A recently developed theory is used to calculate the wave magnetic field amplitude threshold, B-tau, necessary to allow the nonlinear trapping of energetic gyroresonant and Landau resonant electrons by VLF whistler mode waves in the magnetosphere propagating at an arbitrary angle, psi, with respect to the earth's magnetic field. A detailed raytracing study is carried out to establish the variation of psi with position along each magnetic shell in the range L = 2-5 and for frequencies 5-17.8 kHz. It is found that the minimum values of B-tau along each L shell generally occur at points of second-order resonance. In general, for nonducted fixed frequency signals there is only one point of second-order resonance on each L shell, and this is located within a few degrees of the magnetic equator. However, over a narrow range of L, there are as many as three points of second-order resonance and as many as three associated minima in B-tau. At least one of these points is located more than 10 deg from the magnetic equator.

Bell, T. F.↗
NASA NTRS2017

Structure and Properties of the Foreshock at Venus

The interaction of the solar wind with Venus is dominated by the planet's ionosphere that acts as an obstacle to the flow resulting in an induced magnetosphere and bow shock much smaller than their terrestrial counterparts. This study presents a 3-D electromagnetic hybrid (kinetic ions and fluid electrons) simulation of the solar wind interaction with an unmagnetized obstacle to examine the structure and properties of the Cytherean foreshock during periods of near radial IMF, that is, when it lies upstream of the ionosphere. The interaction between the backstreaming ions and the solar wind results in the generation of two classes of ULF waves: (1) parallel propagating sinusoidal waves with periods approx. 20-30 s and (2) highly oblique fast magnetosonic waves. The joint nonlinear evolution of these waves results in the formation of structures called foreshock cavitons with dimensions comparable to the size of the planet. Foreshock cavitons are also present in the terrestrial foreshock. The excavation of plasma and magnetic field from their cores leads to lower average densities and magnetic field strengths in the foreshock. As in the case of Earth, this excavation results in the formation of a fast magnetosonic pulse/shock at the edge of the foreshock named the foreshock compressional boundary. Also similar to Earth is the formation of spontaneous hot flow anomalies (SHFAs) as foreshock cavitons approach the bow shock. The size and properties of SHFAs at Venus are comparable to those at Earth, and their existence has recently been established at Mars and Venus in a companion paper.

Omidi, N.↗
DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗
NASA NTRS1974

Large-amplitude hydromagnetic waves

Several aspects of the theory of large-amplitude hydromagnetic waves and their behavior in the interplanetary medium are examined. The characteristic modes of the full (i.e., nonlinearized) MHD equations and their modification by collisionless and finite-frequency effects are considered. Special attention is paid to the transverse Alfven mode, which is undamped and characterized by strictly constant pressure, density, and B; this seems to be the predominant propagating fluctuation at 1 AU. It is shown that its propagation in the small-wavelength (WKB) approximation is essentially identical to that of the small-amplitude Alfven wave of linearized theory. It is also suggested that its presence at 1 AU may provide a natural explanation of the observed power anisotropy of the fluctuations. A second-order analysis is used to study fluctuations that are not characteristic modes. It is found that for a small range of propagation directions, and subject to third-order effects, a finite-amplitude wave can exist that is linearly polarized with delta B perpendicular to both B sub zero and k; such a wave can damp nonlinearly.

Barnes, A.↗
NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗
NASA NTRS1992

Parametric instabilities of large amplitude Alfven waves with obliquely propagating sidebands

This paper presents a brief report on properties of the parametric decay and modulational, filamentation, and magnetoacoustic instabilities of a large amplitude, circularly polarized Alfven wave. We allow the daughter and sideband waves to propagate at an arbitrary angle to the background magnetic field so that the electrostatic and electromagnetic characteristics of these waves are coupled. We investigate the dependance of these instabilities on dispersion, plasma/beta, pump wave amplitude, and propagation angle. Analytical and numerical results are compared with numerical simulations to investigate the full nonlinear evolution of these instabilities.

Vinas, A. F.↗
NASA NTRS1995

High amplitude waves in the expanding solar wind plasma

We simulated the 1-D nonlinear time-evolution of high-amplitude Alfven, slow and fast magnetoacoustic waves in the solar wind propagating outward at different angles to the mean magnetic (spiral) field, using the expanding box model. The simulation results for Alfven waves and fast magnetoacustic waves fit the observational constraints in the solar wind best, showing decreasing trends for energies and other rms-quantities due to expansion and the appearance of inward propagating waves as minor species in the wind. Inward propagating waves are generated by reflection of Alfven waves propagating at large angles to the magnetic field or they coincide with the occurrence of compressible fluctuations. In our simulations, fast and slow magnetoacoustic waves seem to have a level in the density-fluctuations which is too high when we compare with the observations. Furthermore, the evolution of energies for slow magnetoacoustic waves differs strongly from the evolution of fluctuation energies in situ.

Schmidt, J. M.↗
NASA NTRS1989

Oxygen acoustic solitary waves in a magnetized plasma

Ion-acoustic solitary waves in a magnetized plasma containing an arbitrary mixture of H(+) and O(+) ions are studied. A nonlinear wave equation has been derived from the Poisson-Vlasov equations, including a uniform magnetic field and dissipation due to reflected electrons. When dissipation is ignored, the equation has soliton solutions associated with both oxygen and hydrogen acoustic modes, which can be either rarefactive or compressive depending on the ion concentrations and the electron/ion temperature ratio and, more weakly, on the bulk drifts of the species. If electron reflection is included, the solitary wave can be intensified. Under somewhat restrictive conditions the oxygen solitary wave is rarefactive and propagates with a velocity comparable to that observed by the Viking satellite. The three-dimensional solitons obey a relation of scales parallel to the magnetic field and in the transverse direction. Computer simulations of one-dimensional versions of the nonlinear wave equation are presented.

Qian, S.↗
DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS1992

Observations of short large-amplitude magnetic structures at a quasi-parallel shock

Results of a detailed analysis of short large-amplitude magnetic structures (SLAMS) observed at an encounter of the quasi-parallel blow shock by the AMPTE UKS and IRM satellites are presented. Isolated SLAMS, surrounded by solar wind conditions, and embedded SLAMS, which lie within or form the boundary with regions of significant heating and deceleration, are identified. The duration, polarization, and other characteristics of SLAMS are all consistent with their growth directly out of the ULF wave field, including the common occurrence of an attached whistler as found in ULF shocklets. The plasma rest frame propagation speeds and two-spacecraft time delays for all cases show that the SLAMS attempt to propagate upstream against the oncoming flow, but are convected back downstream. The speeds and delays vary systematically with SLAMS amplitude in the way anticipated from nonlinear wave theory, as do their polarization features. Inter-SLAMS regions and boundary regions with solar wind contain hot deflected ions of lesser density than within the SLAMS.

Schwartz, Steven J.↗
NASA NTRS1972

Convective amplification of Type 1 irregularities in the equatorial electrojet

Wave propagation and refraction of Type 1 irregularities in the equatorial electrojet were investigated. Quantitative calculation of wave refraction in a model electrojet showed that the direction of wave refraction must change sign at one altitude. Waves propagating with the electrons rotate their wave vectors upwards in the upper electrojet and downwards in the lower electrojet during the day, and vice versa at night. Furthermore, the altitude region of largest linear growth rate is also the one with the weakest refraction rate. Consequently, computations of the ray-path integrated wave growth shows that this region would dominate the backscatter spectrum from the electrojet if linear theory were valid, and it is further noted that the maximum amplitude wave should have phase velocities exceeding the ion acoustic speed. It was concluded that propagation alone, without inclusion of nonlinear effects, cannot explain backscatter observations of a constant Doppler frequency shift given by the ion acoustic speed.

Lee, K.↗
NASA NTRS1995

Application of essentially nonoscillatory methods to aeroacoustic flow problems

A finite-difference essentially nonoscillatory (ENO) method has been applied to several of the problems prescribed for the workshop sponsored jointly by the Institute for Computer Applications in Science and Engineering and by NASA Langley Research Center entitled 'Benchmark Problems in Computational Aeroacoustics'. The workshop focused on computational challenges specific to aeroacoustics. Among these are long-distance propagation of a short-wavelength disturbance, propagation of small-amplitude disturbances, and nonreflective boundary conditions. The shock capturing-capability inherent to the ENO method effectively eliminates oscillations near shock waves without the need to add and tune dissipation or filter terms. The method-of-lines approach allows the temporal and spatial operators to be chosen separately in accordance with the demands of a particular problem. The ENO method was robust and accurate for all problems in which the propagating wave was resolved with 8 or more points per wavelength. The finite-wave-model boundary condition, a local nonlinear acoustic boundary condition, performed well for the one-dimensional problems. The buffer-domain approach performed well for the two-dimensional test problem. The amplitude of nonphysical reflections were less than 1 percent of the exiting wave's amplitude.

Atkins, Harold L.↗
NASA NTRS2010

Heat Entrapment Effects Within Liquid Acquisition Devices

We introduce a model problem to address heat entrapment effects or the local accumulation of thermal energy within liquid acquisition devices. We show that the parametric space consists of six parameters, namely the Rayleigh and Prandtl numbers, the aspect ratio, and heat flux ratios for the bottom, side, and top boundaries of the enclosure. For the range of Ra considered 1 to 10(sup 9), beyond Ra on the order of 10(sup 5), convective instability is the dominant mode of convection in comparison to natural convection. The flow field transitions to asymmetric modes at Ra on the order of 10(sup 7). Direct numerical simulation of a large geometric length scale prototype for Ra on the order of 10(sup 9) shows that the flow field evolves from small wavelength instability which gives rise to nonlinear growth of thermals, propagation of the instability occurs via growth of secondary and tertiary modes, and a travelling wave mode occurs prior to asymmetry. The effect of a large aspect ratio is to increase the number of modes in the vertical direction. Due to the slow diffusion of heat in the prototype, asymptotic states are not readily attained, we show that dynamical similarity can be used for a model which allows the attainment of asymptotic states and that transition to a chaotic state occurs for Ra on the order of 10(sup 9) via a broadband power spectrum. These dynamical events show that for the baseline condition in which heat is absorbed from background laboratory environment, higher heat flux is absorbed at the top and bottom boundaries of the enclosure than a nominal value of 34.9 ergs per square centimeter -second.

Duval, W. M. B.↗
NASA NTRS2002

Nonlinear Internal Waves in the South China Sea During ASIAEX

Internal wave distribution maps have been compiled from more than one hundred ERS-1/2, RADARSAT and Space Shuttle SAR images in the South China Sea (SCS) from 1993 to 2000. Based on these distribution maps, most of internal waves in the northeast part of SCS were propagating westward. The wave crest can be as long as 200 km with amplitude of 100 m, due to strong current from the Kuroshio branching out into the SCS. In recent Asian Seas International Acoustics Experiment (ASIAEX), moorings have been deployed in April 2000 and May 2001. Simultaneous RADARSAT ScanSAR images have been collected during the field test to integrate with the model and the in-situ measurements in the SCS. During ASIAEX in May 2001, many large internal waves were observed at the test area and were the major features for acoustic volume interaction. The environmental parameters have been calculated based on extensive CTD castings and mooring data. Nonlinear internal wave models have been applied to simulate the wave evolution on the continental shelf and the results compare reasonably with mooring measurements. The evolution and dissipation of huge internal waves on the shelf break, mode-two waves, elevation waves, and wave-wave interaction are very important issues for acoustic propagation. The implication of internal wave effects on acoustic propagation will also be discussed.

ACOUSTIC PROPAGATION↗