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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of grassland fire dynamics on particulate emission factors

Understanding particulate emission factors in grassland fires is critical for improving air quality assessments and refining emissions inventories. This study explores the dynamic and heterogeneous nature of fire processes and their influence on particulate emissions. Traditional emission factor inventories often rely on static values, overlooking the spatially and temporally variable dynamics of fire behavior. Through field experiments conducted at Konza Prairie Biological Station measuring emissions with uncrewed aerial vehicles, and computational modeling using HIGRAD/FIRETEC, this research highlights the dependency of emission production on localized fire dynamics, including ignition patterns, fuel properties, and atmospheric conditions. The results demonstrate a need for fire dynamic considerations to accurately capture emissions variability.

09 BIOMASS FUELS↗
DOE OSTI2023

On the role of methyl groups in the molecular architectures of mesophase pitches

The role of methyl groups on the liquid–crystal structure of mesophase pitches was investigated by combining experimental characterizations and atomic-scale computational modeling, using three pitches synthesized from different precursors. Of the three pitches, C-9 alkyl benzene and naphthalene-based pitches have 13 and 7 methyl groups per average polyaromatic hydrocarbon, respectively. By contrast, mesophase produced from a coal-tar pitch has about one methyl group. The coal tar–based mesophase pitch is hydrogen deficient or more aromatic compared with C-9 alkyl benzene- and naphthalene-based pitches. Additionally, X-ray diffraction data showed that average coherent domain sizes of C-9 alkyl benzene (3.7 nm) and naphthalene-based (3.6 nm) pitches with more methyl groups are larger than that of coal tar–based mesophase (2.4 nm). Based on the identified features, the influence of the methyl group on the layering structures was investigated via molecular dynamics simulations. The results revealed that methyl groups are critical in mesophase layering in C-9 alkyl benzene- and naphthalene-based pitches, by reducing CH-π interaction. However, similar alignment could be achieved without the same degree of methyl substitutions for the coal tar-based pitch because of stronger π-π interaction than the other precursors. The insights from this study contribute to our understanding of the formation of conventional mesophase pitch and have implications for the processing of coal-derived materials. In conclusion, this knowledge is vital to produce valuable products like carbon fiber and graphite from pitches.

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

MELCOR-TMAP: The integration of MELCOR for fusion and TMAP4 for fusion reactor systems safety analysis and tritium inventory tracking

The design of fusion reactors requires the prediction of radioactive tritium migration during normal operating conditions and accident scenarios. Safety standards such as those from the U.S. Department of Energy (DOE) require that fusion reactor designs limit the release of hazardous material to the environment. In addition to safety requirements, the rarity of tritium supplies for burning D-T fusion plasma requires tritium breeding blankets. Computational modeling of migration of tritium and various reaction products throughout the fusion reactor serves to create designs that fulfill these requirements. MELCOR-TMAP, owned and developed at Idaho National Laboratory (INL), is an integrated system-level engineering code adapted for the analysis of fusion reactor systems. The history of the conception and evolution of the MELCOR-TMAP code, as well as its current capabilities, recent updates including a corrected composite material wall nodalization scheme, and plans for future development, are discussed here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

3D high-fidelity automated neutronics guided optimization of fusion blanket designs

The compact Fusion Pilot Plant (FPP) is defined in the recent National Academies of Sciences, Engineering, and Medicine report as the next step of fusion energy demonstration with a $50$ MWe peak net electricity production, $Q_e$ greater than $1$, and at least $3$ hours of continuous operation. This fusion pilot plant will be a test bed enabling materials, designs, and fuel management assessment, and it will represent an engineering challenge because of its high-fusion power and compact design targets. Previous reactor data is limited to experiments operating in different design space ranges. Therefore, design iterations and assessments should rely on high-fidelity first-principle theoretical and computational models. The high-fidelity integrated modeling of the plasma is a fundamental part of fusion energy research. However, the whole device modeling is often neglected, utilizing low-fidelity, system-level analysis. Recently, the need for high-fidelity multi-physics modeling was recognized, resulting in a selection of integrated tools. Further, autonomous design optimization requires a streamlined framework that perturbs the design point, reruns the analysis, and examines the outputs. However, high-fidelity analysis requires complex geometry specification that is difficult to perturb. This work presents the parametric CAD generation tool TRACER and a new neutronic workflow. TRACER allows the perturbation of the geometry representation, creating geometry files ready for further analysis. The streamlined neutronic workflow allows efficient and accurate calculations. The two new tools coupled together were used to perform a 3D high-fidelity multi-objective, multi-input optimization of an "ARC Class" compact tokamak design. The workflow was driven by an optimization driver for full automation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Investigation of the benefits of diabatic microreactors for process intensification of the water-gas shift reaction within the steam reforming process

Presently, the steam reforming process is the largest industrial source of hydrogen. Improving its efficiency can help to reduce associated carbon emissions and hydrogen production costs. Intensifying the water-gas shift reaction using microreactors with integrated cooling is one way of achieving this. In this study, a 2-D computational model of one of these microreactors is developed, validated with experimental data, and then used to demonstrate how microreactors can enhance the conversion of the water-gas shift reaction beyond what can be achieved using conventional packed bed reactors. These results are then generalized into a full system model of the steam reforming process to demonstrate how microreactors can reduce hydrogen production costs. The results suggest that microreactors can significantly reduce the required reactor volume and catalyst loading for the water-gas shift reaction and can similarly reduce the hydrogen production costs associated with the steam reforming process.

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

Recent advances in rational design of defect-engineered photocatalysts toward sustainable NH 3 synthesis as H 2 carrier: From fundamental and development to machine-learning

In this study, we provide a detailed overview of the fundamental mechanisms underpinning photocatalytic N 2 reduction. We also discuss advances in catalyst design for the synthesis of NH 3 . Particular emphasis is placed on the role of surface defect engineering, which includes the creation of surface defects to enhance the performance of semiconducting photocatalysts for efficient N 2 reduction. In addition, the application of a machine learning-based computational modeling approach is discussed as an important driving force for predicting and regulating NH 3 synthesis efficiency based on catalyst features and reaction conditions. Finally, existing challenges and future perspectives for improving the performance of defect-engineered photocatalysts are outlined to contribute to the ongoing discourse on sustainable ammonia generation. This review aims to clarify recent progress in the rational design of defect-containing photocatalysts for the synthesis of NH 3 and encourages innovative approaches to catalyst optimization rather than solely focusing on new materials.

08 HYDROGEN↗
DOE OSTI2022

Revealing melt flow instabilities in laser powder bed fusion additive manufacturing of aluminum alloy via in-situ high-speed X-ray imaging

Laser metal additive manufacturing technologies enable the fabrication of geometrically and compositionally complex parts unachievable by conventional manufacturing methods. However, the certification and qualification of additively manufactured parts are greatly hindered by the stochastic melt flow instabilities intrinsic to the process, which has not been explicitly revealed by direct observation. Here, we report the mechanisms of the melt flow instabilities in laser powder bed fusion additive manufacturing process revealed by in-situ high-speed high-resolution synchrotron X-ray imaging. Here we identified powder/droplet impact, significant keyhole oscillation, and melting-mode switching as three major mechanisms for causing melt flow instabilities. We demonstrated the detrimental consequences of these instabilities brought to the process, and presented new understanding on the melt flow evolution and keyhole oscillation. This work provides critical insights into process instabilities during laser metal additive manufacturing, which may guide the development of instability mitigation approaches. The results reported here are also important for the development and validation of high-fidelity computational models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Low-stress creep behavior of multiple salt caverns under cyclic operations

Here, the use of multiple salt caverns for large-scale underground hydrogen storage offers a strategic solution to the growing demand for efficient and flexible energy storage, due to the inherent low permeability, high mechanical strength, and self-healing properties of the salt formation. Despite these advantages, the long-term geomechanical response of salt caverns under cyclic injection and withdrawal remains a significant challenge for ensuring operational stability. This mechanistic study highlights the critical role of low-stress creep, also known as pressure solution creep, in governing deformation of salt caverns during cyclic operations. Unlike conventional creep mechanisms that dominate under high-stress conditions, low-stress creep becomes particularly relevant under the moderate stress changes induced by repeated injection and withdrawal, which provides a previously underexplored control on cavern closure. Our three-dimensional computational model integrates the low-stress creep mechanism into the Munson–Dawson creep equation across a spectrum of operational scenarios including variations in cavern pressure cycles, injection-withdrawal frequencies, and caprock mechanical property. The low-stress creep components accelerate the cavern creep closure, which is amplified with shorter cycle frequency and larger pressure difference. Sensitivity analyses capture the complex interactions between multiple caverns and surround formations at the field scale that allows for a quantitative assessment of deformation patterns and closure rates over time and thereby informs the design and operation of sustainable underground hydrogen storage system.

Cyclic operation↗
DOE OSTI2023

Electrochemical Li intercalation in b-As y P 1–y alloys: In-situ Raman spectroscopy study

Structural evolution of b-As y P 1-y alloys with varying As concentration, y was studied during Li-intercalation using in-situ Raman spectroscopy. A monotonic redshift of all the vibrational modes corresponding to P-P bonds (A 1 g , A 2 g , B 2g ), As-As bonds (A 1 g , A 2 g , B 2g ), and As-P bonds of all the samples was observed during the initial stages of the intercalation process due to the softening of each mode caused by the intercalation driven donor-type charge-transfer from Li to b-As y P 1-y . The emergence of a new peak identified as the E g mode of gray As, above an intercalation threshold, is indicative of the presence of an intercalation-driven structural phase segregation process. Further, the A 1g mode of gray As emerging after this intercalation threshold lies in close proximity of the A 2 g Raman mode of b-As y P 1-y . All the peaks beyond the intercalation threshold show an upshift due to the co-existence of gray As with b-As y P 1-y alloys causing strain and thus hardening of the phonon modes. In the sample with the highest As concentration phase segregation takes place during the synthesis process. Computational modeling reveals the co-existence of gray As in the b-As y P 1-y alloys with high As concentrations. It also confirms the existence of a local structural segregation taking place at a critical Li concentration during the intercalation process. Furthermore, this work shows the significance of intercalation on structural changes in the search for new phases of b-As y P 1-y alloys.

2D Materials↗
DOE OSTI2021

Covalent inhibition of hAChE by organophosphates causes homodimer dissociation through long-range allosteric effects

Acetylcholinesterase (EC 3.1.1.7), a key acetylcholine-hydrolyzing enzyme in cholinergic neurotransmission, is present in a variety of states in situ, including monomers, C-terminally disulfide-linked homodimers, homotetramers, and up to three tetramers covalently attached to structural subunits. Could oligomerization that ensures high local concentrations of catalytic sites necessary for efficient neurotransmission be affected by environmental factors? Using small-angle X-ray scattering (SAXS) and cryo-EM, we demonstrate that homodimerization of recombinant monomeric human acetylcholinesterase (hAChE) in solution occurs through a C-terminal four-helix bundle at micromolar concentrations. We show that diethylphosphorylation of the active serine in the catalytic gorge or isopropylmethylphosphonylation by the R P enantiomer of sarin promotes a 10-fold increase in homodimer dissociation. We also demonstrate the dissociation of organophosphate (OP)-conjugated dimers is reversed by structurally diverse oximes 2PAM, HI6, or RS194B, as demonstrated by SAXS of diethylphosphoryl-hAChE. However, binding of oximes to the native ligand-free hAChE, binding of high-affinity reversible ligands, or formation of an S P -sarin-hAChE conjugate had no effect on homodimerization. Dissociation monitored by time-resolved SAXS occurs in milliseconds, consistent with rates of hAChE covalent inhibition. OP-induced dissociation was not observed in the SAXS profiles of the double-mutant Y337A/F338A, where the active center gorge volume is larger than in wildtype hAChE. These observations suggest a key role of the tightly packed acyl pocket in allosterically triggered OP-induced dimer dissociation, with the potential for local reduction of acetylcholine-hydrolytic power in situ. Computational models predict allosteric correlated motions extending from the acyl pocket toward the four-helix bundle dimerization interface 25 Å away.

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

Optimal experimental design for prediction based on push-forward probability measures

Incorporating experimental data is essential for increasing the credibility of simulation-aided decision making and design. This paper presents a method which uses a computational model to guide the optimal acquisition of experimental data to produce data-informed predictions of quantities of interest (QoI). Many strategies for optimal experimental design (OED) select data that maximize some utility that measures the reduction in uncertainty of uncertain model parameters, for example the expected information gain between prior and posterior distributions of these parameters. In this paper, we seek to maximize the expected information gained from the pushforward of an initial (prior) density to the push-forward of the updated (posterior) density through the parameter-to-prediction map. The formulation presented is based upon the solution of a specific class of stochastic inverse problems which seeks a probability density that is consistent with the model and the data in the sense that the push-forward of this density through the parameter-to-observable map matches a target density on the observable data. While this stochastic inverse problem forms the mathematical basis for our approach, we develop a one-step algorithm, focused on push-forward probability measures, that leverages inference-for-prediction to bypass constructing the solution to the stochastic inverse problem. A number of numerical results are presented to demonstrate the utility of this optimal experimental design for prediction and facilitate comparison of our approach with traditional OED.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Analysis of the High-Pressure High-Temperature (HPHT) growth of single crystal diamond

A multi-scale, computational model is developed to describe the growth characteristics of single-crystal diamond in the High-Pressure, High-Temperature (HPHT) process. This model is the first to connect phase-change kinetics governing crystal growth to the continuum transport of carbon through the growth cell. Results show the importance of convective transport driven by buoyant flow in the metallic solvent, which increases the growth rate by nearly an order of magnitude over that obtained under diffusion alone. Parametric studies show how crystal growth may be kinetically-limited or transport-limited, depending on the value of the macroscopic kinetic coefficient. Estimating this kinetic coefficient from growth experiments yields a phase-change Damköhler number of unity, indicating a mixed regime where phase-change kinetics and transport are comparable and strongly coupled in this system. Mechanisms responsible for slowing growth as the crystal size increases are explained. Finally, supersaturation inhomogeneities along the facets of larger crystals are predicted, which may be relevant to solvent inclusion formation during growth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploring the impact of multivalent substitution on high-temperature electrical conductivity in doped lanthanum chromite perovskites: An experimental and ab-initio study

Doped LaCrO 3 perovskites hold promise as robust materials for electrical interconnects and sensor applications in harsh environments. In this work, we investigated the high-temperature behavior of La 1–x Sr x CrO 3 , La 1–x Ca x CrO 3 and La 0.8 Sr 0.2 Cr 1-x Mn x O3 (0.1 ≤ x ≤ 0.4) through a combination of computational modeling and physical characterization up to 1500 °C. Crystalline structural properties were determined and compared with ab-initio calculations, which demonstrated excellent agreement with experimental findings. High-temperature electrical conductivity measurements were performed under different atmospheres. Calcium and strontium/manganese co-doped lanthanum chromites exhibited typical semiconductor exponential trends and conductivity showed proportional correlation with substitutions levels up to 30 %. The DFT modelling was completed up to 1500 °C, including low and high temperature chromite phases and oxygen vacancies insertion. Calculations were correlated with experimental electrical properties. This work expands the understanding of doped lanthanum chromites and paves the way for the development of materials suitable for demanding high-temperature applications.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Structural and Biophysical Properties of a [4Fe-4S] Ferredoxin-Like Protein from Synechocystis sp. PCC 6803 with a Unique Two Domain Structure

Electron carrier proteins (ECPs), binding iron-sulfur clusters, are vital components within the intricate network of metabolic and photosynthetic reactions. They play a crucial role in the distribution of reducing equivalents. In Synechocystis sp. PCC 6803, the ECP network includes at least nine ferredoxins. Previous research, including global expression analyses and protein binding studies, has offered initial insights into the functional roles of individual ferredoxins within this network. This study primarily focuses on Ferredoxin 9 (slr2059). Through sequence analysis and computational modeling, Ferredoxin 9 emerges as a unique ECP with a distinctive two-domain architecture. It consists of a C-terminal iron-sulfur binding domain and an N-terminal domain with homology to Nil-domain proteins, connected by a structurally rigid 4-amino acid linker. Notably, in contrast to canonical [2Fe-2S] ferredoxins exemplified by PetF (ssl0020), which feature highly acidic surfaces facilitating electron transfer with photosystem I reaction centers, models of Ferredoxin 9 reveal a more neutral to basic protein surface. Using a combination of electron paramagnetic resonance spectroscopy and square-wave voltammetry on heterologously produced Ferredoxin 9, this study demonstrates that the protein coordinates 2x[4Fe-4S]2+/1+ redox-active and magnetically interacting clusters, with measured redox potentials of -420 +/- 9 mV and -516 +/- 10 mV vs SHE. A more in-depth analysis of Fdx9's unique structure and protein sequence suggests that this type of Nil-2[4Fe-4S] multi-domain ferredoxin is well conserved in cyanobacteria, bearing structural similarities to proteins involved in homocysteine synthesis in methanogens.

cyanobacteria↗
DOE OSTI2021

Spectroscopy and electronic structure of the low-energy states of ThN

The electronic spectrum of ThN over the range 19,600–21,200 cm -1 is remarkably congested, showing at least twenty vibronic bands originating from the ground state zero-point level. Rotationally resolved spectra for eleven of these bands have been examined to probe the underlying reasons for the dense manifold of states. Dispersed fluorescence spectra and fluorescence decay lifetimes were measured to provide addition insights. Our experimental measurements were complemented by electronic structure calculations for the low-energy doublet and quartet states. The ab initio calculations yielded a density of excited states that was consistent with the congested spectrum. These states were derived from the formal configurations Th 3+ (6d)N 3- , Th 2+ (7s 2 )N 2- (2p 5 ) and Th 2+ (7s6d)N 2 -(2p 5 ). The calculations indicated that the rotational and vibrational constants would be characteristic of the electronic configurations. This prediction was not borne out by the experimental data, suggesting that the full range of inter-configurational state interactions were not captured by the computational model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Revealing the synergistic effects of sequential and simultaneous dual beam irradiations in tungsten via in-situ TEM

Elucidating the synergistic effects of different energetic beams on the radiation response of nuclear materials is critical for developing an improved methodology for their evaluation when exposed to extreme environments. This article describes in-situ sequential (He implantation followed by Kr irradiation and vice-versa) and simultaneous (heavy ion Kr irradiation and He implantation) dual beam irradiations performed on tungsten at 1223 K. Dislocation loop density, average area, and total loop damage as a function of irradiation history and dose/fluence are quantified. The loop Burgers vectors and cavity damage (cavity density, size and total change in volume) are also determined at the final maximum dose for each condition. The loop damage evolution was different in all cases, with the smallest loop damage observed in the simultaneous experiment. Annihilation of Kr generated dislocation loops during He implantation in the Kr + He experiment was an unanticipated observation that may be explained by the dynamic evolution of dislocation loop sink strengths and time-dependent defect fluxes. Dislocation loop raft formation, denuded zones near extended defects, and cavity damage are compared across the different conditions. The phenomena observed and discussed in this work will stimulate further experimental and computational modeling activities leading to improved fundamental understanding of the irradiation response of nuclear materials under reactor-similar environments.

36 MATERIALS SCIENCE↗