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,279 records · Page 71

Physics-constrained, low-dimensional models for magnetohydrodynamics: First-principles and data-driven approaches

Plasmas are highly nonlinear and multiscale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD), although these reduced models hold promise for understanding key physical mechanisms, efficient computation, and real-time optimization and control. Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, and data-driven models, obtained by modern machine learning and system identification, can furnish this gap in the lower levels of the model hierarchy. This work develops a reduced-order modeling framework for compressible plasmas, leveraging decades of progress in projection-based and data-driven modeling of fluids. We begin by formalizing projection-based model reduction for nonlinear MHD systems. To avoid separate modal decompositions for the magnetic, velocity, and pressure fields, we introduce an energy inner product to synthesize all of the fields into a dimensionally consistent, reduced-order basis. Next, we obtain an analytic model by Galerkin projection of the Hall-MHD equations onto these modes. We illustrate how global conservation laws constrain the model parameters, revealing symmetries that can be enforced in data-driven models, directly connecting these models to the underlying physics. We demonstrate the effectiveness of this approach on data from high-fidelity numerical simulations of a three-dimensional spheromak experiment. Finally, this manuscript builds a bridge to the extensive Galerkin literature in fluid mechanics and facilitates future principled development of projection-based and data-driven models for plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developing a Reduced 240-Bus WECC Dynamic Model for Frequency Response Study of High Renewable Integration

The ongoing changes in the generation resource mix, driven by the rapid adoption of inverter-based resources (IBR) as well as the early retirement of synchronous generators, are bringing new challenges to the planning and operation of bulk electric power systems. Consequently, there is an increasing need to understand, design, and quantify the reliability service provision from IBRs by performing integrated scheduling and dynamic simulations. However, test systems that have consistent scheduling and dynamic models rarely exist largely because of the decoupled nature of the two simulations on a synchronous generator-dominated system. This paper develops the dynamic model of a reduced Western Electricity Coordinating Council (WECC) system. In conjunction with the existing scheduling model, it is suitable for integrating scheduling and dynamic simulations. The reduced 240-bus WECC model reflects the generation resource mix of the Western Interconnection as of 2018. Moreover, the developed dynamic model is validated against field frequency events measured by FNET/GridEye and preserves the dominant inter-area oscillation mode in WECC.

dynamic simulation↗

Non-Idealities in Lab-Scale Kinetic Testing: A Theoretical Study of a Modular Temkin Reactor

The Temkin reactor can be applied for industrial relevant catalyst testing with unmodified catalyst particles. It was assumed in the literature that this reactor behaves as a cascade of continuously stirred tank reactors (CSTR). However, this assumption was based only on outlet gas composition or inert residence time distribution measurements. The present work theoretically investigates the catalytic CO2 methanation as a test case on different catalyst geometries, a sphere, and a ring, inside a single Temkin reaction chamber under isothermal conditions. Axial gas-phase species profiles from detailed computational fluid dynamics (CFD) are compared with a CSTR and 1D plug-flow reactor (PFR) model using a sophisticated microkinetic model. In addition, a 1D chemical reactor network (CRN) model was developed, and model parameters were adjusted based on the CFD simulations. Whereas the ideal reactor models overpredict the axial product concentrations, the CRN model results agree well with the CFD simulations, especially under low to medium flow rates. This study shows that complex flow patterns greatly influence species fields inside the Temkin reactor. Although residence time measurements suggest CSTR-like behavior, the reactive flow cannot be described by either a CSTR or PFR model but with the developed CRN model.

Wehinger, Gregor D. (ORCID:0000000217743391)↗

Machine learning for the redox potential prediction of molecules in organic redox flow battery

Here, organic redox flow batteries (ORFB) are recognized as an innovative technology for the large-scale storage of renewable energy. The redox potential of organic redox-active molecules plays a vital role in their performance. Advanced screening techniques like high-throughput experiment and machine learning (ML) have significantly enhanced organic material performance and transformed the field of ORFB. However, the scarcity of experimental data poses a considerable challenge for ML model development in this domain. In our study, we developed lightweight graph-based Gaussian process regression (GPR) models with GPU-accelerated marginalized graph kernel and hybrid kernel to predict the redox potentials of organic redox-active molecules for ORFBs, specifically focusing on small datasets. To evaluate model accuracy, we created a new experimental database of organic redox-active molecules by the data from hundreds of published papers and assembled previous computational datasets. We also considered some key parameters, such as pH conditions and solvent type, to assess their impact on redox potential prediction. Our GPR model predicted redox potentials with high accuracy across all datasets using minimal training data. The study provides powerful tools for molecule screening and design and delivers valuable guidance on designing training datasets for costly experiments.

25 ENERGY STORAGE↗

Evaluation of Reduced Order Model for HT-9 Creep and Modifications to Current HT-9 Creep Model in BISON

In order to leverage existing reduced order models (ROM) for the modeling of LANL-developed HT-9 cladding constitutive behavior, mechanistic-based models have been implemented into BISON. It is posited that the increased fidelity of the ROM will show improvement over engineering-scale models in terms of comparison against experimental measurements. INL will assist LANL in implementing an HT9 ROM in BISON. Once the HT-9 mechanistic constitutive model can be leveraged in nuclear performance simulations, the numerical results of fast reactor models will be analyzed and compared against the Fuels Irradiation & Physics Database (FIPD) and separate effects tests to attempt a validation process. INL and LANL will work together to perform any necessary improvement that is identified during the implementation and use processes. Successful completion of the milestone will enable the general use of validated mechanistic ROMs for HT-9 cladding in BISON. This report provides a brief introduction to the Los Alamos Reduced Order Model Applied to Nonlinear Constitutive Equations (LAROMANCE) code, mechanical testing used to illustrate its creep and plastic deformation predictions for HT-9 cladding, and a comparison with existing models implemented into BISON. Current models have been modified to more appropriately account for primary thermal creep and results for corresponding testing are included in this report.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of a Virtual CFR Engine Model for Knocking Combustion Analysis

Knock is a major bottleneck to achieving higher thermal efficiency in spark ignition (SI) engines. The overall tendency to knock is highly dependent on fuel anti-knock quality as well as engine operating conditions. It is, therefore, critical to gain a better understanding of fuel-engine interactions in order to develop robust knock mitigation strategies. In the present work, a numerical model based on three-dimensional (3-D) computational fluid dynamics (CFD) was developed to capture knock in a Cooperative Fuel Research (CFR) engine. For combustion modeling, a hybrid approach incorporating the G-equation model to track turbulent flame propagation, and a homogeneous reactor multi-zone model to predict end-gas auto-ignition ahead of the flame front and post-flame oxidation in the burned zone, was employed. In addition, a hybrid methodology was implemented wherein a laminar flame speed lookup table generated a priori from a chemical kinetic mechanism could be used to provide flame speed as an input to the G-equation model, instead of using conventional empirical correlations. Multi-cycle Reynolds-averaged Navier-Stokes (RANS) simulations were performed for two different spark timings (STs) corresponding to non-knocking and knocking conditions, with other operating conditions kept the same as those of a standard research octane number (RON) test. Isooctane was considered as the fuel for the numerical study. Two different reduced kinetic mechanisms were employed to describe end-gas auto-ignition chemistry and to generate the flame speed lookup table. Experimental data, including intake/exhaust boundary conditions, was provided by an isooctane ST sweep study conducted in an in-house CFR engine. Moreover, cylinder wall/valve/port surface temperatures and residual gas fraction (RGF) were estimated using a well-calibrated one-dimensional (1-D) model. On the other hand, a novel methodology was also developed to analyze experimental data for the knocking case and identify the most representative cycle. For the non-knocking case, a good agreement was found between experiment and CFD simulation, with respect to cycle-averaged values of 10% burn point (CA10), 50% burn point (CA50), and peak pressure magnitude/location. The virtual CFR engine model was also demonstrated to be capable of predicting average knock characteristics for the knocking case, such as knock point, knock intensity, and energy of resonance, with good accuracy.

42 ENGINEERING↗

SMART Task 6: Evaluation of the Costs of Geologic CO2 Storage for the Illinois Basin Decatur Project Site Using the NRAP/SMART Technoeconomic and Liability Evaluation for Storage (TALES) Model

This is a presentation featuring an analysis related to SMART Task 6 in which CO2 storage costs are presented. The National Energy Technology Laboratory has developed the NRAP/SMART Technoeconomic and Liability Evaluation for Storage (TALES) model to provide quantitative cost-based insights to support developers planning CO2 injection and storage projects. TALES calculates the revenues, costs, and financial performance of candidate CO2 saline storage project based on site-specific activity costs and financial parameters. TALES is being integrated as a module pertaining to storage cost as part of the broader SMART Visualization and Decision Support Platform (SVDSP). In this study, the TALES model was applied using real activity cost data associated with the development and operations at the Illinois Basin Decatur Project (IBDP) CO2 storage project site. Scenario analysis was implemented in which crucial operational and cost attributes were varied and the associated cost implications observed. Key results data and project cost summary metrics like first-year breakeven price of CO2 ($/tonne) and net present value (NPV) are presented in similar fashion to how they will appear in the SVDSP.

Vikara, Derek↗

Reactive molecular dynamics models from ab initio molecular dynamics data using relative entropy minimization

Here, we present two new multiscale molecular dynamics (MS-RMD) models for the hydrated excess proton in water developed directly from ab initio molecular dynamics (AIMD) simulation data of the same system. The potential of mean force along the proton transfer reaction coordinate and radial distribution functions for the MS-RMD models are shown faithfully reproduce those of AIMD. The models are developed using an algorithm based on relative entropy minimization, thus demonstrating the ability of the method to rapidly generate accurate and highly efficient reactive MD force fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Localization landscape theory of disorder in semiconductors. I. Theory and modeling

Here, we present here a model of carrier distribution and transport in semiconductor alloys accounting for quantum localization effects in disordered materials. This model is based on the recent development of a mathematical theory of quantum localization which introduces for each type of carrier a spatial function called localization landscape. These landscapes allow us to predict the localization regions of electron and hole quantum states, their corresponding energies, and the local densities of states. We show how the various outputs of these landscapes can be directly implemented into a drift-diffusion model of carrier transport and into the calculation of absorption/emission transitions. This creates a new computational model which accounts for disorder localization effects while also capturing two major effects of quantum mechanics, namely, the reduction of barrier height (tunneling effect) and the raising of energy ground states (quantum confinement effect), without having to solve the Schrödinger equation. Finally, this model is applied to several one-dimensional structures such as single quantum wells, ordered and disordered superlattices, or multiquantum wells, where comparisons with exact Schrödinger calculations demonstrate the excellent accuracy of the approximation provided by the landscape theory.

36 MATERIALS SCIENCE↗

A Hybrid Framework Combining Model-Based and Data-Driven Methods for Hierarchical Decentralized Robust Dynamic State Estimation: Preprint

This paper combines model-based and data-driven methods to develop a hierarchical decentralized robust dynamic state estimator (DSE). A two-level hierarchy is proposed, where the lower level is comprised of robust model-based decentralized DSEs. The state estimates sent from the lower level are received at the upper level, where they are filtered by a robust data- driven DSE, after a principled sparse selection. This selection allows us to shrink the dimension of the problem at the upper level, and hence to speed up the computational time significantly. The proposed hybrid framework does not depend upon the centralized infrastructure of the control centers, thus it can be completely embedded into the wide-area measurement systems. This feature will ultimately facilitate the placement of hierarchical decentralized control schemes at the phasor data concentrator (PDC) locations. Also, the network model is not necessary, thus a topology processor is not required. Finally, there is no assumption on the dynamics of the electric loads. The proposed framework is tested on a 2000-bus network.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Evaluating SWAT + model uncertainties for human and natural outcomes: Application in a Great Lakes agricultural watershed

Nutrient exports from agricultural lands in the Great Lakes Region pose significant threats to water quality and ecological health through eutrophication, hypoxia, and harmful algal blooms. Climate change and agricultural adaptation practices complicate future nutrient loading due to intensified hydrologic cycles and land use decisions. Our research focuses on evaluating the Soil and Water Assessment Tool (SWAT) plus model parametric uncertainties for human and natural outcomes across different scales. These factors are integral to ensuring a balance between productive agricultural practices and maintaining the health of watershed hydrology. However, uncertainties in modeling such complex interactions pose significant challenges, limiting our ability to precisely determine critical factors that influence crop yield and soil moisture. Our analysis employs Sobol global sensitivity analysis to evaluate first-order, second order, and total-order indices for SWAT crop growth parameter, ensuring comprehensive assessment of individual and interactive effects on model outputs. The objective is to identify the parameters that significantly affect model outputs for crop yield and soil moisture and improve our understanding of their interactions at the basin and hydrological response unit (HRU) scale. Our case study, the Portage River Watershed, which drains into Lake Erie, is chosen to better capture finer scale interactions crucial for predicting nutrient loading under future climate scenarios. This foundational work is aimed at setting the stage for the future development of an agent-based model (ABM). The ABM model would incorporate SWAT outputs to dynamically simulate decision-making processes.

Bunyon, Enock↗

Ion-Exchange Interdiffusion Model with Potential Application to Long-Term Nuclear Waste Glass Performance

Ion exchange is an integral mechanism influencing the corrosion of glasses. Due to the formation of alteration layers in aqueous conditions, it is difficult to conclusively deconvolute the process of ion exchange from other processes, principally dissolution of the glass matrix. Therefore, we have developed a method to isolate alkali diffusion that involves contacting glass coupons with a solution of 6 LiCl dissolved in functionally inert dimethyl sulfoxide. Here, we employ the method at temperatures ranging from 25 to 150 °C with various glass formulations. Glass compositions include simulant nuclear waste glasses, such as SON68 and the international simple glass (ISG), glasses in which the nature of the alkali element was varied, and glasses that contained more than one alkali element. An interdiffusion model based on Fick’s second law was developed and applied to all experiments to extract diffusion coefficients. The model expands established models of interdiffusion to the case where multiple types of alkali sites are present in the glass. Activation energies for alkali ion exchange were calculated and the results are in agreement with those obtained in glass strengthening experiments but are nearly five times higher than values reported for diffusion-controlled processes in nuclear waste glass corrosion experiments. A discussion of the root causes for this apparent discrepancy is provided. The interdiffusion model derived from laboratory experiments is expected to be useful for modeling glass corrosion in a geological repository when the silicon concentration is high.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Learning Interactions in Complex Biological Systems

In this project, we developed novel methods for integrating simulation and data-driven methods to accelerate biomedical discovery. In particular, we exploit agent-based models of disease processes and deep reinforcement learning to identify potential novel multi-drug, patient-specific, and adaptive therapeutic strategies. We use the Innate Immune Response agent-based model (IIRABM) to model sepsis as a demonstration of our approach. We first calibrate the agent-based model using novel stochastic optimization algorithms develop in this project, then use the calibrated model to identify optimal cytokine mediation strategies using deep reinforcement learning (DRL). The learned policy achieves a dramatic reduction in mortality over a set of 500 simulated patients relative to standalone antibiotic therapy. Advantages of our approach are three-fold: 1) the use of simulation allows exploring therapeutic strategies beyond clinical practice and available data, 2) advances in DRL accommodate learning complex therapeutic strategies for complex biological systems, and 3) optimized treatments respond to a patient’s individual disease progression over time, therefore capturing both differences across patients and the inherent randomness of disease progression within a single patient. We hope that this work motivates both considering adaptive, personalized multi-cytokine mediation therapy for sepsis and exploiting simulation with DRL for precision medicine more broadly.

59 BASIC BIOLOGICAL SCIENCES↗

Risk assessment of safety-critical data communication in digital safety feature control system - 205

As one of the safety-critical systems in Advanced Power Reactor-1400 (APR-1400) nuclear power plant (NPP), the Engineered Safety Feature-Component Control System (ESF-CCS) employs safety data link and network communication for transmitting safety component actuation to effectively facilitate various safety-critical field controllers. However, data communication failure risk in the ESF-CCS has yet to be fully quantified, the probability that a safety-critical system in digitalized NPP becomes unsafe due to a network failure must be evaluated to quantify the risk of digital I and C system. In this study, a fault tree model was developed to assess the data link and data network failure-induced unavailability of a safety-critical system function used to generate an automated control signal for actuating design basis accident mitigation equipment. Especially, the risk of data communication failure in a digital safety feature control system was analyzed in consideration of interconnection between safety-critical controllers and the fault-tolerant algorithm implemented in the network configuration. Based on the developed fault tree model, case studies were performed to quantitatively assess the unavailability of ESF-CCS signal generation due to safety data link and network failure and their risk effect on safety signal generation unavailability. This study is expected to provide risk information on the safety-critical data communication in a digitalized NPP instrumentation and control system. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

U(VI) bioreduction with emulsified vegetable oil as the electron donor-Model application to a field test

A one-time 2-hour emulsified vegetable oil (EVO) injection in a fast flowing aquifer decreased U discharge to a stream for over a year. Using a comprehensive biogeochemical model developed in the companion article based on microcosm tests, we approximately matched the observed acetate, nitrate, Fe, U, and sulfate concentrations, and described the major evolution trends of multiple microbial functional groups in the field test. While the lab-determined parameters were generally applicable in the field-scale simulation, the EVO hydrolysis rate constant was estimated to be an order of magnitude greater in the field than in the microcosms. The model predicted substantial biomass (sulfate reducers) and U(IV) accumulation near the injection wells and along the side boundaries of the treatment zone where electron donors (long-chain fatty acids) from the injection wells met electron acceptors (sulfate) from the surrounding environment. While EVO retention and hydrolysis characteristics were expected to control treatment longevity, modeling results indicated that electron acceptors such as sulfate may not only compete for electrons but also play a conducive role in degrading complex substrates and enhancing U(VI) reduction and immobilization. As a result, the spacing of the injection wells could be optimized for effective sustainable bioremediation.

Tang, Guoping↗

Hepatocyte-based in vitro model for assessment of drug-induced cholestasis

Early detection of drug-induced cholestasis remains a challenge during drug development. We have developed and validated a biorelevant sandwich-cultured hepatocytes- (SCH) based model that can identify compounds causing cholestasis by altering bile acid disposition. Human and rat SCH were exposed (24–48 h) to known cholestatic and/or hepatotoxic compounds, in the presence or in the absence of a concentrated mixture of bile acids (BAs). Urea assay was used to assess (compromised) hepatocyte functionality at the end of the incubations. The cholestatic potential of the compounds was expressed by calculating a drug-induced cholestasis index (DICI), reflecting the relative residual urea formation by hepatocytes co-incubated with BAs and test compound as compared to hepatocytes treated with test compound alone. Compounds with clinical reports of cholestasis, including cyclosporin A, troglitazone, chlorpromazine, bosentan, ticlopidine, ritonavir, and midecamycin showed enhanced toxicity in the presence of BAs (DICI ≤ 0.8) for at least one of the tested concentrations. In contrast, the in vitro toxicity of compounds causing hepatotoxicity by other mechanisms (including diclofenac, valproic acid, amiodarone and acetaminophen), remained unchanged in the presence of BAs. A safety margin (SM) for drug-induced cholestasis was calculated as the ratio of lowest in vitro concentration for which was DICI ≤ 0.8, to the reported mean peak therapeutic plasma concentration. SM values obtained in human SCH correlated well with reported % incidence of clinical drug-induced cholestasis, while no correlation was observed in rat SCH. This in vitro model enables early identification of drug candidates causing cholestasis by disturbed BA handling. - Highlights: • Novel in vitro assay to detect drug-induced cholestasis • Rat and human sandwich-cultured hepatocytes (SCH) as in vitro models • Cholestatic compounds sensitize SCH to toxic effects of accumulating bile acids • Drug-induced cholestasis index (DICI) as measure of a drug's cholestatic signature • In vitro findings correlate well with clinical reports on cholestasis.

60 APPLIED LIFE SCIENCES↗

MCNP6 enhancements of delayed-particle production

Over the last decade, there has been an increased interest in the production of delayed-particle signatures from neutron and photon interactions with matter. To address this interest, various radiation transport codes have developed a wide range of delayed-particle physics packages. With the recent merger of the Monte Carlo transport codes MCNP5 and MCNPX, MCNP6 inherited the comprehensive model-based delayed-particle production capabilities developed in MCNPX over the last few years. An integral part of this capability consists of the depletion code CINDER90 which was incorporated into MCNPX in 2004. During this last year, significant improvements have been made to the MCNP6 physics and algorithms associated with delayed-particle production, including the development of a delayed-beta capability, an algorithm enhancement for the delayed-neutron treatment, and a database enhancement for delayed-gamma emission. The delayed-beta feature represents an important component in modeling background signals produced by active interrogation sources. Combined, these improvements provide MCNP6 with a flexible state-of-the-art physics package for generating high-fidelity signatures from fission and activation. This paper provides details of these enhancements and presents results for a variety of fission and activation examples. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗