Search NASA⌕ Search

SEARCH · Search NASA

Results for “Thermal 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 181 records · Page 10

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Modeling Nuclear Thermal Propulsion Startup Transients

A poster for the 2021 intern poster session. This poster details the RELAP-7/Griffin model developed for transient simulations to compare the effects of startup sequencing on the propellant efficiency, startup time, and maximum core temperature of a nuclear thermal propulsion (NTP) system. The sequencing parameters under consideration in this study are the ramping rates of both reactivity insertion and hydrogen propellant mass flow rate. Recommendations are made regarding startup sequencing based on the results produced by this model.

33 ADVANCED PROPULSION SYSTEMS↗

CODE-TO-CODE BENCHMARK STUDY FOR THERMAL STRESS MODELING AND PRELIMINARY ANALYSIS OF THE HIGH-TEMPERATURE SINGLE HEAT PIPE EXPERIMENT

Microreactors are very small nuclear reactors with typical thermal-energy output of up to 20 MWth. The microreactor concept is gaining more and more attention for safe, robust, and reliable supply of electricity for remote locations and to meet industrial process heat needs. Heat-pipe-cooled microreactor is one of the microreactors being investigated at Idaho National Laboratory. In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios for the heat-pipe-cooled microreactor designs. This paper discusses finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory, using two commercial software packages, Abaqus and STAR-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the locations where the largest temperature gradients appeared between heat pipe and electric heater. The temperature fields have good agreements between Abaqus and STAR-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A dynamic thermal algal growth model for pilot-scale open-channel raceways

Here this study describes the makeup of and results from a simulation consisting of a lumped thermal model integrated with a dynamic algae growth model to simulate the microalgae productivity of an open-channel raceway cultivation system. The thermal model considers the dynamic effects of weather, light absorption, convective heat transfer, radiation heat transfer, conductive heat transfer, thermal capacitance, and water control strategies. The dynamic algae growth model solves a set of ordinary differential equations consisting of growth functions dependent on incident radiation, temperature, nutrient availability, basal metabolism, and losses due to dark- and photo-respiration. Relative errors in the predicted ash-free dry weight of Nannochloropsis oceanica are 12.5%, -6.1%, and 4.4% for three separate replicates of cultivation cycles (~4.2 m 2 ) performed at AzCATI during the ATP 3 Unified Field Studies in Fall, Spring, and Summer. This research demonstrates that thermal modeling is an essential contributor to the validation of microalgae growth models.

09 BIOMASS FUELS↗

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Cumulants and correlation functions of net-proton, proton, and antiproton multiplicity distributions in Au+Au collisions at energies available at the BNL Relativistic Heavy Ion Collider

Here, we report a systematic measurement of cumulants, C n , for net-proton, proton, and antiproton multiplicity distributions, and correlation functions, κ n , for proton and antiproton multiplicity distributions up to the fourth order in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.5, 19.6, 27, 39, 54.4, 62.4, and 200 GeV. The C n and κ n are presented as a function of collision energy, centrality and kinematic acceptance in rapidity, y, and transverse momentum, p T . The data were taken during the first phase of the Beam Energy Scan (BES) program (2010–2017) at the BNL Relativistic Heavy Ion Collider (RHIC) facility. The measurements are carried out at midrapidity (|y| < 0.5) and transverse momentum 0.4 < p T < 2.0 GeV/c, using the STAR detector at RHIC. We observe a nonmonotonic energy dependence ($\sqrt{s_{NN}}$ = 7.7–62.4 GeV) of the net-proton C 4 /C 2 with the significance of 3.1σ for the 0–5% central Au+Au collisions. This is consistent with the expectations of critical fluctuations in a QCD-inspired model. Thermal and transport model calculations show a monotonic variation with $\sqrt{s_{NN}}$. For the multiparticle correlation functions, we observe significant negative values for a two-particle correlation function, κ 2 , of protons and antiprotons, which are mainly due to the effects of baryon number conservation. Furthermore, it is found that the four-particle correlation function, κ 4 , of protons plays a role in determining the energy dependence of proton C 4 /C 1 below 19.6 GeV, which cannot be understood by the effect of baryon number conservation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advancing Thermal-Hydraulic Modeling Capabilities for Molten Salt Reactors in Pronghorn

This is an M2 report for the NEAMS program. This work focuses on advancing the modeling and simulation capabilities for Molten Salt Reactors (MSRs) within Pronghorn, addressing key thermal-hydraulic challenges that must be addressed for the safe and efficient operation of these advanced reactors. This report summarizes the development and validation of the following new models in Pronghorn: (i) two-equation RANS turbulence models, (ii) near-wall turbulence corrections due to curvature and roughness, and (iii) multidimensional two-phase drift-flux models. For each of these cases, along with the experimental validation, an application study of interest to MSRs is presented. In addition, this report also presents recent improvements in the domain-overlapping coupling between Pronghorn and SAM to address the cases of weakly compressible flow and flow reversal. Finally, an MSR action that greatly simplifies the setup of MSR models using NEAMS tools is described.

42 ENGINEERING↗

Power Electronics Manufacturing Improvements for Heavy-Duty Fuel Cell Vehicles

The Marel Power Solutions project, funded by the U.S. Department of Energy under Award DE-SC0023801, focused on advancing manufacturing techniques for power electronics in heavy-duty fuel cell vehicles. The research aimed to enhance system efficiency, reduce costs, and support broader adoption of hydrogen fuel cell technology. Key areas of investigation included power topology, thermal modeling, system architecture, and accessibility through software tools. Key Accomplishments: 1. Power Topology: - Developed an interleaved boost converter with optimized phase count, leveraging Marel’s proprietary Power Stacks. - Achieved reduced parasitic inductance and resistance, enabling high efficiency in DC-DC converters. 2. Thermal Modeling: - Integrated innovative cooling systems into compact Silicon Carbide (SiC) modules. - Simulations demonstrated the ability to dissipate significant heat (up to 7.5 kW), ensuring device reliability under heavy loads. 3. System Architecture: - Utilized simulation tools to analyze the impact of various fuel cell and vehicle parameters on efficiency. - Highlighted the role of smaller, modular improvements, such as enhanced DC-DC converters, in achieving system-wide gains. 4. Accessibility: - Evaluated and implemented MATLAB/Simulink code generation tools for real-world hardware applications. - Demonstrated the potential for rapid prototyping of custom power systems with reduced development costs. Impact and Benefits: - Efficiency and Cost Reduction: Marel’s cooling technology enhances SiC die performance, reducing the number of dies required and overall system size. - Scalability and Flexibility: The innovations support tailored solutions for diverse applications, from mass transit to mining vehicles. - Sustainability: The research promotes the integration of electrification technologies, helping meet rising energy demands sustainably. Conclusion: The project’s outcomes advance the state of power electronics for hydrogen fuel cell vehicles, enabling more efficient, compact, and cost-effective solutions. These developments lay a foundation for future innovation, contributing to the broader adoption of clean energy technologies in transportation and other industries.

08 HYDROGEN↗

Modeling of thermal and kinetic processes in non-equilibrium plasma ignition applied to a lean combustion engine

In recent years novel ignition systems have been developed to enable stable and efficient engine operations with lean mixtures. Among them, radio-frequency corona ignition systems create discharges that involve a much wider region compared to traditional spark, and produce non equilibrium plasma with high levels of active radicals and excited species. These devices considerably increase the early flame growth speed and extend stable operating limits. With the aim of expanding the knowledge on high efficiency lean-burn SI engines, this paper investigates and compares the combustion development generated by spark and corona ignitions through computational fluid dynamics, within the Reynolds-Averaged Navier-Stokes framework for turbulence modeling. In order to simultaneously take thermal and chemical effects into account, the Perfectly Stirred Reactor combustion model is used. Experimental data are also collected for validation in an optical access engine, for different mixture levels, from stoichiometric to very lean. Furthermore, the faster burn rate generated by the corona system in the initial stage of the combustion is well predicted by the simulations, in all the relative air-fuel ratio conditions. Remarkably, as the mixture becomes lean, simulations are able to capture the non-linear transition from fast to slow kernel growth, before a self-sustainable flame propagation is established. This correlates very well with the measured engine cyclic variability and the corresponding steep change in the duration of the flame kernel formation. Ultimately, this study highlights the important role of the atomic oxygen, as active radical, in promoting and enhancing the combustion initiated by a corona discharge, in addition to the volumetric ignition effect. By contrast, the validated simulations allow to explain that the high-temperature thermal plasma generated in a traditional spark discharge is insensitive to kinetic aspects.

42 ENGINEERING↗

Thermal coupling-decoupling mechanism of heat transfer across van der Waals interfaces in n-eicosane

Progress toward improving thermal transport properties of n-eicosane, a promising phase change material for thermal energy storage near human body temperature, is hindered by lack of a general theoretical framework to model thermal interfacial conductance (TIC) of the van der Waals (vdW) force bonded molecular interfaces, which have little compositional/structural mismatch. Combining molecular dynamics simulations at temperatures up to the melting point and a time-domain analysis of interfacial heat currents in individual molecules, we unveil that (1) the heat flux across a molecular vdW interface has an alternating coupling-decoupling pattern at the single molecule level due to thermal fluctuations, and (2) the net heat transfer during the thermally coupled periods is orders of magnitude larger than that during the thermally decoupled periods. Furthermore, we demonstrate that a vdW interface's TIC correlates almost linearly with two atomistic-scale coupling-decoupling parameters: αc and σc, that represent the duration and strength of the thermal coupling, respectively. This thermal coupling-decoupling mechanism proposes a new paradigm that explicitly accounts for the effects of structural flexibility and dynamical fluctuation on heat transfer across molecular interfaces, and will likely provide new insights into modeling thermal transport properties in a wide-ranging soft materials, including self-assembled monolayers, paraffins or lipids.

42 ENGINEERING↗

Coupled Hydrological-Thermal-Biogeochemical Modeling for Predicting Arctic Carbon Emissions (CH4PACE)

The Arctic, with nearly 1,700 billion metric tons of carbon stored in its soils, is experiencing significant warming trends that raise concerns about the potential release of greenhouse gases, particularly methane. This report investigates the conversion of carbon stored in Arctic soils and permafrost into methane and its subsequent atmospheric release over the coming decades. Recent data indicate a rise in global methane concentrations, highlighting the importance of understanding the mechanisms of methane generation under anaerobic conditions prevalent in thawing permafrost. Our study employs a multi-faceted approach that integrates fieldwork, benchtop experiments, and numerical simulations to explore the subsurface hydrological-thermal-biogeochemical systems influencing methane emissions. We focus on the relationship between microbial activity in decomposing organic matter under oxygen-limited conditions and the transport of methane out of these saturated wetland soils. Key hydrological parameters, including permeability, porosity, and thermal conductivity, were measured on soil samples from Fairbanks, Alaska, and used to parameterize PFLOTRAN, a coupled hydrological-thermal-reactive transport simulator. Our model enables prediction of the evolution of the active layer and its impact on methane production and release. Our findings underscore the critical need for a mechanistic understanding of methane dynamics in permafrost environments to improve predictions of methane emissions in Earth systems models.

54 ENVIRONMENTAL SCIENCES↗

Heat Generation versus Heat Dissipation Maps for Batteries (Battery Safety Maps) [SWR-21-27]

NREL's Battery Safety Map (Heat Generation versus Heat Dissipation Maps for Batteries) generates a response surface that is grouped into regions where the battery is deemed safe from the propagation of thermal runaway, those where an engineering solution would render the battery safe, and finally instances where thermal runaway is inevitable. The demarcation for the different regions is provided by specifications or constraints imposed by the end user. NREL's Battery Safety Map is a user-defined subroutine that sweeps through a set of parameters describing the operating conditions and short-circuit scenarios for battery safety simulations. For each run, the electrochemical-thermal model implementing reaction kinetics and gas generation are run to evaluate the amount of heat generated. A thermal model then captures the amount of heat dissipation based on boundary conditions defined (e.g., forced convection, constant temperature at the surface, etc.). Based on the rate of heat transfer and the amount of heat generated, each run is classified into "safe," "mitigation possible," or "unsafe" zones. These categories are defined based on constraints or tolerance specifications provided by the end user.

Santhanagopalan, Shriram↗

Testing of an Optical Fiber--Based Gamma Thermometer in the High Flux Isotope Reactor Gamma Irradiation Facility

This report describes the design, thermal modeling, and gamma irradiation testing of an optical fiber–based gamma thermometer (OFBGT), which was irradiated in the High Flux Isotope Reactor (HFIR) Gamma Irradiation Facility (GIF). OFBGTs are a promising technology for application in nuclear reactors because they can provide a distributed measurement of gamma ray heating rate, unlike thermocouple-based gamma thermometers, which are fixed in-core sensors that can be used in boiling water reactors to calibrate local power range monitors. OFBGTs measure gamma ray heating rate by measuring the temperature difference between a pair of optical fibers; one fiber is in thermal contact with a heat sink (usually the reactor coolant), and the other is in thermal contact with a thermally isolated mass. The device can be calibrated with a heating wire within the thermal mass. The OFBGT that was designed and fabricated at Oak Ridge National Laboratory can measure distributed gamma ray heating rate over an effective measurement length of 61 cm, and the outer diameter of the sensor is 12.7 mm, giving the prototypical sensor design a relatively small footprint. The sensor housing is backfilled with Ar to ensure a well-predicted thermal response that is not affected by humidity or chemical interactions during operation. For calibration, the sensor design uses a Ni–Cr wire, which can be supplied with currents from 0 to 1 A to capture the wide range of potential gamma ray heating rates expected in HFIR’s spent fuel elements. The OFBGT was thermally modeled analytically and numerically; both models account for temperature-dependent thermal conductivities of the materials and show good agreement. The thermal response of the sensor inside spent HFIR fuel elements was simulated for times up to 1 year after discharge of the fuel element. Out-of-pile open-air tests indicated that the steady-state response of the sensor matches modeled results within experimental uncertainty. Calibration tests were performed using electrical heating in the HFIR spent fuel pool, above the fuel elements, to establish a relationship between the difference in spectral shift measured by optical fibers located inside and outside the OFBGT and the applied electrical heating. Subsequently, the OFBGT was placed within the the spent fuel element from HFIR cycle 501 to measure the spatial profile of the gamma heating rates. Results showed good agreement between the theoretical and measured gamma dose rate profiles, with maximum deviations of ~10% or less.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

From phase decomposition to evaporation: A multi-modal evaluation of thermally degraded model lightweight high-entropy alloy

Lightweight high-entropy alloys (LHEAs) have the potential to replace conventional lightweight materials due to their superior mechanical properties and thermal stability. However, the thermal degradation pattern of LHEAs from phase decomposition to evaporation is not clear. We develop a new Al-based dual phase (FCC + HCP) LHEA—AlTi 0.45 CuZn, and further investigate its thermal degradation behavior for potential high-temperature structural applications. Using multimodal advanced characterization techniques such as differential scanning calorimetry/thermogravimetric analysis, scanning/transmission electron microscopy, and synchrotron X-ray diffraction/pair distribution function (XRD/PDF), a sequence of thermal degradation events beyond the thermal phase stability limit—between 250 and 360 °C—is observed. These include phase decomposition at ~360 °C, Zn evaporation at ~750 °C, and LHEA melting at 880 °C which results in ~25% cumulative weight loss. The formation of Al-Ti phase off the AlTi 0.45 CuZn matrix is due to the largest negative mixing enthalpy for Al-Ti than other binary pairs. Similarly, Zn evaporation from AlTi 0.45 CuZn LHEA is due to its faster evaporation rate than other constituent elements. The high-resolution synchrotron XRD and PDF results support the aforementioned observations; in addition, they reveal local atomic arrangements, local strain, and sluggish grain growth in the LHEA. Among other LHEAs of close density range (5.55 ≤ ρ ≤ 5.85 g/cc), the investigated LHEA exhibits outstanding nano-indentation hardness values due to the coupled grain size effect and HCP phase strengthening of the FCC matrix. As the search for LHEAs for lightweight applications grows, this study shows the potential use of AlTi 0.45 CuZn LHEA for structural applications even at elevated temperatures.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Seasonal Performance Characterization of a Gen3 Particle-Based Concentrating Solar Plant With a Spatially Resolved Transient Thermal Storage Model

Particle-based Gen3 Concentrating Solar Power (CSP) can be paired with high-temperature power cycles (>700 °C) and can have built-in long duration (≥10 hours) thermal energy storage if the working particles are stored properly in thermal energy storage bins. Although high-temperature and long-duration thermal energy storage can mitigate daily intermittencies in solar irradiation, seasonal variability in local meteorological conditions can still have a significant impact on the overall performance of CSP systems. Aside from daily cloud coverage affecting incident solar radiation, wind speed and ambient temperature are also significant variables regarding system heat attenuation, component efficiency, and overall solar conversion efficiency. In this work, we present simulation results for a Gen3 CSP prototype system in operation over four weeks throughout the year in Albuquerque, NM. The meteorological conditions are taken directly from the TMY3 data at the Albuquerque International Sunport, where the hourly Direct Normal Irradiation (DNI), wind speed, and ambient temperature are of particular interest. An investigation of the sensitivity of individual components like storage and ducting to the local meteorological conditions is provided and extended to the overall performance of the CSP system. The results from this study show that the particle inlet temperature at the particle-to-sCO2 heat exchanger can change as fast as 30 °C/min under standard operation in a passive mode.

Plewe, Kaden (ORCID:0000000286826879)↗

Overview of IMPACT Data Acquisition System and Data Reduction Process

This report documents the development of the data acquisition system (DAS) and data reduction methodologies for the Irradiated Material Property Accelerated Characterization Test (IMPACT) experiment at the Advanced Test Reactor (ATR). The IMPACT experiment is designed to enable in-pile measurement of thermal conductivity in metallic nuclear fuels, specifically U-10Zr, using an instrumented thermal conductivity probe. The DAS supports both passive temperature monitoring and active thermal interrogation of the probe through controlled AC and DC excitation. Significant modifications to laboratory-scale systems were required to accommodate the higher resistance paths associated with the in-pile application. Custom electronics and relay-controlled measurement sequencing were developed to enable the measurement and sufficient power delivery to the sensing region. A reduced-order, axisymmetric thermal model based on the thermal quadrupoles method is presented to support data interpretation. This model enables efficient evaluation of transient heat transfer behavior and facilitates solution of the inverse problem required to extract thermal properties from measured signals. Multiple boundary condition formulations are discussed to address varying experimental time scales and geometries. Additionally, machine learning techniques are introduced to support data reduction and improve confidence in inverse solutions. Convolutional neural networks are applied to identify the presence of gas gaps and other evolving geometric features that significantly impact thermal response during irradiation. These efforts contribute to the broader integration of digital twin frameworks and real-time modeling capabilities within the Advanced Fuels Campaign.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

[Presentation Slides] Code-to-Code Benchmark Study for Thermal Stress Modeling and Preliminary Analysis of the High-temperature Single Heat-Pipe Experiment

In the heat-pipe-cooled microreactor, heat pipes remove heat from the reactor core as a passive heat-transfer device, so the fluid circulation is not required for cooling, which can substantially simplify the overall reactor design. However, given the extremely high temperatures in the core region and potentially large temperature gradients across the structure materials, thermal stresses need to be well-analyzed to ensure structural integrity during normal operations and accident scenarios. This presentation slides discuss finite element method-based thermal-stress analysis for the high-temperature single heat-pipe test article in the Single Primary Heat Extraction and Removal Emulator (SPHERE) facility at Idaho National Laboratory (INL), using two commercial software packages, Abaqus and Star-CCM+. A code-to-code benchmark study was performed to crosscheck the model setup and capability of each code and to gain preliminary insights into the potential thermal stress concerns from the current experimental setup. It is observed the significant thermal stresses happen at the inner surface of the heat pipe hole surrounded by electric heaters where the largest temperature gradients appear. The temperature fields have good agreements between Abaqus and Star-CCM+, while the induced thermal stresses show modest deviations probably due to differences of meshing engines used in these two codes. It is found that the local maximum thermal stresses may reach close to the ultimate tensile strength and yield strength of structural material depending on the heater power. Ultimately, the coupled thermal-structural analysis will help guide the current experimental plan and ensure the facility safety for future experimental study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗