Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “advanced modeling and simulation”

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

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

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Software For Advanced Large-scale Analysis Of Magnetic Confinement For Numerical Design, Engineering & Research (salamander)

As magnetic confinement fusion energy gains traction internationally to enable abundant energy production, designing components for fusion systems is a pressing challenge. During the planned lifetime of a fusion device, components evolve in extreme environments and must withstand large, repeated thermal loads and bombardment by 14 MeV neutrons, plasma ions, and neutral particles (deuterium, tritium, and helium), corrosive conditions, etc. All these physical processes take place simultaneously, interact in intricate ways, and impose important constraints that can affect performance. Experimental data is rare and costly to obtain, making design particularly challenging. Predictive computational frameworks must be an integral part of an accelerated and cost-effective design process by modeling fusion system performance in simulated environments. To better understand component degradation and operational impacts on their performance, the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering & Research (SALAMANDER) is designed as an open-source, fully integrated, multiphysics, multiscale, NQA-1 compliant framework facilitating 3D, high-fidelity fusion system modeling. To that end, SALAMANDER is a MOOSE-based framework, and therefore leverages MOOSE upstream libraries such as PETSc and libMesh to deliver sophisticated finite element, finite volume, and nonlinear solver technology for fusion energy simulations. SALAMANDER couples MOOSE physics module capabilities—such as thermal hydraulics, heat conduction, Navier-Stokes, and thermomechanics—with tritium transport via TMAP8, neutronics via Cardinal, and nascent particle-in-cell capabilities. Direct simulation Monte Carlo methods will be used to address neutral transport near the walls. By coupling all these physics in an integrated application, SALAMANDER will enable high-fidelity modeling of irradiation levels and plasma exposure conditions of plasma facing components and their impact on heat and tritium distributions, as well as the resulting mechanical constraints experienced by the plasma facing components and performance of blanket systems. Furthermore, SALAMANDER will be particularly suited for engineering studies thanks to the stochastic tool module readily available in MOOSE, allowing for extended uncertainty quantification and risk analysis studies. It is also able to use computer-aided design (CAD) meshes to model complex geometries, which is indispensable for fusion systems. SALAMANDER therefore supports design, safety, engineering, and research projects for magnetic confinement fusion systems

Simon, Pierre-Clement [Idaho National Laboratory (↗
NASA NTRSDate not supplied

Particle Interaction Physics Model Formulation for Plume-Surface Interaction Erosion and Cratering

As part of the Game Changing Development (GCD) Program, funded by NASA’s Space Technology Mission Directorate (STMD), the development of simulation capability for the prediction of extra-terrestrial Plume Surface Interaction (PSI) environments has been undertaken by the Fluid Dynamics Branch at NASA/MSFC. The Predictive Simulation Capability (PSC) Element is focused on creating simulation capability for the reliable and accurate prediction of PSI in Martian (~650 Pa) and Lunar (vacuum) ambient environments. In addition to the predictive simulation capability, the GCD Program also contains a companion Ground Testing Element for development of focused datasets for validation of predictive capability as well as a Flight-focused Instrumentation Element. This paper will present the status of implementing and maturing particle-particle interaction constituent physics models essential in simulating the landing surface granular material flow under PSI effects. This gas-particle multi-phase interaction modeling of plume impingement flow on the extra-terrestrial soil material is performed with the Gas-Granular Flow Solver (GGFS) addressed in a companion paper. The response of regolith particle flow induced by lander PSI requires accurate representation of the regolith granular material fluidic behavior and gas-granular interactions. The lunar regolith, as the extreme example, is poorly sorted with broad particle size distributions and large fines content. It has significant cohesion, due to interlocking particle shapes for the very jagged particles. The combination of particle shape and size distribution has been identified as major drivers in the complex particle flow response and resulting crater shape characteristics of extraterrestrial granular material. Constituent models for spherical particles can be formulated directly from particle kinetics theory. Complex particle shapes can be modeled by gluing together elemental spherical shapes into composite particles, requiring a Discrete Element Model (DEM) particle kinetics modeling approach to extract data and formulate constituent models. Mixture constituent models for poly-disperse mixtures (i.e, containing distribution of particle sizes) have recently been developed. The required non-spherical particle mixture granular material response closure models are then obtained through small-scale unit physics DEM simulations for the range of particle shapes, mixtures and packing densities. The granular material response closure models are then implemented in the Eulerian granular flow formulation. This DEM-based constituent model extraction process and formulation of poly-disperse particle mixtures has been successfully developed by small business and academic partners in the development of the Gas-Granular Flow Solver (GGFS) simulation program simulation framework. The currently implemented capabilities have reached the capability level of modeling bi-disperse, non-spherical particle mixtures is being continuously extended towards computational modeling of full range irregular particle mixtures. Under the GCD project, this technology is being further developed, transferred to NASA analysts, and matured towards application readiness. The predictive simulation capability team under the GCD project has acquired the modeling tools and processes of the DEM based constituent model formulation from the GGFS development team and is developing the capability to replicate the existing process. This is the first important step towards the ability of the NASA team to independently perform such model development in a production setting. Further efforts are underway to migrate the DEM based model simulation process performed with the academic based tools to more capable Open Source, highly parallelized simulation tools for efficient operation on NASA HPC assets. Evaluation of the currently implemented (such as mono-disperse and bi-disperse spherical and irregular shape particle constituent model applications) and continuously evolving full-range particle physics models in the GGFS tool is performed by the NASA team to advance application readiness of the simulations. Application testing for complex PSI erosions and cratering scenarios such as the Apollo LM is performed for axi-symmetric and full 3D simulations to aid the tool developers in achieving practical application readiness for NASA projects. Important validation and application testing will further be performed against experimental data generated under the GCD PSI project experimental component.

Peter A Liever↗
NASA NTRS1977

Aerodynamic characteristics of an F-8 aircraft configuration with a variable camber wing at Mach numbers from 0.70 to 1.15

A 0.1-scale model of an F-8 aircraft was tested in the Ames 14-Foot Transonic Wind Tunnel at Mach numbers from 0.7 to 1.15. Angle of attack was varied from -2 deg. to 22 deg. at sideslip angles of 0 deg and -5 deg. Reynolds number, dictated by the atmospheric stagnation pressure, varied with Mach number from 3.4 to 4.0 million based on mean aerodynamic chord. The model was configured with a wing designed to simulate the downward deflection of the leading and trailing edges of an advanced-technology-conformal-variable camber wing. This wing was also equipped with conventional (simple hinge) flaps. In addition, the model was tested with the basic F-8 wing to provide a reference for extrapolating to flight data. In general, at all Mach numbers the use of conformal flap deflections at both the leading edge and trailing edge resulted in slightly higher maximum lift coefficients and lower drag coefficients than with the use of simple hinge flaps. There were also found to be small improvements in the pitching-moment characteristics with the use of conformal flaps.

Boltz, F. W.↗
DOE OSTI2026

Optimization of well design and CO 2 injection strategy for risk reduction in Class VI geological carbon sequestration wells

The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.

25 ENERGY STORAGE↗
DOE OSTI2021

Surrogate Model of Flexible Research Platform EnergyPlus Models to Enable Sensitivity Analysis

This letter report describes the surrogate models developed from the EnergyPlus model of Oak Ridge National Laboratory’s Flexible Research Platform. Two data-driven black-box models were developed, and the outputs of the surrogate models were compared with the EnergyPlus model. The two models developed are a multilayer perceptron deep learning model, and a long short-term memory (LSTM) neural network model. The three factors for selecting the black-box models are scalability, computation time, and accuracy. A total of 107 input variables were the dominant variables in determining the outputs of building energy consumptions and thermal comfort. A total of 54 output variables were identified as the prediction targets, including the system- and zone-level outputs. The large set of the simulation cases were generated by integrating sensor errors into an emulator based on EnergyPlus and Python EMS, which includes advanced control sequences from ASHRAE Guideline 36-2018: High-Performance Sequences of Operation. The surrogate models were developed based on a set of large-scale simulation runs (i.e., 4,000 runs) on a cloud platform. The comparison analysis shows that the two black-box models had good accuracy for predicting new outputs for sensitivity analysis using the root mean square error metric. As a next step, the developed surrogate models will be used to perform sensitivity analysis for different sensor impacts (e.g., sensor types, sensor locations).

42 ENGINEERING↗
DOE OSTI2024

Sienna Modeling Framework [Slides]

NREL's Sienna modeling framework effectively builds, solves, and analyzes the scheduling problems and dynamic simulations of quasi-static infrastructure systems. It uses a modular framework to answer different questions about future energy systems, fundamentally advancing the nation's ability to model individual and integrated infrastructure systems at a range of spatial and temporal scales. This presentation will include NREL power grid researcher Clayton Barrows.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Examining the Impacts of Great Lakes Temperature Perturbations on Simulated Precipitation in the Northeastern United States

Most inland water bodies are not resolved by general circulation models, requiring that lake surface temperatures be estimated. Given the large spatial and temporal variability of the surface temperatures of the North American Great Lakes, such estimations can introduce errors when used as lower boundary conditions for dynamical downscaling. Lake surface temperatures (LSTs) influence moisture and heat fluxes, thus impacting precipitation within the immediate region and potentially in regions downwind of the lakes. For this study, the Advanced Research version of the Weather Research and Forecasting Model (WRF-ARW) was used to simulate precipitation over the six New England states during a 5-yr historical period. The model simulation was repeated with perturbed LSTs, ranging from 10°C below to 10°C above baseline values obtained from reanalysis data, to determine whether the inclusion of erroneous LST values has an impact on simulated precipitation and synoptic-scale features. Results show that simulated precipitation in New England is statistically correlated with LST perturbations, but this region falls on a wet–dry line of a larger bimodal distribution. Wetter conditions occur to the north and drier conditions occur to the south with increasing LSTs, particularly during the warm season. Additionally, the precipitation differences coincide with large-scale anomalous temperature, pressure, and moisture patterns. Care must therefore be taken to ensure reasonably accurate Great Lakes surface temperatures when simulating precipitation, especially in southeastern Canada, Maine, and the mid-Atlantic region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2013

The Digital Astronaut Project Computational Bone Remodeling Model (Beta Version) Bone Summit Summary Report

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur [1-3]. The most commonly used countermeasure against bone loss in microgravity has been prescribed exercise [4]. However, data has shown that existing exercise countermeasures are not as effective as desired for preventing bone loss in long duration, 4 to 6 months, spaceflight [1,3,5,6]. This spaceflight related bone loss may cause early onset of osteoporosis to place the astronauts at greater risk of fracture later in their lives. Consequently, NASA seeks to have improved understanding of the mechanisms of bone demineralization in microgravity in order to appropriately quantify this risk, and to establish appropriate countermeasures [7]. In this light, NASA's Digital Astronaut Project (DAP) is working with the NASA Bone Discipline Lead to implement well-validated computational models to help predict and assess bone loss during spaceflight, and enhance exercise countermeasure development. More specifically, computational modeling is proposed as a way to augment bone research and exercise countermeasure development to target weight-bearing skeletal sites that are most susceptible to bone loss in microgravity, and thus at higher risk for fracture. Given that hip fractures can be debilitating, the initial model development focused on the femoral neck. Future efforts will focus on including other key load bearing bone sites such as the greater trochanter, lower lumbar, proximal femur and calcaneus. The DAP has currently established an initial model (Beta Version) of bone loss due to skeletal unloading in femoral neck region. The model calculates changes in mineralized volume fraction of bone in this segment and relates it to changes in bone mineral density (vBMD) measured by Quantitative Computed Tomography (QCT). The model is governed by equations describing changes in bone volume fraction (BVF), and rates of changes in bone cell populations that remove and replace bone in packets within the bone region. The DAP bone model is unique in several respects. In particular in takes former models of volume fraction changes one step higher in fidelity and separates BVF into separate equations for mineralized and osteoid volume fractions governed by a mineralization rate. This more closely follows the physiology of the remodeling unit cycles where bone is first resorbed and then followed by the action of osteoblasts to lay down collagen matrix which eventually becomes mineralized. In another respect, the modules allow the functional description of the time rate of change of other parameters and variables in the model during a computational simulation. More detailed description of the model, preliminary validation results, current limitation and caveats, and planned advancements are provided in sections 2 through 5. The DAP bone model is being developed primarily as a research tool, and not as a clinical tool like QCT. Even if it transitions to a clinical tool, it is not intended to replace QCT or any other clinical tool. Moreover, the DAP bone model does not predict bone fracture. Its purpose is to provide valuable additional data via "forward prediction" simulations for during and after spaceflight missions to gain insight on, (1) mechanisms of bone demineralization in microgravity, and (2) the volumetric changes at the various bone sites in response to in-flight and post-flight exercise countermeasures. This data can then be used as input to the Keyak [8] (or equivalent) FE analysis method to gain insight on how bone strength may change during and after flight. This information can also be useful to help optimize exercise countermeasure protocols to minimize changes in bone strength during flight, and improve regain of bone strength post-flight. To achieve this goal, the bone model will be integrated with DAP's exercise countermeasure models to simulate the effect of exercise prescriptions on preserving bone. More specifically, the model will accept loading history due to muscle and joint force on bone and produce quantified remodeling within the bone region under influence of the applied stress. Furthermore, because they tend to respond differently, the bone remodeling model includes both trabecular bone and cortical bone.

Pennline, James↗
DOE OSTI2023

Core Modeling & Decision Support Capabilities: FASTSim, RouteE, T3CO & OpenPATH

These are the slides on NREL's Core Modeling and Decision Support Capabilities project for presentation at the 2023 Annual Merit Review for the U.S. Department of Energy, Vehicle Technologies Office. The project is funded by the Energy-Efficient Mobility Systems program to support four core tools: the Future Automotive Systems Technology Simulator (FASTSim), the Route Energy Prediction Model (RouteE), the Transportation Technology Total Cost of Ownership (T3CO) tool, and the Open Platform for Agile Trip Heuristics (OpenPATH).

ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1984

Helicopter model rotor-blade vortex interaction impulsive noise: Scalability and parametric variations

Acoustic data taken in the anechoic Deutsch-Niederlaendischer Windkanal (DNW) have documented the blade vortex interaction (BVI) impulsive noise radiated from a 1/7-scale model main rotor of the AH-1 series helicopter. Averaged model scale data were compared with averaged full scale, inflight acoustic data under similar nondimensional test conditions. At low advance ratios (mu = 0.164 to 0.194), the data scale remarkable well in level and waveform shape, and also duplicate the directivity pattern of BVI impulsive noise. At moderate advance ratios (mu = 0.224 to 0.270), the scaling deteriorates, suggesting that the model scale rotor is not adequately simulating the full scale BVI noise; presently, no proved explanation of this discrepancy exists. Carefully performed parametric variations over a complete matrix of testing conditions have shown that all of the four governing nondimensional parameters - tip Mach number at hover, advance ratio, local inflow ratio, and thrust coefficient - are highly sensitive to BVI noise radiation.

Splettstoesser, W. R.↗
NASA NTRS1985

Helicopter model rotor-blade vortex interaction impulsive noise - Scalability and parametric variations

Acoustic data taken in the anechoic Deutsch-Niederlaendischer Windkanal (DNW) have documented the blade vortex interaction (BVI) impulsive noise radiated from a 1/7-scale model main rotor of the AH-1 series helicopter. Averaged model scale data were compared with averaged full scale, inflight acoustic data under similar nondimensional test conditions. At low advance ratios (mu = 0.164 to 0.194), the data scale remarkable well in level and waveform shape, and also duplicate the directivity pattern of BVI impulsive noise. At moderate advance ratios (mu = 0.224 to 0.270), the scalig deteriorates, suggesting that the model scale rotor is not adequately simulating the full scale BVI noise; presently, no proved explanation of this discrepancy exists. Carefully performed parametric variations over a complete matrix of testing conditions have shown that all of the four governing nondimensional parameters - tip Mach number at hover, advance ratio, local inflow ratio, and thrust coefficient - are highly sensitive to BVI noise radiation.

Boxwell, D. A.↗
NASA NTRS1997

Interannual Variability in Soil Trace Gas (CO2, N2O, NO) Fluxes and Analysis of Controllers

Interannual variability in flux rates of biogenic trace gases must be quantified in order to understand the differences between short-term trends and actual long-term change in biosphere-atmosphere interactions. We simulated interannual patterns (1983-1988) of global trace gas fluxes from soils using the NASA Ames model version of CASA (Carnegie-Ames-Stanford Approach) in a transient simulation mode. This ecosystem model has been recalibrated for simulations driven by satellite vegetation index data from the NOAA Advanced Very High Resolution Radiometer (AVHRR) over the mid-1980s. The predicted interannual pattern of soil heterotropic CO2 emissions indicates that relatively large increases in global carbon flux from soils occurred about three years following the strong El Nino Southern Oscillation (ENSO) event of 1983. Results for the years 1986 and 1987 showed an annual increment of +1 Pg (1015 g) C-CO2 emitted from soils, which tended to dampen the estimated global increase in net ecosystem production with about a two year lag period relative to plant carbon fixation. Zonal discrimination of model results implies that 80-90 percent of the yearly positive increments in soil CO2 emission during 1986-87 were attributable to soil organic matter decomposition in the low-latitudes (between 30 N and 30 S). Soils of the northern middle-latitude zone (between 30 N and 60 N) accounted for the residual of these annual increments. Total annual emissions of nitrogen trace gases (N2O and NO) from soils were estimated to vary from 2-4 percent over the time period modeled, a level of variability which is consistent with predicted interannual fluctuations in global soil CO2 fluxes. Interannual variability of precipitation in tropical and subtropical zones (30 N to 20 S appeared to drive the dynamic inverse relationship between higher annual emissions of NO versus emissions of N2O. Global mean emission rates from natural (heterotrophic) soil sources over the period modeled (1983-1988) were estimated at 57.1 Pg C-CO2yr-1, 9.8Tg (1012 g) N-NO yr-1, and 9.7 Tg N-N2O yr-1. Chemical fertilizer contributions to global soil N gas fluxes were estimated at between 1.3 to 7.3 Tg N-NO yr-1, and 1.2 to 4.0 Tg N-N2O yr-1.

Potter, C.↗
NASA NTRS2005

Anvil Glaciation in a Deep Cumulus Updraught over Florida Simulated with the Explicit Microphysics Model. I: Impact of Various Nucleation Processes

Simulations of a cumulonimbus cloud observed in the Cirrus regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) with an advanced version of the Explicit Microphysics Model (EMM) are presented. The EMM has size-resolved aerosols and predicts the time evolution of sizes, bulk densities and axial ratios of ice particles. Observations by multiple aircraft in the troposphere provide inputs to the model, including observations of the ice nuclei and of the entire size distribution of condensation nuclei. Homogeneous droplet freezing is found to be the source of almost all of the ice crystals in the anvil updraught of this particular model cloud. Most of the simulated droplets that freeze to form anvil crystals appear to be nucleated by activation of aerosols far above cloud base in the interior of the cloud ("secondary" or "in cloud" droplet nucleation). This is partly because primary droplets formed at cloud base are invariably depleted by accretion before they can reach the anvil base in the updraught, which promotes an increase with height of the average supersaturation in the updraught aloft. More than half of these aerosols, activated far above cloud base, are entrained into the updraught of this model cloud from the lateral environment above about 5 km above mean sea level. This confirms the importance of remote sources of atmospheric aerosol for anvil glaciation. Other nucleation processes impinge indirectly upon the anvil glaciation by modifying the concentration of supercooled droplets in the upper levels of the mixed-phase region. For instance, the warm-rain process produces a massive indirect impact on the anvil crystal concentration, because it determines the mass of precipitation forming in the updraught. It competes with homogeneous freezing as a sink for cloud droplets. The effects from turbulent enhancement of the warm-rain process and from the nucleation processes on the anvil ice properties are assessed.

Phillips, Vaughan T. J.↗
DOE OSTI2024

Modeled estimates of permafrost hydrology and related fields for pan-Arctic region over the period 1980-2100

This dataset contains gridded estimates of total runoff (mm/day); subsurface runoff (mm/day); snow water equivalent (SWE, mm/day); vertical upward water flux of combined surface evaporation, evapotranspiration, and snow sublimation (mm/day); and soil temperature (C) across the pan-Arctic basin from simulations of the Permafrost Water Balance Model (PWBM v4). These data were generated to allow researchers to advance understanding of the changes unfolding in the Arctic's terrestrial water cycle under the influences of climate warming, hydrological cycle intensification, and permafrost thaw. The dataset will enable researchers to address important questions spanning multiple disciplines. For example, how will rising river flows exported to the coast impact critical zone ecosystems, near-shore ice dynamics, and ocean biogeochemical processes. Spatial resolution is 25x25 km with N=39926 grid cells spanning the domain. Three forcing datasets were used to drive simulations at an implicit daily time step over the period 1980-2100. Daily air temperature, precipitation, and wind speed data from the WATCH Forcing Data methodology applied to ERA5 reanalysis (W5E5) were used for meteorological forcings over the period 1980-2019. Simulations for the period 1980-2100 used meteorological data from two global climate models; IPSL-CM6A-LR and MPI-ESM1-2-HR. The PWBM hydrological outputs are organized by grid and month in gzipped ASCII text files within each archived tar file. Soil temperature is organized by grid cell and day of year from DOY 205 to DOY 260. No special software is required to work with these data. The package contains spatial domain and river basin information files. These data are most useful for analyses of the spatial and temporal variations in terrestrial hydrology and soil thermal dynamics. A manuscript describing the model simulations and associated analysis has been accepted for publication in The Cryosphere (Rawlins and Karmalkar, 2024).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

Developing Soil Models for Dynamic Impact Simulations

This paper describes fundamental soils characterization work performed at NASA Langley Research Center in support of the Subsonic Rotary Wing (SRW) Aeronautics Program and the Orion Landing System (LS) Advanced Development Program (ADP). LS-DYNA(Registered TradeMark)1 soil impact model development and test-analysis correlation results are presented for: (1) a 38-ft/s vertical drop test of a composite fuselage section, outfitted with four blocks of deployable energy absorbers (DEA), onto sand, and (2) a series of impact tests of a 1/2-scale geometric boilerplate Orion capsule onto soil. In addition, the paper will discuss LS-DYNA contact analysis at the soil/structure interface, methods used to estimate frictional forces, and the sensitivity of the model to density, moisture, and compaction.

Fasanella, Edwin L.↗
NASA NTRSDate not supplied

Digital Lunar Exploration Sites (DLES)

It has been almost 50 years since humans last set foot on the Moon. With NASA’s Artemis program, the United States and its international and commercial partners are embarking on a new endeavor to explore the lunar surface. Before we return, we will have simulated every aspect of these future missions. Many of these models and simulations (M&S) will rely on well-known and commonly-used technologies, some of which trace their origins back to the Apollo program. However, M&S has advanced significantly, as have the underlying computational capabilities. As a result, we are able to model many more aspects of the Artemis vehicles and support systems with significantly improved detail and confidence. Accurately and efficiently modeling the lunar environment will be critical to simulating the Artemis elements and mission activities. This includes characterizing and modeling lunar topography, smaller craters, exposed surface rocks, lunar regolith, surface lighting, and ambient thermal environment. These are all necessary for understanding fundamental behaviors and performance of vehicles and support systems in the lunar environment and are often determining factors in the selection of exploration sites and defining mission profiles. The Astromaterials Research and Exploration Sciences (ARES) team in the Exploration Integration and Science Directorate at NASA’s Johnson Space Center (JSC) and the NASA Exploration Systems Simulations (NExSyS) team in the Simulation and Graphics Branch at JSC are developing and maintaining the Digital Lunar Exploration Sites (DLES) data, documentation, and software packages. This data is being fed directly into a diverse collection of graphics and simulation environments, where it is used to construct the closest known truth for numerous potential Lunar South Pole landing sites. DLES is available to projects across NASA and particularly the Artemis program to support coordinated digital representations of the lunar environment. This paper describes the fundamental need for DLES, the science data sets that are going into DLES, some of the processes used to integrate this data into DLES products, the basic products that constitute DLES, and some examples of DLES in use.

Lunar↗
DOE OSTI2025

Scaling Ultrahigh-Resolution E3SM Land Model for Leadership-Class Supercomputers

This paper presents advancements in scaling the ultrahigh-resolution E3SM Land Model (uELM) for deployment on leadership-class supercomputers, addressing the increased demand for km-scale Earth system modeling. By focusing on km-scale ELM simulations, we enhance predictive capabilities for climate interactions, facilitating improved responses to climate change impacts on energy systems, agriculture, and water resources. Our approach leverages innovative software architecture optimizations, sophisticated data handling techniques, and advanced parallel processing, achieving strong scalability on two leadership supercomputers (2400 nodes (105,600 cores) on Summit, and 1200 nodes (76,800 cores) on Frontier). Results from extensive scalability assessments on the Summit and Frontier also demonstrate outstanding I/O performance (close to 400 GB/s write throughput) and the model's ability to efficiently handle increasing computational demands. This study not only establishes uELM's capability for high-resolution simulations over vast geographical domains, but also sets a foundation for future Earth system modeling breakthroughs.

Wang, Dali [ORNL] (ORCID:0000000168065108)↗
DOE OSTI2024

Celeritas: Accelerating Geant4 with GPUs

Celeritas [1] is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications (specifically, High Lumi- nosity Large Hadron Collider (HL-LHC) simulation) on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a GPU-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.3, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling integration of Celeritas into high energy physics (HEP) experimental frameworks such as CMSSW. On the Summit supercomputer, Celeritas performs EM physics between 6 and 32 faster using the machine’s Nvidia GPUs compared to using only CPUs. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.8–2.3× on GPU and 1.2× on CPU.

Johnson, Seth R.↗