Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer modeling”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

A Computer Model to Study the Response of Energetic Materials to a Range of Dynamic Loads

Over the past decade we developed a model to enable computer simulation of the mechanical and subsequent energetic response of explosives and propellants to mechanical insults such as impacts, perforations, drops, and falls. The model is embedded in computer simulation programs that solve the non-linear, large deformation equations of compressible solid and fluid flow in space and time. It is implemented as a user-defined model, which returns the updated stress tensor and composition that result from the simulation supplied strain tensor change. Although it is multi-phase, in that gas and solid species are present, it is single-velocity, in that the gas does not flow through the porous solid. More than 70 time-dependent variables are made available for additional analyses and plotting. The model encompasses a broad range of possible responses: mechanical damage with no energetic response, and a continuous spectrum of degrees of violence including delayed and prompt detonation. This work describes the basic workings of the model.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Cosmic logic: a computational model

We initiate a formal study of logical inferences in context of the measure problem in cosmology or what we call cosmic logic. We describe a simple computational model of cosmic logic suitable for analysis of, for example, discretized cosmological systems. The construction is based on a particular model of computation, developed by Alan Turing, with cosmic observers (CO), cosmic measures (CM) and cosmic symmetries (CS) described by Turing machines. CO machines always start with a blank tape and CM machines take CO's Turing number (also known as description number or Gödel number) as input and output the corresponding probability. Similarly, CS machines take CO's Turing number as input, but output either one if the CO machines are in the same equivalence class or zero otherwise. We argue that CS machines are more fundamental than CM machines and, thus, should be used as building blocks in constructing CM machines. We prove the non-computability of a CS machine which discriminates between two classes of CO machines: mortal that halts in finite time and immortal that runs forever. In context of eternal inflation this result implies that it is impossible to construct CM machines to compute probabilities on the set of all CO machines using cut-off prescriptions. The cut-off measures can still be used if the set is reduced to include only machines which halt after a finite and predetermined number of steps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

GMRES with embedded ensemble propagation for the efficient solution of parametric linear systems in uncertainty quantification of computational models

In a previous work, embedded ensemble propagation was proposed to improve the efficiency of sampling-based uncertainty quantification methods of computational models on emerging computational architectures. It consists of simultaneously evaluating the model for a subset of samples together, instead of evaluating them individually. A first approach introduced to solve parametric linear systems with ensemble propagation is ensemble reduction. In Krylov methods for example, this reduction consists in coupling the samples together using an inner product that sums the sample contributions. Ensemble reduction has the advantages of being able to use optimized implementations of BLAS functions and having a stopping criterion which involves only one scalar. However, the reduction potentially decreases the rate of convergence due to the gathering of the spectra of the samples. In this paper, we investigate a second approach: ensemble propagation without ensemble reduction in the case of GMRES. This second approach solves each sample simultaneously but independently to improve the convergence compared to ensemble reduction. This raises two new issues which are solved in this paper: the fact that optimized implementations of BLAS functions cannot be used anymore and that ensemble divergence, whereby individual samples within an ensemble must follow different code execution paths, can occur. We tackle those issues by implementing a high-performing ensemble GEMV and by using masks. The proposed ensemble GEMV leads to a similar cost per GMRES iteration for both approaches, i.e. with and without reduction. For illustration, we study the performances of the new linear solver in the context of a mesh tying problem. Furthermore, this example demonstrates improved ensemble propagation speed-up without reduction.

BLAS↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
DOE OSTI2026

A fast and robust computational modeling approach for density and shape predictions in powder metallurgy hot isostatic pressing

Powder metallurgy hot isostatic pressing (PM-HIP) is an advanced manufacturing process that produces near-net-shape parts with high material utilization and uniform microstructures. PM-HIP is frequently used for producing small-scale parts with complicated geometries and is potentially economical for producing large-scale parts. However, excessive post-HIP shape distortions can reduce its effectiveness and economic advantage, especially for larger parts. A PM-HIP computational model can predict and help mitigate these distortions. However, due to complex deformation mechanisms and thermo-mechanical coupling present in PM-HIP processes, these non-linear computational models sometimes become numerically unstable. The numerical instabilities in these models can lead to very slow convergence or no convergence at all, which often translates to slow and unreliable models. These limitations are more pronounced in large models with complicated geometries. Hence, in this work, an alternative modeling approach is presented that improves numerical stability and computational performance. The presented approach achieves these improvements through approximating the fully coupled thermo-mechanical PM-HIP model as a decoupled model and adding inertial damping to the model’s mechanical part. In conclusion, a comparison with the fully coupled model indicated a slight dip in prediction accuracy (<5% error) but significant improvements in numerical stability (>20 times larger time step size) and computational performance (5-10 times speed-up with less computational resource usage) when using the presented approach.

Hot isostatic pressing↗
DOE OSTI2017

Web-Based Computational Model Builder For Nuclear Energy Advanced Modeling and Simulation (Final Report)

According to the World Nuclear Association, the U.S. is the largest producer of nuclear power worldwide. In fact, the U.S. provides over 30 percent of the world's nuclear generation of electricity. Although few new reactors have been built in the past 30 years, the association anticipates that four to six reactors may be built by 2020. The first of the reactors will be built in response to 16 license applications, which have been completed since 2007, to build 24 additional nuclear reactors. Currently, both government and industry are working to accelerate approval for the construction and design of new nuclear energy plants. The U.S. nuclear industry has already achieved significant developments in the implementation of nuclear power plants due to advancements in refueling, maintenance, and safety systems at current power plants. Meanwhile, changes in government policy, which have occurred since the late 1990s, have provided for noteworthy expansion in nuclear capacity. For example, the Energy Policy Act of 2005 evoked investment in electricity infrastructure such as nuclear power. To design future systems and continue operational improvements at existing U.S. plants, advanced modeling and simulation of nuclear reactors are crucial. We developed a Cloud/Web-based tool, HPCCloud, that enables nuclear energy firms, or other manufacturing and engineering companies, to innovate using advanced modeling and simulation. Through the Web-Based Computational Model Builder for Nuclear Energy Advanced Modeling and Simulation grant, we develop HPCCloud to simplify and democratize advanced modeling and simulation in the nuclear energy industry on a range of nuclear engineering applications. It leverages millions of investment dollars from the Department of Energy's Consortium for Advanced Simulation of Light Water Reactors (CASL) and the Office of Nuclear Energy's research and development in Nuclear Energy Advanced Modeling and Simulation (NEAMS). Our approach simplifies the workflow, alleviates the need for in-house computational science and engineering experts, and lowers the capital investments required for advanced modeling and simulation. The Cloud/Web-based simulation environment guides end-users and developers through the advanced modeling and simulation lifecycle. In addition to providing significantly improved, intuitive software, the environment offers reproducible workflows where one can save, share, and even publish the full path of data from input to final analyzed results. HPCCloud addresses identified challenges to the nuclear energy industry's use of advanced modeling and simulation. With this unique Cloud/Web-based environment, we enable an ecosystem of synergistic activities, fueling new research and development for nuclear energy. HPCCloud lays the groundwork for commercialization, through customization, technology integration, and advanced research and development services.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Assessment of radiation dose in nuclear cardiovascular imaging using realistic computational models

Purpose: Nuclear cardiology plays an important role in clinical assessment and has enormous impact on the management of a variety of cardiovascular diseases. Pediatric patients at different age groups are exposed to a spectrum of radiation dose levels and associated cancer risks different from those of adults in diagnostic nuclear medicine procedures. Therefore, comprehensive radiation dosimetry evaluations for commonly used myocardial perfusion imaging (MPI) and viability radiotracers in target population (children and adults) at different age groups are highly desired. Methods: Using Monte Carlo calculations and biological effects of ionizing radiation VII model, we calculate the S-values for a number of radionuclides (Tl-201, Tc-99m, I-123, C-11, N-13, O-15, F-18, and Rb-82) and estimate the absorbed dose and effective dose for 12 MPI radiotracers in computational models including the newborn, 1-, 5-, 10-, 15-yr-old, and adult male and female computational phantoms. Results: For most organs, {sup 201}Tl produces the highest absorbed dose whereas {sup 82}Rb and {sup 15}O-water produce the lowest absorbed dose. For the newborn baby and adult patient, the effective dose of {sup 82}Rb is 48% and 77% lower than that of {sup 99m}Tc-tetrofosmin (rest), respectively. Conclusions: {sup 82}Rb results in lower effective dose in adults compared to {sup 99m}Tc-labeled tracers. However, this advantage is less apparent in children. The produced dosimetric databases for various radiotracers used in cardiovascular imaging, using new generation of computational models, can be used for risk-benefit assessment of a spectrum of patient population in clinical nuclear cardiology practice.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

FY-24 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the computational modeling of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY24. This year’s modeling campaign focuses on the continual testing of the RELAP5 model against experimental data. Several fault cases were simulated with the RELAP5 model in FY23 to study the accuracy of the model under complex flow conditions. Similarly, in FY24, a case is simulated where the chimney outlet is throttled with a valve while the flow is undergoing two-phase oscillations. The purpose of this case is to study the effects of increasing pressure drop in the two-phase region on the system behavior. Even though the model is able to capture the experimental data qualitatively, it overpredicts the pressure drops experienced by the flow in the two-phase regime. Nevertheless, the model and the experiment show that oscillations are stabilized with increased pressure drop in the two-phase region.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Computational modeling of proton acceleration with multi-picosecond and high energy, kilojoule, lasers

Here, we use computational modeling to investigate proton beam generation from kilojoule, multi-picosecond laser pulses pertinent to several recently commissioned, large-scale laser facilities. The dependencies of proton acceleration on electron source parameters including pulse duration, temperature, and flux are independently and systematically evaluated. Proton acceleration is found to depend not only on the source size and peak temperature of the injected electrons but also on the rate of increase for a more physical time-varying temperature. Simulations of a 10 ps, sub-relativistic intensity (8 × 10 17 W/cm 2 ) at 1 μm wavelength laser pulse show that energetic electrons generated within the expanding under-dense laser-produced plasma sustain the proton acceleration for ~20 ps. This results in 15 MeV energy gain of the protons, well above what would be predicted based on conventional intensity scalings or what has been observed with shorter pulses. Using this prolonged acceleration, a scheme consisting of a 1 ps and 10 ps double pulse is shown to further boost proton maximum energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Multi-scale computational modeling of residual stress in selective laser melting with uncertainty quantification

Selective laser melting (SLM) is a powder-based additive manufacturing technique which creates parts by fusing together successive layers of powder with a laser. The quality of produced parts is highly dependent on the proper selection of processing parameters, requiring significant testing and experimentation to determine parameters for a given machine and material. Computational modeling could potentially be used to shorten this process by identifying parameters through simulation. However, simulating complete SLM builds is challenging due to the difference in scale between the size of the particles and laser used in the build and the size of the part produced. Often, continuum models are employed which approximate the powder as a continuous medium to avoid the need to model powder particles individually. While computationally expedient, continuum models require as inputs effective material properties for the powder which are often difficult to obtain experimentally. Building on previous works which have developed methods for estimating these effective properties along with their uncertainties through the use of detailed models, this work presents a part scale continuum model capable of predicting residual thermal stresses in an SLM build with uncertainty estimates. Model predictions are compared to experimental measurements from the literature.

42 ENGINEERING↗
DOE OSTI2022

People who inject drugs in metropolitan Chicago: A meta-analysis of data from 1997-2017 to inform interventions and computational modeling toward hepatitis C microelimination

Progress toward hepatitis C virus (HCV) elimination in the United States is not on track to meet targets set by the World Health Organization, as the opioid crisis continues to drive both injection drug use and increasing HCV incidence. A pragmatic approach to achieving this is using a microelimination approach of focusing on high-risk populations such as people who inject drugs (PWID). Computational models are useful in understanding the complex interplay of individual, social, and structural level factors that might alter HCV incidence, prevalence, transmission, and treatment uptake to achieve HCV microelimination. However, these models need to be informed with realistic sociodemographic, risk behavior and network estimates on PWID. We conducted a meta-analysis of research studies spanning 20 years of research and interventions with PWID in metropolitan Chicago to produce parameters for a synthetic population for realistic computational models (e.g., agent-based models). We then fit an exponential random graph model (ERGM) using the network estimates from the meta-analysis in order to develop the network component of the synthetic population.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Mechanistic insights into metal, nitrogen doped carbon catalysts for oxygen reduction: progress in computational modeling

We report metal and nitrogen doped carbon materials (denoted as M–N–C) synthesized through high-temperature pyrolysis have been found to exhibit activity for oxygen reduction reaction (ORR) approaching that of Pt and electrochemical stability higher than previous MN 4 -containing macrocyclic molecular catalysts. Tremendous efforts have thus been devoted to the advancement of M–N–C catalysts as an economical alternative to Pt-based catalysts for proton exchange membrane fuel cell cathodes with a focus on simultaneously improving activity and stability. To this end, novel computational modeling techniques have been developed and applied to acquire knowledge crucial for accelerating the pace of M–N–C catalyst development. In this review, recent progress in computational method development, as well as the predictions of chemical structure of active sites, reaction pathways, ORR kinetics, and catalyst stability in electrochemical environments, are critically surveyed. Moreover, the crucial role of computational modeling to elucidate the functional mechanism of M–N–C catalysts for ORR in acid media and enable rational design of M–N–C catalysts is discussed with a visionary outlook for the field.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Computational modeling of nanoscale and microscale particle deposition, retention and dosimetry in the mouse respiratory tract

Comparing effects of inhaled particles across rodent test systems and between rodent test systems and humans is a key obstacle to the interpretation of common toxicological test systems for human risk assessment. These comparisons, correlation with effects and prediction of effects, are best conducted using measures of tissue dose in the respiratory tract. Differences in lung geometry, physiology and the characteristics of ventilation can give rise to differences in the regional deposition of particles in the lung in these species. Differences in regional lung tissue doses cannot currently be measured experimentally. Regional lung tissue dosimetry can however be predicted using models developed for rats, monkeys, and humans. A computational model of particle respiratory tract deposition and clearance was developed for BALB/c and B6C3F1 mice, creating a cross-species suite of available models for particle dosimetry in the lung. Airflow and particle transport equations were solved throughout the respiratory tract of these mice strains to obtain temporal and spatial concentration of inhaled particles from which deposition fractions were determined. Particle inhalability (Inhalable fraction, IF) and upper respiratory tract (URT) deposition were directly related to particle diffusive and inertial properties. Measurements of the retained mass at several post-exposure times following exposure to iron oxide nanoparticles, micro- and nanoscale C60 fullerene, and nanoscale silver particles were used to calibrate and verify model predictions of total lung dose. Interstrain (mice) and interspecies (mouse, rat and human) differences in particle inhalability, fractional deposition and tissue dosimetry are described for ultrafine, fine and coarse particles.

Lung deposition, nanosized particles, particle↗
DOE OSTI2019

A computational model for exploring particle acceleration during reconnection in macroscale systems

A new computational model suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems is presented. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfvén waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Review of Computational Models for the Flow of Milled Biomass Part I: Discrete-Particle Models

Biomass is a renewable and sustainable energy resource. Current design of biomass handling and feeding equipment leverage both experiments and numerical modeling. This paper reviews the state-of-the-art discrete element methods (DEM) for the flow of milled biomass (Part I), accompanied by a comprehensive review on continuum-based computational models (Part II). The present review on DEM is primarily focused on the features and suitability of various particle shape models for different types of milled biomass because particle shape is the predominant attribute controlling the flow behavior of complex-shaped granular material. The general strengths and weaknesses in the applicability of those models for the milled biomass modeling are summarized. In particular, comments are provided to balance the numerical model capabilities and the computational cost for the development of DEM models. To our best knowledge, this is the first-of-its-kind review on DEM specifically for biomass. Our study indicates that the current DEM models require further development, calibration, and validation based on a deep understanding of biomass particle contact mechanics and experimental data support before they can be reliably used for predictive simulations in handling and feeding systems.

42 ENGINEERING↗