Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational fluid dynamics simulation”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Blood flow imaging by optimal matching of computational fluid dynamics to 4D-flow data

Three-dimensional, time-resolved blood flow measurement (4D-flow) is a powerful research and clinical tool, but improved resolution and scan times are needed. Therefore, this study aims to (1) present a postprocessing framework for optimization-driven simulation-based flow imaging, called 4D-flow High-resolution Imaging with a priori Knowledge Incorporating the Navier-Stokes equations and the discontinuous Galerkin method (4D-flow HIKING), (2) investigate the framework in synthetic tests, (3) perform phantom validation using laser particle imaging velocimetry, and (4) demonstrate the use of the framework in vivo. An optimizing computational fluid dynamics solver including adjoint-based optimization was developed to fit computational fluid dynamics solutions to 4D-flow data. Synthetic tests were performed in 2D, and phantom validation was performed with pulsatile flow. Reference velocity data were acquired using particle imaging velocimetry, and 4D-flow data were acquired at 1.5 T. In vivo testing was performed on intracranial arteries in a healthy volunteer at 7 T, with 2D flow as the reference. Results Synthetic tests showed low error (0.4%-0.7%). Phantom validation showed improved agreement with laser particle imaging velocimetry compared with input 4D-flow in the horizontal (mean -0.05 vs -1.11 cm/s, P < .001; SD 1.86 vs 4.26 cm/s, P < .001) and vertical directions (mean 0.05 vs -0.04 cm/s, P = .29; SD 1.36 vs 3.95 cm/s, P < .001). In vivo data show a reduction in flow rate error from 14% to 3.5%. Phantom and in vivo results from 4D-flow HIKING show promise for future applications with higher resolution, shorter scan times, and accurate quantification of physiological parameters.

4D-flow MRI↗
DOE OSTI2021

Computational Characterization and Boundary Condition Models of a Jet Interaction Fluidic Oscillator

Fluidically oscillating jet actuators that generate sweeping jets when supplied with a pressurized fluid have been used to mitigate separation and reduce drag in a range of flow control applications. Their implementation in future flight vehicles would require fundamental understanding and accurate predictive techniques of the physics of their internal flow and jet formation. The present investigations focus on high-fidelity, time-accurate simulations to characterize the flow physics of the actuator in quiescent conditions. An important element of the present simulations is to demonstrate the ability of the computational fluid dynamics (CFD) solver to predict the jet characteristics and provide a basis for the development of improved boundary conditions (BC) without entirely resolving the geometrical features of the fluidic device. The CFD-predicted oscillation frequencies of the engendered jets were found to be in excellent agreement with experiments, even on two-dimensional meshes. The study revealed that three-dimensional simulations are required to capture some of the flow features of the sweeping jet such as the double peak in time-averaged velocity distributions downstream of the actuator’s orifice that were measured in experiments. Furthermore, several approaches for modeling the actuator were implemented and assessed in quiescent conditions. The evaluation of a boundary condition at the device throat, based on the phase-averaged flow variables, provides the basis for devising surface-based boundary conditions. The influence and necessity of including turbulent characteristics as part of the boundary conditions have also been identified.

42 ENGINEERING↗
DOE OSTI2024

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗
DOE OSTI2021

Simulations of the shock-driven Kelvin–Helmholtz instability in inclined gas curtains

In this paper, we present simulation results for the two-dimensional, shock-driven Kelvin–Helmholtz instability. Simulations are performed with a Mach 2.0 shock propagating through a finite-thickness curtain of gas inclined at an angle α0=30° with respect to the shock plane. After the passage of the shock, the gas curtain is accelerated along its axis. A perturbation develops due to shock reflection near the lower wall, and a Kelvin–Helmholtz instability forms near the vertical center of the curtain. This is the first known numerical reproduction of these phenomena that have previously been observed in experiments with an inclined cylindrical gas column. The effects of varying Mach number and column width were explored in detail to complement experimental data. Additionally, the dependence of the Kelvin–Helmholtz wavelength on Mach number closely matches the relationship observed in experiments. This supports the notion that the observed instability is effectively two-dimensional and inviscid (like classical Kelvin–Helmholtz). The growth rate of the perturbations in the gas curtain was also found to be similar for different Mach numbers. The perturbation at the curtain foot, previously unreported in experiments, was found to have a similar relationship to Mach number as the Kelvin–Helmholtz instability. Both perturbation wavelengths are found to be proportional to layer width. Simulations were performed with the fast interfaces and transport in the atmosphere, an exascale ready, graphics processing unit-accelerated compressible flow solver developed at the University of New Mexico.

42 ENGINEERING↗
DOE OSTI2023

Coupled Monte Carlo Transport and Conjugate Heat Transfer for Wire-Wrapped Bundles Within the MOOSE Framework

Cardinal is an open-source application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics (CFD) to the Multiphysics Object-Oriented Simulation Environment (MOOSE), closing neutronics and thermal-fluid gaps in conducting high-resolution multiscale and multiphysics analyses of nuclear systems. We first provide a brief introduction to Cardinal's software design, data mapping, and coupling strategy to highlight our approach to overcoming common challenges in high-fidelity multiphysics simulations. Here we then present two Cardinal simulations for hexagonal pin bundles. The first is a validation of Cardinal's conjugate heat transfer coupling of NekRS's Reynolds-Averaged Navier Stokes model with MOOSE's heat conduction physics for a bare seven-pin Freon-12 bundle flow experiment. Predictions for pin surface temperatures under three different heating modes agree reasonably well with experimental data and similar CFD modeling from the literature. The second simulation is a multiphysics coupling of OpenMC, NekRS, and BISON for a reduced-scale, seven-pin wire-wrapped version of an Advanced Burner Reactor bundle. Wire wraps are approximated using a momentum source model, and coupled predictions are provided for velocity, temperature, and power distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Subject-specific modeling framework for particle deposition using computational fluid dynamics

Quantifying particle deposition and dose in the respiratory tract requires a physiologically realistic representation and reproducible computational workflows. However, existing modeling frameworks, such as the International Commission on Radiological Protection (ICRP) compartmental models and the Multiple Path Particle Dosimetry (MPPD) tool, lack detailed deposition profiles and subject-specific capabilities. The combination of advances in computer vision algorithms applied to the respiratory tract and Computational Fluid and Particle Dynamics (CFPD) allows high-fidelity simulations of particle behavior in anatomically accurate geometries derived from individual CT scans. The segmentation, preprocessing, and file preparation task for a CFPD simulation was often time-consuming, and no prior studies to-date have yet presented a fully automated framework. This work presents a fully automated workflow to obtain individualized particle deposition profiles in the human respiratory tract. The pipeline starts with segmenting upper and lower airway geometries using morphological and deep learning-based methods, generating three-dimensional (3D) models from CT imaging data. Next, a series of algorithms are presented to quality check and prepare the 3D geometry for a CFD or CFPD simulation. The preprocessing step includes correcting geometric artifacts, enforcing a physically consistent mesh, and automatically identifying and capping multiple outlets, which is required for CFD/CFPD simulations. These processed models are then input into open-source (OpenFOAM) or commercial (StarCCM+) CFD solvers, where flow and transient particle transport equations — including turbulence and particle–wall interactions are solved under realistic breathing conditions. Finally, the resulting particle deposition profiles can be integrated with Monte Carlo radiation transport codes and state-of-the-art computational phantoms to assess organ-specific absorbed doses in scenarios of radioactive aerosol inhalation. The presented work streamlines respiratory tract segmentation, preprocessing for CFD/CFPD simulations, and integration with dose assessment workflows, reducing manual intervention and improving access to high-fidelity, subject-specific modeling. The high precision in predicted particle deposition and dose distributions can improve personalized treatment strategies in respiratory medicine and refine dose estimates for radiation protection.

AI↗
DOE OSTI2023

Near-continuum, hypersonic oxygen flow over a double cone simulated by direct simulation Monte Carlo informed from quantum chemistry

A large-scale, fully resolved direct simulation Monte Carlo (DSMC) computation of a non-equilibrium, reactive flow of pure oxygen over a double cone is presented. Under the simulated near-continuum conditions, the computational demands are shown to be significant because of the wide range of length scales that must be resolved. Therefore, robust grid adaption capabilities and efficient parallelization of the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code that is utilized in this work are essential. The thermochemical and transport collision models were selected for efficiency and simplicity. First-principles data, obtained from the highly accurate direct molecular simulation method, were used to inform the collision models’ parameters. Importantly, because SPARTA implements molecular collision models using collision-specific energies, the resulting macroscopic relaxation rates were evaluated a posteriori via zero-dimensional heat bath simulations. The comparisons of surface properties, namely heat flux and pressure, show very close agreement with previous computational fluid dynamics (CFD) results. Differences with the measurements were found to be similar to the CFD simulations. The unresolved discrepancy with the measurements could be due to inconsistent free stream conditions with the actual experimental data or missing physical phenomena altogether, for example atomic and molecular oxygen electronically excited states, three-dimensional effects, or more complex gas–surface interactions. As shown in this work, the advantages of obtaining a DSMC particle solution for these flows reside in the method's ability to be directly informed from first principles and to seamlessly describe internal energy non-equilibrium for all modes. With the advent of exascale computing and beyond, particle methods will be an increasingly important tool to verify the validity of physical assumptions in reduced-order models via fully resolved, experimental-scale simulations, down to the level of molecular-level distributions.

Mechanics↗
DOE OSTI2025

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

CFD modeling of non-catalytic, partial-oxidation engine reformer for flare mitigation

Flaring associated natural gas is commonly employed in the oil and gas industry to reduce methane (CH 4 ) emissions but generates carbon dioxide (CO 2 ) and harmful pollutants, significantly contributing to air pollution and posing risks to public health. To mitigate this impact, M2X Energy Inc. has developed a small-scale, modular gas-to-methanol system. This system features an engine reformer that performs fuel-rich partial oxidation of wellhead gas to produce syngas—a mixture of carbon monoxide (CO) and hydrogen (H 2 )—followed by a downstream reactor for methanol synthesis. This study focused on computational fluid dynamics (CFD) modeling of the engine reformer to simulate partial oxidation chemistry, predict the rich-burn operating limit, and assess syngas quality, ultimately aiding in design and operational optimization. The CFD model, developed within a Reynolds-Averaged Navier-Stokes (RANS) turbulence framework, incorporated sub-models for turbulent combustion, a chemical mechanism with polycyclic aromatic hydrocarbon (PAH) pathways, and soot emissions to accurately capture the fuel-rich, turbulent jet ignition and combustion processes. Model validation against experimental data showed good agreement across pre- and main-chamber pressures, apparent heat release rates, and exhaust gas concentrations of key species (H 2 , CO, CO 2 , CH 4 ) for varying intake equivalence ratios. Here, the model identified a rich-burn operating limit near a fuel-air equivalence ratio of 2.35, consistent with experimental observations. Furthermore, syngas quality analysis revealed that extending the rich-burn limit through engine reformer optimization could enhance syngas production, contributing to higher methanol synthesis efficiency.

Computational Fluid Dynamics↗
DOE OSTI2023

Predicting Wind Loading and Instability in Solar Tracking PV Arrays

Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.

aerodynamics↗
DOE OSTI2025

CFD Simulation of Aerobic Gas Fermentation to Enable Commercial Conversion of CO 2 into Aquaculture and Animal Feed: Cooperative Research and Development (Final Report)

NovoNutrients’ fermentation technology uses energy from hydrogen to transform industrial CO2 emissions into premium animal feed ingredients and other valuable products. A single NovoNutrients’ commercial manufacturing plant will capture and convert over 200,000 tons/yr of CO2 into over 100,000 tons/yr of high-protein feed. Key to the rapid and widespread deployment of the technology is maximization of its productivity and energy efficiency. Robust, physically based computational models of the technology will significantly increase productivity and efficiency, accelerating NovoNutrients’ technology to manufacturing scale. NREL has unique capabilities for creating and running such computational models. NREL's existing aerobic bioreaction computational fluid dynamics (CFD) models will be adapted to NovoNutrients’ gas fermentation (CO2, H2, O2) technology. The multiphysics CFD simulations require thousands of high-performance computing (HPC) node hours to simulate the complex geometries and contents of NovoNutrients’ industrial bioreactors. The experimentally validated CFD models were used to identify optimally efficient and productive bioreactor designs and operating conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modeling and simulation of gas separations with spiral–wound membranes

Models for gas separations with spiral–wound membranes are developed and found to exhibit good agreement with experiments performed on N 2 /O 2 mixtures. The two–dimensional (2D) model can be accurately approximated by a one–dimensional (1D) surrogate model when the spacer widths are chosen to make the channel pressure drops small. Subsequently, the separation of propane/propylene mixtures from the recycle purge stream of a polypropylene reactor is investigated. Assuming ideal gas is found to lead to significant overestimations in membrane stage cuts (sometimes more than 10%), an extent comparable to that associated with extrapolating constant olefin permeance from a low–pressure condition. While olefin permeance can change significantly with pressure, using a constant–permeance formulation can result in a small (< 2.5%) underprediction in stage cut if the value for the permeance is taken from the feed condition. Finally, membrane properties and costs necessary for a viable separation process are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Preparing an Incompressible-Flow Fluid Dynamics Code for Exascale-Class Wind Energy Simulations

The U.S. Department of Energy has identified exascale-class wind farm simulation as critical to wind energy scientific discovery. A primary objective of the ExaWind project is to build high-performance, predictive computational fluid dynamics (CFD) tools that satisfy these modeling needs. GPU accelerators will serve as the computational thoroughbreds of next-generation, exascale-class supercomputers. Here, we report on our efforts in preparing the ExaWind unstructured mesh solver, Nalu-Wind, for exascale-class machines. For computing at this scale, a simple port of the incompressible-flow algorithms to GPUs is insufficient. To achieve high performance, one needs novel algorithms that are application aware, memory efficient, and optimized for the latest-generation GPU devices. The result of our efforts are unstructured-mesh simulations of wind turbines that can effectively leverage thousands of GPUs. In particular, we demonstrate a first-of-its-kind, incompressible-flow simulation using Algebraic Multigrid solvers that strong scales to more than 4000 GPUs on the Summit supercomputer.

algebraic multigrid↗
DOE OSTI2025

CFD modeling of natural circulation in LiCl-KCl molten salt closed loop

Characterizing flow within a molten salt closed-loop system is crucial for assessing system requirements, evaluating performance, and identifying potential flaws. Direct flow measurement using instrumentation is challenging due to extreme environmental conditions and the limitations associated with measuring molten salt flow under natural convection. Here, this study aims to provide comprehensive insights into the thermal-hydraulic behavior of a closed loop, with a particular focus on temperature distribution and velocity prediction. The Computational Fluid Dynamics (CFD) model demonstrated the capability to effectively simulate and predict both temperature distributions and flow velocities within the molten salt loop. The CFD model's predictive capability was validated by its ability to replicate temperature measurements under varying boundary conditions. The analysis revealed that the CFD model tends to underpredict temperatures in the cold leg and overpredict them in the hot leg, highlighting the need for continuous model refinement and acknowledging the limitations of using a steady-state approach. Furthermore, the potential of using external temperature measurements to estimate internal molten salt temperatures and predict flow velocity was explored, revealing that this approach could introduce up to a 5.5% error in flow velocity calculations. Line probes mapping temperature distributions across the tube's cross-section and molten salt provided valuable insights into temperature gradients, emphasizing the need for a thermal conductivity equation for molten salt with lower uncertainty to achieve more accurate temperature predictions of the system.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗