SEARCH · Engineering Papers

Results for “Calculational Margin”

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 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2026

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

California Price Response Potential Study

California's energy landscape is undergoing a significant transformation, driven by the increasing integration of renewable energy sources, the increased adoption of distributed energy resources, the electrification of end-use loads, and the growing need for grid efficiency. To address these challenges, recent revisions to the State’s Load Management Standards (LMS) require all of California’s large utilities and community choice aggregators (CCAs) to offer dynamic electricity pricing options to customers by 2027. Dynamic pricing, which involves varying electricity rates based on real-time supply and demand conditions, offers a promising solution for optimizing grid operations, reducing costs, and incentivizing efficient use of grid capacity. Effective implementation of dynamic pricing requires understanding the potential impacts on customer bills, system load, and the cost-effectiveness of automation technologies. This study aims to evaluate the load response of various end-use devices to hourly dynamic prices. The end-uses studied here are space cooling, space heating, water heating, crop irrigation, pool and spa pumps, and electric vehicle (EV) charging, all for both residential and commercial applications, except for crop irrigation. In 2030, these end uses are forecasted to account for 18% of annual electricity demand in the state, but 40% of demand in the peak net load hour. By modeling possible price-responsive load dispatch algorithms and assessing the resulting impacts on both individual bills and the overall grid, we seek to inform policymakers and utilities about the potential benefits and challenges associated with dynamic pricing, and considerations for the design of dynamic pricing tariffs. Additionally, we will explore the cost effectiveness of adopting automation technologies to enable devices to respond more effectively to real-time price signals. This study considers a range of price profiles, accounting for differences across utilities and customer classes, and presents scenarios for dynamic price design via variation in the percentage of total customer electric costs that are allocated dynamically (versus constituting a fixed portion of the hourly volumetric price). We present results focused primarily on 2030, forecasting electricity prices under both low and high-cost scenarios, to inform longer-term tariff design considerations. We design tariffs by starting with 2019 prices that were calculated according to CalFUSE guidance (CPUC, 2022) and that have been used in recent studies; these prices are all-in volumetric rates that vary by utility and are revenue-neutral to each customer class. They are developed by considering six electricity cost components that are allocated hourly based on system load indicators (gross and net load, and wholesale prices). These prices are forecasted to 2030 for low and high cost scenarios, considering recent trends in total electricity costs with and without years of substantial wildfire mitigation investments. These tariffs, which allocate all costs on an hourly basis, are considered our “Full” dynamic tariff design scenario, while two additional scenarios explore allocating a portion of costs as a flat volumetric charge: the “Medium” scenario allocates 50% of revenue dynamically (and keeps 50% flat), while the “Mild” scenario allocates 20% of revenue dynamically. The 20% dynamic allocation on the Mild scenario aims to represent a case where only the marginal operating costs of the grid are included in the dynamic price.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2021

Origin and resummation of threshold logarithms in the lattice QCD calculations of PDFs

Many present lattice QCD approaches to calculate the parton distribution functions (PDFs) rely on a factorization formula or effective theory expansion of certain Euclidean matrix elements in boosted hadron states. In the quasi- and pseudo-PDF methods, the matching coefficient in the factorization or expansion formula includes large logarithms near the threshold, which arise from the subtle interplay of collinear and soft divergences of an underlying 3D momentum distribution. We use the standard prescription to resum such logarithms in the Mellin-moment space at next-to-leading logarithmic accuracy, which also accounts for the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution, and we show that it can suppress the PDF at large x . Unlike the deep inelastic scattering and Drell-Yan cross sections, the resummation formula is away from the Landau pole. We then apply our formulation to reanalyze the recent lattice results for the pion valence PDF, and find that within the current data sensitivity, the effect of threshold resummation is marginal for the accessible moments and the PDF at large x.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

The focus of this work was to characterize the fundamental flow physics and the overall performance effects due to increased rotor tip clearance heights in axial compressors. Data have been collected in the three-stage axial research compressor at Purdue University with a specific focus on analyzing the multistage effects resulting from the tip leakage flow. Three separate rotor tip clearance heights were studied with nominal tip clearance heights of 1.5%, 3.0%, and 4.0% based on a constant annulus height. Overall compressor performance was investigated at four corrected speedlines (100%, 90%, 80%, and 68%) for each of the three tip clearance configurations using total pressure and total temperature rakes distributed throughout the compressor. The results have confirmed results from previous authors showing a decrease of total pressure rise, isentropic efficiency, and stall margin which is approximately linear with increasing tip clearance height. The stall inception mechanisms have also been evaluated at the same corrected speeds for each of the tip clearance configurations. Detailed flow field measurements have been collected at two loading conditions, nominal loading (NL) and high loading (HL), on the 100% corrected speedline for the smallest and largest tip clearance heights (1.5% and 4.0%). Steady detailed radial traverses of total pressure at the exit of each stator row have been supported by flow visualization techniques to identify regions of flow recirculation and separation. Furthermore, detailed radial traverses of time-resolved total pressures at the exit of each rotor row have been measured with a fast-response pressure probe. These data have helped to quantify the size of the leakage flow at the exit of each rotor. Thermal anemometry has also been implemented to evaluate the time-resolved three-dimensional components of velocity throughout the compressor and calculate blockage due to the rotor tip leakage flow throughout the compressor. These measurements have also been used to calculate streamwise vorticity. Time-resolved static pressure measurements have been collected over the rotor tips for all rotors with each of the three tip clearance configurations for up to five loading conditions along the 100% corrected speedline using fast-response piezoresistive pressure sensors. These time-resolved static pressure measurements, as well as the time-resolved total pressures and velocities have helped to reveal a profound influence of the upstream stator vane on the size and shape of the rotor tip leakage flow. Finally, a novel particle image velocimetry (PIV) technique has been developed as a proof-of- concept. In contrast to PIV methods that have been typically been utilized for turbomachinery applications in the past, the method used for this study introduced the laser light through the same access window that was also used to image the flow. This new method addresses potential concerns related to the intrusive laser-introducing techniques that have typically been utilized by other authors in the past. Ultimately, the data collected for this project represent a unique data set which contributes to build a better understanding of the tip leakage flow field and its associated loss mechanisms. These data will facilitate future engine design goals leading to small blade heights in the rear stages of high pressure compressors and aid in the development of new blade designs which are desensitized to the performance penalties attributed to rotor tip leakage flows.

turbomachinery↗

NASA NTRS2015

An Experimental Investigation of the Flow Physics Associated With End Wall Losses and Large Rotor Tip Clearances as Found in the Rear Stages of a High Pressure Compressor

Berdanier, Reid A.↗

DOE OSTI2020

TDCOSMO - IV. Hierarchical time-delay cosmography – joint inference of the Hubble constant and galaxy density profiles

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data.Key words: gravitational lensing: strong / galaxies: general / galaxies: kinematics and dynamics / distance scale / cosmological parameters / cosmology: observations⋆ The full analysis is available at https://github.com/TDCOSMO/hierarchy_analysis_2020_public.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2013

Constraining Landscape History and Glacial Erosivity Using Paired Cosmogenic Nuclides in Upernavik, Northwest Greenland

High-latitude landscape evolution processes have the potential to preserve old, relict surfaces through burial by cold-based, nonerosive glacial ice. To investigate landscape history and age in the high Arctic, we analyzed in situ cosmogenic Be(sup 10) and Al (sup 26) in 33 rocks from Upernavik, northwest Greenland. We sampled adjacent bedrock-boulder pairs along a 100 km transect at elevations up to 1000 m above sea level. Bedrock samples gave significantly older apparent exposure ages than corresponding boulder samples, and minimum limiting ages increased with elevation. Two-isotope calculations Al(sup26)/B(sup 10) on 20 of the 33 samples yielded minimum limiting exposure durations up to 112 k.y., minimum limiting burial durations up to 900 k.y., and minimum limiting total histories up to 990 k.y. The prevalence of BE(sup 10) and Al(sup 26) inherited from previous periods of exposure, especially in bedrock samples at high elevation, indicates that these areas record long and complex surface exposure histories, including significant periods of burial with little subglacial erosion. The long total histories suggest that these high elevation surfaces were largely preserved beneath cold-based, nonerosive ice or snowfields for at least the latter half of the Quaternary. Because of high concentrations of inherited nuclides, only the six youngest boulder samples appear to record the timing of ice retreat. These six samples suggest deglaciation of the Upernavik coast at 11.3 +/- 0.5 ka (average +/- 1 standard deviation). There is no difference in deglaciation age along the 100 km sample transect, indicating that the ice-marginal position retreated rapidly at rates of approx.120 m yr(sup−1).

glacial↗

DOE OSTI2022

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Comparative Analysis of Standard and Advanced USL Methodologies for Nuclear Criticality Safety

The American National Standards Institute/American Nuclear Society national standards 8.1 and 8.24 provide guidance on the requirements and recommendations for establishing confidence in the results of the computerized models used to support operation with fissionable materials. By design, the guidance is not prescriptive, leaving freedom to the analysts to determine how the various sources of uncertainties are to be statistically aggregated. Due to the involved use of statistics entangled with heuristic recipes, the resulting safety margins are often difficult to interpret. Also, these technical margins are augmented by additional administrative margins, which are required to ensure compliance with safety standards or regulations, eliminating the incentive to understand their differences. With the new resurgent wave of advanced nuclear systems, e.g., advanced reactors, fuel cycles, and fuel concepts, focused on economizing operation, there is a strong need to develop a clear understanding of the uncertainties and their consolidation methods to reduce them in manners that can be scientifically defended. In response, the current studies compare the analyses behind four notable methodologies for upper subcriticality limit estimation that have been documented in the nuclear criticality safety literature: the parametric, nonparametric, Whisper, and TSURFER methodologies. Specifically, the work offers a deep dive into the various assumptions of the noted methodologies, their adequacies, and their limitations to provide guidance on developing confidence for the emergent nuclear systems that are expected to be challenged by the scarcity of experimental data. Here, to limit the scope, the current work focuses on the application of these methodologies to criticality safety experiments, where the goal is to calculate a bias, a bias uncertainty, and a tolerance limit for k eff in support of determining an upper subcriticality limit for nuclear criticality safety.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

DOE OSTI2021

Feasibility and Prospects of the High Resolution Rapid Refresh Model For Dynamic Line Rating

The ampacity of transmission lines is defined as the maximum amount of current the conductor can safely carry. It is necessary for transmission line operators to apply ampacity limits due to the thermal properties of the conductor. Dynamic Line Rating (DLR) is a technology and technique that uses the environmental conditions or a set of the conditions to calculate the ampacity of the conductor. The way the DLR is calculated has depended upon some amount of physical technology to implement the solution. The High Resolution Rapid Refresh (HRRR) Model is an atmospheric model that may provide enough resolution to perform dynamic line rating without any hardware at all. The objective is to use this model to perform Dynamic Line Rating. This will solve two of the most pervasive problems in Dynamic Line Rating, how to instrument a transmission line hundreds of miles long and then how to predict near future ampacity sufficient for transmission operators to effectively use, This technique will use the Idaho National Laboratory site as the test location for analysis. The High Resolution Rapid Refresh model is a reliable, freely available, and well-maintained weather forecasting model for the entire United States of America. The model provides enough special resolution that additional local weather monitoring devices at location are not necessary to conduct dynamic line rating. The HRRR ratings are lower on average than the locally measured ratings which allows for an additional margin of safety. The adoption of using purely digital methods for dynamic line rating could allow for wide scale adoption of the technology with very little overhead expenditure of implementation.

17 WIND ENERGY↗

NASA NTRS2016

Assessment of Aerosol Optical Property and Radiative Effect for the Layer Decoupling Cases over the Northern South China Sea During the 7-SEAS Dongsha Experiment

The aerosol radiative effect can be modulated by the vertical distribution and optical properties of aerosols, particularly when aerosol layers are decoupled. Direct aerosol radiative effects over the northern South China Sea (SCS) were assessed by incorporating an observed data set of aerosol optical properties obtained from the Seven South East Asian Studies (7-SEAS)/Dongsha Experiment into a radiative transfer model. Aerosol optical properties for a two-layer structure of aerosol transport were estimated. In the radiative transfer calculations, aerosol variability (i.e., diversity of source region, aerosol type, and vertical distribution) for the complex aerosol environment was also carefully quantified. The column-integrated aerosol optical depth (AOD) at 500nm was 0.1-0.3 for near-surface aerosols and increased 1-5 times in presence of upper layer biomass-burning aerosols. A case study showed the strong aerosol absorption (single-scattering albedo (omega) approx. = 0.92 at 440nm wavelength) exhibited by the upper layer when associated with predominantly biomass-burning aerosols, and the omega (approx. = 0.95) of near-surface aerosols was greater than that of the upper layer aerosols because of the presence of mixed type aerosols. The presence of upper level aerosol transport could enhance the radiative efficiency at the surface (i.e., cooling) and lower atmosphere (i.e., heating) by up to -13.7 and +9.6W/sq m2 per AOD, respectively. Such enhancement could potentially modify atmospheric stability, can influence atmospheric circulation, as well as the hydrological cycle over the tropical and low-latitude marginal northern SCS.

biomass burning↗

DOE OSTI2020

Ultrafast dynamics in the high-symmetry and in the charge density wave phase of 2 H –NbSe 2

Here we investigate carrier and collective mode dynamics in 2H-NbSe 2 using time-resolved optical pump-probe spectroscopy and compare the results with first-principles calculations. Broadband ultrafast reflectivity studies of 2H-NbSe 2 in a wide temperature interval covering the normal, charge density wave (CDW) and superconducting phase were performed. Spectral features observed in the transient reflectivity experiment were associated with specific optical transitions obtained from band structure calculations. Displacive excitation of coherent phonons showed CDW-associated coherent oscillations of the soft phonon mode across the whole spectral range. Temperature evolution of this coherent phonon mode in the low-excitation linear regime shows softening of the mode down to the CDW transition temperature T CDW with subsequent hardening below T CDW . The global fit of the broadband probe data reveals four different relaxation times associated with characteristic electron-electron, electron-phonon, and phonon-phonon relaxation processes. From first-principles calculations of electron-phonon coupling we associate the few picosecond electron-phonon relaxation time $\tau_2$ with a specific group of phonons with frequencies around 20 meV. On the other hand, the anomalously long relaxation time of $\tau_3 \sim$ to 100 ps is associated with anharmonicity-driven phonon-phonon scattering. All relaxation processes result from anomalies near the second order CDW phase transition that are reflected in the temperature dependencies of the characteristic relaxation times and amplitudes of coherent oscillations. At the highest fluences we observe electronic melting of the CDW and disappearance of the mode hardening below T CDW .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Modeling gas, hydrates, and slope stability on the U.S. Atlantic margin during Pleistocene glacial cycles

Changes in temperature and sea level can cause dissociation of methane hydrates in shallow marine sediments, leading to seafloor destabilization. Along the U.S. Atlantic margin, there exists a well-documented history of slope failure and numerous recorded occurrences of gas seeps. Several studies have linked slope failure in the region to gas seepage and hydrate dissociation driven by glacial-interglacial transitions, but this linkage has not been quantitatively demonstrated. Along the shelf edge, in an area where shallow methane gas seeps have been identified, we modeled methane gas and hydrate formation using a one-dimensional fluid flow model. Methane gas formation was modeled over the last 120,000 years to simulate a glacial-interglacial cycle. Here, we ran this model at 16,044 individual locations in the region between 29° N – 45° N and 82° W – 66° W at a resolution of 1 x 1 arcminutes, focusing specifically on water depths between 200 and 1000 m that bracket the seafloor outcrop of the base of the hydrate stability zone. Using historic temperature and pressure records from the last 120,000 years, sediment properties in the area, and factor of safety calculations, we found that hydrate dissociation alone is unlikely to cause sediment failure in the region, implying that an additional driving force would be necessary for failure to occur.

Hydrate dissociation↗

DOE OSTI2023

Retrieval of temperature and humidity profiles from ground-based high-resolution infrared observations using an adaptive fast iterative algorithm

Various retrieval algorithms have been developed for retrieving temperature and water vapor profiles from Atmospheric Emitted Radiance Interferometer (AERI) observations. The physical retrieval algorithm, named AERI Optimal Estimation (AERIoe), outperforms other retrieval algorithms in many aspects except the retrieval time, which is significantly increased due to the complex radiative transfer process. The calculation of the Jacobian matrix is the most computationally intensive step of the physical retrieval algorithm. Interestingly, an analysis of the change in AERI observations' information content with respect to Jacobians revealed that the AERIoe algorithm's performance presents negligible dependence on these metrics. Thus, the Jacobian matrix could remain unchanged when the variation in the atmospheric state is small in the retrieval process to reduce the most time-consuming computation. On the basis of the above findings, a fast physical–iterative retrieval algorithm was proposed by adaptively recalculating Jacobians in keeping with the changes in the atmospheric state. Experiments with synthetic observations demonstrate that the proposed method experiences an average reduction in retrieval time by an impressive 59 % compared to the original AERIoe algorithm while achieving maximum root-mean-square errors of less than 0.95 K and 0.22 log(ppmv) for heights below 3 km for the temperature and water vapor profile, respectively. Further analyses revealed that the fast-retrieval algorithm reached an acceptable convergence rate of 98.7 %, marginally lower than AERIoe's 99.9 % convergence rate for the 826 cases used in this study.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗

NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗

NASA NTRSDate not supplied

Geochemical Insights Into Volcanic and Lithospheric Evolution of Mars

Introduction: The Martian lithosphere plays a pivotal role in volcanic processes, altering surface features, climate shifts, and the planet's early habitability. Crucially, the rigid lithosphere significantly influences Mars' long-term thermal conditions. Studying various geological eons’ volcanic compositions unveils how the lithosphere has evolved, even though understanding the early Noachian volcanic chemistry remains challenging due to weathering and resurfacing. In this study, newly discovered Noachian volcanic terranes[1-2] along with Hesperian and Amazonian volcanic terranes[3] are considered(Figure 1), to infer the evolution of the Martian lithosphere. Figure 1:Mars topographic map derived from MOLA-HRSC data [4] highlights the study regions including the Noachian volcanic terranes, Arabia Terra (AT), Thaumasia Minor (TM), and the Southwest Margin of Icaria Planum (SWIP)[5].Additionally, the Hesperian volcanic terranes such as Syrtis Major (SM), Hesperia Planum (HP), and, the Amazonian volcanic terranes encompass Elysium Mons (EM), Alba Patera (AP), and Olympus Mons (OM) are marked. Methods and Datasets: We use the most recent elemental mass fraction maps derived from gamma spectroscopy data from the Mars Odyssey 2001mission [6]to ascertain the bulk composition of the Martian landscape. Utilizing the updated geochemical provinces dataset of Mars[7], we consider bulk chemistry variations within our study regions. The maps, at a resolution of 50×50, encompass elements like Al, Ca, Fe, Si, K, Th, H2O, Cl, and S, detected through Gamma-Ray Spectroscopy (GRS) with decimeter scale depth sensitivity[6].Additional major and minor elements like Mg, Na, Ti, P, and Mn, using mass balance methods [8]serve as inputs for melting and crystallization simulations. We employ the pMELTS model with FMQ-3to FMQ oxygen fugacity to estimate the Pressure(P)and melting degree(F)in various eons of volcanic terranes. We have also used conventional thermobarometric calculations to determine their formation P-T conditions, using silica activity and olivine-melt Mg-exchange thermometry[9-10]. Results: Chemical Composition of Studied Volcanic Terranes: To assess formation conditions, we have investigated if GRS-measured compositions signify primary igneous processes by evaluating weathering through the Chemical Index of Alteration (CIA), K/Th ratios, and ternary plots for aqueous alteration. Consistent K/Th ratios and low CIA (<50) across volcanic regions suggest minimal large-scale weathering, affirming well-preserved igneous material extending at decimeter depths. In addition, the calculated bulk compositions align with Martian mantle equilibrium, confirming studied volcanic terranes-Arabia Terra, Thaumasia Minor, and SWIP -as primary or less altered compositions, reinforcing their representativeness at a regional GRS scale on Mars. We have also estimated the bulk compositions of Hesperian and Amazonian volcanic terranes shown in Figure 1. All volcanic terranes represent the primary or minimally altered composition. Temporal Evolution of Martian Volcanic Terranes: Our results show P-T conditions ranging from 1.3-1.6 GPa and 1350-1390°C for Noachian terranes, while Hesperian volcanic terranes vary from 1.6-1.7 GPa and 1370-1395°C, and Amazonian volcanic terranes range between 1.9-2.8 GPa and 1380-1415°C. Our analyses correspond to lithospheric thicknesses/ depth of melting for Noachian which ranges from110-135 km, 120-235 km for Hesperian terranes, and 160-235 km for Amazonian terranes. Partial melting percentages (F) range from 8-12% for Noachian, 10-11% for Hesperian, and 10-12% for Amazonian. We have also calculated mantle potential temperatures (Tp) to show variations across epochs(Figure 2). Heat flux estimation shows variations from 51-65 mW/m2 for Noachian, 38-45 mW/m2for Hesperian, and 27-39 mW/m2for Amazonian terranes, implying temporal changes in lithospheric thickness and heat flow. Comparison of Noachian volcanic terranes with Hesperian and Amazonian reveals variations in lithospheric thickness and heat flux, indicating potential implications for Martian surface conditions, including volcanic activity, climate evolution, and early habitability. Figure 2: P-T diagram displays formation conditions of studied Martian volcanic regions using GRS data, alongside Noachian-age surface basalt samples from Gusev, Gale, and Jezero. Dashed lines represent Martian mantle potential temperatures, calculated by adding the latent heat of fusion to the equilibrium temperature at which melt and solid mantle coexist and subtracting the Martian adiabatic gradient corresponding to the depth of melt. The solid line shows the Martian mantle solidus [11].Discussion& Implication: The study explores Martian volcanic activity across eons. Uniform lithospheric thickness from the Noachian to Hesperian eons suggests persistent weak plumes until Hesperian, influencing heat flow variations. Noachian's high heat flux, likely due to concentrated radioactive elements, shows spatial heterogeneity, indicating early Mars' differentiated mantle. Volcanic outgassing plays avital role in the formation of the Martian climate, likely sustained water on the surface, influencing warm-wet or cold-dry conditions. Hesperian's favorable lithospheric heat flow (45 mW/m2) supports extended water existence, resembling Earth's conditions. Whereas, Amazonian's diminished heat flux (25 mW/m2) corresponds to water loss and arid climate. Hesperian's similar lithospheric thickness to mid-Noachian implies prolonged water retention, influencing climate and potential for astrobiological research, expanding perspectives for future Martian missions. Future work: As volcanism was active throughout its history. The emergence of Mars' extensive volcanic regions, known as large igneous provinces, which are linked to intermittent volcanic activity driven by mantle plumes, poses challenges to understanding the scale and persistence of these phenomena in Martian convection models. Therefore, understanding the regional-scale changes in eruptive processes within Martian volcanic provinces over time remains unclear and crucial for deciphering the evolution of the Martian interior without Earth-like plate tectonics. We will continue this work with petrological and thermoelastic analyses to study the compositional and spatiotemporal evolution of the Elysium Volcanic Province.

A Rani↗

NASA NTRS2020

The interactive global fire module pyrE (v1.0)

Fires affect the composition of the atmosphere and Earth’s radiation balance by emitting a suite of reactive gases and particles. An interactive fire module in an Earth system model (ESM) allows us to study the natural and anthropogenic drivers, feedbacks, and interactions of open fires. To do so, we have developed pyrE, the NASA GISS (Goddard Institute for Space Studies) interactive fire emissions *module. The pyrE module is driven by environmental variables like flammability and cloud-to-ground lightning, calculated by the GISS ModelE ESM, and parameterized by anthropogenic impacts based on population density data. Fire emissions are generated from the flaming phase in pyrE (active fires). Using pyrE, we examine fire occurrence, regional fire suppression, burned area, fire emissions, and how it all affects atmospheric composition. To do so, we evaluate pyrE by comparing it to satellite-based datasets of fire count, burned area, fire emissions, and aerosol optical depth (AOD). We demonstrate pyrE’s ability to simulate the daily and seasonal cycles of open fires and resulting emissions. Our results indicate that interactive fire emissions are biased low by 32 %–42 %, depending on emitted species, compared to the GFED4s (Global Fire Emissions Database) inventory. The bias in emissions drives underestimation in column densities, which is diluted by natural and anthropogenic emissions sources and production and loss mechanisms. Regionally, the resulting AOD of a simulation with interactive fire emissions is underestimated mostly over Indonesia compared to a simulation with GFED4s emissions and to MODIS AOD. In other parts of the world pyrE’s performance in terms of AOD is marginal to a simulation with prescribed fire emissions.

Fires↗