Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “propagation database models”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Analyzing Hardware and Software Common Cause Failures in Digital Instrumentation and Control Systems using Dual Error Propagation Method

This paper develops a methodology for quantifying software common cause failures (CCFs) in digital instrumentation and control (I&C) systems of nuclear power plants. To support the transition of analog I&C systems to digital in nuclear power plants, probabilistic risk assessment (PRA) techniques are used. The hardware components of the I&C systems have reliability databases that can be used in the PRA studies. However, the failure data for redundant software components of the systems is sparse. Failure of components constitutes a CCF, wherein two or more components or systems fail due to a single shared cause and coupling mechanism. This paper proposes a quantification approach that can simultaneously model hardware and software components, incorporate the CCFs of software systems in the models, and bridge the gap between the failure quantification of models and the development of CCF parametric databases. We demonstrate the dual error propagation method (DEPM) by developing I&C systems failure models for a representative digital reactor trip system. The DEPM models are built to simulate the control and data flows within the systems and can accommodate failure states. By expanding DEPM to software CCFs, we generated alpha factor parameter estimates for each of the modeled error propagation mechanisms.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

ROADRUNNER uranium nitride MiniFuel: Experimental design, fabrication and pre-irradiation baseline characterization for accelerated burnup testing

Uranium nitride (UN) is a promising fuel candidate for advanced reactor systems owing to its high uranium density and thermal conductivity; however, its qualification remains constrained by the scarcity of well-controlled irradiation performance data. Here, to address this limitation, the ROADRUNNER (Research On ADvancing the peRformance of UraNium Nitrides in Extreme enviRonments) campaign employs the MiniFuel platform in the High Flux Isotope Reactor (HFIR) to enable accelerated burnup irradiation testing under tightly controlled and largely isothermal conditions. This paper presents the experimental design, fuel fabrication, and pre-irradiation baseline characterization of the ROADRUNNER UN MiniFuel campaign. Thirty-six UN minidisc specimens were fabricated with systematically varied as-fabricated density (86–96% of theoretical density), carbon impurity content (961–5240 ppm), oxygen content (≤ ∼2000 ppm), and grain size (2.5–24 μm). The irradiation matrix spans nominal fuel temperatures of 873 K, 1173 K, and 1473 K and target burnups of 3.75%, 6.0%, and 7.5% fissions per initial metal atom (FIMA). Neutronic and thermal analyses were performed to define specimen-specific burnup accumulation and temperature histories, establishing the boundary conditions for subsequent in-pile behavior. Comprehensive pre-irradiation characterization—including dimensional metrology, density verification, impurity analysis, X-ray diffraction, Raman spectroscopy, scanning electron microscopy, X-ray computed tomography, and confocal profilometry—provides a detailed baseline for post-irradiation examination. Pre-irradiation data were further used to generate predictive estimates of fission gas release and swelling using existing empirical correlations. This quantitative comparison reveals substantial inter-model divergence at intermediate and elevated temperatures that exceeds propagated input uncertainties, highlighting structural gaps in the historical irradiation database. The ROADRUNNER irradiation campaign is currently underway in HFIR, with initial firs cycle completed in late 2025 and remaining targets scheduled through 2027. The experimental design and baseline dataset presented here establish the framework needed to interpret forthcoming post-irradiation measurements and to provide discriminating data for the validation and refinement of physics-based UN fuel performance models.

Lopes, Denise Adorno [Oak Ridge National Laborator↗
DOE OSTI2015

Parallel database search and prime factorization with magnonic holographic memory devices

In this work, we describe the capabilities of Magnonic Holographic Memory (MHM) for parallel database search and prime factorization. MHM is a type of holographic device, which utilizes spin waves for data transfer and processing. Its operation is based on the correlation between the phases and the amplitudes of the input spin waves and the output inductive voltage. The input of MHM is provided by the phased array of spin wave generating elements allowing the producing of phase patterns of an arbitrary form. The latter makes it possible to code logic states into the phases of propagating waves and exploit wave superposition for parallel data processing. We present the results of numerical modeling illustrating parallel database search and prime factorization. The results of numerical simulations on the database search are in agreement with the available experimental data. The use of classical wave interference may results in a significant speedup over the conventional digital logic circuits in special task data processing (e.g., √n in database search). Potentially, magnonic holographic devices can be implemented as complementary logic units to digital processors. Physical limitations and technological constrains of the spin wave approach are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Impact of fission yield covariance matrices on decay heat uncertainty quantification with the DARWIN2 package

Although fission yields are strongly correlated, correlations between them are usually not taken into account when performing decay heat uncertainty calculations with the CEA DARWIN2 package due to the lack of reference covariance matrices associated to the JEFF-3.1.1 evaluation. However, covariance matrices for {sup 235}U and {sup 239}Pu thermal fission yields have been produced recently by the subgroup 37 of OECD/NEA Working Party on International Nuclear Data Evaluation Cooperation (WPEC). The study presented in this paper evaluates the effect of those covariances on the decay heat uncertainty calculation, both on fission burst experiments and on integral decay heat measurements that are part of the DARWIN2 experimental validation database. Although some differences are observed, partly due to the variety of models used to produce the matrices, the propagation of fission yield covariances always leads to a reduction of the decay heat uncertainty for cooling time above 10 seconds, indicating that not taking them into account is a conservative hypothesis from a safety point of view. Nevertheless, the strong impact of those covariances on the decay heat uncertainty also highlights the need of documented and consistent fission yield uncertainties and correlation data. On top of that, all the covariance matrices used for this study represent the correlation between the physical model parameters, but none of them includes the experimental correlations. An evaluation of the experimental correlations would also be of strong interest for decay heat uncertainty calculations. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Thermodynamic modeling with uncertainty quantification using the modified quasichemical model in quadruplet approximation: Implementation into PyCalphad and ESPEI

The modified quasichemical model in the quadruplet approximation (MQMQA) considers the first- and the second-nearest-neighbor coordination and interactions, particularly useful in describing short-range ordering (SRO) in complex liquids such as molten salts, slag in metal processing, and electrolytic solutions. Here, the present work implements the MQMQA into the Python based open-source software PyCalphad for thermodynamic calculations. This endeavor facilitates the development of MQMQA-based thermodynamic database with uncertainty quantification (UQ) and propagation (UP) using the open-source software ESPEI. A new database structure based on Extensible Markup Language (XML) is proposed for ESPEI evaluation of MQMQA model parameters. Using the KF-NiF 2 , KCl-NaCl-MgCl 2 , and CaCl 2 -CaF 2 -LiCl-LiF salt systems as examples, we demonstrate the successful implementation of MQMQA in PyCalphad through thermodynamic calculations of Gibbs energy, equilibrium quadruplet fractions, and phase diagram, as well as database development with UQ and UP using ESPEI. Furthermore, as an application of the present implementation, both the LiF–TbF 3 and LiF-HoF 3 systems have been modeled by MQMQA for the first time, which are in good agreement with experiments. The present implementation hence offers an open-source capability for performing CALPHAD modeling for complex liquids with SRO using MQMQA plus a new XML database structure.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Prefabricated High-Strength Rebar Systems with High-Performance Concrete for Accelerated Construction of Nuclear Concrete Structures

This project aimed to reduce the field construction times and fabrication costs of reinforced concrete (RC) nuclear structures through the use of: 1) high-strength steel deformed reinforcing bar (rebar); 2) prefabricated high-strength rebar (HSR) assemblies with headed anchorages; and 3) high-strength concrete (HSC). The focus of the project was not to develop new innovations on materials, but rather to fill a major knowledge gap on the effectiveness, code conformity, and viability (i.e., practicality, commercial availability) of existing high-strength materials in nuclear structures, especially in shear walls, their connections/joints, and around penetrations. These advances promote the Nuclear Energy Enabling Technologies Advanced Methods for Manufacturing (NEET-1) goals of accelerating deployment schedules by at least 6 months and reducing component fabrication costs by 20% or more. Toward this vision, the following objectives were identified and completed: Develop limit-benefit and cost-benefit analysis frameworks: The structural and economic limits and benefits of HSR and HSC in nuclear walls were evaluated. Importantly, it was found that modest increases in concrete strength up to 10.0 ksi and rebar yield strength up to 100 ksi result in the most economic designs and improved structural performance (Barbachyn et al., 2017a and Barbachyn et al., 2017b). Also, an industry survey was distributed to determine the cost benefits associated with prefabrication. Importantly, on-site worker-hours are significantly reduced, allowing for significantly accelerated construction schedules (Devine et al. 2018b). Develop an optimization methodology for structural design: An optimization study was deemed not necessary since the trivial solution with maximized use of high-strength materials and prefabrication always resulted in the largest cost and time savings. Instead, the benefits (and limits) of high-strength materials and prefabrication on construction costs and on-site construction time were investigated through a parametric study. Cost-reductions through the use of HSR and HSC were found to be greater than 40%. Also, on-site construction time of nuclear shear walls can be reduced by up to 83% through the use of prefabricated rebar assemblies. Conduct experimental evaluations of nuclear shear walls and wall-foundation joints: A large scale experimental program was completed on shear-critical RC cantilever deep beams. From these experiments, it was demonstrated that increasing rebar yield strength had a greater effect on increasing lateral strength than increasing the concrete compressive strength. Further, using HSR together with HSC resulted in the highest lateral strength and displacement ductility (Devine et al. 2018a). Additionally, large-scale reversed-cyclic lateral load tests were conducted on RC shear walls. It was demonstrated that through the utilization of HSR and HSC, a wall with 55% reduced vertical and horizontal rebar area was able to achieve 91% of the peak lateral strength of a conventional wall using state-of-practice materials with the same geometry. Further, trim reinforcement with headed anchorages around the penetrations proved to be effective. Develop validated numerical simulation models: Numerical modeling was completed in two phases during this project using 3 different finite element software packages. The first phase used existing measured data from previous experimental studies on stocky RC walls to validate simplified predictive models. Models using two different programs were validated based on a database of 38 stocky walls (Barbachyn et al., 2017a). The second phase occurred after experimental testing using a third and more detailed finite element package. Prediction improvements for the initial stiffness, lateral deformations, damage propagation, and failure mechanisms were targeted. The resulting detailed models more accurately predicted the peak lateral strength, initial and secant stiffnesses, lateral deformations, and cracking patterns of the wall test specimens subjected to reversed-cyclic lateral loading as compared to the simplified pretest models. Develop validated design procedures, tools, and criteria: This project evaluated the accuracy of existing closed-form design equations from U.S. building design standards and previous research (Barbachyn et al. 2017a and Devine et al. 2018a). It was shown that the current predictive equations for shear strength capacity need to be revisited as they provide significant scatter and gross overestimates, indicating the potential for unconservative designs. Further, for both the deep beam and shear wall specimens, the numerical models and flexural capacity predictions provided the most accurate strength predictions. Develop cost-effective field construction procedures: To promote the use of prefabricated rebar assemblies, the following construction conditions where prefabrication would be used instead of in-place rebar tying were identified from an industry survey: 1) to save on construction schedule; 2) to improve safety and/or quality control; 3) for areas with heavy rebar congestion; and 4) for structures with significant repetition in rebar layout. All of these conditions directly apply and benefit nuclear construction where construction costs and schedules are extensive, safety and quality are of utmost importance, the reinforcement is extremely dense, and the reinforcement layout is repetitive in many regions of the shear walls. Additionally, an experimental study was conducted to investigate the changes in bar spacing (with respect to code-required rebar spacing and placement tolerances) when a horizontally assembled prefabricated two-dimensional (2D) rebar mat or three-dimensional (3D) rebar cage is lifted (or tripped) to a vertical position. It was found that the bars involved in both the lateral translation of the rebar assemblies (i.e., movement of the assembly in its horizontal position) and the tripping process were most susceptible to large spacing changes, and thus, their positions within the assembly would need to be checked after the placement of the assembly in place within the formwork. Importantly, the spacings of all of the other bars in the assembly were not affected significantly by the lateral translation and tripping process, thereby eliminating the need to recheck their positions within the assembly (Devine et al. 2018b).

36 MATERIALS SCIENCE↗
DOE OSTI2015

Large Eddy Simulation Modeling of Flashback and Flame Stabilization in Hydrogen-Rich Gas Turbines Using a Hierarchical Validation Approach

This project was a combined computational and experimental effort to improve predictive capability for boundary layer flashback of premixed swirl flames relevant to gas-turbine power plants operating with high-hydrogen-content fuels. During the course of this project, significant progress in modeling was made on four major fronts: 1) use of direct numerical simulation of turbulent flames to understand the coupling between the flame and the turbulent boundary layer; 2) improved modeling capability for flame propagation in stratified pre-mixtures; 3) improved portability of computer codes using the OpenFOAM platform to facilitate transfer to industry and other researchers; and 4) application of LES to flashback in swirl combustors, and a detailed assessment of its capabilities and limitations for predictive purposes. A major component of the project was an experimental program that focused on developing a rich experimental database of boundary layer flashback in swirl flames. Both methane and high-hydrogen fuels, including effects of elevated pressure (1 to 5 atm), were explored. For this project, a new model swirl combustor was developed. Kilohertz-rate stereoscopic PIV and chemiluminescence imaging were used to investigate the flame propagation dynamics. In addition to the planar measurements, a technique capable of detecting the instantaneous, time-resolved 3D flame front topography was developed and applied successfully to investigate the flow-flame interaction. The UT measurements and legacy data were used in a hierarchical validation approach where flows with increasingly complex physics were used for validation. First component models were validated with DNS and literature data in simplified configurations, and this was followed by validation with the UT 1-atm flashback cases, and then the UT high-pressure flashback cases. The new models and portable code represent a major improvement over what was available before this project was initiated.

42 ENGINEERING↗
DOE OSTI2025

Machine learning enhanced predictions of ICRF heating: Overcoming numerical limitations via data curation

In this work, we present the development of robust surrogate models for Ion Cyclotron Range of Frequencies (ICRF) and High-Harmonic Fast Wave (HHFW) heating predictions in fusion plasmas. Building upon our previous efforts to achieve real-time capable models, we identify the cause of the outliers found using TORIC in certain HHFW heating scenarios. The outliers are observed to be spurious ion Bernstein wave (IBW)-like modes caused by a wavelength control algorithm designed to address challenging scenarios with high perpendicular wavenumbers. The effect arises from the modulation in the perpendicular susceptibility, which can induce sign reversal and IBW-like propagation for scenarios featuring normalized ion Larmor radius λ i ≫ 1. We use TORIC with this algorithm disabled to generate a novel HHFW-NSTX database that is free of outliers. Surrogate models trained on this database, including Random Forest Regressor (RFR), Multi-Layer Perceptrons, and Gaussian Process Regressors (GPR), demonstrate the ability to accurately predict HHFW heating profiles, with regression scores of R 2 ∈[0.93−0.99]. Additionally we demonstrate that it is possible to generalize predictions beyond training data by the use of both RFR and GPR models, enabling the prediction of scenarios previously limited to the original model. GPR models also provide uncertainty quantification, offering insights into model confidence. This work introduces a comprehensive Verification, Validation, and Uncertainty Quantification methodology for surrogate modeling, applicable not only to ICRF heating but also to other RF heating challenges and fusion physics problems. Beyond accelerated inference, these models show effective extrapolation capabilities, providing an alternative for addressing numerical challenges.

Artificial neural networks↗
DOE OSTI2019

An artificial neural network framework for reduced order modeling of transient flows

Here, this paper proposes a supervised machine learning framework for the non-intrusive model order reduction of unsteady fluid flows to provide accurate predictions of non-stationary state variables when the control parameter values vary. Our approach utilizes a training process from full-order scale direct numerical simulation data projected on proper orthogonal decomposition (POD) modes to achieve an artificial neural network (ANN) model with reduced memory requirements. This data-driven ANN framework allows for a nonlinear time evolution of the modal coefficients without performing a Galerkin projection. Our POD-ANN framework can thus be considered an equation-free approach for latent space dynamics evolution of nonlinear transient systems and can be applied to a wide range of physical and engineering applications. Within this framework we introduce two architectures, namely sequential network (SN) and residual network (RN), to train the trajectory of modal coefficients. We perform a systematic analysis of the performance of the proposed reduced order modeling approaches on prediction of a nonlinear wave-propagation problem governed by the viscous Burgers equation, a simplified prototype setting for transient flows. We find that the POD-ANN-RN yields stable and accurate results for test problems assessed both within inside and outside of the database range and performs significantly better than the standard intrusive Galerkin projection model. Finally, our results show that the proposed framework provides a non-intrusive alternative to the evolution of transient physics in a POD basis spanned space, and can be used as a robust predictive model order reduction tool for nonlinear dynamical systems.

42 ENGINEERING↗
DOE OSTI2023

Evaluation of a CFD Commercial Package to the Modeling of Acoustic Wave Propagation in Bubbly-Liquid Column

This document reports on the considerations and development of an initial computational model based on the first principles for predicting acoustic propagation in a wellbore based on a logging-while-drilling (LWD) tool. This tool could be used to create a large database of simulated signals mimicking gas-influx-while-drilling events (so-called “kicks”). While other logging techniques, such as electrical resistivity and bulk density, are available in LWD tools, the acoustic logging method was identified in a complementary experimental endeavor as the most appropriate for early kick detection due to its high sensitivity to the presence of gas in liquids.

03 NATURAL GAS↗
DOE OSTI2017

Alfvén Wave Turbulence as a Coronal Heating Mechanism: Simultaneously Predicting the Heating Rate and the Wave-induced Emission Line Broadening

We test the predictions of the Alfvén Wave Solar Model (AWSoM), a global wave-driven magnetohydrodynamic (MHD) model of the solar atmosphere, against high-resolution spectra emitted by the quiescent off-disk solar corona. AWSoM incorporates Alfvén wave propagation and dissipation in both closed and open magnetic field lines; turbulent dissipation is the only heating mechanism. We examine whether this mechanism is consistent with observations of coronal EUV emission by combining model results with the CHIANTI atomic database to create synthetic line-of-sight spectra, where spectral line widths depend on thermal and wave-related ion motions. This is the first time wave-induced line broadening is calculated from a global model with a realistic magnetic field. We used high-resolution SUMER observations above the solar west limb between 1.04 and 1.34 R {sub ⊙} at the equator, taken in 1996 November. We obtained an AWSoM steady-state solution for the corresponding period using a synoptic magnetogram. The 3D solution revealed a pseudo-streamer structure transversing the SUMER line of sight, which contributes significantly to the emission; the modeled electron temperature and density in the pseudo-streamer are consistent with those observed. The synthetic line widths and the total line fluxes are consistent with the observations for five different ions. Further, line widths that include the contribution from the wave-induced ion motions improve the correspondence with observed spectra for all ions. We conclude that the turbulent dissipation assumed in the AWSoM model is a viable candidate for explaining coronal heating, as it is consistent with several independent measured quantities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

How does uncertainty of soil organic carbon stock affect the calculation of carbon budgets and soil carbon credits for croplands in the U.S. Midwest?

Cropland carbon budget depicts the amount of carbon flowing in and out of agroecosystems and the changes in carbon stocks of soil and living biomass during the same period. Soil carbon credit is the additional change in soil carbon stock under certain farming practices compared with the business-as-usual practices. Accurately calculating cropland carbon budget and soil carbon credit is critical to assessing climate change mitigation potential in agroecosystems. The calculation of cropland carbon budget and soil carbon credit is sensitive to local soil and climatic conditions, especially initial soil organic carbon (SOC) stock, which is determined by both SOC concentration (SOC%) and bulk density (Bulk_Density). SOC stock data are either from soil sampling or gridded public survey data. In agroecosystem models, SOC stock data are a key model input for quantifying cropland carbon budget and soil carbon credit. However, various types and degrees of uncertainties exist in SOC stock datasets, which propagate to the quantification of SOC stock change. In particular, a large discrepancy is found in two widely used SOC stock datasets — Rapid Carbon Assessment dataset (RaCA) and Gridded Soil Survey Geographic Database (gSSURGO) — in the U.S. Midwest, with a relative difference (quantified using Normalized Root Mean Square Error, NRMSE) of 48.0% for 0–30 cm SOC stock between the two datasets. It remains largely unclear how uncertainty in SOC stocks affects the calculation of cropland carbon budget and soil carbon credit. To address this question, we used a well-validated process-based agroecosystem model, ecosys, to assess the impacts of SOC stock uncertainty on carbon budget and soil carbon credit calculation in the U.S. Midwestern corn-soybean rotation systems. Our results reveal the following findings: (1) A sizable discrepancy exists in simulated cropland carbon budget between using gSSURGO and using RaCA for their SOC% and Bulk_Density as model inputs, with a Pearson correlation coefficient (r) of only 0.4 for simulated change of SOC stock (ΔSOC) using these two different soil datasets. (2) Simulated cropland carbon budget components were more sensitive to initial SOC% than to Bulk_Density. For example, the upper and lower quartiles of multi-year averaged ΔSOC were –29.8 and 4.8 gC/m 2 /year for the selected counties respectively, with an uncertainty of 13.7 and 0.7 gC/m 2 /year induced by uncertainties in initial SOC% and Bulk_Density, respectively. (3) Both simulated ΔSOC and its uncertainty were negatively correlated with initial SOC%, whereas ΔSOC was negatively correlated with air temperature, and ΔSOC uncertainty was positively correlated with air temperature. (4) The uncertainty of calculated soil carbon credits was much smaller compared with the uncertainty of calculated absolute carbon budgets assuming the same SOC stock uncertainty level in the inputs. Specifically, in our assessment comparing planting cover crops vs no cover crop, the uncertainty of calculated soil carbon credits induced by initial SOC% uncertainty was less than 4% (relative to the quantified value of the soil carbon credits) for 90% of the cases. Our analysis highlights that high accuracy measurement of SOC% as inputs is needed for the calculation of cropland carbon budgets; however, soil carbon credit quantification is much less sensitive to the initial SOC% inputs, and the current publicly available soil datasets (e.g., gSSURGO) are largely suitable for the calculation of soil

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

InfraPy: Python-Based Signal Analysis Tools for Infrasound

InfraPy is a Python-based analysis toolkit being development at LANL. The algorithms are intended for ground-based nuclear detonation detection applications to detect, locate, and characterize explosive sources using infrasonic observations. The implementation is usable as a stand-alone Python library or as a command line driven tool operating directly on a database. With multiple scientists working on the project, we've begun using a LANL git repository for collaborative development and version control. Current and planned work on InfraPy focuses on the development of new algorithms and propagation models. Collaboration with Southern Methodist University (SMU) has helped identify bugs and limitations of the algorithms. The current focus of usage development is focused on library imports and CLI.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

MarpoDB: An open registry for Marchantia polymorpha genetic parts

Marchantia polymorpha is an extant relative of the earliest terrestrial plants and has attracted a substantial interest as a model organism for evolutionary and developmental studies. Given its relatively simple genome, compact gene families, simple morphology, ease of propagation and transformation, M. polymorpha is becoming a promising platform for plant synthetic biology. Modular genetic parts have been essential for development of synthetic biology approaches, so we sought to design an engineering oriented database for M. polymorpha genetic parts where each gene is a standalone functional unit. MarpoDB is a database of M. polymorpha genes and genetic parts, which is tailored to become an integral tool for a synthetic biology workflow. Among its features are precompiled cross-database querying to InterPro, Pfam signatures and non-redundant Viridiplantae BLAST annotations; BLAST querying to M. polymorpha genes; sequence export in GenBank format; recoding of sequences to the common syntax for type IIS assembly and exchange of DNA parts; and a minimalistic, intuitive and interactive user interface for gene models and sequence exploration. Furthermore, we have implemented user input to encourage feedback, collaboration and exchange between the MarpoDB community. MarpoDB source-code is released on GitHub to promote development of computational tools for synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Predictive LES Modeling and Validation of High-Pressure Turbulent Flames and Flashback in Hydrogen-enriched Gas Turbines

The overall objective of this project is to make advances in large-eddy simulation modeling of swirl-flame flashback at elevated pressure, where the reactants can be premixed or stratified, and composed of methane enriched with hydrogen. The modeling effort was conducted with a companion experimental validation effort that focused on high-speed measurements of flashback in a model swirl-flame combustor. The modeling effort focused on the development of manifold-based combustion model for flashback simulations. The main complexity in this configuration comes from the flames’ propagation through a stratified fuel-air mixture along with heat loss near the walls. For this purpose, a non-adiabatic flamelet model was developed and tested using canonical flow configurations. Further, extensive validation in premixed configurations is used to establish the model’s accuracy in capturing heat loss and flame propagation. Finally, the model was tested using the UT flame flashback experiment, and was shown to accurately capture the flame flashback process. The end products from this work are a) a full validated flameout-based approach for modeling flame flashback; b) a complete implementation of the model in the open source OpenFOAM framework, which has been made available to both industry and academia; c) an extensive suite of validation studies performed using the OpenFOAM solver. In addition, new techniques for flame modeling, subfilter modeling of scalars, and simulation of low Mach number flows were developed and tested. In the experimental validation effort, flashback was studied in a premixed and stratified swirl flame by using high-speed imaging diagnostics. Fuels tested included CH 4 enriched with H 2 , with H 2 enrichment up to 84% by volume. These flashback experiments are conducted at pressures ranging from 1 to 5 atm. A swirler-based fuel-injection system was used to create the fuel-air stratification in the radial direction. For the high-pressure measurements, an optically accessible elevated pressure chamber was used. The spatial distribution of the equivalence ratio under non-reacting conditions was characterized using planar laser-induced fluorescence with acetone as the fuel tracer. It was observed that fuel-air pockets were distributed across the mixing tube width, although in an average sense, the fuel-air mixture was radially stratified. The global behavior of upstream flame propagation was investigated for different levels of hydrogen-enrichment. For stratified hydrogen-rich flashback, the propagation path of the flame changed from the inner wall to outer wall, which was induced by the faster chemistry of stoichiometric mixtures that are frequently present near the outer wall. This behavior of hydrogen-rich flashback persists even at elevated pressures up to 5 atm, although the propagation of the flame occurs as a wide flame tongue as opposed to the acute-tipped flame structures present in the atmospheric cases. The end product of the experimental effort is greatly improved knowledge of flashback physics and a unique database that can be used for model validation.

08 HYDROGEN↗
DOE OSTI2022

Uncertainty propagation in pore water chemical composition calculation using surrogate models

Performance assessment in deep geological nuclear waste repository systems necessitates an extended knowledge of the pore water chemical conditions prevailing in host-rock formations. In the last two decades, important progress has been made in the experimental characterization and thermodynamic modeling of pore water speciation, but the influence of experimental artifacts and uncertainties of thermodynamic input parameters are seldom evaluated. In this respect, we conducted an uncertainty propagation study in a reference geochemical model describing the pore water chemistry of the Callovian-Oxfordian clay formation. Nineteen model input parameters were perturbed, including those associated to experimental characterization (leached anions, exchanged cations, cation exchange selectivity coefficients) and those associated to generic thermodynamic databases (solubilities). A set of 13 quantities of interest were studied by the use of polynomial chaos expansions built non-intrusively with a least-squares forward stepwise regression approach. Training and validation sets of simulations were carried out using the geochemical speciation code PHREEQC. The statistical results explored the marginal distribution of each quantity of interest, their bivariate correlations as well as their global sensitivity indices. The influence of the assumed distributions for input parameters uncertainties was evaluated by considering two parametric domain sizes.

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

Global adjoint tomography—model GLAD-M25

SUMMARY Building on global adjoint tomography model GLAD-M15, we present transversely isotropic global model GLAD-M25, which is the result of 10 quasi-Newton tomographic iterations with an earthquake database consisting of 1480 events in the magnitude range 5.5 ≤ Mw ≤ 7.2, an almost sixfold increase over the first-generation model. We calculated fully 3-D synthetic seismograms with a shortest period of 17 s based on a GPU-accelerated spectral-element wave propagation solver which accommodates effects due to 3-D anelastic crust and mantle structure, topography and bathymetry, the ocean load, ellipticity, rotation and self-gravitation. We used an adjoint-state method to calculate Fréchet derivatives in 3-D anelastic Earth models facilitated by a parsimonious storage algorithm. The simulations were performed on the Cray XK7 ‘Titan’ and the IBM Power 9 ‘Summit’ at the Oak Ridge Leadership Computing Facility. We quantitatively evaluated GLAD-M25 by assessing misfit reductions and traveltime anomaly histograms in 12 measurement categories. We performed similar assessments for a held-out data set consisting of 360 earthquakes, with results comparable to the actual inversion. We highlight the new model for a variety of plumes and subduction zones.

58 GEOSCIENCES↗