Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

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,135 records · Page 63

One Dimensional Analysis Model of a Condensing Spray Chamber Including Rocket Exhaust Using SINDA/FLUINT and CEA

Modeling droplet condensation via CFD codes can be very tedious, time consuming, and inaccurate. CFD codes may be tedious and time consuming in terms of using Lagrangian particle tracking approaches or particle sizing bins. Also since many codes ignore conduction through the droplet and or the degradating effect of heat and mass transfer if noncondensible species are present, the solutions may be inaccurate. The modeling of a condensing spray chamber where the significant size of the water droplets and the time and distance these droplets take to fall, can make the effect of droplet conduction a physical factor that needs to be considered in the model. Furthermore the presence of even a relatively small amount of noncondensible has been shown to reduce the amount of condensation. It is desirable then to create a modeling tool that addresses these issues. The path taken to create such a tool is illustrated. The application of this tool and subsequent results are based on the spray chamber in the Spacecraft Propulsion Research Facility (B2) located at NASA's Plum Brook Station that tested an RL-10 engine. The platform upon which the condensation physics is modeled is SINDAFLUINT. The use of SINDAFLUINT enables the ability to model various aspects of the entire testing facility, including the rocket exhaust duct flow and heat transfer to the exhaust duct wall. The ejector pumping system of the spray chamber is also easily implemented via SINDAFLUINT. The goal is to create a transient one dimensional flow and heat transfer model beginning at the rocket, continuing through the condensing spray chamber, and finally ending with the ejector pumping system. However the model of the condensing spray chamber may be run independently of the rocket and ejector systems detail, with only appropriate mass flow boundary conditions placed at the entrance and exit of the condensing spray chamber model. The model of the condensing spray chamber takes into account droplet conduction as well as the degrading effect of mass and heat transfer due to the presence of noncondensibles. The one dimension model of the condensing spray chamber makes no presupposition on the pressure profile within the chamber, allowing the implemented droplet physics of heat and mass transfer coupled to the SINDAFLUINT solver to determine a transient pressure profile of the condensing spray chamber. Model results compare well to the RL-10 engine pressure test data.

Chemical Equilibrium↗

One Dimensional Analysis Model of a Condensing Spray Chamber Including Rocket Exhaust Using SINDA/FLUINT and CEA

Modeling droplet condensation via CFD codes can be very tedious, time consuming, and inaccurate. CFD codes may be tedious and time consuming in terms of using Lagrangian particle tracking approaches or particle sizing bins. Also since many codes ignore conduction through the droplet and or the degradating effect of heat and mass transfer if noncondensible species are present, the solutions may be inaccurate. The modeling of a condensing spray chamber where the significant size of the water droplets and the time and distance these droplets take to fall, can make the effect of droplet conduction a physical factor that needs to be considered in the model. Furthermore the presence of even a relatively small amount of noncondensible has been shown to reduce the amount of condensation [Ref 1]. It is desirable then to create a modeling tool that addresses these issues. The path taken to create such a tool is illustrated. The application of this tool and subsequent results are based on the spray chamber in the Spacecraft Propulsion Research Facility (B2) located at NASA's Plum Brook Station that tested an RL-10 engine. The platform upon which the condensation physics is modeled is SINDAFLUINT. The use of SINDAFLUINT enables the ability to model various aspects of the entire testing facility, including the rocket exhaust duct flow and heat transfer to the exhaust duct wall. The ejector pumping system of the spray chamber is also easily implemented via SINDAFLUINT. The goal is to create a transient one dimensional flow and heat transfer model beginning at the rocket, continuing through the condensing spray chamber, and finally ending with the ejector pumping system. However the model of the condensing spray chamber may be run independently of the rocket and ejector systems detail, with only appropriate mass flow boundary conditions placed at the entrance and exit of the condensing spray chamber model. The model of the condensing spray chamber takes into account droplet conduction as well as the degrading effect of mass and heat transfer due to the presence of noncondensibles. The one dimension model of the condensing spray chamber makes no presupposition on the pressure profile within the chamber, allowing the implemented droplet physics of heat and mass transfer coupled to the SINDAFLUINT solver to determine a transient pressure profile of the condensing spray chamber. Model results compare well to the RL-10 engine pressure test data.

Condensing Spray Model↗

Uncertainty Quantification of Capacitor Switching Transient Location using Machine Learning

Identification of capacitor switching transient location provides valuable insight into the state of the associated equipment. Machine learning (ML) models, and convolutional neural networks (CNNs) in particular, have demonstrated remarkable performance in signal location. However, ML models are data driven whose predictions are affected by noise in data and may also suffer from large extrapolation errors when applied to new conditions. Uncertainty quantification (UQ) is necessary to ensure model trustworthiness and avoid overconfident predictions in extrapolation. Here, in this work, we propose a novel UQ method, called PI3NN, to quantify prediction uncertainty of ML models and integrate the method with CNNs for transient source location. PI3NN calculates Prediction Intervals by training 3 Neural Networks and uses root-finding methods to determine the interval precisely. Additionally, PI3NN can identify out-of-distribution (OOD) data in a nonstationary condition to avoid overconfident prediction. Results indicate that with PI3NN, transient signals are not only correctly identified, but when said signals are subject to corruptions characteristic of an actual power monitoring system (e.g. non-ideal sensors), the model recognizes when it is uncertain about its predictions, effectively letting the user know when to accept or discard the results.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Climatic factors and human population changes in Eurasia between the Last Glacial Maximum and the early Holocene

Archaeological records document a significant expansion of populations from the Last Glacial Maximum (LGM, ~23–19 ka) to the early Holocene (EH, ~9 ka) in Eurasia, which is often attributed to the influence of orbital-scale climate changes. Yet, information remains limited concerning the climatic factor(s) which were responsible for conditioning demographic patterns. Here, in this work, we present results from an improved Minimalist Terrestrial Resource Model (MTRM), forced by a transient climate simulation from the LGM to the EH. Simulated potential hunter-gatherer population densities and spatial distributions across Eurasia are supported by observed archaeological sites in Europe and China. In the low latitudes, potential population size change was predominantly controlled by precipitation and its strong influence on plant and animal resources. In the middle-high latitudes, temperature was the dominant driver in influencing potential population size change and animal resource availability. Different regional responses of potential populations to climate change across Eurasia - owing to variations in available food resources between the LGM and EH - provide a better understanding of human dispersal during the Late Pleistocene.

54 ENVIRONMENTAL SCIENCES↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Uncertainty Quantification in High-Low Dynamic System Coupling using RAVEN and TRANSFORM

This work demonstrates new functionality and applications stemming from the development of high-fidelity to low-fidelity (high-low) coupling for system simulations and to further explore the capabilities of the Risk Analysis Virtual Environment (RAVEN) in the performance of uncertainty quantification in this kind of high-low coupled system models. The work builds from previous work on high-low coupling that utilized COBRA-TF (CTF), the high-fidelity subchannel analysis code, with a low fidelity model built in ORNL’s TRANSFORM, the system analysis code, utilizing the Functional Mock-Up Interface (FMI). Steady-state and transient analysis examples using the high/low coupled models generated from CTF and TRANSFORM/FMI are investigated. The workflows for both steady-state and transient coupled simulations are described. A steady-state parameter sweep and uncertainty analysis of the primary flow rates and reactor power are demonstrated. Likewise, a transient pump trip and power ramp sensitivity studies are also demonstrated. This work elucidates some of the potential benefits and future needs of using RAVEN for high/low system coupling analysis of energy systems. It also shows some of the difficulties that can be encountered in coupling system simulations.

Williams, Wesley↗

Comparison of observed rheological properties of hard wheat flour dough with predictions of the Giesekus-Leonov, White-Metzner and Phan-Thien Tanner models

The measured rheological behavior of hard wheat flour dough was predicted using three nonlinear differential viscoelastic models. The Phan-Thien Tanner model gave good zero shear viscosity prediction, but overpredicted the shear viscosity at higher shear rates and the transient and extensional properties. The Giesekus-Leonov model gave similar predictions to the Phan-Thien Tanner model, but the extensional viscosity prediction showed extension thickening. Using high values of the mobility factor, extension thinning behavior was observed but the predictions were not satisfactory. The White-Metzner model gave good predictions of the steady shear viscosity and the first normal stress coefficient but it was unable to predict the uniaxial extensional viscosity as it exhibited asymptotic behavior in the tested extensional rates. It also predicted the transient shear properties with moderate accuracy in the transient phase, but very well at higher times, compared to the Phan-Thien Tanner model and the Giesekus-Leonov model. None of the models predicted all observed data consistently well. Overall the White-Metzner model appeared to make the best predictions of all the observed data.

Non-NASA Center↗

Analysis of AP1000 Small-Break Loss-of-Coolant Accident Using Reactor Transient Simulator

The Westinghouse Electric Company’s Advanced Passive Reactor (AP1000) is characterized by the incorporation of passive safety systems (PSSs) designed to ensure core cooling during transient events. The assessment of PSSs requires evaluation of their performance through a combination of experiments and simulations employing various thermal-hydraulic codes. In addition, detailed evaluation of PSSs for a specific reactor system transient analysis such as loss-of-coolant-accident analysis supports understanding representative integral effects test facility development and the further evolution model development and assessment process. Developing a reactor system code is a complex and time-consuming process that requires significant engineering expertise and effort. It can take several months to even years to complete in the early stages of reactor system design and analysis. However, this process can be expedited through the use of transient simulator models for similar reactor systems, which can be used for lesson learning and training purposes. This study uses the Personal Computer Transient Analyzer (PCTRAN) code. The main advantage of PCTRAN is its ease of use and ability to run faster than real time. This study presents the results obtained for a small-break loss-of-coolant accident (SBLOCA) for two breaks using the full version (licensed) of PCTRAN. The purpose of this investigation is to evaluate the overall system behavior during the postulated SBLOCA event as well as assess the capability of the PCTRAN code to reproduce the system response during transient events. The obtained results were compared with the Westinghouse NOTRUMP system code. The PCTRAN code proved to be reliable in predicting the qualitative behavior of the system in both transient cases. As for the system response, it was found that it is contingent on the activation time of the PSSs. The differences in reactor coolant system pressure between the two codes were attributed to the critical flow model and simplification of mass and energy balance. Despite PCTRAN’s limitations, it can still provide a reasonable prediction of various reactor parameters such as pressure, mass flow rate, and void fraction during a SBLOCA scenario. It is worth noting that PCTRAN currently employs a bulk approach similar to that of the Modular Accident Analysis Program (MAAP) and MELCOR codes. However, the upcoming version of PCTRAN will include an artificial intelligence–based detection and accident prevention system, as well as different models for different reactor components. Consequently, PCTRAN has the potential to be upgraded to match the system thermal-hydraulic codes of the U.S. Nuclear Regulatory Commission and become more widely used in cybersecurity to safeguard nuclear power plants from cyberattacks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Autonomous Inverter Controls for Resilient and Secure Grid Operation: Vector Control Design for Grid Forming

The project addresses both fundamental and practical challenges of GFM/GFL inverter control for the power grids with high inverter based resources (IBRs) penetration. A data- driven modeling technique is applied to accurately model dynamics of PWM inverters, including electromagnetic-transient (EMT). Systematic and integrative designs of grid- forming (GFM) and grid-following (GFL) primary controls are developed to guarantee system performance under either normal or abnormal operating conditions without violating constraints. This modeling and control framework provides black-start capability in case of an outage without relying on rotating generators, and its secondary control is also shown to enhance resilience against cyber-physical attacks.

14 SOLAR ENERGY↗

Assessing Dynamic Behaviors in Converter- Dominated Power Systems via RMS and EMT Simulations: A Study of Hawaii’s NELHA Microgrid

The transition from conventional power systems to converter-based microgrids has significantly advanced sustainability, clean energy integration, and operational reliability. However, this paradigm shift introduces operational challenges due to the intermittent nature of renewable energy sources and the non-linear characteristics of power electronic loads, inducing voltage fluctuations and harmonic distortions that complicate voltage and frequency regulation. Accurate dynamic modeling is hypothesized to be critical for capturing such effects, enabling reliable simulation and control strategy development. This study introduces an innovative dynamic modeling framework for a real-world converter-based microgrid, utilizing both root mean square (RMS) and electromagnetic transient (EMT) simulation methods. The microgrid was modeled in DIgSILENT PowerFactory, with simulations calibrated against high-resolution field measurements from SEL-735 power quality meters. Results show that RMS simulations effectively characterize steady-state dynamics, while EMT simulations are essential for capturing high-frequency transients and non-linear effects from photovoltaic inverters and variable frequency drives (VFDs). This complementary approach provides a comprehensive understanding of microgrid behavior, providing critical insights for improving simulation accuracy, advancing protection schemes, and improving resilience in future low-inertia power networks.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Resident and transient coyotes exhibit differential patterns of movement behavior across heterogeneous landscapes in the southeastern United States

Coyotes (Canis latrans) are a highly adaptable canid species whose behavioral plasticity has allowed them to persist in a wide array of habitats throughout North America. As generalists, coyotes can alter movement patterns and change territorial strategies between residency (high site fidelity) and transiency (low site fidelity) to maximize fitness. Uncertainty remains about resident and transient coyote movement patterns and habitat use because research has reached conflicting conclusions regarding patterns of habitat use by both groups. We quantified effects of habitat on resident and transient coyote movement behavior using first passage time (FPT) analysis, which assesses recursive movement along an individual's movement path to delineate where they exhibit area-restricted search (ARS) behaviors relative to habitat attributes. We quantified monthly movement rates for 171 coyotes (76 residents and 53 transients) and then used estimated FPT values in generalized linear mixed models to quantify monthly habitat use for resident and transient coyotes. Transients had greater movement rates than residents across all months except January. Resident FPT values were positively correlated with agricultural land cover during fall and winter, but negatively correlated with agriculture during spring. Resident FPT values were also negatively correlated with developed habitats during May–August, deciduous land cover during June–August, and wetlands during September–January except November. FPT values of transient coyotes were positively correlated with developed areas throughout much of the year and near wetlands during July–September. Transient FPT values were negatively correlated with agriculture during all months except June and July. High FPT values (ARS behavior) of residents and transients were generally correlated with greater densities of edge habitat. Although we observed high individual variation in space use, our study found substantive differences in habitat use between residents and transients, providing further evidence that complexity and plasticity of coyote habitat use is influenced by territorial strategy.

54 ENVIRONMENTAL SCIENCES↗

Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models

Results from the fully and biogeochemically coupled simulations in which CO 2 increases at a rate of 1 % yr –1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO 2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm –1 ) and ocean (0.79 ± 0.07 PgC ppm –1 ), while the carbon–climate feedback over land (–45.1 ± 50.6 PgC °C –1 ) is about 3 times larger than over ocean (–17.2 ± 5.0 PgC °C –1 ). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 °C EgC –1 and is similar to that found in CMIP5 models (1.63 ± 0.48 °C EgC –1 ) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

54 ENVIRONMENTAL SCIENCES↗

Simplified aeroelastic modeling of horizontal axis wind turbines

Certain aspects of the aeroelastic modeling and behavior of the horizontal axis wind turbine (HAWT) are examined. Two simple three degree of freedom models are described in this report, and tools are developed which allow other simple models to be derived. The first simple model developed is an equivalent hinge model to study the flap-lag-torsion aeroelastic stability of an isolated rotor blade. The model includes nonlinear effects, preconing, and noncoincident elastic axis, center of gravity, and aerodynamic center. A stability study is presented which examines the influence of key parameters on aeroelastic stability. Next, two general tools are developed to study the aeroelastic stability and response of a teetering rotor coupled to a flexible tower. The first of these tools is an aeroelastic model of a two-bladed rotor on a general flexible support. The second general tool is a harmonic balance solution method for the resulting second order system with periodic coefficients. The second simple model developed is a rotor-tower model which serves to demonstrate the general tools. This model includes nacelle yawing, nacelle pitching, and rotor teetering. Transient response time histories are calculated and compared to a similar model in the literature. Agreement between the two is very good, especially considering how few harmonics are used. Finally, a stability study is presented which examines the effects of support stiffness and damping, inflow angle, and preconing.

Wendell, J. H.↗

Coupled SAM/Griffin Model of a Reference Pebble Bed High-Temperature Gas Cooled Reactor for Multi-Physics Simulations

An effort has been dedicated to developing a reference model for multi-physics coupled simulation of the pebble bed high-temperature gas-cooled reactor (PB-HTGR) with SAM and Griffin computer codes for safety analysis and licensing purpose. The reference problem of a prototypical reactor design serves as the foundation for the U.S. NRC (Nuclear Regulatory Commission) to verify the adequacy of computer codes and evaluation models for specific reactor types. In this work, a SAM model of the HTR-PM reactor has been developed based on publicly available design information and the multi-dimensional Pronghorn model developed by Idaho National Laboratory. The SAM HTR-PM model includes a multi-dimensional core region and 0- D/1-D fluid components. Additionally, a simplified air RCCS loop is modeled for decay heat removal. The Griffin model is based on a recent work by Idaho National Laboratory. The coupling between the models is achieved through the MOOSE MultiApp system. Both steady-state and transient scenarios were simulated to demonstrate the coupled model’s capability for multi-physics simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Liquid Core Detection and Strand Condition Monitoring in a Continuous Caster Using Optical Fiber

Real-time monitoring of the liquid core position during the continuous casting of steel has been demonstrated using low-cost distributed optical-fiber-based strain sensors. These sensors were installed on the containment roll support structures in the segments of a production continuous caster to detect the position of the solid–liquid interface and monitor the strand condition during the continuous casting. Distributed Fiber Bragg Grating sensors (FBGs) were used in this work to monitor strain at six roll positions in the caster. The sensor performance was first validated by comparing optical strain measurements with conventional strain gauge measurements in the lab. Next, optical strain measurements were performed on an isolated caster segment in a segment maintenance facility using hydraulic jacks to simulate the presence of a liquid core under the roll. Finally, the sensors were evaluated during caster operation. The sensors successfully detected the load increase associated with the presence of a liquid core under each instrumented roll location. Incidents of bulging and roll eccentricity were also detected using frequency analysis of the optical strain signal. The liquid core position measurements were compared using predictions from computer models (digital twins) in use at the production site. The measurements were in good agreement with the model predictions, with a few exceptions. Under certain transient caster operating conditions, such as spraying practice changes and SEN exchanges, the model predictions deviated slightly from the liquid core position determined from strain measurements.

47 OTHER INSTRUMENTATION↗

Lessons Learned from RELAP5-3D Modeling of HTTF

Modeling of HTTF using the RELAP5-3D ring model has shown an inability to predict transient temperature rise during HTTF experiments. This presentation provides an overview of the RELAP5-3D ring model, the HTTF modeling activities using the ring model, and reasons why the models failed to reproduce measured temperature rise. It also presents an overview of the path forward for RELAP5-3D modeling of HTTF to better predict the behavior of HTTF experiments

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Temporal Evolution of the Plasma Sheath Surrounding Solar Cells in Low Earth Orbit

High voltage solar array interactions with the space environment can have a significant impact on array performance and spacecraft charging. Over the past 10 years, data from the International Space Station has allowed for detailed observations of these interactions over long periods of time. Some of the surprising observations have been floating potential transients, which were not expected and are not reproduced by existing models. In order to understand the underlying processes producing these transients, the temporal evolution of the plasma sheath surrounding the solar cells in low Earth orbit is being investigated. This study includes lumped element modeling and particle-in-cell simulation methods. This presentation will focus on recent results from the on-going investigations.

Willis, Emily M.↗