Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density”

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,567 records · Page 87

Modeling Variable-Impedance, Magnetically Insulated, Transmission Lines

Coupling generator current to electron-beam loads and dynamic-inductance loads (z pinches) is greatly improved with low-inductance vacuum transmission lines. Typically, magnetically insulated transmission lines (MITLs) utilize nearly constant impedance geometries where the minimum impedance (anode-cathode gap) is limited by electron loss during the setup of magnetic insulation. Constant-impedance MITL designs are not optimized for low inductance (Zτ) when given a constraint on vacuum electron flow. We describe the design and the modeling of MITLs that have a nonconstant impedance. The baseline variable-impedance MITL designs are based on equilibrium vacuum-electron flow models [1] that are incorporated into the Screamer circuit code [2]. Such designs satisfy the MITL limits in the areal density of the loss current and have magnetically insulated flow that is stable everywhere along the MITL. Different MITL designs are then modeled analytically and with a 2-D particle-in-cell computer code in order to quantify subtle electron losses driven by changes to the vacuum impedance as a function of distance along the transmission line. We present a variable-impedance MITL design that delivers more current with a lower inductance than found in constant-impedance MITL designs.

47 OTHER INSTRUMENTATION↗

Experimental Electron Density Distribution in Actinide Compounds – an Experimental Atoms in Molecules (AIM) Approach

Understanding the basic chemistry of the actinides, in particular of uranium, is essential to address issues related to nuclear fuels, nuclear waste and environmental concerns. Atoms in general are chemically bound to their neighbors via the distribution of electrons within the intervening space. However, this only represents a small proportion of the total electrons, the others are tightly bound to the atomic nuclei and do not participate in chemical bonding. The distribution of electrons in crystalline materials may be examined using the scattering of X-rays. As all electrons scatter X-rays, obtaining the information about the bonding electrons in the presence of a large excess of “irrelevant” electrons requires extremely accurate scattering measurements. Indeed, until recently this was considered impossible for elements as heavy as uranium. In our current project we have made significant improvements in the data quality obtainable for uranium containing compounds. In order to determine the electron density distribution we fit a mathematical model describing the scattering effect to the observed scattering data. This mathematical model has a number of variable parameters that describe the radial and angular distribution of electrons around an atomic center. We have developed such a model to be particularly useful for elements as heavy as uranium. Given the best fit parameters from this model, we are able to map the electron distribution around all atoms in a chemical compound. The topology of this distribution provides us with important characteristics with regard to the chemical bonds. We have now applied these methods to several uranium compounds. In particular we have studied the bonds made by uranium to oxygen, fluorine and chlorine in a number of different environments. The experimental protocols and data analysis protocols developed during this project will permit others in the field to explore the chemical bonding of other heavy element systems. From the few compounds studied to date, we are already discovering the parameters that are characteristic of bonding in such systems, and providing essential experimental benchmarks for theoreticians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Variation in supercurrent connectivity and vortex pinning of state-of-the-art pulsed laser deposited RE Ba 2 Cu 3 O 7- x coated conductors

REBa 2 Cu 3 O 7-x (REBCO, RE: rare earth) coated conductors (CCs) suffer from great critical current I c differences between different manufacturers, I c variations within individual manufacturers, and often significant lengthwise fluctuations. The understanding of such variations is complicated by the lack of a direct correlation between I c and the critical current density J c . In fact, although J c is the fundamental property determined by the local vortex pinning landscape, I c is often limited by variable current blocking mechanisms. An important practical complexity is commercial practices of performing J c and I c evaluations based on I c at 77 K and self-field (sf), where connectivity variation dominates over vortex pinning variations. However, at higher fields and lower temperatures, vortex pinning becomes more complex and highly variable, making predictions of I c and J c at arbitrary temperature T, magnetic field H, and field orientation θ, quite uncertain. To address some aspects of this problem, we conducted detailed spool-to-spool performance characterization on recently manufactured REBCO CCs. Despite J c (77 K, sf) varies by only ∼11%, J c (77 K, 1 T) of its minimum and maximum (for H//ab-plane) show variations of ∼21% and ∼32%, respectively. This emphasizes the shortcoming in using J c (77 K, sf) as parameter for evaluating even the low field performance. An even more remarkable spool-to-spool J c variation of ∼68% was observed at 20 K and 15 T for H//ab-plane. To identify the origin of such lack of reproducibility we performed microstructural characterizations, which revealed, within the REBCO layer, large variation in the density of copper oxide (CuO x ) particles ranging from 0.1 to 2 μm in size. We believe that they play a decisive role in reducing the effective cross-section of the REBCO layer by not simply blocking current themselves, but also by nucleating off-axis REBCO grains, whose misoriented grain boundaries adversely impact REBCO grain-to-grain connectivity. The REBCO growth associated with high density CuO x particles also leads to the more disordered spatial arrangement of BaHfO 3 precipitate arrays, which, when self-aligning along the ab-planes, generate stronger pinning enhancing J c (H//ab) at all temperatures. In this way, we established that variations of both connectivity and vortex pinning are thus directly coupled. Our results also explain why the so-called ‘lift-factor’, typically defined by the ratio between I c (T, H) and I c (77 K, sf), frequently turns out to be unreliable.

36 MATERIALS SCIENCE↗

Hard X-Ray Variability of the Black Hole Candidate GRO J0422+32 During its 1992 Outburst

We have studied the hard X-ray variability of the soft X-ray transient GRO J0422+32 with BATSE (Burst and Transient Source Experiment) in the 20-100 keV energy band. Our analysis covers 180 days following the first X-ray detection of the source on August 5, 1992, fully covering its primary and secondary X-ray outbursts. We compute power density spectra (PDSs) in the 20-50, 50-100, and 20-100 keV energy bands, in the frequency interval 0.002-0.488 Hz. The PDSs of GRO J0422+32 are approximately flat up to a break frequency, and decay as a power law above it, with index approximately 1. During the first 70 days of the X-ray outburst, the PDSs of GRO J0422+32 show a significant quasi-periodic oscillation (QPO) peak near approximately 0.2 Hz, superposed on the power-law tail. The break frequency of the PDSs obtained during the primary X-ray outburst of GRO J0422+32 occurs at 0.041 +/- 0.006 Hz; during the secondary outburst, the break is at 0.081 +/- 0.015 Hz. The power density at the break ranged between 44% and 89% / square root of Hz (20-100 keV). The canonical anticorrelation between the break frequency and the power density at the break, observed in Cyg X-1 and other BHCs in the low state, is not observed in the PDSs of GRO J0422+32. We compare our results with those of similar variability studies of Cyg X-1. The relation between the spectral slope and the amplitude of the X-ray variations of GRO J0422+32 is similar to that of Cyg X-1; however, the relation between the hard X-ray flux and the amplitude of its variation is opposite to what has been found in Cyg X-1. Phase lags between the X-ray flux variations of GRO J0422+32 at high and low photon energies could only be derived during the first 30 days of its outburst. During this period, the variations in the 50-100 keV band lag those in the 20-50 keV energy band by an approximately constant phase difference of 0.039(3) rad in the frequency interval 0.02-0.20 Hz. The time lags of GRO J0422+32 during the first 30 days of the outburst decrease with frequency as a power law with index 0.9 for v > 0.01 Hz.

vanderHooft, F.↗

Linear Temporal Stability Analysis of a Low-Density Round Gas Jet Injected into a High-Density Gas

It has been observed in previous experimental studies that round helium jets injected into air display a repetitive structure for a long distance, somewhat similar to the buoyancy-induced flickering observed in diffusion flames. In order to investigate the influence of gravity on the near-injector development of the flow, a linear temporal stability analysis of a round helium jet injected into air was performed. The flow was assumed to be isothermal and locally parallel; viscous and diffusive effects were ignored. The variables were represented as the sum of the mean value and a normal-mode small disturbance. An ordinary differential equation governing the amplitude of the pressure disturbance was derived. The velocity and density profiles in the shear layer, and the Froude number (signifying the effects of gravity) were the three important parameters in this equation. Together with the boundary conditions, an eigenvalue problem was formulated. Assuming that the velocity and density profiles in the shear layer to be represented by hyperbolic tangent functions, the eigenvalue problem was solved for various values of Froude number. The temporal growth rates and the phase velocity of the disturbances were obtained. The temporal growth rates of the disturbances increased as the Froude number was reduced (i.e. gravitational effects increased), indicating the destabilizing role played by gravity.

Lawson, Anthony L.↗

Noninteracting electrons in a prototypical one-dimensional sinusoidal potential

A prototypical model of a one-dimensional metallic monatomic solid containing noninteracting electrons is studied, where the argument of the cosine potential energy, periodic with the lattice, contains the first reciprocal lattice vector G 1 =2π/a, where a is the lattice constant. The time-independent Schrödinger equation can be written in reduced variables as a Mathieu equation for which numerically exact solutions for the band structure and wave functions are obtained. The band structure has band gaps that increase with increasing amplitude q of the cosine potential. In the extended-zone scheme, the energy gaps decrease with increasing index n of the Brillouin-zone boundary ka=nπ, where k is the crystal momentum of the electron. The wave functions of the band electron are derived for various combinations of k and q as complex combinations of the real Mathieu functions with even and odd parity, and the normalization factor is discussed. The wave functions at the bottoms and tops of the bands are found to be real or imaginary, respectively, corresponding to standing waves at these energies. Irrespective of the wave vector k within the first Brillouin zone, the electron probability density is found to be periodic with the lattice. The Fourier components of the wave functions are derived versus q, which reveal multiple reciprocal-lattice-vector components with variable amplitudes in the wave functions unless q = 0. The magnitudes of the Fourier components are found to decrease exponentially as a power of n for n~3 to 45 for ka = π/2 and q = 2, and a precise fit is obtained to the data. The probability densities and probability currents obtained from the wave functions are also discussed. The probability currents are found to be zero for crystal momenta at the tops and bottoms of the energy bands, because the wave functions for these crystal momenta are standing waves. To conlcude, the band structure is calculated from the central equation and compared to the numerically exact band structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hard X-Ray Variability of the Black Hole Candidate GRO J0422+32 During its 1992 Outburst

We have studied the hard X-ray variability of the soft X-ray transient GRO J0422+32 with the Burst and Transient Source Experiment (BATSE) in the 20-100 keV energy band. Our analysis covers 180 days following the first X-ray detection of the source on 1992 August 3, fully covering its primary and secondary X-ray outburst. We computed power density spectra (PDSS) in the 20-50, 50-100, and 20-100 keV energy bands., in the frequency interval 0.002-0.488 Hz. The PDSs of GRO J0422+32 are approximately flat up to a break frequency and decay as a power law above, with index about 1. During the first 70 days of the X-ray outburst, the PDSs of GRO J0422+32 show a significant QPO Peak near about 0.2 Hz, superposed on the power-law tail. The break frequency of the PDSs obtained during the primary X-ray outburst of GRO J0422+32 occurs at 0.041+/-0.006 Hz; during the secondary outburst the break is at 0.081q0.01.5 Hz. The power density at the break ranged between 44 and 89% Hz(exp -1/2) 20-100 keV). The canonical anticorrelation between the break frequency and the power density at the break, observed in Cyg X-1 and other BHCs in the low state, is not observed in the PDSs of GRO J0422+32. We compare our results with those of similar variability studies of Cyg X-1. The relation between the spectral slope and the amplitude of the X-ray variations of GRO J0422+32 is similar to that of Cyg X-1.

VanderHooft, F.↗

Understanding Temporal and Spatial Variability of the Lunar Helium Atmosphere Using Simultaneous Observations from LRO, LADEE, and ARTEMIS

Simultaneous measurements of helium in the exosphere of the Moon are made from the Lunar Reconnaissance Orbiter (LRO) Lyman Alpha Mapping Project (LAMP) and the Lunar Atmosphere and Dust Environment Explorer (LADEE) Neutral Mass Spectrometer (NMS) through the entire 5-month span of the LADEE mission. In addition, the ARTEMIS mission monitored the solar wind alpha particle flux to the Moon. Modeling the lunar helium exosphere, we relate the LAMP polar observations to the LADEE equatorial observations. Further, using the ARTEMIS alpha flux in the Monte Carlo model reproduces the temporal variations in helium density. Comparing the LAMP data to the LADEE data shows excellent agreement. Comparing those with the ARTEMIS data reveals that the solar wind alpha flux is the primary driver to variability in the helium exosphere throughout the LADEE mission. Using a decay time for exospheric helium of 5 days, we determine that the solar wind contributes 64 +/- 5% of the helium to the lunar exosphere. The remaining 36 +/- 5% is presumed to come from outgassing of radiogenic helium from the interior of the Moon. Furthermore, the model reproduces the measurements if 63 +/- 6% of the incident alpha particles are converted to thermalized helium atoms through the interaction between the alphas and the lunar surface. However, these values are dependent on both inferred source rates from LAMP and LADEE observations and on the assumed time constant of the exospheric decay rate.

Moon↗

The Multi-Layer Variable Absorbers in NGC 1365 Revealed by XMM-Newton and NuSTAR

Between 2012 July and 2013 February, NuSTAR and XMM-Newton performed four long-look joint observations of the type 1.8 Seyfert, NGC 1365. We have analyzed the variable absorption seen in these observations in order to characterize the geometry of the absorbing material. Two of the observations caught NGC 1365 in an unusually low absorption state, revealing complexity in the multi-layer absorber that had previously been hidden. We find the need for three distinct zones of neutral absorption in addition to the two zones of ionized absorption and the Compton-thick torus previously seen in this source. The most prominent absorber is likely associated with broad-line region clouds with column densities of around approximately 10 (sup 23) per square centimeter and a highly clumpy nature as evidenced by an occultation event in 2013 February. We also find evidence of a patchy absorber with a variable column around approximately 10 (sup 22) per square centimeter and a line-of-sight covering fraction of 0.3-0.9, which responds directly to the intrinsic source flux, possibly due to a wind geometry. A full-covering, constant absorber with a low column density of approximately 1 by 10 (sup 22) per square centimeter is also present, though the location of this low density haze is unknown.

galaxies: active aEuro" galaxies: individual (NGC ↗

Seasonal variability of wake impacts on US mid-Atlantic offshore wind plant power production

Abstract. The mid-Atlantic will experience rapid wind plant development due to its promising wind resource located near large population centers. Wind turbines and wind plants create wakes, or regions of reduced wind speed, that may negatively affect downwind turbines and plants. We evaluate wake variability and annual energy production with the first yearlong modeling assessment using the Weather Research and Forecasting model, deploying 12 MW turbines across the domain at a density of 3.14 MW km−2, matching the planned density of 3 MW km−2. Using a series of simulations with no wind plants, one wind plant, and complete build-out of lease areas, we calculate wake effects and distinguish the effect of wakes generated internally within one plant from those generated externally between plants. We also provide a first step towards uncertainty quantification by testing the amount of added turbulence kinetic energy (TKE) by 0 % and 100 %. We provide a sensitivity analysis by additionally comparing 25 % and 50 % for a short case study period. The strongest wakes, propagating 55 km, occur in summertime stable stratification, just when New England's grid demand peaks in summer. The seasonal variability of wakes in this offshore region is much stronger than the diurnal variability of wakes. Overall, yearlong simulated wake impacts reduce power output by a range between 38.2 % and 34.1 % (for 0 %–100 % added TKE). Internal wakes cause greater yearlong power losses, from 29.2 % to 25.7 %, compared to external wakes, from 14.7 % to 13.4 %. The overall impact is different from the linear sum of internal wakes and external wakes due to non-linear processes. Additional simulations quantify wake uncertainty by modifying the added amount of turbulent kinetic energy from wind turbines, introducing power output variability of 3.8 %. Finally, we compare annual energy production to New England grid demand and find that the lease areas can supply 58.8 % to 61.2 % of annual load. We note that the results of this assessment are not intended to make nor are they suitable to make commercial judgments about specific wind projects.

17 WIND ENERGY↗

Investigation of a Light Gas Helicon Plasma Source for the VASIMR Space Propulsion System

An efficient plasma source producing a high-density (approx.10(exp 19/cu m) light gas (e.g. H, D, or He) flowing plasma with a high degree of ionization is a critical component of the Variable Specific Impulse Magnetoplasma Rocket (VASIMR) concept. We are developing an antenna to apply ICRF power near the fundamental ion cyclotron resonance to further accelerate the plasma ions to velocities appropriate for space propulsion applications. The high degree of ionization and a low vacuum background pressure are important to eliminate the problem of radial losses due to charge exchange. We have performed parametric (e.g. gas flow, power (0.5 - 3 kW), magnetic field , frequency (25 and 50 MHz)) studies of a helicon operating with gas (H2 D2, He, N2 and Ar) injected at one end with a high magnetic mirror downstream of the antenna. We have explored operation with a cusp and a mirror field upstream. Plasma flows into a low background vacuum (<10(exp -4) torr) at velocities higher than the ion sound speed. High densities (approx. 10(exp 19/cu m) have been achieved at the location where ICRF will be applied, just downstream of the magnetic mirror.

Squire, J. P.↗

How Far Is Quasar UV/Optical Variability from a Damped Random Walk at Low Frequency?

Studies have shown that UV/optical light curves of quasars can be described using the prevalent damped random walk (DRW) model, also known as the Ornstein–Uhlenbeck process. A white noise power spectral density (PSD) is expected at low frequency in this model; however, a direct observational constraint to the low-frequency PSD slope is difficult due to the limited lengths of the light curves available. Meanwhile, quasars show scatter in their DRW parameters that is too large to be attributed to uncertainties in the measurements and dependence on the variation of known physical factors. In this work we present simulations showing that, if the low-frequency PSD deviates from the DRW, the red noise leakage can naturally produce large scatter in the variation parameters measured from simulated light curves. The steeper the low-frequency PSD slope, the larger scatter we expect. Based on observations of SDSS Stripe 82 quasars, we find that the low-frequency PSD slope should be no steeper than −1.3. The actual slope could be flatter, which consequently requires that the quasar variabilities should be influenced by other unknown factors. We speculate that the magnetic field and/or metallicity could be such additional factors.

79 ASTRONOMY AND ASTROPHYSICS↗

A combined machine learning and density functional theory study of binary Ti-Nb and Ti-Zr alloys: Stability and Young’s modulus

The multicomponent Ti alloys, specifically the -phase, have experienced a strong growth over the last decades, due to their outstanding properties of ultra-high strength and low Young’s modulus. These properties play a significant role in many aerospace and biomedical applications. Selection and optimization of multicomponent alloys is challenging due to the vast chemical and compositional space. Here we investigate the use of machine learning techniques informed by density functional calculations to guide the selection of Nb- and Zr-based Ti binary alloys. From the cubic structures obtained from high throughput calculations and literature, we identify several structures with Young’s moduli below 40 GPa. The multivariant decision tree methods provide efficient surrogate models to identify structure variables have high influences on the energetic stability and Young’s modulus. Furthermore, we implement a workflow of incorporating DFT provided results and machine learning method to explore the chemical and composition space of other binary and multicomponent alloys, to eventually accelerate the material design via taking advantages of identified key variables.

36 MATERIALS SCIENCE↗

High temperature gasification of high heating-rate chars using a flat-flame reactor

The increasing interest in gasification and oxy-fuel combustion of biomass has heightened the need for a detailed understanding of char gasification in industrially relevant environments (i.e., high temperature and high-heating rate). Despite innumerable studies previously conducted on gasification of biomass, very few have focused on such conditions. Consequently, in this study the high-temperature gasification behaviors of biomass-derived chars were investigated using non-intrusive techniques. Two biomass chars produced at a heating rate of approximately 10 4 K/s were subjected to two gasification environments and one oxidation environment in an entrained flow reactor equipped with an optical particle-sizing pyrometer. A coal char produced from a common U.S. low sulfur subbituminous coal was also studied for comparison. Both char and surrounding gas temperatures were precisely measured along the centerline of the furnace. Despite differences in the physical and chemical properties of the biomass chars, they exhibited rather similar reaction temperatures under all investigated conditions. On the other hand, a slightly lower particle temperature was observed in the case of coal char gasification, suggesting a higher gasification reactivity for the coal char. A comprehensive numerical model was applied to aid the understanding of the conversion of the investigated chars under gasification atmospheres. In addition, a sensitivity analysis was performed on the influence of four parameters (gas temperature, char diameter, char density, and steam concentration) on the carbon conversion rate. Here, the results demonstrate that the gas temperature is the most important single variable influencing the gasification rate.

09 BIOMASS FUELS↗

Atomistically-informed modeling of point defect clustering and evolution in irradiated ThO 2

A cluster dynamics (CD) model has been developed to investigate the nucleation and growth of point defect clusters, i.e., interstitial prismatic loops and nanoscale and sub-nanoscale voids, in ThO 2 during irradiation by energetic particles. The model considers cluster off-stoichiometry due to the asymmetry of point defect generation on the O and Th sublattices under irradiation, as well as the point defect diffusivities and the defect binding energies to clusters. The energies were established using detailed molecular dynamics simulations considering the statistical variability of cluster configuration. A high-order adaptive time-integration has been used to solve the model. The predicted loop density and their average size is in good agreement with reported experimental observations for proton irradiated ThO 2 at 600°C. In conclusion, the model did not predict void evolution due to the sluggish kinetics of cation vacancies, explaining the absence of voids in proton irradiated ThO 2 (and other oxides) at relatively low temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The relative influence of climate and housing development on current and projected future fire patterns and structure loss across three California landscapes

Climate and land use patterns are expected to change dramatically in the coming century, raising concern about their effects on wildfire patterns and subsequent impacts to human communities. The relative influence of climate versus land use on fires and their impacts, however, remains unclear, particularly given the substantial geographical variability in fire-prone places like California. Here, we developed a modeling framework to compare the importance of climatic and human variables for explaining fire patterns and structure loss for three diverse California landscapes, then projected future large fire and structure loss probability under two different climate (hot-dry or warm-wet) and two different land use (rural or urban residential growth) scenarios. The relative importance of climate and housing pattern varied across regions and according to fire size or whether the model was for large fires or structure loss. The differing strengths of these relationships, in addition to differences in the nature and magnitude of projected climate or land use change, dictated the extent to which large fires or structure loss were projected to change in the future. Despite this variability, housing and human infrastructure were consistently more responsible for explaining fire ignitions and structure loss probability, whereas climate, topography, and fuel variables were more important for explaining large fire patterns. For all study areas, most structure loss occurred in areas with low housing density (from 0.08 to 2.01 units/ha), and expansion of rural residential land use increased structure loss probability in the future. Regardless of future climate scenario, large fire probability was only projected to increase in the northern and interior parts of the state, whereas climate change had no projected impact on fire probability in southern California. Given the variation in fire-climate relationships and land use effects, policy and management decision-making should be customized for specific geographical regions.

54 ENVIRONMENTAL SCIENCES↗

Drivers of Decadal Carbon Fluxes Across Temperate Ecosystems

Long-running eddy covariance flux towers provide insights into how the terrestrial carbon cycle operates over multiple timescales. Here, we evaluated variation in net ecosystem exchange (NEE) of carbon dioxide (CO 2 ) across the Chequamegon Ecosystem-Atmosphere Study AmeriFlux core site cluster in the upper Great Lakes region of the USA from 1997 to 2020. The tower network included two mature hardwood forests with differing management regimes (US-WCr and US-Syv), two fen wetlands with varying levels of canopy sheltering and vegetation (US-Los and US-ALQ), and a very tall (400 m) landscape-level tower (US-PFa). Together, they provided over 70 site-years of observations. The 19-tower Chequamegon Heterogenous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 campaign centered around US-PFa provided additional information on the spatial variation of NEE. Decadal variability was present in all long-term sites, but cross-site coherence in interannual NEE in the earlier part of the record became weaker with time as non-climatic factors such as local disturbances likely dominated flux time series. Average decadal NEE at the tall tower transitioned from carbon source to sink to near neutral over 24 years. Respiration had a greater effect than photosynthesis on driving variations in NEE at all sites. Declining snowfall offset potential increases in assimilation from warmer springs, as less-insulated soils delayed start of spring green-up. Higher CO 2 increased maximum net assimilation parameters but not total gross primary productivity. Stand-scale sites were larger net sinks than the landscape tower. Clustered, long-term carbon flux observations provide value for understanding the diverse links between carbon and climate and the challenges of upscaling these responses across space.

54 ENVIRONMENTAL SCIENCES↗

Extended MHD modeling of tearing-driven magnetic relaxation

Discrete relaxation events in reversed-field pinch relevant configurations are investigated numerically with nonlinear extended magnetohydrodynamic modeling, including the Hall term in Ohm’s law and first-order ion finite Larmor radius effects. Our results show variability among relaxation events, where the Hall dynamo effect may help or impede the MHD dynamo effect in relaxing the parallel current density profile. The competitive behavior arises from multi-helicity conditions where the dominant magnetic fluctuation is relatively small. The resulting changes in parallel current density and parallel flow are aligned in the core, consistent with experimental observations. Analysis of simulation results also confirms that force density from fluctuation-induced Reynolds stress arises subsequent to the drive from fluctuation-induced Lorentz force density. Transport of momentum density is found to be dominated by the fluctuation-induced Maxwell stress over most of the cross section with viscous and gyroviscous contributions being large in the edge region. The findings resolve a discrepancy with respect to the relative orientation of current density and flow relaxation, which had not been realized or investigated in Ref. [King et. al. Phys. Plasmas 19, 055905 (2012)], where only the magnitude of flow relaxation is actually consistent with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗