Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient problem”

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 217 records · Page 12

Nuclear fuels for transient test reactors

Transient test reactors with the ability to test fissile specimens under extreme conditions have been crucial tools in the development of nuclear technologies. Less than 10 unique facility designs have ever been constructed, most of which remain operational today and still use the original nuclear fuel constructed for them more than 40 years ago. Historic fuel systems for transient test reactors vary in significant ways which have marked influences on reactor capabilities. Eventually, new fuel will be needed to support the longevity of transient test reactor missions. This paper reviews precedent transient reactor fuel systems in the context of their unique requirements. A few key conclusions are illustrated by comparing and contrasting these transient test reactors. Fuel composites which are mostly graphite can enable transient reactors with very high neutron fluence capability (>2E16 n/cm 2 ) and are amenable to longer “shaped” transients but cannot achieve pulses <10 ms in duration. Reducing the graphite-to-uranium ratio can yield a very narrow pulse capability but delivers less fluence and requires cores with considerably more fissile material. Designs based on uranium dioxide (UO 2 ) make use of readily available materials to create compact cores with narrow pulse width capabilities but with moderate neutron fluence capabilities (~2E15 n/cm 2 ). Uranium zirconium hydride (U-ZrHx) is a well-established fuel system for pulsing reactors which has been intermittently manufactured throughout the decades. U-ZrHx offers similar capabilities to UO 2 designs in terms of nuclear kinetics, but with about half the fluence capability (~1E15 n/cm 2 ). An evolution of the UO 2 system, termed “ternary ceramic” fuel, shows that dispersing UO 2 in zirconium oxide and calcium oxide can increase fluence capability greatly (~7E15 n/cm 2 ), but is not presently a commonly available fuel form. A unique composite of UO 2 and beryllium oxide (UO 2 -BeO) can be used to create a core with similar kinetics and compact core geometry as U-ZrHx designs, but with significantly higher fluence capability (~6E15 n/cm 2 ). Like ternary ceramic fuel, newly fabricated UO 2 -BeO would require reestablishing its historic manufacturing process which would be further complicated by the health hazards associated with beryllium. In conclusion, like most engineering problems, there is no perfect solution, but this paper outlines the advantages and disadvantages of candidate fuel options to help guide detailed evaluations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Growth-induced Donnan exclusion influences swelling kinetics in highly charged dynamic polymerization hydrogels

Polymeric gels crosslinked by DNA sequences can exploit DNA strand-displacement reactions to promote swelling through dynamic polymerization. The degree of swelling and the rate of swelling must be directly tunable to achieve the promise of programmable soft matter. Though the kinetics of the strand-displacement reaction provide insertion rates up to 10 4 /Molar/second as measured in bulk solution, DNA hydrogel swelling can take upwards of 30 h to complete. Computational modeling of the reaction-induced swelling of these gels with our recently-developed reactive electrochemomechanical theory (Zimmerman et al., 2024) suggests that their extraordinarily slow swelling is partly due to a scaling mismatch between the addition of charge and the addition of fluid volume, leading to a large transient increase in the fixed charge density. The significant increase in the gel’s fixed charge density, due to the binding of negatively charged DNA, sharply restricts the concentration of mobile hairpins through the phenomenon of Donnan charge exclusion, an effect commonly exploited in nanofiltration applications using polymeric membranes. The scaling problem is overcome when the mean additional swelling provided to the hydrogel by addition of a crosslink is above a critical value, thus the swelling outpaces the charge accumulation, leading the fixed charge density to drop and significantly accelerating the swelling process. This study shows that Donnan exclusion can explain the kinetics of DNA hydrogel swelling, and studies ways to modulate the reaction speed by either modifying the salt concentration or increasing or decreasing the number of base pairs in each DNA sequence.

42 ENGINEERING↗

Feasibility Study on Implementing a Staggered-Grid Finite Volume Method for System Analysis Code Development Under the MOOSE Framework

Here, this work summarizes a feasibility study on testing numerical algorithms that are suitable and efficient for advanced system analysis code development under the mutli-physics framework, MOOSE. The key to the test bed is the implementation of high-order one-dimensional staggered-grid finite volume method (SG-FVM), and its direct interaction with the linear/nonlinear solver, PETSc. The test bed utilized a more flexible code structure to enable the finite volume method implementation and direct interacting with the solver package, instead of using the natively supported finite element method by the framework. Using a suite of selected test problems with different problem sizes and levels of complexity, the implemented SG-FVM demonstrated superior performance improvement against a direct finite element method implementation through MOOSE. On two computer systems, the speedup was observed to be significant, with at least one order of magnitude of solving time reduction. For a complex reactor model, transient simulation was performed using the newly developed finite volume method code, the results of which agree very well with the reference results from the finite element method code. Overall, this study demonstrates a successful feasibility study on the proposed numerical algorithms and software structure to support advanced system analysis tool development.

MOOSE↗

Nonequilibrium dynamics of two-level systems directly after cryogenic alternating bias

Two-level systems (TLSs) are a dominant decoherence mechanism in superconducting qubits, microwave kinetic inductance detectors, and quantum dots. TLSs are tunneling states commonly found in amorphous materials and interfaces, which interact electromechanically, leading to energy loss. Fluctuations in the frequency of the TLSs are a significant problem for qubit operation and force recalibration every few hours. In this study, we probe the effect of alternating bias at cryogenic temperatures on TLSs in a purpose-built LC oscillator. When an in situ alternating bias is applied, the steady-state spectra of the TLSs disappear. Spectroscopy reveals transient behavior, in which the TLS frequency fluctuates on the order of minutes. Thermal cycling above 10 K reverses these effects, restoring the TLS spectrum to its original state. Importantly, the intrinsic loss tangent of the LC oscillator remains unchanged before and after the application of the alternating bias. We propose that the disappearance of the steady-state spectrum is caused by nonequilibrium energy buildup from strain in the oxide film introduced by the pulsed voltage-bias sequence. Understanding this nonequilibrium energy could shed light on TLS fluctuations and reduce the calibration requirements for qubits.

general physics↗

Demonstration of CFD to support the coupled analysis of a reactor pressure vessel subjected to pressurized thermal shock

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended duration or under accident condition. As one of the most significant components of a reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Aging degradation mechanisms (e.g. thermo-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The prediction of the temperature field is achieved by using computational fluid dynamics (CFD) simulation. In this report, a demonstration CFD standalone simulation is performed to support coupled analysis for Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS). The study use a simplified computational domain to represents a real RPV. The purpose of the study is to demonstrate the transient temperature response of RPV to ECC injection. The CFD model is built in a robust and efficient way for further coupled calculation. The next steps of this work, including the coupled thermal and tensor mechanics capabilities using Cardinal are expected to be complete by the end of FY21 for the demo problem. After this, into FY22, the capability will be demonstrated for a realistic RPV.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DECOVALEX-2019 (Task D: INBEB Final Report)

The DECOVALEX Project is an on-going international research collaboration, established in 1992, to advance the understanding and modeling of coupled Thermal (T), Hydrological (H), Mechanical (M) and Chemical (C) processes in geological in geological systems. This document is the final report of Task D which was proposed by the DECOVALEX Project and coordinated by CIMNE_UPC, presenting the technical definitions of the problems studied, approaches applied, achievements made and outstanding issues for future research. Task D of the DECOVALEX 2019 project is devoted to the study of the hydro-mechanical (HM) and thermo-hydro-mechanical (THM) Interactions in Bentonite Engineered Barriers (INBEB). Special attention is paid to the evolution of barrier heterogeneity under transient conditions. A proper understanding of the processes occurring during the transient period requires the appropriate simulation of the behavior of the engineered barrier by means of coupled numerical analyses. In this context, the main objective of the Task is to assess the capabilities of numerical formulations and codes to simulate and interpret the HM and THM behavior of the bentonite barrier, with a particular focus on the final state reached on saturation.

58 GEOSCIENCES↗

Transient Efficiency Flexibility and Reliability Optimization of Coal-Fired Power Plants: Model-Predictive Control Library Development (Report)

This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.

01 COAL, LIGNITE, AND PEAT↗

Development of Explainable, Knowledge-Guided AI Models to Enhance the E3SM Land Model Development and Uncertainty Quantification

Focal Area(s): (2)Predictive modeling using AI techniques and AI-derived model components; use of AI and other tools to design a prediction system comprising of a hierarchy of models. (3) Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge- guided AI. Science Challenge: The Energy Exascale Earth System Model (E3SM) is a fully coupled, state-of-the-science Earth system model that uses code optimized for DOE's advanced computers to address the most critical scientific questions facing our nation and society (Golaz et al., 2019). The E3SM Land model (ELM) is designed to understand how the changes in terrestrial land surfaces will interact with other Earth system components and has been used to understand hydrologic cycles, biogeophysics, and ecosystem dynamics. In spite of great successes, the ELM has several known issues that restrain rapid improvements. For example, the ELM uses equilibrium models to simulate dynamic land-climate interactions and it requires long model spin-up time to identify suitable initial conditions for transient simulations. The ELM lacks built-in uncertainty mechanisms that can improve the robustness of model predictions. The ELM is a holistic, deterministic model system with a rigid design, and in many situations, it is hard to modify the ELM system to incorporate new theory/hypothesis and new data across scales to address emerging science problems (such as predicting the impacts of water cycle extremes). In addition, The ELM is technically optimized for traditional CPU-centric computers and it cannot fully utilize the current and incoming leadership computers for model simulations and uncertainty quantification (UQ). The success of artificial intelligence (AI) has inspired scientists to use AI models to discover intrinsic features from simulation data (Chattopadhyay et al., 2020) and observational data (Reichstein et al., 2019) to gain further process understanding of Earth science problems. However, autonomous AI model training through deep learning usually requires a huge amount of annotated data. To overcome the limitations from the data and computing resources, knowledge-guided AI models are necessary where human-knowledge is ingested in model construction (Banino et al., 2018) and training process (Silver et al., 2016) for efficient learning. Herein, we present a new way that leverages the process understanding from the ELM to guide AI model development for the ELM enhancement and UQ. We hope this study can inspire further Earth and environmental system model developments and transformations.

54 ENVIRONMENTAL SCIENCES↗

PETSc DMNetwork: A Library for Scalable Network PDE-Based Multiphysics Simulations

We present DMNetwork, a high-level package included in the PETSc library for the simulation of multiphysics phenomena over large-scale networked systems. The library aims at applications that have networked structures such as those in electrical, gas, and water distribution systems. DMNetwork provides data and topology management, parallelization for multiphysics systems over a network, and hierarchical and composable solvers to exploit the problem structure. DMNetwork eases the simulation development cycle by providing the necessary infrastructure through simple abstractions to define and query the network components. This article presents the design of DMNetwork, illustrates its user interface, and demonstrates its ability to solve multiphysics systems, such as an electric circuit, a network of power grid and water subnetworks, and transient hydraulic systems over large networks with more than 2 billion variables on extreme-scale computers using up to 30,000 processors.

extreme-scale computing↗

Analytical Large-Signal Modeling of Inverter-Based Microgrids With Koopman Operator Theory for Autonomous Control

The microgrid (MG) plays a crucial role in the energy transition, but its nonlinearity presents a significant challenge for large-signal power systems studies in the electromagnetic transient (EMT) time scale. In this paper, we develop a large-signal linear MG model that considers the detailed dynamics of the primary and zero-control levels based on the Koopman operator (KO) theory. Firstly, a set of observable functions is carefully designed to capture the nonlinear dynamics of the MG. The corresponding linear KO is then analytically derived based on these observables, resulting in the linear representation of the original nonlinear MG with observables as the new coordinate. The influence of external input on the system dynamics is also considered during the derivation, enabling control of the MG. We solve the voltage control problem using the traditional linear quadratic integrator (LQI) method to demonstrate that textbook linear control techniques can accurately control the original nonlinear MG via the developed KO-linearized MG model. Our proposed KO linearization method is generic and can be easily extended for different control objectives and MG structures using our analytical derivation procedure. Here, we validate the effectiveness of our methodology through various case studies.

25 ENERGY STORAGE↗

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗

Multiphysics Analysis of Load Following and Safety Transients for MicroReactors

The tools developed within the U.S. Department of Energy (DOE) Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aim at providing high fidelity multiphysics modeling capabilities to support design and licensing of various types of advanced nuclear reactors, including the technologies being developed by U.S. microreactor vendors relying on heat pipe and gas-cooled technologies. In FY-2022, the NEAMS Multiphysics Applications team made significant progress both in demonstrating capabilities applied to microreactor problems, and in supporting NEAMS developers.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Improvement and Verification of Online Cross Section Generation Capability of Griffin for TRISO-fueled Reactors

Griffin, a MOOSE-based reactor multiphysics code jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE Office of Nuclear Energy’s NEAMS program, has pursued the development of an online multigroup cross section generation capability for a few years to enable high-fidelity, problem-dependent neutronics analyses of advanced thermal reactors. Recent advancements in Griffin’s online multigroup cross section generation capability have significantly improved the accuracy, robustness, and efficiency of self-shielding calculations for both prismatic and pebble-bed TRISO-fueled reactor applications. Key developments include a unified fuel self-shielding method applicable to both TRISO and annular compact/spherical shell fuel zone geometries; an advanced Dancoff Category-based Equivalence Theory using a bell function for non-fuel resonance treatment, achieving more than an order-of-magnitude speedup compared to the Tone method; an on-the-fly multigroup equivalence approach to mitigate group condensation errors; and a streaming correction method for pebble-bed homogenization. A proof-of-concept demonstration of on-the-fly group condensation with consistent P0 transport correction was also achieved. The method reproduced direct fine-group solutions with excellent accuracy (eigenvalue errors within 10 pcm and pin-power differences within 0.5%), but due to performance limitations of the current fixed-source solver, improvements to solver efficiency will be addressed in future work. Verification tests were performed on graphite-moderated TRISO-fueled two-dimensional core benchmark problems representing gas-cooled microreactors, heat pipe-cooled microreactors, gas-cooled pebble-bed reactors, and fluoride salt-cooled high-temperature reactors. Across all cases, Griffin showed excellent agreement with Serpent2 continuous energy Monte Carlo solutions: eigenvalue errors within 200 pcm, pin-power root-mean-square errors within 2%, and control rod and drum worth errors less than 2%. It should be noted that, for the benchmark problem, cross section generation contributed less than 3% of the total simulation times. These results demonstrate that Griffin’s online cross section generation capability delivers accurate and efficient reactor physics solutions across a wide spectrum of TRISO-fueled advanced reactor designs. With further improvements to the fine-group fixed-source solver and planned extensions to depletion, transients, and coupled neutron–gamma transport, Griffin will be well-positioned to become a powerful and comprehensive tool for advanced reactor analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

COUPLED GRIFFIN AND PRONGHORN SIMULATION OF THE MOLTEN SALT FAST REACTOR (MSFR) FOR THE VIRTUAL TEST BED

The Virtual Test Bed (VTB) repository hosts a wide range of challenge problems to showcase modeling and simulation capabilities to support advanced reactor demonstration. An overview of a coupled multiphysics model for Molten Salt Reactors (MSR) is presented here. The analysis leverages MOOSE-based tools (notably Griffin and Pronghorn) for neutronic and thermal hydraulic simulation. Neutron precursor drifting capability is showcased, along with some initial turbulence models. Both steady-state and transient coupled multiphysics results for an MSR concept are discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis of Molten Salt Thermal Energy Storage Tank Sizing and Preheating

Molten Salt Thermal energy storage (TES) is essential for concentrating solar power (CSP) plants, enhancing their capacity factor and dispatchability to ultimately reduce electricity costs. This technology is becoming more important for other technologies including nuclear power, thermal-electrical storage, and concentrating solar thermal (CST). While molten salt TES tanks have been successfully deployed in lower-temperature parabolic-trough plants, the current state-of-the-art central receiver systems operating at temperatures up to 565 degrees Celsius have experienced costly and frequent failures. These failures are associated with the immaturity of the technology and problems in tank design, fabrication, commissioning and, sometimes, aggressive operation. A primary contributing factor is the absence of a dedicated technical standard for high-temperature molten nitrate TES tanks; designers must currently rely on existing codes like API 650 and ASME Section II, which are insufficient because they fail to account for the unique demands of high temperatures, thermal cycling, and transient conditions.

25 ENERGY STORAGE↗

Code Benchmark of the HTTF Pressurized Conduction Cooldown Test Using SAM

The High Temperature Test Facility (HTTF) at Oregon State University is an integral system test facility to simulate postulated reactor transients of prismatic high-temperature gas-cooled reactors(HTGRs). A series of test campaigns was launched, providing abundant test data that could be used to benchmark reactor system analysis codes like the System Analysis Module (SAM). In this study, a SAM model of the facility is developed based on the two-dimensional (2D) ring model approach. All components including the ceramic matrix, graphite heaters, helium coolant channels, core barrel, upcomer, pressure vessel, and reactor cavity cooling system are modeled as concentric cylindrical rings. The model is used to simulate one of the benchmark problems-Pressurized Conduction Cooldown (PCC)-within the scope of the Organisation for Economic Co-operation and Development Nuclear Energy Agency International HTTF Benchmark. The simulations consist of two parts. In the first part, operating and boundary conditions as well as thermophysical properties of materials are specified for the benchmark problem. In this work, results from the first part will be used in code-to-code comparison. In the second part, the SAM model is used to simulate Test PG-27, which is the first PCC test carried out in the HTTF, with only two of the ten heater banks activated. The results in the second part are used for code-to-data comparison. Because the helium coolant flow rate is not measured in this facility, it is estimated using the input power and inlet/outlet coolant temperatures. Additionally, radial heat flow in the ceramic blocks is complicated by hundreds of cylindrical coolant channels and heater rods embedded in them. As such, it is necessary to deduce an effective thermal conductivity for the ceramic to analyze the core thermal behavior. SAM predictions of the helium coolant and ceramic temperatures are compared with test data measured in three equivalent sectors. Overall, the SAM results agree reasonably well with test data within the variation of data among the three sectors, which demonstrates SAM's capability in capturing transient effects in HTGR using the simplified 2D ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Multiobjective Dynamic Optimization for Optimal Load-Following of Natural Gas Combined Cycle Power Plants under Stress Constraints

With the increasing penetration of intermittent renewable energy sources into the electric grid, there is an associated need for conventional thermal power plants that are designed to operate at base-load conditions to cycle their load rapidly and frequently and operate under low-load conditions. Fast load-following operations lead to plant efficiency loss and the reduction of equipment life. A high-fidelity dynamic model of a natural gas combined cycle power plant, with rigorous equipment level submodels, is developed to capture the plant transient performance and the thermomechanical stress evolution at the high-pressure drum is computed to assess the drum life consumption. Stress evolution is modeled for the location of the edge at the drum/downcomer junction that experiences high circumferential stress amplitude during the fast load-following operation, leading to higher fatigue damage than locations considered under existing design standards. A dynamic optimization problem is solved for optimal load-following operation. Depending on the value of the stress constraint and the desired ramp rate, satisfying the desired stress constraint may be infeasible without relaxing the ramp rate. A multiobjective dynamic optimization problem is solved using a lexicographic approach that minimizes ramp rate relaxation and maximizes efficiency while satisfying stress constraints, avoiding spraying to saturation, and maintaining main steam and reheat steam temperature within bounds. It was observed that the optimal ramp rate can be highly nonlinear as opposed to the industry-standard linear ramp rate. Nonlinear optimal ramp-rate profiles not only help to avoid stress constraints that may be unavoidable by using the linear profile but also result in higher efficiency than the linear profile. Furthermore, the study shows that there is a strong tradeoff between the relaxation in the plant ramp rate and the time-average thermal efficiency.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗