Search NASA⌕ Search

SEARCH · Search NASA

Results for “lumped thermal model”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

A title-gap flow model for use in aerodynamic loads assessment of space shuttle thermal protection system: Parallel gap faces

The problem of predicting aerodynamic loads on the insulating tiles of the space shuttle thermal protection system (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed, based upon a slow viscous (Stokes) flow analysis, and is verified against experimental data. The tile-gap pressure field is derived from a solution of the two-dimensional Laplace equation; the mass-flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculation. Analytical and experimental work to improve TPS flow predictions and a possible shuttle TPS hardware modification are recommended.

Dwoyer, D. L.↗
NASA NTRS1982

Flow through the tile gaps in the Space Shuttle Thermal Protection System

The problem of predicting aerodynamic loads on the insulating tiles of the Space Shuttle Thermal Protection System (TPS) is discussed and seen to require a method for predicting pressure and mass flux in the gaps between tiles. A mathematical model of the tile-gap flow is developed based upon a slow viscous (Stokes) flow analysis and is verified against available experimental data. This model derives the tile-gap pressure field from a solution of the two-dimensional Laplace equation; the mass flux vector is then calculated from the pressure gradient. The means for incorporating this model into a lumped-parameter network analogy for porous-media flow is also given. The flow model shows tile-gap mass flux to be very sensitive to the gap width indicating a need for coupling the TPS flow and tile displacement calculations. Finally recommendations are made concerning additional analytical and experimental work to improve TPS flow predictions.

Dwoyer, D. L.↗
NASA NTRS2004

Ventless pressure control of two-phase propellant tanks in microgravity

This work studies pressurization and pressure control of a large liquid hydrogen storage tank. A finite element model is developed that couples a lumped thermodynamic formulation for the vapor region with a complete solution of the Navier-Stokes and energy equations for the flow and temperature fields in the liquid. Numerical results show that buoyancy effects are strong, even in microgravity, and can reposition a vapor bubble that is initially at the center of the tank to a region near the tank wall in a relatively short time. Long-term tank pressurization with the vapor bubble at the tank wall shows that after an initial transient lasting about a week, the final rate of pressure increase agrees with a purely thermodynamic analysis of the entire tank. However, the final pressure levels are quite different from thermodynamic predictions. Numerical results also show that there is significant thermal stratification in the liquid due to the effects of natural convection. A subcooled jet is used to provide simultaneous cooling and mixing in order to bring the tank pressure back down to its initial value. Three different jet speeds are examined. Although the lowest jet speed is ineffective at controlling the pressure because of insufficient penetration into the liquid region, the highest jet speed is shown to be quite effective at disrupting thermal stratification and reducing the tank pressure in reasonable time.

Weightlessness↗
NASA NTRS1992

Galactic disks, infall, and the global value of Omega

Stringent limits on the current rate of infall of satellite systems onto spiral galaxies are set on the basis of the thinness and coldness of Galactic disks. For infalling satellites on isotropically oriented circular orbits, it is shown that, due to scattering, the thermal energy gain of the disk exceeds the satellite energy loss from dynamical friction by a factor of 1.6, with 25 percent deposited in z motion and 75 percent in planar motions. It is found that no more than 4 percent of the Galactic mass inside the solar radius can have accreted within the last 5 billion years, or else its scale and its Toomre Q-parameter would exceed observed values. In standard cold-dark-matter-dominated models for the growth of structure with Omega sub tot of 1, the mass accreted in dark matter lumps rises faster than t exp 2/3 and would exceed 28 percent in the last 5 Gyr. It is proposed that heating from satellite infall accounts for a substantial fraction of the increase of velocity dispersion and scale height with age that is observed in the Galaxy.

Toth, G.↗
NASA NTRS1989

The transient performance of a two-phase fluid reservoir

Thermal control of future large, high power spacecraft will require a two-phase fluid central bus. The two-phase fluid reservoir is a critical component in the two-phase fluid bus. It both controls the saturation temperature and provides a space for volumetric changes. A dynamic reservoir simulation model does not currently exist, but it is needed to expedite efforts and reduce risk. During 1989 an effort was made to develop a simulation model of the transient performance of a two-phase fluid reservoir. As a beginning, a preliminary model was developed. It is based upon component mathematical models in lumped parametric form and build upon five component mathematical models for calculating dynamic responses of two-phase fluid reservoirs, primary feedback elements, controller commands, heater actuators, and reservoir heaters. As much as possible, the model took advantage of the available SINDA'85/FLUINT thermal/fluid integrator. Additional calculation logic and computer subroutines were developed to complete implementation of the model. The model is capable of simulating dynamic response of an equilibrium two-phase fluid reservoir. Modification of the model to include the liquid/vapor nonequilibrium is required for applications of the model to simulate performance of reservoir in which the liquid and vapor phases of the reservoir fluid are not in equilibrium. In addition, the model in its present form, needs to be refined in several respects. More empirical data are needed to guide the model development. The model may then be used to conduct a full parametric study of two-phase fluid reservoirs. More complexities in two-phaes flow regions in laboratory and flight conditions may have to be considered eventually if empirical data cannot be simulated satisfactorily. System with other components arrangement also need to be simulated if optimization is ever to be attained. The present model does, however, preliminarily demonstrates that such analyses are quite possible and offers a far less expensive method to understand the transient of a two-phase fluid reservoir system than a totally headware approach.

Chi, Joseph↗
NASA NTRS2006

A Relict-Grain-Bearing Porphyritic Olivine Compound Chondrule from LL3.0 Semarkona that Experienced Limited Remelting

Chondrule D8n in LL3.0 Semarkona is a porphyritic olivine (PO) chondrule, 1300 x 1900 microns in size, with a complicated thermal history. The oldest recognizable portion of D8n is a moderately high-FeO, PO chondrule that is modeled as having become enmeshed in a dust ball containing a small, intact, low-FeO porphyritic chondrule and fine-grained material consisting of forsterite, kamacite, troilite, and possibly reduced C. The final chondrule melting event may have been a heat pulse that preferentially melted the low-FeO material and produced a low-FeO, opaque-rich, exterior region, 45-140 microns in thickness, around the original chondrule. Ai one end of the exterior region, a kamacite- and troilite-rich lump 960 pm in length formed. During the final melting event, the coarse, moderately ferroan olivine phenocrysts within the original chondrule appear to have been partly resorbed (These relict phenocrysts have the highest concentrations of FeO, MnO, and Cr2O3-7.5, 0.20, and 0.61 wt%, respectively-in D8n.). Narrow olivine overgrowths crystallized around the phenocrysts following final chondrule melting; their compositions seem to reflect mixing between melt derived from the exterior region and the resorbed margins of the phenocrysts. During the melting event, FeO in the relict phenocrysts was reduced, producing numerous small blebs of Ni-poor metallic Fe along preexisting curvilinear fractures. The reduced olivine flanking the trails of metal blebs has lower FeO than the phenocrysts but virtually identical MnO and Cr2O3 contents. Subsequent parent-body aqueous alteration in the exterior region of the chondrule formed pentlandite and abundant magnetite.

Rubin, Alan E.↗
NASA NTRS2017

Heat Loads Due to Small Penetrations in Multilayer Insulation Blankets

The main penetrations (supports and piping) through multilayer insulation systems for cryogenic tanks have been previously addressed by heat flow measurements. Smaller penetrations due to fasteners and attachments are now experimentally investigated. The use of small pins or plastic garment tag fasteners to each the handling and construction of multilayer insulation (MLI) blankets goes back many years. While it has long been understood that penetrations and other discontinuities degrade the performance of the MLI blanket, quantification of this degradation has generally been lumped into gross performance multipliers (often called degradation factors or scale factors). Small penetrations contribute both solid conduction and radiation heat transfer paths through the blanket. The conduction is down the stem of the structural element itself while the radiation is through the hole formed during installation of the pin or fastener. Analytical models were developed in conjunction with MLI perforation theory and Fouriers Law. Results of the analytical models are compared to experimental testing performed on a 10 layer MLI blanket with approximately 50 small plastic pins penetrating the test specimen. The pins were installed at 76-mm spacing inches in both directions to minimize the compounding of thermal effects due to localized compression or lateral heat transfer. The testing was performed using a liquid nitrogen boil-off calorimeter (Cryostat-100) with the standard boundary temperatures of 293 K and 78 K. Results show that the added radiation through the holes is much more significant than the conduction down the fastener. The results are shown to be in agreement with radiation theory for perforated films.

Structural attachments↗
NASA NTRS1993

A novel pulse technique for thermal diffusivity determination of high-temperature levitated materials - Inverse analysis and experimental set-up

In order to determine the thermal diffusivity of materials at high temperatures, a two-step extended flash technique, which is applicable to levitated spherical samples, is proposed. The containerless flash method is modeled as an axisymmetric transient conduction heat transfer problem within the sphere subjected to a radiative boundary condition on its surface. The problem of nonlinearity which arises from the radiative heat transfer boundary condition is solved by replacing it with the measured time-dependent surface temperature data, thus giving rise to an 'Inverse Analysis'. Upon obtaining the analytic solution for the temperature field, the determination of the thermal diffusivity turns into a minimization problem, whereby the difference between the actual surface boundary condition and the analytic solution is minimized. In performing the proposed experiments, there is a need to undertake a cool-down experiment which should be governed by the lumped parameter analysis.

Shen, F.↗
DOE OSTI2023

Model-based quantification of margins and uncertainties in metal additive manufacturing for process design and qualification

Laser powder bed fusion (LPBF) Additive Manufacturing (AM) has the potential to enable the production of components with novel designs and material properties unachievable otherwise. However, process repeatability is a challenge, making qualification ill-defined and greatly reducing the utility of what could be an important manufacturing technology. In this work, a combination of modeling, uncertainty quantification (UQ), and experimentation are used in an effort to predict and bound the range of possible outcomes of the LPBF process. Quantities of interest predicted are melt pool dimensions, microstructure features, and mechanical distortions. A combination of high fidelity thermal-fluid models, microstructure growth models, and reduced fidelity, rapid thermal and mechanical models are used. Uncertainty propagation techniques are used to predict probability distributions of quantities of interest from estimates of process uncertainties. Repeated experiments are done to quantify observed probability distributions and compared to predicted distributions to determine if predictions are precise and accurate. Novel modeling methods are microstrucutre characterization techniques are also discussed. It is found that high fidelity models do a generally good job bounding experimentally observed melt pool morphologies for both bead-on-plate and powder bed cases. Microstructure models are able to bound a number of experimentally observed microstructure statistics, but with low precision due to challenges with calibrating the microstructure growth model parameters. A developed modified inherent strain distortion model does not accurately predict observed distortions. A lumped laser distortion model shows promise in being both accurately and precisely bounding observed outcomes from the deflection comb build, but requires further evaluation on more builds and geometries.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Modeling a Transient Pressurization with Active Cooling Sizing Tool

As interest in the area of in-space zero boil-off cryogenic propellant storage develops, the need to visualize and quantify cryogen behavior during ventless tank self-pressurization and subsequent cool-down with active thermal control has become apparent. During the course of a mission, such as the launch ascent phase, there are periods that power to the active cooling system will be unavailable. In addition, because it is not feasible to install vacuum jackets on large propellant tanks, as is typically done for in-space cryogenic applications for science payloads, instances like the launch ascent heating phase are important to study. Numerous efforts have been made to characterize cryogenic tank pressurization during ventless cryogen storage without active cooling, but few tools exist to model this behavior in a user-friendly environment for general use, and none exist that quantify the marginal active cooling system size needed for power down periods to manage tank pressure response once active cooling is resumed. This paper describes the Transient pressurization with Active Cooling Tool (TACT), which is based on a ventless three-lump homogeneous thermodynamic self-pressurization model1 coupled with an active cooling system estimator. TACT has been designed to estimate the pressurization of a heated but unvented cryogenic tank, assuming an unavailable power period followed by a given cryocooler heat removal rate. By receiving input data on the tank material and geometry, propellant initial conditions, and passive and transient heating rates, a pressurization and recovery profile can be found, which establishes the time needed to return to a designated pressure. This provides the ability to understand the effect that launch ascent and unpowered mission segments have on the size of an active cooling system. A sample of the trends found show that an active cooling system sized for twice the steady state heating rate would results in a reasonable time for tank pressure recovery with ZBO of a liquid oxygen propellant tank.

Guzik, Monica C.↗
DOE OSTI2026

Dielectric Spectroscopy Cable NDE for Unequal Aging Over Different Lengths

This Pacific Northwest National Laboratory milestone report assesses the effect of unequal aging over different lengths of cable on dielectric spectroscopy (DS) bulk impedance measurements from the cable end. DS measurements can indicate the remaining useful life of cables; however, most benchmark tests are based on accelerated aging of the entire cable. If only a portion of the cable is exposed to a harsh environment (as is frequently the case), the DS measurement will indicate a significantly less-aged cable than if the entire cable were exposed to a uniform environmental stress. The DS measurement is primarily related to the intrinsic insulation material relative permittivity—typically between 1.5 and 3.0. Insulation aging increases permittivity and, correspondingly, the cable’s overall capacitance measured from the cable end. To judge the health of the overall cable system, operators are primarily interested in the condition of the most severely aged section since that is where a failure is most likely to occur. This can be expressed as the dimensionless relative permittivity of the severely aged section divided by the permittivity of a pristine cable. This relative permittivity is an intrinsic material parameter independent of cable length and can serve as a quantitative indication of the insulation condition. If the relative percentages of cable exposed to a harsh and benign environment are known, the DS measurement of the entire cable length, including both pristine and aged sections, can be compensated to predict the relative permittivity (relative to pristine permittivity) of the most severly aged segment. A test was performed on three 50 ft (16 m) cables with 10%, 50%, and 90% of the cable inside a thermal aging oven. Predictably, the aging effect on the DS response was greatest in the 90% sample, followed by the 50% sample, and then the 10% sample. The more interesting result is the ratio of the aged insulation permittivity to the pristine insulation permittivity. Based on prior work, an end-of-life cable has a value of approximately 1.35. If the relative percentages of pristine and severely aged cable lengths are known, and the total cable-length capacitance of the initial pristine cable and the combined pristine plus aged cable are known, the relative aged permittivity (ratio of aged section permittivity / pristine section permittivity) can be estimated. The relative cable lengths exposed to a harsh environment may be known or estimated based on plant layout drawings or cable reflectometry tests that can locate the oven entry and exit points. The compensation was verified using a lumped-parameter model to predict the aged-segment permittivity e2/e1, where e2 is the aged-segment permittivity and e1 is the pristine-segment permittivity. For all three percentages of aged cable lengths, the ratio was quite similar, as would be expected since all cable segments were exposed to the same aging environment. This ratio of aged permittivity to pristine permittivity can be compared to other damage indicating tests, including elongation at break, to assess cable condition.

ARENA Test Bed↗
DOE OSTI2025

Thermo-hydraulic steam pipe models for district heating simulations: Simplifications to balance accuracy and simulation speed

Steam piping networks are essential for optimizing performance in industrial processes and district heating systems. However, dynamic models that balance thermo-hydraulic accuracy with computational efficiency remain limited. In response, this paper presents a new discretized steam pipe model based on the plug flow approach, capturing key thermo-hydraulic behaviors while simplifying steam phase change processes. Implemented in Modelica, the model accurately calculates temperature and pressure distributions along steam pipelines. To improve computational efficiency for district-scale simulations, five model simplifications are introduced: lumped thermo-hydraulic functions, empirical correlations, fluid state approximations, steady-state dynamics and inclusion of flow derivatives. These simplified models achieve 85%-98% accuracy in predicting pressure drop and condensation losses, including dynamic condensate behavior during pipe warm-up—a factor often overlooked in existing models. The models support diverse network configurations, scaling effectively to systems with multiple distribution pipes and connected building loads. Discrete models provide detailed insights but exhibit a cubic increase in simulation time as the network scales by N connected building O(N 2.42 ). In contrast, lumped models simulate 10–28 times faster than discrete, offering quadratic scaling of simulation time O(N 1.73 ). However, they still require 6 times more computation time than a lossless network, highlighting the inherent computational challenges of modeling compressible fluid flow. In conclusion, the steady-state lumped variant, with its near-linear scalability in computational time O(N 1.01 ), emerges as an efficient solution for preliminary design evaluations and extensive parametric studies.

15 GEOTHERMAL ENERGY↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization

NASA’s near-future long-term space missions necessitate advancements in cryogenic fluid management (CFM), which includes safe and reliable long-term propellant storage. Consequently, NASA STMD (Space Technology Mission Directorate) has established the CFM Portfolio Project to improve CFM technologies for upcoming missions. As a part of CFM Modeling Portfolio, NASA Marshall Space Flight Center’s (MSFC) Fluid Dynamics branch within the Propulsion Systems Department is tasked with assessing and improving computational tools used to support flight projects such as Human Lander System and Commercial Lunar Payload Services. One of the challenging modeling problems is that of self-pressurization of propellant tank due to heat leakage over long time periods. Reduced order and nodal tools find it extremely difficult to accurately predict self-pressurization under transient conditions or where complex flow patterns or thermal gradients exist, and application of 3-D CFD (computational fluid dynamics) simulations is necessary to characterize these problems. Until recently, CFD simulations for these long-term processes (order of hours or days) have been too impractical to conduct due to prohibitive wall time and computational resource requirements. The requisite CFD tool need to be efficient, computationally scalable, modular with ability to incorporate various physics models, and robust enough to not accumulate conservation errors over several hours of simulated time. NASA MSFC’s Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. In this paper, we validate Loci-Stream for predicting self-pressurization of a flight scale propellant tank so it can serve as a reliable design and analysis tool for NASA’s CFM application needs. Liquid hydrogen tank pressurization tests carried out at the K-site testing facility provide a reliable data set for this purpose. These tests were simulated using Loci-Stream solver with VOF module as well as a hybrid approach which uses a lumped model for the ullage gas domain and CFD simulation of the liquid propellant. Both are shown to have very good predictive capabilities over multiple K-site experiments.

Hong Q. Yang↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization

NASA's near-future long-term space missions necessitate advancements in cryogenic fluid management (CFM), which includes safe and reliable long-term propellant storage. Consequently, NASA STMD (Space Technology Mission Directorate) has established the CFM Portfolio Project to improve CFM technologies for upcoming missions. As a part of CFM Modeling Portfolio, NASA Marshall Space Flight Center's (MSFC) Fluid Dynamics branch within the Propulsion Systems Department is tasked with assessing and improving computational tools used to support flight projects such as Human Lander System and Commercial Lunar Payload Services. One of the challenging modeling problems is that of self-pressurization of propellant tank due to heat leakage over long time periods. Reduced order and nodal tools find it extremely difficult to accurately predict self-pressurization under transient conditions or where complex flow patterns or thermal gradients exist, and application of 3-D CFD (computational fluid dynamics) simulations is necessary to characterize these problems. Until recently, CFD simulations for these long-term processes (order of hours or days) have been too impractical to conduct due to prohibitive wall time and computational resource requirements. The requisite CFD tool need to be efficient, computationally scalable, modular with ability to incorporate various physics models, and robust enough to not accumulate conservation errors over several hours of simulated time. NASA MSFC's Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. In this paper, we validate Loci-Stream for predicting self-pressurization of a flight scale propellant tank so it can serve as a reliable design and analysis tool for NASA's CFM application needs. Liquid hydrogen tank pressurization tests carried out at the K-site testing facility provide a reliable data set for this purpose. These tests were simulated using Loci-Stream solver with VOF module as well as a hybrid approach which uses a lumped model for the ullage gas domain and CFD simulation of the liquid propellant. Both are shown to have very good predictive capabilities over multiple K-site experiments.

cryogenic fluid management↗
DOE OSTI2021

Full-Loop Simulation of the Combustion of Biomass in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency, and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. This report summarizes the development of a comprehensive computational fluid dynamics (CFD) model for biomass combustion in a 50 kW th CFB combustor employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite at the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL). The experimental combustor was designed, built, and operated at CanmetENERGY, Natural Resources Canada in Ottawa, Canada. This work is the culmination of a multi-year collaboration between NETL and CanmetENERGY. The hydrodynamics in the full loop are validated against cold-flow experiments conducted at CanmetENERGY. To mimic the dense packing in the standpipe in the absence of a collisional model in the MFiX-PIC approach, the ram valve in the return leg is partially closed to obstruct the flow. A combustion reaction scheme was developed that considers the volatiles released by the fuel as a lumped species that fully converts to CO 2 and H 2 O. The so-called enhanced simplified combustion scheme produces an excellent match with experimental data in terms of the species concentrations of CO 2 and O 2 at the outlet and the freeboard, demonstrating the utility of the comprehensive model for a range of flow conditions and operating temperatures. Furthermore, the model can be readily adapted to simulate combustion of a variety of biomass fuels as well as co-fired systems under air and oxy-fuel conditions.

09 BIOMASS FUELS↗
DOE OSTI2024

A study on the impact of using a subchannel resolution for modeling of large break loss of coolant accidents

The nuclear industry is investigating the feasibility of transitioning from 18- to 24-month fuel cycles because of the positive impact it would have on the operational costs for the current fleet of light-water reactors. A challenge to making this change is the increased risk of fuel fragmentation, relocation, and dispersal (FFRD) due to the known potential for ceramic fuel to pulverize into fine particles at the higher discharge burnups. Previous work has been performed by the Nuclear Energy Advanced Modeling and Simulation program to assess FFRD risk in high-burnup cores using the BISON fuel performance code and a coarse mesh thermal hydraulics (T/H) solution for a loss-of-coolant accident (LOCA) using the TRACE system T/H code. Because of the importance of the T/H solution for FFRD assessment, this study seeks to investigate the impact of using higher-fidelity subchannel techniques for modeling of the LOCA transient. CTF was used to model a subregion of a high-burnup core that was depleted by the Virtual Environment for Reactor Applications (VERA) multiphysics core simulator. Both coarse-mesh and pin-resolved models were created in CTF, and a consistent coarse-mesh TRACE model was also developed to allow for benchmarking the code results. Further, a large-break loss-of-coolant accident (LBLOCA) reflood transient was simulated using these three models, and results were compared. Results showed some consistent differences between the CTF and TRACE coarse models, including a higher peak cladding temperature (PCT) prediction in CTF and later quenching in CTF; however, the transient clad temperature behavior was similar, and these differences are likely due to post-critical heat flux heat transfer modeling differences and minimum film boiling temperature model differences. The pin-resolved results indicate that the PCT in the lumped model is often under-predicted by as much as 70 °C and that PCT occurs at a different location than the high-power pin in the assembly. The lumped model predicts a difference of 10 °C or less between the average and hot pins in the assembly, whereas the pin-resolved model predicts a range of over 100 °C. These results indicate that higher-fidelity T/H results may have an impact on predicted core behavior during LOCA, which may be important to consider when assessing FFRD risk.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗