Search NASA⌕ Search

SEARCH · Search NASA

Results for “process 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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cooperative Research and Development Agreement among National Energy Technology Laboratory, The Board of Regents of the University System of Georgia by and on behalf of Georgia Institute of Technology, and University of Texas at El Paso [Abstract]

Georgia Institute of Technology (GIT), University of Texas at El Paso (UTEP), and the National Energy Technology Laboratory (NETL) will collaborate and develop cyber-physical systems for advanced power systems. GIT will test virtual models on hybrid power system hardware and cyber physical components. UTEP will test novel actuators on hybrid power system hardware. NETL will access GIT’s expertise in the engineering and real time simulation of power system processes for cyber-physical systems and simulations to control hardware. NETL will establish a relationship to promote UTEP’s capabilities in additive manufacturing and high-temperature piezo-electric devices with the potential for controlling high-temperature flows with a light and inexpensive device at high response rates. This project is the result of an EERE Advanced Manufacturing Office open Multi-Topic awarded proposal led by GIT titled, “Robust Combined Heat and Hybrid Power (CHHP) for High Electrical Efficiency Cogeneration”.

42 ENGINEERING↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

DOE ASR Final Report on “Use of ARM Observations to Investigate the Role of Tropical Radiative Processes and Cloud Radiative Effects in Climate Simulations”

The overall objective of this ASR funded project is to investigate the role of cloud radiative effects, especially those associated with tropical thin cirrus clouds in the tropical tropopause layer, by analyzing the ARM observations combined with numerical models. In particular, we have processed and analyzed the observations from the Raman lidar at the ARM SGP and TWP sites. In the tenure of the project (8/15/2013 – 8/14/2016 and with a no-cost extension to 8/14/2017), we have been concentrating on (i) developing an automated feature detection scheme of clouds and aerosols for the ARM Raman lidar; (ii) developing an automated retrieval of cloud and aerosol extinctions for the ARM Raman lidar; (iii) investigating cloud radiative effects based on the observations on the simulated temperatures in the tropical tropopause layer using a radiative-convective model; and (iv) examining the effect of changes of atmospheric composition on the tropical lower-stratospheric temperatures. In addition, we have examined the biases in the CALIPSO-inferred aerosol direct radiative effects using ground-based Raman lidars at the ARM SGP and TWP sites, and estimated the impact of lidar detection sensitivity on assessing global aerosol direct radiative effects. We have also investigated the diurnal cycle of clouds and precipitation at the ARM site using the cloud radar observations along with simulations from the multiscale modeling framework. The main results of our research efforts are reported in the six referred journal publications that acknowledge the DOE Grant DE-SC0010557.

47 OTHER INSTRUMENTATION↗
DOE OSTI2011

Characterization of Settler Tank, KW Container and KE Container Sludge Simulants

The Sludge Treatment Project (STP), managed by CH2M Hill Plateau Remediation Company (CHPRC) has specified base formulations for non-radioactive sludge simulants for use in the development and testing of equipment for sludge sampling, retrieval, transport, and processing. In general, the simulant formulations are based on the average or design-basis physical and chemical properties obtained by characterizing sludge samples. The simulants include surrogates for uranium metal, uranium oxides (agglomerates and fine particulate), and the predominant chemical phases (iron and aluminum hydroxides, sand). Specific surrogate components were selected to match the nominal particle-size distribution and particle-density data obtained from sludge sample analysis. Under contract to CHPRC, Pacific Northwest National Laboratory (PNNL) has performed physical and rheological characterization of simulants, and the results are reported here. Two base simulant types (dry) were prepared by STP staff at the Maintenance and Storage Facility and received by PNNL in February 2009: Settler Tank Simulant and KW Container Sludge Simulant. A third simulant, KE Container Sludge Simulant was received by PNNL in December 2010. The objectives of this simulant characterization effort were to provide baseline characterization data on simulants being used by STP for process development and equipment testing and provide a high-level comparison of the simulant characteristics to the targets used to formulate the simulants.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Comparing quantile regression forest and mixture density long short-term memory models for probabilistic post-processing of satellite precipitation-driven streamflow simulations

Abstract. Deep learning (DL) and machine learning (ML) are widely used in hydrological modelling, which plays a critical role in improving the accuracy of hydrological predictions. However, the trade-off between model performance and computational cost has always been a challenge for hydrologists when selecting a suitable model, particularly for probabilistic post-processing with large ensemble members. This study aims to systematically compare the quantile regression forest (QRF) model and countable mixtures of asymmetric Laplacians long short-term memory (CMAL-LSTM) model as hydrological probabilistic post-processors. Specifically, we evaluate their ability in dealing with biased streamflow simulations driven by three satellite precipitation products across 522 nested sub-basins of the Yalong River basin in China. Model performance is comprehensively assessed using a series of scoring metrics from both probabilistic and deterministic perspectives. Our results show that the QRF model and the CMAL-LSTM model are comparable in terms of probabilistic prediction, and their performances are closely related to the flow accumulation area (FAA) of the sub-basin. The QRF model outperforms the CMAL-LSTM model in most sub-basins with smaller FAA, while the CMAL-LSTM model has an undebatable advantage in sub-basins with FAA larger than 60 000 km2 in the Yalong River basin. In terms of deterministic predictions, the CMAL-LSTM model is preferred, especially when the raw streamflow is poorly simulated and used as input. However, setting aside the differences in model performance, the QRF model with 100-member quantiles demonstrates a noteworthy advantage by exhibiting a 50 % reduction in computation time compared to the CMAL-LSTM model with the same ensemble members in all experiments. As a result, this study provides insights into model selection in hydrological post-processing and the trade-offs between model performance and computational efficiency. The findings highlight the importance of considering the specific application scenario, such as the catchment size and the required accuracy level, when selecting a suitable model for hydrological post-processing.

Geology↗
DOE OSTI2025

Flexible User-Defined Domain Decomposition in Kilometer-Scale E3SM Land Model Simulation

The Energy Exascale Earth System Model (E3SM) Land Model (ELM) has been extended to kilometer-scale (km-ELM) resolutions, enabling high-fidelity simulations of terrestrial processes at 1 km x 1 km grid spacing. In ELM, domain decomposition partitions the computational domain across processors, ensuring efficient parallel execution. Currently, round-robin decomposition is applied, providing a straightforward way to distribute computational workload. As ELM continues evolving at the kilometer-scale (km-scale), particularly with integrating lateral flow modeling, decomposition strategies must also account for the increased workload and data movement. This paper introduces a flexible user-defined domain decomposition framework, allowing users to customize domain partitioning based on application requirements. The impact of different decomposition strategies is evaluated across various applications concerning computation, communication, and I/O. Results demonstrate that while 1D partitioning yields superior I/O performance, k-nearest neighbors (KNN) clustering effectively reduces inter-process communication overhead. This study lays the groundwork for scalable partitioning in large-scale land surface simulations, enhancing next-generation Earth system modeling.

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

Compression behavior of microcrystalline cellulose spheres: Single particle compression and confined bulk compression across regimes

Particle characteristics can drastically influence the process-structure-property-performance aspects of granular materials in compression. In this study, we aim to computationally simulate the mechanical processes of stress redistribution in compacts including the kinematics of particle rearrangement during densification and particle deformation leading to fragmentation. Confined compression experiments are conducted with three sets of commercial microcrystalline cellulose particles nearly spherical in shape with different mean particle size. Experimentally measured compression curves from tall powder columns are fitted with the Kenkre et al. (J. of American Chemical Society, Vol. 79, No. 12) model. This model provides a basis to derive several common two-parameter literature models and as a framework to incorporate statistical representations of critical particle behaviors. We focus on the low-stress compression data and the model comparisons typically not discussed in the literature. Furthermore, additional single particle compressions report fracture strength with particle size for comparison to the apparent particle strength extracted from bulk compression data.

42 ENGINEERING↗
DOE OSTI2025

On-lattice kinetic Monte Carlo approaches for modeling molecular anisotropy in resveratrol crystallization

Stilbenes are a class of organic compounds with broad-ranging pharmaceutical and agricultural applications, which are typically isolated and purified through recrystallization. We are motivated by reducing experimental waste and optimizing yield via developing predictive simulations for processing-dependent crystal morphologies. Using resveratrol as a model stilbene system, we have developed an approach for simulating crystallization with molecular resolution using on-lattice kinetic Monte Carlo. In this work, we highlight modifications to the Stochastic Parallel PARticle Kinetic Simulator (SPPARKS) software package, which were essential to this application. Key enhancements include the incorporation of non-orthogonal cell shapes and monomer anisotropy approximations using bound hard spheres. This new SPPARKS application has been applied to resveratrol with attachment energy libraries obtained from density functional theory, resulting in excellent agreement with experimental morphology prediction.

crystallization↗
DOE OSTI2013

Experimental Studies of Elementary Particle Interactions at High Energies

This is the final report of a program of research on "Experimental Studies of Elementary Particle Interactions at High Energies'' of the High Energy Physics (HEP) group of The Rockefeller University. The research was carried out using the Collider Detector at Fermilab (CDF) and the Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC) at CERN. Three faculty members, two research associates, and two postdoctoral associates participated in this project. At CDF, we studied proton-antiproton collisions at an energy of 1.96 TeV. We focused on diffractive interactions, in which the colliding antiproton loses a small fraction of its momentum, typically less than 1%, while the proton is excited into a high mass state retaining its quantum numbers. The study of such collisions provides insight into the nature of the diffractive exchange, conventionally referred to as Pomeron exchange. In studies of W and Z production, we found results that point to a QCD-based interpretation of the diffractive exchange, as predicted in a data-driven phenomenology developed within the Rockefeller HEP group. At CMS, we worked on diffraction, supersymmetry (SUSY), dark matter, large extra dimensions, and statistical applications to data analysis projects. In diffraction, we extended our CDF studies to higher energies working on two fronts: measurement of the single/double diffraction and of the rapidity gap cross sections at 7 TeV, and development of a simulation of diffractive processes along the lines of our successful model used at CDF. Working with the PYTHIA8 Monte Carlo simulation authors, we implemented our model as a PYTHIA8-MBR option in PYTHIA8 and used it in our data analysis. Preliminary results indicate good agreement. We searched for SUSY by measuring parameters in the Constrained Minimal Supersymmetric extension of the Standard Model (CMSSM) and found results which, combined with other experimental constraints and theoretical considerations, indicate that the CMSSM is not a viable model. Expressing our results in terms of simple topologies, we exclude squark masses below 0.75 TeV and gluino masses below 1.1 TeV. Astrophysical measurements suggest that about 80% of the matter density of the Universe is non-luminous. One of the theories on dark matter attributes it to Weakly Interacting Massive Particles (WIMPs). We searched for WIMPs in 7 TeV and 8 TeV collisions at CMS and set limits on WIMP production rates, which are competitive and complementary to those of direct detection experiments. Searching for monojets (events with only one jet), which in a popular model could be produced by a jet paired by a gravitino that escapes into extra dimensions, we significantly improved the previously set limit. Our results have been used to set limits on Higgs decay to invisible particles and on production of top squarks in compressed SUSY scenarios. Statistics. We computed Bayesian reference priors for several types of measurement and used them in the analysis of CMS data; investigated the applicability of bootstrap methods to HEP measurements; studied several issues associated with simple-versus-simple hypothesis testing and applied the resulting methods to the measurement of some properties of the top quark and Higgs boson.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Co-scheduling Ensembles of In Situ Workflows

Molecular dynamics (MD) simulations are widely used to study large-scale molecular systems. HPC systems are ideal platforms to run these studies, however, reaching the necessary simulation timescale to detect rare processes is challenging, even with modern supercomputers. To overcome the timescale limitation, the simulation of a long MD trajectory is replaced by multiple short-range simulations that are executed simultaneously in an ensemble of simulations. Analyses are usually co-scheduled with these simulations to efficiently process large volumes of data generated by the simulations at runtime, thanks to in situ techniques. Executing a workflow ensemble of simulations and their in situ analyses requires efficient co- scheduling strategies and sophisticated management of computational resources so that they are not slowing down each other. In this paper, we propose an efficient method to co-schedule simulations and in situ analyses such that the makespan of the workflow ensemble is minimized. We present a novel approach to allocate resources for a workflow ensemble under resource constraints by using a theoretical framework modeling the workflow ensemble’s execution. We evaluate the proposed approach using an accurate simulator based on the WRENCH simulation framework on various workflow ensemble configurations. Results demonstrate the significance of co-scheduling simulations and in situ analyses that couple data together to benefit from data locality, in which inefficient scheduling decisions can lead to slowdown in makespan up to a factor of 30.

Do, Tu Mai Anh↗
DOE OSTI2016

Investigation of boundary-layer wind predictions during nocturnal low-level jet events using the Weather Research and Forecasting model

Simulations of two periods featuring three consecutive low level jet (LLJ) events in the US Upper Great Plains during the autumn of 2011 were conducted to explore the impacts of various setup configurations and physical process models on simulated flow parameters within the lowest 200 m above the surface, using the Weather Research and Forecasting (WRF) model. Sensitivities of simulated flow parameters to the horizontal and vertical grid spacing, planetary boundary layer (PBL) and land surface model (LSM) physics options, were assessed. Data from a Light Detection and Ranging (lidar) system, deployed to the Weather Forecast Improvement Project (WFIP; Finley et al. 2013) were used to evaluate the accuracy of simulated wind speed and direction at 80 m above the surface, as well as their vertical distributions between 120 and 40 m, covering the typical span of contemporary tall wind turbines. All of the simulations qualitatively captured the overall diurnal cycle of wind speed and stratification, producing LLJs during each overnight period, however large discrepancies occurred at certain times for each simulation in relation to the observations. 54-member ensembles encompassing changes of the above discussed configuration parameters displayed a wide range of simulated vertical distributions of wind speed and direction, and potential temperature, reflecting highly variable representations of stratification during the weakly stable overnight conditions. Root mean square error (RMSE) statistics show that different ensemble members performed better and worse in various simulated parameters at different times, with no clearly superior configuration . Simulations using a PBL parameterization designed specifically for the stable conditions investigated herein provided superior overall simulations of wind speed at 80 m, demonstrating the efficacy of targeting improvements of physical process models in areas of known deficiencies. However, the considerable magnitudes of the RMSE values of even the best performing simulations indicate ample opportunities for further improvements.

Mirocha, Jeff D.↗
DOE OSTI2020

Coupling solidification kinetics with phase-behavior computations in hydrodynamic simulations of high-pressure, dynamic-compression processes

In this study, we report a numerical scheme to integrate models for the kinetics of solidification processes together with phase-behavior computations in the context of continuum-scale hydrodynamic simulations. The objective of the phase-behavior computations is to determine the pressure and temperature, given the following three sets of inputs: (1) an appropriate equation of state to describe our system, (2) the phase fraction(s) produced by the kinetic models, (3) and the volume and internal energy obtained by solving the conservation equations that govern the hydrodynamic behavior. The kinetics are assumed to be governed by the Kolmogorov–Johnson–Mehl–Avrami equation, and the nucleation and growth rates that enter into that equation are functions of the pressure and temperature produced by the phase-behavior computations. Our formulation allows for the fluid and solid phases to be at different temperatures (thermal nonequilibrium) and pressures (arising from surface-tension-induced Laplace contributions). The formulation is presented in a fairly general setting that is independent of any particular material, although we demonstrate it in some examples related to high-energy-density science applications where materials are rapidly compressed to pressures exceeding several gigapascals in less than a microsecond. We conclude with a critical evaluation of our approach and provide suggestions for future work to improve the predictive capabilities and generality of the models.

42 ENGINEERING↗
DOE OSTI2011

Chromospheric anemone jets and magnetic reconnection in partially ionized solar atmosphere

The solar optical telescope onboard Hinode with temporal resolution of less than 5 s and spatial resolution of 150 km has observed the lower solar atmosphere with an unprecedented detail. This has led to many important findings, one of them is the discovery of chromospheric anemone jets in the solar chromosphere. The chromospheric anemone jets are ubiquitous in solar chromosphere and statistical studies show that the typical length, life time and energy of the chromospheric anemone jets are much smaller than the coronal events (e.g., jets/flares/CMEs). Among various observational parameters, the apparent length and maximum velocity shows good correlation. The velocity of chromospheric anemone jets is comparable to the local Alfven speed in the lower solar chromosphere. Since the discovery of chromospheric anemone jets by Hinode, several evidences of magnetic reconnection in chromospheric anemone jets have been found and these observations are summarized in this paper. These observations clearly suggest that reconnection occurs quite rapidly as well as intermittently in the solar chromosphere. In the solar corona ({lambda}{sub i} > {delta}{sub SP}), anomalous resistivity arises due to various collisionless processes. Previous MHD simulations show that reconnection becomes fast as well as strongly time-dependent due to anomalous resistivity. Such processes would not arise in the solar chromosphere which is fully collisional and partially-ionized. So, it is unclear how the rapid and strongly time-dependent reconnection would occur in the solar chromosphere. It is quite likely that the Hall and ambipolar diffusion are present in the solar chromosphere and they could play an important role in driving such rapid, strongly time-dependent reconnection in the solar chromosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Outcomes from a Collaborative Approach to a Code Comparison Study for Enhanced Geothermal Systems

A code comparison study, supported by the U.S. Department of Energy, Geothermal Technologies Office, has recently been completed whose stated purpose was to test, diagnose differences, and demonstrate the state of modeling capabilities developed at national laboratories and academic institutes within the United States for modeling enhanced geothermal systems. This three-year study comprised two sets of problems; those classified as benchmark problems and those classified as challenge problems. The benchmark problems were structured to test the ability of the collection of numerical simulators to solve various combinations of coupled thermal, hydrologic, geomechanical, and geochemical processes. This class of problems was strictly defined in terms of properties, driving forces, initial conditions, and boundary conditions. Study participants submitted solutions to problems for which their simulation tools were deemed capable or nearly capable. Some participating codes were originally developed for EGS applications whereas some others were designed for different applications but can simulate processes similar to those in EGS. Solution submissions from both were encouraged. In some cases participants made small incremental changes to their numerical simulation codes to address specific elements of the problem, and in other cases participants submitted solutions with existing simulation tools, acknowledging the limitations of the code. The challenge problems were based on the enhanced geothermal systems research conducted in hot dry rock at Fenton Hill, near Los Alamos, New Mexico, between 1974 and 1995, covering two research stimulation, development and circulation phases in two separate reservoirs. Both challenge problems have specific questions to be answered via numerical simulation in three topical areas: 1) reservoir creation/stimulation, 2) reactive and passive transport, and 3) thermal recovery. Whereas the benchmark class of problems were designed to test capabilities for modeling coupled processes under strictly specified conditions, the stated objective for the challenge class of problems was to demonstrate what new understanding of the Fenton Hill experiments could be realized via the application of modern numerical simulation tools by recognized expert practitioners. Critical observations and data from the Fenton Hill experiments were varied and scattered among a number of individual experiments, and numerical simulation solutions were sought that satisfied as many observations at possible. Achieving agreement in multiple experimental observations from Fenton Hill often required participants to critically think about stimulation mechanisms involving natural and hydraulic fractures and re-evaluate their conceptual models and numerical solution approaches. This process has yielded new insights to Fenton Hill reservoirs and direction for future EGS research. The keys to the success of this study were the high quality numerical simulation tools available, the expertise of the developers and practitioners, and the collaborative approach to reporting preliminary results among the study participants. This paper describes results submitted to the challenge problems and the outcomes of a collaborative approach to conducting a code comparison study.

Enhanced Geothermal Systems, numerical simulation,↗
DOE OSTI2026

Evaluating Energy Savings Measures for Building Retrofits in an Interactive Workshop Environment

As part of the Better Buildings Summit 2025, a workshop was developed to simulate the decision process for implementing a set of energy savings measures on an existing building. As part of the workshop a "game" was created around energy savings strategies, and workshop participants selected renovation cards that provided data around the impacts of that renovation choice. Participants then were able to calculate overall savings and payback from the selected retrofit opportunities. This publication has been produced to package the resources used in the activity in an organized document so that the materials can be distributed and reused more broadly. It provides a comprehensive framework for replicating the workshop, including step-by-step instructions, retrofit "game cards" with data points, and facilitation guidance. By using these materials, readers can conduct their own version of this interactive workshop with stakeholders in their networks, such as building owners, energy managers, or policy makers. The goal is to empower users to simulate the decision-making processes for energy projects, foster collaboration, and drive meaningful discussions around cost-effective retrofit opportunities and goal-oriented outcomes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The LSST DESC DC2 Simulated Sky Survey

Here, we describe the simulated sky survey underlying the second data challenge (DC2) carried out in preparation for analysis of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) by the LSST Dark Energy Science Collaboration (LSST DESC). Significant connections across multiple science domains will be a hallmark of LSST; the DC2 program represents a unique modeling effort that stresses this interconnectivity in a way that has not been attempted before. This effort encompasses a full end-to-end approach: starting from a large N-body simulation, through setting up LSST-like observations including realistic cadences, through image simulations, and finally processing with Rubin's LSST Science Pipelines. This last step ensures that we generate data products resembling those to be delivered by the Rubin Observatory as closely as is currently possible. The simulated DC2 sky survey covers six optical bands in a wide-fast-deep area of approximately 300 deg 2 , as well as a deep drilling field of approximately 1 deg 2 . We simulate 5 yr of the planned 10 yr survey. The DC2 sky survey has multiple purposes. First, the LSST DESC working groups can use the data set to develop a range of DESC analysis pipelines to prepare for the advent of actual data. Second, it serves as a realistic test bed for the image processing software under development for LSST by the Rubin Observatory. In particular, simulated data provide a controlled way to investigate certain image-level systematic effects. Finally, the DC2 sky survey enables the exploration of new scientific ideas in both static and time domain cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A GPU Accelerated Mixed‐Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems

In nature, many complex multi‐physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU‐accelerated, mixed‐precision configuration for the FDiRW solver to maximize efficiency through GPU multi‐threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU‐accelerated mixed‐precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU‐accelerated mixed‐precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium‐sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Moist thermodynamics of the Madden-Julian Oscillation in a cloud resolving simulation

The moist thermodynamic processes that determine the scale and energy of the Madden-Julian Oscillation are investigated using moisture and eddy available potential energy (EAPE) budget analyses on a high resolution regional model simulation. The two MJO episodes observed during the winter of 2007-2008 are realistically simulated. In the model, small differences among the timescales of convective vertical transport, mixing and condensation of moisture determine the timescale of MJO moistening. Furthermore, various cloud types play a damping or destabilizing contribution role in the EAPE budget of the MJO, depending on their characteristic latent heating profile and its covariance with the temperature fluctuations. The results are used identify possible sources of the difficulties in simulating MJO in low resolution models that rely on cumulus parameterizations.

54 ENVIRONMENTAL SCIENCES↗