The Darwin model as a tool for electromagnetic plasma simulation
Self consistent Darwin model for slow wave electromagnetic plasma simulation
Kaufman, A. N.↗
SEARCH · Engineering Papers
Results for “simulated 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.
Study of Ardmore, Oklahoma storm clouds. I - Convective storm cloud initiation and development based on the remote sensing gravity-wave-induced convection. II - Satellite infrared remote sensing and numerical simulation
The role of gravity waves is discussed with respect to the vertical velocity of convection. Specific attention is given to wave-induced convection which contributes to the fractions of formation and the development of severe convective storms. Large-amplitude gravity waves and convective instability were investigated in storm clouds above Ardmore, Oklahoma. Rapid-scan satellite imagery and radar summaries provide evidence of water-vapor condensation related to convection which is introduced by gravity waves. Gravity wave periods of 35 minutes are found to initiate weak convection, which can be intensified by gravity waves with periods of 20 minutes. The convective motion reaches a maximum about one hour before funnel clouds develop. Other mechanisms which contribute to convective motion are considered, but gravity waves are the major contributor to the initiation, formation, and development of mesoscale storm clouds. Cloud modeling based on satellite imagery and sounding data showed that by overshooting cloud tops that penetrated the tropopause, storm clouds mature; that these clouds collapse about 9 minutes before the touchdown of tornadoes; and that cloud tops collapse at a high rate about 6 minutes before tornadoes lift off.
Hung, R. J.↗
An Incidence Loss Model for Wave Rotors with Axially Aligned Passages
A simple mathematical model is described to account for the losses incurred when the flow in the duct (port) of a wave rotor is not aligned with the passages. The model, specifically for wave rotors with axially aligned passages, describes a loss mechanism which is sensitive to incident flow angle and Mach number. Implementation of the model in a one-dimensional CFD based wave rotor simulation is presented. Comparisons with limited experimental results are consistent with the model. Sensitivity studies are presented which highlight the significance of the incidence loss relative to other loss mechanisms in the wave rotor.
Paxson, Daniel E.↗
Global simulation of auroral arcs
Numerical simulation of global formation of auroral arcs is carried out for a three-dimensional coupled ionosphere-magnetosphere system. With the presence of the stationary convection electric field, the field-aligned current driven by the polarization and drift of the ionospheric conductivity enhancement is shown to be responsible for the feedback instability of the shear Alfven wave. The simulation result demonstrates the dynamic process of auroral arc formation, where an ionospheric perturbation nearly resonant with the magnetospheric Alfven wave (toroidal mode) can eventually grow into auroral arcs. Many of the important characteristics of quiet auroral arcs are demonstrated, and it is shown that the overall distribution of auroras depends critically on the stationary state of the auroral oval characterized by the large-scale Birkeland current, the electric field, and the electron density distribution.
Miura, A.↗
Computation of acoustic waves in a jet
A numerical treatment of acoustic waves in a jet is described. The full time dependent Euler equations are used in both linear and nonlinear formulations. The computational region of integration is artificially bounded and boundary conditions are developed to simulate outgoing waves and to enable the computational domain to be substantially restricted. Higher order methods and coordinate transformations are introduced to further reduce the number of grid points as well as to increase the efficiency of the program. Numerical results are presented for time harmonic sources as well as for sources with more complicated time dependence.
Bayliss, A.↗
Long-Time Numerical Integration of the Three-Dimensional Wave Equation in the Vicinity of a Moving Source
We propose a family of algorithms for solving numerically a Cauchy problem for the three-dimensional wave equation. The sources that drive the equation (i.e., the right-hand side) are compactly supported in space for any given time; they, however, may actually move in space with a subsonic speed. The solution is calculated inside a finite domain (e.g., sphere) that also moves with a subsonic speed and always contains the support of the right-hand side. The algorithms employ a standard consistent and stable explicit finite-difference scheme for the wave equation. They allow one to calculate tile solution for arbitrarily long time intervals without error accumulation and with the fixed non-growing amount of tile CPU time and memory required for advancing one time step. The algorithms are inherently three-dimensional; they rely on the presence of lacunae in the solutions of the wave equation in oddly dimensional spaces. The methodology presented in the paper is, in fact, a building block for constructing the nonlocal highly accurate unsteady artificial boundary conditions to be used for the numerical simulation of waves propagating with finite speed over unbounded domains.
Ryabenkii, V. S.↗
Assessment of the Atmospheric Channel for Short (Ka-Band and Optical) Wavelengths
Atmospheric turbulence under clear sky conditions is an impairment of the atmospheric channel that greatly affects propagation of optical signal in the troposphere. The turbulence manifests itself in a number of forms within the optical domain, from the twinkling of a star in a clear night, to resolution degradation in a large aperture telescope. Therefore, a body of analytical, numerical, and experimental tools has been developed in optics to study, simulate, and control effects of atmospheric turbulence on an optical signal. Incidentally, there has been an increasing demand for high data rate returns from NASA missions which has led to envision utilizing a carrier signal in the Ka-Band range. The impact of atmospheric turbulence effects must be evaluated and considered for this frequency domain. The purpose of this work is to show that when the turbulence strength from the optical case to the KaBand ease is properly scaled, one can apply the same mathematical simulation developed for optical to predict turbulence effects within the Ka-Band domain. As a demonstration of this principle, we present how the scintillations of a Ka-Band downlink return of a deep space signal was successfully reproduced through wave-optics simulation.
optical communications↗
Density response of the mesospheric sodium layer to gravity wave perturbations
Lidar observations of the mesospheric sodium layer often reveal wavelike features moving through the layer. It is often assumed that these features are a layer density response to gravity waves. Chiu and Ching (1978) described the approximate form of the linear response of atmospheric layers to gravity waves. In this paper, their results are used to predict the response of the sodium layer to gravity waves. These simulations are compared with experimental observations and a good correlation is found between the two. Because of the thickness of the sodium layer and the density gradients found in it, a linear model of the layer response is not always adequate to describe gravity wave-sodium layer interactions. Inclusion of nonlinearities in the layer response is briefly discussed. Experimental data is seen to contain features consistent with the predicted nonlinearities.
Shelton, J. D.↗
Asymptotic analysis of dissipative waves with applications to their numerical simulation
Various problems involving the interplay of asymptotics and numerics in the analysis of wave propagation in dissipative systems are studied. A general approach to the asymptotic analysis of linear, dissipative waves is developed. It was applied to the derivation of asymptotic boundary conditions for numerical solutions on unbounded domains. Applications include the Navier-Stokes equations. Multidimensional traveling wave solutions to reaction-diffusion equations are also considered. A preliminary numerical investigation of a thermo-diffusive model of flame propagation in a channel with heat loss at the walls is presented.
Hagstrom, Thomas↗
Impact of Ring Current Ions on Electromagnetic Ion Cyclotron Wave Dispersion Relation
Effect of the ring current ions in the real part of electromagnetic ion Cyclotron wave dispersion relation is studied on global scale. Recent Cluster observations by Engebretson et al. showed that although the temperature anisotropy of is energetic (> 10 keV) ring current protons was high during the entire 22 November 2003 perigee pass, electromagnetic ion cyclotron waves were observed only in conjunction with intensification of the ion fluxes below 1 keV by over an order of magnitude. To study the effect of the ring current ions on the wave dispersive properties and the corresponding global wave redistribution, we use a self-consistent model of interacting ring current and electromagnetic ion cyclotron waves, and simulate the May 1998 storm. The main findings of our simulation can be summarized as follows: First, the plasma density enhancement in the night MLT sector during the main and recovery storm phases is mostly caused by injection of suprathermal plasma sheet H + (approximately < 1 keV), which dominate the thermal plasma density. Second, during the recovery storm phases, the ring current modification of the wave dispersion relation leads to a qualitative change of the wave patterns in the postmidnight-dawn sector for L > 4.75. This "new" wave activity is well organized by outward edges of dense suprathermal ring current spots, and the waves are not observed if the ring current ions are not included in the real part of dispersion relation. Third, the most intense wave-induced ring current precipitation is located in the night MLT sector and caused by modification of the wave dispersion relation. The strongest precipitating fluxes of about 8 X 10(exp 6)/ (cm(exp 2) - s X st) are found near L=5.75, MLT=2 during the early recovery phase on 4 May. Finally, the nightside precipitation is more intense than the dayside fluxes, even if there are less intense waves, because the convection field moves ring current ions into the loss cone on the nightside, but drives them out of the loss cone on the dayside. So convection and wave scattering reinforce each other in the nightside, but interfere in the dayside sector.
Khazanov, G. V.↗
Finite-difference theory for sound propagation in a lined duct with uniform flow using the wave envelope concept
Finite difference equations are derived for sound propagation in a two dimensional, straight, soft wall duct with a uniform flow by using the wave envelope concept. This concept reduces the required number of finite difference grid points by one to two orders of magnitude depending on the length of the duct and the frequency of the sound. The governing acoustic difference equations in complex notation are derived. An exit condition is developed that allows a duct of finite length to simulate the wave propagation in an infinitely long duct. Sample calculations presented for a plane wave incident upon the acoustic liner show the numerical theory to be in good agreement with closed form analytical theory. Complete pressure and velocity printouts are given to some sample problems and can be used to debug and check future computer programs.
Baumeister, K. J.↗
Ionospheric Response to Electron Precipitation Driven by Whistler-Mode Chorus Waves Measured by Arase Satellite and Simulated by STET and SPEAH-RIT
The origin of pulsating auroras has long been an open question. Arase satellite observation reported the direct evidence that precipitated electrons in pulsating aurora are driven by whistler-mode chorus wave activity. Our study presents the theoretical analysis of this observational event based on SuperThermal Electron Transport (STET) code that simulates the highly dynamic environment of measured wave intensities and particle fluxes. Specifically, the STET code simulated results confirm the delicate loss-cone observation results of this mission and reveal the broader energy range of precipitated electron fluxes that were not measurable by Arase satellite. These highly dynamic precipitating electron fluxes were applied to a newly developed code, the Superthermal Proton, Electron and Atomic Hydrogen tRansport in the Ionosphere and Thermosphere (SPEAH-RIT), to simulate ionospheric response of this phenomena, in particular, the enhancement of electron density, plasma temperatures, and ionospheric electric conductance at the footprints of Arase satellite.
George V. Khazanov↗
Application of PML Absorbing Boundary Conditions to the Benchmark Problems of Computational Aeroacoustics
Accurate numerical non-reflecting boundary conditions are important in all the proposed benchmark problems of the Second Workshop. Recently, a new absorbing boundary condition has been developed using Perfectly Matched Layer (PML) equations for the Euler equations. In this approach, a region with a width of a few grid points is introduced adjacent to the non-reflecting boundaries. In the added region, Perfectly Matched Layer equations are constructed and applied so that the out-going waves are absorbed inside the layer with little reflection to the interior domain. It will be demonstrated in the present paper that the proposed absorbing boundary condition is quite general and versatile, applicable to radiation boundaries as well as inflow and outflow boundaries. It is also easy to implement. The emphasis of the paper will be on the application of the PML absorbing boundary condition to problems in Categories 1, 2, and 3. In Category 1, solutions of problems 1 and 2 are presented. Both problems are solved using a multi-domain polar grid system. Perfectly Matched Layer equations for a circular boundary are constructed and their effectiveness assessed. In Category 2, solutions of problem 2 are presented. Here, in addition to the radiation boundary conditions at the far field in the axisymmetric coordinate system, the inflow boundary condition at the duct inlet is also dealt with using the proposed Perfectly Match Layer equations. At the inlet, a PML domain is introduced in which the incident duct mode is simulated while the waves reflected from the open end of the duct are absorbed at the same time. In Category 3, solutions of all three problems are presented. Again, the PML absorbing boundary condition is used at the inflow boundary so that the incoming vorticity wave is simulated while the outgoing acoustic waves are absorbed with very little numerical reflection. All the problems are solved using central difference schemes for spatial discretizations and the optimized Low-Dissipation and Low-Dispersion Runge-Kutta scheme for the time integration. Issues of numerical accuracy and efficiency are also addressed.
Hu, Fang Q.↗
Constraints on cosmic ray propagation in the galaxy
The goal was to derive a more detailed picture of magnetohydrodynamic turbulence in the interstellar medium and its effects on cosmic ray propagation. To do so, radio astronomical observations (scattering and Faraday rotation) were combined with knowledge of solar system spacecraft observations of MHD turbulence, simulations of wave propagation, and modeling of the galactic distribution to improve the knowledge. A more sophisticated model was developed for the galactic distribution of electron density turbulence. Faraday rotation measure data was analyzed to constrain magnetic field fluctuations in the ISM. VLBI observations were acquired of compact sources behind the supernova remnant CTA1. Simple calculations were made about the energies of the turbulence assuming a direct link between electron density and magnetic field variations. A simulation is outlined of cosmic ray propagation through the galaxy using the above results.
Cordes, James M.↗
Evaluation of Intranasal Scopolamine for the Prevention of Wave Motion Induced Motion Sickness While Maintaining Performance on Operationally Relevant Tasks
Motion sickness represents one of the greatest clinical challenges impacting crew activities during and following g-transitions. Our overall goal is to characterize the effectiveness of motion sickness countermeasures during controlled laboratory experiments using capsule wave motion simulation and during field testing in operational environments. This laboratory study focused on prevention of motion sickness using intranasal scopolamine using a double-blinded repeated measures design in 30 subjects (19M, 11F). Intranasal scopolamine was provided by Defender Pharmaceuticals, Inc. (DPI-386 Nasal Gel, referred to as Inscop) self-administered by a nasal pump (Aptar Pharma) that delivers 0.4 mg dose (0.2 mg / nostril). During each session, subjects were exposed to complex wave motion on a six degree-of-freedom platform that included pitch, roll and heave at provocative stimulus frequencies (0.1-0.25 Hz) while seated in an illuminated cabin deprived of external visual cues. Motion sickness symptoms were compared across treatment and placebo control sessions counterbalanced across subjects and separated by at least one week. The time-to-motion sickness endpoint was based on severe malaise, defined as symptom score of ≥8 points using the Pensacola Diagnostic Index [1]. The bioavailability of scopolamine for each session was estimated from plasma concentrations obtained every 15 min[2]. Side effects during the treatment session were minimal, and performance was not impaired on a test battery including motion perception tracking, tablet-based eye-hand coordination and psychomotor vigilance testing. The plasma concentration remained near-peak levels throughout the 45 min motion sickness testing for most subjects. The percentile ranking on the Motion Sickness Susceptibility Questionnaire [MSSQ, 3]was moderately correlated with the motion sickness time-to-endpoint during the placebo control session (rho = –0.3, p=0.056). Seventeen subjects did not reach an endpoint during their placebo session and were eliminated from subsequent analysis. Another subject was excluded due to insufficient plasma concentration during the treatment session. For the remaining 12 subjects, the change in time-to-motion sickness endpoint between placebo and treatment sessions was moderately correlated with plasma concentration (rho =0.48, p = 0.056), improving on average 4.4 ± 18.4 min, mean ± std with Inscop versus placebo. Our results are consistent with previous findings that intranasal delivery of scopolamine can be effective at reducing motion sickness symptoms with minimal cognitive or sedative side effects. Future work is needed to optimize the delivery of Inscop for rescue (treatment) of symptoms following g-transitions as one of the key advantages of this formulation is self-administration in a suited environment. We will also explore how the combination of Inscop with non-pharmaceutical sensory aids, e.g., vibrotactile feedback of Earth vertical, may further mitigate motion sickness and improve task performance.
S J Wood↗
Simulation of the nonlinear evolution of electron plasma waves
Electrostatic waves driven by an electron beam in an ambient magnetized plasma were studied using a quasi-1D PIC simulation of electron plasma waves (i.e., Langmuir waves). The results disclose the presence of a process for moving wave energy from frequencies and wavenumbers predicted by linear theory to the Langmuir-like frequencies during saturation of the instability. A decay process for producing backward propagating Langmuir-like waves, along with low-frequency waves, is observed. The simulation results, however, indicate that the backscattering process is not the conventional Langmuir wave decay. Electrostatic waves near multiples of the electron plasma frequency are generated by wave-wave coupling during the nonlinear stage of the simulations, confirming the suggestion of Klimas (1983).
Nishikawa, K.-I.↗
Numerical solutions for the complete shock wave structure behind supersonic edge delta wings
Computerized simulation of shock waves and flow fields behind supersonic edge delta wings and wing-body combinations
Kutler, P.↗
Long-lived spiral waves in N-body simulations
Results of N-body simulations of disc galaxies using a two-dimensional Cartesian N-body code are presented. Both trailing arm spirals (TAS) and leading arm spirals (LAS) were used with varieties of pitch angles and pattern speeds. LAS perturbations transferred their energy to TAS via swing amplification; TAS perturbations led to TAS arms. In both cases the spiral arms persisted for more than 5 rotation periods, but the maximum amplitude 2-armed spirals were generated by LAS perturbations. The persistence of the trailing arm spiral waves is thought to be caused by the kinematic spiral arm mechanism described by Kalnajs (1973).
Comins, Neil F.↗