Search NASA⌕ Search

SEARCH · Search NASA

Results for “active turbulence”

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

Origin of the ring ellipticity in the black hole images of M87*

We investigate the origin of the elliptical ring structure observed in the images of the supermassive black hole M87*, aiming to disentangle contributions from gravitational, astrophysical, and imaging effects. Leveraging the enhanced capabilities of the Event Horizon Telescope (EHT)'s 2018 array, including improved (u,v)-coverage from the Greenland Telescope, we measured the ring's ellipticity using five independent imaging methods, obtaining a consistent average value of τ = 0.08−0.02+0.03 with a position angle of ξ = 50.1−7.6+6.2 degrees. To interpret this measurement, we compared it to general relativistic magnetohydrodynamic (GRMHD) simulations spanning a wide range of physical parameters including the thermal or nonthermal electron distribution function, spins, and ion-to-electron temperature ratios in both low- and high-density regions. We find no statistically significant correlation between spin and ellipticity in GRMHD images. Instead, we identify a correlation between ellipticity and the fraction of non-ring emission, particularly in nonthermal models and models with higher jet emission. These results indicate that the ellipticity measured from the M87* emission structure is consistent with that expected from simulations of turbulent accretion flows around black holes, where it is dominated by astrophysical effects rather than gravitational ones. Future high-resolution imaging, including space very long baseline interferometry and long-term monitoring, will be essential to isolate gravitational signatures from astrophysical effects.Key words: accretion, accretion disks / black hole physics / gravitation / galaxies: active** These authors have contributed equally to this work.*** NASA Hubble Fellowship Program, Einstein Fellow.† Deceased.

Dahale, Rohan [Granada U., Theor. Phys. Astrophys.↗
DOE OSTI2025

Optimization and active stabilization of a far-infrared laser for NSTX-U high poloidal wavenumber scattering diagnostics

The far-infrared (FIR) laser output beam power and profile are important parameters in the laser-aided diagnostics, directly influencing the spatial resolution and signal-to-noise ratio of measurements. Here, this work focuses on developing a systematic control method to enhance FIR laser beam quality through optimized mirror alignment and real-time feedback-based precision cavity length tuning. A 150 W CO 2 laser, aligned with the waveguide axis using a HeNe reference laser, serves as the pump source. The sensitivity of FIR beam intensity to pump gas pressure and thermal expansion is investigated, revealing that even a 1 µm cavity expansion can significantly degrade output power stability to about two-thirds of its original value. To address this, a feedback control module has been designed and implemented for active cavity length adjustment, stabilizing the output power at ∼30 mW. In addition, maintaining a high formic acid gas pressure ($>$190 mTorr) within the cavity ensures reliable operation. The optimized FIR laser will be deployed on the National Spherical Torus eXperiment-U high poloidal wavenumber scattering system for studying electron-scale turbulence in tokamak plasmas.

Xu, Xinhang [Univ. of California, Davis, CA (Unite↗
DOE OSTI2020

Impact of resonant magnetic perturbations on zonal flows and microturbulence

Results from a joint experimental and computational effort studying the effect of resonant magnetic perturbations (RMPs) on microturbulence levels and their connection to zonal flows in the DIII-D tokamak L-mode are presented. Beam emission spectroscopy measurements show a direct increase in density fluctuations at microturbulent scales with increasing RMP amplitude, suggesting that magnetic activity introduced by the RMP affects the regulation of microturbulence on DIII-D. This is analogous to how MHD-scale magnetic fluctuations arising from tearing modes have been observed in simulations to increase microturbulence levels in the reversed-field pinch (RFP). In the RFP, this is attributed to magnetic fluctuations eroding turbulence-limiting zonal flows; this work examines if a similar mechanism is present for DIII-D microturbulence. Gyrokinetic simulations find that the application of an RMP corresponds directly to a decrease in zonal flow levels, producing a similar increase of turbulent fluctuation levels over a range of RMP amplitudes as observed in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

CFD analysis of a solid pin-fueled small modular fluoride salt-cooled reactor

A computational fluid dynamics (CFD) study is performed for the pre-conceptual solid pin-fueled small modular fluoride salt-cooled reactor developed by the Oak Ridge National Laboratory (ORNL). The solid and fluid regions in a 1/12th section of a fuel assembly are modeled. The solid region includes fuel and non-fuel pins, with heat generation in the active region of the fuel pins. Two different power profiles, uniform and center-peaked, are considered. The FLiBe coolant flows from the bottom to the top of the core, parallel to the bank of fuel and non-fuel pins. The k-ω shear-stress transport (SST) model is chosen as the baseline turbulence model. The effects of grid refinement, inlet turbulence specification, and turbulence models on the temperature and pressure drop predictions are studied. Turbulence model sensitivity is investigated by comparing the results from the k-ω SST model with other two-equation models (k-ω baseline or BSL and k-ε realizable) as well as with anisotropic Reynolds stress transport models (linear pressure-strain and stress-BSL). Here the results show that the choice of turbulence model has a significant impact on the pin temperature and bundle pressure drop predictions. The results from the baseline turbulence model show good agreement between the fuel pin temperatures and bundle pressure drop values predicted using a subchannel model previously developed by the authors of the present study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Mallat Scattering Transformation based surrogate for Magnetohydrodynamics

Abstract A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x , y , and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Outlook for active plasma spectroscopy with entangled two-photon absorption

Entangled two-photon absorption (ETPA) exhibits a quantum-enhanced two-photon absorption cross-section and a linear scaling with incident power that exceeds classical two-photon absorption (TPA) at low photon flux. Here, we assess the feasibility and applications of ETPA for active plasma spectroscopy. The ETPA-absorbed power for an Ar + species in a helicon plasma scenario is calculated to be as high as 0.6 nW cm −1 along the incident beam path for nondegenerate-entangled photon pairs that are near-resonant with an intermediate level. The ETPA absorption power exceeds classical TPA and Bremsstrahlung emission, a necessary condition for a feasible fluorescence measurement. Applications of ETPA plasma spectroscopy may include trace impurity detection, ground-state excitation, and continuous-wave pumping for high-bandwidth turbulence measurements. Possible benefits of plasma measurements are less harm to optical components with lower incident flux and reduced stray light that can contaminate optical diagnostics.

active plasma spectroscopy↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Intermittent Frontogenesis in the Alboran Sea

Abstract We present a phenomenological description and dynamical analysis of the Alboran fronts using a realistic simulation at submesoscale resolution. The study is focused on east Alboran fronts emerging within relatively strong flows that separate from the Spanish coast into the basin interior. Despite modest lateral shifting associated with the position of the Alboran anticyclonic gyres and variations in intensity, these fronts present a similar structure and dynamical configuration as the climatological Almeria–Oran front. The statistical analysis of our solution shows that strained-induced frontogenesis is a recurrent submesoscale mechanism associated with these fronts, and the process is assessed in terms of the advective Lagrangian frontogenetic tendencies associated with buoyancy and velocity horizontal gradients. Intermittency in their strength and patterns is indicative of high variability in the occurrence of active frontogenesis in association with the secondary (overturning) circulation across the frontal gradient. As a result, we find many episodes with strong surface fronts that do not have much associated downwelling. Frontogenesis and the associated secondary circulation are further explored during two particular frontal events, both showing strong downwelling of (1) cm s −1 extending down into the pycnocline. A frontogenetic contribution of turbulent vertical momentum mixing to the secondary circulation is identified in the easternmost region during the cold season, when the dynamics are strongly influenced by the intrusion of the salty Northern Current. The background vertical velocity fields observed during the analyzed events indicate other currents in the submesoscale range, including tidal and topographic internal waves.

Capó, Esther↗
DOE OSTI2019

Predicting the impact of particle-particle collisions on turbophoresis with a reduced number of computational particles

A common feature of wall-bounded turbulent particle-laden flows is enhanced particle concentrations in a thin layer near the wall due to a phenomenon known as turbophoresis. Even at relatively low bulk volume fractions, particle-particle collisions regulate turbophoresis in a critical way, making simulations sensitive to collisional effects. Lagrangian tracking of every particle in the flow can become computationally expensive when the physical number of particles in the system is large. Artificially reducing the number of particles in the simulation can mitigate the computational cost. When particle-particle collisions are an important aspect determining the simulation outcome, as in the case when turbophoresis plays an active role, simply reducing the number of particles in the simulation significantly alters the computed particle statistics. This paper introduces a computational particle treatment for particle-particle collisions which reproduces the results of a full simulation with a reduced number of particles. This is accomplished by artificially enhancing the particle collision radius based on scaling laws for the collision rates. The proposed method retains the use of deterministic collision models and is applicable for both low and high Stokes number regimes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Involute Working Group – Validation of CFD Turbulence Models for Steady-State Safety Analysis of Corner Geometry

Worldwide, three research reactors have fuel plates curved as a circle-involute (a spiral generated around a circle): Oak Ridge National Laboratory High Flux Isotope Reactor located in Tennessee, U.S.A. [1], Institut Laue-Langevin High Flux Reactor located in Grenoble, France [2], and Technical University of Munich Research Neutron Source Heinz Maier-Leibnitz located in Garching, Germany [3]. All three reactors are currently using highly enriched uranium ( 235 U/U ≥ 20 wt%) as fuel, and all three are actively engaged in activities to convert to low-enriched uranium ( 235 U/U < 20 wt%) fuel.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

A Model Intercomparison Study of Aerosol‐Cloud‐Turbulence Interactions in a Cloud Chamber: 1. Model Results

This study presents the first model intercomparison of aerosol‐cloud‐turbulence interactions in a controlled cloudy Rayleigh‐Bénard Convection chamber environment, utilizing the Pi Chamber at Michigan Technological University. We analyzed simulated cloud chamber‐averaged statistics of microphysics and thermodynamics in a warm‐phase, cloudy environment under steady‐state conditions at varying aerosol injection rates. Simulation results from seven distinct models (DNS, LES, and a 1D turbulence model) were compared. Our findings demonstrate that while all models qualitatively capture observed trends in droplet number concentration, mean radius, and droplet size distributions at both high and low aerosol injection rates, significant quantitative differences were observed. Notably, droplet number concentrations varied by over two orders of magnitude between models for the same injection rates, indicating sensitivities to the model treatments in droplet activation and removal and wall fluxes. Furthermore, inconsistencies in vertical relative humidity profiles and in achieving steady‐state liquid water content suggest the need for further investigation into the mechanisms driving these variations. Despite these discrepancies, the models generally reproduced consistent power‐law relationships between the microphysical variables. This model intercomparison underscores the importance of controlled cloud chamber experiments for validating and improving cloud microphysical parameterizations. Recommendations for future modeling studies are also highlighted, including constraining wall conditions and processes, investigating droplet/aerosol removal (including sidewall losses), and conducting simplified experiments to isolate specific processes contributing to model divergence and reduce model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Modeling the effects of active wake mixing on wake behavior through large-scale coherent structures

The use of active wake mixing (AWM) to mitigate downstream turbine wakes has created new opportunities for reducing power losses in wind farms. However, many current analytical or semi-empirical wake models do not capture the flow instabilities that are excited through the blade pitch actuation. In this work, we develop a framework, which accounts for the impacts of the large-scale coherent structures and turbulence on the mean flow, for modeling AWM. The framework uses a triple-decomposition approach for the unsteady flow field and models the mean flow and fine-scale turbulence with a parabolized Reynolds-averaged Navier–Stokes (RANS) system. The wave components are modeled using a simplified spatial linear stability formulation that captures the growth and evolution of the coherent structures. Comparisons with high-fidelity large eddy simulations (LESs) of the turbine wakes showed that this framework was able to capture the additional wake mixing and faster wake recovery in the far-wake regions for both the pulse and helix AWM strategies with minimal computational expense. In the near-wake region, some differences are observed in both the RANS velocity profiles and initial growth of the large-scale structures, which may be due to some simplifying assumptions used in the model.

17 WIND ENERGY↗
DOE OSTI2022

EAST steady-state long pulse H-mode with core-edge integration for CFETR

Recent EAST experiment has successfully demonstrated long pulse steady-state high plasma performance scenario and core-edge integration since the last IAEA in 2018. A discharge with a duration over 60s with β P ~2.0, β N ~1.6, H 98y2 ~1.3 and internal transport barrier on electron temperature channel is obtained with multi-RF power heating and current drive. A higher β N (β N ~1.8, β p ~2.0, H 98y2 ~1.3, n e /n GW ~0.75) with a duration of 20s is achieved by using the modulated neutral beam and multi-RF power, where several normalized parameters are close or even higher than the phase III 1GW scenario of CFETR steady-state. High-Z impurity accumulation in the plasma core is well controlled in a low level by using the on-axis ECH. Modelling shows that the strong diffusion of TEM turbulence in the central region prevents tungsten impurity to accumulate. More recently, EAST has demonstrated compatible core-edge integration discharges in the high βp scenario: high confinement H 98y2 >1.2 with high β P ~2.5/β N ~2.0 and f bs ~50% is sustained with reduced divertor heat flux at high density n e /n GW ~0.7 and moderate q 95 ~6.7. By combining active impurity seeding through radiative divertor feedback control and strike point splitting induced by resonant perturbation coil, the peak heat flux is reduced by 20-30% on the ITER-like tungsten divertor, here a mixture of 50% neon and 50% D2 is applied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical coupling of aerosol emissions, dry removal, and turbulent mixing in the E3SM Atmosphere Model version 1 (EAMv1) – Part 1: Dust budget analyses and the impacts of a revised coupling scheme

Abstract. An earlier study evaluating dust life cycle in the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1) has revealed that the simulated global mean dust lifetime is substantially shorter when higher vertical resolution is used, primarily due to significant strengthening of dust dry removal in source regions. This paper demonstrates that the sequential splitting of aerosol emissions, dry removal, and turbulent mixing in the model's time integration loop, especially the calculation of dry removal after surface emissions and before turbulent mixing, is the primary reason for the vertical resolution sensitivity reported in that earlier study. Based on this reasoning, we propose a revised numerical process coupling scheme that requires the least amount of code changes, in which the surface emissions are applied before turbulent mixing instead of before dry removal. The revised scheme allows newly emitted particles to be transported aloft by turbulence before being removed from the atmosphere, and hence better resembles the dust life cycle in the real world. Sensitivity experiments show that the revised process coupling substantially weakens dry removal and strengthens vertical mixing in dust source regions. It also strengthens the large-scale transport from source to non-source regions, strengthens dry removal outside the source regions, and strengthens wet removal and activation globally. In transient simulations of the years 2000–2009 conducted using 1∘ horizontal grid spacing, 72 vertical layers, and unchanged tuning parameters of emission strength, the revised process coupling leads to a 40 % increase in the global total dust burden and an increase of dust lifetime from 1.8 to 2.5 d in terms of 10-year averages. Weakened dry removal and increased mixing ratios are also seen for other aerosol species that have substantial surface emissions, although the changes in mixing ratio are considerably smaller for the submicron species than for dust and sea salt. Numerical experiments confirm that the revised coupling scheme significantly reduces the strong and non-physical sensitivities of model results to vertical resolution in the original EAMv1. This provides a motivation for adopting the revised scheme in EAM as well as for further improvements on the simple revision presented in this paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NO x = NO + NO 2 ). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NO x is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO 2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

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

The NASA ACTIVATE Mission

The NASA Aerosol Cloud Meteorology Interactions over the Western Atlantic Experiment (ACTIVATE) conducted 162 joint flights with two aircraft over the northwest Atlantic to study aerosol–cloud interactions (ACIs), which represent the largest uncertainty in estimating total anthropogenic radiative forcing. The combination of a high-flying King Air and low-flying HU-25 Falcon, equipped with remote sensing and in situ instruments, characterized trace gases, aerosol particles, clouds, and meteorological variables with data collected nearly simultaneously below, within, and above marine boundary layer (MBL) clouds. Flights spanning warm and cold seasons across 3 years (2020–22) provided a broad range of conditions associated with aerosol particles, cloud properties (including particle size and phase), and meteorology, ideally suited for robust ACI calculations and assessing how well models simulate a wide range of MBL clouds from stratiform to cumulus. ACTIVATE data suggest that drivers of cloud droplet number concentration N d , including aerosol particles and MBL dynamics, vary between winter and summer months with a stronger potential to convert aerosol particles into cloud droplets in winter. Models of varying complexity not only highlight some skills in simulating winter and summer cloud types but also identify challenges that still need to be addressed such as treatment of turbulence, wet scavenging, and mesoscale organization. Remote sensing advances range from new retrieval methods for N d , cloud phase classification, vertically resolved aerosol and cloud condensation nuclei number concentration, and ocean surface wind speed. This work describes these scientific and technological advances along with efforts in outreach and open data science.

aerosol indirect effect↗
DOE OSTI2023

Azimuthal structures and turbulent transport in Penning discharge

Azimuthal structures in cylindrical Penning discharge are studied with a 2D3V radial–azimuthal PIC/MCC model with the axial magnetic field. The discharge is self-consistently supported by ionization due to the axial injection of electrons. It is shown that the steady-state discharge can be supported in two different regimes with different types of observed azimuthal structures. The transition between the regimes is controlled by the mechanism of the energy input to the discharge. In the first regime (low energy of the injected electrons), with the pronounced m = 1 spoke activity, the power input is dominated by the energy absorption due to the radial current and self-consistent electric field. In the other regime (higher energy of the injected electrons), with prevalent small-scale m > 1 spiral structures, and the lower values of the anomalous transport, the total energy deposited to the discharge is lower and is mostly due to the direct input of the kinetic energy from the axial electron beam. Here we show that the large (m = 1) spoke and small-scale structures occur as a result of Simon–Hoh and lower hybrid instabilities driven by the electric field, density gradient, and collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗