Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized gradient approximations”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Multifractal characterizations of nonstationarity and intermittency in geophysical fields: Observed, retrieved, or simulated

Geophysical data rarely show any smoothness at any scale, and this often makes comparison with theoretical model output difficult. However, highly fluctuating signals and fractal structures are typical of open dissipative systems with nonlinear dynamics, the focus of most geophysical research. High levels of variability are excited over a large range of scales by the combined actions of external forcing and internal instability. At very small scales we expect geophysical fields to be smooth, but these are rarely resolved with available instrumentation or simulation tools; nondifferentiable and even discontinuous models are therefore in order. We need methods of statistically analyzing geophysical data, whether measured in situ, remotely sensed or even generated by a computer model, that are adapted to these characteristics. An important preliminary task is to define statistically stationary features in generally nonstationary signals. We first discuss a simple criterion for stationarity in finite data streams that exhibit power law energy spectra and then, guided by developments in turbulence studies, we advocate the use of two ways of analyzing the scale dependence of statistical information: singular measures and qth order structure functions. In nonstationary situations, the approach based on singular measures seeks power law behavior in integrals over all possible scales of a nonnegative stationary field derived from the data, leading to a characterization of the intermittency in this (gradient-related) field. In contrast, the approach based on structure functions uses the signal itself, seeking power laws for the statistical moments of absolute increments over arbitrarily large scales, leading to a characterization of the prevailing nonstationarity in both quantitative and qualitative terms. We explain graphically, step by step, both multifractal statistics which are largely complementary to each other. The geometrical manifestations of nonstationarity and intermittency, 'roughness' and 'sparseness', respectively, are illustrated and the associated analytical (differentiability and continuity) properties are discussed. As an example, the two techniques are applied to a series of recent measurements of liquid water distributions inside marine stratocumulus decks; these are found to be multifractal over scales ranging from approximately 60 m to approximately 60 km. Finally, we define the 'mean multifractal plane' and show it to be a simple yet comprehensive tool with many applications including data intercomparison, (dynamical or stochastic) model and retrieval validations.

Davis, Anthony↗
DOE OSTI2024

An extended variational method for the resistive wall mode in toroidal plasma confinement devices

The external-kink stability of a toroidal plasma surrounded by a rigid resistive wall is investigated. The well-known analysis of Haney and Freidberg is rigorously extended to allow for a wall that is sufficiently thick that the thin-shell approximation does not necessarily hold. A generalized Haney–Freidberg formula for the growth-rate of the resistive wall mode is obtained. Thick-wall effects do not change the marginal stability point of the mode but introduce an interesting asymmetry between growing and decaying modes. Growing modes have growth-rates that exceed those predicted by the original Haney–Freidberg formula. On the other hand, decaying modes have decay-rates that are less than those predicted by the original formula. The well-known Hu–Betti formula for the rotational stabilization of the resistive wall mode is also generalized to take thick-wall effects into account. Increasing wall thickness facilitates the rotational stabilization of the mode, because it decreases the critical toroidal electromagnetic torque that the wall must exert on the plasma. On the other hand, the real frequency of the mode at the marginal stability point increases with increasing wall thickness.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

'Downward control' of the mean meridional circulation and temperature distribution of the polar winter stratosphere

According to the 'downward control' principle, the extratropical mean vertical velocity on a given pressure level is approximately proportional to the meridional gradient of the vertically integrated zonal force per unit mass exerted by waves above that level. In this paper, a simple numerical model that includes parameterizations of both planetary and gravity wave breaking is used to explore the influence of gravity wave breaking in the mesosphere on the mean meridional circulation and temperature distribution at lower levels in the polar winter stratosphere. The results of these calculations suggest that gravity wave drag in the mesosphere can affect the state of the polar winter stratosphere down to altitudes below 30 km. The effect is most important when planetary wave driving is relatively weak: that is, during southern winter and in early northern winter. In southern winter, downwelling weakens by a factor of 2 near the stratospause and by 20% at 30 km when gravity wave drag is not included in the calculations. As a consequence, temperatures decrease considerably throughout the polar winter stratosphere (over 20 K above 40 km and as much as 8 K at 30 km, where the effect is enhanced by the long radiative relaxation timescale). The polar winter states obtained when gravity wave drag is omitted in this simple model resemble the results of simulations with some general circulation models and suggest that some of the shortcomings of the latter may be due to a deficit in mesospheric momentum deposition by small-scale gravity waves.

Garcia, Rolando R.↗
DOE OSTI2021

Non-conformal attractor in boost-invariant plasmas

We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modeled and observed properties related to the direct aerosol radiative effect of biomass burning aerosol over the southeastern Atlantic

Biomass burning smoke is advected over the southeastern Atlantic Ocean between July and October of each year. This smoke plume overlies and mixes into a region of persistent low marine clouds. Model calculations of climate forcing by this plume vary significantly in both magnitude and sign. NASA EVS-2 (Earth Venture Suborbital-2) ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) had deployments for field campaigns off the west coast of Africa in 3 consecutive years (September 2016, August 2017, and October 2018) with the goal of better characterizing this plume as a function of the monthly evolution by measuring the parameters necessary to calculate the direct aerosol radiative effect. Here, this dataset and satellite retrievals of cloud properties are used to test the representation of the smoke plume and the underlying cloud layer in two regional models (WRF-CAM5 and CNRM-ALADIN) and two global models (GEOS and UM-UKCA). The focus is on the comparisons of those aerosol and cloud properties that are the primary determinants of the direct aerosol radiative effect and on the vertical distribution of the plume and its properties. The representativeness of the observations to monthly averages are tested for each field campaign, with the sampled mean aerosol light extinction generally found to be within 20 % of the monthly mean at plume altitudes. When compared to the observations, in all models, the simulated plume is too vertically diffuse and has smaller vertical gradients, and in two of the models (GEOS and UM-UKCA), the plume core is displaced lower than in the observations. Plume carbon monoxide, black carbon, and organic aerosol masses indicate underestimates in modeled plume concentrations, leading, in general, to underestimates in mid-visible aerosol extinction and optical depth. Biases in mid-visible single scatter albedo are both positive and negative across the models. Observed vertical gradients in single scatter albedo are not captured by the models, but the models do capture the coarse temporal evolution, correctly simulating higher values in October (2018) than in August (2017) and September (2016). Uncertainties in the measured absorption Ångstrom exponent were large but propagate into a negligible (<4 %) uncertainty in integrated solar absorption by the aerosol and, therefore, in the aerosol direct radiative effect. Model biases in cloud fraction, and, therefore, the scene albedo below the plume, vary significantly across the four models. The optical thickness of clouds is, on average, well simulated in the WRF-CAM5 and ALADIN models in the stratocumulus region and is underestimated in the GEOS model; UM-UKCA simulates cloud optical thickness that is significantly too high. Overall, the study demonstrates the utility of repeated, semi-random sampling across multiple years that can give insights into model biases and how these biases affect modeled climate forcing. The combined impact of these aerosol and cloud biases on the direct aerosol radiative effect (DARE) is estimated using a first-order approximation for a subset of five comparison grid boxes. A significant finding is that the observed grid box average aerosol and cloud properties yield a positive (warming) aerosol direct radiative effect for all five grid boxes, whereas DARE using the grid-box-averaged modeled properties ranges from much larger positive values to small, negative values. It is shown quantitatively how model biases can offset each other, so that model improvements that reduce biases in only one property (e.g., single scatter albedo but not cloud fraction) would lead to even greater biases in DARE. Across the models, biases in aerosol extinction and in cloud fraction and optical depth contribute the largest biases in DARE, with aerosol single scatter albedo also making a significant contribution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1995

Navier-Stokes Aerodynamic Simulation of the V-22 Osprey on the Intel Paragon MPP

The paper will describe the Development of a general three-dimensional multiple grid zone Navier-Stokes flowfield simulation program (ENS3D-MPP) designed for efficient execution on the Intel Paragon Massively Parallel Processor (MPP) supercomputer, and the subsequent application of this method to the prediction of the viscous flowfield about the V-22 Osprey tiltrotor vehicle. The flowfield simulation code solves the thin Layer or full Navier-Stoke's equation - for viscous flow modeling, or the Euler equations for inviscid flow modeling on a structured multi-zone mesh. In the present paper only viscous simulations will be shown. The governing difference equations are solved using a time marching implicit approximate factorization method with either TVD upwind or central differencing used for the convective terms and central differencing used for the viscous diffusion terms. Steady state or Lime accurate solutions can be calculated. The present paper will focus on steady state applications, although time accurate solution analysis is the ultimate goal of this effort. Laminar viscosity is calculated using Sutherland's law and the Baldwin-Lomax two layer algebraic turbulence model is used to compute the eddy viscosity. The Simulation method uses an arbitrary block, curvilinear grid topology. An automatic grid adaption scheme is incorporated which concentrates grid points in high density gradient regions. A variety of user-specified boundary conditions are available. This paper will present the application of the scalable and superscalable versions to the steady state viscous flow analysis of the V-22 Osprey using a multiple zone global mesh. The mesh consists of a series of sheared cartesian grid blocks with polar grids embedded within to better simulate the wing tip mounted nacelle. MPP solutions will be shown in comparison to equivalent Cray C-90 results and also in comparison to experimental data. Discussions on meshing considerations, wall clock execution time, load balancing, and scalability will be provided.

Vadyak, Joseph↗
NASA NTRS2009

High Temperature VARTM of Phenylethynyl Terminated Imides

Fabrication of composite structures using vacuum assisted resin transfer molding (VARTM) is generally more affordable than conventional autoclave techniques. Recent efforts have focused on adapting VARTM for the fabrication of high temperature composites. Due to their low melt viscosity and long melt stability, certain phenylethynyl terminated imides (PETI) can be processed into composites using high temperature VARTM (HT-VARTM). However, one of the disadvantages of the current HT-VARTM resin systems has been the high porosity of the resultant composites. For aerospace applications, the desired void fraction of less than 2% has not yet been achieved. In the current study, two PETI resins, LaRC PETI-330 and LaRC PETI-8 have been used to make test specimens using HT-VARTM. The resins were infused into ten layers of IM7-6K carbon fiber 5-harness satin fabric at 260 C or 280 C and cured at 371 C. Initial runs yielded composites with high void content, typically greater than 7% by weight. A thermogravimetric-mass spectroscopic study was conducted to determine the source of volatiles leading to high porosity. It was determined that under the thermal cycle used for laminate fabrication, the phenylethynyl endcap was undergoing degradation leading to volatile evolution. By modifying the thermal cycle used in laminate fabrication, the void content was reduced significantly (typically approximately 3%). Densities of the composites were determined using a density gradient column and the glass transition temperatures of the cured composites were measured by dynamic mechanical analysis. Photomicrographs of the panels were taken and void contents were determined by acid digestion. The results of this work are presented herein.

Cano, Roberto J.↗
DOE OSTI2022

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

On the temperature and velocity through the photosphere of a sunspot penumbra

We investigate the structure in depth of a sunspot penumbra by means of the inversion code of the radiative transfer equation proposed by Ruiz Cobo & del Toro Iniesta (1992), applied to a set of filtergrams of a sunspot, scanning the Fe I line at 5576.1 A, with a sampling interval of 30 mA, from -120 to 120 mA from line center (data previously analyzed by Title et al. 1993). The temperature structure of this penumbra is obtained for each of the 801 pixels selected (0.32 sec x 0.32 sec). On the average, the temperatures seem to decrease as we move inward, but the differences are of the order of the rms values (approximately equal 100-200 K) at a given distance to sunspot center. The outer parts of the penumbra have also a bigger curvature in the T versus log tau(sub 5) relation than the inner parts. We realize, however, that these differences might be influenced by possible stray light effects. Compared to the quiet Sun, penumbral temperatures are cooler at deep layers and hotter at high layers. A mean penumbral model atmosphere is presented. The asymmetries observed in the intensity profile (the line is magnetically insensitive) are deduced to be produced by strong gradients of the line-of-sight velocity that sharply vary spatially along slices of almost constant distance to sunspot center. These variations suggest that such gradients are not only needed to explain the broadband circular polarization observed in sunspots (see Sanchez Almeida & Lites 1992) but are a main characteristic of the fine-scale penumbra. The results are compatible with an Evershed flow present everywhere, but its gradient with depth turns out to vary so that the flow seems to be mainly concentrated in some penumbral fibrils when studied through Dopplergrams. Finally, as by-products of this study, we put constraints to the practical usefulness of the Eddington-Barbier relation, and we explain the values of the Fourier Dopplergrams to be carrying information of layers around the centroid of the generalized response function of Dopplergrams to velocity fluctuations.

Del Toro Iniesta, J. C.↗
DOE OSTI2023

Adaptive, problem-tailored variational quantum eigensolver mitigates rough parameter landscapes and barren plateaus

Abstract Variational quantum eigensolvers (VQEs) represent a powerful class of hybrid quantum-classical algorithms for computing molecular energies. Various numerical issues exist for these methods, however, including barren plateaus and large numbers of local minima. In this work, we consider the Adaptive, Problem-Tailored Variational Quantum Eiegensolver (ADAPT-VQE) ansätze, and examine how they are impacted by these local minima. We find that while ADAPT-VQE does not remove local minima, the gradient-informed, one-operator-at-a-time circuit construction accomplishes two things: First, it provides an initialization strategy that can yield solutions with over an order of magnitude smaller error compared to random initialization, and which is applicable in situations where chemical intuition cannot help with initialization, i.e., when Hartree-Fock is a poor approximation to the ground state. Second, even if an ADAPT-VQE iteration converges to a local trap at one step, it can still “burrow” toward the exact solution by adding more operators, which preferentially deepens the occupied trap. This same mechanism helps highlight a surprising feature of ADAPT-VQE: It should not suffer optimization problems due to barren plateaus and random initialization. Even if such barren plateaus appear in the parameter landscape, our analysis suggests that ADAPT-VQE avoids such regions by design.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1996

Regulation of Methane Oxidation in a Freshwater Wetland by Water Table Changes and Anoxia

The effects of water table fluctuations and anoxia on methane emission and methane oxidation were studied in a freshwater marsh. Seasonal aerobic methane oxidation rates varied between 15% and 76% of the potential diffusive methane flux (diffusive flux in the absence of aerobic oxidation). On an annual basis, approximately 43% of the methane diffusing into the oxic zone was oxidized before reaching the atmosphere. The highest methane oxidation was observed when the water table was below the peat surface. This was confirmed in laboratory experiments where short-term decreases in water table levels increased methane oxidation but also net methane emission. Although methane emission was generally not observed during the winter, stems of soft rush (Juncus effusus) emitted methane when the marsh was ice covered. Indigenous methanotrophic bacteria from the wetiand studied were relatively anoxia tolerant. Surface peat incubated under anoxic conditions maintained 30% of the initial methane oxidation capacity after 32 days of anoxia. Methanotrophs from anoxic peat initiated aerobic methane oxidation relatively quickly after oxygen addition (1-7 hours). These results were supported by culture experiments with the methanotroph Methylosinus trichosporium OB3b. This organism maintained a greater capacity for aerobic methane oxidation when starved under anoxic compared to oxic conditions. Anoxic incubation of M. trichosporium OB3b in the presence of sulfide (2 mM) and a low redox potential (-110 mV) did not decrease the capacity for methane oxidation relative to anoxic cultures incubated without sulfide. The results suggest that aerobic methane oxidation was a major regulator of seasonal methane emission front the investigated wetland. The observed water table fluctuations affected net methane oxidation presumably due to associated changes in oxygen gradients. However, changes from oxic to anoxic conditions in situ had relatively little effect on survival of the methanotrophic bacteria and thus on methane oxidation potential per se.

Roslev, Peter↗
NASA NTRS2007

Effect of Head-Down Bed Rest and Artificial Gravity Countermeasure on Cardiac Autonomic and Advanced Electrocardiographic Function

To study the effects of 21 days of head-down bed rest (HDBR), with versus without an artificial gravity (AG) countermeasure, on cardiac autonomic and advanced electrocardiographic function. Fourteen healthy men participated in the study: seven experienced 21 days of HDBR alone ("HDBR controls") and seven the same degree and duration of HDBR but with approximately 1hr daily short-arm centrifugation as an AG countermeasure ("AG-treated"). Five minute supine high-fidelity 12-lead ECGs were obtained in all subjects: 1) 4 days before HDBR; 2) on the last day of HDBR; and 3) 7 days after HDBR. Besides conventional 12-lead ECG intervals and voltages, all of the following advanced ECG parameters were studied: 1) both stochastic (time and frequency domain) and deterministic heart rate variability (HRV); 2) beat-to-beat QT interval variability (QTV); 3) T-wave morphology, including signal-averaged T-wave residua (TWR) and principal component analysis ratios; 4) other SAECG-related parameters including high frequency QRS ECG and late potentials; and 5) several advanced ECG estimates of left ventricular (LV) mass. The most important results by repeated measures ANOVA were that: 1) Heart rates, Bazett-corrected QTc intervals, TWR, LF/HF power and the alpha 1 of HRV were significantly increased in both groups (i.e., by HDBR), but with no relevant HDBR*group differences; 2) All purely "vagally-mediated" parameters of HRV (e.g., RMSSD, HF power, Poincare SD1, etc.), PR intervals, and also several parameters of LV mass (Cornell and Sokolow-Lyon voltages, spatial ventricular activation times, ventricular gradients) were all significantly decreased in both groups (i.e., by HDBR), but again with no relevant HDBR*group differences); 3) All "generalized" or "vagal plus sympathetic" parameters of stochastic HRV (i.e., SDNN, total power, LF power) were significantly more decreased in the AG-treated group than in the HDBR-only group (i.e., here there was a relevant HDBR*group difference); and 4) QTV index was also significantly more changed (increased) in the AG-treated group than in the HDBR-only group, although this was clearly due to a greater decrease in generalized HRV and not to a greater increase in QTV proper because there was no relevant HDBR*group effect for either the SDNN or the RMSSD of QTV. Brief daily AG treatment by short-arm centrifuge during each of 21 days of HDBR does not appear to protect against HDBR-related losses of cardiac autonomic function or of LV mass as estimated by ECG.

Schlegel, T. T.↗
NASA NTRS2024

Global Optimization via Quadratic Disjunctive Programming for Water Networks Design with Energy Recovery

Generalized disjunctive programming (GDP) models with bilinear and concave constraints, often seen in water network design, are challenging optimization problems. This work proposes quadratic and piecewise linear approximations for nonlinear terms to reformulate GDP models into quadratic GDP (QGDP) models that suitable solvers may solve more efficiently. We illustrate the benefits of the quadratic reformulation with a water treatment network design problem in which nonconvexities arise from bilinear terms in the mixers’ mass balances and concave investment cost functions of treatment units. Given the similarities with water network design problems, we suggest quadratic approximation for the GDP model for the optimal design of a large-scale reverse electrodialysis (RED) process. This power technology can recover energy from salinity differences between by-product streams of the water sector, such as desalination brine mixed with regenerated wastewater effluents. The solver Gurobi excels in handling QGDP problems, but weighing the problem’s precision and tractability balance is crucial. The piecewise linear approximation yields more accurate, yet larger QGDP models that may require longer optimization times in large-scale process synthesis problems.

Water Networks↗
DOE OSTI2023

Generalized moving least squares vs. radial basis function finite difference methods for approximating surface derivatives

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. Further, the second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

Galileo Direction Finding of Jovian Radio Emissions

The Galileo spacecraft, in orbit about Jupiter, has observed distinct spin modulation of plasma wave emissions near the Ganymede (G1 and G2) encounters in the frequency range from about 100 kHz to approximately 6 MHz. Assuming circularly polarized, transverse electromagnetic radiation, we have used the spin modulation of the sweep-frequency receivers of the electric dipole antenna over many spins to estimate the source location in the spin plane of the spacecraft. Hectometric (HOM) and decametric (DAM) emission is observed by Galileo as a general and continuous background with frequent bursts that last tens of minutes and can be separated by minutes or hours. We have analyzed HOM and DAM emissions observed near Jupiter just after the GI and G2 encounters, including two HOM/DAM "arc" signatures observed after the G2 encounter. These latter appear to be low-frequency extensions of DAM arcs, with source regions along either the Io or the Ganymede flux tube. While the uncertainties associated with the data analysis do not allow a precise source location, the HOM/DAM emission observed near the G1 and G2 encounters is consistent with a gyroresonant source region, but it is necessary to require refraction due to the Io torus to understand the results. To explain emission from apparent source regions above a gyroresonant source region, wave refraction from asymmetries in the Io plasma torus that extend along magnetic field lines is postulated. Alternatively, if such torus density asymmetries do not exist, emission with sources above a gyroresonant source region would require another free-energy source such as energetic plasma beams in the presence of density gradients or temperature anisotropies.

Menietti, J. D.↗
NASA NTRS1999

Mars Global Surveyor Thermal Emission Spectrometer (TES) Observations: Atmospheric Temperatures During Aerobraking and Science Phasing

Between September 1997, when the Mars Global Surveyor spacecraft arrived at Mars, and September 1998 when the final aerobraking phase of the mission began, the Thermal Emission Spectrometer (TES) has acquired an extensive data set spanning approximately half of a Martian year. Nadir-viewing spectral measurements from this data set within the 15-micrometers CO2 absorption band are inverted to obtain atmospheric temperature profiles from the surface up to about the 0.1 mbar level. The computational procedure used to retrieve the temperatures is presented. Mean meridional cross sections of thermal structure are calculated for periods of time near northern hemisphere fall equinox, winter solstice, and spring equinox, as well as for a time interval immediately following the onset of the Noachis Terra dust storm. Gradient thermal wind cross sections are calculated from the thermal structure. Regions of possible wave activity are identified using cross sections of rms temperature deviations from the mean. Results from both near-equinox periods show some hemispheric asymmetry with peak eastward thermal winds in the north about twice the magnitude of those in the south. The results near solstice show an intense circumpolar vortex at high northern latitudes and waves associated with the vortex jet core. Warming of the atmosphere aloft at mid-northern latitudes suggests the presence of a strong cross-equatorial Hadley circulation. Although the Noachis dust storm did not become global in scale, strong perturbations to the atmospheric structure are found, including an enhanced temperature maximum aloft at high northern latitudes resulting from intensification of the Hadley circulation. TES results for the various seasonal conditions are compared with published results from Mars general circulation models, and generally good qualitative agreement is found.

Conrath, Barney J.↗
DOE OSTI2025

Dependence of ion-cyclotron range of frequencies instabilities on species mix and fast-ion distribution: I. Stability

Ion cyclotron emission (ICE) driven by fast ions has the potential to be a reactor relevant diagnostic of confined and escaping fast ions; moreover, controlled experiments in magnetic fusion devices can shed light on similar magnetospheric instabilities below and above the ion cyclotron frequency f ci . Dedicated DIII-D experiments studied instabilities below f ci and at cyclotron harmonics l f ci by injecting thirteen different neutral beam populations into plasmas with different thermal compositions of H, D, and 3He at five different values of toroidal field. High-frequency magnetic loops diagnose mode properties. Each beam type is individually injected for an approximate slowing-down time (~100 ms), resulting in a database of 2529 conditions that are analyzed in frequency bands corresponding to ICE harmonics and < f ci global Alfvén eigenmodes. The sub-cyclotron modes are more unstable at low field and in plasmas with large H concentration. For D injection, first harmonic ICE at f cD is more unstable in hydrogen plasmas but, for H injection, the converse is true (ICE at f cH is more unstable in D), suggesting that thermal same-species ions damp fundamental ICE. ICE stability also depends upon the pitch-angle anisotropy gradient and upon the ‘bump-on-tail’ gradient ∂f/∂v. Although the frequency of unstable ICE in these plasmas always occurs near the central cyclotron frequency, the precise value also depends on the Alfvén speed and beam direction, consistent with the idea that instability occurs for modes that simultaneously satisfy ω ≃ 2πlf ci , the magnetosonic dispersion relation, and the Doppler-shifted resonance condition.

electromagnetic ion cyclotron instabilities↗
NASA NTRS1996

Theoretical Basis for Finite Difference Extrapolation of Sonic Boom Signatures

Calculation of sonic boom signatures for aircraft has traditionally followed the methods of Whitham' and Walkden. The wave disturbance generated by the vehicle is obtained by area rule linearized supersonic flow methods, which yield a locally axisymmetric asymptotic solution. This solution is acoustic in nature, i.e., first order in disturbance quantities, and corresponds to ray acoustics. Cumulative nonlinear distortion of the signature is incorporated by using this solution to adjust propagation speed to first order, thus yielding a solution second order in disturbance quantities. The effects of atmospheric gradients are treated by Blokhintzov's method of geometrical acoustics. Both nonlinear signature evolution and ray tracing are applied as if the pressure field very close to the vehicle were actually that given by the source term (the 'F-function') of the asymptotic linearized flow solution. The viewpoint is thus that the flow solution exists at a small radius near the vehicle, and may be treated as an input to an extrapolation procedure consisting of ray tracing and nonlinear aging. The F-function is often regarded as a representation of a near-field pressure signature, and it is common for computational implementations to treat it interchangeably with the pressure signature. There is a 'matching radius' between the source function and the subsequent propagation extrapolation. This viewpoint has been supported by wind tunnel tests of simple models, and very typically yields correct results for actual flight vehicles. The assumption that the F-function and near-field signature are interchangeable is generally not correct. The flowfield of a vehicle which is not axisymmetric contains crossflow components which are very significant at small radii and less so at larger distances. From an acoustical viewpoint, the crossflow is equivalent to source diffraction portions of the wave field. Use of the F-function as a near field signature effectively assumes that the diminution of the crossflow/diffraction component may be applied all at once at the matching radius noted above. This approximation, though not rigorously validated, is responsible for the usual correct far-field results. On the other hand, if an actual near-field signature (either from wind tunnel or CFD data) is used at a starting point rather than one based on th effective source distribution, the predicted far-field signature is generally wrong.

Plotkin, Kenneth J.↗