Results for “Average model” — Page 53

NASA NTRS1982

A model study of Fuego volcanic aerosol dispersion in the lower stratosphere

A zonally averaged time-dependent primitive equation model is used to simulate the dispersion of both a carbon 14 injection and the volcanic aerosol from the 1974 Fuego eruption. It is noted that both injections occurred at low latitudes to midlatitudes in the Northern Hemisphere. The eddy flux terms, which account for the major portion of the transport in the lower stratosphere of this model, are specified in a manner similar to that of Harwood and Pyle (1975). Comparisons with data underline the ability of the model to simulate the vertical character of the tracer while maintaining reasonable meridional transport times. For the aerosol study, the simulated 1/e decay time at 37 deg N and for the 16- to 21-km altitude region is 9 months, whereas lidar measurements at the same latitude give a decay time of 8 months. The simulated vertical width at half-maximum for the aerosol tracer at 37 deg N and 19 deg N and for 6 months after the event possesses values of 5.0 km and 3.6 km, respectively, whereas the observed lidar values are 4.4 km and 3.0 km, respectively. The tracer transport to the Southern Hemisphere also is in qualitative agreement with the limited data that are available.

Remsberg, E. E.↗

NASA NTRS1984

Simple performance model for ring and line cusp ion thrusters

A model is developed for high magnetic field strength cusped thrusters which results in a single equation describing the thruster performance. Such a simple result is made possible by formulating the model in terms of the average energy expended in producing ions in the discharge plasma and the fraction of these ions that are extracted into the beam. A key feature of the model is that it allows one to calculate the average plasma ion energy cost as a function of the propellant utilization knowing only thruster geometric design parameters and the propellant gas. The model is exercised to indicate the effects on performance resulting from changes in the extracted ion fraction, the quality of the containment of primary electrons, the propellant mass flow rate, the propellant gas and the transparency of the accelerator system to neutral atoms.

Brophy, J. R.↗

DOE OSTI2022

Reionization time of the Local Group and Local-Group-like halo pairs

ABSTRACT Patchy cosmic reionization resulted in the ionizing UV background asynchronous rise across the Universe. The latter might have left imprints visible in present-day observations. Several numerical simulation-based studies show correlations between the reionization time and overdensities and object masses today. To remove the mass from the study, as it may not be the sole important parameter, this paper focuses solely on the properties of paired haloes within the same mass range as the Milky Way. For this purpose, it uses CoDaII, a fully coupled radiation hydrodynamics reionization simulation of the local Universe. This simulation holds a halo pair representing the Local Group, in addition to other pairs, sharing similar mass, mass ratio, distance separation, and isolation criteria but in other environments, alongside isolated haloes within the same mass range. Investigations of the paired halo reionization histories reveal a wide diversity although always inside-out, given our reionization model. Within this model, haloes in a close pair tend to be reionized at the same time but being in a pair does not bring to an earlier time their mean reionization. The only significant trend is found between the total energy at z = 0 of the pairs and their mean reionization time: Pairs with the smallest total energy (bound) are reionized up to 50 Myr earlier than others (unbound). Above all, this study reveals the variety of reionization histories undergone by halo pairs similar to the Local Group, that of the Local Group being far from an average one. In our model, its reionization time is ∼625 Myr against 660 ± 4 Myr (z ∼ 8.25 against 7.87 ± 0.02) on average.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Computationally Efficient Multiscale Neural Networks Applied to Fluid Flow in Complex 3D Porous Media

Abstract The permeability of complex porous materials is of interest to many engineering disciplines. This quantity can be obtained via direct flow simulation, which provides the most accurate results, but is very computationally expensive. In particular, the simulation convergence time scales poorly as the simulation domains become less porous or more heterogeneous. Semi-analytical models that rely on averaged structural properties (i.e., porosity and tortuosity) have been proposed, but these features only partly summarize the domain, resulting in limited applicability. On the other hand, data-driven machine learning approaches have shown great promise for building more general models by virtue of accounting for the spatial arrangement of the domains’ solid boundaries. However, prior approaches building on the convolutional neural network (ConvNet) literature concerning 2D image recognition problems do not scale well to the large 3D domains required to obtain a representative elementary volume (REV). As such, most prior work focused on homogeneous samples, where a small REV entails that the global nature of fluid flow could be mostly neglected, and accordingly, the memory bottleneck of addressing 3D domains with ConvNets was side-stepped. Therefore, important geometries such as fractures and vuggy domains could not be modeled properly. In this work, we address this limitation with a general multiscale deep learning model that is able to learn from porous media simulation data. By using a coupled set of neural networks that view the domain on different scales, we enable the evaluation of large ( $$>512^3$$ > 512 3 ) images in approximately one second on a single graphics processing unit. This model architecture opens up the possibility of modeling domain sizes that would not be feasible using traditional direct simulation tools on a desktop computer. We validate our method with a laminar fluid flow case using vuggy samples and fractures. As a result of viewing the entire domain at once, our model is able to perform accurate prediction on domains exhibiting a large degree of heterogeneity. We expect the methodology to be applicable to many other transport problems where complex geometries play a central role.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Validation of RANS Models for Predicting Bypass Flow in Prismatic VHTRs

In prismatic Very-High Temperature Reactors (VHTRs), gaps occur between neighboring fuel assemblies referred to as bypass gaps. This work aims to measure the flow distribution between the bypass gaps between neighboring assemblies and the coolant channels within the assemblies. Computational fluid dynamics (CFD) was done using Reynolds-Averaged Navier Stokes (RANS) modeling in STAR-CCM+ to assess the ability of RANS models to accurately capture this phenomena. The RANS results were compared to experimental worked performed by Seoul National University (SNU) and Korea Atomic Energy Research Institute (KAERI) for validation. The CFD results showed good comparison with the experimental data available from SNU and flow split was characterized for multiple configurations, inlet mass flow rates, and bypass gap sizes.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1996

Enhanced Water Vapor Absorption within Tropospheric Clouds: A Partial Explanation for Anomalous Absorption

Comparisons between solar flux measurements and predictions obtained from theoretical radiative transfer models indicate that most of these models underestimate the globally averaged solar energy absorbed by cloudy atmospheres by up to 25Wm&sup-2;.The origin of this anomalous absorption has not yet been established, but it has been attributed to a variety of sources including oversimplified or missing physical processes in the existing models, uncertainties in the input data, and even measurement errors. We used a sophisticated atmospheric radiative transfer model to provide improved constraints on the physical processes that contribute to the absorption of solar radiation by Earth's atmosphere. The results are described herein.

solar flux measurement water vapor absorption atmo↗

NASA NTRS2015

Uncertainty Quantification of Turbulence Model Closure Coefficients for Transonic Wall-Bounded Flows

The goal of this work was to quantify the uncertainty and sensitivity of commonly used turbulence models in Reynolds-Averaged Navier-Stokes codes due to uncertainty in the values of closure coefficients for transonic, wall-bounded flows and to rank the contribution of each coefficient to uncertainty in various output flow quantities of interest. Specifically, uncertainty quantification of turbulence model closure coefficients was performed for transonic flow over an axisymmetric bump at zero degrees angle of attack and the RAE 2822 transonic airfoil at a lift coefficient of 0.744. Three turbulence models were considered: the Spalart-Allmaras Model, Wilcox (2006) k-w Model, and the Menter Shear-Stress Trans- port Model. The FUN3D code developed by NASA Langley Research Center was used as the flow solver. The uncertainty quantification analysis employed stochastic expansions based on non-intrusive polynomial chaos as an efficient means of uncertainty propagation. Several integrated and point-quantities are considered as uncertain outputs for both CFD problems. All closure coefficients were treated as epistemic uncertain variables represented with intervals. Sobol indices were used to rank the relative contributions of each closure coefficient to the total uncertainty in the output quantities of interest. This study identified a number of closure coefficients for each turbulence model for which more information will reduce the amount of uncertainty in the output significantly for transonic, wall-bounded flows.

Schaefer, John↗

NASA NTRS2018

Validation of a CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in a Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

NASA NTRS2018

Validation of A CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in A Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011 [1]. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. [2]. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

DOE OSTI2021

Mitigating the Sign Problem through Basis Rotations

Quantum Monte Carlo simulations of quantum many-body systems are plagued by the Fermion sign problem. The computational complexity of simulating Fermions scales exponentially in the projection time β and system size. The sign problem is basis dependent and an improved basis, for fixed errors, leads to exponentially quicker simulations. Herein, we show how to use sign-free quantum Monte Carlo simulations to optimize over the choice of basis on large two-dimensional systems. We numerically illustrate these techniques decreasing the “badness” of the sign problem by optimizing over single-particle basis rotations on one- and two-dimensional Hubbard systems. We find a generic rotation which improves the average sign of the Hubbard model for a wide range of U and densities for L×4 systems. In one example improvement, the average sign (and hence simulation cost at fixed accuracy) for the 16×4 Hubbard model at U/t=4 and n=0.75 increases by exp[8.64(6)β]. For typical projection times of β ≳ 100, this accelerates such simulation by many orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRSDate not supplied

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗

NASA NTRSDate not supplied

Computational Modeling of Two Mars Powered Descent Vehicle Concepts Tested in the Langley Unitary Plan Wind Tunnel

Future human Mars missions will require powered descent starting at supersonic conditions, something which has never been done before at Mars. Computational powered descent flowfield simulations have been completed at full-scale Mars conditions, but the available ground test data are not suitable for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles. Testing was conducted in the NASA Langley Unitary Plan Wind Tunnel in order to investigate the aerodynamic interference of sub-scale versions of two Mars powered descent vehicle concepts at supersonic Mach numbers (2.4 and 3.5): a model based on a blunt hypersonic inflatable aerodynamic decelerator (HIAD) and the second representing a more slender rigid vehicle with body flaps (CobraMRV). This paper covers computational flowfield predictions completed at wind tunnel conditions and comparisons to the test data. On the blowing HIAD models, the time-averaged pressure inboard of the nozzles was generally well-predicted, especially if the nozzles are canted outward, when the nozzles are located further from the nose. At intermediate CobraMRV thrust coefficients, CFD pressures are more accurately predicted than they are for the HIAD models, largely due to the nozzle locations and pointing directions. Overall, the CFD pressure coefficients were predicted within 0.2 of the steady pressure measurements for all blowing models, with smaller discrepancies at higher HIAD thrust, and larger discrepancies at lower HIAD CobraMRV thrust. All HIAD models were predicted to have a gradually decreasing axial force coefficient as the total thrust increases, in agreement with available pressure sensitive paint data. On models with canted nozzles or with nozzles further from the nose, the force coefficient was slightly higher for a given thrust. On the CobraMRV model, the CFD also shows consistent results between solvers and follows trends revealed in the data; the aerodynamic force coefficient remains near the non-blowing value at a tunnel Mach number of 2.4 regardless of thrust, and increases above that level at a Mach number of 3.5, consistent with the discrete pressure data. CFD analysis at tunnel and flight conditions will continue as flight system designs concepts mature.

Supersonic Retropropulsion↗

DOE OSTI2025

Active Learning Surrogates for Integrating Electron Microscopy and Computational Insights from Simulations in Autonomous Experiments

Artificial Intelligence (AI) combined with simulations and experiments has great potential to accelerate scientific discovery across technology and pharmaceuticals. However, the gap between simulations and experiments is challenging due to disparities in time and scale, making it difficult to estimate properties like energy and electronic states from experiments, and to provide feedback based on theoretical insights.Our research addresses the challenge by developing unique deep kernel based surrogate models that learns from microscopic images, mapping structural features to energy differences from defect formation. We start with full-training using simulated images to determine optimal settings, establishing a baseline for active learning. Using these settings from the baseline, active learning is trained, and predicts structures along simulation trajectories based on uncertainty and energetic stability, thus reducing data requirements, simulation time and computational costs. The results demonstrate that the model achieves a low average error margin of approximately 0.03 meV, indicating good performance. To enhance feature extraction and reconstruction capabilities, we developed an autoencoder-decoder as additional surrogate to create latent space to capture essential features, enabling precise comparisons between simulations and experiments. The results from this model achieved a reconstruction loss of around 0.2 and accurately reconstructed molecular structures.Overall, this work advances the steering of experiments through computational simulations by employing a surrogate models that actively predicts the trajectories of structural evolution, achieving time-to-solution comparable to experimental measurements.

Saranathan, Gayathri [Hewlett-Packard]↗

DOE OSTI2022

Use Attainability Analysis Aquatic Life Uses for Perennial Reach of Sandia Canyon [Slides]

The following information will be presented during the February 23, 2022 Accord Pueblos Technical Meeting. The purpose of the study is to determine if natural thermal conditions are preventing the attainment of Coldwater Aquatic Life Use in the perennial reach of the upper Sandia Canyon Assessment Unit (Sandia AU). The New Mexico Water Quality Standards allow for a change in the designated use if a Use Attainability Analysis (UAA) demonstrates that the use is not attainable due to one or more of six factors listed in 40 CFR 131.10(g), including naturally occurring pollutant concentrations. The Sandia AU is located in a perennial reach of upper Sandia Canyon between Sigma Canyon and NPDES Outfall 001. The classified Segment 20.6.4.126 NMAC comprises perennial waters within Los Alamos National Laboratory boundaries and includes the Sandia AU. Persistent surface flows originate from NPDES permitted effluent releases. These releases have occurred since the early 1950’s and continue today. The UAA examines several lines of evidence. NMED’s Air-Water Temperature Correlation (AWTC) model is used for identifying appropriate stream classifications and attainable aquatic life use subcategories. The model correlates between July average air temperatures (ATEMP) and maximum weekly average stream temperatures. Air temperatures are obtained from PRISM and LANL Meteorological Towers and used to derive ATEMP. Thermographs were placed in the Sandia AU to obtain measured stream temperatures during summer months when stream temperatures are the highest. Measured data is used with AWTC modeled data to determine if the Sandia AU is meeting its natural air temperature driven thermal condition.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Assessment of the Impact of Realistic Sensor Physics and the Integration of Ex-Core Sensors on Reactor Power Synthesis

In the work documented in this report, a weighting function–based core power synthesis method was applied to multiple Monte Carlo N-Particle (MCNP) reactor models, which are informed based on simulated self-powered neutron detector (SPND) responses. The weighting function method used has been coined the point-based iterative (PBI) method. The goal of this application is to assess the impact of considering realistic sensor physics in the generation of the simulated SPND outputs as well as to consider how the synthesis is impacted based on the inclusion of ex-core detectors in the model. The NuScale small modular reactor (SMR) and Westinghouse AP1000 pressurized water reactor (PWR) are the models that served as the testbeds for the assessment of realistic sensor physics; this was achieved by using Geant4 SPND models in comparison with analytical models, such that the effect of electron transport in realistic SPND geometries in the Geant4 model can be understood in terms of synthesis error and convergence time. The comparison was considered for fuel burnup–induced perturbations, for a range of sensor string densities and synthesized power distribution axial fidelities. The Texas A&M Testing, Research, Isotopes, General Atomics Reactor (TAMU TRIGA) reactor MCNP model was used to assess the impact of ex-core sensors; this was done by performing synthesis with and without the ex-core detectors and by quantifying the synthesis error and number of iterations associated with Gaussian-type perturbations in many locations in the core. The TAMU TRIGA model was particularly pertinent for this study because of the interest in future experimental tests with SPNDs in this reactor, as well as the ease of modifying the MCNP model to include ex-core detectors with heterogeneously described response functions. Results from the comparison between the Geant4 and analytical SPND models indicate that similar average and maximum synthesis errors were obtained for burnup-induced perturbations in both the NuScale SMR and the AP1000. This was true for a range of sensor string densities and axial fidelities. However, there were marked differences between both the Geant4 and analytically informed models in terms of the iterations required to converge on the synthesized power distribution. Namely, the Geant4-informed models tended to lead to fewer iterations, except for a few sensor–core configurations that had particularly numerous iterations. Results from the ex-core sensor assessment with the TAMU TRIGA model indicate that the inclusion of ex-core sensors drastically reduces the synthesis error of Gaussian-type perturbations close to the edge of the core, and it slightly reduces synthesis errors for perturbations closer to the center of the core. This was achieved with a minimal increase in computational cost—that is, the number of iterations required for convergence. The errors were identified to be in the same location as the perturbation in the core, indicating that the methodology remains robust for unperturbed regions of the core. A secondary result from this study with the TAMU TRIGA was yielded by analysis of the neutron flux levels in the in-core and ex-core sensor locations of the core; these flux levels indicate that SPNDs could be used as both in-core and ex-core sensors, so long as the emitter material is sensitive to thermal neutrons. The results from these studies provide a quantitative understanding of the importance of considering realistic sensor physics and including ex-core sensors to perform accurate and timely power distribution synthesis of a reactor core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗