Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffuse 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 1,063 records · Page 59

Simulating energetic ions and enhanced fusion rates from ion-cyclotron resonance heating with a full-wave/Fokker–Planck model

Reproducing fast-ion enhanced fusion rates from ion-cyclotron resonance heating (ICRH) in tokamaks requires the self-consistent coupling of a full-wave solver and a Fokker–Planck solver, which evolves multiple simultaneously resonant ion species. We introduce a new self-consistent model that iterates the TORIC full-wave solver with the CQL3D Fokker–Planck solver using the integrated plasma simulator (IPS). This model evolves the bounce-averaged ion distribution functions in both parallel and perpendicular velocity-space with a quasilinear radio frequency (RF) diffusion operator valid in the ion finite Larmor radius (FLR) limit and the RF electric fields with the resultant non-Maxwellian FLR dielectric tensor. This produces non-Maxwellian ICRH simulations that are fully self-consistent, fast, and interoperable with integrated modeling frameworks, such as TRANSP/GACODE/IPS-FASTRAN. We demonstrate our model's capabilities by validating it against experimental data in Alcator C-Mod. We then perform the first RF heating simulations of SPARC using self-consistent non-Maxwellian ion distributions to investigate the potential to enhance fusion rates using ion cyclotron resonance heating generated fast ions.

Physics↗

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗

Results of a zonally truncated three-dimensional model of the Venus middle atmosphere

Although the equatorial rotational speed of the solid surface of Venus is only 4 m s(exp-1), the atmospheric rotational speed reaches a maximum of approximately 100 m s(exp-1) near the equatorial cloud top level (65 to 70 km). This phenomenon, known as superrotation, is the central dynamical problem of the Venus atmosphere. We report here the results of numerical simulations aimed at clarifying the mechanism for maintaining the equatorial cloud top rotation. Maintenance of an equatorial rotational speed maximum above the surface requires waves or eddies that systematically transport angular momentum against its zonal mean gradient. The zonally symmetric Hadley circulation is driven thermally and acts to reduce the rotational speed at the equatorial cloud top level; thus wave or eddy transport must counter this tendency as well as friction. Planetary waves arising from horizontal shear instability of the zonal flow (barotropic instability) could maintain the equatorial rotation by transporting angular momentum horizontally from midlatitudes toward the equator. Alternatively, vertically propagating waves could provide the required momentum source. The relative motion between the rotating atmosphere and the pattern of solar heating, which as a maximum where solar radiation is absorbed near the cloud tops, drives diurnal and semidiurnal thermal tides that propagate vertically away from the cloud top level. The effect of this wave propagation is to transport momentum toward the cloud top level at low latitudes and accelerate the mean zonal flow there. We employ a semispectral primitive equation model with a zonal mean flow and zonal wavenumbers 1 and 2. These waves correspond to the diurnal and semidiurnal tides, but they can also be excited by barotropic or baroclinic instability. Waves of higher wavenumbers and interactions between the waves are neglected. Symmetry about the equator is assumed, so the model applies to one hemisphere and covers the altitude range 30 to 110 km. Horizontal resolution is 1.5 deg latitude, and vertical resolution is 1.5 km. Solar and thermal infrared heating, based on Venus observations and calculations drive the model flow. Dissipation is accomplished mainly by Rayleigh friction, chosen to produce strong dissipation above 85 km in order to absorb upward propagating waves and limit extreme flow velocities there, yet to give very weak Rayleigh friction below 70 km; results in the cloud layer do not appear to be sensitive to the Rayleigh friction. The model also has weak vertical diffusion, and very weak horizontal diffusion, which has a smoothing effect on the flow only at the two grid points nearest the pole.

Newman, M.↗

Collisional broadening of nonlinear resonant wave–particle interactions

A general procedure for understanding plasma behaviour when resonant wave–particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principal results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas, narrow collisional boundary layers must be retained for resonant wave–particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear processes.

Physics↗

Direct Simulation of Extinction in a Slab of Spherical Particles

The exact multiple sphere superposition method is used to calculate the coherent and incoherent contributions to the ensemble-averaged electric field amplitude and Poynting vector in systems of randomly positioned nonabsorbing spherical particles. The target systems consist of cylindrical volumes, with radius several times larger than length, containing spheres with positional configurations generated by a Monte Carlo sampling method. Spatially dependent values for coherent electric field amplitude, coherent energy flux, and diffuse energy flux, are calculated by averaging of exact local field and flux values over multiple configurations and over spatially independent directions for fixed target geometry, sphere properties, and sphere volume fraction. Our results reveal exponential attenuation of the coherent field and the coherent energy flux inside the particulate layer and thereby further corroborate the general methodology of the microphysical radiative transfer theory. An effective medium model based on plane wave transmission and reflection by a plane layer is used to model the dependence of the coherent electric field on particle packing density. The effective attenuation coefficient of the random medium, computed from the direct simulations, is found to agree closely with effective medium theories and with measurements. In addition, the simulation results reveal the presence of a counter-propagating component to the coherent field, which arises due to the internal reflection of the main coherent field component by the target boundary. The characteristics of the diffuse flux are compared to, and found to be consistent with, a model based on the diffusion approximation of the radiative transfer theory.

spheres↗

Collisional broadening of nonlinear resonant wave-particle interactions

A general procedure for understanding plasma behavior when resonant wave-particle interactions are the sole destabilizing and transport mechanism or only heating and/or current drive source is highlighted without recourse to involved numerical or analytical treatments. These phenomena are characterized by transport that appears to be collisionless even though collisions play a central role in narrow collisional boundary layers. The order of magnitude estimates, which include nonlinear effects, are shown to provide expressions in agreement with the principle results of recent toroidal Alfvén eigenmode (TAE), toroidal magnetic field ripple, and heating and current drive treatments. More importantly, the retention of nonlinearities leads to new estimates of the alpha particle energy diffusivity at saturation for TAE modes, and the ripple threshold at which superbanana plateau evaluations of alpha particle transport are modified by nonlinear radial drift effects. In addition, the estimates indicate when quasilinear descriptions for heating and current drive will begin to fail. The phenomenological procedure demonstrates that in magnetic fusion relevant plasmas narrow collisional boundary layers must be retained for resonant wave-particle interactions as they enhance the role of collisions, and make stochastic particle motion unlikely to be more important than other nonlinear process.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A block-spectral adaptive H-/$p$-refinement strategy for shock-dominated problems

An adaptive H-/p-refinement strategy using a novel sensor is devised and tested in a block-spectral compressible Euler code equipped with adaptive-mesh refinement (AMR) and high-order flux-reconstruction numerics. At each Gauss quadrature point (or solution point) within each spectral block (or mesh element) the discrete velocity jump ΔU = ∂U/∂y 1 Δy 1 + ∂V/∂y 2 Δy 2 + ∂W/∂y 3 Δy 3 is calculated and normalized by the local speed of sound, a. Here, the grid spacing, Δx i , is calculated in each direction as the distance between auxiliary Gauss-Lobatto points, staggered relative to the solution points. The polynomial order is increased from p = 0 to p = p max in regions of weak compression, (ΔU/a) crit < ΔU/a < 0 and kept at p = p max in regions of flow expansion ΔU/a ≥ 0, while staying at the H = 0 base mesh level. Regions experiencing strong compressions, i.e. ΔU/a < (ΔU/a) crit , are H-refined up to H = H max where H max is applied at the location of maximum compression, ΔU/a = min(ΔU/a) in the domain, while keeping p = 0 to guarantee robustness and monotonicity of the solution in the H refined region. The critical value of (ΔU/a) crit = -0.06 is found to effectively separate smooth and non-smooth solution regions, supported by a 1D detonation initiation test case in ideal gas and a shock-to-detonation transition in high explosives. Using this value, the Sod shock tube, Shu-Osher problem, double Mach reflection and a 2D detonation in a high-explosive are simulated with the proposed adaptive H-/p-refinement. In the Sod shock tube case, p-refinement resolves the (weak) contact discontinuity while H-refinement enhances the grid resolution in the shock exploiting the monotonicity of the p = 0 reconstruction. For the Shu-Osher problem, p-refinement captures the small-scale oscillations trailing the shock that would be otherwise attenuated, while H-refinement triggered by the ΔU-sensor appropriately tracks the shock. In the double Mach reflection problem, H-refinement confines the numerical diffusion around the reflected shock while p-refinement recaptures many physical features trailing the shock. Finally, in the 2D high-explosive detonation case, H-refinement follows the leading shock and resolves the curvature of the detonation wave, while p-refinement adds resolution to the trailing reaction zone. Finally, the proposed methodology is tested in a detonation-wave propagation test case in high-explosives with numerical predictions comparing favorably against experiments.

97 MATHEMATICS AND COMPUTING↗

Thermal conductivity of earth materials at high temperatures.

The total thermal conductivity (lattice plus radiative) of several important earth materials is measured in the temperature range from 500 to 1900 K. A new technique is used in which a CO2 laser generates a low-frequency temperature wave at one face of a small disk-shaped sample, and an infrared detector views the opposite face to detect the phase of the emerging radiation. Phase data at several frequencies yield the simultaneous determination of the thermal diffusivity and the mean extinction coefficient of the material. The lattice, radiative, and total thermal conductivities are then calculated. Results for single-crystal and polycrystalline forsterite-rich olivines and an enstatite indicate that, even in relatively pure large-grained material, the radiative conductivity does not increase rapidly with temperature. The predicted maximum total thermal conductivity at a depth of 400 km in an olivine mantle is 0.020 cal/cm/sec/deg C, which is less than twice the surface value.

Schatz, J. F.↗

The relative diffusion of strong magnetic fields and tenuous gases

General equations are given for treating the diffusion of a magnetic field in a fluid with electrical conductivity and a resistive diffusion coefficient. Consequences of nonnegligible gas motion are investigated by considering the one-dimensional problem where a tenuous isothermal gas with pressure expressed as a function of x-coordinate and time is enmeshed in a magnetic field in the y-direction. Approximate equations are obtained for the case of weak fluctuations in gas or field pressure, and the properties of a nonlinear diffusion equation are illustrated with the example of an initial thin slab of gas confined by a uniform field on either side. The advance of a gas front into a nonvanishing gas density is illustrated by using the progressive wave solution. Application of the results to the evolution of streamers and filaments observed in the solar corona is briefly discussed

Parker, E. N.↗

On the origin of the MeV energy nucleon flux associated with CIRs

Measurements of recurrent enhancements of interplanetary nucleon flux in the MeV energy range are presented and interpreted. The features recur at the solar rotation period in association with stream-stream plasma interaction regions corotating with the sun. At distances from the sun less than 1 AU, the maximum intensities of the hydrogen and helium components increase with increasing distance between 300%/AU and 600%/AU. A model is proposed which predicts the acceleration of nuclei from keV plasma energies to MeV energies by means of transit time damping of magnetosonic waves as solar wind plasma flows from the sun. Numerical solutions of the transport equation are derived to demonstrate that the model does not reproduce radial variations of hydrogen and helium fluxes, and estimates are made of the radial diffusive mean free path. Finally, the observations are found to be consistent with nucleon acceleration at the CIR shocks beyond 1 AU with subsequent diffusion toward the sun.

Christon, S. P.↗

The Shuttle Orbiter high resolution accelerometer package experiment - Preliminary flight results

A description of the High Resolution Accelerometer Package (HiRAP) experiment, designed to measure rarefied flow aerodynamic accelerations, is given. The ground test calibration factors, as well as post-flight data processing techniques to extract aerodynamic accelerations, are discussed and applied to the recorded reentry data of the recent STS-6 Shuttle Orbiter flight. The ratio of the measured normal-to-axial aerodynamic accelerations during reentry is used to obtain the first flight measurement of the lift-to-drag ratio (L/D) of a winged entry vehicle in the rarefied flow flight regime. The preliminary result for the free-molecule flow L/D is 0.10 + or - 0.03. The initial flight results on L/D are compared with current predictions which are based upon theoretical and empirical considerations. The measured free-molecule flow L/D value is higher than predicted by a factor of about three, indicating that surface reflection is not completely diffuse as currently assumed. In the rarefied flow transition regime, the Shuttle data book-bridging formula fits the flight data adequately. Upper altitude density profiles are also deduced from the measurements and presented. The density profiles show a wave phenomena with amplitude of about 60 percent, relative to a standard model.

Blanchard, R. C.↗

The generation of low-frequency electrostatic waves in association with auroral kilometric radiation

Observations of the electron distribution in the source region of auroral kilometric radiation (AKR) indicate that it can have positive gradients with respect to both v-perpendicular and v-parallel. It is shown that this type of distribution is unstable to both the electron cyclotron maser instability, which is responsible for AKR, and the bump-in-tail instability which generates electron acoustic waves in the hiss band. Simulations indicate that the two instabilities do not develop independently but compete for the available free energy due to the reduction in the positive gradients of the distribution caused by the quasi-linear diffusion associated with each instability. The dominant instability depends on the ratio of the density of the energetic electrons n(E) to that of the background electrons n(b). For n(E)/n(b) greater than about 1, the maser instability dominates, and the bump-in-tail instability is suppressed, whereas the reverse is true if n(E)/n(b) is smaller than about 1. Comparison with observations indicates that probably n(E)/n(b) is greater than about 1 in the source region of AKR.

Winglee, R. M.↗

Modulational instability of lower hybrid waves at the magnetopause

The role of lower bybrid waesat the magnetopause is reexamined. It is found that for the maximum observed wave power, the lower hybrid waves are unstable to a modulational instability on the magnetosheath side of the magnetopause. The modulational instabitlity leads to localized field structures oriented predominantly along the magnetic field. Such patchy lower hybrid turbulence has been observed by some spacecraft. As a result of the large T(sub i)/T(sub e) ratio, the waves saturate by ion heating; as a result, unlike other settings (e.g. comets, critical ionization phenomena) energetic electrons are not expected. The stochasitc electron transport in the presence of such turbulence is analyzed and results in strongly anistropic electron diffusion, with the dominant direction across the magnetic field. The diffusion rate exceeds significantly that expected from quasi-linear considerations and, for magnetospause parameters, also exceeds the rate discussed by Sonnerup (1980).

Shapiro, V. D.↗

Kinetic Alfven waves as a source of plasma transport at the dayside magnetopause

As the shocked solar wind with variable plasma density and magnetic field impinges on the dayside magnetopause, it is likely to generate large-scale Alfven waves at the solar wind magnetosphere interface. However, large gradients in the density and magnetic field at the magnetopause boundary effectively couple large-scale Alfven waves with kinetic Alfven waves may play an important role in plasma transport at the dayside magnetopause and in electron acceleration along field lines. The transport can occur because, unlike the magnetohydrodynamic (MHD) shear Alfven wave, the kinetic Alfven wave has an associated parallel electric field which breaks down the 'frozen-in' condition and decouples the plasma from field lines. We calculate the average deviation of the plasma from the field line from which we estimate the diffusion coefficient associated with these 'bundles' of lead to acceleration of electrons along field lines in the magnetopause boundary and may possibly provide an explanation for observed counterstreaming electron beams characterized by energies of 50-200 eV.

Lee, L. C.↗

A Study of Supernova Remnants with Center-Filled X-Ray Morphology

CTA 1 is a center-filled supernova remnant (SNR) whose morphology and spectrum indicate the presence of a central pulsar, a synchrotron nebula, and a thermal component associated with the expansion of the blast wave into the interstellar medium. The centrally bright emission surrounds the position of a faint point source of x-rays observed with the ROSAT PSPC. Here we report on ASCA observations that confirm the nonthermal nature of the diffuse emission from the central regions of the remnant. We also present evidence for weak thermal emission that appears to increase in strength toward the outer boundary of the SNR. Thus, CTA 1 appears to be an x-ray composite remnant. Both the aftermath of the explosive supernova event and the energetic compact core are observable.

Slane, Patrick O.↗

Stable and low diffusive hybrid upwind splitting methods

A new concept for upwinding is introduced, named the hybrid upwind splitting (HUS), which is achieved by combining the basically distinct flux vector splitting (FVS) and the flux difference splitting (FDS) approaches. The HUS approach yields upwind methods which share the robustness of the FVS schemes in the capture of nonlinear waves and the accuracy of some of the FDS schemes. Numerical illustrations are presented proving the relevance of the HUS methods for viscous calculations.

Coquel, Frederic↗

Vertical constituent transport in the mesosphere

Ground-based microwave spectroscopy measurements of mesospheric CO and H2O vertical mixing ratio profiles are used to infer vertical mixing rates in the upper mesosphere. The CO and H2O data consistently imply vertical eddy diffusion coefficients in the 70- to 85-km region of 100,000-200,000 sq cm/s during spring through summer at midlatidues. Although chemical acceleration of vertical transport is substantial for O and O3, below the mesopause, the divergences of their associated fluxes are modest, with at most a factor of 2 effect on the concentrations of O and O3 for measured variability in gravity wave activity. Comparison of Solar Mesosphere Explorer (SME) O3 data with model results reinforces the conclusions of slow vertical mixing in the upper mesosphere as a consequence of the reduced HO(x) catalytic loss of odd oxygen. The changes in chemical rate constants recommended by Rusch and Eckman (1985), in conjunction with slow vertical mixing, yield good agreement with SME O3 data. The slow vertical mixing deduced in this study is consistent with upper limits obtained from studies of the mesospheric heat budget and could be construed as evidence for an advectively controlled mesosphere. A comparison of the vertical eddy diffusion coefficients for momentum stresses, constituent transport, and heat transport suggests that the eddy Prandtl number must be of order 10.

Strobel, Darrell F.↗

Computational diagnostics and characterization of combustion recession in diesel sprays

While low-temperature combustion (LTC) strategies have been found to mitigate nitrogen oxides and particulate matter emissions in diesel engines, studies have also associated LTC with an increase in unburnt hydrocarbons. With more recent studies on diesel after end-of-injection (AEOI), combustion recession is identified as a phenomenon where at near nozzle region, high-temperature ignition (HTI) combustion can propagate back to the nozzle tip consuming the unburnt hydrocarbons AEOI. Current literature has suggested that combustion recession is controlled by auto-ignition. However, high-fidelity simulations and detailed analysis of such a mechanism are missing. In this study, comprehensive Large Eddy Simulations of a reacting spray at “Spray A” conditions are performed, where detailed analysis of combustion recession concerning flame morphology and propagation modes are included. In particular, this study demonstrated for the first time that while combustion recession is mainly auto-ignition dominated (consistent with the literature), a cool flame was found to deflagrate towards the richer regions of the mixture, promoting mixing and increasing the mixture temperature. This leads to HTI kernels, which then grow and develop as deflagrative waves, therefore sustaining the combustion recession process. The study also detailed the extinction mechanism of combustion: the entrainment wave will overlean the near-nozzle mixtures, rendering it unable to support HTI, which leads to the extinction of the upstream flame AEOI in lower reactivity mixtures. Combustion recession is also observed to be contingent on the chemical and diffusion processes, even at low scalar dissipation rates. Finally, a new criterion for combustion recession based on chemical explosive mode is proposed and validated with previous combustion recession index to quantify the extent of HTI in near-nozzle mixtures AEOI. The newly developed metric combined with a previous experimentally-based metric can provide simple but valuable measurements of the degree and propensity of the upstream flame AEOI.

CEMA↗