SEARCH · Search NASA

Results for “variable density flow”

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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗

NASA NTRS1993

Observational evidence of preferred flow regimes in the Northern Hemisphere winter stratosphere

Ten years of stratospheric geopotential height data are analyzed in an attempt to determine whether there are preferred flow regimes in the Northern Hemisphere winter stratosphere. The data are taken from Stratospheric Sounding Units on board NOAA satellites. The probability density estimate of the amplitude of the wavenumber 1 10-mb height is found to be bimodal. The density distribution is composed of a dominant large-amplitude mode and a less frequent low-amplitude mode. When the wavenumber 1 10-mb height data are projected onto the phase plane defined by the 10-mb zonal-mean winds and wavenumber 1 100-mb heights, three preferred regimes are evident. The small-amplitude mode separates into a strong zonal wind-weak wave regime and a weak zonal wind-weak wave regime. The large-amplitude mode is an intermediate zonal wind-strong wave regime. Transitions between the large-amplitude regime and the weak zonal wind-weak wave regime are found to be associated with major stratospheric warmings. The clustering of the stratospheric data into the preferred flow regimes is interpreted in light of the bifurcation properties of the Holton and Mass model. The interannual variability of the Northern Hemisphere winter stratosphere is interpreted in terms of the relative frequency of the observed preferred regimes.

Pierce, R. B.↗

NASA NTRS2015

Stage-by-Stage and Parallel Flow Path Compressor Modeling for a Variable Cycle Engine, NASA Advanced Air Vehicles Program - Commercial Supersonic Technology Project - AeroServoElasticity

This paper covers the development of stage-by-stage and parallel flow path compressor modeling approaches for a Variable Cycle Engine. The stage-by-stage compressor modeling approach is an extension of a technique for lumped volume dynamics and performance characteristic modeling. It was developed to improve the accuracy of axial compressor dynamics over lumped volume dynamics modeling. The stage-by-stage compressor model presented here is formulated into a parallel flow path model that includes both axial and rotational dynamics. This is done to enable the study of compressor and propulsion system dynamic performance under flow distortion conditions. The approaches utilized here are generic and should be applicable for the modeling of any axial flow compressor design accurate time domain simulations. The objective of this work is as follows. Given the parameters describing the conditions of atmospheric disturbances, and utilizing the derived formulations, directly compute the transfer function poles and zeros describing these disturbances for acoustic velocity, temperature, pressure, and density. Time domain simulations of representative atmospheric turbulence can then be developed by utilizing these computed transfer functions together with the disturbance frequencies of interest.

Compressor Modeling Modeling↗

NASA NTRS2006

NASA Tech Briefs, July 2006

Topics covered include: Airport Remote Tower Sensor Systems; Implantable Wireless MEMS Sensors for Medical Uses; Embedded Sensors for Measuring Surface Regression; Coordinating an Autonomous Earth-Observing Sensorweb; Range-Measuring Video Sensors; Stability Enhancement of Polymeric Sensing Films Using Fillers; Sensors for Using Times of Flight to Measure Flow Velocities; Receiver Would Control Phasing of a Phased-Array Antenna; Modern Design of Resonant Edge-Slot Array Antennas; Carbon-Nanotube Schottky Diodes; Simplified Optics and Controls for Laser Communications; Coherent Detection of High-Rate Optical PPM Signals; Multichannel Phase and Power Detector; Using Satellite Data in Weather Forecasting: I; Using Dissimilarity Metrics to Identify Interesting Designs; X-Windows PVT Widget Class; Shuttle Data Center File-Processing Tool in Java; Statistical Evaluation of Utilization of the ISS; Nanotube Dispersions Made With Charged Surfactant; Aerogels for Thermal Insulation of Thermoelectric Devices; Low-Density, Creep-Resistant Single-Crystal Superalloys; Excitations for Rapidly Estimating Flight-Control Parameters; Estimation of Stability and Control Derivatives of an F-15; Tool for Coupling a Torque Wrench to a Round Cable Connector; Ultrasonically Actuated Tools for Abrading Rock Surfaces; Active Struts With Variable Spring Stiffness and Damping; Multiaxis, Lightweight, Computer-Controlled Exercise System; Dehydrating and Sterilizing Wastes Using Supercritical CO2; Alpha-Voltaic Sources Using Liquid Ga as Conversion Medium; Ice-Borehole Probe; Alpha-Voltaic Sources Using Diamond as Conversion Medium; White-Light Whispering-Gallery-Mode Optical Resonators; Controlling Attitude of a Solar-Sail Spacecraft Using Vanes; and Wire-Mesh-Based Sorber for Removing Contaminants from Air.

Source record↗

NASA NTRS2021

Predicting near-saturated hydraulic conductivity in urban soils

Pedotransfer functions (PTFs) provide point predictions of soil hydraulic properties from more readily measured soil characteristics, yet uncertainties and biases in measurement methods, sampling distributions, and boundary conditions can limit accuracy when estimating near-saturated hydraulic conductivity (K(n)). These limitations may be particularly problematic in understudied urban landscapes that often contain altered hydraulic properties. To better treat deficiencies in PTF performance, we addressed three objectives, which were to: 1) develop PTFs to predict urban K(n), 2) assess bulk density and coarse fragments as explanatory variables; and 3) evaluate the predictive capability of these PTFs by comparing their output to measured hydraulic conductivity values from three other studies of urban soil hydraulics. We used artificial neural networks (ANN) and random forest (RF) approaches to predict urban K(n), with the training dataset including 307 tension infiltrometer tests and other measurements drawn from urban soil assessments in 11 U.S. cities. The PTFs utilized a hierarchy of inputs, starting with percentage sand, silt, clay, and then adding percentage coarse fragments and bulk density. The ANN models performed similar to the RF models, and all models exhibited similar or better predictive performance as models results collected from published articles. The inclusion of bulk density or coarse fragments did not improve accuracy over soil texture alone. Possible reasons for this result include low correlation between K(n) and bulk density and the exclusion of large voids during flow measurements with tension infiltrometers. The models have been made available as an open-source software package to encourage adoption by users working in urban systems.

Jinshi Jian↗

NASA NTRS1982

High resolution measurements of nightside ion troughs at Venus - Evidence of electrodynamic perturbations

The Bennett rf ion mass spectrometer of the Pioneer Venus Orbiter was expressly designed to provide variable temporal resolution for measurements of thermal ion composition and density. The Explore-Adapt mode is used to obtain priority for measuring the most prominent ion species; in the 2/16 configuration, the two dominant ions within the available range of 16 species are selectively sampled at the highest rate of 0.2 sec/sample. The high-resolution measurements are combined with independent observations from the magnetic field, neutral mass spectrometer, and electron temperature experiments in investigating sharply structured troughs in the low-altitude nightside ion concentrations. The results suggest a close correlation between the structure in the ion distributions and the structured configuration of the magnetic field that is draped about the planet. In the regions of the ion depletions, sharp fluctuations in electron temperature and anomalous increases in the density of neutral gases suggest that the ion depletion may be associated with dynamic perturbation in the ion and neutral flows and/or local joule heating.

Taylor, H. A., Jr.↗

NASA NTRS2003

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗

NASA NTRS2010

Internal Acoustics of a Pintle Valve with Supercritical Helium Flow

Large amplitude flow unsteadiness is a common phenomenon within the high flow rate ducts and valves associated with propulsion systems. Boundary layer noise, shear layers and vortex shedding are a few of the many sources of flow oscillations. The presence of lightly damped acoustic modes can organize and amplify these sources of flow perturbation, causing undesirable loading of internal parts. The present study investigates the self-induced acoustic environment within a pintle valve subject to high Reynolds Number flow of helium gas. Experiments were conducted to measure the internal pressure oscillations of the Ares I Launch Abort System (LAS) Attitude Control Motor (ACM) valve. The AGM consists of a solid propellant gas generator with eight pintle valves attached to the aft end. The pintle valve is designed to deliver variable upstream conditions to an attache( converging diverging nozzle. In order to investigate the full range of operating conditions 28 separate tests were conducted with varying pintle position and upstream pressure. Helium gas was utilized in order to closely mimic the speed of sound of the gas generator exhaust, minimizing required scaling during data analysis. The recordec pressure measurements were interrogated to multiple ends. The development of root mean square (RMS) value! versus Reynolds Number and Pintle position are important to creating bounding unsteady load curves for valve internal parts. Spectral analysis was also performed, helping to identify power spectral densities (PSD) of acoustic natural frequencies and boundary layer noise. An interesting and unexpected result was the identification of an acoustic mode within the valve which does not respond until the valve was over 60% open. Further, the response amplitude around this mode can be as large or larger than those associated with lower frequency modes.

Fishbach, Sean R.↗

NASA NTRS2010

Solar Flare Impulsive Phase Observations from SDO and Other Observatories

With the start of normal operations of the Solar Dynamics Observatory in May 2010, the Extreme ultraviolet Variability Experiment (EVE) and the Atmospheric Imaging Assembly (AIA) have been returning the most accurate solar XUV and EUV measurements every 10 and 12 seconds, respectively, at almost 100% duty cycle. The focus of the presentation will be the solar flare impulsive phase observations provided by EVE and AIA and what these observations can tell us about the evolution of the initial phase of solar flares. Also emphasized throughout is how simultaneous observations with other instruments, such as RHESSI, SOHO-CDS, and HINODE-EIS, will help provide a more complete characterization of the solar flares and the evolution and energetics during the impulsive phase. These co-temporal observations from the other solar instruments can provide information such as extending the high temperature range spectra and images beyond that provided by the EUV and XUV wavelengths, provide electron density input into the lower atmosphere at the footpoints, and provide plasma flows of chromospheric evaporation, among other characteristics.

Chamberlin, Phillip C.↗

NASA NTRS1996

Microstructural Development during Directional Solidification of Peritectic Alloys

A thorough understanding of the microstructures produced through solidification in peritectic systems has yet to be achieved, even though a large number of industrially and scientifically significant materials are in this class. One type of microstructure frequently observed during directional solidification consists of alternating layers of primary solid and peritectic solid oriented perpendicular to the growth direction. This layer formation is usually reported for alloy compositions within the two-phase region of the peritectic isotherm and for temperature gradient and growth rate conditions that result in a planar solid-liquid interface. Layered growth in peritectic alloys has not previously been characterized on a quantitative basis, nor has a mechanism for its formation been verified. The mechanisms that have been proposed for layer formation can be categorized as either extrinsic or intrinsic to the alloy system. The extrinsic mechanisms rely on externally induced perturbations to the system for layer formation, such as temperature oscillations, growth velocity variations, or vibrations. The intrinsic mechanisms approach layer formation as an alternative type of two phase growth that is inherent for certain peritectic systems and solidification conditions. Convective mixing of the liquid is an additional variable which can strongly influence the development and appearance of layers due to the requisite slow growth rate. The first quantitative description of layer formation is a model recently developed by Trivedi based on the intrinsic mechanism of cyclic accumulation and depiction of solute in the liquid ahead of the interface, linked to repeated nucleation events in the absence of convection. The objective of this research is to characterize the layered microstructures developed during ground-based experiments in which external influences have been minimized as much as possible and to compare these results to the current the model. Also, the differences between intrinsic and externally influenced layer formation were explored. The choice of alloy system is critical to a study of the formation of layered microstructures. The ideal system would have a well-characterized phase diagram, equal densities of both elements in the liquid state to minimize compositionally-driven convective flows, a low peritectic temperature to simplify directional solidification and the achievement of a high temperature gradient in the liquid, a broad composition range for the peritectic reaction, and a reasonable hardness at room temperature to facilitate handling and metallographic preparation. The In-Sn system was selected initially due to a very low peritectic temperature and the nearly equal densities of In and Sn in the liquid state. Since the In-rich peritectic reaction had apparently not been utilized previously for solidification research, experiments were conducted to check the phase diagram in the region of interest. The alloys in this system proved to be difficult to handle and prepare in bulk form with the equipment available, so experiments were initiated with the Sn-Cd system. Layered microstructures had been observed previously in Sn-Cd.

Lograsso, Thomas A.↗

NASA NTRS1996

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.↗

NASA NTRS2006

Global Aspects of Heliosphere-Geosphere Coupling

The magnetosphere serves both to hold off the solar wind and to couple it selectively to the Earth through the auroral zones. It has long been understood that the plasmasphere consists of geogenic plasmas that expand out of the sunlit low latitude ionosphere, but it was initially assumed that the hot low density plasma beyond the plasmasphere is largely of solar wind origin with a minor admixture of ionospheric plasmas, discovered via mass spectrometric observations in the early 70's. Since then, elaborate simulation models have been developed that have taken us from an era of cartoon physics to a new era of quantitative global comparisons between observations and theory. In most current global circulation models of the magnetosphere, the ionospheric load on the system is taken to lie exclusively in the thin F layer of the ionosphere. This layer is coupled with solar wind and magnetospheric plasmas via Maxwell stresses communicated by field aligned current systems, and with the thermosphere via ion-neutral charge exchange and Coulomb collisions. However, recent observations have shown us that ionospheric plasmas flow sporadically in various forms into the plasmasphere and up into the high latitude circulation cells of the outer magnetosphere, as driven by the solar wind and its variable intensity and magnetic field. Under some conditions, ionospheric material is observed to be the principal component of plasmas at the dayside magnetopause and in the plasma sheet and ring current regions. Given a global model of magnetospheric circulation, it is relatively straightforward to investigate the behaviors of ionospheric plasmas in response to solar wind drivers, and we report the results of such efforts here. We find that ionospheric plasmas dominate the plasma pressure in the magnetosphere in some regions, particularly when the solar wind is especially intense and-or southward directed. This result violates the assumption that the ionospheric load is confined to the F layer, and shows that the ionosphere is often an important dynamic element of the solar wind-ionosphere system throughout the magnetosphere, especially during larger geospace storms. This means that future global circulation models must account for ionospheric plasma inertia, heat capacity, and pressure, to be quantitatively credible.

Moore, Thomas E.↗

NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗

NASA NTRS2007

Piloted Ignition to Flaming in Smoldering Fire-Retarded Polyurethane Foam

Experimental results are presented on the piloted transition from smoldering to flaming in the fire-retarded polyurethane foam Pyrell . The samples are small rectangular blocks with a square cross section, vertically placed in the wall of a vertical wind tunnel. Three of the vertical sample sides are insulated and the fourth side is exposed to an upward oxidizer flow of variable oxygen concentration and to a variable radiant heat flux. The gases emitted from the smoldering reaction pass upwards through a pilot, which consists of a coiled resistance heating wire. In order to compensate for the solid-phase and gas-phase effects of the fire retardants on the piloted transition from smoldering to flaming in Pyrell, it was necessary to assist the process by increasing the power supplied to the smolder igniter and the pilot (compared to that used for non-fire retarded foam). The experiments indicate that the piloted transition from smoldering to flaming occurs when the gaseous mixture at the pilot passes the lean flammability limit. It was found that increasing the oxygen concentration or the external heat flux increases the likelihood of a piloted transition from smoldering to flaming, and generally decreases the time delay to transition. The piloted transition to flaming is observed in oxygen concentrations of 23% and above in both low-density and high-density Pyrell. Comparisons with previous experiments show that the piloted transition from smoldering to flaming is possible under a wider range of external conditions (i.e. lower oxygen concentration) than the spontaneous transition from smoldering to flaming. The results show that the fire retardants in Pyrell are very effective in preventing the piloted transition to flaming in normal air, but Pyrell is susceptible to smoldering and the piloted transition to flaming in oxygen-enriched environments. Therefore, precautions should be taken in the design of applications of Pyrell in oxygen-enriched environments to reduce to the risk of a piloted transition to flaming.

Putzeys, O.↗

DOE OSTI2026

Comprehensive Characterization of Water Group Ion Composition and Distributions in Saturn's Magnetosphere With Cassini Plasma Spectrometer Data

Saturn's magnetosphere is continuously supplied with neutrals from the Enceladus plume and the icy rings, which undergo ionization and charge-exchange to form a complex water-group plasma environment. While the Cassini Plasma Spectrometer (CAPS) instrument has provided extensive compositional information, detailed separation of individual water-group ion species in time-of-flight (TOF) data has not previously been achieved. In this study, we perform forward modeling of CAPS-IMS energy-per-charge (E/Q) and TOF spectra obtained between 2004 and 2012 to resolve O + , OH + , H 2 O + , and H 3 O + and to characterize their plasma properties, including number density, temperature, and thermodynamic κ. Our results demonstrate that O + is the dominant thermal ion species throughout Saturn's magnetosphere, comprising up to ∼70% of the total ion population beyond ∼5 Saturn radii (R S ). In contrast, molecular ions such as OH + , H 2 O + , and H 3 O + dominate closer to Enceladus but rapidly dissociate into atomic ions between ∼5 and 10 R S . This radial region is also characterized by the steepest increase in plasma flow speed, which rises from ∼40% to ∼80% of rigid corotation. Simultaneously, ion velocity distributions approach Maxwell–Boltzmann equilibrium, as indicated by high kappa values. These findings provide new constraints on the ion–neutral chemistry that regulates the balance between molecular and atomic ions in Saturn's magnetosphere. They also emphasize the critical role of the 5–10 R S region as a transition zone for both plasma composition and dynamics. Our results refine previous CAPS-based studies and underscore the need to incorporate seasonal variability and ionospheric coupling into future global models of Saturn's plasma environment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1999

Nonlinear Effects in Single-Pass ICRF Heating

The Variable Specific Impulse Magnetoplasma Rocket (VASIMR) concept employs Ion Cyclotron Resonant Frequency (ICRF) heating as the main power deposition mechanism. Since the ions accelerate to the full energy in a single pass through the cyclotron resonance, their response to the RF-field will be essentially nonlinear - hence the motivation to amend the commonly used linear approach to the problem. In a collisionless plasma, the energy gain of an accelerated ion is limited by the time the particle spends at the resonance. This time is affected by: (1) incident flow velocity, (2) longitudinal grad B force, (3) ambipolar electric field, and (4) ponderomotive force of the RF-field. Our analysis shows that the grad B force is the dominant factor at low to moderate levels of RF-power. We present nonlinear scaling for the energy gain and the absorption efficiency with RF-power and plasma parameters. We also demonstrate that the nonlinear regime exhibits a steep decrease in the plasma density at the resonance.

Arefiev, A. V.↗

DOE OSTI2024

Multizone Modeling of Black Hole Accretion and Feedback in 3D GRMHD: Bridging Vast Spatial and Temporal Scales

Simulating accretion and feedback from the horizon scale of supermassive black holes (SMBHs) out to galactic scales is challenging because of the vast range of scales involved. Elaborating on H. Cho et al., we describe and test a "multizone" technique, which is designed to tackle this difficult problem in three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. While short-timescale variability should be interpreted with caution, the method is demonstrated to be well-suited for finding dynamical steady states over a wide dynamic range. We simulate accretion on a nonspinning SMBH ($a\ast$ = 0) using initial conditions and the external galactic potential from a large-scale galaxy simulation and achieve a steady state over eight decades in radius. As found in H. Cho et al., the density scales with radius as ρ ∝ r –1 inside the Bondi radius R B , which is located at R B = 2 × 10 5 r g (≈60 pc for M87), where r g is the gravitational radius of the SMBH; the plasma-β is ~ unity, indicating an extended magnetically arrested state; the mass accretion rate $\dot{M}$ is ≈1% of the analytical Bondi accretion rate ${\dot{M}}_{{\rm{B}}};$ and there is continuous energy feedback out to ≈100R B (or beyond > kpc) at a rate $\approx 0.02\dot{M}{c}^{2}$. Surprisingly, no ordered rotation in the external medium survives as the magnetized gas flows to smaller radii, and the final steady solution is very similar to when the exterior has no rotation. Using the multizone method, we simulate GRMHD accretion over a wide range of Bondi radii, R B ~ 10 2 –10 7 r g , and find that $\dot{M}/{\dot{M}}_{{B}}\approx {({R}_{{B}}/6\,{r}_{g})}^{-0.5}$.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2011

Eucrite Impact Melt NWA 5218 - Evidence for a Large Crater on Vesta

Northwest Africa (NWA) 5218 is a 76 g achondrite that is classified as a eucrite [1]. However, an initial classification [2] describes it as a "eucrite shock-melt breccia...(in which) large, partially melted cumulate basalt clasts are set in a shock melt flow...". We explore the petrology of this clast-bearing impact melt rock (Fig. 1), which could be a characteristic lithology at large impact craters on asteroid Vesta [3]. Methods: Optical microscopy, scanning electronmicroscopy, and Raman spectroscopy were used on a thin section (Fig. 1) for petrographic characterization. The impact melt composition was determined by 20 m diameter defocused-beam analyses with a Cameca SX-100 electron microprobe. The data from 97 spots were corrected for mineral density effects [4]. Constituent mineral phases were analyzed with a focusedbeam. Bidirectonal visible and near-infrared (VNIR) and biconical FT-IR reflectance spectra were measured on the surface of a sample slab on its central melt area and on an eucrite clast, and from 125-500 m and <125 m powders of melt. Results: General petrography: The sample specimen is a coherent, medium dark-grey (N4), melt rock. The thin section captures a central, subophitic-textured melt that contains ~1 cm to tens of m-size subangular to rounded, variably-shocked eucrite clasts. Clasts >100 m are coarse-grained with equigranular ~1 mm size plagioclase, quartz, and clinopyroxene (Fig. 1). Single crystals of chromite, ilmenite, zircon, Ca-Mg phosphate, Fe-metal, and troilite are embedded in the melt. Polymineralic clasts are mostly compositionally similar to the above mentioned larger clasts but scarce granulitic fragments are observed as well.

Wittmann, Axel↗