Search NASA⌕ Search

SEARCH · Search NASA

Results for “stochastic control”

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 811 records · Page 45

A Physics-Informed Reinforcement Learning Framework for Economic-Thermal Co-Optimization of Crypto Mining Data Centers: Preprint

The rapid expansion of cryptocurrency mining has created a new class of high-density data centers characterized by extreme thermal flux and high sensitivity to volatile economic markets. Traditional thermal management strategies, typically reliant on rule-based control, maintain static setpoints that fail to account for fluctuating electricity prices and cryptocurrency values - factors critical to mining profitability. To address this, we present a physics-informed reinforcement learning (PIRL) framework for economic-thermal co-optimization in crypto mining data centers. This framework consists of a proximal policy optimization (PPO) agent, a virtual testbed powered by high-fidelity physics-based models, and an interactive frontend dashboard. The PPO agent is trained using the virtual testbed and strict hardware safety limits. This physics-informed approach allows the agent to learn a stochastic policy that dynamically balances mining revenue against operational costs by co-optimizing HVAC cooling setpoints and IT computational hashrate. The simulation results demonstrate that the integrated framework achieved an 8.62% increase in net operational profit compared to traditional baseline strategies while strictly adhering to safety-critical temperature constraints (coolant supply temperature < 32 degrees C). This work provides a scalable template for the deployment of reinforcement learning in mission critical facilities where economic volatility and physical safety must be managed simultaneously.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Operational Resilience Metric to Evaluate Inertia and Inverter-based Generation on the Grid

In an effort to reduce carbon emissions and curtail the effects of climate change there has been considerable effort to increase the penetration of inverter-based renewable energy sources. The adoption of renewable generation over conventional inertia-based generation sources is forming considerable challenges for the operation and stability of the power system. The power system has been designed around generation units characterized by high inertia and primary frequency response (PFR), allowing them to respond to disturbances such as faults or generators tripping off line. In contrast, the modern inverter-based assets are characterized by no contribution to inertia and they typically provide stochastic generation at their maximum output, thus having no contribution to PFR during low-frequency disturbance events. Because of this, inertia in power systems is reducing, resulting in a faster rate of frequency change after a disturbance occurs. As more inverter-based generation units are added to the grid it is important to understand the stability of the system and the size of disturbance a system is capable of withstanding. This paper presents a resilience metric that evaluates the maximum size of disturbance a systems can withstand based on the system inertia and the primary frequency control of inverter and inertia-based generation. The results are shown visually and are based on the real-time operation of generation units and their characteristics such as latency, ramp rates, and energy constraints. It is demonstrated that the real-time positioning or bias of the generating units has an effect of the size of disturbance that a system can withstand, i.e. its resilience. It is expected that this type of analysis can help operators increase the resilience of power systems in the future.

13 HYDRO ENERGY↗

Ice Supersaturation Variability in Cirrus Clouds: Role of Vertical Wind Speeds and Deposition Coefficients

Abstract Aircraft measurements reveal ice supersaturation statistics in cirrus (ISSs) with broad maxima around ice saturation and pronounced variance. In this study, processes shaping ISSs in midlatitude and tropical upper tropospheric conditions are systematically investigated. Water vapor deposition and sublimation of size‐resolved ice crystal populations are simulated in an air parcel framework. Mesoscale temperature fluctuations (MTFs) due to gravity waves force the temporal evolution of supersaturation. Various levels of background wave forcing and cirrus thickness are distinguished in stochastic ensemble simulations. Kinetic limitations to ice mass growth are brought about by supersaturation‐dependent deposition coefficients that represent efficient and inefficient growth modes as a function of ice crystal size. The simulations identify a wide range of deposition coefficients in cirrus, but most values stay above 0.01 such that kinetic limitations to water uptake remain moderate. Supersaturation quenching times are long, typically 0.5–2 hr. The wave forcing thus causes a remarkably large variability in ISSs and cirrus microphysical properties except in the thickest cirrus, producing ensemble‐mean ISSs in line with in‐situ measurements. ISS variance is controlled by MTFs and increases with decreasing cirrus integral radii. In comparison, the impact of ice crystal growth rates on ISSs is small. These results contribute to efforts directed at identifying and solving issues associated with ice‐supersaturated areas and non‐equilibrium cirrus physics in global models.

Meteorology & Atmospheric Sciences↗

NASA Tech Briefs, March 2014

Topics include: Data Fusion for Global Estimation of Forest Characteristics From Sparse Lidar Data; Debris and Ice Mapping Analysis Tool - Database; Data Acquisition and Processing Software - DAPS; Metal-Assisted Fabrication of Biodegradable Porous Silicon Nanostructures; Post-Growth, In Situ Adhesion of Carbon Nanotubes to a Substrate for Robust CNT Cathodes; Integrated PEMFC Flow Field Design for Gravity-Independent Passive Water Removal; Thermal Mechanical Preparation of Glass Spheres; Mechanistic-Based Multiaxial-Stochastic-Strength Model for Transversely-Isotropic Brittle Materials; Methods for Mitigating Space Radiation Effects, Fault Detection and Correction, and Processing Sensor Data; Compact Ka-Band Antenna Feed with Double Circularly Polarized Capability; Dual-Leadframe Transient Liquid Phase Bonded Power Semiconductor Module Assembly and Bonding Process; Quad First Stage Processor: A Four-Channel Digitizer and Digital Beam-Forming Processor; Protective Sleeve for a Pyrotechnic Reefing Line Cutter; Metabolic Heat Regenerated Temperature Swing Adsorption; CubeSat Deployable Log Periodic Dipole Array; Re-entry Vehicle Shape for Enhanced Performance; NanoRacks-Scale MEMS Gas Chromatograph System; Variable Camber Aerodynamic Control Surfaces and Active Wing Shaping Control; Spacecraft Line-of-Sight Stabilization Using LWIR Earth Signature; Technique for Finding Retro-Reflectors in Flash LIDAR Imagery; Novel Hemispherical Dynamic Camera for EVAs; 360 deg Visual Detection and Object Tracking on an Autonomous Surface Vehicle; Simulation of Charge Carrier Mobility in Conducting Polymers; Observational Data Formatter Using CMOR for CMIP5; Propellant Loading Physics Model for Fault Detection Isolation and Recovery; Probabilistic Guidance for Swarms of Autonomous Agents; Reducing Drift in Stereo Visual Odometry; Future Air-Traffic Management Concepts Evaluation Tool; Examination and A Priori Analysis of a Direct Numerical Simulation Database for High-Pressure Turbulent Flows; and Resource-Constrained Application of Support Vector Machines to Imagery.

Source record↗

Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes

Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. Here, to address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.

organic↗

A General Framework for Bounding Approximate Dynamic Programming Schemes

For years, there has been interest in approximation methods for solving dynamic programming problems, because of the inherent complexity in computing optimal solutions characterized by Bellman’s principle of optimality. A wide range of approximate dynamic programming (ADP) methods now exists. It is of great interest to guarantee that the performance of an ADP scheme be at least some known fraction, say ß , of optimal. This letter introduces a general approach to bounding the performance of ADP methods, in this sense, in the stochastic setting. The approach is based on new results for bounding greedy solutions in string optimization problems, where one has to choose a string (ordered set) of actions to maximize an objective function. This bounding technique is inspired by submodularity theory, but submodularity is not required for establishing bounds. Instead, the bounding is based on quantifying certain notions of curvature of string functions; the smaller the curvatures the better the bound. The key insight is that any ADP scheme is a greedy scheme for some surrogate string objective function that coincides in its optimal solution and value with those of the original optimal control problem. The ADP scheme then yields to the bounding technique mentioned above, and the curvatures of the surrogate objective determine the value ß of the bound. The surrogate objective and its curvatures depend on the specific ADP.

discrete event systems↗

Megadroughts in Southwestern North America in ECHO-G Millennial Simulations and Their Comparison to Proxy Drought Reconstructions

Simulated hydroclimate variability in millennium-length forced transient and control simulations from the ECHAM and the global Hamburg Ocean Primitive Equation (ECHO-G) coupled atmosphere-ocean general circulation model (AOGCM) is analyzed and compared to 1000 years of reconstructed Palmer drought severity index (PDSI) variability from the North American Drought Atlas (NADA). The ability of the model to simulate megadroughts in the North American southwest is evaluated. (NASW: 25deg42.5degN, 125deg-105degW). Megadroughts in the ECHO-G AOGCM are found to be similar in duration and magnitude to those estimated from the NADA. The droughts in the forced simulation are not, however, temporally synchronous with those in the paleoclimate record, nor are there significant differences between the drought features simulated in the forced and control runs. These results indicate that model-simulated megadroughts can result from internal variability of the modeled climate system rather than as a response to changes in exogenous forcings. Although the ECHO-G AOGCM is capable of simulating megadroughts through persistent La Nina-like conditions in the tropical Pacific, other mechanisms can produce similarly extreme NASW moisture anomalies in the model. In particular, the lack of low-frequency coherence between NASW soil moisture and simulated modes of climate variability like the El Nino-Southern Oscillation, Pacific decadal oscillation, and Atlantic multidecadal oscillation during identified drought periods suggests that stochastic atmospheric variability can contribute significantly to the occurrence of simulated megadroughts in the NASW. These findings indicate that either an expanded paradigm is needed to understand multidecadal hydroclimate variability in the NASW or AOGCMs may incorrectly simulate the strength and/or dynamics of the connection between NASW hydroclimate variability and the tropical Pacific.

variability↗

The role of the droplet interface in controlling the multiphase oxidation of thiosulfate by ozone

Predicting reaction kinetics in aqueous microdroplets, including aerosols and cloud droplets, is challenging due to the probability that the underlying reaction mechanism can occur both at the surface and in the interior of the droplet. Additionally, few studies directly measure the surface activities of doubly charged anions, despite their prevalence in the atmosphere. Here, deep-UV second harmonic generation spectroscopy is used to probe surface affinities of the doubly charged anions thiosulfate, sulfate, and sulfite, key species in the thiosulfate ozonation reaction mechanism. Thiosulfate has an appreciable surface affinity with a measured Gibbs free energy of adsorption of -7.3 ± 2.5 kJ mol -1 in neutral solution, while sulfate and sulfite exhibit negligible surface propensity. The Gibbs free energy is combined with data from liquid flat jet ambient pressure X-ray photoelectron spectroscopy to constrain the concentration of thiosulfate at the surface in our model. Stochastic kinetic simulations leveraging these novel measurements show that the primary reaction between thiosulfate and ozone occurs at the interface and in the bulk, with the contribution of the interface decreasing from ~65% at pH 5 to ~45% at pH 13. Additionally, sulfate, the major product of thiosulfate ozonation and an important species in atmospheric processes, can be produced by two different pathways at pH 5, one with a contribution from the interface of >70% and the other occurring predominantly in the bulk (>98%). The observations in this work have implications for mining wastewater remediation, atmospheric chemistry, and understanding other complex reaction mechanisms in multiphase environments. Future interfacial or microdroplet/aerosol chemistry studies should carefully consider the role of both surface and bulk chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Derivation of heterogeneous material distributions and their sensitivity to HM-coupled two-phase flow models exemplified with the LASGIT experiment

Abstract Advective gas transport in bentonite, a possible buffer material in repositories for radioactive materials, is difficult to simulate in numerical continuum models, partly due to the complicated microstructure of bentonite. To generate reliable models of repositories nevertheless, spatially distributed heterogeneous material properties can be used to allow localization of gas flow. In this study, a pore-size-dependent stochastic approach of the gas entry pressure is derived from Mercury Intrusion Porosimetry, which is used to replicate measurements from the LASGIT experiment. In addition, three benchmark tests are simulated to investigate the dependence of heterogeneous distributions of material properties on the mesh discretization, the temporal dependence, and the coupling between the processes influenced by the heterogeneous parameters. The numerical modeling results of the LASGIT experiment show that the onset of gas flow into the system and the subsequent increase in pressure and stress can be well reproduced using heterogeneous distributions. Compared to a model with homogeneous material properties, heterogeneous distributions may allow the generation of dilatancy-controlled microfractures—an important feature with regard to the advective gas flow in bentonites. However, it can be observed that the heterogeneous distributions in LASGIT are less significant, as technical gaps or differences in material types could have a greater impact.

Environmental Sciences & Ecology↗

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗

Visual Performance Challenges to Low-Frequency Perturbations After Long-Duration Space Flight, and Countermeasure Development

Astronauts experience sensorimotor disturbances after long-duration space flight. After a water landing, crewmembers may need to egress the vehicle within a few minutes for safety and operational reasons in various sea state conditions. Exposure to even low-frequency motions induced by sea conditions surrounding a vessel can cause significant motor control problems affecting critical functions. The first objective of this study was to document human visual performance during simulated wave motion below 2.0 Hz. We examined the changes in accuracy and reaction time when subjects performed a visual target acquisition task in which the location of the target was offset vertically during horizontal rotation at an oscillating frequency of 0.8 Hz. The main finding was that both accuracy and reaction time varied as a function of target location, with greater performance decrements occurring when vertical targets were acquired at perturbing frequencies of 0.8 Hz in the horizontal plane. A second objective was to develop a countermeasure, base d on stochastic resonance (SR), to enhance sensorimotor capabilities with the aim of facilitating rapid adaptation to gravitational transitions after long-duration space flight. SR is a mechanism by which noise can enhance the response of neural systems to relevant sensory signals. Recent studies have shown that applying imperceptible stochastic electrical stimulation to the vestibular system (SVS) significantly improved balance and oculomotor responses. This study examined the effectiveness of SVS on improving balance performance. Subjects performed a standard balance task while bipolar SVS was applied to the vestibular system using constant current stimulation through electrodes placed over the mastoid process. The main finding of this study was that balance performance with the application of SR showed significant improvement in the range of 10%-25%. Ultimately an SR-based countermeasure might be fielded either as preflight training to enhance adaptability, or as a miniature patch-type stimulator worn post flight.

Mulavara, Ajitkumar↗

A control-inspired approach for energy transition planning under uncertainty

As the global carbon footprint continues to grow, many countries are implementing carbon emission reduction policies which have incentivized the expansion of low-carbon and renewable technologies. However, the speed and scale of deployment falls short of that needed to meet climate goals. Energy system models serve as key tools for guiding investment decisions and helping policymakers evaluate the effects of various policies on the development of an energy system. This study focuses on the energy system of the United States and builds upon prior work by incorporating more geographic granularity to account for the trade of commodities and addresses transmission congestion through electricity price adjustments. Furthermore, real-world characteristics, such as delays in constructing new liquid fuel production and electricity generation facilities, are integrated using a sequential decision-making approach that better reflects how decisions can be updated as uncertainties unfold. Results demonstrate that stochastic programming combined with sequential decision-making produces energy transition pathways that are robust to multiple uncertain futures. Additionally, considering real-world characteristics significantly impacts the deployment of renewable technologies and the ability to meet carbon emission reduction goals while also reliably meeting demand. These findings highlight the importance of accounting for uncertainty and real-world characteristics to avoid overly optimistic projections in energy system planning.

energy systems↗

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗

Rapid Optimization of Total Variation with Applications in Imaging, Additive Manufacturing, and Qualification

Total Variation optimization penalizes the gradient of a control variable or state. While this work focuses on image processing in particular, it has also found applications in inverse problems and topology optimization. In image processing, the goal is to maintain faithfulness to the original image while denoising and/or deblurring. Additionally, bilevel optimization over the spatially varying regularization weights can illuminate interfaces such as damage regions and other anomalies. We will address two fundamental challenges with TV-optimization: (i) the typical slow convergence of existing TV-optimization methods, and (ii) the selection of spatially varying TV parameters to promote interface detection. Additionally, we will apply such techniques to image data collected in additive manufacturing. In said context, stochasticity in build events induces flaws in the manufactured piece, compromising the integrity of said part. There is a critical need for in-situ monitoring to spot anomalies once they form, and in this setting we apply our total variation and hyperparameter solvers. We will develop a customized algorithm based on for extreme-scale TV-optimization that achieves super-linear or quadratic-convergence, a critical property for real-time, image-by-image analysis. A worst-case outcome is a preprocessing step that enhances image quality in-situ, specifically for out-of-focus and noisy images.

36 MATERIALS SCIENCE↗

Integrating Multi-Source Data for Bi-Level Traffic Simulator Calibration: A Literature Review and Highway Case Study

Traffic simulation serves as a powerful tool for pre-evaluating policies and technologies. In this context, simulation-based Dynamic traffic assignment (DTA) models are capable of capturing traffic dynamics. They are well-known as critical tools in controlling and predicting traffic situations. The reliability of simulation results heavily depends on the calibration process. Most studies in the literature formulate and calibrate simulators based on a single source of collected data or multiple data sets with the same spatiotemporal characteristics. However, in practice, traffic data is collected by various tools with usually different spatial and temporal resolutions. This study introduces a novel approach to taking into account diverse input data from a variety of sources. An iterative bi-level solution is proposed. to equally treat traffic flow and speed data. The upper level solves flow calibration with the exact solution method, and the lower level calibrates the speed with the simultaneous perturbation stochastic approximation (SPSA) algorithm. Subsequently, the effectiveness of the proposed model is investigated using data from a six-mile section of Nashville's I-24 highway in Tennessee. The results demonstrate that our proposed model creates an effective feedback loop between the optimizer and the simulator for calibrating flow and speed to reduce the error between simulated and real data.

42 ENGINEERING↗

Uncertainty Error Modeling for Non-Linear State Estimation With Unsynchronized SCADA and µPMU Measurements

Distribution systems of the future smart grid require enhancements to the reliability of distribution system state estimation (DSSE) in the face of low measurement redundancy, unsynchronized measurements, and dynamic load profiles. Micro phasor measurement units (µPMUs) facilitate co-synchronized measurements with high granularity, albeit at an often prohibitively expensive installation cost. Supervisory control and data acquisition (SCADA) measurements can supplement µPMU data, although they are received at a slower sampling rate. Further complicating matters is the uncertainty associated with load dynamics and unsynchronized measurements–not only are the SCADA and µPMU measurements not synchronized with each other, but the SCADA measurements themselves are received at different time intervals with respect to one another. This paper proposes a non-linear state estimation framework which models dynamic load uncertainty error by updating the variances of the unsynchronized measurements, leading to a time-varying system of weights in the weighted least squares state estimator. Case studies are performed on the 33-Bus Distribution System in MATPOWER, using Ornstein–Uhlenbeck stochastic processes to simulate dynamic load conditions.

Cooper, Austin↗

Assessing the potential of composite confining systems for secure and long-term CO 2 retention in geosequestration

A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.

58 GEOSCIENCES↗

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗