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 1,747 records · Page 97

A composite position independent monitor of reactor fuel irradiation using Pu, Cs, and Ba isotope ratios

When post-irradiation materials from the nuclear fuel cycle are released to the environment, certain isotopes of actinides and fission products carry signatures of irradiation history that can potentially aid a nuclear forensic investigation into the material's provenance. In this study, combinations of Pu, Cs, and Ba isotope ratios that produce position (in the reactor core) independent monitors of irradiation history in spent light water reactor fuel are identified and explored. These position independent monitors (PIMs) are modeled for various irradiation scenarios using automated depletion codes as well as ordinary differential equation solutions to approximate nuclear physics models. Experimental validation was performed using irradiated low enriched uranium oxide fuel from a light water reactor, which was sampled at 8 axial positions from a single rod. Plutonium, barium and cesium were chemically separated and isotope ratio measurements of the separated solutions were made by quadrupole and multi-collector inductively coupled mass spectrometry (Cs and Pu, respectively) and thermal ionization mass spectrometry (Ba). The effect of axial variations in neutron fluence and energy spectrum are evident in the measured isotope ratios. Furthermore, two versions of a combined Pu and Cs based PIM are developed. A linear PIM model, which can be used to solve for irradiation time is found to work well for natural U fuel with <10% 240 Pu and known or short cooling times. A non-linear PIM model, which cannot be solved explicitly for irradiation time without additional information, can nonetheless still group samples by irradiation history, including high burnup LEU fuel with unknown cooling time. 137 Ba/ 138 Ba is also observed to act as a position independent monitor; it is nearly single valued across the sampled fuel rod, indicating that samples sharing an irradiation history (same irradiation time and cooling time) in a reactor despite experiencing different neutron fluxes will have a common 137 Ba/ 138 Ba ratio. Modeling of this Ba PIM shows it increases monotonically with irradiation and cooling time, and a confirmatory first order analytical solution is also presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Carbon-enhanced metal-poor star frequencies in the galaxy: corrections for the effect of evolutionary status on carbon abundances

We revisit the observed frequencies of carbon-enhanced metal-poor (CEMP) stars as a function of the metallicity in the Galaxy, using data from the literature with available high-resolution spectroscopy. Our analysis excludes stars exhibiting clear overabundances of neutron-capture elements and takes into account the expected depletion of surface carbon abundance that occurs due to CN processing on the upper red giant branch. This allows for the recovery of the initial carbon abundance of these stars, and thus for an accurate assessment of the frequencies of carbon-enhanced stars. The correction procedure we develop is based on stellar-evolution models and depends on the surface gravity, log g, of a given star. Our analysis indicates that for stars with [Fe/H] ≤–2.0, 20% exhibit [C/Fe] ≥+0.7. This fraction increases to 43% for [Fe/H] ≤–3.0 and 81% for [Fe/H] ≤–4.0, which is higher than have been previously inferred without taking the carbon abundance correction into account. These CEMP star frequencies provide important inputs for Galactic and stellar chemical evolution models, as they constrain the evolution of carbon at early times and the possible formation channels for the CEMP-no stars. We also have developed a public online tool with which carbon corrections using our procedure can be easily obtained.

79 ASTRONOMY AND ASTROPHYSICS↗

Study of Data-Driven Mesh-Model Optimization in System Thermal-Hydraulic Simulation

Quantification of nuclear power plant safety risk requires a systematic and yet practical approach to identification of accident scenarios, assessment of their likelihood and consequences. Instrumental to this goal is risk-informed safety margin characterization (RISMC) framework, whose realization requires computationally robust and affordable methods for sufficiently accurate simulation of complex multi-dimensional physical phenomena, such as turbulent and multi-phase flow. The CFD-like correlation-based codes with 3D simulation capability (e.g., GOTHIC) ensure computational efficiency using coarse mesh size and the sub-grid phenomena in the boundary layer that can be captured by adequate constitutive correlations. However, the error sources and user effects on the selection of mesh size and models lead to unpredictable simulation error, while rich High-Fidelity (HF) data from experiments and numerical simulation using validated code or DNS are not fully explored. It would be useful to have a “smart” data-driven multi-scale framework in which the low-resolution models can be “taught” to emulate high-resolution models. The objective of this work is to develop and evaluate a physical-based data-driven mesh-model optimization approach (Optimized Mesh/Model Information System, OMIS) to estimate the simulation error and give advice on the optimized selection of coarse mesh size and Low-Fidelity (LF) models for System Thermal Hydraulic (STH) simulation to achieve accuracy comparable to that of HF models. This approach takes advantages of computational efficiency of coarse-mesh simulation and application of Machine Learning (ML) algorithms.

97 MATHEMATICS AND COMPUTING↗

WE-AB-BRA-02: Development of Biomechanical Models to Describe Dose-Volume Response to Liver Stereotactic Body Radiation Therapy (SBRT) Patients

Purpose: Previous studies have shown that radiotherapy treatment for liver metastases causes marked liver hypertrophy in areas receiving low dose and atrophy/fibrosis in areas receiving high dose. The purpose of this work is to develop and evaluate a biomechanical model-based dose-response model to describe these liver responses to SBRT. Methods: In this retrospective study, a biomechanical model-based deformable registration algorithm, Morfeus, was expanded to include dose-based boundary conditions. Liver and tumor volumes were contoured on the planning images and CT/MR images three months post-RT and converted to finite element models. A thermal expansion-based relationship correlating the delivered dose and volume response was generated from 22 patients previously treated. This coefficient, combined with the planned dose, was applied as an additional boundary condition to describe the volumetric response of the liver of an additional cohort of metastatic liver patients treated with SBRT. The accuracy of the model was evaluated based on overall volumetric liver comparisons and the target registration error (TRE) using the average deviations in positions of identified vascular bifurcations on each set of registered images, with a target accuracy of the 2.5mm isotropic dose grid (vector dimension 4.3mm). Results: The thermal expansion coefficient models the volumetric change of the liver to within 3%. The accuracy of Morfeus with dose-expansion boundary conditions a TRE of 5.7±2.8mm compared to 11.2±3.7mm using rigid registration and 8.9±0.28mm using Morfeus with only spatial boundary conditions. Conclusion: A biomechanical model has been developed to describe the volumetric and spatial response of the liver to SBRT. This work will enable the improvement of correlating functional imaging with delivered dose, the mapping of the delivered dose from one treatment onto the planning images for a subsequent treatment, and will further provide information to assist with the biological characterization of patients’ response to radiation.

60 APPLIED LIFE SCIENCES↗

Nuclear Technology Research and Development (NTRD) Technical Monthly November FY18

Single crystals of U 3 Si 2 were grown at various crystal rotation and pull rates to establish the best parameters for growth of crystals with diameters ≥5 mm (Figure 1). Crystals were pulled from the melt using a modified Czochralski method. U 3 Si 2 single crystals will enable more fundamental property studies for this proposed high density fuel providing data to develop and inform physics-based fuel performance models. (D. Byler)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Semi–Analytical Modeling of Transient Stream Drawdown and Depletion in Response to Aquifer Pumping

Analytical and semi–analytical models for stream depletion with transient stream stage drawdown induced by groundwater pumping are developed to address a deficiency in existing models, namely, the use of a fixed stream stage condition at the stream–aquifer interface. Here field data are presented to demonstrate that stream stage drawdown does indeed occur in response to groundwater pumping near aquifer–connected streams. A model that predicts stream depletion with transient stream drawdown is developed based on stream channel mass conservation and finite stream channel storage. The resulting models are shown to reduce to existing fixed–stage models in the limit as stream channel storage becomes infinitely large, and to the confined aquifer flow with a no–flow boundary at the streambed in the limit as stream storage becomes vanishingly small. The model is applied to field measurements of aquifer and stream drawdown, giving estimates of aquifer hydraulic parameters, streambed conductance, and a measure of stream channel storage. The results of the modeling and data analysis presented herein have implications for sustainable groundwater management.

54 ENVIRONMENTAL SCIENCES↗

Percheron Power Archimedes Screw Turbine Analysis (Abstract)

The development of the Percheron Power (Percheron) composite Archimedes screw turbine (AST) is based upon employing newly developed theoretical mathematical models for optimizing the water volume carried through the AST for a given diameter, combined with minimizing efficiency losses due to bypass flow and friction. The key feasibility issues associated with producing the composite blades will be to understand the strength, fatigue, and wear performance of the AST blades and associated components, and how these factors change along the length of the turbine. To get to this understanding of material property requirements, Percheron needs free-surface computational fluid dynamics (CFD) models of the AST designs and to perform interactive modeling between the turbine blade optimization and component design based on finite element method stress analysis (FEA). The workscope and analyses requested by Percheron in the Small Business Voucher (SBV) program are well aligned with PNNL capability and recent work experience in performing CFD and FEA analysis for hydropower and marine energy turbines. The CFD analysis will provide estimates of power, efficiency, and surface forces on the blades. The forces on the blade will be transferred to the FEA analysis that will be used to evaluate the material stresses. The CFD and FEA analyses will be repeated for different AST designs provided to PNNL by Percheron. The work conducted in this SBV project will complement the existing Percheron FOA project by providing better defined parameters for design and material selection for the turbine unit prior to manufacture of initial test units.

42 ENGINEERING↗

A Study of Propane Combustion in a Spark-Ignited Cooperative Fuel Research (CFR) Engine

Liquefied petroleum gas (LPG), whose primary composition is propane, is a promising candidate for heavy-duty vehicle applications as a diesel fuel alternative due to its CO 2 reduction potential and high knock resistance. To realize diesel-like efficiencies, spark-ignited LPG engines are proposed to operate near knock-limit over a wide range of operating conditions, which necessitates an investigation of fuel-engine interactions that leads to end-gas autoignition with propane combustion. This work presents both experimental and numerical studies of stoichiometric propane combustion in a sparkignited (SI) cooperative fuel research (CFR) engine. Engine experiments are initially conducted at different compression ratio (CR) values, and the effects of CR on engine combustion are characterized. A three-pressure analysis (TPA) model based on the two-zone combustion concept is developed in GT-Power and validated using test results to estimate in-cylinder wall temperatures, residual gas fraction, etc. This model is further utilized to examine end-gas chemistry by enabling the SI turbulent flame combustion and unburned gas chemical kinetics modules. Finally, a three-dimensional (3D) computational fluid dynamic (CFD) model of the CFR engine is developed in CONVERGE, where the G-equation and SAGE detailed chemical kinetics models are implemented for combustion modeling. Here, a 153 species reduced chemical kinetics mechanism derived from the detailed NUIGMech1.1 mechanism based on the ignition delay and laminar flame speed (LFS) studies is used to generate an LFS lookup table and to describe end-gas autoignition chemistry. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations are then performed for the tested CRs, and the numerical model is shown to be capable of predicting the propane combustion characteristics, particularly the end-gas autoignition behavior.

20 FOSSIL-FUELED POWER PLANTS↗

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Experimental characterization and constitutive modeling of the mechanical behavior of molybdenum under electromagnetically applied compression-shear ramp loading

Magnetically applied pressure-shear (MAPS) is a new experimental technique that provides a platform for direct measurement of material strength at extreme pressures. The technique employs an imposed quasi-static magnetic field and a pulsed power generator that produces an intense current on a planar driver panel, which in turn generates high amplitude magnetically induced longitudinal compression and transverse shear waves into a planar sample mounted on the drive panel. In order to apply sufficiently high shear traction to the test sample, a high strength material must be used for the drive panel. Molybdenum is a potential driver material for the MAPS experiment because of its high yield strength and sufficient electrical conductivity. To properly interpret the results and gain useful information from the experiments, it is critical to have a good understanding and a predictive capability of the mechanical response of the driver. In this work, the inelastic behavior of molybdenum under uniaxial compression and biaxial compression-shear ramp loading conditions is experimentally characterized. It is observed that an imposed uniaxial magnetic field ramped to approximately 10 T through a period of approximately 2500 μs and held near the peak for about 250 μs before being tested appears to anneal the molybdenum panel. In order to provide a physical basis for model development, a general theoretical framework that incorporates electromagnetic loading and the coupling between the imposed field and the inelasticity of molybdenum was developed. Based on this framework, a multi-axial continuum model for molybdenum under electromagnetic loading is presented. The model reasonably captures all of the material characteristics displayed by the experimental data obtained from various experimental configurations. Additionally, data generated from shear loading provide invaluable information not only for validating but also for guiding the development of the material model for multiaxial loadings.

36 MATERIALS SCIENCE↗

Advanced High-Fidelity Lumped EMT Grid Modelling & Comparison

There is increased interest in performing electromagnetic transient (EMT) studies of large power grids for stability purposes. The large power grids may not represent the entire grid, but portions of interest. For example, studies are being performed to identify the impact of extreme events on solar power plants in California that requires the EMT model of the grid in California, but not the complete western interconnection (WI) grid model to which the grid in California is connected. The EMT models of grids are necessary to perform integration studies of large power electronics systems like utility-scale photovoltaic (PV) systems, high-voltage direct current (HVdc) systems, flexible alternating current transmission systems (FACTs), and others. Traditionally, models of grids available to utilities are the transient stability (TS) models. Hence, it assumes significance to identify methods of developing high-fidelity EMT models for portions of the grid (of interest) from the information available in the existing TS models of the entire grid. In this paper, a conversion mechanism along with upgradation based on data-driven methods is applied to develop EMT models of grids. The conversion mechanism utilizes information from TS models of grids to generate EMT models of grids. These models can be extracted for the required portion of the grid. While these models adequately represent the local voltage behaviour of the grid, they are inadequate in representing the frequency behaviour of the grid. The data-driven methods are applied to identify frequency dynamics of the grid, and the EMT models are upgraded with the data-driven methods. The resulting upgraded EMT models of grids can perform the studies mentioned above. These models are tested under contingencies of interest to validate the developed models.

Debnath, Suman↗

Modeling Long-term Creep Performance for Welded Nickel-base Superalloy Structures for Power Generation Systems

The goal of this project is to model long-term creep performance for nickel-base superalloy weldments in high temperature power generation systems. The project uses physics-based modeling methodologies and algorithms for predicting alloy properties in heterogeneous material structures. The modeling methodology will be demonstrated on a gas turbine combustor liner weldment of Haynes 282 precipitate-strengthened nickel-base superalloy. The major developments are: (1) microstructure-property relationships under creep conditions and microstructure characterization (2) modeling inhomogeneous microstructure in superalloy weld (3) modeling mesoscale plastic deformation in superalloy weld and (4) a constitutive creep model that accounts for weld and base metal microstructure and their long term evolution. The developed modeling technology is aimed to provide a more efficient and accurate assessment of a material’s long-term performance compared with current testing and extrapolation methods. This modeling technology will also accelerate development and qualification of new materials in advanced power generation systems. This document is a final technical report for the project, covering efforts conducted from October 2014 to December 2016.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Generic and Multifunctional Electromagnetic Transient Model for Grid-Following Inverters

This article presents a generic and multi-functional electromagnetic transient (EMT) dynamic model of grid following (GFL) inverter-based resource (IBR) using the PSCAD software platform. The features of the developed model includes the flexibility in selecting various types and combinations of DC sources covering PV modules, battery modules, ideal DC source module, as well as flexibility in selecting either switched or averaged model of inverter. This model also covers exhaustive lists of controller logic covering open-loop/closed-loop PQ dispatch control, DC voltage and AC terminal voltage control along with the conventional current control designed in dq-domain, ate- domain and positive-negative sequence domain. Moreover, this model is equipped with the flexibility in selecting various types of current limiting schemes that includes saturation-based as well as latching-based current limiter, anti-windup protection. Moreover, the EMT model is agnostic to the MVA rating and is suitable for interfacing transmission systems by being complaint with the IEEE Std. 2800. The generality in the power circuits and the multi-functional options in operation and control of the developed EMT model makes it suitable for both academia and industry to study various power system aspects not limited to but such as fault behavior of GFL IBR and impacts on protection system, transient stability of a system interfaced with large number of GFL IBRs etc.

14 SOLAR ENERGY↗

Modeling the effect of habitat availability and quality on endangered winter-run Chinook salmon ( Oncorhynchus tshawytscha ) production in the Sacramento Valley

Reconciliation between water uses and sustaining aquatic species populations requires an effort to identify and quantify essential habitat characteristics for ecosystem health and valuation methods to predict ecosystem response to restoration actions. This process is particularly challenging for anadromous fish species such as California's Sacramento River winter-run Chinook salmon, due to their limited geographic range and diverse life history habitat requirements. Tools, such as life-cycle models, are needed to manage population dynamics and quantify the composite effects of processes across space and time. Nevertheless, complex institutions can hinder result interpretation and communication, and limit model use in decision-making. This paper focuses on the federally endangered and endemic Sacramento River winter-run Chinook (Oncorhynchus tshawytscha) by developing a Winter-Run Habitat-based Population Model (WRHAP). Additionally, WRHAP is a conceptual, freshwater rearing stage model that includes alternative rearing habitats reported in the literature (e.g., floodplains, off-channel and tributaries), defines rules of habitat use based on instream conditions, and incorporates a juvenile growth module that combines bioenergetics modeling with empirical growth rates. Model outputs reasonably follow observed out-migration patterns and provide a realistic smolt size distribution arriving at the San Francisco Bay. This effort demonstrates the importance of currently “non-critical” habitats (as defined by the Endangered Species Act) for juvenile development, with floodplain habitat contributing to a quarter of out-migrating biomass (despite < 18% availability and <10% of total rearing days), and off-channel growth being one of the most sensitive parameters (explaining ~13% of average juvenile weight variance). The model shows the utility of a simple population model to explore relationships between habitat quality/quantity and juvenile development, and to assist water/environmental management and decision-making processes focused on species recovery.

54 ENVIRONMENTAL SCIENCES↗

Sensor enabled data-driven predictive analytics for modeling and control with high penetration of DERs in distribution systems

The electric power grid is undergoing a tremendous transformation due to the increasing penetration of renewable energy resources beginning with wind and more recently with the distributed energy resources (DERs) such as solar and battery storage. DERs have dramatically changed the role of the distribution systems in the overall power grid, and they are expected to contribute a significant portion of power generation in the future. If current trends for DERs continue, system operation and control will need to change dramatically for improved grid reliability and resiliency. As renewable resources increase in penetration, new and challenging operational, planning, and design problems are expected to emerge. Some of the key challenges that arise in the planning and operation of the future grid are: 1) Quantifying the impact of high DER penetration in distribution systems on bulk grid behavior over multiple time scales. 2) Identifying whether a particular DER configuration/settings have a large impact on the overall grid behavior. These challenges can be addressed in an offline manner using detailed T&D grid models and they can also be addressed in an online manner using sensor measurements. In particular, the advancement and planned growth in sensor technology in power grid over various voltage levels provide us with a unique opportunity to tackle these challenges from a data analytic perspective without needing detailed T&D grid models. A few questions that naturally arise when addressing the challenges from DERs using sensor data are: 1) How can we use limited sensor measurements to monitor & control voltage stability and small signal stability of the bulk system? 2) How can we ensure that the developed data analytic methods are robust to data availability and quality issues? 3) How can we compute the developed analytics in a scalable manner using streaming measurements? In this project, we addressed the aforementioned challenges arising from DERs and answered the questions raised above on how to effectively use the sensor measurements to enhance the reliability and performance of the electric grid. Thus, the overarching goal of this project is to develop effective reduced/representative system models from data that make the computational complexity sufficiently manageable so as to be useful to simulate, analyze, and even control complex non-linear power systems dynamics with large penetrations of DERs. In order to achieve the objective, the project team established a four-fold technical approach 1) Formulated a combined transmission-distribution co-simulation framework for data generation and validation, 2) Derived reduced/representative models of power systems based on data-driven methods for efficient computation and appropriate representation of system behavior, 3) Developed data driven characterization of power system behavior based on transfer operator theory, machine learning and optimization for model estimation, 4) Incorporated a scalable data management and processing architecture using distributed Kafka streaming applications that coordinate input data streams to the developed data analytics. The key accomplishments of the project are: 1) Development of a scalable multi-timescale T&D co-simulation framework (both for steady state and for dynamic co-simulation) using commercial solvers (PSSE and GridLAB-D). The steady-state T&D co-simulation interface is shared with our industry partner (PJM). 2) A structured reduced order dynamic model of distribution systems that can represent partial motor stalling along with a systematic procedure to derive the model parameters. 3) A PMU based online method to monitor, localize and mitigate fault-induced delayed voltage recovery using DER reactive support and load control in distribution systems. 4) Development of linear operator based robust methodologies for dynamic state estimation, uncertainty quantification, system identification and trajectory prediction for power system dynamics. 5) An adaptive damping control for utilizing wind energy resources to provide oscillation damping and system stability. 6) Implementation of Kafka-based framework for efficient processing of streaming data using Linux-based local virtual environment.

DER integration↗

Fast Adaptive Neural Control of Resonant Extraction at Fermilab

We present progress on the development of a machine learning (ML) regulation system for third-order resonant extraction of the beam delivered to the Mu2e experiment at Fermilab. We consider classical and ML-based controllers optimized on semi-analytic simulations and provide performance comparisons for several models. Additionally, we discuss the efficiency of each model in training, which has implications for future work on adaptive control. We also discuss progress on developing optimized implementations of ML models for edge-based inference.

Whitbeck, A. [Fermilab] (ORCID:0000000342245164)↗

Report on the LLNL Global Full-waveform Inversion Workflow and Progress

LLNL has developed the SPiRaL global seismic tomography model based upon millions of body wave travel times and surface wave dispersion curves in the period range of 25-200 seconds (Simmons et al., 2021). The motivation to construct the SPiRaL model (which stands for S-, P-, Rayleigh, and Love waves) is based upon the guidance provided by a 3-D model workshop held in Berkeley in 2007 and documented in a paper by Zucca et al. (2009). It was concluded in Zucca et al. (2009) that a “practical and obtainable” model is one that predicts body wave travel times for event location while also predicting surface wave velocities for moment tensor estimation. We have largely followed that guidance which has led to the construction of SPiRaL. Looking forward, we now seek to improve the global model through full waveform inversion (FWI). This necessitates the development/customization of an efficient and flexible modeling workflow to automate the process using LLNL’s supercomputing platforms. This report briefly outlines the work that being performed at LLNL to develop the necessary tools and framework to update the SPiRaL global model with full waveform information.

58 GEOSCIENCES↗