Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Scale down”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2025

Water intensity of photovoltaic module manufacturing at the terawatt scale

As the U.S. ramps photovoltaic (PV) manufacturing to the terawatt scale and emphasizes re-shoring manufacturing, potential regional impacts on the U.S. water supply should be considered, particularly since many PV companies rely almost exclusively on public water supplies for manufacturing. This work surveys the academic literature and PV manufacturer reports to estimate the water intensity of monocrystalline silicon, multicrystalline silicon, and cadmium telluride modules manufactured at the terawatt scale, determining that on average, cadmium telluride manufacturing is less water intensive on a per megawatt scale – this is anticipated to be true for all thin film PV manufacturing. While much lower than the water intensity of thermoelectric (e.g., coal) energy generation, significant issues and gaps with PV manufacturing data quality in academic studies are identified which cause estimates to vary by over 1000x (0.04 – 49 trillion liters/terawatt). Data issues are discussed and the need for accurate accounting of water resources (e.g., via continuous, updated information during PV manufacturing) is highlighted. The opportunity to reconfigure decommissioned thermoelectric sites to PV manufacturing is also explored. Finally, factors that influence PV manufacturing water intensity, from individual manufacturing steps to trends across the PV industry, are examined and water conservation opportunities are presented.

14 SOLAR ENERGY↗
DOE OSTI2023

Developing and tuning a community scale energy model for a disadvantaged community

This work describes the development of a community-scale energy model for a mixed-use low-income community located in Huntington Beach, CA. An accurate community-scale energy model is useful for evaluating the use of limited capital resources used to invest in clean energy technologies. This work lays out the process of developing such a model while relying primarily on publicly available data and highlighting critical partnerships necessary for model development success. The primary contribution of this work is the demonstration of the process used to develop an accurate energy model for a disadvantaged community when minimal building and energy use data is available. The heart of the model is the physics-based community scale energy modeling platform URBANopt. Using a bottom-up load modeling approach, energy simulated energy use falls within 3% or less of aggregate annual utility data, and within 10% or less aggregate monthly utility data. The demonstrated model development and tuning process can be used by others to characterize other atypical communities, which may differ significantly from prototypical models.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Towards intelligent emergency control for large-scale power systems: Convergence of learning, physics, computing and control

Here, this paper has delved into the pressing need for intelligent emergency control in large-scale power systems, which are experiencing significant transformations and are operating closer to their limits with more uncertainties. Learning-based control methods are promising and have shown effectiveness for intelligent power system control. However, when they are applied to large-scale power systems, there are multifaceted challenges such as scalability, adaptiveness, and security posed by the complex power system landscape, which demand comprehensive solutions. The paper first proposes and instantiates a convergence framework for integrating power systems physics, machine learning, advanced computing, and grid control to realize intelligent grid control at a large scale. Our developed methods and platform based on the convergence framework have been applied to a large (more than 3000 buses) Texas power system, and tested with 56 000 scenarios. Our work achieved a 26% reduction in load shedding on average and outperformed existing rule-based control in 99.7% of the test scenarios. The results demonstrated the potential of the proposed convergence framework and DRL-based intelligent control for the future grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Go slow to go fast? A review of the impacts of permitting on large-scale solar project development

State and local permitting challenges could impede the ability of large-scale solar (LSS) to meet growing electricity demand in the United States. Here, we review research that explores LSS permitting and its impacts on the pace and scale of LSS project development. Research on LSS permitting is relatively scarce, such that we support our review with research in the context of wind permitting, where appropriate. Further, few studies attempt to rigorously quantify the effects of permitting on the pace and scale of LSS project development. The available evidence allows us to identify various hypotheses and identify gaps for further research. Our review suggests that differences in permitting policies, regulations, and ordinances explain relatively little variation in LSS permitting and development outcomes across jurisdictions, except where jurisdictions implement rules designed to impede LSS. The evidence suggests LSS permitting challenges largely accrue during the implementation of permitting processes. Recent research suggests that community opposition to project development is a key driver of LSS permitting challenges, given that project opponents often use permitting processes to translate opposition into legal action. We call on future researchers to more concretely describe the LSS permitting challenge and to identify the specific actors responsible for implementing solutions.

Community opposition↗
DOE OSTI2023

Bridging adsorption behavior of confined CH 4 -CO 2 binary mixtures across scales

An accurate understanding of the competitive adsorption of CH 4 -CO 2 binary mixtures in nano-confined systems is critical for engineering CO 2 storage in shale gas reservoirs. Due to difficulties in making reliable experimental observations in nano-scale, atomistic simulations (ASs), such as the Grand Canonical Monte Carlo (GCMC) method, provide a viable approach to studying the adsorption behavior of confined fluids. ASs are, however, limited in the size of the compositional domain due to the high computational cost. This work proposes a framework that combines AS and the lattice Boltzmann (LB) method to bridge the physics of confined fluids across scales. The Peng–Robinson equation of state (PR-EoS) produces fugacity coefficients, which serve as input for conducting multi-component GCMC simulations. These GCMC simulations explore the competitive adsorption behavior of CH 4 -CO 2 in nano-slits at various composition, pressure, and channel-width conditions. Both components generate adsorption layers with high densities near the walls with CO 2 preferentially adsorbing compared to CH 4 on the organic walls of carbon sheets. At the mesoscale, a pseudopotential model represents the intermolecular forces in multi-component, multiple-relaxation-time LB simulations. The LB simulations are in good agreement with the GCMC results, allowing us to obtain values for tunable LB parameters. We then extend the use of LB to simulate adsorption behavior in complex networks with nano-sized channels. The phase behavior and fluid properties in the complex geometries of nano-channels differ from nano-slits and bulk systems. Furthermore, the bridging of physics from GCMC (microscale) to LB (mesoscale) via the macroscale PR-EoS connects the adsorption behavior of binary systems across scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Extrapolation of thermal conductivity in non-equilibrium molecular dynamics simulations to bulk scale

Predictions of the bulk scale thermal conductivity of solids using non-equilibrium molecular dynamics (MD) simulations have relied on the linear extrapolation of the thermal resistivity versus the reciprocal of the system length in the simulations. Several studies have reported deviation of the extrapolation from linearity near the micro-scale, raising a concern of its applicability to large systems. To investigate this issue, present work conducted extensive MD simulations of silicon with two different potentials (EDIP and Tersoff-II) for unprecedented length scales up to 10.3 μm and simulation times up to 530 ns. For large systems ≥0.35 μm in size the non-linearity of the extrapolation of the reciprocal of the thermal conductivity is mostly due to ignoring the dependence of the thermal conductivity on temperature. To account for such dependence, the present analysis fixes the temperature range for determining the gradient for calculating the thermal conductivity values. However, short systems ≤0.23 μm in size show significant non-linearity in the calculated thermal conductivity values using a temperature window of 500 ± 10 K from the simulations results with the EDIP potential. Since these system sizes are shorter than the mean phonon free path in EDIP (~0.22 μm), the nonlinearity may be attributed to phonon transport. For the MD simulations with the Tersoff-II potential there is no significant non-linearity in the calculated thermal conductivity values for systems ranging in size from 0.05 to 5.4 μm.

42 ENGINEERING↗
DOE OSTI2026

Representativity error scaling of models of the high temperature test facility

Error scaling is a critical step toward the validation of modeling capabilities for high-temperature gas-cooled reactors (HTGRs). Extensive effort is being made to bring HTGRs into the validation basis of numerous thermal-hydraulics codes. Here, this paper demonstrates how systems-level codes can be leveraged to perform error scaling analyses between an experimental facility and a plant-to-be facility. Specifically, we focus on two conduction cooldown experiments from the High Temperature Test Facility (HTTF) and the General Atomics 350 MW th modular high-temperature gas-cooled reactor (MHTGR-350). The error scaling methodology employed in this study is representativity, which in the context of this work was used to quantify how well experiments captured the physics of the plant facility by comparing sensitivity vectors between the two facilities. In addition, two different RELAP5–3D models of the HTTF were included in the comparison to determine if modeling methodologies notably impact the representativity results. The key figures of merit are the maximum block temperature and the coolant outlet temperature. The time-dependent maximum block temperature had a low representativity of below 0.1 for both models across both transients. The coolant outlet temperature had a higher representativity of around 0.6 for both models during the pressurized conduction cooldown, but it was below 0.2 for the depressurized conduction cooldown. Overall, the HTTF experiments were not representative of the transients in the MHTGR-350. However, these results can play a significant role in informing future potential experiments for HTGR systems that iterate on what the HTTF accomplished.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Multi-scale model based design of membrane reactor/separator processes for intensified hydrogen production through the water gas shift reaction

This paper aims to first quantify the impact of various diffusion models (Maxwell-Stefan, Wilke, Dusty-Gas) on the predictions of a multi-scale membrane reactor/separator mathematical model, and to then demonstrate this model's use for the design and process intensification of membrane reactor/separator systems for hydrogen production. This multi-scale model captures velocity, temperature and species' concentration profiles along the catalyst pellet's radial direction, and along the reactor's axial direction, by solving the momentum, energy, and species transport equations, accounting for convection, conduction, reaction, and diffusion mechanisms. In the first part of work, the effect of pellet-scale design parameters (mean pore diameter, volumetric porosity, tortuosity factor, etc.) and various species' flux models on the model predictions is studied. In the second part, the study focuses on the comparison, in terms of their process intensification characteristics, of various hydrogen production processes. These include a conventional high-temperature shift reactor (HTSR)/low-temperature shift reactor (LTSR) sequence, a novel HTSR/membrane separator (MS)/LTSR/MS sequence, and a process that involves low-temperature shift membrane reactors-LTSMR in a series.

08 HYDROGEN↗
DOE OSTI2022

Detection and attribution of long-term and fine-scale changes in spring phenology over urban areas: A case study in New York State

Spring phenology plays an essential role in climate change, terrestrial ecosystem, and public health. Field-based monitoring and understanding of changes in spring phenology for long periods and in large regions are challenging due to the limited in-site observations. Space-based remotely sensed observations offer great potentials for monitoring decadal spring phenology changes from regional to global scales. However, the coarse-scale remotely sensed observations are insufficient to capture fine-scale spring phenology dynamics, especially in urban areas, and this makes it challenging for understanding the combined effects of climate change and urbanization on spring phenology. We derived the start of phenology season (SOS) in New York State using 30 m Landsat observations from 1990 to 2015 to understand the impact of the environment and urbanization on SOS. The results show that SOS for different years reveals heterogeneous spatial distribution. Most regions of New York State have been experiencing significant spring phenology changes in form of earlier onset of vegetation greening, ranging from 0.2 to 0.6 day/year during 1990 to 2015, and this trend varies slightly with latitudes and urbanization levels. Further, spatial correlation analysis shows that the increase in temperature and urbanization could both promote the advancement of SOS. However, the effect of urbanization (partial correlation coefficient (R) ranges from −0.289 to −0.542) on SOS is greater than the effect of temperature (R ranges from 0.006 to −0.192). The study generates a high spatio-temporal resolution spring phenology dataset for ecological, environmental and public health studies, especially in urban areas, and reveals the importance of better accounting for the urbanization effects when quantifying the SOS dynamics in phenology models.

Landsat↗
DOE OSTI2024

Convergence of variational Monte Carlo simulation and scale-invariant pre-training

We provide theoretical convergence bounds for the variational Monte Carlo (VMC) method as applied to optimize neural network wave functions for the electronic structure problem. Here, we study both the energy minimization phase and the supervised pre-training phase that is commonly used prior to energy minimization. For the energy minimization phase, the standard algorithm is scale-invariant by design, and we provide a proof of convergence for this algorithm without modifications. The pre-training stage typically does not feature such scale-invariance. We propose using a scale-invariant loss for the pretraining phase and demonstrate empirically that it leads to faster pre-training.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2023

Core-scale numerical simulation and comparison of breakdown of shale and resulting fractures using sc-CO 2 and water as injectants

Supercritical carbon dioxide (sc-CO 2 ) is an alternative to water for stimulation of low permeability systems such as shale gas and geothermal resources. Previously core-scale experimental studies have compared the behavior of CO 2 to water injection for sample breakdown. Due to differences in experimental setup and core sample preparation, inconsistent or even apparently contradictory conclusions have resulted. To reconcile this contradiction, a phase-field numerical model is applied to understand hydraulic fracturing experiments using Green River shale found in the literature. The finite element numerical model incorporates a rate-dependent phase-field fracture model developed separately to describe fracture initiation and growth. We investigate the impact of various material and fluid properties on the resulting fractures. Most importantly, we study the effect of fluid properties and boundary conditions on the breakdown pressure, including the direction of the resulting fracture plane. Model results predict that (1) sc-CO 2 injection in the laboratory may result in greater breakdown pressure than that of water under no-flow boundary conditions because lower viscosity sc-CO 2 may result in pressure build up at the core boundary that opposes fracture initiation and (2) lower viscosity sc-CO 2 also produces fast-propagating fractures that are less influenced by the bedding plane on their resulting fracture topology. Here our model offers a straightforward explanation and reconciliation of existing experimental observations, as well as a means to extrapolate to new conditions. Exploration of field-scale conditions suggests less pronounced or no elevation in breakdown pressure when sc-CO 2 is injected because the pressure build up effect at the system boundary is significantly less or absent at field length scales.

42 ENGINEERING↗
DOE OSTI2023

Dynamics of streamflow permanence in a headwater network: Insights from catchment-scale model simulations

The hillslope and channel dynamics that govern streamflow permanence in headwater systems have important implications for ecosystem functioning and downstream water quality. Recent advancements in process-based, semi-distributed hydrologic models that build upon empirical studies of streamflow permanence in well-monitored headwater catchments show promise for characterizing the dynamics of streamflow permanence in headwater systems. However, few process-based models consider the continuum of hillslope-stream network connectivity as a control on streamflow permanence in headwater systems. The objective of this study was to expand a process-based, catchment-scale hydrologic model to better understand the spatiotemporal dynamics of headwater streamflow permanence and to identify controls of streamflow expansion and contraction in a headwater network. Further, we aimed to develop an approach that enhanced the fidelity of model simulations, yet required little additional data, with the intent that the model might be later transferred to catchments with limited long-term and spatially explicit measurements. This approach facilitated network-scale estimates of the controls of streamflow expansion and contraction, albeit with higher degrees of uncertainty in individual reaches due to data constraints. Our model simulated that streamflow permanence was highly dynamic in first-order reaches with steep slopes and variable contributing areas. The simulated stream network length ranged from nearly 98±2% of the geomorphic channel extent during wet periods to nearly 50±10% during dry periods. The model identified a discharge threshold of approximately 1 mm d –1 , above which the rate of streamflow expansion decreases by nearly an order of magnitude, indicating a lack of sensitivity of streamflow expansion to hydrologic forcing during high-flow periods. Overall, we demonstrate that process-based, catchment-scale models offer important insights on the controls of streamflow permanence, despite uncertainties and limitations of the model. We encourage researchers to increase data collection efforts and develop benchmarks to better evaluate such models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessing mechanical response of CO 2 storage into a depleted carbonate reef using a site-scale geomechanical model calibrated with field tests and InSAR monitoring data

Geomechanical risks of injection have raised concerns regarding secure CO 2 storage. In this work, a combined monitoring and modeling approach is used to assess the stress changes and surface uplift associated with CO 2 injection into a depleted carbonate reef of the Michigan basin. A site-scale geomechanical model is built by assigning mechanical properties of formations using well-log and experimental data. Gravity load is applied to the model to estimate the vertical component of stress as well as different lateral boundary displacement scenarios to estimate horizontal stresses. We used a poroelastic pressure-dependent model (instead of a linear elastic mechanical earth model) to calibrate initial stresses using hydraulic fracture test data measured at depleted reservoir status. Multi-phase fluid flow-geomechanical simulations are performed to estimate the poroelastic response during (1) primary depletion (2) field-scale CO 2 injection phase (3) a hypothetical forecast scenario in which well bottom hole pressure (BHP) reach 45000 KPa. The predicted surface uplift is less than 1 mm at the end of the field-scale CO 2 injection phase which is in good agreement with Interferometric Synthetic Aperture Radar (InSAR) uplift measurement. Although the InSAR data shows an insignificant uplift, hydromechanical modeling of injection shows that CO 2 injection still causes reservoir deformation emphasizing the role of carbonate overburden and reservoir formation mechanical properties and limited size of reef on diminishing the surface deformation. Modeling indicates poroelastic response of caprock matters to estimate uplift. The lower permeability of the top two layers provides additional barrier to large uplift. Also, history of subsidence due to production should be accounted to predict uplift due to a follow up injection correctly. This report shows the significance of combining a calibrated geomechanical model with field measured stresses and monitoring data to be used as a tool to ensure the safety of CO 2 storage.

42 ENGINEERING↗
DOE OSTI2021

Development of small-scale plutonium electrorefining in molten CaCl 2

Molten salt extraction and electrorefining of plutonium were successfully conducted with only a 200-g initial charge. Due to this small scale, the electrorefining cell components were redesigned relative to conventional, multi-kg standards. The efficiency and yield of this small-scale process are comparable to those of larger-scale exemplars. Chemical analysis of the plutonium product and the anode heel is presented. Here, the product was 99.8% pure. The anode heel was able to be depleted beyond the expected freezing point without significant impurity carryover.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

A linked-scale coupled model of mass erosion and redistribution in plasma-exposed micro-foam surfaces

Surface evolution due to exposure to harsh environments is of importance in many scientific and technological applications. In plasma-exposed materials, the surface receives charged particles from the plasma, leading to a series of processes that drive the system far from equilibrium and may lead to the severe deterioration of the surface properties. Although surface morphological changes are driven by atomic collisions taking place over picoseconds and nanometers, these processes result in mass loss and redistribution of matter over much larger length and time scales. This necessitates a multi-scale approach capable of capturing the range of processes linking primary atomic collision events with engineering-level surface geometry changes. In this paper, we develop a computational model to simulate the morphological evolution and effective erosion rate of micro-architected tungsten foams during low-energy plasma ion bombardment. Furthermore, the model acts on several length scales, with the energy and angular dependence of the sputtering yield of flat tungsten surfaces determined using the SRIM code based on the binary collision approximation. This information is introduced into a low-fluence, short-term Monte Carlo raytracing model, and further into a high-fluence, long-term particle transport model. In the latter, material particles representing billions of atoms are described in a digitized 3-D representation of the foam structure from X-ray tomography data, and are sputtered off and redeposited using the atomistic information. We show that the redeposition of sputtered atoms leads to partial self-healing in the bottom layers with a sharp reduction in the sputtering coefficient of low density foams, roughly 25% of the solid W value. This is in qualitative agreement with recent experiments on low porosity W structures. At high fluence, the foam structure degrades considerably as there are fewer ligaments available to recapture sputtered atoms and, consequently, the sputtering rate increases again.

36 MATERIALS SCIENCE↗