Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model based development”

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 OSTI2021

A Scoping Study to Develop a Computational Fluid Dynamics Based Model to Predict Radiological Materials Packaging Temperatures within a Generic Staging Building

The objective of this work is to perform computational fluid dynamics simulations of a ventilated radiological-material-package staging building to determine the effect of including or excluding a variety of physical effects and computational methods on both the predicted package temperatures and required computational resources. The generic building contains 640 drum-packages containing heat-generating, radiological material supported on four racks that are eight levels tall. Additionally, there is a forced ventilation system, lighting, and insulated walls. Computational models were constructed that included or excluded (a) shelving, (b) effects of unsteadiness, and (c) radiation heat transfer. Simulations with each drum modeled separately were compared to simpler simulations with sets of four drums represented by an equivalent box-package. These models employed between 106 to 107 elements. Steady state simulations predicted package temperatures that were within 0.1°C of simulations that included transient effects and required only one-eighth the computational resources. Calculations that excluded shelving predicted temperatures within 0.6°C of simulations that included shelving and required one-fourth the computational resources. Excluding radiation heat transfer systematically increased temperatures by around 1.5°C but reduced computational resources by a factor of four. The box-package model reduced the computational resources by a factor of 3, but systematically predicted higher temperatures by around 1.1°C. These results will be used to develop an efficient computational fluid dynamics model to assess the ability of different staging building designs to prevent the temperature of package components from exceeding specified limits.

42 ENGINEERING↗
DOE OSTI2023

Validation of IMEP on Alcator C-Mod and JET-ILW ELMy H-mode plasmas

Abstract The recently developed integrated model based on engineering parameters (IMEP) (Luda et al 2020 Nucl. Fusion 61 126048; Luda et al 2021 Nucl. Fusion 60 036023), so far validated on ASDEX Upgrade, has been tested on a database of 3 Alcator C-Mod and 55 JET-ILW ELMy (type I) H-mode stationary phases. The empirical pedestal transport model included in IMEP, consisting now of imposing a fixed value of R < ∇ T e > / T e , t o p = − 82.5 , allows an accurate prediction of the pedestal top temperature (when the pedestal top density is fixed to the experimental measurements) across these three machines with different sizes, when the pedestal is peeling–ballooning (PB) limited. Cases far from the ideal PB boundary, corresponding to high edge Spitzer resistivity, are instead strongly overpredicted by IMEP. A comparison between the predictions of Europed and IMEP for a subset of JET-ILW cases shows that IMEP can more accurately reproduce the experimental pedestal width. This allows IMEP to better capture profile effects on the pedestal stability, and therefore to correctly describe the negative effect of fueling on the pedestal pressure for PB limited cases. A strong correlation between the separatrix density and the fueling rate has been identified for a subset of JET-ILW cases, when taking into account different divertor configurations. Overall, these promising results encourage further developments of integrated models to obtain reliable predictions of pedestal and global confinement using only engineering parameters for present and future machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Design and operational considerations of catalytic membrane reactors for ammonia synthesis

Production of ammonia using hydrogen derived from renewable electricity instead of hydrocarbon reforming would dramatically reduce the carbon footprint of this commodity chemical. Novel technologies such as catalytic membrane reactors may potentially be more compatible with distributed ammonia production than the conventional Haber-Bosch process. Here, a reactor model is developed based on integrating a standard industrial iron catalyst into a catalytic membrane reactor (CMR) equipped with an inorganic membrane that is selective to NH 3 over N 2 /H 2 . CMR performance is studied as functions of wide ranges of membrane properties and operating conditions. Conversion and ammonia recovery are dictated principally by the ammonia permeance, and the benefits by using membranes become significant above 100 GPU = 3.4 × 10 –8 mol m –2 s –1 Pa –1 . To be effective, the CMR requires a minimum selectivity for ammonia of 10 over both nitrogen and hydrogen, and purity scales with the effective selectivity. Increasing the pressure of operation significantly improves all metrics, and at P = 30 bar with a quality membrane ammonia is almost completely recovered, enabling direct recycle of un-reacted hydrogen and nitrogen without need for recompression. Temperature drives conversion and scales monotonically without thermodynamic limitations in a CMR. Alternatively, the temperature may be reduced as low as 300°C while achieving conversion levels surpassing equilibrium limits at T = 400°C in a conventional reactor.

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

Leveraging operational information from wastewater injection wells to evaluate CO 2 injection performance for carbon storage applications in the Appalachian Basin

Abstract Geologic parameters, geophysical logging, injection testing, and operational metrics from wastewater injection wells were integrated to develop a preliminary design of a carbon storage facility in the Appalachian Basin. A scattered group of 10–20 commercial wastewater injection wells dispose off produced water from oil and gas wells in the region, utilizing a sequence of stacked deep saline formations for injection zones. These wastewater injection wells provide practical benchmarks for understanding the feasibility of carbon dioxide (CO 2 ) storage. Geologic models were developed based on characterization data from the wastewater injection wells. Reservoir simulations were calibrated according to injection testing and operational data from the wastewater injection wells. Long‐term operational data on injection flow rates and pressures measured in the wastewater injection wells were especially useful to evaluate the performance of carbon storage applications. The simulations were used to estimate injection pressures, radius of CO 2 saturation, and pressure response for industrial scale CO 2 storage applications. Results were also used to provide a design basis in terms of number of injection wells, well spacing, area of review, injection system components, monitoring plan, and CO 2 pipeline distribution system. The analysis demonstrates that there is sufficient injectivity in the deep saline formations in the west‐central Appalachian Basin to store commercial volumes of anthropogenic CO 2 . The geologic system appears suitable for supporting CO 2 injection rates of 0.5–1.0 million metric tons per year at injection pressures below formation fracture pressure in a single well. The long‐term operational data of wastewater injection wells within the study area suggested a lower permeability‐thickness values than indicated by initial reservoir tests. A workflow for developing realistic permeability values for input into reservoir simulations is presented. © 2020 Society of Chemical Industry and John Wiley & Sons, Ltd.

Valluri, Manoj↗
DOE OSTI2024

TATB thermal decomposition: Expanding the molecular profile with cryo-focused pyrolysis GC-MS

Understanding the molecular composition of high explosives during thermal decomposition is vital for predicting the sensitivity, safety, and performance of explosive materials. The thermal decomposition of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) has been linked to the formation of furazans through a series of dehydration reactions of the NO2 and NH2 groups on the phenyl ring, along with breakdown into small molecules (≤120 amu). Molecular identification of compounds formed in this transformation of the furazans to light gases has been lacking. To address this, we have applied a pseudo-confined sampling system in a cryo-focused pyrolysis gas chromatography-mass spectrometry (pyGC-MS) system to molecularly identify these intermediates. By design, sublimation of TATB, which has complicated MS analyses of thermal degradation, was significantly reduced and additional compounds were identified with potential structural information. In addition to the known furazan compounds, one of these compounds forms from the loss of oxygen from benzo-trifurazan (F3) and produces an open ring structure that may be the first step in the formation of lower molecular weight furazan breakdown products. The loss of a nitro group from benzo-monofurazan (F1) was also discovered and implicates the formation of oxidizing NO2 gas in the thermal decomposition mechanism. So these findings are vital for understanding the proper heat flow from energetic materials on a molecular level, necessary when measuring enthalpy and developing decomposition models based on kinetic parameters.

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

Development of a conduction-based model for analyzing frozen startup of alkali-metal heat pipes

One key area of interest in heat pipe modeling/simulation is to analyze the startup behavior of the liquid-metal heat pipes (LMHPs) from a frozen state. This so-called ‘frozen startup’ process involves a complex set of nonlinear mass and heat transport phenomena, including phase transitions from solid to liquid and vapor, multiphase interactions, microporous wick flow, and compressible vapor dynamics. The complexity of these processes makes it challenging to simulate LMHP’s frozen startup using conventional numerical methods or commercial computational fluid dynamics (CFD) software. This paper presents a simplified conduction-based modeling approach that can provide practical insights into the entire LMHP frozen startup process, while alleviating the challenges of modeling its complex physics. The theoretical foundation and physical assumptions of the proposed model are based solely on heat-conduction equation, allowing for a more tractable simulation without sacrificing essential physical accuracy. The proposed model was implemented in a commercial CFD software, and its prediction was compared with the experimental data obtained from sodium heat-pipe startup experiments. The comparison highlights the proposed model's ability to capture the transient thermal behavior of LMHP during frozen startup. This study not only validates the conduction-based frozen startup modeling method but also shows its potential as a practical and efficient tool for understanding the startup performance of the LMHP systems.

Microreactor↗
DOE OSTI2016

Development of a Low-Level Ar-37 Calibration Standard

Argon-37 is an important environmental signature of an underground nuclear explosion. Producing and quantifying low-level 37Ar standards is an important step in the development of sensitive field measurement instruments for use during an On-Site Inspection, a key provision of the Comprehensive Nuclear-Test-Ban Treaty. This paper describes progress at Pacific Northwest National Laboratory (PNNL) in the development of a process to generate and quantify low-level 37Ar standard material, which can then be used to calibrate sensitive field systems at activities consistent with soil background levels. The 37Ar used for our work was generated using a laboratory-scale, high-energy neutron source to irradiate powdered samples of calcium carbonate. Small aliquots of 37Ar were then extracted from the head space of the irradiated samples. The specific activity of the head space samples, mixed with P10 (90% stable argon:10% methane by mole fraction) count gas, is then derived using the accepted Length-Compensated Internal-Source Proportional Counting method. Due to the low activity of the samples, a set of three Ultra-Low Background Proportional-Counters designed and fabricated at PNNL from radio-pure electroformed copper was used to make the measurements in PNNL’s shallow underground counting laboratory. Very low background levels (<10 counts/day) have been observed in the spectral region near the 37Ar emission feature at 2.8 keV. Two separate samples from the same irradiation were measured. The first sample was counted for 12 days beginning 28 days after irradiation, the second sample was counted for 24 days beginning 70 days after irradiation (the half-life of 37Ar is 35.0 days). Both sets of measurements were analyzed and yielded very similar results for the starting activity (~0.1 Bq) and activity concentration (0.15 mBq/ccSTP argon) after P10 count gas was added. A detailed uncertainty model was developed based on the ISO Guide to the Expression of Uncertainty in Measurement. This paper presents a discussion of the measurement analysis, along with assumptions and uncertainty estimates.

Absolute radiation measurement↗
DOE OSTI2022

Machine learning for surrogate process models of bioproduction pathways

Technoeconomic analysis and life-cycle assessment are critical to guiding and prioritizing bench-scale experiments and to evaluating economic and environmental performance of biofuel or biochemical production processes at scale. Traditionally, commercial process simulation tools have been used to develop detailed models for these purposes. However, developing and running such models can be costly and computationally intensive, which limits the degree to which they can be shared and reproduced in the broader research community. This study evaluates the potential of an automated machine learning approach to develop surrogate models based on conventional process simulation models. The analysis focuses on several high-value biofuels and bioproducts for which pathways of production from biomass feedstocks have been well-established. The results demonstrate that surrogate models can be an accurate and effective tool for approximating the cost, mass and energy balance outputs of more complex process simulations at a fraction of the computational expense.

09 BIOMASS FUELS↗
DOE OSTI2022

Magnetic anisotropy in permalloy antidot square lattice

Here, magnetic anisotropy of Permalloy (Py) antidot square lattice was investigated by torquemetry method using Rotation Magneto-Optic Kerr Effect (ROTMOKE). We find that there exists a field-dependent 4-fold magnetic anisotropy with the easy magnetization axis along the axis of the antidot square lattice. In addition, there also exists an artifact of a uniaxial magnetic anisotropy in ROTMOKE result. We show that both results are due to the period wiggling of the magnetization in space which was confirmed by magnetic imaging using magnetic transmission soft x-ray microscopy (MTXM). Micromagnetic simulation from MuMax3 supports the wiggling structure of the magnetization, as well as reproduces ROTMOKE result. A simplified model was developed based on the periodic wiggling of the magnetization and successfully explored the physical origin of the field-dependent 4-fold anisotropy and the artifact of the uniaxial anisotropy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Applicability of discrete element method with spherical and clumped particles for constitutive study of granular materials

Discrete element method (DEM) has been intensively used to study the constitutive behaviour of granular materials. However, to what extent a real granular material can be reproduced by virtual DEM simulations remains unclear. This study attempts to answer this question by comparing DEM simulations with typical features of experimental granular materials. Three groups of models with spherical and clumped particles are investigated from four perspectives: (i) deviatoric stress and volumetric behaviour; (ii) critical state behaviour; (iii) stress-dilatancy relationship; and (iv) the evolution of principal stress ratio against axial strain. The results demonstrate that DEM with spherical or clumped particles is capable of qualitatively describing macroscopic deviatoric stress responses, volumetric behaviour, and critical state behaviour observed in experiments for granular materials. On the other hand, some qualitative deviations between experiments and the investigated DEM simulations are also observed, in terms of the stress-dilatancy behaviour and principal stress ratio against axial strain, which are proven to be critical for constitutive modelling. The results demonstrate that DEM with spherical or clumped particles may not necessarily fully capture experimental features of granular materials even from a qualitative perspective. It is thus encouraged to thoroughly validate DEM with experiments when developing constitutive models based on DEM observations.

58 GEOSCIENCES↗
DOE OSTI2015

Investigation of precipitate refinement in Mg alloys by an analytical composite failure model

An analytical model is developed to simulate precipitate refinement in second phase strengthened magnesium alloys. The model is developed based on determination of the stress fields inside elliptical precipitates embedded in a rate dependent inelastic matrix. The stress fields are utilized to determine the failure mode that governs the refinement behavior. Using an AZ31 Mg alloy as an example, the effects the applied load, aspect ratio and orientation of the particle is studied on the macroscopic failure of a single α-Mg17Al12 precipitate. Additionally, a temperature dependent version of the corresponding constitutive law is used to incorporate the effects of temperature. In plane strain compression, an extensional failure mode always fragments the precipitates. The critical strain rate at which the precipitates start to fail strongly depends on the orientation of the precipitate with respect to loading direction. The results show that the higher the aspect ratio is, the easier the precipitate fractures. Precipitate shape is another factor influencing the failure response. In contrast to elliptical precipitates with high aspect ratio, spherical precipitates are strongly resistant to sectioning. In pure shear loading, in addition to the extensional mode of precipitate failure, a shearing mode may get activated depending on orientation and aspect ratio of the precipitate. The effect of temperature in relation to strain rate was also verified for plane strain compression and pure shear loading cases.

Tabei, Ali↗
DOE OSTI2025

Recrystallization driven softening and heating rate dependencies of FeCrAl nuclear fuel cladding during accident transients

A refined understanding of FeCrAl cladding behavior during rapid transients is critical for its potential deployment in light-water reactors. Current assessments focus on transient burst testing metrics such as balloon geometry, burst temperature, and hoop stress, often used as proxies for simpler conventional tensile properties. However, directly correlating isothermal tensile and creep data with accident transient scenarios remains a challenge, although it is essential for high fidelity model development. Recent modeling based on tensile tests up to 800 °C, conducted with both immediate loading and a 10-minute soak, showed that immediate loading better predicts experimental burst temperatures, indicating a thermal softening effect. Building upon this observation, the current study connects transient performance, microstructural evolution, and high-temperature tensile properties by leveraging results from C26M claddings burst tests performed at heating rates of 1–50 °C/s and hoop stresses from 25 to 100 MPa. At 25 MPa, rupture temperatures varied by only 6 °C, but at 100 MPa, the difference reached 116 °C, with faster heating yielding higher burst temperatures. Microstructural analysis identified recrystallization as the primary cause of heating rate-dependent softening, eliminating prior cold-working. In-situ thermomechanical data linked ballooning onset to localized instabilities, similar to ultimate tensile strength behavior in conventional tensile tests. High heating rates correlated with immediate loading tensile data, while lower rates matched soaked data. Furthermore, by linking burst performance to microstructural evolution and tensile properties, this work provides a foundation for more accurate modeling of FeCrAl claddings and potentially other Fe-based materials under accident conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2021

Pressure Effects on the Relaxation of an Excited Ethane Molecule in High-Pressure Bath Gases

Here, we use molecular dynamics to calculate the rotational and vibrational energy relaxation of C 2 H 6 in Ar, Kr, and Xe bath gases over a pressure range of 10 to 400 atm and at temperatures of 300 K and 800 K. The C 2 H 6 is instantaneously excited by 80 kcal/mol randomly distributed into both vibrational and rotational modes. The computed relaxation rates show little sensitivity to the identity of the noble gas in the bath. Vibrational relaxation rates show a non-linear pressure dependence at 300 K. At 800 K the reduced range of bath gas densities covered by the range of pressures do not yet show any non-linearity in the pressure dependence. Rotational relaxation is characterized with two relaxation rates. The slower rate is comparable to the vibrational relaxation rate. The faster rate has a linear pressure dependence at 300 K but an irregular, nonlinear pressure dependence at 800 K. To understand this, a model was developed based on approximating the periodic box used in the molecular dynamics simulations by an equal-volume collection of cubes where each cube is sized to allow only single occupancy by the noble gas or the molecule. Combinatorial statistics then leads to a pressure and temperature dependent analytic distribution of the bath gas species the molecule encounters in a collision. This distribution, the dissociation energy of molecule/bath gas complexes and bath gas clusters, and the computed energy release per collision combine to show that only at 300 K is the energy release sufficient to dissociate likely complexes and clusters. This suggests that persistent and pressure-dependent clusters and complexes at 800 K may be responsible for the non-linear pressure dependence of rotational relaxation.

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

Observation and rationalization of nitrogen oxidation enabled only by coupled plasma and catalyst

Abstract Heterogeneous catalysts coupled with non-thermal plasmas (NTP) are known to achieve reaction yields that exceed the contributions of the individual components. Rationalization of the enhancing potential of catalysts, however, remains challenging because the background contributions from NTP or catalysts are often non-negligible. Here, we first demonstrate platinum (Pt)-catalyzed nitrogen (N 2 ) oxidation in a radio frequency plasma afterglow at conditions at which neither catalyst nor plasma alone produces significant concentrations of nitric oxide (NO). We then develop reactor models based on reduced NTP- and surface-microkinetic mechanisms to identify the features of each that lead to the synergy between NTP and Pt. At experimental conditions, NTP and thermal catalytic NO production are suppressed by radical reactions and high N 2 dissociation barrier, respectively. Pt catalyzes NTP-generated radicals and vibrationally excited molecules to produce NO. The model construction further illustrates that the optimization of productivity and energy efficiency involves tuning of plasma species, catalysts properties, and the reactor configurations to couple plasma and catalysts. These results provide unambiguous evidence of synergism between plasma and catalyst, the origins of that synergy for N 2 oxidation, and a modeling approach to guide material selection and system optimization.

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

Highest transmittance and high-mobility amorphous indium gallium zinc oxide films on flexible substrate by room-temperature deposition and post-deposition anneals

Amorphous indium gallium zinc oxide (a-IGZO) thin films of the highest transmittance reported in literature were initially deposited onto flexible polymer substrates at room temperature. The films were annealed in vacuum, air, and oxygen to enhance their electrical and optical performances. Electrical and optical characterizations were done before and after anneals. A partial reversal of the degradation in electrical properties upon annealing in oxygen was achieved by subjecting the films to subsequent vacuum anneals. A model was developed based on film texture and structural defects which showed close agreement between the measured and calculated carrier mobility values at low carrier concentrations (2-6 x 10{sup 19} cm{sup -3}).

36 MATERIALS SCIENCE↗
DOE OSTI2013

Phonon thermal transport outside of local equilibrium in nanowires via molecular dynamics

We study thermal transport through Pt nanowires that bridge planar contacts as a function of wire length and vibrational frequency of the contacts. When phonons in the contacts have lower average frequencies than those in the wires thermal transport occurs under conditions away from local equilibrium with low-frequency phonons experiencing a higher thermal gradient than high-frequency ones. This results in a size-dependent increase in the effective thermal conductivity of the wire with decreasing vibrational frequencies of the contacts. The interfacial resistivity when heat flows from the wire to the contact is also size-dependent and has the same physical origin in the lack of full equilibration in short nanowires. We develop a model based on a 1D atomic chain that captures the salient physics of the MD results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗