SEARCH · Engineering Papers

Results for “Model initialization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

DOE OSTI2016

CONVECTIVE PROPERTIES OF ROTATING TWO-DIMENSIONAL CORE-COLLAPSE SUPERNOVA PROGENITORS

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20 M {sub ⊙} star. Using the Modules for Experiments in Stellar Astrophysics we construct one-dimensional (1D) stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in 1D stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of the resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

Effects of electron-phonon coupling and electronic thermal conductivity in high energy molecular dynamics simulations of irradiation cascades in nickel

The two-temperature model has been applied to investigate the effects of the electronic subsystem on 150 keV Ni ion cascades in nickel using molecular dynamics simulation. We explore the effects of the magnitude of the electron-phonon coupling and the electronic thermal conductivity on defect production and cluster formation. It has been found that stronger electron-phonon coupling allows larger and more rapid energy feedback to the atomic subsystem, leading to reduction of number of point defects and suppression of the formation of larger defect clusters. It was observed that larger electronic thermal conductivity results in slightly increased number of point defects and larger size vacancy clusters. The latter takes place because of suppression of point defects recombination in faster cooling areas of initial damage.

Zarkadoula, Eva↗

DOE OSTI2017

Local and global Λ polarization in a vortical fluid

We compute the fermion spin distribution in the vortical fluid created in off-central high energy heavy-ion collisions. We employ the event-by-event (3+1)D viscous hydrodynamic model. The spin polarization density is proportional to the local fluid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial fluid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Robustness of the Ensemble Score Filter to the Type of Assimilated Observation Networks

Recent advances in data assimilation (DA) have focused on developing more flexible approaches that can better accommodate nonlinearities in models and observations. However, it remains unclear how the performance of these advanced methods depends on the observation network characteristics. In this study, we present initial experiments with the surface quasi‐geostrophic model, in which we compare a recently developed ensemble filter using score‐based diffusion models with the standard Local Ensemble Transform Kalman Filter (LETKF). Our results show that the analysis solutions respond differently to the number, spatial distribution, and nonlinear fraction of assimilated observations. We also find notable changes in the multiscale characteristics of the analysis errors. Given that standard DA techniques will eventually be replaced by more advanced methods, we hope this study sets the ground for future efforts to reassess the value of Earth observing systems in the context of newly emerging algorithms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2020

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗

DOE OSTI2019

The DOE E3SM Coupled Model Version 1: Description and Results at High Resolution

This study provides an overview of the coupled high-resolution Version 1 of the Energy Exascale Earth System Model (E3SMv1) and documents the characteristics of a 50-year-long high-resolution control simulation with time-invariant 1950 forcings following the HighResMIP protocol. In terms of global root-mean-squared error metrics, this high-resolution simulation is generally superior to results from the low-resolution configuration of E3SMv1 (due to resolution, tuning changes, and possibly initialization procedure) and compares favorably to models in the CMIP5 ensemble. Ocean and sea ice simulation is particularly improved, due to better resolution of bathymetry, the ability to capture more variability and extremes in winds and currents, and the ability to resolve mesoscale ocean eddies. The largest improvement in this regard is an ice-free Labrador Sea, which is a major problem at low resolution. Interestingly, several features found to improve with resolution in previous studies are insensitive to resolution or even degrade in E3SMv1. Most notable in this regard are warm bias and associated stratocumulus deficiency in eastern subtropical oceans and lack of improvement in El Niño. Another major finding of this study is that resolution increase had negligible impact on climate sensitivity (measured by net feedback determined through uniform +4K prescribed sea surface temperature increase) and aerosol sensitivity. Cloud response to resolution increase consisted of very minor decrease at all levels. Large-scale patterns of precipitation bias were also relatively unaffected by grid spacing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Basin and Petroleum System Modeling of Gas Hydrate Deposits in the Walker Ridge Area, Northern Gulf of Mexico - Poster

Poster presented by NETL at ICGH, 06/17 - Abstract: We are currently constructing an Earth model of the Terrebonne Basin for comprehensive 3D modeling. Well data and a high-quality 3D seismic survey underpin the geologic conditions for the model. Additional model constraints include basal heat flow, seafloor temperature, and paleowater depth through time. Our initial investigation comprises the development of a “pseudo” 2D model centered on the deep WR313 #001 well. We are using this test case to establish the workflow for integrating BPSM with gas hydrate analyses; this will eventually be applied to the entire mini-basin. A key aspect of investigating the gas hydrate system with BPSM is integrating it with the other active petroleum system(s) in the region, because the methane in hydrate deposits can derive from both biogenic and thermogenic sources. Some of the questions we may consider are: What is the history of these deposits with respect to the gas hydrate stability zone, recognizing that this area is not one characterized by simple continuous sedimentation and subsidence, but instead by periodic uplift, tilting, plate rotation, sea-level change, salt tectonics, and evolution of temperature gradients. In conclusion, BPSM has been called the ‘great integrator’ in petroleum exploration. Development of a BPSM model of gas hydrates in the Terrebonne Basin of the northern GoM will provide a vehicle within which to integrate other early exploration and assessment research being conducted on gas hydrates, such as in frontier areas like Japan, as a potential resource likely to provide many decades of energy if proven to be commercially producible

Dafov, Laura↗

DOE OSTI2025

Insights Into Thermal Runaway Mechanisms: Fast Tomography Analysis of Metal Agglomerates in Lithium-Ion Batteries

Thermal Runaway (TR) in lithium-ion batteries (LIB) is a critical technological and social concern. Whilst such events are rare, TR is characterized by uncontrollable heating leading to catastrophic failures. To deepen the understanding of the failure process and subsequently develop more accurate TR prediction models and as a result safer battery systems, we present in this work high-speed X-ray tomography for in-depth investigations of the copper current collector melting and agglomeration during TR. The melting process presents valuable real-time internal information about heat evolution during TR, previously challenging to access but crucially important for validating TR models. In this work, controlled failure studies combined with high-speed X-ray tomography were performed on two different commercial LIB models, subjecting them to both external heating and nail penetration to induce TR. Through real-time observation via high-speed tomography, followed by segmentation, rendering, and analysis, the formation of copper agglomerates was qualitatively and quantitatively characterized and visualized for the first time. Agglomerates tended to form either from the battery's outermost layers or centrally, depending on the method of TR initiation, and gives an indirect insight into the internal temperature evolution and distribution. Moreover, an initial comparative analysis between the battery models also revealed differences in agglomerate size, which has been linked to the thicker copper current collectors of one of the cell models. We further discuss the impact of larger copper agglomerates on heat distribution and safety. This study not only sheds light on the intricate dynamics of TR in LIBs but also underscores the pivotal role of 'gold-standard' imaging techniques in advancing battery safety, crucial for the robust modeling of TR and the future design of electric vehicle safety systems.

25 ENERGY STORAGE↗

DOE OSTI2024

Utilizing waste heat in wastewater treatment plants for water desalination: Modeling and Multi-Objective optimization of a Multi-Effect desalination system using Decision Tree Regression and Pelican optimization algorithm

This paper examines the feasibility of using waste heat from wastewater treatment plants (WWTPs) for water desalination. A model was developed to utilize waste heat from the gensets at As Samra WWTP in Jordan, using real data and TRNSYS® software to calculate available waste heat. The desalination process was then modeled with ASPEN PLUS® software, focusing on multi-effect desalination (MED). Both series and parallel configurations for the MED system were compared. The study investigated the effects of system feeding flow rate, feeding pressure, and heat input on productivity, performance ratio, and recovery ratio. The study also introduces a novel optimization technique combining machine learning and modern optimization algorithms to maximize system productivity and performance. Initially, a decision tree regression (DTR) model is developed to establish relationships between key independent variables (flow rate, feed pressure, and heat input) and dependent variables (productivity, performance ratio, and recovery ratio). The Pelican Optimization Algorithm (POA) is then used to identify the optimal values of the independent variables for maximum productivity and performance. The results show that using a series configuration yields a system productivity of 3984.2 kg/hr, a performance ratio of 3.78, and a recovery ratio of 0.991 at a feed flow rate of 4000 kg/hr, feed pressure of 3 bars, and heat input of 719 kW. Optimal productivity (4421 kg/hr), performance ratio (3.81), and recovery ratio (0.851) are achieved at a feed flow rate of 5166 kg/hr, feed pressure of 3.2 bars, and heat input of 794 kW. In conclusion, the techno-economic assessment indicates a levelized cost of water of 1.63 USD/m 3 for parallel configurations and 1.65 USD/m 3 for series configurations, with a payback period of less than two years.

42 ENGINEERING↗

DOE OSTI2024

Emulator of PR‐DNS: Accelerating Dynamical Fields With Neural Operators in Particle‐Resolved Direct Numerical Simulation

Abstract Particle‐resolved direct numerical simulations (PR‐DNS) play an increasing role in investigating aerosol‐cloud‐turbulence interactions at the most fundamental level of processes. However, the high computational cost associated with high resolution simulations poses considerable challenges for large domain or long duration simulation using PR‐DNS. To address these issues, here we present an emulator of the complex physics‐based PR‐DNS developed by use of the data‐driven Fourier Neural Operator (FNO) method. The effectiveness of the method is showcased by presenting turbulence and temperature fields in a two‐dimensional space. The results demonstrate high accuracy at various resolutions and the emulator is two orders of magnitude cheaper in terms of computational demand compared to the physics‐based PR‐DNS model. Furthermore, the FNO emulator exhibits strong generalization capabilities for different initial conditions and ultra‐high‐resolution without the need to retrain models. These findings highlight the potential of the FNO method as a promising tool to simulate complex fluid dynamics problems with high accuracy, computational efficiency, and generalization capabilities, enhancing our understanding of the aerosol‐cloud‐precipitation system.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2016

RECOVERY FROM GIANT ERUPTIONS IN VERY MASSIVE STARS

We use a hydro-and-radiative-transfer code to explore the behavior of a very massive star (VMS) after a giant eruption—i.e., following a supernova impostor event. Beginning with reasonable models for evolved VMSs with masses of 80 M{sub ⊙} and 120 M{sub ⊙}, we simulate the change of state caused by a giant eruption via two methods that explicitly conserve total energy. (1) Synthetically removing outer layers of mass of a few M{sub ⊙} while reducing the energy of the inner layers. (2) Synthetically transferring energy from the core to the outer layers, an operation that automatically causes mass ejection. Our focus is on the aftermath, not the poorly understood eruption itself. Then, using a radiation-hydrodynamic code in 1D with realistic opacities and convection, the interior disequilibrium state is followed for about 200 years. Typically the star develops a ∼400 km s{sup −1} wind with a mass loss rate that begins around 0.1 M{sub ⊙} yr{sup −1} and gradually decreases. This outflow is driven by κ-mechanism radial pulsations. The 1D models have regular pulsations but 3D models will probably be more chaotic. In some cases a plateau in the mass-loss rate may persist about 200 years, while other cases are more like η Car which lost >10 M{sub ⊙} and then had an abnormal mass loss rate for more than a century after its eruption. In our model, the post-eruption outflow carried more mass than the initial eruption. These simulations constitute a useful preliminary reconnaissance for 3D models which will be far more difficult.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

SAM Theory Manual

The System Analysis Module (SAM) is an advanced and modern system analysis tool being developed at Argonne National Laboratory under the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. SAM development aims for advances in physical modeling, numerical methods, and software engineering to enhance its user experience and usability for reactor transient analyses. To facilitate the code development, SAM utilizes an object-oriented application framework (MOOSE), and its underlying meshing and finite-element library (libMesh) and linear and non-linear solvers (PETSc), to leverage modern advanced software environments and numerical methods. SAM focuses on modeling advanced reactor concepts such as SFRs (sodium fast reactors), LFRs (lead-cooled fast reactors), and FHRs (fluoride-salt-cooled high temperature reactors) or MSRs (molten salt reactors). These advanced concepts are distinguished from light-water reactors in their use of single-phase, low-pressure, high-temperature, and low Prandtl number (sodium and lead) coolants. As a new code development, the initial effort has been focused on modeling and simulation capabilities of heat transfer and single-phase fluid dynamics responses in Sodium-cooled Fast Reactor (SFR) systems. The system-level simulation capabilities of fluid flow and heat transfer in general engineering systems and typical SFRs have been verified and validated. This document provides the theoretical and technical basis of the code to help users understand the underlying physical models (such as governing equations, closure models, and component models), system modeling approaches, numerical discretization and solution methods, and the overall capabilities in SAM. As the code is still under ongoing development, this SAM Theory Manual will be updated periodically to keep it consistent with the state of the development.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Peridynamic Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. Here, the computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

42 ENGINEERING↗

DOE OSTI2015

THE DYNAMICAL EVOLUTION OF STELLAR BLACK HOLES IN GLOBULAR CLUSTERS

Our current understanding of the stellar initial mass function and massive star evolution suggests that young globular clusters (GCs) may have formed hundreds to thousands of stellar-mass black holes (BHs), the remnants of stars with initial masses from ∼20-100 M {sub ☉}. Birth kicks from supernova explosions may eject some BHs from their birth clusters, but most should be retained. Using a Monte Carlo method we investigate the long-term dynamical evolution of GCs containing large numbers of stellar BHs. We describe numerical results for 42 models, covering a broad range of realistic initial conditions, including up to 1.6 × 10{sup 6} stars. In almost all models we find that significant numbers of BHs (up to ∼10{sup 3}) are retained all the way to the present. This is in contrast to previous theoretical expectations that most BHs should be ejected dynamically within a few gigayears The main reason for this difference is that core collapse driven by BHs (through the Spitzer {sup m}ass segregation instability{sup )} is easily reverted through three-body processes, and involves only a small number of the most massive BHs, while lower-mass BHs remain well-mixed with ordinary stars far from the central cusp. Thus the rapid segregation of stellar BHs does not lead to a long-term physical separation of most BHs into a dynamically decoupled inner core, as often assumed previously. Combined with the recent detections of several BH X-ray binary candidates in Galactic GCs, our results suggest that stellar BHs could still be present in large numbers in many GCs today, and that they may play a significant role in shaping the long-term dynamical evolution and the present-day dynamical structure of many clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

3D glasma initial state for relativistic heavy ion collisions

We extend the impact-parameter-dependent Glasma model to three dimensions using explicit small-x evolution of the two incoming nuclear gluon distributions. We compute rapidity distributions of produced gluons and the early-time energy momentum tensor as a function of space-time rapidity and transverse coordinates. Finally, we study rapidity correlations and fluctuations of the initial geometry and multiplicity distributions and make comparisons to existing models for the three-dimensional initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2016

Mean-Field-Game Model for Botnet Defense in Cyber-Security

We initiate the analysis of the response of computer owners to various offers of defence systems against a cyber-hacker (for instance, a botnet attack), as a stochastic game of a large number of interacting agents. We introduce a simple mean-field game that models their behavior. It takes into account both the random process of the propagation of the infection (controlled by the botner herder) and the decision making process of customers. Its stationary version turns out to be exactly solvable (but not at all trivial) under an additional natural assumption that the execution time of the decisions of the customers (say, switch on or out the defence system) is much faster that the infection rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗