Search NASA⌕ Search

SEARCH · Search NASA

Results for “drop”

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

The molten salt tritium transport experiment: A pumped fluoride salt loop for hydrogen isotope experimentation

Molten salt reactors (MSRs) and fusion reactors propose to use molten salts as coolants and breeder blanket materials, respectively. Tritium, however, poses safety concerns in both reactor types due to its ability to permeate through reactor materials and potential for environmental release. This manuscript addresses the tritium transport phenomena in molten salts and presents the design and analysis of the Molten Salt Tritium Transport Experiment (MSTTE). MSTTE is a forced-convection fluoride salt loop intended to measure hydrogen isotope permeation through structural materials in a flowing salt system. In the first phase, MSTTE will use FLiNaK salt and deuterium as surrogates for FLiBe and tritium, with future plans to utilize tritium and FLiBe. MSTTE couples a Copenhagen Atomics pumped salt loop with an external test section that introduces hydrogen isotopes into the loop and measures transport phenomena. The Hydrogen Injection System (HIS) controls hydrogen isotope introduction into the molten salt loop. Here, the permeation test section measures the permeation rate through stainless steel tubing in contact with flowing salt. Computational fluid dynamics (CFD) analysis ensures fully developed salt flow in the permeation test section. MSTTE is modeled with MELCOR-TMAP to predict the permeation rate as a function of experimental variables such as source term, salt flow rate, and salt temperature. Results indicate that the source term is the only parameter with a significant effect on the permeation rate. Pressure drop analysis suggests that the loop should operate below 200 LPM to maintain a pressure drop below 200 kPa. Additionally, finite-element analysis assesses thermal stress during loop operation to ensure the experiment's safe design. MSTTE will provide semi-integral data on tritium transport phenomena in molten salts and serve as a testbed for advancing molten salt technology.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Metal additively manufactured wavy fin cold-plate architecture for improved thermal-hydraulic performance

Rapid growth in artificial intelligence and data center workloads demands high-performance liquid cooling to manage increasing chip power. This study presents two metal-additive-manufactured cold plates with sinusoidal fins, constant-amplitude wavy fins and linearly variable-amplitude wavy fins and compares them against metal-additive-manufactured straight fins using experiments conducted at 1 kW heat dissipation as well as high-fidelity 3D conjugate computational fluid dynamic simulations. The cold plates were printed in AlSi10Mg material and underwent design using a Python-automated workflow prior to manufacture and testing. The experiments show that wavy fins reduce the normalized thermal resistance by 35 to 45 % at water flow rates from 1 to 4 LPM. At a fixed 20 kPa pressure drop, the variable-waviness design lowered peak surface temperature by 9 °C and thermal resistance by 51 %, while edge-channel maldistribution in the constant wavy fin design limited gains. A thermal resistance breakdown revealed that 55–63 % of the total thermal resistance in wavy designs comes from base heat conduction, 27–33 % from fin heat conduction, and 9–13 % from fin heat convection, indicating the need to address conduction bottlenecks. Parametric sweeps identify a 3 mm fin pitch as optimal, and that horizontal inlet/outlet manifolds further reduce pressure drop by 30–60 % and thermal resistance by 9–16 % relative to vertical inlet-outlet manifolds. The results yield comprehensive guidelines for fin geometry, manifold alignment, material selection and additive-manufacturing constraints to realize high-performance liquid-cooled cold plates for power-dense electronics.

3d printing↗
DOE OSTI2026

Optimization of well design and CO 2 injection strategy for risk reduction in Class VI geological carbon sequestration wells

The safety and durability of Class VI wells are critical for geological carbon sequestration (GCS). However, current GCS operations face unique challenges: unlike traditional Class II wells, Class VI CO 2 injection wells operate at rates up to 100 times higher, dramatically increasing the risk of wellbore leakage and structural compromise due to severe temperature drops and associated mechanical stresses. Despite existing guidelines on material selection, there remains a substantial gap in understanding how rapid CO 2 injection rates, low surface temperatures, and variable reservoir conditions interact to threaten long-term well integrity. This study presents a comprehensive, original workflow integrating advanced analytical and numerical models for both well flow and well integrity analysis. By systematically simulating a wide range of field-relevant scenarios—including variations in injection rate, CO 2 temperature, and reservoir pressure—this work provides the first cross-validated assessment of cooling effects on wellbore. The results reveal that extreme temperature drops, up to 60 °C, can occur under high injection rates, particularly in depleted reservoirs, significantly increasing the risk of cement failure. Building on these insights, the study proposes innovative, practical well design and operational strategies, including ductile cement formulations, pre-stressing techniques, advanced insulation coatings, and proactive management of injection rates. The safety of Class VI well extends beyond simply using CO 2 resistant materials. Cement materials should possess optimal thermo-hydraulic-mechanical-chemical properties for effective performance. This work provides a scientific basis for optimizing Class VI well designs, with direct benefits for minimizing environmental risk, lowering operational costs, and enhancing the long-term reliability of GCS.

25 ENERGY STORAGE↗
DOE OSTI2024

2D analysis of tokamak divertor-plasma detachment-bifurcation with operational parameters and geometries

UEDGE simulations with density scans for various input power, transport coefficients and outer poloidal leg length are performed to study the conditions for the existence of a bifurcation-like drop of T e at the outer strike point, commonly referred to as a detachment cliff, when transitioning to a detached plasma from an attached plasma in the outer divertor as the upstream density increases (McLean et al., 2015). The simulation results show that a detachment cliff tends to occur with a higher power input regardless of diffusivities and leg length. Further analysis of change of plasma profiles at a cliff indicate that, in addition to the sharp reduction of the E x B drift fluxes in the outer divertor studied in Jaervinen et al., (2018), the substantial change of the Mach number in the outer divertor and the decrease of the outer mid-plane T e due to the radiation front moving across the separatrix into the confinement region above the X-point consistently occur for all UEDGE density scans that have a detachment cliff. UEDGE time-dependent simulation of the evolution of a detachment cliff shows that the rapid increase of radiation above the X-point occurs in a time scale of ~0.3–0.5, which could possibly be the trigger for the formation of a detachment cliff, quicker than the Mach number change in a time scale of ~1 ms and the drop of T e in a time scale of ~2–3 ms in the outer divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Flow reversal benchmark of a one-sided heated narrow rectangular channel with CATHARE and RELAP5

Flow reversal in narrow coolant channels can be a crucial phenomenon for the safety of research reactors with a downward nominal flow direction. During a loss of forced flow accident, the downward flow stagnates briefly before transitioning into an upward natural circulation flow. The fuel may be damaged if dryout occurs and threshold fuel and/or cladding temperatures are exceeded. A comprehensive study is provided for flow reversal in narrow rectangular channels by examining experimental data and conducting software model analyses. The literature on flow reversal was reviewed, and selected experimental datasets were used to benchmark against CATHARE and RELAP5 models and also compare the code calculations with each other. The experimental data comes from flow reversal tests conducted with a narrow rectangular channel with one-sided heating. The results were compared with experimental data for successful flow reversal tests and predicted dryout power for dryout conditions. Also, the study examined the effects of the pump coastdown period, inlet liquid temperature, system pressure, and localized pressure drops. The experimental results showed that shorter coastdown periods, reduced pressure drops, and lower coolant inlet temperatures increased the dryout power. However, the system pressure did not noticeably affect the results. The simulation results showed that both CATHARE and RELAP5 agreed with experimental data, capturing the trends of the experimental results. Slight differences between each code calculation, as well as the predicted and measured dryout powers, were attributed to experimental uncertainties and the modeling of physical phenomena such as wall nucleation, interfacial heat transfer, drag coefficients, and critical heat flux. Overall, this study provides an understanding of flow reversal and the prediction capabilities of thermal-hydraulics software models. In conclusion, a future study of the flow reversal benchmark of a narrow rectangular channel with two-sided heating may provide additional valuable insights.

CATHARE↗
DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Kelvin-Helmholtz instability as one of the key features for fast and efficient emulsification by hydrodynamic cavitation

The paper investigates the oil–water emulsification process inside a micro-venturi channel. More specifically, the possible influence of Kelvin-Helmholtz instability on the emulsification process. High-speed visualizations were conducted inside a square venturi constriction with throat dimensions of 450 µm by 450 µm, both under visible light and X-Rays. We show that cavity shedding caused by the instability results in the formation of several cavity vortices. Their rotation causes the deformation of the oil stream into a distinct wave-like shape, combined with fragmentation into larger drops due to cavitation bubble collapse. Later on, the cavity collapse further disperses the larger drops into a finer emulsion. Thus, it turns out that the Kelvin-Helmholtz instability is similarly characteristic for hydrodynamic cavitation emulsification inside a microchannel as is the Rayleigh-Taylor instability for acoustically driven emulsion formation.

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

Direct Air Capture Using Aqueous Amino Acid Solvents in a Crossflow Absorber

Carbon dioxide (CO 2 ) is the most abundant of all greenhouse gases (GHGs). CO 2 levels in the atmosphere are 50% higher than in the preindustrial era, trapping heat. CO 2 removal from the atmosphere by direct air capture (DAC) is needed to achieve the internationally agreed global temperature goals. The most common CO 2 capture technology is absorption by aminebased solvents in packed columns. Amino acid solutions have recently gained attention due to their advantages over traditional amine solvents. To be implemented effectively, DAC industrial processes need to handle large airflow rates in separation absorbers. The large air and solvent flow rates preclude the use of countercurrent columns due to high-pressure drops and the occurrence of flooding. Crossflow air−liquid absorbers are used to handle large air and liquid volumes due to their lower-pressure drop. The objective of this work is to study the influence of crossflow absorber geometric parameters and operating conditions on product formation and process efficiency. An already derived theoretical model for countercurrent absorbers has been modified to simulate the operation of a crossflow DAC absorber. The predictive model was implemented into a computer code that was used to study the efficiency of the processes as geometric equipment dimensions and operating parameters vary. Practical suggestions are made to design more efficient DAC processes.

Absorption↗
DOE OSTI2025

Morphological Degradation of Oxygen Evolution Reaction-Electrocatalyzing Nickel Selenides at Industrially Relevant Current Densities

We investigated electrodeposited nanoparticulate nickel selenide (pre)catalysts that transform into nickel oxides/ oxyhydroxides under oxygen evolution reaction conditions in alkaline solutions. Previous studies of this transformation were conducted at lower current densities than those of industrial relevance (≥1 A cm −2 ). We used ultramicroelectrodes (UMEs) to achieve such current densities, benefiting from their small size, ensuring low absolute currents and low ohmic drop but high current densities. Morphological degradation of the catalyst material was only observed at current densities exceeding 1 A cm −2 but not for smaller ones. Using X-ray absorption, Xray photoemission spectroscopy, and X-ray diffraction, we confirmed that the degradation was accompanied by the literature-known transformation of nanoparticulate Ni 3 Se 2 (bulk)/NiSe (surface) into nickel oxyhydroxide. The transformation of the precatalyst goes along with a significant improvement in the charge transfer kinetics observed by decreasing Tafel slopes with ongoing experimental time extracted from cyclic voltammetry (CV) experiments and electrochemical impedance spectroscopy (EIS) in the high-frequency range. However, these kinetic improvements are accompanied by limitations in mass transport concluded from decreasing current responses at high overpotentials in CVs and increasing impedance in the low-frequency range of the EIS spectra after extended CV cycling. These mass transport limitations originated from morphological degradations at the UME exceeding 1 A cm −2 which we proved by applying identical location scanning electron microscopy. This has not been reported in studies that have been limited to lower current densities before. Our findings showcase how UMEs can be used to study (pre)catalysts (herein nickel selenides) under current densities of industrial relevance in the absence of ohmic drop-related ambiguities, combined with in-depth materials characterization studies, e.g., identical location microscopy and advanced spectroscopic methods. This approach enables direct evaluation and comparison of catalyst materials and thus demonstrates how to overcome long-standing limitations of electrocatalyst design and testing.

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

A Versatile Method for Creating Ultrathin Films of Polyzwitterions with Antifouling Properties

Protein repellent (i.e., antifouling) properties of ultrathin polymeric films have been ascribed to the formation of a dense layer of adsorbed polymers, which causes steric hindrance to protein adsorption on underlying substrates. The antifouling properties of the films are hypothesized to be independent of polymer hydrophilicity based on previous studies involving uncharged polymers. Here, to test the hypothesis and extend it to the realm of charged polymers, we have developed a grafting method for creating ultrathin layers of polyzwitterions on solid surfaces. To demonstrate the efficacy of this method, poly[1-(3-sulfopropyl-1)-2-vinylpyridinium betaine] (P2VPPS) was synthesized via free-radical polymerization and immobilized onto surfaces by using a photoreactive benzophenone derivative as the cross-linker, which was covalently bonded to either quartz or Si/SiOx substrates through a silane anchor. We found that ultrathin films with thicknesses of less than 10 nm can be obtained in dry states using 365 nm ultraviolet (UV) cured/cross-linked surface coatings. These films were found to be hydrophilic and undergo structural rearrangements in the presence of water. The antifouling performance of the films was evaluated using Pseudomonas aeruginosa strain PAO1 bacterium as the fouling agent and exhibited antifouling properties with a kinetic control via the coating method (i.e., spin-coating and drop-casting). In particular, the surface area coverages due to the PAO1 cell attachment on the quartz substrates coated with P2VPPS were found to be ∼5% and ∼0% for the spin-coated and drop-casted films, respectively, after 4 h. These results show that the grafting method can be used to generate surfaces with antifouling properties based on ultrathin films of polyzwitterions.

Pseudomonas aeruginosa↗
DOE OSTI2024

Incorrect computation of Madden-Julian oscillation prediction skill

The Madden–Julian oscillation (MJO) is a major tropical weather system and one of the largest sources of predictability for subseasonal-to-seasonal weather forecasts. Skillful prediction of the MJO has been a highly active area of research due to its large socio-economic impacts. Silini et al., herein S21, developed a machine learning model to predict the MJO, which they claimed to have an MJO prediction skill of 26–27 days over all seasons and 45 days for December–February (DJF) winter. If true, this would make the skill of their model competitive with that of the state-of-the-art dynamical MJO prediction systems at 20–35 days. However, here we show that the MJO prediction was calculated incorrectly in S21, which spuriously increased the performance of their model. Correctly computed skill of their model was substantially lower than that reported in S21; the skill for all seasons drops to 11–12 days and the skill for forecasts initialized during DJF drops to 15 days. Our findings clarify that the S21 machine learning model is not competitive with state-of-the-art numerical weather prediction models in predicting the MJO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Insights into the broadband emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Aims.We present the first multiwavelength study of Mrk 501 that contains simultaneous very-high-energy (VHE)γ-ray observations and X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). Methods.We used radio-to-VHE data from a multiwavelength campaign carried out between March 1, 2022, and July 19, 2022 (MJD 59639 to MJD 59779). The observations were performed by MAGIC,Fermi-LAT,NuSTAR,Swift(XRT and UVOT), and several other instruments that cover the optical and radio bands to complement the IXPE pointings. We characterized the dynamics of the broadband emission around the X-ray polarization measurements through its multiband fractional variability and correlations, and compared changes observed in the polarization degree to changes seen in the broadband emission using a multi-zone leptonic scenario. Results.During the IXPE pointings, the VHE state is close to the average behavior, with a 0.2–1 TeV flux of 20%–50% of the emission of the Crab Nebula. Additionally, it shows low variability and a hint of correlation between VHEγ-rays and X-rays. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings, taken in March 2022 (MJD 59646 to 59648 and MJD 59665 to 59667), with a synchrotron peak frequency > 1 keV. For the third IXPE pointing, in July 2022 (MJD 59769 to 59772), the synchrotron peak shifts toward lower energies and the optical/X-ray polarization degrees drop. All three IXPE epochs show an atypically low Compton dominance in theγ-rays. The X-ray polarization is systematically higher than at lower energies, suggesting an energy stratification of the jet. While during the IXPE epochs the polarization angles in the X-ray, optical, and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. Such results further support the hypothesis of an energy-stratified jet. We modeled broadband spectra taken simultaneous to the IXPE pointings, assuming a compact zone that dominates in the X-rays and the VHE band, and an extended zone stretching farther downstream in the jet that dominates the emission at lower energies.NuSTARdata allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change in the magnetization and/or the emission region size, which directly connects the shift in the synchrotron peak to lower energies with the drop in the polarization degree.

Astronomy & Astrophysics↗
DOE OSTI2025

The Three Hundred Project: Modeling baryon and hot-gas fraction evolution in simulated clusters

The baryon fraction of galaxy clusters, expressed as the ratio between the mass in baryons (including both stars and cold or hot gas) and the total mass, is a powerful tool to provide information on the cosmological parameters, while the hot-gas fraction provides indications on the physics of the intracluster plasma and its interplay with the processes that drive galaxy formation. Using cosmological hydrodynamical simulations of about 300 simulated massive galaxy clusters with a median mass M 500 ≈ 7 × 10 14 M ⊙ at z = 0, we model the relations between total mass and either baryon fraction or the hot gas fractions at overdensities Δ = 2500, 500, and 200 with respect to the cosmic critical density, and their evolution from z ∼ 0 to z ∼ 1.3. We utilized the simulated galaxy clusters from the Three Hundred project, which include star formation and feedback from both supernovae and active galactic nuclei. We fit the simulation results for such scaling relations against three analytic forms (linear, quadratic, and logarithmic in a logarithmic plane) and three forms for the redshift dependence, and we considered as a variable both the inverse of the cosmic scale factor, (1 + z), and the Hubble expansion rate, E(z). We show that power-law dependencies on cluster mass poorly describe the investigated relations. A power law fails to simultaneously capture the flattening of the total baryon and gas fractions at high masses, their drop at low masses, and the transition between these two regimes. The other two functional forms provide a more accurate description of the curvature in mass scaling. The fractions measured within smaller radii exhibit a stronger evolution than those measured within larger radii. From the analysis of these simulations, we evince that as long as we include systems in the mass range herein investigated, the baryon or gas fraction can be accurately related to the total mass through either a parabola or a logarithm in the logarithmic plane. The trends are common to all modern hydro simulations, although the amplitude of the drop at low masses might differ. Being able to observationally determine the gas fraction in groups will thus provide constraints on the baryonic physics.

galaxy clusters↗
DOE OSTI2025

Neutron generation dynamics inside a MA-class dense plasma focus Z-pinch

Dense plasma focii (DPFs) are appealing as energy efficient sources of short pulses of ions, neutrons, and x rays. The output of these sources is expected to scale with input current (I4), but has been shown to drop at the MA level [S. K. H. Auluck, “On the failure of neutron yield scaling in the dense plasma focus,” Phys. Plasmas 30, 080701 (2023)]. New results on the MegaJOuLe Neutron Imaging Radiography DPF showed neutron yield production in agreement with the input current scaling beyond the previously observed drop. This work provides insight into the pinch formation on a DPF and reports on the two different mechanisms leading to neutron generation inside a DPF using a combination of kinetic simulations and experimental data. A combination of particle-in-cell (PIC) and 1D shock theory results are used to describe the pinch formation and disassembly and the corresponding thermonuclear and beam-target mechanisms. The temporal evolution of the pinch column predicted by the PIC simulations shows qualitative agreement with the experimental data from plasma photon emission as well as temporal neutron pulse shapes. In MJ-class DPFs, both thermonuclear and beam-target mechanisms can occur over the course of the implosion and contribute to the total neutron production. Hence the neutron source size of a DPF will change throughout the implosion. Experimental neutron radiographs show the increase in source size as the pinch breaks apart, in agreement with simulation's prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Superlattice reflection signatures of the insulator–metal transition in V 3 O 5 thin films

Nanoelectronic systems that are inspired by the brain are increasingly looking to insulator–metal transition (IMT) materials as they can mimic the response characteristics of neurons to temperature changes so that these can be used in robotic and computational applications. V 3 O 5 has an insulator–metal transition at ∼430 or ∼80 K higher than VO 2 and provides a unique high-temperature opportunity for these types of applications. Here, in this work, we track the structural evolution of V 3 O 5 thin films across the IMT through conventional selected-area electron diffraction (SAED) and four-dimensional scanning TEM (4D-STEM), correlated with temperature-dependent resistance measurements. SAED patterns show reversible evidence of superlattice reflections associated with the IMT—present below T IMT and absent above it—consistent with the accompanying drop in resistance. At room temperature, nanobeam electron diffraction patterns further reveal three local configurations: (i) type I regions with clean patterns lacking superlattice reflections and spot splitting; (ii) type II regions exhibiting rows of superlattice reflections and split spots indicative of crystallographic variants; and (iii) type III regions with negligible superlattice reflections but larger spot splitting suggestive of overlapping domains of insulating and conducting phases likely driven by local lattice distortions. Upon heating, the superlattice reflections disappear between 413 and 453 K, concurrent with the resistance drop at T IMT , consistent with the emergence of a conducting phase. The overall diffraction geometry remains essentially unchanged up to 573 K, implying that relative domain orientations persist through the transition. These observations reveal nanoscale structural heterogeneity in V 3 O 5 thin films across the IMT and inform operation in regimes where mixed-phase textures are expected. A plausible indexing framework rationalizing the observed geometries is presented in the Discussion section, alongside its limitations and alternative interpretations.

25 ENERGY STORAGE↗
DOE OSTI2024

A RELAP5-3D Model of the Lobo Lead Loop

The goal of our research is to build upon the capability of RELAP5-3D to model molten lead systems. Molten lead has several potential uses in future advanced reactors like the lead fast reactor or fusion reactors that utilize dual-coolant lead lithium blankets. This potential for use in future generations of reactors highlights the necessity to develop molten lead models to ensure that they can accurately predict the thermohydraulic behavior. We have developed a RELAP5-3D model of the Lobo Lead Loop facility, located at the University of New Mexico, to verify the accuracy of RELAP5-3D via comparison to existing computational fluid dynamics results and analytical calculations. It was found that RELAP5-3D accurately calculated radiative heat transfer (within < 1%) when compared to theoretical calculations. In addition, pressure drop calculations done in RELAP5-3D demonstrated reasonable agreement within 20 kPa, mostly within ~7-15%, when compared to the computational fluid dynamics model of the facility developed by the University of New Mexico, and captured the dependence of pressure drop on flow velocity accurately. Finally, a hypothetical loss of flow transient was imposed on the RELAP5-3D model to determine the feasibility of performing a similar experiment with the Lobo Lead Loop. It was found that such an experiment could be possible as the RELAP5-3D model indicated that the temperatures of the fluid would not exceed limiting temperatures of the structure (1658 K) nor the maximum temperature of the electromagnetic pump inlet (823 K). Although there is not experimental data to begin validation, the model will be readily available for future validation studies when the experimental data is generated, especially as the model continues to evolve over time. Furthermore, the results so far demonstrate a promising first step in the verification/validation of the RELAP5-3D model of the Lobo Lead Loop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Large-scale Multiphysics Simulations of Small Modular Reactors Operating in Natural Circulation

Thanks to the advancements in high-performance computing, advanced modeling and simulation have become crucial in driving the development and deployment of next-generation nuclear reactors, such as small modular reactors (SMRs). SMRs offer the promise of cost-effective baseload electricity production and improved safety, while addressing some of the challenges associated with large reactor designs, such as high capital costs and extended construction timelines. As part of the Exascale Computing Project, the large-scale multiphysics simulation of an entire SMR primary system has been achieved by combining computational fluid dynamics and neutronics. In addition to the successful demonstration of full-core SMR simulations, the current study integrated the impact of natural circulation into the system. Natural circulation is the primary mechanism driving coolant circulation in SMRs. The mass flow rate in the core depends on the core power, and a numerical model has been developed to predict it. The pressure drop caused by the helical coil steam generator was also accounted for by developing a pressure drop correlation based on high-fidelity large eddy simulation results, further improving prediction accuracy. In conclusion, the results of the study demonstrate that the implemented natural circulation model is effective in predicting the responses of SMR full-core multiphysics simulations.

ECP↗
DOE OSTI2024

Velocity-space distribution function of fast ions in a sawtoothing plasma

This study explores the influence of sawtooth oscillations on the velocity space distribution of fast ions in tokamak plasma discharges. The relevant Fokker–Planck equation for fast ions is solved analytically. Two distinct effects arising from the temperature drop associated with a sawtooth crash and their impact on the distribution function of fast ions are considered. The first effect involves the modulation of the fusion alpha particle source on the timescale of the sawtooth period, linked to the drop in fusion yield resulting from the sawtooth temperature relaxations. The second effect is tied to the increase of the slowing-down time during the sawtooth ramp, causing particles born later in the sawtooth cycle to experience reduced slowing down compared to those born right after the crash, creating an accumulation-like mechanism at higher energies. In regimes where the sawtooth period is shorter than the fast ion slowing-down time, the combined influence of these effects gives rise to fast ion distribution functions that transiently exhibit positive slopes in velocity space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗