SEARCH · Engineering Papers

Results for “Dynamic equivalent”

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 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS2005

Growth from Solutions: Kink dynamics, Stoichiometry, Face Kinetics and stability in turbulent flow

1. Kink dynamics. The first segment of a polygomized dislocation spiral step measured by AFM demonstrates up to 60% scattering in the critical length l*- the length when the segment starts to propagate. On orthorhombic lysozyme, this length is shorter than that the observed interkink distance. Step energy from the critical segment length based on the Gibbs-Thomson law (GTL), l* = 20(omega)alpha/(Delta)mu is several times larger than the energy from 2D nucleation rate. Here o is tine building block specific voiume, a is the step riser specific free energy, Delta(mu) is the crystallization driving force. These new data support our earlier assumption that the classical Frenkel, Burton -Cabrera-Frank concept of the abundant kink supply by fluctuations is not applicable for strongly polygonized steps. Step rate measurements on brushite confirms that statement. This is the1D nucleation of kinks that control step propagation. The GTL is valid only if l* <Dk/vk, the diffusion path of a kink that has diffusivity Dk and average growth velocity vk. This is equivalent to supersaturations sigma less than approx. alpha/2l*, where alpha is the building block size. For lysozyme, sigma much less than (1%). Conventionally used interstep distance generated by screw dislocation, 19(omega)alpha/Delta(mu) should be replaced by the very different real one, approx.4l*. 2. Stoichiometry. Kink, and thus step and face rates of a non-Kossel complex molecular monocomponent or any binary, AB, lattice was found theoretically to be proportional to 1/(zeta(sup 1/2) + zeta(sup - 1/2)), where zeta = [B]/[A] is the stoichiometry ratio in solution. The velocities reach maxima at zeta = 1. AFM studies of step rates on CaOxalate monohydrate (kidney stones) from aqueous solution was found to obey the law mentioned above. Generalization for more complex lattice will be discussed. 3. Turbulence. In agreement with theory, high precision in-situ laser interferometry of the (101) KDP crystal face shows step bunching if solution flows parallel to the step flow. The bunch height increases with the distance the bunch travels, i.e. with the face size. However, when the flow rate, u, increases, at u greater than approx. 1 m / s , the average step bunch height decreases as 1/u. The pheonomenon is attributed to the turbulent rather than laminar viscous boundary layer where diffusivity Dt = 0.5u(sub tau),y, i.e. increases linearly with the distance y from the solid face. Friction velocity, u(sub tau) approx. u(sup 7/8). Dramatically larger rate of the mass/heat transport within the turbulent, as compared to the laminar, viscous layer will be discussed.

Chernov, A. A.↗

DOE OSTI2023

Turboexpander for Direct Cooling in Hydrogen Vehicle Fueling Infrastructure

Hydrogen fuel cell electric vehicles (FCEVs) have been identified as one of a few options for zero carbon emissions transportation. A major advantage of FCEVs is that they can fuel quickly and follow a familiar fueling behavior to hydrocarbon-fueled vehicles. Whether light duty or heavy duty, the goal for a hydrogen dispenser is to fuel a vehicle in the same amount of time as the fossil fuel equivalent. When hydrogen is dispensed into the vehicle storage system, however, the temperature rises due to the Joule-Thomson effect and the heat of compression. Typically, vehicles store the compressed hydrogen in composite overwrapped pressure vessels that have a polymer liner with an operational temperature limit of 85°C. This temperature limit can be exceeded during fast fueling if hydrogen is not precooled. Precooling allows for a dispenser to fuel a vehicle at a faster flow rate by preventing the storage tank on the vehicle from overheating. Fueling protocols and requirements are presented in SAE J2601 Fueling Protocols for Light Duty Gaseous Hydrogen Surface Vehicles [1]. A heavy-duty equivalent is under development with similar requirements for precooling. Currently, conventional precooling for light-duty vehicle refueling uses a heat exchanger and chiller to cool the hydrogen gas to -40°C before entering the vehicle. The precooling system represents a significant part of the station capital and operating costs, so if the cost of the precooling system can be reduced by improving its efficiency, the overall station capital and operating cost can be reduced. In this project, National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (SNL) researchers teamed up to investigate the turboexpander precooling application. A turboexpander is a device that places a turbine in a flow path where a pressure differential can be attained. This expansion device will extract work and lower the temperature of the fluid as the pressure reduces. While initial calculations based on established principles showed potential for a turboexpander to generate cooled gas, much work needs to be done to prove the concept. Turboexpanders typically work best under steady state conditions, while the dispenser is a very dynamic flow system. Dynamic turboexpander systems have been proven, such as a turbocharger on a gasoline vehicle. The inlet pressure at a dispenser is also much higher than any other known turboexpander system but should behave similarly to higher density fluids at lower pressures. Having both performed initial calculations, NREL and SNL researchers teamed up to investigate the turboexpander precooling application further. A project was soon built around the idea with SNL performing system modeling using previously proven capabilities and NREL performing hardware characterization with established station capabilities. Creare LLC was contracted as the turbomachinery expert to design and build the concept device. Part way through the project, however, contracting issues with the funding partner caused the project to terminate early before building and characterizing the concept device. While the project could not continue, many key findings were already learned. This paper is a summary of those findings.

08 HYDROGEN↗

DOE OSTI2026

SuperLab 2.0 Showcase: Connecting Five Labs to Tackle Grid Complexity and Unlock Unique Grid Asset Potential

SuperLab 2.0 (5-Lab Demo) is a collaborative, national-scale experiment showcasing the coordination of geographically distributed energy assets in real time. The demonstration integrates 25 physical and digital assets, spanning wind, PV, batteries, electrolyzers, DC fast chargers, microgrid controllers, building automation systems, small modular reactor (SMR), control centers, and gas turbines, across five DOE national laboratories-NLR, INL, NETL, LBNL, and SNL. These assets are unified using Energy Sciences Network (ESnet), a low-latency, high-performance U.S. Department of Energy's (DOE) network, and controlled via a centralized energy controller hosted at NLR's ARIES facility. The demonstration validates the ability to stress-test hybrid energy systems under dynamic scenarios to de-risk advanced control strategies for greater resilience and flexibility. SuperLab 2.0 (5-Lab Demo) showcased a major advancement in federated national laboratory collaboration, enabling real-time, cross-laboratory experimentation to coordinate geographically dispersed distributed energy resources (DERs) using various communication protocols and networks. SuperLab 2.0 (5-Lab Demo) built on previous demonstrations conducted between NLR-PNNL and NLR-INL connecting diverse assets including distant protection devices, a SMR simulator, and a high temperature electrolyzer (HTE). Previous demos were based on a single connection between two labs with minimal coordination challenges. The 5-Lab demo with a centralized controller, distributed testbeds across different geographical locations, and use of protocols-based communication represents a scenario closer to real-world grid operations that coordinate resources across a region to meet system needs. This experiment studied how local DER controllers interact with a centralized energy controller during normal and abnormal events to maintain reliability. The SuperLab team across the five labs implemented a notional power system model equivalent of transmission and distribution lines, represented by the data networks interconnecting the labs. Each lab continuously exchanged local parameters (such as P and Q) from its Hardware-In-Loop (CHIL) and Power Hardware-In-Loop (PHIL) assets through centralized energy controller at NLR, enabling real-time interaction and coordination across sites. By leveraging ESnet as the communication backbone, the team successfully operated the distributed assets as a unified power system, with each bus represented by a different laboratory. This setup mirrors how assets interact in real-world power systems across dispersed locations with various protocols and latencies. At each lab site, assets were operated using their own local controllers which were coordinated through an overarching operation and control layer of centralized energy controller, equivalent to how an energy management system (EMS) orchestrates assets across a regional or national grid. SuperLab's federated connectivity utilized a Digital Real-Time Simulators (DRTS)-type gateway to connect Controller Hardware-In-Loop (CHIL) and PHIL assets between labs. To enable this federated connection through ESnet, a deterministic network was established where latency variations were consistent. This consistency allowed the development of digital filters for the power system assets across CHIL and PHIL interfaces to avoid unstable and unreliable grid conditions. This report provides an overview of the cross-laboratory configuration and offers insights into interconnecting geographically distributed research assets to test them as if they were co-located. This experiment represents a step toward linking nine DOE national laboratories, enabling nation-wide simulations that can address utility-driven challenges with grid resilience, flexibility, and modernization.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Pilot Testing of a Modular Oxygen Production System Using Oxygen Binding Adsorbents (Final Report)

In this project, RTI and Air Liquide have focused on the development of innovative oxygen separation materials and technologies based on reversibly binding of oxygen to enable smaller and cheaper air separations. In alignment with that object the team has had multiple achievements to advance the technology including: Synthesized and characterized novel M-CoSalen (RTIO2Sorb-1) material which showed a dynamic oxygen adsorption capacity in excess of 1wt%. Developed synthesis routes for RTIO2Sorb-1 with commercially relevant techniques and produced batches of 0.25kg and over 4 kg of total synthesis. Synthesized and characterized extrudate forms of the RTIO2Sorb-1 with relevant crush strength and maintained dynamic oxygen adsorption capacity in near 1wt%. Developed techniques for forming structured beds in fiber shapes with conventional materials. Developed a model for O 2 sorption processes to analyze various bed configurations and cycle parameters. Studied classic O 2 solid sorbents reported in the literature to learn oxygen binding mechanisms. Design a 10 kg/day O 2 VPSA pilot system with 2 or 4-bed operation. Completed a techno-economic analysis based on the experimental results and modeling. The overall objective of the project was to design, fabricate, and test a modular O 2 production system and perform a techno-economic analysis (TEA) after testing to determine the cost-benefit of the advanced modular air separation system. The goals of this technology development project were to achieve a bed-size factor (BSF) of less than 600 lb-adsorbent/TPD O 2 (ton/day O 2 ) (as compared with the state-of-the-art BSF of 850), and O 2 -purity greater than 95% at a cost that is projected to be equivalent or lower than the current state of the art (SOTA), commercially available large-scale cryogenic air separation systems. To achieve these goals and objective, the project team executed on (1) oxygen binding adsorbent optimization and scale up, (2) adsorbent material formation process studies to form the adsorbent material into structured beds for rapid pressure swing adsorption (PSA) cycles with low pressure drop, fast mass transfer, and low attrition, (3) cycle development studies to optimize the PSA process, and (4) develop simulation tools for rapid cycle modeling and numerical evaluation/optimization.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2026

Non-Invasive Electrode Impedance Estimation for Optimized Charge Profile Parameterization of Lithium-Ion Batteries

This work presents a non-invasive method for parameterizing a physically motivated equivalent circuit model of lithium-ion batteries using operando electrochemical impedance spectroscopy and time-domain data. The proposed model consists exclusively of linear circuit elements, enabling computationally efficient simulation and real-time implementation on battery management system chips. By integrating frequency- and time-domain measurements, the model accurately estimates internal states such as the negative electrode potential, achieving a root mean square error of 12.3 mV during fast charging validation. Parameterization requires only rate tests with sinusoidal perturbations at three different ambient temperatures, making the approach experimentally accessible. The model reveals key insights into battery behavior, including rate-dependent overpotentials primarily governed by charge transfer kinetics at the positive electrode, and temperature-dependent impedance contributions from both charge transfer and solid-state diffusion processes. Validation using reference electrodes confirms the model’s ability to detect lithium plating onset and reproduce impedance behavior across a wide range of operating conditions. The approach enables in situ optimization of fast charging profiles and lays the foundation for future extensions incorporating aging effects and plating dynamics.

25 ENERGY STORAGE↗

DOE OSTI2023

How thermal fluctuations influence the function of the FeMo cofactor in nitrogenase enzymes

The catalytic mechanism of N 2 fixation by nitrogenase remains unresolved in how the strong N≡N bond is activated and why the reductive elimination of H 2 is required. Here, we use density functional theory and physiologically relevant thermal simulations to elucidate the mechanism of the complete nitrogenase catalytic cycle. Over the accumulation of four reducing equivalents, we find that protons and electrons transfer to the FeMo cofactor to weaken and break its bridge Fe–S bond, leading to temporary H 2 S formation that exposes the Fe sites to weakly bind N 2 . Remarkably, we find that subsequent H 2 formation is responsible for chemical activation to an N=N double bond accompanied by a low barrier for H 2 release. We emphasize that finite temperature effects smooth out mechanistic differences between DFT functionals observed at 0 K, thus leading to a consistent understanding as to why H formation is an obligatory step in N 2 adsorption and activation.

DFT↗

DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

The 1st AIAA Stability and Control Prediction Workshop was created to establish best practices for the prediction of stability & control derivatives using computational fluid dynamics and assess the limitations of these computational methods when those best practices are applied. The inaugural workshop considers the ONERA version of the NASA/Boeing Common Research Model(CRM), which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel tests have been conducted for this configuration with longitudinal tests having been previously published, in addition to unpublished data at small sideslip angles that will serve as ‘blind’ data for workshop data comparisons. Participants were provided a ‘family’ of unstructured grids for the full-span ONERA CRM model with the wind-tunnel sting included. This family of mixed-element grids consists of 5 levels of refinement (tiny, coarse, medium, fine, and extra fine) with surface and volume mesh scaling, resulting in a size range of 14.6 to 53.4 million nodes. In addition, a medium refinement mesh has been provided for the ONERA CRM configuration without a sting to evaluate the sting’s impact on static longitudinal stability characteristics. In addition to these workshop-provided grids, the present work also considers an equivalent ‘family’ of computational grids generated using Heldenmesh™, a rapid grid generation software by Helden Aerospace Corporation for creating high-quality, three-dimensional, mixed-element unstructured meshes. Because of the authors’ familiarity with this software, these additional grids were generated as a comparison to the workshop-provided grids and to better understand the implications of using volume-mirrored grids for stability and control predictions. The present work will contribute to the workshop with test case data generated using the NASA FUN3D code, a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Numerical simulations will be conducted using the Quadratic Constitutive Relationship (QCR) version of the Spalart-Allmaras (SA) turbulence model with negative turbulence variable provisions. Both steady and2nd-order, time-accurate simulation results are to be generated and compared for select test cases, as time permits, to investigate their impact on FUN3D predictions. The present work will consider the three primary workshop test cases: (1) grid convergence study, (2) Mach number effect on static stability, and (3) wind tunnel sting increments. Additionally, data will be provided for the two optional test cases, which include:(1) static stability derivative calculations and (2)sideslip angle sweeps. In each of the test cases, the vehicle is stationary, and the body is assumed to be rigid, where vehicle deformation has been accounted for in the model configuration geometry. For all test cases, longitudinal and lateral force and moment aerodynamic coefficients will be provided for the total configuration, in addition to a component-level breakdown that includes the port wing, starboard wing, fuselage, and tail.

CFD↗

DOE OSTI2022

Seabed bathymetry and friction modeling in MoorDyn

This paper presents the implementation and verification of a seabed bathymetry feature and seabed friction feature to the open-source, lumped-mass mooring system dynamics modeler, MoorDyn, which is part of the National Renewable Energy Laboratory’s aero-hydro-servo-elastic simulation tool, OpenFAST. Variations in seabed slope, as well as the frictional effects of mooring lines moving along the seabed, will affect the mooring line tensions of a floating platform and the consequent platform response. These new features are especially relevant for modeling mooring systems in deep-water coastal areas where seabed depth can change significantly over an entire mooring footprint. The bathymetry feature models the seabed as a rectangular grid of variable water depths in place of the existing, uniform water depth in MoorDyn. The friction force is primarily represented as a Coulombic friction force, or the product of a kinetic friction coefficient and the seabed contact normal force, with the ability to differentiate between transverse and axial motion of a line node on the seabed in any bathymetry grid. These capabilities were tested by running MoorDyn and OpenFAST simulations over a variety of seabed and environmental conditions; the resulting fairlead tensions, node tensions, and mooring line kinematics were verified against equivalent OrcaFlex simulations. The results match closely, meaning the features are verified, which will increase the overall fidelity of OpenFAST and FAST.Farm simulations.

17 WIND ENERGY↗

NASA NTRS2011

Receiver Gain Modulation Circuit

A receiver gain modulation circuit (RGMC) was developed that modulates the power gain of the output of a radiometer receiver with a test signal. As the radiometer receiver switches between calibration noise references, the test signal is mixed with the calibrated noise and thus produces an ensemble set of measurements from which ensemble statistical analysis can be used to extract statistical information about the test signal. The RGMC is an enabling technology of the ensemble detector. As a key component for achieving ensemble detection and analysis, the RGMC has broad aeronautical and space applications. The RGMC can be used to test and develop new calibration algorithms, for example, to detect gain anomalies, and/or correct for slow drifts that affect climate-quality measurements over an accelerated time scale. A generalized approach to analyzing radiometer system designs yields a mathematical treatment of noise reference measurements in calibration algorithms. By treating the measurements from the different noise references as ensemble samples of the receiver state, i.e. receiver gain, a quantitative description of the non-stationary properties of the underlying receiver fluctuations can be derived. Excellent agreement has been obtained between model calculations and radiometric measurements. The mathematical formulation is equivalent to modulating the gain of a stable receiver with an externally generated signal and is the basis for ensemble detection and analysis (EDA). The concept of generating ensemble data sets using an ensemble detector is similar to the ensemble data sets generated as part of ensemble empirical mode decomposition (EEMD) with exception of a key distinguishing factor. EEMD adds noise to the signal under study whereas EDA mixes the signal with calibrated noise. It is mixing with calibrated noise that permits the measurement of temporal-functional variability of uncertainty in the underlying process. The RGMC permits the evaluation of EDA by modulating the receiver gain using an external signal. Without the RGMC, samples of calibrated references from radiometers form an ensemble data set of the natural occurring fluctuations within a receiver. By driving the gain of an otherwise stable receiver with an external signal, the conceptual framework and generalization of the mathematics of EDA can be tested. A series of measurements was conducted to evaluate and characterize the performance of the RGMC. Test signals stepped the RGMC across its dynamic range of performance using a radiometer that sampled four noise references; analysis indicates that the RGMC successfully modulated the receiver gain with an external signal. Calibration algorithms applied to four noise references demonstrate the RGMC produced ensemble data sets of the external signal.

Jones, Hollis↗

DOE OSTI2022

Enabling Production of Algal Biofuels by Techno-Economic Optimization of Co-Product Suites

Recent techno-economic analysis (TEA) has underscored that for algal biofuels to be cost competitive with petroleum fuels, co-products are necessary to offset the cost of fuel production. The co-product suite must scale with fuel production while also maximizing value from the non-fuel precursor components. The co-product suite also depends on algal biomass composition, which is highly dynamic and depends on environmental conditions during cultivation. Intentional shifts in composition during cultivation are often associated with reduced biomass productivity, which can increase feedstock production costs for the algae-based biorefinery. The optimal algae-based biorefinery configuration is thus a function of many factors. We have found that comprehensive TEA, which requires the construction of process models with detailed mass and energy balances, along with a complete accounting of capital and operating expenditures for a commercial-scale production facility, provides invaluable insight into the viability of a proposed biorefinery configuration. This insight is reflected in improved viability for one biorefining approach that we have developed over the last 10 years, namely, the Combined Algal Processing (CAP) approach. This approach fractionates algal biomass into carbohydrate-, lipid-, and protein-rich fractions, and tailors upgrading chemistry to the composition of each fraction. In particular, transitioning from valorization of only the lipids to a co-product suite from multiple components of high-carbohydrate algal biomass can reduce the minimum fuel selling price (MFSP) from more than $8/gallon of gasoline equivalent (GGE) to $2.50/GGE. This paper summarizes that progress and discusses several surprising implications in this optimization approach.

09 BIOMASS FUELS↗

DOE OSTI2023

Modification of Lie's transform perturbation theory for charged particle motion in a magnetic field

It is pointed out that the conventional Lie transform perturbation theory for the guiding center motion of charged particles in a magnetic field needs to be modified for ordering inconsistency. There are two reasons. First, the ordering difference between the temporal variation of gyrophase and that of the other phase space coordinates needs to be taken into account. Second, it is also important to note that the parametric limit of the derivative of a function is not equivalent to the derivative of the limit function. When these facts are taken into account, the near identity transformation rule for one form related to the Lagrangian is modified. With the modified near identity transformation rule, the drift motion of charged particles can be described in the first order, instead of the second order and beyond through a tedious expansion process as in the conventional formulation. Here, this resolves the discrepancy between the direct and Lie transform treatments in the Lagrangian perturbation theory for charged particle motion in a magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

ISMIP6 Antarctica: A Multi-Model Ensemble of the Antarctic Ice Sheet Evolution Over the 21st Century

Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between -7:8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica ass change varies between -6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

Hélène Seroussi↗

DOE OSTI2020

ISMIP6 Antarctica: a multi-model ensemble of the Antarctic ice sheet evolution over the 21st century

Abstract. Ice flow models of the Antarctic ice sheet are commonly used to simulate its future evolution in response to different climate scenarios and assess the mass loss that would contribute to future sea level rise. However, there is currently no consensus on estimates of the future mass balance of the ice sheet, primarily because of differences in the representation of physical processes, forcings employed and initial states of ice sheet models. This study presents results from ice flow model simulations from 13 international groups focusing on the evolution of the Antarctic ice sheet during the period 2015–2100 as part of the Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). They are forced with outputs from a subset of models from the Coupled Model Intercomparison Project Phase 5 (CMIP5), representative of the spread in climate model results. Simulations of the Antarctic ice sheet contribution to sea level rise in response to increased warming during this period varies between −7.8 and 30.0 cm of sea level equivalent (SLE) under Representative Concentration Pathway (RCP) 8.5 scenario forcing. These numbers are relative to a control experiment with constant climate conditions and should therefore be added to the mass loss contribution under climate conditions similar to present-day conditions over the same period. The simulated evolution of the West Antarctic ice sheet varies widely among models, with an overall mass loss, up to 18.0 cm SLE, in response to changes in oceanic conditions. East Antarctica mass change varies between −6.1 and 8.3 cm SLE in the simulations, with a significant increase in surface mass balance outweighing the increased ice discharge under most RCP 8.5 scenario forcings. The inclusion of ice shelf collapse, here assumed to be caused by large amounts of liquid water ponding at the surface of ice shelves, yields an additional simulated mass loss of 28 mm compared to simulations without ice shelf collapse. The largest sources of uncertainty come from the climate forcing, the ocean-induced melt rates, the calibration of these melt rates based on oceanic conditions taken outside of ice shelf cavities and the ice sheet dynamic response to these oceanic changes. Results under RCP 2.6 scenario based on two CMIP5 climate models show an additional mass loss of 0 and 3 cm of SLE on average compared to simulations done under present-day conditions for the two CMIP5 forcings used and display limited mass gain in East Antarctica.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Effects of elastic softening and helium accumulation kinetics on surface morphological evolution of plasma-facing tungsten

Based on a continuous-domain model, capable of accessing the spatiotemporal scales relevant to fuzz formation on the surface of plasma-facing component (PFC) tungsten, we report self-consistent simulation results that elucidate the effects of elastic softening and helium (He) accumulation kinetics on the surface morphological response of PFC tungsten. The model accounts for the softening of the elastic moduli in the near-surface region of PFC tungsten, including both thermal softening at high temperature and softening due to He accumulation upon He implantation. The dependence of the elastic moduli on the He content follows an exponential scaling relation predicted by molecular-dynamics simulations, while the He content in the near-surface region of PFC tungsten evolves according to a first-order saturation kinetics, consistent with experimental and simulation results reported in the literature. We establish that this elastic softening accelerates both nanotendril growth on the PFC surface and the onset of fuzz formation. We also explore the role of the rate of He accumulation to a saturation level in the near-surface region of irradiated tungsten in the onset of fuzz formation. For PFC tungsten surfaces such as W(110) where, under typical irradiation conditions, the characteristic time scale for stress-driven surface diffusion is comparable to the characteristic time scale for He accumulation, we find that accelerating the rate of He accumulation accelerates the growth rate of nanotendrils emanating from the surface. Additionally, we present a systematic parametric study of the PFC surface morphological response to explore its dependence on the He accumulation kinetics that is controlled by the irradiation conditions for low-energy implantation. Finally, we introduce an incubation time for nanotendril growth on the PFC surface, a concept equivalent to that of incubation fluence discussed in the literature, to predict and explain the minimum exposure time required to observe fuzz formation on PFC tungsten surfaces.

36 MATERIALS SCIENCE↗

NASA NTRS1998

Central Plasma Sheet Ion Properties as Inferred from Ionospheric Observations

A method of inferring central plasma sheet (CPS) temperature, density, and pressure from ionospheric observations is developed. The advantage of this method over in situ measurements is that the CPS can be studied in its entirely, rather than only in fragments. As a result, for the first time, comprehensive two-dimensional equatorial maps of CPS pressure, density, and temperature within the isotropic plasma sheet are produced. These particle properties are calculated from data taken by the Special Sensor for Precipitating Particles, version 4 (SSJ4) particle instruments onboard DMSP F8, F9, F10, and F11 satellites during the entire year of 1992. Ion spectra occurring in conjunction with electron acceleration events are specifically excluded. Because of the variability of magnetotail stretching, the mapping to the plasma sheet is done using a modified Tsyganenko [1989] magnetic field model (T89) adjusted to agree with the actual magnetotail stretch at observation time. The latter is inferred with a high degree of accuracy (correlation coefficient -0.9) from the latitude of the DMSP b2i boundary (equivalent to the ion isotropy boundary). The results show that temperature, pressure, and density all exhibit dawn-dusk asymmetries unresolved with previous measurements. The ion temperature peaks near the midnight meridian. This peak, which has been associated with bursty bulk flow events, widens in the Y direction with increased activity. The temperature is higher at dusk than at dawn, and this asymmetry increases with decreasing distance from the Earth. In contrast, the density is higher at dawn than at dusk, and there appears to be a density enhancement in the low-latitude boundary layer regions which increases with decreasing magnetic activity. In the near-Earth regions, the pressure is higher at dusk than at dawn, but this asymmetry weakens with increasing distance from the Earth and may even reverse so that at distances X less than approx. 10 to -12 R(sub E), depending on magnetic activity, the dawn sector has slightly higher pressure. The temperature and density asymmetries in the near-Earth region are consistent with the ion westward gradient/curvature drift as the ions ExB convect earthward. When the solar wind dynamic pressure increases, CPS density and pressure appear to increase, but the temperature remains relatively constant. Comparison with previously published work indicates good agreement between the inferred pressure, temperature, and density and those obtained from in situ data. This new method should provide a continuous mechanism to monitor the pressure, temperature, and density in the magnetotail with unprecedented comprehensiveness.

Wing, Simon↗

NASA NTRS2009

Assessment of the Short-Term Radiometric Stability between Terra MODIS and Landsat 7 ETM+ Sensors

The Landsat 7 (L7) Enhanced Thematic Mapper (ETM+) sensor was launched on April 15th, 1999 and has been in operation for over nine years. It has six reflective solar spectral bands located in the visible and shortwave infrared part of the electromagnetic spectrum (0.5 - 2.5 micron) at a spatial resolution of 30 m. The on-board calibrators are used to monitor the on-orbit sensor system changes. The ETM+ performs solar calibrations using on-board Full Aperture Solar Calibrator (FASC) and the Partial Aperture Solar Calibrator (PASC). The Internal Calibrator Lamp (IC) lamps, a blackbody and shutter optics constitute the on-orbit calibration mechanism for ETM+. On 31 May 2003, a malfunction of the scan-line corrector (SLC) mirror assembly resulted in the loss of approximately 22% of the normal scene area. The missing data affects most of the image with scan gaps varying in width from one pixel or less near the centre of the image to 14 pixels along the east and west edges of the image, creating a wedge-shaped pattern. However, the SLC failure has no impacts on the radiometric performance of the valid pixels. On December 18, 1999, the Moderate Resolution Imaging Spectroradiometer (MODIS) Proto-Flight Model (PFM) was launched on-board the NASA's EOS Terra spacecraft. Terra MODIS has 36 spectral bands with wavelengths ranging from 0.41 to 14.5 micron and collects data over a wide field of view angle (+/-55 deg) at three nadir spatial resolutions of 250 m, 500 in 1 km for bands 1 to 2, 3 to 7, and 8 to 36, respectively. It has 20 reflective solar bands (RSB) with spectral wavelengths from 0.41 to 2.1 micron. The RSB radiometric calibration is performed by using on-board solar diffuser (SD), solar diffuser stability monitor (SDSM), space-view (SV), and spectro-radiometric calibration assembly (SRCA). Through the SV port, periodic lunar observations are used to track radiometric response changes at different angles of incidence (AOI) of the scan mirror. As a part of the AM Constellation satellites, Terra MODIS flies approximately 30 minutes behind L7 ETM+ in the same orbit. The orbit of L7 is repetitive, circular, sunsynchronous, and near polar at a nominal altitude of 705 km (438 miles) at the Equator. The spacecraft crosses the Equator from north to south on a descending node between 10:00 AM and 10:15 AM. Circling the Earth at 7.5 km/sec, each orbit takes nearly 99 minutes. The spacecraft completes just over 14 orbits per day, covering the entire Earth between 81 degrees north and south latitude every 16 days. The longest continuous imaging swath that L7 sensor can collect is for a 14-minute subinterval contact period which is equivalent to 35 full WRS-2 scenes. On the other hand, Terra can provide the entire corresponding orbit with wider swath at any given ETM+ collection without contact time limitation. There are six spectral matching band pairs between MODIS (bands 3, 4, 1, 2, 6, 7) and ETM+ (bands 1, 2, 3, 4, 5, 7) sensor. MODIS has narrower spectral responses than ETM+ in all the bands. A short-term radiometric stability was evaluated using continuous ETM+ scenes within the contact period and the corresponding half orbit MODIS scenes. The near simultaneous earth observations (SNO) were limited by the smaller swath size of ETM+ (187 km) as compared to MODIS (2330 km). Two sets of continuous granules for MODIS and ETM+ were selected and mosaiced based on pixel geolocation information for non cloudy pixels over the North American continent. The Top-of- Atmosphere (TOA) reflectances were computed for the spectrally matching bands between ETM+ and MODIS over the regions of interest (ROI). The matching pixel pairs were aggregated from a finer to a coarser pixel resolution and the TOA reflectance values covering a wide dynamic range of the sensors were compared and analyzed. Considering the uncertainties of the absolute calibration of the both sensors, radiometric stability was verified for the band pairs. The Railroad Valley Playa, Nada (RVPN) was included in the path of this continuous orbit, which served as a verification point between the shortterm and the long-term trending results from previous studies. This work focuses on monitoring the short-term on-orbit stability of MODIS and the ETM+ RSB. It also provides an assessment of the absolute calibration differences between the two sensors over their wide dynamic ranges.

Choi, Taeyoung↗

DOE OSTI2023

Antarctic permafrost processes and antiphase dynamics of cold-based glaciers in the McMurdo Dry Valleys inferred from 10 Be and 26 Al cosmogenic nuclides

Abstract. Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10Be and 26Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10Be and 26Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10Be and 26Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10Be and 26Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea.

58 GEOSCIENCES↗