Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “equations of fluid dynamics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Uncovering turbulent plasma dynamics via deep learning from partial observations

One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. Furthermore, this technique presents a novel paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Uncovering turbulent plasma dynamics via deep learning from partial observations

One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. This technique presents a novel paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Sneezing and asymptomatic virus transmission

The novel coronavirus disease (COVID-19) spread pattern continues to show that geographical barriers alone cannot contain a virus. Asymptomatic carriers play a critical role in the nature of this virus quickly escalating into a global pandemic. Asymptomatic carriers may transmit the virus unintentionally through sporadic sneezing. A novel Computational Fluid Dynamics (CFD) approach has been proposed with a realistic modeling of a human sneeze achieved by the combination of state-of-the-art experimental and numerical methods. This modeling approach may be suitable for future engineering analyses aimed at reshaping public spaces and common areas, with the main objective to accurately predict the spread of aerosol and droplets that may contain pathogens. This study shows that the biomechanics of a human sneeze, including complex muscle contractions and relaxations, can be accurately modeled by the angular head motion and the dynamic pressure response during sneezing. These have been considered as the human factors and were implemented in the CFD simulation by imposing a momentum source term to the coupled Eulerian–Lagrangian momentum equations. The momentum source was modeled by the measured dynamic pressure response in conjunction with the angular head motion. This approach eliminated the need to create an ad hoc set of inlet boundary conditions. With this proposed technique, it is easier to add multiple fixed and/or moving sources of sneezes in complex computational domains. Additionally, extensive sensitivity analyses based on different environmental conditions were performed, and their impact was described in terms of potential virus spread.

42 ENGINEERING↗
DOE OSTI2021

A comparative study of machine learning models for predicting the state of reactive mixing

Mixing phenomena are important mechanisms controlling flow, species transport, and reaction processes in fluids and porous media. Accurate predictions of reactive mixing are critical for many Earth and environmental science problems such as contaminant fate and remediation, macroalgae growth, and plankton biomass evolution. Here, to investigate the evolution of mixing dynamics under different scenarios (e.g., anisotropy, fluctuating velocity fields), a finite-element-based numerical model was built to solve the fast, irreversible bimolecular reaction-diffusion equations to simulate a range of reactive-mixing scenarios. A total of 2,315 simulations were performed using different sets of model input parameters comprising various spatial scales of vortex structures in the velocity field, time-scales associated with velocity oscillations, the perturbation parameter for the vortex-based velocity, anisotropic dispersion contrast (i.e., ratio of longitudinal-to-transverse dispersion), and molecular diffusion. The outputs comprised concentration profiles of reactants and products. The inputs to and outputs from these simulations were concatenated into feature and label matrices, respectively, to train 20 different machine learning (ML) models intended to emulate system behavior. These 20 ML emulators, based on linear methods, Bayesian methods, ensemble learning methods, and multilayer perceptrons (MLPs), were trained to classify the state of mixing and predict three quantities of interest (QoIs) characterizing species production, decay (i.e., average concentration, square of average concentration), and degree of mixing (i.e., variances of species concentration). Unsurprisingly, linear classifiers and regressors failed to reproduce the QoIs; however, ensemble methods (classifiers and regressors) and the MLP model accurately classified the state of reactive mixing and the QoIs. Among ensemble methods, random forest and decision-tree-based AdaBoost faithfully predicted the QoIs. At run time, trained ML emulators produced results times faster than the finite-element simulations. Due to their low computational expense and high accuracy, ensemble and MLP models are excellent emulators for these numerical simulations and great utilities in uncertainty quantification exercises, which can require 1,000s of forward model runs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

PDC Modifications for Analysis of Gas-Cooled Reactors with Direct Helium Brayton Cycle

Capabilities of the Plant Dynamics Code (PDC) have been extended to allow steady-state and transient simulation of graphite-cooled reactors with direct helium Brayton cycle. On the cycle side, the most significant code modification is the addition of helium properties, in the format required by the code’s equations. Since the code was already formulated to handle more complex real gas properties, adding helium as working fluid that behaves like ideal gas was fairly straightforward. A reactor module was added to PDC to simulate a reactor cooled by the working fluid of the Brayton cycle. Two options are supported: channel type, typical for graphite gas-cooled reactors, and pin type, typical for light-water and liquid metal-cooled reactors. The reactor module is an extension of the electrical heater model and simulates heat deposition in the fuel and transfer of this heat from the fuel to the coolant through the matrix and tube materials. The new reactor module becomes the third option in PDC for modeling heat addition to the cycle, besides previously modeled heat addition heat exchanger and electrical heater. In addition to those changes, other minor code modifications and improvements were introduced during the work of expanding PDC to modeling of gas-cooled reactors. These modifications are summarized in the last chapter of this report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Coupled Induction Machine and HVAC Models for Simulating HVAC Performance Considering Grid Dynamics in Buildings

This paper presents the development of novel models that integrate induction machines with HVAC equipment, such as pumps, heat pumps, and chillers, to analyze the impact of electrical parameters on the operational performance of thermo-fluid systems. The proposed model employs a coupling technique that captures the dynamic interactions between induction machines and HVAC systems. By integrating electrical, thermal, and mechanical dynamics, the models provide a comprehensive framework for simulating real-world scenarios, including interactions with the electrical grid. This achievement was made possible through the development of a Computationally Efficient and Accurate Induction Machine (CEAIM) model. Implemented using the equation-based Modelica language, the CEAIM model has been validated against experimental results, manufacturer data sheets, and various operating conditions. Its performance has been compared with existing induction machine models in the Modelica Standard Library (MSL), demonstrating superior accuracy and computational efficiency. The CEAIM model predicts torque, speed, and power consumption with a coefficient of determination (R 2 ) ranging from 0.98 to 1 and a coefficient of variation of root mean square error (CVRMSE) between 0.27% and 6.67%. Additionally, CEAIM scales more efficiently than conventional MSL models, with a slower computational growth rate in large-scale simulations. After thorough validation of the CEAIM model, it was coupled with HVAC equipment as this approach provides a detailed multi-dimensional view of capturing electrical transients and mechanical performance. To support this, a case study was conducted to showcase its capabilities.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Local correlation-based transition models for high-Reynolds-number wind-turbine airfoils

Abstract. Modern wind-turbine airfoil design requires robust performance predictions for varying thicknesses, shapes, and appropriate Reynolds numbers. The airfoils of current large offshore wind turbines operate with chord-based Reynolds numbers in the range of 3–15 million. Turbulence transition in the airfoil boundary layer is known to play an important role in the aerodynamics of these airfoils near the design operating point. While the lack of prediction of lift stall through Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) is well known, airfoil design using CFD requires the accurate prediction of the glide ratio (L/D) in the linear portion of the lift polar. The prediction of the drag bucket and the glide ratio is greatly affected by the choice of the transition model in RANS CFD of airfoils. We present the performance of two existing local correlation-based transition models – one-equation model (γ− SA) and two-equation model (γ-Reθt‾- SA) coupled with the Spalart–Allmaras (SA) RANS turbulence model – for offshore wind-turbine airfoils operating at a high Reynolds number. We compare the predictions of the two transition models with available experimental and CFD data in the literature in the Reynolds number range of 3–15 million including the AVATAR project measurements of the DU00-W-212 airfoil. Both transition models predict a larger L/D compared to fully turbulent results at all Reynolds numbers. The two models exhibit similar behavior at Reynolds numbers around 3 million. However, at higher Reynolds numbers, the one-equation model fails to predict the natural transition behavior due to early transition onset. The two-equation transition model predicts the aerodynamic coefficients for airfoils of various thickness at higher Reynolds numbers up to 15 million more accurately compared to the one-equation model. As a result, the two-equation model predictions are more comparable to the predictions from eN transition model. However, a limitation of this model is observed at very high Reynolds numbers of around 12–15 million where the predictions are very sensitive to the inflow turbulent intensity. The combination of the two-equation transition model coupled with the Spalart–Allmaras (SA) RANS turbulence model is a good method for performance prediction of modern wind-turbine airfoils using CFD.

17 WIND ENERGY↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Validation of Oregon State University High Temperature Test Facility Experiments Using Pronghorn

The OSU High Temperature Test Facility is a quarter-scale diameter, 1/64 scale volume test facility meant to replicate thermophysical phenomena in the prototypical General Atomics Modular High Temperature Gas Reactor. Tests pertaining to conduction cooldown events were performed from 2016-2019, providing a large database by which computational methods that are applicable to different length scales can be validated. One of these codes is Pronghorn, which is a coarse-meshed, porous-based subchannel thermal hydraulics code based on the MOOSE application, with the intention of better capturing the physics of both conduction and convection heat transfer within the OSU HTTF core. The goal of this summer project is to develop the framework by which Pronghorn can perform validation exercises of the HTTF core for benchmarking, by generating a mesh appropriate to the geometry of the HTTF core, developing input decks that accurately capture the initial and boundary conditions, materials, and relevant equations to the physics seen in the HTTF core, and using a postprocessor to compare simulation results to various experimental data. While validation of codes is a multi-year project, a mesh has been generated and tested in Pronghorn that meets mass conservation and basic heat transfer principles. The next step is to accurate depict the fluid inlet and outlet boundary conditions, which will be performed using computational fluid dynamics software.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2021

Parameter space for magnetization effects in high-energy-density plasmas

Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Demonstration and performance testing of extreme-resolution simulations with static meshes on Summit (CPU & GPU) for a parked-turbine configuration and an actuator-line (mid-fidelity model) wind farm configuration (ECP-Q4 FY2020 Milestone Report)

The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources. The primary physics codes in the ExaWind simulation environment are Nalu-Wind, an unstructured-grid solver for the acoustically incompressible Navier-Stokes equations, AMR-Wind, a block-structured-grid solver with adaptive mesh refinement capabilities, and OpenFAST, a wind-turbine structural dynamics solver. The Nalu-Wind model consists of the mass-continuity Poisson-type equation for pressure and Helmholtz-type equations for transport of momentum and other scalars. For such modeling approaches, simulation times are dominated by linear-system setup and solution for the continuity and momentum systems. For the ExaWind challenge problem, the moving meshes greatly affect overall solver costs as reinitialization of matrices and recomputation of preconditioners is required at every time step. The choice of overset-mesh methodology to model the moving and non-moving parts of the computational domain introduces constraint equations in the elliptic pressure-Poisson solver. The presence of constraints greatly affects the performance of algebraic multigrid preconditioners.

17 WIND ENERGY↗
DOE OSTI2025

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗