Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “building 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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Development of a metamodelling framework for building energy models with application to fifth-generation district heating and cooling networks

Fully defined physics-based building energy models can accurately represent building systems; however, generating models based on high-level parameters is time consuming and simulation time of complex models can be slow. This article discusses the development of a Metamodelling Framework to create metamodels from a building energy modelling dataset. The framework generates metamodels using either linear regression, random forests, or support vector regressions. A fifth-generation district heating and cooling system analysis use case was used to motivate the development of the framework. The use case required quick and accurate representations of annual building loads reported hourly. Typical annual building modelling approaches can result in a runtime of 10 min. The metamodels runtime was reduced to less than 10 s to load and run an annual simulation with user-defined covariates. The results of the metamodel performance and an abbreviated topology analysis based on the motivating use case will be presented.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Gaussian process-based surrogate modeling framework for process planning in laser powder-bed fusion additive manufacturing of 316L stainless steel

Here, Laser Powder-Bed Fusion (L-PBF) metal-based additive manufacturing (AM) is complex and not fully understood. Successful processing for one material, might not necessarily apply to a different material. This paper describes a workflow process that aims at creating a material data sheet standard that describes regimes where the process can be expected to be robust. The procedure consists of building a Gaussian process-based surrogate model of the L-PBF process that predicts melt pool depth in single-track experiments given a laser power, scan speed, and laser beam size combination. The predictions are then mapped onto a power versus scan speed diagram delimiting the conduction from the keyhole melting controlled regimes. This statistical framework is shown to be robust even for cases where experimental training data might be suboptimal in quality, if appropriate physics-based filters are applied. Additionally, it is demonstrated that a high-fidelity simulation model of L-PBF can equally be successfully used for building a surrogate model, which is beneficial since simulations are getting more efficient and are more practical to study the response of different materials, than to re-tool an AM machine for new material powder.

316L stainless steel↗
DOE OSTI2022

Modeling rarefied gas chemistry with QuiPS, a novel quasi-particle method

The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Overall, variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.

42 ENGINEERING↗
DOE OSTI2023

Residual estimation for grid modification in wall-modeled large eddy simulation using unstructured high-order methods

Here, the accuracy and computational cost of a large eddy simulation are highly dependent on the computational grid. Building optimal grids manually from a priori knowledge is not feasible in most practical use cases; instead, solution-adaptive strategies can provide a robust and cost-efficient method to generate a grid with the desired accuracy. We adapt the residual estimation algorithm developed by Toosi and Larsson for Discontinuous Galerkin Spectral Elements Methods (DGSEM) to guide the grid-adaptation process. The core of the method is the computation of the estimated modeling residual using the polynomial basis functions used in DGSEM, and the averaging of the estimated residual over each element. The final method is assessed in multiple channel flow test cases and for the transonic flow over an airfoil, in both cases making use of mortar interfaces between elements with hanging nodes. The method is found to be robust and reliable, and to provide solutions on grids with significantly fewer elements at comparable accuracy compared to when using human-generated grids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Neural-based time series forecasting of loss of coolant accidents in nuclear power plants

During the last few years, deep learning in neural networks has demonstrated impressive successes in the areas of computer vision, speech and image recognition, text generation, and many others. However, sensitive engineering areas such as nuclear engineering benefited less from these efficient techniques. In this work, deep learning expert systems are utilized to model and predict time series progression of a design-basis nuclear accident, featuring a loss of coolant accident. Two major findings are accomplished in this work. First, the ability to train expert systems with high accuracy, which could help nuclear power plant operators to figure out plant responses during the accident. Second, building fast, efficient, and accurate deep models to simulate nuclear phenomena, which could be valuable to nuclear computational science. In this work, large amount of time series data is obtained from simulation tools by simulating different conditions of the base-case/nominal accident scenario. Four critical outputs/responses are monitored during the accident (e.g. temperature, pressure, break flow rate, water level). Two approaches are adopted in this work. The first approach is to use feedforward deep neural networks (DNN) to fit all time steps and outputs in a single model. The second approach is to use long short-term memory (LSTM) to fit all time steps together for each reactor response separately. Both DNN and LSTM demonstrate very good performance in predicting the test and base-case scenarios, with accuracy as low as 92% and as high as 99%, where these test scenarios are unknown to the expert systems and are not included in the model training. In addition, both approaches demonstrate a significant reduction in computational costs, as the deep expert system is able to accurately predict the accident 100,000 times faster than the original simulation tool. Given sufficient data, the methodology adopted in this study demonstrates that DNN/LSTM expert systems can be used as a decision support system to model advanced time series phenomena within nuclear power plants with high accuracy and negligible computational costs.

42 ENGINEERING↗
DOE OSTI2025

MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8) for advanced open-source tritium transport and fuel cycle modeling

Tritium management is critical for the safety, sustainability, and economics of fusion energy systems, and advanced and reliable modeling tools help accelerate the development of tritium technologies. This paper presents the Tritium Migration Analysis Program, Version 8 (TMAP8), an open-source, MOOSE-based application developed to provide state-of-the-art tritium transport and fuel cycle modeling capabilities. TMAP8 aims to expand the capabilities of previous versions (i.e., TMAP4 and TMAP7) by leveraging modern computational techniques, ensuring high software quality assurance standards (key to building trust), and enabling multispecies, multiscale, and multiphysics simulations for integrated tritium transport modeling in complex geometries. This paper outlines TMAP8’s scope and rigorous development practices, emphasizing its transparency, accessibility, modularity, and reliability. We present the current suite of verification and validation cases based on those from TMAP4, demonstrating TMAP8’s accuracy and reliability against analytical solutions and experimental data. Additionally, the paper showcases TMAP8’s integrated fuel cycle modeling capabilities, highlighting its applicability at various scales and levels. The TMAP8 code and documentation are openly available, promoting collaborative development and widespread adoption within the fusion community. Future work will soon expand TMAP8’s verification and validation suite to include those from TMAP7 and other recent experimental studies for validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Toward Lithium Recovery Using Modular (and Membraneless) Phase Separation and Extraction (MPSE) Technology with Ionic Liquid (IL) Solvents: Effect of Coatings

In this work, a single-channel slope-plate Modular and Membrane-less Phase Separation and Extraction (MPSE) device was fabricated using 3D printing. The slope-plate was designed with an insertable glass slide to enable surface modification with coatings. Self-assembled monolayers (SAMs) and perfluoropolyether (Zdol) coatings were applied to tailor surface wettability and enhance phase separation. Using a model biphasic system of water and hexadecane, both coatings significantly improved the separation efficiency compared to the uncoated device by promoting selective wetting. Building on these results, lithium extraction from a simulated saline solution was investigated using 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][NTf₂]) as the ionic liquid (IL) extractant. The IL effectively extracted Li⁺ from the aqueous phase in bulk extraction, demonstrating its potential for lithium recovery. However, the low interfacial tension between the IL and aqueous phases posed challenges for phase separation. The application of SAM and Zdol coatings effectively mitigated this issue. Overall, integrating tailored surface chemistry with the slope-plate MPSE design shows great promise as an efficient and scalable platform for studying liquid–liquid separation and optimizing ionic liquid–based extraction processes for lithium recovery from saline sources.

3D printing↗
DOE OSTI2025

A Universal Model of Cation Effects in Electrocatalysis

Electrolyte cations are conventionally viewed as inert spectators in electrocatalysis. However, a wealth of observations show that catalytic rates are often highly sensitive to cation identity. Despite their prevalence, these cation effects have resisted a unified mechanistic explanation, with different physical phenomena implicated across reaction chemistries, catalyst compositions, and choice of solvent. In this perspective, we describe a general framework for understanding cation effects in electrocatalysis based on electrostatics. We argue that cations influence reaction rates by modifying the strength of the electric field present at the catalyst surface, which alters the energetics of adsorbed intermediates and transition states according to their dipole moments and polarizabilities. The magnitude of this field depends on how cations arrange at the electrode surface, controlled by their size, shape, solvation, and packing efficiency. Cations that can arrange more densely result in a steeper potential drop at the electrode surface and consequently a stronger electric field. Our model further identifies two criteria for observing cation effects: (1) the operating potential must be negative of the electrode’s potential of zero total charge, ensuring that cations accumulate at the interface, and (2) the energetics of the kinetically relevant elementary step must be field sensitive. This framework reconciles previously inconsistent trends, including why cation effects appear only for some catalysts, why reaction selectivity is sensitive to cation identity, and why activity can increase with cation size on certain metals but decrease on others. Supported by kinetic measurements, spectroscopy, and atomistic simulations, the model provides both conceptual value for building intuition about catalysis at charged interfaces and predictive value for anticipating trends for new reactions, catalysts, and electrolytes. We conclude by highlighting the importance of electric fields across electrochemical, thermochemical, and biological catalysis and propose that considering the electrostatic environment around active sites offers new opportunities for improving activity and selectivity.

cation effects↗
DOE OSTI2021

The Average Shape of Sea Ice Ridge Keels

Through analysis of over 64,000 ridge profiles identified from moored upward-looking sonars, we identify a well-defined average shape of pressure ridge keels that is concave or cusped, not triangular as widely assumed in other literature. On the basis of this average shape, we put forward a new, dimensional-definition, of a pressure ridge cross-section that follows a negative exponential form and allows an average pressure ridge cross section to be constructed with knowledge or choice of a single parameter. The horizonal asymptote of the profile represents the draft of the “background” ice in which the ridge is embedded. The draft of this background ice scales with keel depth and is typically greater than can be accounted through thermodynamic growth, indicating that ridges tend to be embedded in fields of mechanically thickened rubble. Using a variational ridge model we simulated keel shapes for a range of ridge building conditions. Here, the model results agree well with the observations and indicate the cusped shape of an average ridge profile arises from the varying angle of horizontal shear in the ice cover when ridges form. The modeling results also explain the elongated tail of the area draft distribution of the pack.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Non-adiabatic excited state molecular dynamics of phenylene ethynylene dendrimer using a multiconfigurational Ehrenfest approach

Photoinduced dynamics of electronic and vibrational unidirectional energy transfer between meta-linked building blocks in a phenylene ethynylene dendrimer is simulated using a multiconfigurational Ehrenfest in time-dependent diabatic basis (MCE-TDDB) method, a new variant of the MCE approach developed by us for dynamics involving multiple electronic states with numerous abrupt crossings. Excited-state energies, gradients and non-adiabatic coupling terms needed for dynamics simulation are calculated on-the-fly using the Collective Electron Oscillator (CEO) approach. In conclusion, a comparative analysis of our results obtained using MCE-TDDB, the conventional Ehrenfest method and the surface-hopping approach with and without decoherence corrections is presented.

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

Modeling the nonlinear hysteretic response in DAE experiments of Berea sandstone: A case-study

Dynamic acousto-elasticity (DAE) allows probing the instantaneous state of a material while the latter slowly and periodically is changed by an external, dynamic source. In DAE investigations of geo-materials, hysteresis of the material's modulus defect displays intriguing features which have not yet been interpreted in terms of any specific mechanism occurring at atomic or mesoscale. Here, experimental results on dry Berea sandstone, which is the rock type best investigated by means of a DAE technique, are analyzed in terms of three rheological models providing simplified representations of mechanisms involving dislocations interacting with point defects which are distributed along the dislocations' core or glide planes, and microcracks with finite stiffness in compression. Constitutive relations linking macroscopic strain and stress are derived. From the latter, the modulus defect associated to each mechanism is recovered. These models are employed to construct a composite one which is capable of reproducing several of the main features observed in the experimental data. The limitations of the present approach and, possibly, of the current implementation of DAE are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Study of shock waves and related phenomena motivated by astrophysics

This study discusses the recent research in High-Energy-Density Physics at our Center. Our work in complex hydrodynamics is now focused on mode coupling in the Richtmyer-Meshkov process and on the supersonic Kelvin-Helmholtz instability. These processes are believed to occur in a wide range of astrophysical circumstances. In radiation hydrodynamics, we are studying radiative reverse shocks relevant to cataclysmic variable stars. Our work on magnetized flows seeks to produce magnetized jets and study their interactions. We build the targets for all these experiments, and simulate them using our CRASH code. We also conduct diagnostic research, focused primarily on imaging x-ray spectroscopy and its applications to scattering and fluorescence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Beyond mass: detecting secondary halo properties with galaxy-galaxy lensing

ABSTRACT Secondary halo properties beyond mass, such as the mass accretion rate (MAR), concentration, and the half mass scale, are essential in understanding the formation of large-scale structure and dark matter haloes. In this paper, we study the impact of secondary halo properties on the galaxy-galaxy lensing observable, ΔΣ. We build an emulator trained on N-body simulations to model ΔΣ and quantify the impact of different secondary parameters on the ΔΣ profile. We focus on the impact of MAR on ΔΣ. We show that a 3σ detection of variations in MAR at fixed halo mass could be achieved with the Hyper Suprime Cam survey assuming no baryonic effects and a proxy for MAR with scatter <1.5. We show that the full radial profile of ΔΣ depends on secondary properties at fixed halo mass. Consequently, an emulator that can perform full shape fitting yields better than two times improvement upon the constraints on MAR than only using the outer part of the halo. Finally, we highlight that miscentring and MAR impact the radial profile of ΔΣ in a similar fashion, implying that miscentring and MAR need to be modelled jointly for unbiased estimates of both effects. We show that present-day lensing data sets have the statistical capability to place constraints on halo MAR within our assumptions. Our analysis opens up new possibilities for observationally measuring the assembly history of the dark matter haloes that host galaxies and clusters.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Predictably missing satellites: subhalo abundances in Milky Way-like haloes

ABSTRACT On small scales there have been a number of claims of discrepancies between the standard cold dark matter (CDM) model and observations. The ‘missing satellites problem’ infamously describes the overabundance of subhaloes from CDM simulations compared to the number of satellites observed in the Milky Way. A variety of solutions to this discrepancy have been proposed; however, the impact of the specific properties of the Milky Way halo relative to the typical halo of its mass has yet to be explored. Motivated by recent studies that identified ways in which the Milky Way is atypical, we investigate how the properties of dark matter haloes with mass comparable to our Galaxy’s – including concentration, spin, shape, and scale factor of the last major merger – correlate with the subhalo abundance. Using zoom-in simulations of Milky Way-like haloes, we build two models of subhalo abundance as functions of host halo properties. From these models we conclude that the Milky Way most likely has fewer subhaloes than the average halo of the same mass. We expect up to 30 per cent fewer subhaloes with low maximum rotation velocities ($V_{\rm max}^{\rm sat} \sim 10$ km s−1) at the 68 per cent confidence level and up to 52 per cent fewer than average subhaloes with high rotation velocities ($V_{\rm max}^{\rm sat} \gtrsim 30$ km s−1, comparable to the Magellanic Clouds) than would be expected for a typical halo of the Milky Way’s mass. Concentration is the most informative single parameter for predicting subhalo abundance. Our results imply that models tuned to explain the missing satellites problem assuming typical subhalo abundances for our Galaxy may be overcorrecting.

Fielder, Catherine E.↗
DOE OSTI2020

Dirty bosons on the Cayley tree: Bose-Einstein condensation versus ergodicity breaking

Building on large-scale quantum Monte Carlo simulations, we investigate the zero-temperature phase diagram of hard-core bosons in a random potential on site-centered Cayley trees with branching number K=2. In order to follow how the Bose-Einstein condensate (BEC) is affected by the disorder, we focus on both the zero-momentum density, probing the quantum coherence, and the one-body density matrix (1BDM) whose largest eigenvalue monitors the off-diagonal long-range order. We further study its associated eigenstate which brings useful information about the real-space properties of this leading eigenmode. Upon increasing randomness, we find that the system undergoes a quantum phase transition at finite disorder strength between a long-range ordered BEC state, fully ergodic at large scale, and a new disordered Bose glass regime showing conventional localization for the coherence fraction while the 1BDM displays a nontrivial algebraic vanishing BEC density together with a nonergodic occupation in real space. These peculiar properties can be analytically captured by a simple phenomenological description on the Cayley tree which provides a physical picture of the Bose glass regime.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Monolayer CrCl 3 as an Ideal Test Bed for the Universality Classes of 2D Magnetism

The monolayer halides CrX 3 (X=Cl, Br, I) attract significant attention for realizing 2D magnets with genuine long-range order (LRO), challenging the Mermin-Wagner theorem. In this work, we show that monolayer CrCl 3 has the unique benefit of exhibiting tunable magnetic anisotropy upon applying a compressive strain. This opens the possibility to use CrCl 3 for producing and studying both ferromagnetic and antiferromagnetic 2D Ising-type LRO as well as the Berezinskii-Kosterlitz-Thouless (BKT) regime of 2D magnetism with quasi-LRO. Using state-of-the-art density functional theory, we explain how realistic compressive strain could be used to tune the monolayer's magnetic properties so that it could exhibit any of these phases. Building on large-scale quantum Monte Carlo simulations, we compute the phase diagram of strained CrCl 3 , as well as the magnon spectrum with spin-wave theory. Our results highlight the eminent suitability of monolayer CrCl 3 to achieve very high BKT transition temperatures, around 50 K, due to their singular dependence on the weak easy-plane anisotropy of the material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Spatiotemporal Crossover between Low- and High-Temperature Dynamical Regimes in the Quantum Heisenberg Magnet

The stranglehold of low temperatures on fascinating quantum phenomena in one-dimensional quantum magnets has been challenged recently by the discovery of anomalous spin transport at high temperatures. Whereas both regimes have been investigated separately, no study has attempted to reconcile them. For instance, the paradigmatic quantum Heisenberg spin-$1/2$ chain falls at low-temperature within the Tomonaga-Luttinger liquid framework, while its high-temperature dynamics is superdiffusive and relates to the Kardar-Parisi-Zhang universality class in $1+1$ dimensions. This work aims at reconciling the two regimes. Building on large-scale matrix product state simulations, we find that they are connected by a temperature-dependent spatiotemporal crossover. In this work, as the temperature $T$ is reduced, we show that the onset of superdiffusion takes place at longer length and time scales $\propto 1/T$. This prediction has direct consequences for experiments including nuclear magnetic resonance: it is consistent with earlier measurements on the nearly ideal Heisenberg $S=1/2$ chain compound Sr$_2$CuO$_3$ yet calls for new and dedicated experiments.

1-dimensional spin chains↗
DOE OSTI2016

Poster - 47: A parametrized prediction model of rectal toxicity in focal SBRT of low risk prostate cancer

There has been a recent trend towards watchful waiting in place of intervention for early stage prostate cancer (CaP). However, this approach can allow for disease progression, and subsequent whole-gland therapies such as prostatectomy and whole gland irradiation can result in functional deficits or rectal toxicities or both. A controversial alternative approach for this patient cohort is the use of focal therapy, where the treatment is focussed on an identified dominant index lesion (DIL). This work aims to investigate the treatment parameters for focal SBRT of the prostate under which clinically acceptable rectal NTCP levels can be achieved. For each of 25 low risk CaP patients, a hypothetical 2 cc DIL was modeled in the right-posterior quadrant of the prostate, and was used to build a PTV as the target for SBRT simulation. An SBRT prescriptions of 41 Gy and 37 Gy in 5 fractions were chosen, corresponding to the boost levels used in previous CaP dose escalation studies. DVH data were exported and used to calculate rectal NTCP values based on the Lyman-Kutcher-Burman (LKB) model using the QUANTEC reccommended model parameters. Rectal NTCP dependence on DIL-to-rectum separation, dose level, and DIL volume were investigated. The final goal of this ongoing work is to create a map of the maximum allowable prescription dose for a given patient geometry that achieves a clinically acceptable rectal NTCP level.

60 APPLIED LIFE SCIENCES↗