Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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

A hybrid smoothed-particle hydrodynamics model of oxide skins on molten aluminum

A computational model of aluminum melting is proposed which captures both the thermal fluid-solid phase transition and the mechanical effects of oxidation. The model hybridizes ideas from smoothed particle hydrodynamics and bonded particle models to simulate both hydrodynamic flows and solid elasticity. Oxidation is represented by dynamically adding and deleting spring-like bonds between surface fluid particles to represent the formation and rupture of the oxide skin. Various complex systems are simulated to demonstrate the adaptability of the method and to illustrate the significant impact of skin properties on material flow. As a result, initial comparison to experiments of a melting aluminum cantilever highlights that the computational model can reproduce key qualitative features of aluminum relocation.

Clemmer, Joel Thomas↗
DOE OSTI2021

Deterioration of concrete due to ASR: Experiments and multiscale modeling

The process of ASR (Alkali-Silica Reaction) induced expansion and damage in pavement concrete specimens is investigated using laboratory experiments and computational modeling. In the experimental program, the concrete specimens are subject to CS-CPT (climate simulation concrete prism test) to obtain ASR induced expansion with and without external supply of alkali. The dissolution rates of the granodiorite used in the concrete mix and the gel formation rates are determined under concrete-like conditions (pH 13.8, with/without Ca(OH){sub 2} and NaCl) at different temperatures. A micromechanics based computational model with aggregate-scale diffusion and reaction kinetics coupled to an Eigenstrain based micromechanics damage model is developed for the simulation of ASR induced expansion and damage. Data from the experimental program are used to calibrate and validate the computational model. Model predictions show that for the given concrete mixture, ASR induced expansion in the specimen exposed to water is predominantly governed by microcracking in the aggregate, while the expansion in the specimen subjected to external alkali supply is governed by microcracking in both the aggregates and the cement paste.

36 MATERIALS SCIENCE↗
DOE OSTI2020

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗
DOE OSTI2022

Quantum Markov chain Monte Carlo with digital dissipative dynamics on quantum computers

Modeling the dynamics of a quantum system connected to the environment is critical for advancing our understanding of complex quantum processes, as most quantum processes in nature are affected by an environment. Modeling a macroscopic environment on a quantum simulator may be achieved by coupling independent ancilla qubits that facilitate energy exchange in an appropriate manner with the system and mimic an environment. This approach requires a large, and possibly exponential number of ancillary degrees of freedom which is impractical. In contrast, we develop a digital quantum algorithm that simulates interaction with an environment using a small number of ancilla qubits. By combining periodic modulation of the ancilla energies, or spectral combing, with periodic reset operations, we are able to mimic interaction with a large environment and generate thermal states of interacting many-body systems. We evaluate the algorithm by simulating preparation of thermal states of the transverse Ising model. Our algorithm can also be viewed as a quantum Markov chain Monte Carlo process that allows sampling of the Gibbs distribution of a multivariate model. To illustrate this we evaluate the accuracy of sampling Gibbs distributions of simple probabilistic graphical models using the algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Validation Metrics for Fixed Effects and Mixed-Effects Calibration

The modern scientific process often involves the development of a predictive computational model. To improve its accuracy, a computational model can be calibrated to a set of experimental data. A variety of validation metrics can be used to quantify this process. Some of these metrics have direct physical interpretations and a history of use, while others, especially those for probabilistic data, are more difficult to interpret. In this work, a variety of validation metrics are used to quantify the accuracy of different calibration methods. Frequentist and Bayesian perspectives are used with both fixed effects and mixed-effects statistical models. Through a quantitative comparison of the resulting distributions, the most accurate calibration method can be selected. Two examples are included which compare the results of various validation metrics for different calibration methods. It is quantitatively shown that, in the presence of significant laboratory biases, a fixed effects calibration is significantly less accurate than a mixed-effects calibration. This is because the mixed-effects statistical model better characterizes the underlying parameter distributions than the fixed effects model. The results suggest that validation metrics can be used to select the most accurate calibration model for a particular empirical model with corresponding experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Robust Solution Verification Experiments on Nonuniform Meshes

The activities of verification, validation, and uncertainty quantification (VVUQ) provide a comprehensive means to assess the credibility of computational models. Within VVUQ, solution verification assesses numerical errors and evaluates whether the simulation is sufficiently accurate for its intended applications. As computational modeling gains traction in the development of complex, high-consequence systems, the need for robust solution verification intensifies, particularly because experimental data for these systems are often limited. This work examines improvements in the robustness of Richardson extrapolation (RE), a method commonly used in solution verification to study the discretization error of computational models using a power law. Nonuniform mesh refinement is discussed alongside other pollutants that affect the robustness of the power law model. Maximum likelihood estimation (MLE) is proposed as a robust strategy to address the uncertainty generated by nonuniform mesh refinement. An exploratory computational fluid dynamics (CFD) study of a 2D planar Poiseuille flow is conducted to determine if nonuniform mesh noise can be modeled with this MLE approach for more robust RE.

Weinmeister, Justin [ORNL] (ORCID:0000000160090237↗
DOE OSTI2020

Computational micromechanics model based failure criteria for chopped carbon fiber sheet molding compound composites

Chopped carbon fiber sheet molding compound has a great potential in lightweight automotive, marine, and aerospace applications. One of the most challenging tasks is to predict the failure strength of the material due to its anisotropy and heterogeneity, as well as complex stress states in real-world working conditions. In this work, a novel constitutive model of carbon fiber chip is proposed to capture the pre- and post-failure behaviors under different loading modes. On this basis, we propose a new computational micromechanics model, which is calibrated and validated by uniaxial tensile, compressive, and in-plane shear experiments. Furthermore, a set of microstructures representative volume element (RVE) models under complex loading conditions are reconstructed to understand the relationship between the microstructure characteristics and the failure envelopes. Finally, several modified versions of classical failure criteria are proposed for anisotropic materials with consideration of the fiber orientation tensor. The modified Tsai-Wu failure criterion, which shows the best accuracy among all failure criteria, is highlighted in the comparative study.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A computational fluid dynamics model to estimate local quantities in firebrand char oxidation

Firebrand burning is a complex phenomenon that is influenced by several parameters which are difficult to fully explore experimentally. Computational fluid dynamics models capable of predicting local quantities are essential for accurate prediction of char oxidation in firebrands. This article presents a computational fluid dynamics model to estimate firebrand mass loss, diameter change, and surface temperature during char oxidation. The model was validated using previously conducted wind tunnel experiments. These experiments were conducted for firebrands of two different aspect ratios, which were arranged in three different configurations (single, horizontal array, and vertical array), and for four different wind speeds (0.5, 1, 1.5, and 2 m/s). The computational fluid dynamics results were compared with a previous 1 D model. In all the test cases, the computational fluid dynamics model predicted the physical phenomena with significantly improved accuracy compared to a 1 D model. The char oxidation model presented in this article can be coupled with other models to study firebrand generation and trajectory, biomass pyrolysis, fluidized bed reactors, and coal combustion.

Engineering↗
DOE OSTI2024

Innovative Natural Gas Technologies For Efficiency Gain In Reliable And Affordable Thermochemical Electricity Generation (INTEGRATE)

Hyper is a cutting-edge system that pairs computational models that simulate the response of different systems with physical components such as sensors and actuators to transfer said calculated response to the HyPer facility. Computational models must operate in real-time in order to resolve the system response when coupled to physical hardware. Real-time is defined by the minimal HyPer system response time of Δt = 0.080 s. The integration of these models serves as a foundation for performance characterization based on real-world SOFC hardware and paves the way for comprehensive fuel flexibility studies. This report underscores the innovative approach to cyber-physical system deployment at NETL’s HyPer Facility, highlighting the critical role of computational models in enhancing the performance and adaptability of energy technologies.

03 NATURAL GAS↗
DOE OSTI2024

Robust Combined Heat and Hybrid Power (CHHP) for High Electrical Efficiency Cogeneration

Georgia Tech (Prime Recipient), the University of Texas at El Paso (UTEP) and the National Energy Technology Laboratory (NETL) investigated a hybrid fuel cell/ gas turbine system concept as a combined heat and hybrid power (CHHP) system for both robust and high power-to-process heat ratio cogeneration. The novelty of the proposed system entailed the distinct, elevated electrical efficiencies it maintains while simultaneously supporting a broad span of heating needs (e.g., supply temperatures) demanded across variable heat loads. The scope included: 1) leveraging a pioneering national lab facility configured for dynamic system operability development of hybrid fuel cell/gas turbine cycles; 2) enabling technology development to adjust and modulate the quality and quantity of thermal supply to bottoming heat loads via novel extreme temperature gas bypass valves. Hybrid fuel cell/gas turbine systems have primarily been reduced-to-practice in a constrained (e.g., initial proof-of-concept) manner and have still demonstrated considerable electrical efficiencies. However, these pre-pilot systems have focused upon electrical efficiencies and electrical power generation as the exclusive energy demand. Such hybrid systems had not been extensively researched or developed for flexible and variable operation consisting of both power and heat demands; however, these variable combined power and heat demands are characteristic of many types of manufacturers such as animal/poultry processing, bakeries and milk/flour/pastry manufacturing, textile mills, and electrochemical processing. Commercially, developing the system into a working combined heat and power system benefits these types of manufacturers by allowing them to meet their power and heat demands at a lower cost, higher efficiency, and/or through onsite generation. Therefore, the technical scope of this project was largely to study and facilitate these hybrid systems as combined heat and hybrid power (CHHP) systems that include dynamic operability for variable heat and power loads and/or grid dynamics for various types of manufacturers. Simulation results were used to predict the performance of the CHHP system and conceptually develop it to achieve desired dynamic operability. Experimentally, the primary goal was to design, manufacture, and experiment upon a high-temperature bypass valve. Experimental data included air mass flow rates through the valve orifice when the valve was changed to variable extent between fully closed and fully open. The experimental data was then used to create a semi-empirical computational model of the bypass valve. Concluded simulation goals for the research included developing computational heat exchanger models for the hybrid system inclusive of the bottoming heat exchanger and the recuperative heat exchanger, and then combining the computational recuperator model with the computational valve model. Afterwards, the computational models were then integrated to predict the dynamic operation of hybrid fuel cell/gas turbine cycles throughout a design space and reporting such. The scope stated in the preceding paragraph was packaged into five specific goals: 1) enabling the simulation of dynamic combined heat and power through the creation of computational, modular heat exchanger models; 2) simulation and exploration of the CHHP system’s performance by integrating the heat exchanger models with the national lab’s pre-existing hybrid system (computational) simulation, but without the recuperator bypass valve concept in order to initially determine how the (baseline) system behaves and can be controlled in order to meet variable heat and power demands; 3) development and initial deployment of the high-temperature recuperator bypass valve technology in order to confirm and characterize the approach; 4) usage of the experimental data for the valve to create a semi-empirical computational model for the bypass valve which could then be combined with the heat exchanger computational models; 5) repeat of the second task of simulating and exploring the system’s performance, but this time including the bypass valve to resolve its efficacy. Tasks were successfully completed, and the general notion of flexibly operating, high electrical efficiency CHHP was further corroborated. Supportive details are provided in the report.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Understanding the Computing and Analysis Needs for Resiliency of Power Systems from Severe Weather Impacts

As the frequency and intensity of severe weather has increased, its effect on the electric grid has manifested in the form of significantly more and larger outages in the United States. This has become especially true for regions that were previously isolated from weather extremes. In this paper, we analyze the weather impacts on the electric power grid across a variety of weather conditions, draw correlations, and provide practical insights into the operational state of these systems. High resolution computational modeling of specific meteorological variables, computational approaches to solving power system models under these conditions, and the types of resiliency needs are highlighted as goal-oriented computing approaches are being built to address grid resiliency needs. An example analysis correlating outages to 1km day-ahead weather from two historical winter storms, calculated on a large cluster using a combination of interpolated and extrapolated inputs from multiple instrumented sites to workflows that produce primary meteorological outputs, is shown as initial proof of concept.

analysis↗
DOE OSTI2022

Gamma Radiation-Induced Degradation of Acetohydroxamic Acid (AHA) in Aqueous Nitrate and Nitric Acid Solutions Evaluated by Multiscale Modelling

Acetohydroxamic acid (AHA) has been proposed for inclusion in advanced, single-cycle, used nuclear fuel reprocessing solvent systems for the reduction and complexation of plutonium and neptunium ions. For this application, a detailed description of the degradation of AHA in dilute aqueous nitric acid is required. To this end, we present a comprehensive, multiscale computer model for the coupled radiolytic and hydrolytic degradation of AHA in aqueous sodium nitrate and nitric acid solutions. Rate coefficients for the reactions of AHA and HA with the oxidizing nitrate radical were measured for the first time using pulsed electron radiolysis and used as inputs in the kinetic model. The computer model results are validated by comparison to experimental data from steady state gamma irradiations, for which the agreement is excellent. Here, the presented model accurately predicts the yields of the major degradation products of AHA: acetic acid, hydroxylamine, nitrous oxide, and molecular hydrogen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

UQpy: A general purpose Python package and development environment for uncertainty quantification

In this paper, we present the UQpy software toolbox, an open-source Python package for general uncertainty quantification (UQ) in mathematical and physical systems. The software serves as both a user-ready toolbox that includes many of the latest methods for UQ in computational modeling and a convenient development environment for Python programmers advancing the field of UQ. The paper presents an introduction to the software's architecture and existing capabilities, divided in the code in a set of modules centered around different UQ tasks such as sampling methods, generation of random processes and random fields, probabilistic inverse modeling, reliability analysis, surrogate modeling, and active learning. The paper also highlights the importance of the RunModel module, which is used to drive simulations in the uncertainty analyses performed in UQpy. This module conveniently allows the user to define computational models directly in Python, or to run simulations from a third-party software in serial or in parallel. To illustrate the various capabilities, two examples are tracked throughout the paper and analyzed repeatedly for various UQ tasks. The first is a Python model solving a nonlinear structural dynamics problem, used to illustrate UQpy's capabilities in sampling and forward propagation of high dimensional random vectors (stochastic processes), and probabilistic inference. The second model is a third-party Abaqus finite element model solving the thermomechanical response of a beam structure. This example is used to illustrate UQpy's capabilities in variance reduction sampling techniques, reliability analysis, surrogate modeling and active learning techniques.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Flow instabilities in helical-coil steam generators for small modular reactors: A review

Here, this study covers the research and discoveries in two-phase flow-boiling instabilities available in the literature—specifically for a helical-coil steam generator (HCSG), including experimental findings, theoretical research, computational models, and system code analyses—supporting research and development of representative small modular reactors (SMRs). Like other new and advanced reactor systems, water-cooled SMRs require experimental data from both integral and separate thermal-hydraulics test facilities for the verification and validation (V&V) of the computational models and computer codes in order to design and obtain regulatory approval. The complex dynamics of two-phase flow-boiling instabilities includes flow regimes physics phenomena, flow-channel geometries, heat-transfer behavior, and interactions among the solid–liquid-gas within the system boundary, all of which are pivotal for understanding the design and operational challenges of SMRs. This study focuses on identifying the relevant knowledge gaps on boiling instabilities—specifically for a HCSG—and provides insights about future research direction optimizing the transport of thermal energy, mass-flow rates, and boundary conditions that ensure the adequate heat-transfer performance, operational stability, and safety associated with SMR systems.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

Methods and devices for preventing computationally explosive calculations in a computer for model parameters distributed on a hierarchy of geometric simplices

A computer-implemented method of preventing computationally explosive calculations. The method includes obtaining, by a processor of the computer, measured data of one of a physical process or a physical object; performing hierarchical numerical modeling of a physical process inclusive of an Earth model containing at least one of (a) infrastructure in the ground and (b) a formation feature in the ground, wherein predicted data is generated; comparing the measured data to the predicted data to calculate an estimated error; analyzing the estimated error via an inversion process to update the at least one of the Earth model and infrastructure model so as to reduce the estimated error and to determine a final composite Earth model of at least one of the infrastructure and the feature; and using the final composite Earth model to characterize at least one of the process and the physical object.

Source record↗
DOE OSTI2023

RCSB Protein Data Bank: visualizing groups of experimentally determined PDB structures alongside computed structure models of proteins

Recent advances in Artificial Intelligence and Machine Learning (e.g., AlphaFold, RosettaFold, and ESMFold) enable prediction of three-dimensional (3D) protein structures from amino acid sequences alone at accuracies comparable to lower-resolution experimental methods. These tools have been employed to predict structures across entire proteomes and the results of large-scale metagenomic sequence studies, yielding an exponential increase in available biomolecular 3D structural information. Given the enormous volume of this newly computed biostructure data, there is an urgent need for robust tools to manage, search, cluster, and visualize large collections of structures. Equally important is the capability to efficiently summarize and visualize metadata, biological/biochemical annotations, and structural features, particularly when working with vast numbers of protein structures of both experimental origin from the Protein Data Bank (PDB) and computationally-predicted models. Moreover, researchers require advanced visualization techniques that support interactive exploration of multiple sequences and structural alignments. This paper introduces a suite of tools provided on the RCSB PDB research-focused web portal RCSB. org, tailor-made for efficient management, search, organization, and visualization of this burgeoning corpus of 3D macromolecular structure data.

3D visualization↗