SEARCH · Engineering Papers

Results for “mesh data”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Experimental Testing of Particle Erosion and Attrition in CSP Systems

In CSP systems using solid particles for heat transfer, erosion from the falling particles can damage the system components upon impact, and from sliding against the component walls. Furthermore, the resulting attrition of particles can change the thermal performance due to alterations in particle morphology. The research on these issues has been limited to industrial applications, and don’t necessarily coincide with the operating conditions in CSP systems. In this work, we investigated the durability of the particles selected as the heat transfer media, as well as special substrate materials that are of interest as lining materials to the particle storage tanks and hoppers. Three different erosion phenomena were studied: impact, abrasion, and attrition erosion. The low-temperature experiments revealed negligible losses from attrition and abrasion erosion. Although the substrates indicated visible surface damage from impact erosion, the actual material mass loss was small even after a throughput of 500kg of particles 63 days of 1MW CSP plant operation assuming 14 hours per day of operation. Erosion is typically classified as impact, abrasion, or attrition erosion. For example, the CSP systems utilizing solid particles as heat transfer media will likely experience impact erosion on the receiver hopper and heat exchanger walls as the solid particles fall from the reservoir. Additionally, there will be sliding erosion along the receiver walls as particles move through the hopper, and lastly, there will also be attrition erosion as the particles rub against each other. The extent of erosion depends on many factors including but not limited to particle shape and size, material hardness/brittleness and ductility, and particle impact speed and angle. Mechanistic erosion equations that are available in the literature are mostly derived from data pertaining to specific experiments and therefore can’t be applied to all applications without suitable adjustments. Survey of the literature relating to particle-based receivers revealed that the most prevalent approach has been to choose low velocities to limit erosion, but to our knowledge, little effort has been put in to fully characterize erosion at the expected operating conditions. A relatively recent work investigated the erosion of wire mesh proposed for use in a particle receiver, where the results actually showed an increase in mass as oxidation occurred on the mesh surface but the work did not investigate the particle attrition rate, and the erosion rate of a solid substrate. The end goal of this work is to develop a comprehensive particle and substrate durability model that will allow a comprehensive understanding of particle-based CSP operating conditions on component durability.

14 SOLAR ENERGY↗

DOE OSTI2017

ASSESSMENT OF MULTIPHASE CFD WITH ZERO CLOSURE MODEL FOR BOILING WATER REACTOR FUEL ASSEMBLIES

Extension of existing boiling closures to the high void fraction conditions of BWRs is one of challenges for M-CFD applications. This paper presents the assessment of an Eulerian M-CFD model with Zero closure model that has been specifically designed for two-phase flow and heat transfer applications in BWR fuel assemblies. The baseline Zero closure aims at including the simplest and most robust representation of the key physical mechanisms, rather than aiming at absolute accuracy of the predictions. The capabilities of M-CFD with Zero were evaluated against the international OECD/NRC BFBT benchmark data through the comparative analyses of local void fraction, sub-channel void fraction error and exit quality. INL high performance computing system has been employed to demonstrate the scalability of the M-CFD application to nuclear reactor fuel assemblies. The M-CFD simulation results show a good agreement with the experimental data for both quantitative and qualitative results, with the highest uncertainty being observed in the prediction of the local void fraction profile at low void fraction condition. This paper presents a recent view of the best practices for M-CFD of BWR fuel assemblies to provide a useful reference to identify and overcome common challenges in the applications of M-CFD tools. We provide recommendations for geometric modeling, mesh generation and solver configurations to optimize the methodology application.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2019

Stress field analysis in a stony meteorite under thermal fatigue and mechanical loadings

A micromechanical approach relying on 3D image-based finite element (FE) models, together with experimental testing, is employed to investigate the stress distribution in an L6 ordinary chondrite's microstructure subjected to thermal and mechanical loadings. An in-situ X-ray micro-computed tomography (micro-CT) setup capable of applying thermal and mechanical loads shed light on crack initiation mechanisms in chondrite samples. Micro-CT images were further used to reconstruct realistic microstructural models of the samples, which were transformed into FE models using an automated non-iterative mesh generation algorithm. FE approximation of stress fields showed that under mechanical loads, the stress often concentrates in matrix-particle interfacial regions with higher curvatures. Furthermore, this could result in the nucleation and eventually propagation of cracks in such regions. However, when the thermal cycling takes place, the stress concentrates more uniformly along with particle interfaces, resulting in their debonding from the surrounding matrix. Both observations were verified by experimental data obtained from in-situ micro-CT imaging.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Compatible, energy conserving, bounds preserving remap of hydrodynamic fields for an extended ALE scheme

From the very origins of numerical hydrodynamics in the Lagrangian work of von Neumann and Richtmyer [83], the issue of total energy conservation as well as entropy production has been problematic. Because of well known problems with mesh deformation, Lagrangian schemes have evolved into Arbitrary Lagrangian–Eulerian (ALE) methods [39] that combine the best properties of Lagrangian and Eulerian methods. Energy issues have persisted for this class of methods. We believe that fundamental issues of energy conservation and entropy production in ALE require further examination. The context of the paper is an ALE scheme that is extended in the sense that it permits cyclic or periodic remap of data between grids of the same or differing connectivity. The principal design goals for a remap method then consist of total energy conservation, bounded internal energy, and compatibility of kinetic energy and momentum. We also have secondary objectives of limiting velocity and stress in a non-directional manner, keeping primitive variables monotone, and providing a higher than second order reconstruction of remapped variables. Particularly, the new contributions fall into three categories associated with: energy conservation and entropy production, reconstruction and bounds preservation of scalar and tensor fields, and conservative remap of nonlinear fields. Our paper presents a derivation of the methods, details of implementation, and numerical results for a number of test problems. The methods requires volume integration of polynomial functions in polytopal cells with planar facets, and the requisite expressions are derived for arbitrary order.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

Assessing the geomechanical stability of interbedded hydrate-bearing sediments under gas production by depressurization at NGHP-02 Site 16

Establishing the geomechanical stability of marine sediments in the vicinity of a production well is one of the key design considerations in planning offshore gas production from marine hydrate reservoirs. This work presents an assessment of the sediment stability at India's National Gas Hydrate Program, Expedition 2 (NGHP-02) Site 16 Area B offshore eastern India, for which gas production is to be carried out by depressurization. One important feature of the study is that extensive calibration of constitutive model parameters has been conducted based on laboratory test data from pressured core samples. From analysis perspective, the site is challenging because the hydrate reservoir consists of thin layers of hydrate-bearing sands interbedded with mud. Moreover, depressurization at the depth of a reservoir more than 2750 m below sea surface will lead to a pore pressure drop, and accordingly an effective confining stress increase as high as 25 MPa. In dealing with thin interbedded hydrate-bearing strata, meshing requirements for flow and geomechanical analysis are quite different from those for reservoirs with thicker massive layers, An axisymmetric model and one-way coupling simulations were thus adopted for this study, in which the geomechanical study utilizes pore pressure and hydrate saturation output from the flow study, but the flow study does not takes the porosity changes from the geomechanical analysis. Instead, the reduction of porosity due to sediment deformation in the flow study is based on a pressure-dependent pore compressibility relationship derived from geomechanical modeling. The rationality is validated through back computing the pore compressibility from the geomechanical deformation results. Finally, the study shows that large compression in the reservoir will result in movement of the sediments from above and below, as well as laterally in smaller magnitudes; and the sediment is deemed stable during the gas production period.

58 GEOSCIENCES↗

DOE OSTI2019

Evaluating Large-Eddy Simulations of Boundary-Layer Cloud Response to Climate Change using MAGIC Observations

Over the cool waters off California, a persistent and extensive deck of low-lying clouds forms atop a shallow, turbulent layer of moist air extending 500-2000 meters above the ocean, above which the air is much warmer and drier. Its west edge varies with weather and season and is typically half way to Hawaii. This cloud deck and similar cloud decks elsewhere reflect sunlight back to space, cooling the underlying ocean and indeed the entire planet. They must be accurately represented by a climate model to correctly simulate our current climate or future global warming. Most climate models suggest these cloud decks will thin and break up in a warmer CO 2 -enriched climate, so extra sunlight would reach the ocean surface and amplify warming worldwide. But these models have difficulty accurately simulating the cloud thickness and extent because the turbulent cloud-forming air currents are much too small to be simulated by their coarse computational grid. In 2012-2013, the DOE Atmospheric Radiation Measurement MAGIC program instrumented a container ship transiting between California and Hawaii to observe clouds and rain, small particles (aerosols), and weather conditions. Data were obtained from 20 4-6 day cruises in each direction. This study tested whether a sophisticated ‘large-eddy simulation’ (LES) of the turbulent air motions that create the clouds in this region accurately predicts the variations in cloud properties seen across the different cruises. If so, that model is more believable for predicting cloud changes in this region in a future warmer climate. Our LES used a numerical mesh over a 6.4x6.4 km region with 50 m horizontal and 5 m vertical grid spacing to simulate the cloud-forming eddies and the sharp overlying transition between moist and dry air. To better compare with MAGIC observations, we configured our model to simulate a region that moved along with the ship. We used a weather forecast model to specify effects of the large-scale air flow. We compared our results with balloon ascents, a radar and other instruments pointed upward at the clouds, measurements of downward sunlight and infrared radiation, aerosols, and weather variables for all suitable cruise periods. Our LES was remarkably skillful in predicting the cloudiness variations measured by the ship. MAGIC also sampled mesoscale cells, regions of cloud thickening and thinning over distances of 10-50 km. An LES run over a larger region 100 km on a side shows these cells form spontaneously, enhanced by strong radiative heat loss at the cloud tops.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1993

Micromechanics of crenulated fibers in carbon/carbon composites

The influence of crenulated noncircular fibers on the micromechanical stress states due to a transverse strain and to a temperature change in carbon/carbon composites is examined using the finite element method. Stresses at the interface of both fully bonded and fully disbonded fibers having two crenulation amplitudes and with two fiber volume fractions are presented. In each case, these interface stresses are compared to stresses at the interface of circular fibers which have the same degree of disbond and fiber volume fraction and are under the same loading conditions. For the disbonded cases, deformed meshes showing locations of fiber/matrix contact are also included. In addition to the interface stress states, selected composite properties are also computed and compared in each case examined. Interest in studying noncircular fibers stems from a desire to increase the transverse properties of carbon/carbon by introducing a mechanical interlocking between the fiber and the matrix. Results presented here indicate that this interlocking does in fact occur. Evidence from the interface stress data suggests, however, that any possible advantage of this interlocking may be outweighed by the disadvantage of stress concentrations which arise at the interface due to the crenulated geometry of the fibers.

Carapella, E. E.↗

NASA NTRS1998

Numerical Modeling Studies of Wake Vortex Transport and Evolution Within the Planetary Boundary Layer

In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Lin, Yuh-Lang↗

NASA NTRS2006

Thermal Dispersion Within a Porous Medium Near a Solid Wall

The regenerator is a key component to Stirling cycle machine efficiency. Typical regenerators are of sintered fine wires or layers of fine-wire screens. Such porous materials are contained within solid-waH casings. Thermal energy exchange between the regenerator and the casing is important to cycle performance for the matrix and casing would not have the same axial temperature profile in an actual machine. Exchange from one to the other may allow shunting of thermal energy, reducing cycle efficiency. In this paper, temperature profiles within the near-wall region of the matrix are measured and thermal energy transport, termed thermal dispersion, is inferred. The data show how the wall affects thermal transport. Transport normal to the mean flow direction is by conduction within the solid and fluid and by advective transport within the matrix. In the near-wall region, both may be interrupted from their normal in-core pattern. Solid conduction paths are broken and scales of advective transport are damped. An equation is presented which describes this change for a wire screen mesh. The near-wall layer typically acts as an insulating layer. This should be considered in design or analysis. Effective thermal conductivity within the core is uniform. In-core transverse thermal effective conductivity values are compared to direct and indirect measurements reported elsewhere and to 3D numerical simulation results, computed previously and reported elsewhere. The 3-D CFD model is composed of six cylinders in cross flow, staggered in arrangement to match the dimensions and porosity of the matrix used in the experiments. The commercial code FLUENT is used to obtain the flow and thermal fields. The thermal dispersion and effective thermal conductivities for the matrix are computed from the results.

Simon, T.↗

NASA NTRS2015

Parametric Design of Injectors for LDI-3 Combustors

Application of a partially calibrated National Combustion Code (NCC) for providing guidance in the design of the 3rd generation of the Lean-Direct Injection (LDI) multi-element combustion configuration (LDI-3) is summarized. NCC was used to perform non-reacting and two-phase reacting flow computations on several LDI-3 injector configurations in a single-element and a five-element injector array. All computations were performed with a consistent approach for mesh-generation, turbulence, spray simulations, ignition and chemical kinetics-modeling. Both qualitative and quantitative assessment of the computed flowfield characteristics of the several design options led to selection of an optimal injector LDI- 3 design that met all the requirements including effective area, aerodynamics and fuel-air mixing criteria. Computed LDI-3 emissions (namely, NOx, CO and UHC) will be compared with the prior generation LDI- 2 combustor experimental data at relevant engine cycle conditions.

injectors↗

NASA NTRS2015

Impedance Eduction in Sound Fields With Peripherally Varying Liners and Flow

A two-dimensional impedance eduction theory is extended to three-dimensional sound fields and peripherally varying duct liners. The approach is to first measure the acoustic pressure field at a series of flush-mounted wall microphones located around the periphery of the flow duct. The numerical solution for the acoustic pressure field at these microphones is also obtained by solving the three-dimensional convected Helmholtz equation using the finite element method. A quadratic objective function based on the difference between the measured and finite element solution is constructed and the unknown impedance function is obtained by minimizing this objective function. Impedance spectra educed for two uniform-structure liners (a wire-mesh and a conventional liner) and a hard-soft-hard peripherally varying liner (for which the soft segment is that of the conventional liner) are presented. Results are presented at three mean flow Mach numbers and fourteen sound source frequencies. The impedance spectra of the uniform-structure liners are also computed using a two-dimensional impedance eduction theory. The primary conclusions of the study are: 1) when measured data is used with the uniform-structure liners, the three-dimensional theory reproduces the same impedance spectra as the two-dimensional theory except for frequencies corresponding to very low or very high liner attenuation; and 2) good agreement between the educed impedance spectra of the uniform structure conventional liner and the soft segment of the peripherally varying liner is obtained.

Watson, W. R.↗

NASA NTRSDate not supplied

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

This paper presents the results of a computational study on porous-plug nozzles for their eventual application in reducing takeoff noise levels for supersonic civil transport. The study explores a wide range of porous-plug geometries, encompassing variations in hole size, plug length and porosity. Computational fluid dynamics (CFD) analyses are employed in concert with an advanced automated mesh refinement (AMR) scheme, Sketch-to-Solution (S2S), to investigate porous plug nozzle design strategies that may diminish shock strength, broad-band shock noise (BBSN), flow separation, and other undesirable characteristics associated with high takeoff noise levels. A selected subset of the porous-plug nozzles that were analyzed computationally were manufactured at subscale using stereolithography 3D printing. These nozzles were tested experimentally, obtaining far-field noise spectral data and schlieren flow visualization pictures. Detailed comparisons between CFD predictions and experimental results were conducted. These comparisons focus on gaining a better understanding of the underlying physical mechanisms responsible for noise reduction with the porous plugs relative to corresponding solid plugs.

computational fluid dynamics↗

NASA NTRSDate not supplied

Investigation of Physical Mechanisms for Jet Noise Reduction by Plug Nozzle Porosity

Nozzles↗

NASA NTRS1978

Monthly mean simulation experiments with a course-mesh global atmospheric model

Substitution of observed monthly mean sea-surface temperatures (SSTs) as lower boundary conditions, in place of climatological SSTs, failed to improve the model simulations. While the impact of SST anomalies on the model output is greater at sea level than at upper levels the impact on the monthly mean simulations is not beneficial at any level. Shifts of one and two days in initialization time produced small, but non-trivial, changes in the model-generated monthly mean synoptic fields. No improvements in the mean simulations resulted from the use of either time-averaged initial data or re-initialization with time-averaged early model output. The noise level of the model, as determined from a multiple initial state perturbation experiment, was found to be generally low, but with a noisier response to initial state errors in high latitudes than the tropics.

Spar, J.↗

DOE OSTI2022

Numerical and experimental investigation of the flame kernel growth in a methane/air mixture near the lean flammability limit

Lean combustion has the potential to improve the thermal efficiency of spark-ignition engines, but it faces the significant challenge of increased cycle-to-cycle variation due to low mixture reactivity and unstable flame dynamics. Computational fluid dynamics (CFD) employing predictive models can guide engine design and optimize operating strategies for lean combustion. However, ignition and combustion models have rarely been validated at fuel-lean conditions, and a fundamental understanding of the early flame kernel growth process is also lacking for a successful sub-model development. Here, the present study develops a numerical simulation framework used to investigate early flame kernel growth in methane/air mixtures. A nanosecond-pulsed discharge (NPD) approach is employed to effectively decouple the flame kernel growth from the electrical discharge due to their difference in timescales, and equivalence ratios near the experimentally measured lean flammability limit (LFL) are selected to focus on challenging mixture conditions. Three numerical investigations, such as the choice of turbulence modeling, grid size, and grid control strategies, are examined to match both LFL and flame kernel structure measured from experiments. It is demonstrated that a quasi-direct numerical simulation (QDNS) with a fixed grid embedding of 10 μm can predict the LFL as φ CFD =0.61 and match the displacement speed of the kernel’s boundary marked in schlieren images. To predict the LFL and flame kernel shape, a fine grid (Δ≤12.5 μm) is needed to capture the consumption of formaldehyde (CH 2 O) in kernel’s reaction branches attached to the anode, and adaptive mesh refinement is replaced with the fixed embedding due to loss of simulation accuracy. Also, it is found that a large-eddy simulation (LES) using the Dynamic Structure model is not suitable for the NPD-induced flame kernel simulation because artificial sub-grid turbulent kinetic energy induced by shock dynamics alters the flow velocity calculation, resulting in divergence of LES from QDNS. Lastly, the simulation well matches the experimental data for the flame kernel evolution in three mixture conditions (φ = 0.7, 0.61, 0.55), showing toroidal flame kernel expansion and flame kernel growth/extinction.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

CASL VMA Milestone Report FY16 (L3:VMA.VUQ.P13.08): Westinghouse Mixing with STAR-CCM+

STAR-CCM+ (STAR) is a high-resolution computational fluid dynamics (CFD) code developed by CD-adapco. STAR includes validated physics models and a full suite of turbulence models including ones from the k-ε and k-ω families. STAR is currently being developed to be able to do two phase flows, but the current focus of the software is single phase flow. STAR can use imported meshes or use the built in meshing software to create computation domains for CFD. Since the solvers generally require a fine mesh for good computational results, the meshes used with STAR tend to number in the millions of cells, with that number growing with simulation and geometry complexity. The time required to model the flow of a full 5x5 Mixing Vane Grid Assembly (5x5MVG) in the current STAR configuration is on the order of hours, and can be very computationally expensive. COBRA-TF (CTF) is a low-resolution subchannel code that can be trained using high fidelity data from STAR. CTF does not have turbulence models and instead uses a turbulent mixing coefficient β. With a properly calibrated β, CTF can be used a low-computational cost alternative to expensive full CFD calculations performed with STAR. During the Hi2Lo work with CTF and STAR, STAR-CCM+ will be used to calibrate β and to provide high-resolution results that can be used in the place of and in addition to experimental results to reduce the uncertainty in the CTF results.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗