Parametric wave amplification and mixing in the ionosphere
Ionospheric parametric amplification and frequency mixing due to polarized TEM waves nonlinear interaction with longitudinal electrostatic waves
SEARCH · Engineering Papers
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.
Ionospheric parametric amplification and frequency mixing due to polarized TEM waves nonlinear interaction with longitudinal electrostatic waves
Harmonic radiation from electrical power transmission lines in the range of a few kHz leaks into the magnetosphere and stimulates a coherent wave instability, resulting in strong amplification of the input waves and the generation of free-running emissions. A description is given of some recent observational results that provide new information on the power line radiation (PLR) phenomenon. It is pointed out that PLR stimulates many subtle and complex wave-particle interactions in the magnetosphere that are similar to those simulated by controlled transmitter signals. These interactions undoubtedly affect both wave and particle environments in the magnetosphere. However, a quantitative assessment of their importance is not possible until further information becomes available.
The structure of an artificially generated sinusoidal water wave train of fixed frequency under the influence of wind is analyzed. Artificially generated waves of 1 Hz are studied at seven wind speeds in the range of 140-400 cm/s. It is observed that the water wave train deformed by wind consists of two components at both the fundamental mode and the harmonics. The amplitude and phase of the wave components are derived, and the dispersion relation and component phase speeds are examined. The data reveal that the amplitude of the forced and free-traveling second harmonics correlate with previous theories, and that the deviation of the measured phase speed from the linear theory is caused by the nonlinearity of the primary wave, the interaction between short gravity waves and the primary wave, and the advection effects of wind drift.
A two-dimensional computational code, PRLUS2D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent-boundary-layer interactions. The problem of compression corners at supersonic speeds was solved using the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were used. Computational results showed that the RPLUS2D code compared very well with experimentally obtained data for supersonic compression corner flows, except in the case of large separated flows resulting from the interactions between the shock wave and turbulent boundary layer. The computational results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien two-equation turbulence model performing better than the Baldwin-Lomax model.
Electromagnetic ion cyclotron (EMIC) waves are an important magnetospheric emission, which is excited near the magnetic equator with frequencies below the proton gyro-frequency. The source of bee energy for wave growth is provided by temperature anisotropy of ring current (RC) ions, which develops naturally during inward convection from the plasma sheet These waves strongly affect the dynamic s of resonant RC ions, thermal electrons and ions, and the outer radiation belt relativistic electrons, leading to non-adiabatic particle heating and/or pitch-angle scattering and loss to the atmosphere. The rate of ion and electron scattering/heating is strongly controlled by the Wave power spectral and spatial distributions, but unfortunately, the currently available observational information regarding EMIC wave power spectral density is poor. So combinations of reliable data and theoretical models should be utilized in order to obtain the power spectral density of EMIC waves over the entire magnetosphere throughout the different storm phases. In this study, we present the simulation results, which are based on two coupled RC models that our group has developed. The first model deals with the large-scale magnetosphere-ionosphere electrodynamic coupling, and provides a self-consistent description of RC ions/electrons and the magnetospheric electric field. The second model is based on a coupled system of two kinetic equations, one equation describes the RC ion dynamics and another equation describes the power spectral density evolution of EMIC waves, and self-consistently treats a micro-scale electrodynamic coupling of RC and EMIC waves. So far, these two models have been applied independently. However, the large-scale magnetosphere-ionosphere electrodynamics controls the convective patterns of both the RC ions and plasmasphere altering conditions for EMIC wave-particle interaction. In turn, the wave induced RC precipitation Changes the local field-aligned current distributions and the ionospheric conductances, which are crucial for a large-scale electrodynamics. The initial results from this new self-consistent model of the magnetospheric electric field, RC and EMIC waves will be shown in this presentation.
Using a dynamical-system approach, we have investigated the efficiency of large-amplitude whistler waves for causing microburst precipitation in planetary radiation belts by modeling the microburst energy and particle fluxes produced as a result of nonlinear wave-particle interactions. We show that wave parameters, consistent with large amplitude oblique whistlers, can commonly generate microbursts of electrons with hundreds of keV-energies as a result of Landau trapping. Relativistic microbursts (greater than 1 MeV) can also be generated by a similar mechanism, but require waves with large propagation angles Theta (sub k)B greater than 50 degrees and phase-speeds v(sub phi) greater than or equal to c/9. Using our result for precipitating density and energy fluxes, we argue that holes in the distribution function of electrons near the magnetic mirror point can result in the generation of double layers and electron solitary holes consistent in scales (of the order of Debye lengths) to nonlinear structures observed in the radiation belts by the Van Allen Probes. Our results indicate a relationship between nonlinear electrostatic and electromagnetic structures in the dynamics of planetary radiation belts and their role in the cyclical production of energetic electrons (E greater than or equal to 100 keV) on kinetic timescales, which is much faster than previously inferred.
Abstract In the diffuse aurora, magnetospheric electrons, initially precipitated from the inner plasma sheet via wave-particle interaction processes, degrade in the atmosphere toward lower energies, and produce secondary electrons via impact ionization of the neutral atmosphere. These initially precipitating electrons of magnetospheric origin can also be additionally reflected back into the magnetosphere, leading to a series of multiple reactions by the two magnetically conjugate atmospheres that can greatly impact the initially precipitating flux at the upper ionospheric boundary (700-800 km). The resultant population of secondary and primary electrons cascades toward lower energies and escape back to the magnetosphere. Escaping upward electrons traveling from the ionosphere can be trapped in the magnetosphere, as they travel inside the loss cone, via Coulomb collisions with the cold plasma, or by interactions with various plasma waves. Even though this scenario is intuitively transparent, this magnetosphere-ionosphere coupling element is not considered in any of the existing diffuse aurora research. Nevertheless, as we demonstrate in this letter, this process has the potential to dramatically affect the formation of electron precipitated fluxes in the regions of diffuse auroras.
During the last six months, we have continued our study of the excitation of electromagnetic waves in the whistler frequency range and the role that these waves will play in the acceleration of electrons and ions in the auroral region. A paper entitled 'Electron Beam Excitation of Upstream Waves in the Whistler Mode Frequency Range' was listed in the Journal of Geophysical Research. In this paper, we have shown that an anisotropic electron beam (or gyrating electron beam) is capable of generating both left-hand and right-hand polarized electromagnetic waves in the whistler frequency range. Since right-hand polarized electromagnetic waves can interact with background electrons and left-hand polarized waves can interact with background ions through cyclotron resonance, it is possible that these beam generated left-hand and right-hand polarized electromagnetic waves can accelerate either ions or electrons (or both), depending on the physical parameters under consideration. We are currently carrying out a comprehensive study of the electromagnetic whistler and lower hybrid like waves observed in the auroral zone using both wave and particle data. Our first task is to identify these wave modes and compare it with particle observations. Using both the DE-1 particle and wave measurements, we can positively identify those electromagnetics lower hybrid like waves as fast magnetosonic waves and the upper cutoff of these waves is the local lower hybrid frequency. From the upper cutoff of the frequency spectrum, one can infer the particle density and the result is in very good agreement with the particle data. Since these electromagnetic lower hybrid like waves can have frequencies extended down to the local ion cyclotron frequency, it practically confirms that they are not whistler waves.
Attention is directed to the acoustics research of the 1950s and 1960s for guidance in understanding and quantizing the turbulence amplification that can occur in regions of shock-wave boundary-layer interaction. Three primary turbulence amplifier-generator mechanisms are identified and shown, by linear analysis, to be responsible for turbulence amplification across a shock wave in excess of 100% of the incident turbulence intensity.
Interactions of long-wave solitons propagating in shear flows are described by a coupled pair of Korteweg-de Vries equations. The basic equation of motion for the analysis is the quasi-geostrophic forecast equation, and the interaction of two wave modes is studied. The solution for mode 1-mode 2 interaction of solitary waves in an asymmetric shear flow of a barotropic atmosphere with divergence is constructed. Streamline patterns for certain flows are obtained. An unsteady solitary wave solution for a modified Korteweg-de Vries equation is derived.
This paper describes a combined experimental and theoretical study of three-dimensional swept shock wave-turbulent boundary layer interactions at Mach 4 generated by a sharp fin of angles alpha equals 16 and 20 degrees. The theoretical model is the three-dimensional compressible Reynolds-averaged Navier-Stokes equations with turbulence incorporated through the algebraic eddy viscosity model of Baldwin and Lomax. Previous computations have been performed by Horstman using the Baldwin-Lomax, Cebeci-Smith and Jones Launder models. Computed results for the surface pressure, skin friction and streamline angles are compared with experiment and previous numerical results. The present results display good agreement with experimental data for surface pressure and surface flow direction. All turbulence models fail to accurately predict the peak skin friction. The computed flowfields are in agreement with many of the features of the quasi-conical flowfield model of Settles.
The problem of a normal shock wave impinging on a flat plate, turbulent boundary layer is considered for the case where the external flow is transonic. Asymptotic methods are employed. It is shown that there are two outer regions, including the outer part of the boundary layer and the external flow, in which inviscid flow governing equations hold, and two regions near the wall, in which Reynolds and/or viscous stresses need be taken into account. The solutions in the outer regions lead to the pressure distribution on the wall, for which an analytical expression is presented, valid under those conditions when separation is imminent but has not yet occurred. The solutions valid in the inner regions lead to a separation criterion.
The nature of the 3-s ultralow frequency (ULF) wave in the Earth's foreshock region and the associated wave-particle interaction are not yet well understood. In this work, we investigate the 3-s ULF waves using Magnetospheric Multiscale (MMS) observations. By combining the plasma rest frame wave properties obtained from multiple methods with the instability analysis based on the velocity distribution in the linear wave stage, the ULF wave is determined to be due to the ion/ion nonresonant mode instability. The interaction between the wave and ions is analyzed using the phase relationship between the transverse wave fields and ion velocities and using the longitudinal momentum equation. During the stage when ULF waves have sinusoidal waveforms up to |dB|/|B 0 | ~ 3, where dB is the wave magnetic field and B0 is the background magnetic field, the wave electric fields perpendicular to B 0 do negative work to solar wind ions; along B 0 , a longitudinal electric field develops, but the V × B force is stronger and leads to solar wind ion deceleration. During the same wave stage, the backstreaming beam ions gain energy from the transverse wave fields and get deceleration along B 0 by the longitudinal electric field. The ULF wave leads to electron heating, preferentially in the direction perpendicular to the local magnetic field. Secondary waves are generated within the ULF waveforms, including whistler waves near half of the electron cyclotron frequency, high-frequency electrostatic waves, and magnetosonic whistler waves. The work improves the understanding of the nature of 3-s ULF waves and the associated wave-particle interaction.
Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.
The transition Reynolds number for shear layers produced by interactions between weak and strong shock waves is determined on the basis of experiments performed in a 20-in. (Mach 6) and an 11-in. (Mach 6.9) hypersonic tunnel. A variable angle wedge was used to generate a planar shock wave which interacted with the bow wave of a blunt body. An average value of the transition length (defined as the length along the shear layer from the shock interaction to the point where turbulence became visible on schlieren photographs) was used to determine the transition Reynolds number.
Fluctuating surface pressure measurements have been made to investigate the effectiveness of boundary layer separators (BLS's) in reducing the fluctuating pressure loads produced by separated shock wave turbulent boundary layer interactions. Measurements have been made under unswept and swept compression corner interactions in a Mach 5 flow. BLS's fix the separation location and eliminate the large-amplitude, low-frequency fluctuating pressure loads upstream of the compression corners. The loads on the unswept compression corner face are reduced by as much as 59%. The BLS's also shift the mean pressure distribution on the unswept corner face in the streamwise direction. Results show that the loads on the corner face vary with the BLS height and the distance between the BLS and the compression corner. Suggestions for the optimum placement and the use of the BLS's are also made.
A combined experimental and numerical study has been conducted on the hypersonic shock-wave turbulent-boundary layer interaction at Mach 8.2 generated by a single fin of angles alpha = 10 and 15 deg. Three models are considered: (1) the 3D compressible Reynolds-averaged Navier-Stokes (RANS) equations using the k-epsilon turbulence model, (2) the 3D RANS using the Rodi turbulence model, and (3) the conical RANS using the Baldwin-Lomax algebraic turbulence model. The computations are compared with various experimental data. The computations using models (1) and (2) show quantitatively very similar results and very good agreement with experimental data for surface pressure and skin friction. Comparison with boundary layer profiles of pitot pressure and yaw angle are also generally good, but the peak surface heat transfer is overestimated by up to 48 percent. The effect of the laminar boundary layer on the fin is restricted to the immediate vicinity of the fin surface. Conical calculations using model (3) show substantially poorer agreement with experiment.
Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually involve turbo-ramjet configurations. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. In the first phase of our study, an over/under nozzle configuration was analyzed. The two plumes from the turbojet and ramjet interact at the end of a common 2-D cowl, where they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. For the problem analyzed, the turbojet engine operates at a higher nozzle pressure ratio than the ramjet, causes the turbojet plume overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data were used to compare with a computational fluid dynamics (CFD) study using the PARC2D code. The CFD results were in good agreement with both static pressure distributions on the cowl separator and on nozzle walls. The thrust coefficients were also in reasonable agreement. In addition, inviscid relationships were developed around the confluence point, where the two exhaust jets meet, and these results compared favorably with the CFD results. In the second phase of our study, a 3-D CFD solution was generated to compare with the 2-D solution. The major difference between the 2-D and 3-D solutions was the interaction of the shock waves, generated by the plume interactions, on the sidewall. When a shock wave interacts with a sidewall and sidewall boundary layer, it is called a glancing shock sidewall interaction. These interactions entrain boundary layer flow down the shockline into a vortical flow pattern. The 3-D plots show the streamlines being entrained down the shockline. The pressure of the flow also decreases slightly as the sidewall is approached. Other difference between the 2-D and 3-D solutions were a lowering of the nozzle thrust coefficient value from 0.9850 (2-D) to 0.9807 (3-D), where the experimental value was 0.9790. In the third phase of our study, a different turbo-ramjet configuration was analyzed. The confluence of a supersonic turbojet and a subsonic ramjet in the turbine based combined-cycle (TBCC) propulsion system was studied by a 2-D CFD code. In the analysis, Mach 1.4 primary turbojet was mixed with the subsonic ramjet secondary flow in an ejector mode operation. Reasonable agreements were obtained with the supplied I-D TBCC solutions. For low downstream backpressure, the Fabri choke condition (Break-Point condition) was observed in the secondary flow within mixing zone. For sufficient high downstream backpressure, the Fabri choke no longer exist, the ramjet flow was reduced and the ejector flow became backpressure dependent. Highly non-uniform flow at ejector exit were observed, indicated that for smooth downstream combustion, the mixing of the two streams probably required some physical devices.