Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time step”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2023

The role of stiffness in training and generalization of ResNets

Neural ordinary differential equations (NODEs) have recently regained popularity as large-depth limits of a large class of neural networks. In particular, residual neural networks (ResNets) are equivalent to an explicit Euler discretization of an underlying NODE, where the transition from one layer to the next is one time step of the discretization. The relationship between continuous and discrete neural networks has been of particular interest. Notably, analysis from the ordinary differential equation viewpoint can potentially lead to new insights for understanding the behavior of neural networks in general. In this work, we take inspiration from differential equations to define the concept of stiffness for a ResNet via the interpretation of a ResNet as the discretization of a NODE. Here, we then examine the effects of stiffness on the ability of a ResNet to generalize, via computational studies on example problems coming from climate and chemistry models. We find that penalizing stiffness does have a unique regularizing effect, but we see no benefit to penalizing stiffness over L 2 regularization (penalization of network parameter norms) in terms of predictive performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MLMOD: Machine Learning Methods for Data-Driven Modeling in LAMMPS

MLMOD is a software package for incorporating machine learning approaches and models into simulations of microscale mechanics and molecular dynamics in LAMMPS. Recent machine learning approaches provide promising data-driven approaches for learning representations for system behaviors from experimental data and high fidelity simulations. The package facilitates learning and using data-driven models for (i) dynamics of the system at larger spatial-temporal scales (ii) interactions between system components, (iii) features yielding coarser degrees of freedom, and (iv) features for new quantities of interest characterizing system behaviors. MLMOD provides hooks in LAMMPS for (i) modeling dynamics and time-step integration, (ii) modeling interactions, and (iii) computing quantities of interest characterizing system states. The package allows for use of machine learning methods with general model classes including Neural Networks, Gaussian Process Regression, Kernel Models, and other approaches. Here we discuss our prototype C++/Python package, aims, and example usage. For related papers, examples, updates, and additional information see https://github.com/atzberg/mlmod and http://atzberger.org/.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Hybrid Natural Gas Geothermal Combined Cycle Power Plant Analysis

Low temperature geothermal resources, including those associated with oil and gas production, are an underutilized source of low carbon energy. The present work investigates coupling of low-temperature geothermal resources with concentrating solar and/or natural gas energy sources to increase the number of locations at which power generation from low temperature geothermal resources would be technically and economically viable. Stand-alone and hybrid geothermal power cycles are simulated using SimTech IPSEpro process modeling software. Design point strategies for a hybrid power cycle that may operate with either a single heat source or two simultaneous heat input sources are considered. Additionally, off-design power plant operation is investigated to consider the impacts of the heat source availability and ambient temperature variations. The off-design analysis utilizes a modeling tool that predicts power plant performance at each time step as a function of topping cycle heat input (from gas or solar), bottoming cycle heat input (from geothermal), and ambient temperature. Addition of a steam topping cycle to an organic Rankine cycle geothermal power plant provides opportunities to increase the efficiency and power output relative to a stand-alone geothermal power plant. Additionally, use of the waste heat from gas turbine power generation in a geothermal bottoming cycle provides opportunities to increase the amount of power generation associated with each unit of carbon dioxide emitted. This paper will describe the hybrid plant configuration evaluated, discuss the predicted power cycle performance, and compare with stand-alone natural gas and geothermal power generation cases. The power cycle investigated is expected to be applicable for use with conventional hydrothermal resources as well as with geologic thermal energy storage applications and/or enhanced geothermal systems. The steam topping cycle could use a concentrated solar heat source for fully renewable hybrid plant configuration. A plant initially constructed with a natural gas topping cycle heat source could be converted to a solar heat source part way though the power cycle operational life to achieve life cycle carbon emission reductions.

concentrating solar power↗
DOE OSTI2025

Vidyut3d: A Gpu Accelerated Fluid Solver for Non-Equilibrium Plasmas on Adaptive Grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

Sitaraman, Hariswaran↗
DOE OSTI2020

Analysis of the performance of a hybrid CPU/GPU 1D2D coupled model for real flood cases

Coupled 1D2D models emerged as an efficient solution for a two-dimensional (2D) representation of the floodplain combined with a fast one-dimensional (1D) schematization of the main channel. At the same time, high-performance computing (HPC) has appeared as an efficient tool for model acceleration. In this work, a previously validated 1D2D Central Processing Unit (CPU) model is combined with an HPC technique for fast and accurate flood simulation. Due to the speed of 1D schemes, a hybrid CPU/GPU model that runs the 1D main channel on CPU and accelerates the 2D floodplain with a Graphics Processing Unit (GPU) is presented. Since the data transfer between sub-domains and devices (CPU/GPU) may be the main potential drawback of this architecture, the test cases are selected to carry out a careful time analysis. Here, the results reveal the speed-up dependency on the 2D mesh, the event to be solved and the 1D discretization of the main channel. Additionally, special attention must be paid to the time step size computation shared between sub-models. In spite of the use of a hybrid CPU/GPU implementation, high speed-ups are accomplished in some cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Sierra/SolidMechanics 4.56 Capabilities In Development

This user's guide documents capabilities in Sierra/SolidMechanics which remain "in-development" and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.56 User's Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Update of ASHRAE Standard 140 Section 5.2 and Related Sections (BESTEST Building Thermal Fabric Test Cases)

This report documents new work intended to update the current building thermal fabric software-to-software comparative test specifications and example results of ANSI/ASHRAE Standard 140-2017, Sections 5.2.1, 5.2.2, 5.2.3, and related supporting sections. The Building Energy Simulation Test and Diagnostic Method (BESTEST) test specification and example results that provide the basis for these Standard 140 sections were originally published in 1995. The state of the art in building energy simulation modeling has advanced since then. Additionally, comments on the test specification over the years highlighted ambiguities, further necessitating the update to the original work. This new work was conducted for the U.S. Department of Energy (DOE) by Argonne National Laboratory (Argonne) and National Renewable Energy Laboratory (NREL), in collaboration with a working group of international experts associated with the American Society of Heating, Refrigerating, and Air-Conditioning Engineers (ASHRAE) Standing Standard Project Committee (SSPC) 140. SSPC 140 is responsible for ANSI/ASHRAE Standard 140. The new work, presented in this report, defined the test cases and specifications such that building energy simulation programs with time steps of one hour or less can perform the tests. Revisions to the test specification and update of example results follow the vetting process of iterative field trials established for previous BESTEST work that has been adapted for ASHRAE Standard 140.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Multigrid-Reduction-in-Time for the Rotating Shallow Water Equations

We consider multilevel time-parallel methods for the numerical solution of the rotating shallow water equations. In particular, the multigrid-reduction-in-time (MGRIT) algorithm is used for the parallel time integration. An asymptotic model is used at the coarse levels while the full model is employed at the finer levels. The asymptotic model is well-suited for highly oscillatory partial differential equations like the rotating shallow water equations because it can accurately and stably take the required large time-steps on coarse levels. Our work exploits the flexibility of the MGRIT algorithm in terms of the number of levels and relaxation schemes to show some computational benefits, especially with respect to FCF-relaxation and data reuse.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A File Format and API for Dynamic Radar Cross Section Data

Often the Radar Cross-Section (RCS) of a target is incorrectly assumed to be a single number by those unfamiliar with electromagnetic scattering. In actuality, a target's RCS depends on many factors. These factors include radar signal frequency, radar observation angle, as well as target orientation. Another possible parameter (often not considered) is time. The RCS of targets may change over time due to movement, environmental changes, etc. In order to accurately represent the dynamic RCS of a target in a time-stepped analysis, the ability to interface with large RCS datasets efficiently is desired. To this end, a file format and API (written in C++) were developed and are described in this report.

42 ENGINEERING↗
DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗
DOE OSTI2020

Cavitating Acoustics Modeling in Sierra/SD

Motivated by the need to simulate the effects of underwater explosion on ship structures, we develop a new cavitating acoustics formulation. The proposed approach is consistent with existing methods where the cavitation phenomenon is captured with a bilinear constitutive law. However, the new formulation is in terms of velocity potential, as opposed to the existing displacement-potential and pressure formulations. Also unique to the proposed formulation is a new generalized time-stepping procedure specific to cavitating acoustics, which has the ability to introduce numerical damping to control frothing. Numerical examples of varying complexity are presented to illustrate the effectiveness of the proposed approach and the ability to use velocity potential as a primary field variable for cavitating acoustics simulations.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2020

Modeling the Functional Forms of Grid Disturbances

This report introduces a functional form that may be used to quantitatively predict the impacts of new grid tools and changing system qualities on the likelihoods, durations, and depths of various grid disturbances. Each disturbance scenario is modeled to have three component stages—avoidance, reactance and recovery, which together parametrically estimate one disturbance’s impacts. The modeled scenario is then placed and replicated within an analysis period to represent the likelihood or frequency of the scenario and its consequent impacts. Whereas analysts have struggled to define and apply metrics for grid resilience, the functional form introduced by this report shares units of measurement with accepted grid-status measures (e.g., numbers of customers currently experiencing a service outage). Furthermore, the integrated and averaged functional form over an analysis period provides a meaningful normalized performance metric (e.g., customer outage minutes per year) that is ultimately independent of the duration of the period of. The approach may be applied similarly regardless of the severity or frequency of the disturbances that are being analyzed. Because metrics can be chosen to be identical in both the hypothetical future and the actual historical past, the historical past eventually becomes the test of the future predictions, at least in a statistical sense. The authors originally developed this approach to facilitate analysis of the effects of transactive energy (TE) systems effects on electric power grid resilience. TE systems invite energy suppliers and consumers to actively collaborate toward the discovery of, and their responses to, the locational value of energy. The findings from this process are often embodied as energy prices, the dynamics of which indicate the locational value of energy and can further represent important grid service needs. While some academic papers claim to quantify the value of a specific TE system design toward grid resilience, the answer, in general, has been elusive. Not only do multiple and conflicting definitions of resilience and reliability exist, but countless TE systems are being invented. We conclude the following: (1) The ideal analysis should harmonize rather than differentiate resilience and reliability. Therefore, this report uses the more general term disturbance whenever the overloaded terms resilience and reliability can be avoided. (2) The effectiveness of TE systems must be mapped to underlying qualities of a TE system, thereby avoiding presumptions that every TE design offers similar advantages. The authors seek to evaluate the parametric effects of TE system qualities (e.g., spatial granularity, granularity of time steps, length of future prediction horizon) on avoiding, reacting to, and recovering from grid disturbances. Furthermore, any advantages (or disadvantages) must be fairly compared with the many alternative tools, systems, and strategies that might offer comparable benefits.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Sierra/SolidMechanics 5.0 Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.0 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and $\textit{J}$-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Stability and Conservation properties of Hermite-based approximations of the Vlasov-Poisson System

Spectral approximation based on Hermite-Fourier expansion of the Vlasov-Poisson model for a collisionless plasma in the electro-static limit is provided by adding high-order artificial collision operators of Lenard-Bernstein type. These differential operators are suitably designed in order to preserve the physically-meaningful invariants (number of particles, momentum, energy). In view of time-discretization, stability results in appropriate norms are presented. In this study, necessary conditions link the magnitude of the artificial collision term, the number of spectral modes of the discretization, as well as the time-step. The analysis, carried out in full for the Hermite discretization of a simple linear problem in one-dimension, is then partly extended to cover the complete nonlinear Vlasov-Poisson model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Sierra/SolidMechanics 4.58. Capabilities In Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 4.58 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as peridynamics and the reproducing kernel particle method (RKPM), numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and /-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Sierra/SolidMechanics 5.2 Capabilities in Development

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.2 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

42 ENGINEERING↗
DOE OSTI2020

Local-nonlocal coupling in Emu/PDMS

A technique called the splice method for coupling local to peridynamic subregions of a body is described. The method relies on ghost nodes, whose values of displacement are interpolated from nearby physical nodes, to make each subregion visible to the other. In each time step, the nodes in each subregion treat the nodes in the other subregion as boundary conditions. Adaptively changing the subregions is possible through the creation and deletion of ghost nodes. Example problems in 2D and 3D illustrate how the method is used to perform multiscale modeling of fracture and impact events within a larger structure.

42 ENGINEERING↗