Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reference Model 2”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Adaptative Site Management for a 115 Acre Chlorinated Solvent Plume with Two Separate Source Areas at Kennedy Space Center, Florida

Background/Objectives. During Resource Conservation and Recovery Act (RCRA) Facility Investigation (RFI) activities, Geosyntec delineated a chlorinated volatile organic compound (CVOC) plume at the National Aeronautics and Space Administration’s (NASA’s) Vehicle Assembly Building (VAB) area located at KSC, Florida. The RFI activities identified an approximate 115-acre dissolved plume (primarily vinyl chloride) and a trichloroethene (TCE) source area in an active aerospace complex that is surrounded by sensitive wetland/waterbodies. Due to the size of the impacted area, the Corrective Measure Design included a multi-component strategy: (i) address the source area via bioremediation; (ii) protect sensitive wetlands from impacted groundwater discharge via biosparging; and (iii) Long Term Monitoring (LTM) of the remaining dissolved plume. After the Corrective Measures implementation (CMI), NASA and Geosyntec worked with Florida Department of Environmental Protection (FDEP) to implement an adaptive site management for the complex, 115-acre site outside of the traditional RCRA process. The adaptive site management approach relied on performing supplemental assessments and implementing Interim Measures (IMs) to further assess and implement remedies over time while working within site and budget constraints, with an overall goal of achieving enough mass reduction to transition the entire site to LTM and eventually achieve site closure. Approach/Activities. After the biosparge barrier was operational and bioremediation within the source area (referred to as Hot Spot 1) achieved the Corrective Action Objective (CAO), supplemental assessment of the area between Hot Spot 1 and the biosparge barrier was performed. The conceptual site model was updated using the supplemental assessment results and an air sparge system IM was designed to treat an approximate 1.2 acre area (referred to as Hot Spot 2). After installation of the air sparge system, supplemental assessment within the remainder of the 115-acre dissolved plume was performed and a second TCE source area was identified. The TCE source area and associated areas with elevated CVOC concentrations (referred to as Hot Spot 3) were delineated and a bioremediation IM was implemented. Also, the downgradient impacts from Hot Spot 3 were adjacent to a sensitive waterbody, and negotiations with the FDEP allowed the area to be monitored using LTM. Results/Lessons Learned. The performance of supplemental assessment activities and implementation of remedial alternatives as IMs allowed NASA to successfully address groundwater impacts over time, while working within the FDEP regulatory framework. The implementation of the CMI and multiple IMs has achieved the following goals: (i) the biosparge barrier has mitigated the potential discharge of impacted groundwater to an adjacent wetland; (ii) enhanced bioremediation within Hot Spot 1 achieved the CAO within 2 years and transitioned the area into LTM; (iii) operation of an air sparge system within Hot Spot 2 removed TCE as a constituent of concern and contributed to a reduction (approximately 43%) in the impacted groundwater area outside the air sparge treatment area (plume collapse); and (iv) bioremediation within Hot Spot 3 removed approximately 80% of the CVOC mass and contributed to a reduction (approximately 47%) in the impacted groundwater area outside the bioremediation IM treatment area. Overall, the adaptive approach is protecting the sensitive water bodies surrounding the complex site and reducing the area of impacted groundwater, which is moving the entire site towards LTM.

Rebecca C Daprato↗
DOE OSTI2025

High-Resolution Fire Weather Index Data for the Conterminous US (1980–2099), Version 1

This dataset presents a suite of high-resolution fire weather index datasets calculated from observation (gridMet, Livneh, Daymet V4), reanalysis (AgERA5), downscaled hydro-climate projections over the conterminous United States (CONUS) based on multiple selected Global Climate Models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). Aside from the daily FWI datasets, we also include a set of FWI extreme indicators at annual, seasonal, and monthly scales, including 1) fwixx: maximum FWI; 2) fwisa: mean FWI. For annual fwisa, it refers to the season with the maximum seasonal average; 3) fwils: Length of fire season over a specified period, where fire season is defined as the days exceeding the median value of the normalized FWI during the reference period (1980-1984); 4) fwixd: Number of extreme fire weather days over a specified period, where extreme day is defined as the day with FWI > the 95th percentile of the FWI during the reference period (1980-1984). All FWI datasets cover 1980-2020 baseline and the model simulated products including the downscaled products additionally include 2021-2099 near-future periods under the high-end (SSP585) emission scenario.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-Sourced Collaboration for the Production and Refining of Rare Elements and Critical Metals (Final Technical Report)

The project objective was to develop a feasible and cost-effective method for recovering rare earth elements (REEs) and critical materials (CMs) from coal and coal byproducts, resulting in high-purity individually separated REEs and CMs. The targeted REEs included Y, Pr, Nd, Gd, Dy, and Sm, with a purity of over 99.5%, while the CMs included Co, Mn, Ga, Sr, Li, Ni, Zn, and Ge, with a purity of over 90%. The project aimed to design a prototype facility capable of producing 1-3 tonnes/day of high-purity REO mixes. The work was divided into four designated circuits: 1) REE extraction and concentration, 2) REE separation and purification, 3) RE metal production, and 4) CM production. To achieve these goals, the project involved 11 tasks, including technology reviews, research, process flow diagram development, mass balance estimation, and preliminary technical-economic analysis. The project team included researchers from the University of Kentucky, University of Alabama and Virginia Tech as well as process specialists from Argonne National Laboratory. MP Materials provided technical support regarding rare earth markets and processing while Alliance Coal performed resource assessment. The project included a market analysis for Nd/Pr, Tb, Dy, Gd, Y, Co, Mn, Li, Sr, Ga, Ni, Zn, and Ge. These analyses provided insights into the supply and demand trends as well as historic and future projections of market price relative to purity requirements for these elements. Two coal resources were selected for the project: the West Kentucky No. 13 (Baker) Seam and an undisclosed lignite resource in the Illinois coal basin. The estimated quantities of REEs in these resources were calculated based on production samples and drilling data. It was estimated that there is adequate supply for an operation producing one metric ton daily of higher purity mixed rare earth oxides (MREO) for approximately 20 years at a site located in western Kentucky. In Circuit 1, project data was obtained from a pilot heap leach and REE concentration facility. It was concluded that the existing circuit, which generated a MREO concentrate, two types of CM mixed products, and Li- and Sr-containing waters, would be suitable feed for circuits 2-4. Data from the first-of-its-kind coal coarse refuse heap leach pilot pad played a crucial role in estimating reliable elemental concentrations of the pregnant leaching solution (PLS). The average total REE concentration in the PLS was found to be 28.6 ppm. In Circuit 2, several concepts were explored including a novel process referred to as solvent-assisted chromatography (SAC). This concept involved a novel columnar reactor that incorporated multiple mixer/settlers, thereby enabling the operation of counter-flowing aqueous and organic phases. Unfortunately, due to project time constraints, a complete fundamental modeling analysis could not be completed to fully evaluate the technology. Molten salt electrowinning was considered as an alternative for circuit 3 following circuit 2 purification circuit utilizing the novel SAC process. A mass and energy balance of Nd reduction to metal in a fluoride containing molten salt electrolyte was conducted. Comparisons were made with the current state of Asian molten salt electrorefining, and potential improvements in siphoning rare earth metals (REM) from the reactor were presented. A cost estimate was performed for the production of 1 tonne per day, which yielded a total of $2.29 million for the nine electrowinning (EW) cells required. The selected option for circuits 2 and 3 was a plasma distillation process, which initially separates rare earth elements (REEs) from other elements. This is followed by selective electrowinning in various ionic liquids. The selection was made on the basis of thermodynamic modeling and experimental data previously published by a project partner. The combination offers an innovative approach to integrated refining and RE metal production. For Circuit 4, an extensive literature review was conducted for the processing of the CMs. The ultimate decision was to utilize a combined plasma and ionic liquid process as well to produce individual high-purity concentrates of Zn, Ni, Co, Mn, and Mg. A separate flowsheet for Li and Sr was recommended, which would yield carbonates of these elements. Due to the lack of suitable experimental data at this time, a process recommendation could not be provided but several methods have been proposed for consideration. Lastly, a techno-economic analysis (TEA) was conducted to assess the effectiveness of the proposed process for further investigation. The TEA results revealed a capital expense (CapEx) of $737 million and an annual operational expense (OpEx) of $220 million. Due to the selected elements, the hypothetical heap leach pad can produce 1 metric tonne per day of REO equivalent, but a conscious decision was made to only treat targeted REEs, resulting in the production of 0.4 metric tonne of REM. An estimated annual revenue of $90.87 million was projected based on standard market pricing information provided by the funding agency. During the TEA, ten different modules were evaluated for costing purposes. The precipitation circuit was identified as the largest single operational expense, followed by the Mg/Mn process due to the amount of treated metal. In terms of capital expenditures, the heap leach process incurred the highest cost, followed by the Mg/Mn process. The scalability of the plasma process is a crucial consideration since the reactors cannot be scaled beyond the largest demonstrated size due to their reliance on surface area of the slag and vapor phase. The purity estimate for the REEs are generally 98%±2% to produce a metal. The purity level being lower than the project objective was due to the lack of specific experimental data needed to tighten the tolerance of the estimates. Based on literature and previous experience, the CMs are estimated as follows; Ga (95%+, metal), Sr (95%+, carbonate), Li (95%+, carbonate), Ni (98%±2%, metal), Zn (95%+, metal sponge), Ge (95%+, metal), Co (98%±2%, metal), and Mn (98%±2%, metal).

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1997

Global Plate Velocities from the Global Positioning System

We have analyzed 204 days of Global Positioning System (GPS) data from the global GPS network spanning January 1991 through March 1996. On the basis of these GPS coordinate solutions, we have estimated velocities for 38 sites, mostly located on the interiors of the Africa, Antarctica, Australia, Eurasia, Nazca, North America, Pacific, and South America plates. The uncertainties of the horizontal velocity components range from 1.2 to 5.0 mm/yr. With the exception of sites on the Pacific and Nazca plates, the GPS velocities agree with absolute plate model predictions within 95% confidence. For most of the sites in North America, Antarctica, and Eurasia, the agreement is better than 2 mm/yr. We find no persuasive evidence for significant vertical motions (less than 3 standard deviations), except at four sites. Three of these four were sites constrained to geodetic reference frame velocities. The GPS velocities were then used to estimate angular velocities for eight tectonic plates. Absolute angular velocities derived from the GPS data agree with the no net rotation (NNR) NUVEL-1A model within 95% confidence except for the Pacific plate. Our pole of rotation for the Pacific plate lies 11.5 deg west of the NNR NUVEL-1A pole, with an angular speed 10% faster. Our relative angular velocities agree with NUVEL-1A except for some involving the Pacific plate. While our Pacific-North America angular velocity differs significantly from NUVEL-1A, our model and NUVEL-1A predict very small differences in relative motion along the Pacific-North America plate boundary itself. Our Pacific-Australia and Pacific- Eurasia angular velocities are significantly faster than NUVEL-1A, predicting more rapid convergence at these two plate boundaries. Along the East Pacific Pise, our Pacific-Nazca angular velocity agrees in both rate and azimuth with NUVFL-1A.

Larson, Kristine M.↗
DOE OSTI2016

Maternal body burden of cadmium and offspring size at birth

Increasing evidence suggests an inverse association between cadmium (Cd) and size at birth, potentially greatest among female neonates. We evaluated whether greater maternal body burden of Cd is associated with reduced neonatal anthropometry (birthweight, birth length, head circumference, and ponderal index) and assessed whether these associations differ by infant sex. The analytic sample for the present study (n=396) was derived from a subcohort of 750 women randomly drawn from among all participants (N=4344) in the Omega Study, a prospective pregnancy cohort. Creatinine-corrected Cd in maternal clean-catch spot urine samples (U-Cd) was quantified by inductively coupled plasma mass spectrometry. Continuous log{sub 2}-transformed Cd (log{sub 2}-Cd) and U-Cd tertiles (low<0.29 μg/g creatinine, middle 0.29–0.42 μg/g creatinine, high≥0.43 μg/g creatinine) were used in multivariable linear regression models. Females had reduced birth length with greater U-Cd tertile, whereas males birth length marginally increased [β(95% CI) females: low=reference, middle=−0.59 cm (−1.37, 0.19), high=−0.83 cm (−1.69, 0.02), p-trend=0.08; males: low=reference, middle=0.18 cm (−0.59, 0.95), high=0.78 cm (–0.04, 1.60), p-trend=0.07; p for interaction=0.03]. The log{sub 2}-Cd by infant sex interaction was statistically significant for ponderal index [p=0.003; β(95% CI): female=0.25 kg/m{sup 3} (−0.20, 0.70); male=−0.63 kg/m{sup 3} (−1.01, −0.24)] and birth length [p<0.001; β(95% CI): female=−0.47 cm (−0.74, −0.20), male=0.32 cm (0.00, 0.65)]. Our findings suggest potential sex-specific reversal of Cd’s associations on birth length and contribute to the evidence suggesting Cd impairs fetal growth. - Highlights: • Cadmium levels in the general population potentially adversely affect size at birth. • Maternal urinary cadmium was inversely related to birth length among female infants. • For male infants, maternal cadmium was positively associated with birth length. • Maternal cadmium was also inversely related to ponderal index among male infants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Feeling the Pulse of the Stratosphere: An Emerging Opportunity for Predicting Continental-Scale Cold-Air Outbreaks 1 Month in Advance

Extreme weather events such as cold-air outbreaks (CAOs) pose great threats to human life and the socioeconomic well-being of modern society. In the past, our capability to predict their occurrences has been constrained by the 2-week predictability limit for weather. Here, we demonstrate here for the first time that a rapid increase of air mass transported into the polar stratosphere, referred to as the pulse of the stratosphere (PULSE), can often be predicted with a useful degree of skill 4–6 weeks in advance by operational forecast models. We further show that the probability of the occurrence of continental-scale CAOs in midlatitudes increases substantially above normal conditions within a short time period from 1 week before to 1–2 weeks after the peak day of a PULSE event. In particular, we reveal that the three massive CAOs over North America in January and February of 2014 were preceded by three episodes of extreme mass transport into the polar stratosphere with peak intensities reaching a trillion tons per day, twice that on an average winter day. Therefore, our capability to predict the PULSEs with operational forecast models, in conjunction with its linkage to continental-scale CAOs, opens up a new opportunity for 30-day forecasts of continental-scale CAOs, such as those occurring over North America during the 2013/14 winter. A real-time forecast experiment inaugurated in the winter of 2014/15 has given support to the idea that it is feasible to forecast CAOs 1 month in advance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

From hard spheres to hard-core spins

A system of hard spheres exhibits physics that is controlled only by their density. This comes about because the interaction energy is either infinite or zero, so all allowed configurations have exactly the same energy. The low-density phase is liquid, while the high-density phase is crystalline, an example of “order by disorder” as it is driven purely by entropic considerations. Here we study a family of hard spin models, which we call hard-core spin models, where we replace the translational degrees of freedom of hard spheres with the orientational degrees of freedom of lattice spins. Their hard-core interaction serves analogously to divide configurations of the many spin system into allowed and disallowed sectors. We present detailed results on the square lattice in d = 2 for a set of models with $\mathbb{Z}_n$ symmetry, which generalize Potts models, and their U(1) limits, for ferromagnetic and antiferromagnetic senses of the interaction, which we refer to as exclusion and inclusion models. As the exclusion and inclusion angles are varied, we find a Kosterlitz-Thouless phase transition between a disordered phase and an ordered phase with quasi-long-ranged order, which is the form order by disorder takes in these systems. These results follow from a set of height representations, an ergodic cluster algorithm, and transfer matrix calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
DOE OSTI2016

MO-FG-202-07: Real-Time EPID-Based Detection Metric For VMAT Delivery Errors

Purpose: To create and test an accurate EPID-frame-based VMAT QA metric to detect gross dose errors in real-time and to provide information about the source of error. Methods: A Swiss cheese model was created for an EPID-based real-time QA process. The system compares a treatmentplan- based reference set of EPID images with images acquired over each 2° gantry angle interval. The metric utilizes a sequence of independent consecutively executed error detection Methods: a masking technique that verifies infield radiation delivery and ensures no out-of-field radiation; output normalization checks at two different stages; global image alignment to quantify rotation, scaling and translation; standard gamma evaluation (3%, 3 mm) and pixel intensity deviation checks including and excluding high dose gradient regions. Tolerances for each test were determined. For algorithm testing, twelve different types of errors were selected to modify the original plan. Corresponding predictions for each test case were generated, which included measurement-based noise. Each test case was run multiple times (with different noise per run) to assess the ability to detect introduced errors. Results: Averaged over five test runs, 99.1% of all plan variations that resulted in patient dose errors were detected within 2° and 100% within 4° (∼1% of patient dose delivery). Including cases that led to slightly modified but clinically equivalent plans, 91.5% were detected by the system within 2°. Based on the type of method that detected the error, determination of error sources was achieved. Conclusion: An EPID-based during-treatment error detection system for VMAT deliveries was successfully designed and tested. The system utilizes a sequence of methods to identify and prevent gross treatment delivery errors. The system was inspected for robustness with realistic noise variations, demonstrating that it has the potential to detect a large majority of errors in real-time and indicate the error source. J. V. Siebers receives funding support from Varian Medical Systems.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

Basic research for the geodynamics program

Work performed and data obtained in geodynamic research is reported. The purpose was to obtain utilization of: (1) laser and very long baseline interferometry (VLBI); (2) range difference observation in geodynamics; (3) development of models for ice sheet and crustal deformations. The effects of adopting new precession, nutation and equinox corrections on the terrestrial reference frame are investigated.

Mueller, I. I.↗
NASA NTRS1990

Determination of optimum fin profile for a zero-G capillary drained condenser

This paper presents the analytical formulation and numerical results for heat transfer in a high heat flux condenser that relies on capillary flow along shaped fins (Gregorig surfaces) and a drainage network embedded in the condenser walls. Results are shown for a variety of fin profile shapes in order to show the geometric trade-offs involved in seeking a maximum effective heat transfer coefficient for the fin. Predictions of the model show excellent agreement with previously reported measurements for steam. Based on this work, a profile has been selected for a 2 kW ammonia condenser currently under development for use in space. In that design the fin half width is 0.5 mm and the model predicts a heat transfer coefficient referred to the base of the fin of 9 W/sq cm deg C for a heat flux of 10/W sq cm at the base.

Mccormick, John A.↗
NASA NTRS2014

Expansion of the Real-Time SPoRT-Land Information System for NOAA/National Weather Service Situational Awareness and Local Modeling Applications

The NASA Short-term Prediction Research and Transition (SPoRT) Center in Huntsville, AL is running a real-time configuration of the Noah land surface model (LSM) within the NASA Land Information System (LIS) framework (hereafter referred to as the "SPoRT-LIS"). Output from the real-time SPoRT-LIS is used for (1) initializing land surface variables for local modeling applications, and (2) displaying in decision support systems for situational awareness and drought monitoring at select NOAA/National Weather Service (NWS) partner offices. The experimental CONUS run incorporates hourly quantitative precipitation estimation (QPE) from the National Severe Storms Laboratory Multi- Radar Multi-Sensor (MRMS) which will be transitioned into operations at the National Centers for Environmental Prediction (NCEP) in Fall 2014.This paper describes the current and experimental SPoRT-LIS configurations, and documents some of the limitations still remaining through the advent of MRMS precipitation analyses in the SPoRT-LIS land surface model (LSM) simulations.

Case, Jonathan L↗
DOE OSTI2022

Discrepancy quantification between experimental and simulated data of CO 2 adsorption isotherm using hierarchical Bayesian estimation

Here, to quantitatively analyze the inconsistencies commonly observed between experimental and simulated adsorption isotherms, parameter estimation of adsorption isotherm models was conducted by hierarchical Bayesian estimation with parameter uncertainties being quantified as probability distributions. The estimation method was implemented using Markov Chain Monte Carlo (MCMC) to analyze multiple data sets obtained from different sources, including a publicly available database. To describe the discrepancies of experimental and simulated adsorption data, the simulation data was set as the reference to which experimental measurements were compared. We applied the proposed approach to analyze CO 2 adsorption isotherms that are measured and simulated on zeolite 13X and MIL-101(Cr). In these case studies, the discrepancy of CO 2 adsorption isotherm was successfully quantified between experimental measurements and predictions given by molecular simulations using Grand Canonical Monte Carlo (GCMC), where uncertainties were quantified as probability distributions. Furthermore, experimental data sets that agree well with the GCMC simulation have been identified, providing insights into experimental and measurement methods as well as choosing the right assumptions in the molecular simulation.

42 ENGINEERING↗
DOE OSTI2025

Multiphysics Meshfree Degradation Modeling of Energy Storage Materials with Kernel Enrichment

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Renewable Energy Laboratory (NREL), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. Reference: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

25 ENERGY STORAGE↗
DOE OSTI2026

Kernel Enriched Meshfree Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Laboratory of the Rockies (NLR), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. References: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Extension and Demonstration of NEAMS Multiphysics Tools to Lead-Cooled, Sodium-Cooled, and Molten Salt Fast Reactor Applications

The SHARP toolkit is a high-fidelity reactor simulation tool developed under the U.S. Department of Energy, Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Campaign. SHARP toolkit is comprised of the neutronics module PROTEUS thermal hydraulics module Nek5000, and structural mechanics module Diablo. During FY17 and FY18, the PROTEUS and Nek5000 components of SHARP were applied to solve challenging sodium-cooled fast reactor (SFR) problems. In particular, selected hot channel factors (HCF) for a prototype metal-fueled SFR design (the AFR-100) were analyzed in high fidelity, and the “SHARP zooming capability” for SFRs was developed and demonstrated to reduce computational expense for full core problems in cases where detailed data is needed in selected fuel assemblies. After the previous success applying SHARP to challenging SFR problems, the focus in FY19 expanded to additional fast reactor applications: lead cooled fast reactors (LFR), sodium cooled fast reactors (SFR), and fast molten salt reactors (MSR). The specific technical tasks were (1) assessment of hot channel factors for LFR, for which no data currently exists, (2) implementation of a gamma transport capability in the PROTEUS solvers to extend the zooming capability to non-fueled assemblies, (3) demonstration of zooming capability in assemblies of the Versatile Test Reactor (VTR), and (4) development of a new coupled capability for simulating fast MSRs, specifically modeling precursor transport inside and outside the core. The progress on tasks (1) and (3) are discussed primarily in this document. Tasks (2) and (4) are briefly mentioned in this document, and the interested reader is referred to separate reports for detailed information. First-of-a-kind hot channel factor (HCF) estimation for LFR with high fidelity codes (PROTEUS/Nek5000) was successfully demonstrated in this study. Selected HCF were computed and compared with SFR data (AFR-100, EBR-II). The findings confirm that different reactor types, design parameters and uncertainties lead to different HCFs. Careful estimation of HCF for a specific design is necessary to obtain appropriate HCFs. In addition to improvement in HCF accuracy, high fidelity tools generate data to help the designer better understand the mechanism of the impact from these uncertainties. For example, the impact of cladding thickness manufacturing tolerance resulted in non-intuitive effects in the corner pins of the LFR assembly. This procedure of computing HCF using high fidelity models shows promise and flexibility for being repeated for any arbitrary reactor of choice, since these tools are targeted to solve many reactor types and geometries. Progress was made towards extending the previously demonstrated SHARP zooming capability to non-fueled SFR assemblies. In particular, a gamma transport capability was implemented in both high fidelity PROTEUS solvers in order to accurately account for heat deposition caused by gamma particles, which accounts for ~10% of total core power. Neutronics verification cases were carried out for a candidate Versatile Test Reactor (VTR) design using the new gamma transport capability in PROTEUS. Comparisons were made with continuous energy MCNP calculations and shown to agree well. The models for the full core design with heterogeneous control and fuel assemblies is in progress and will be completed in FY20. In the fast molten salt reactor (MSR) domain, a new capability was developed to model the flow of neutron precursors both inside and outside of the core, which impacts the effective delayed neutron fraction of the system. The PROTEUS-SN and PROTEUS-NODAL solvers were coupled to Nek5000 in order to simulate precursor flow. The workflows were demonstrated on test problems. The fast MSR modeling capability was supported under two other work packages, and comprehensive details are included in a separate report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2017

OVERFLOW Contribution to HiLiftPW-3

We plan to perform the following sets of computations: For all our contributions (except where stated) Code: OVERFLOW, Turbulence model: SAnegRCQCR2000. - 1. Results will be submitted for both the full chord flap gap (Case 1a) and partially-sealed Chord Flap gap (Case 1c): 1. Grid Refinement Study; 2. Grids: structured overset grids supplied by HiLiftPW committee; 3. Connectivity: Domain Connectivity Framework, DCF. - 2. Results will be submitted for JAXA Standard Model and Nacelle/Pylon Off (Case 2a), Nacelle/Pylon On (Case 2c): 1. Alpha Study; 2. Grids: structured overset grids supplied by HiLiftPW committee; 3. Connectivity: Pegasus 5 (Peg5). - 3. A study of the effects of different connectivity paradigms: 1. DCF vs Peg5 for HLCRM cases; 2. DCF vs. C3P (NASA Ames) vs. Peg5 for JSM cases; 3. JSM grids will be the focus where we will hopefully see some type of trends with reference to wind tunnel data. - 4. Adaption cases will be attempted for (and submitted where appropriate): 1. Cases 1c,1d: HLCRM; 2. Cases 2c and 2d: JSM; 3. Grid: Near Body grids provided by committee, OffBody grids Cartesian; 4. AMR NearBody and OffBody Adaption. - 5. Case 3 Turbulence model verification study: 1. Grid: Series of 3 finest grids as defined on http://turbmodels.larc.nasa.gov/airfoilwakeverif.html; 2. Turbulence models: SAneg and SAneg— RCQCR2000. OVERFLOW 2.2 is a Reynolds-averaged Navier-Stokes (RANS) code developed by NASA...

OVERFLOW↗
DOE OSTI2021

Effects of organized convection parameterization on the MJO and precipitation in E3SMv1. Part I: Mesoscale heating

Mesoscale organization of convection is typically not represented in global circulation models, and hence its influence on the global circulation is not accounted for. The heating component of a parameterization that represents the dynamical and physical effects of circulations associated with organized convection, referred to as the multiscale coherent structure parameterization (MCSP), is implemented in the Energy Exascale Earth System Model version 1 (E3SMv1). Numerical simulations are conducted to assess its impact on the simulated climate. Besides E3SMv1 simulations, we performed high-resolution (1 km) simulations using the Weather Research and Forecasting (WRF) Model to determine the temperature tendencies induced by mesoscale convective systems embedded in deep convection. We tuned the free parameters of the MCSP based on the WRF simulations. MCSP heating enhances Kevin wave spectra in E3SMv1, improves the representation of the Madden-Julian Oscillation, and impacts on mean precipitation that vary with locations and seasons, including a reduction of biases over the tropical Pacific.

Climate Modeling↗