Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical modeling”

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 OSTI2026

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗
DOE OSTI2024

VTS Dynamic Model and Experimental Validation for Vibration Suppression

High Q-factor superconducting radio-frequency (SRF) cavities are extremely sensitive to mechanical vibrations, which can cause an uncontrolled shift of the electromagnetic resonance frequency, with detrimental effects on the cavity performance. To avoid this phenomenon (i.e. microphonics) the dynamic behavior of the assembly must be studied, to design vibration suppression systems. In this paper, the numerical model of the Vertical Test Stand (VTS) insert adopted at the Fermi National Accelerator Laboratory for cold tests was studied, in the framework of the experimental characterization of Dark SRF cavities conducted by the SQMS (Superconducting Quantum Materials and Systems) center. Due to the complexity of the assembly, experimental activity is of utmost importance to validate the model. Thus, vibration testing was performed on the actual system, and a good correlation between numerical and experimental results was found. Additionally, two different setups were compared, with and without implementing a preliminary passive vibration suppression system, to assess its effectiveness in reducing the microphonics phenomenon. The obtained results allowed for the design of an improved isolation system, which will be the subject of future experimental testing.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗
DOE OSTI2024

Summary Report of the FY24 DOE Contributions to the GIF VHTR CMVB

The Generation-IV Forum (GIF) Very-High-Temperature Reactor-Computational Methods Validation and Benchmark (VHTR-CMVB) initiative, involving organizations from Korea Atomic Energy Research Institute (KAERI) (South Korea), Institute of Nuclear and New Energy Technology of Tsinghua University (INET) (China), U.S. Department of Energy (DOE) (U.S.), Joint Research Centre (JRC) (Europe), and Japan Atomic Energy Agency (JAEA) (Japan), is dedicated to the verification and validation of tools for High-Temperature Gas-Cooled Reactors (HTGRs) analysis, using data shared by Computational Methods Validation and Benchmark (CMVB) signatories. For FY24, the US DOE CMVB has committed to several critical activities. Under WP1, led by the US, the integration of the High Temperature Gas Cooled Reactor - Pebble-Bed Module (HTR-PM) Phenomena Identification and Ranking Table (PIRT) into the comparative PIRT is progressing, with a new draft of the comparison tables issued earlier this year and currently being utilized by INET for their contribution. Neutronic validation efforts under WP3 include the preparation of the burnup analysis benchmark, preliminary calculations, and the development of reference models and results. In WP2, a validation exercise for hot gas mixing in the lower plenum of HTR-PM is in progress, using experimental data from INET (China) to validate modeling approaches. A model of the experimental facility has been developed using StarCCM+, with initial calculations slated for presentation at the GIF CMVB meeting this fall. Another WP2 activity focuses on validating numerical models for air-cooled Reactor Cavity Cooling System (RCCS) with experimental data from the Wisconsin Madison RCCS facility. A high-fidelity model, developed using NEK-RS, is currently being validated with available data from a low power forced convection test. These efforts are aimed at enhancing and confirming the accuracy of HTGR analysis tools, ensuring their alignment with experimental data and regulatory requirements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Analysis of heat transfer and AuNPs-mediated photo-thermal inactivation of E. coli at varying laser powers using single-phase CFD modeling

In the wake of the COVID-19 pandemics, the demand for innovative and effective methods of bacterial inactivation has become a critical area of research, providing the impetus for this study. The purpose of this research is to analyze the AuNPs-mediated photothermal inactivation of E. coli. Gold nanoparticles irradiated by laser represent a promising technique for combating bacterial infection that combines high-tech and scientific progress. The intermediate aim of the work was to present the calibration of the model with respect to the gold nanorods experiment. The purpose of this work is to study the effect of initial concentration of E. coli bacteria, the design of the chamber and the laser power on heat transfer and inactivation of E. coli bacteria. Using the CFD simulation, the work combines three main concepts. 1. The conversion of laser light to heat has been described by a combination of three distinctive approximations: a- Discrete particle integration to take into account every nanoparticle within the system, b- Rayleigh-Drude approximation to determine the scattering and extinction coefficients and c- Lambert–Beer–Bourger law to describe the decrease in laser intensity across the AuNPs. 2. The contribution of the presence of E. coli bacteria to the thermal and fluid-dynamic fields in the microdevice was modeled by single-phase approach by determining the effective thermophysical properties of the water-bacteria mixture. 3. An approach based on a temperature threshold attained at which bacteria will be inactivated, has been used to predict bacterial response to temperature increases. The comparison of the thermal fields and temporal temperature changes obtained by the CFD simulation with those obtained experimentally confirms the accuracy of the light-heat conversion model derived from the aforementioned approximations. The results show a linear relationship between maximum temperature and variation in laser power over the range studied, which is in line with previous experimental results. It was also found that the temperature inside the microchamber can exceed 55 °C only when a laser power higher than 0.8 W is used, so bacterial inactivation begins. The experimental data allows to determinate the concentration of nanoparticles. This parameter is introduced into the mathematical model obtaining the same number of AuNPs. However, this assumption introduces a certain simplification, as in the mathematical model the distribution of nanoparticles is uniform. This work is directly connected to the use of gold nanoparticles for energy conversion, as well as the field of bacterial inactivation in microfluidic systems such as lab-on-a-chip. Presented mathematical and numerical models can be extended to the entire spectrum of wavelengths with particular use of white light in the inactivation of bacteria. This work represents a significant advancement in the field, as to the best of the authors’ knowledge, it is the first to employ a single-phase computational fluid dynamics (CFD) approach specifically combined with the thermal inactivation of bacteria. Moreover, this research pioneers the use of a numerical simulation to analyze the temperature threshold of photothermal inactivation of E. coli mediated by gold nanorods (AuNRs). The integration of these methodologies offers a new perspective on optimizing bacterial inactivation techniques, making this study a valuable contribution to both computational modeling and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On measuring stimulated photon–photon scattering using multiple ultraintense lasers

Stimulated photon–photon scattering is a predicted consequence of quantum electrodynamics that has yet to be measured directly. Measuring the cross section for stimulated photon–photon scattering is the aim of a flagship experiment for NSF OPAL, a proposed laser user facility with two, 25-PW beamlines. We present optimized experimental designs for achieving this challenging and canonical measurement. A family of experimental geometries is identified that satisfies the momentum- and energy-matching conditions for two selected laser frequency options. Numerical models predict a maximum signal exceeding 1000 scattered photons per shot at the experimental conditions envisaged at NSF OPAL. Experimental requirements on collision geometry, polarization, cotiming and copointing, background suppression, and diagnostic technologies are investigated numerically. These results confirm that a beam cotiming shorter than the pulse duration and control of the copointing on a scale smaller than the shortest laser wavelength are needed to robustly scatter photons on a per-shot basis. Finally, we assess the bounds that a successful execution of this experiment may place on the mass scale of Born–Infeld nonlinear electrodynamics beyond the standard model of physics.

Beyond the Standard Model↗
DOE OSTI2026

Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor

The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.

Geomechanics↗
DOE OSTI2024

Using Parameter Sweep in WaterTAP to Analyze New Water Treatment Technologies

We describe a powerful and generalized parameter sweep tool in this report that was originally developed to analyze the performance of existing and novel water treatment models being developed in WaterTAP. Since WaterTAP is built upon IDAES and Pyomo, the parameter sweep tool can be used to systematically explore and debug the behavior of most Pyomo and IDAES numerical models. In order to enable meaningful analyses, the parameter sweep tool has been designed with the following features: 1) Model flexibility: The parameter sweep tool does not enforce any restrictions on the types of models that can be used with it. As long as a Pyomo model can be solved and the parameter is active and mutable, the tool only needs functions that describe how to run the model, the sweep parameters, and the output quantities of interest. 2) Flexible sampling: The parameter sweep tool has inbuilt functions to generate samples from a random distribution or a multidimensional Euclidean space. Furthermore, the users have to ability to supply samples generated from a tool of their choice. 3) Multiple sweep types: A user can choose from one of 3 types of parameter sweeps depending on their needs. 4) Detailed outputs: Outputs generated by the parameter sweep tool can be stored in detailed H5 file or user-friendly CSV files for post processing. 5) Parallel computing: The parameter sweep supports shared and distributed memory parallel computing to enable the use of high performance computers (HPC) for large-scale analyses. 6) Modular: The parameter sweep tool is self-contained and can easily be integrated within an outer-loop analysis or as desired by the user. 7) Ease of use: The tool is well documented and a simple sweep can be easily executed by following the online documentation in a few lines of code. We demonstrate the use of the parameter sweep tool on a simple water treatment system from the WaterTAP repository and show its parallel scaling performance on an Apple laptop and NREL's Eagle HPC. The parameter sweep tool is actively being used with models currently being developed within WaterTAP and we expect its use to grow beyond it to other IDAES and Pyomo models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Understanding Climate and Extreme Weather Events in the Greater New York Area

The goal of this project is to create a framework for modeling extreme events using climate information for risk assessment in the Greater New York City (NYC) region. The central research question we addressed was how best can we use information from historical climate data and numerical models of the climate system to estimate changing risks of rainfall extremes at specific locations – starting with the NYC region as a testbed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic Modelling of HiveWind Floating Wind Substructure in OpenFAST

Current aero-elastic modes tend to consider the wind turbine floating substructures as rigid. However, the new generation of offshore wind platforms are leaner, less rigid. The OpenFAST code has been updated to consider the flexibility of the floating substructure. This presentation shows how a numerical model of the HiveWind design was created within OpenFAST as well as some verifications performed. When comparing the flexible model to a simplified rigid model (previous state-of-the-art), there are some significant differences in key performance indicators. The new flexibility capability within OpenFAST allows to calculate ultimate and fatigue loads within the floating substructure, allowing for the development of more economical designs.

ENGINEERING,WIND ENERGY↗
DOE OSTI2024

Nanoscale structural correlations in a model cuprate superconductor

Understanding the extent and role of inhomogeneity is a pivotal challenge in the physics of cuprate superconductors. While it is known that structural and electronic inhomogeneity is prevalent in the cuprates, it has proven difficult to disentangle compound-specific features from universally relevant effects. Here, in this study, we combine advanced neutron and x-ray diffuse scattering with numerical modeling to obtain insight into bulk structural correlations in HgBa 2 ⁢ CuO 4+δ . This cuprate exhibits a high optimal transition temperature of nearly 100 K, pristine charge-transport behavior, and a simple average crystal structure without long-range structural instabilities, and is therefore uniquely suited for investigations of intrinsic inhomogeneity. We uncover diffuse reciprocal-space patterns that correspond to prominent nanoscale correlations of atomic displacements perpendicular to the CuO 2 planes. The real-space nature of the correlations is revealed through three-dimensional pair distribution function analysis and complementary numerical refinement. We find that relative displacements of ionic and CuO 2 layers play a crucial role, and that the structural inhomogeneity is not directly caused by the presence of conventional point defects. The observed correlations are therefore intrinsic to HgBa 2 ⁢ CuO 4+δ , and thus likely important for the physics of cuprates more broadly. It is possible that the structural correlations are closely related to the unusual superconducting fluctuations and Mott-localization in these complex oxides. As advances in scattering techniques yield increasingly comprehensive data, the experimental and analysis tools developed here for large volumes of diffuse scattering data can be expected to aid future investigations of a wide range of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Debunking common myths in coastal circulation modeling

Despite tremendous progress in algorithm development, computational efficiency and transition into operations over the past two decades, coastal modeling still lacks scientific rigor due to proliferation of many ‘gray’ areas related to various modeling choices made by modelers. Here, in this paper, we propose some guiding principles for the modeling community to improve performance, and we also debunk commonly held myths that make the coastal modeling lack rigor. Using our own experience in developing seamless cross-scale unstructured-grid based models for the past two decades, we describe in unprecedented detail the end-to-end modeling process (i.e., from digital elevation models (DEMs) to mesh generation to post analysis), and demonstrate that defensible modeling is within reach for any end user by following three guiding principles: (1) Bathymetry is a first order forcing in coastal domains and thus should be respected in all aspects of modeling; (2) Oceanographic processes are driven across multiple spatial scales and so models should enable appropriate resolution as needed; and (3) Model assessment should focus on physical processes. Through qualitative and quantitative model assessments, we demonstrate the fundamental role played by bathymetry/topography as embedded in DEMs in making the results defensible, which is unfortunately glossed over in many modeling studies. Focusing on process-based assessment simplifies the calibration process. A major conclusion of this work is that model developers and operators should maximize the scientific rigor for in silico oceanography by avoiding some common pitfalls that rely on error compensation at the expense of representation of physical system processes. We present some best practice procedures for defensive and trustworthy numerical modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measurement of Close-in Ground Motion from an Underground Chemical Explosion

Understanding the geophysical response near an underground explosion is crucial for generating insights into the source and emplacement conditions that produce distinct observations in monitoring scenarios occurring at greater distances. Recently, Shot A of the Low Yield Nuclear Monitoring (LYNM) Physics Experiment 1 (PE1) series was conducted at the Nevada National Security Site to provide ground truth for subsurface explosion signal models. This experiment resulted in measuring near-source ground motion at distances ranging from 70 to 1000 m/kt with a 99% success rate, yielding high-fidelity knowledge of the near-field response that can serve as benchmarks for future numerical modeling and experiment planning. However, technical challenges exist in observing near-source phenomena while safeguarding sensitive data acquisition components from the detrimental effects of ground motion in the subsurface. This report outlines tools and techniques to address challenges associated with observing near-source accelerations and within the tunnel drift of the PE1 test bed. Additionally, we describe key systems designed with both modern advancements and legacy guidance to maximize the collection of high-quality ground motion data, which may be applied to constitutive and computational models, leading to new or improved understanding of the near- and far-field signals produced by underground explosions.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗
DOE OSTI2024

Multiobjective Constrained Symbolic Regression for Predictive Modeling of Material Creep Behavior

When creep testing is repeated on samples of the same alloy under the same parametric conditions (i.e., stress and temperature), the resulting strain/time curves can vary from each other considerably as shown in Figure 1 [1]. The time required to creep test a material to rupture can extend to the order of years. Because of this, a numerical model that can quickly analyze the incomplete results of an ongoing experiment to predict 1) the incomplete portion of the strain/time curve leading up to the rupture point and 2) the rupture point itself would be of great utility to the materials community. Such a model has the potential to save 1) the time required to finish running the experiment to rupture 2) the associated monetary cost of finishing said experiment. Furthermore, it would be advantageous if the predictive model could give a parametric function modeling strain/time curves for material scientists to investigate the impact of the temperature and stress parameters on the resulting creep behavior. This work introduces a piecewise symbolic regression algorithm to predict the remainder of the strain/time curve. Preliminary results show good model performance.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗
DOE OSTI2024

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗