Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Thermomechanical Relaxation of Vibrationally Excited NO in a Hypersonic Turbulent Boundary Layer

The effects of thermal nonequilibrium are studied using a DC glow discharge plasma to selectively modify the vibrational temperature of an otherwise cold flow. These experiments were performed at 𝑴=5.7and 𝑹𝒆=6×10 6 /m in the Actively Controlled Expansion tunnel at Texas A&M University on a canonical 2.75° half-angle wedge test article. A∼47 W plasma was generated across spanwise electrodes, and removable trips produced either a laminar or turbulent boundary layer. NO planar laser-induced fluorescence (PLIF) quantified the extent of the NO 𝑿 2 𝚷 (𝒗=1) vibrational perturbation due to the plasma and its subsequent interplay with the boundary layer, the first such data of which the authors are aware. NOPLIF measurements were taken at three locations along the test article, and quasi-temperature measurements are reported at each. Despite not quantifiably perturbing the boundary layer, the plasma was shown to alter the NO vibrational temperature profile therein. Independent analysis of the thermodynamics and fluid mechanics of the system suggests each set of physics could explain the results observed near and far from the wall, respectively.

Casey J. Broslawski↗
NASA NTRSDate not supplied

Thermomechanical Relaxation of Vibrationally Excited NO in a Hypersonic Turbulent Boundary Layer

The effects of thermal nonequilibrium are studied using a DC glow discharge plasma to selectively modify the vibrational temperature of an otherwise cold flow. These experiments were performed at 𝑴=5.7and 𝑹𝒆=6×10 6 /m in the Actively Controlled Expansion tunnel at Texas A&M University on a canonical 2.75° half-angle wedge test article. A∼47 W plasma was generated across spanwise electrodes, and removable trips produced either a laminar or turbulent boundary layer. NO planar laser-induced fluorescence (PLIF) quantified the extent of the NO 𝑿 2 𝚷 (𝒗=1) vibrational perturbation due to the plasma and its subsequent interplay with the boundary layer, the first such data of which the authors are aware. NOPLIF measurements were taken at three locations along the test article, and quasi-temperature measurements are reported at each. Despite not quantifiably perturbing the boundary layer, the plasma was shown to alter the NO vibrational temperature profile therein. Independent analysis of the thermodynamics and fluid mechanics of the system suggests each set of physics could explain the results observed near and far from the wall, respectively.

Casey J. Broslawski↗
DOE OSTI2025

Measurement of directed flow in Au + Au collisions at $\sqrt{s{NN}}$ = 19.6 and 27 GeV with the STAR event plane detector

In heavy-ion collision experiments, the global collectivity of final-state particles can be quantified by anisotropic flow coefficients (𝑣 𝑛 ) . The first-order flow coefficient, also referred to as the directed flow (𝑣 1 ), describes the collective sideward motion of produced particles and nuclear fragments in heavy-ion collisions. It carries information on the very early stage of the collision, especially at large pseudorapidity (𝜂), where it is believed to be generated during the nuclear passage time. Directed flow therefore probes the onset of bulk collective dynamics during thermalization, providing valuable experimental guidance to models of the pre-equilibrium stage. In 2018, the Event Plane Detector (EPD) was installed in STAR and used for the Beam Energy Scan phase-II (BES-II) data taking. The combination of EPD (2.1 < |𝜂| < 5.1) and high-statistics BES-II data enables us to extend the 𝑣 1 measurement to the forward and backward 𝜂 regions. In this paper, we present the measurement of 𝑣 1 over a wide 𝜂 range in Au + Au collisions at $\sqrt{s_{NN}}$ = 19.6 and 27 GeV using the STAR EPD. The results of the analysis at $\sqrt{s_{NN}}$ = 19.6 GeV exhibit excellent consistency with the previous PHOBOS measurement, while elevating the precision of the overall measurement. The increased precision of the measurement also revealed finer structures in heavy-ion collisions, including a potential observation of the first-order event-plane decorrelation. Multiple physics models were compared to the experimental results. Only a transport model and a three-fluid hybrid model can reproduce a sizable 𝑣 1 at large 𝜂 as was observed experimentally. The model comparison also indicates 𝑣 1 at large 𝜂 might be sensitive to the QGP phase transition.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Heterogeneous Cluster Energetics and Nonlinear Thermodynamic Response in Supercritical Fluids

Microstructural heterogeneities arising from molecular clusters directly affect the nonlinear thermodynamic properties of supercritical fluids. We present a physical model to elucidate the relation between energy exchange and heterogeneous cluster dynamics during the transition from liquidlike to gaslike conditions. By analyzing molecular-dynamics data and employing physical principles, the model considers contributions from three key processes, namely, changing cluster density, cluster separation, and transfer of molecules between clusters. We show that the proposed model is consistent with the energetics at subcritical conditions and can be used to explain the nonlinear behavior of thermodynamic response functions, including the peak in the isobaric heat capacity.

42 ENGINEERING↗
NASA NTRS2002

Characteristics of Pool Boiling on Graphite-Copper Composite Surfaces

Nucleate pool boiling performance of different liquids on graphite-copper composite (Gr-Cu) surfaces has been experimentally studied and modeled. Both highly wetting fluids, such as freon-113 and pentane, and a moderately wetting fluid (water) were tested on the Gr-Cu surfaces with different graphite-fiber volume fractions to reveal the enhancement effects of the composite surfaces on the nucleate pool boiling. Results of the experiments show that the graphite-fiber volume fraction has an optimum value. The Gr-Cu composite surface with 25 percent graphite-fiber volume (f=0.25) has a maximum enhancement effect on the nucleate boiling heat transfer comparing to the pure copper surface. For the highly wetting fluid, the nucleate boiling heat transfer is generally enhanced on the Gr- Cu composite surfaces by 3 to 6 times shown. In the low heat flux region, the enhancement is over 6 times, but in the high heat flux region, the enhancement is reduced to about 40%. For the moderately wetting fluid (water), stronger enhancement of nucleate boiling heat transfer is achieved on the composite surface. It shown the experimental results in which one observes the nucleate boiling heat transfer enhancement of 5 to 10 times in the low heat flux region and an enhancement of 3 to 5 times in the high heat flux region. Photographs of bubble departure during the initial stage of nucleate boiling indicate that the bubbles detached from the composite surface are much smaller in diameter than those detached from the pure copper surface. Typical photographs are presented.It shows that the bubbles departed from the composite surface have diameters of only O(0.1) mm, while those departed from the pure copper surface have diameters of O(1) mm. It is also found that the bubbles depart from the composite surface at a much higher frequency, thus forming vapor columns. These two phenomena combined with high thermal conductivity of the graphite fiber are considered the mechanisms for such a significant augmentation in nucleate boiling heat transfer on the composite surfaces. A physical model is developed to describe the phenomenon of bubble departure from the composite surface: The preferred site of bubble nucleation is the fiber tip because of higher tip temperature than the surrounding copper base and poor wettability of the graphite tip compared with that of the base material (copper). The high evaporation rate near the contact line produces the vapor cutback due to the vapor recoil pushing the three-phase line outwards from the fiber tip, and so a neck of the bubble is formed near the bubble bottom. Evaporation and surface tension accelerate the necking process and finally result in the bubble departure while a new small bubble is formed at the tip when the surface tension pushes the three-phase line back to the tip. The process is schematically shown. The proposed model is based on and confirmed by experimental results.

Zhang, Nengli↗
DOE OSTI2023

Prediction of hemiwicking dynamics in micropillar arrays

Dynamic hemiwicking behavior is observable in both nature and a wide range of industrial applications ranging from biomedical devices to thermal management. We present a semi-analytical modeling framework (without empirical fitting coefficients) to predict transient capillary-driven hemiwicking behavior of a liquid through a nano/microstructured surface, specifically a micropillar array. In our model framework, the liquid domain is discretized into micropillar unit cells to enable the time marching of the hemiwicking front. A simplified linear pressure drop is assumed along the hemiwicking length such that the local meniscus curvature, contact angle, and effective liquid height are determined at each time step in our transient model. This semi-analytical model is validated with experimental data from our own experiments and from published literature for different fluids. Our model predicts hemiwicking dynamics with <20% error over a broad range of micropillar geometries with height-to-pitch ratio ranging between ≈0.34 and 6.7 and diameter-to-pitch ratio in the range of ≈0.25–0.7 and without any fitting parameters. For lower diameter-to-pitch ratio data points related to sparse micropillar array arrangements, we suggest modifications to the semi-analytical model. This work sheds light on complex and dynamic solid–liquid–vapor interfacial interactions which could serve as a guide for the design of textured surfaces for wicking enhancement in multi-phase thermal and mass transport technologies and applications.

Mechanics↗
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Thermal photon production in Gubser inviscid relativistic fluid dynamics

The Gubser solution to inviscid relativistic fluid dynamics is used to examine the role of transverse expansion on the energy spectrum of photons radiated by quark-gluon plasma. Transverse flow is shown to be a modest effect on the energy spectrum of photons as a whole, despite its large effect on rare high-energy photons produced at low temperatures. An exact expression is derived for the volume of the plasma as a function of its temperature. Furthermore, a simple formula is obtained for the energy spectrum of high-energy thermal photons, which is used to relate the inverse slope T eff of the photon spectrum at energy E to the maximum temperature of the plasma T 0 , finding T eff ≈ T 0 / (1 + $\frac{5}{2}$ T 0 /E).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2023

Thermodynamic behavior of correlated electron-hole fluids in van der Waals heterostructures

Abstract Coupled two-dimensional electron-hole bilayers provide a unique platform to study strongly correlated Bose-Fermi mixtures in condensed matter. Electrons and holes in spatially separated layers can bind to form interlayer excitons, composite Bosons expected to support high-temperature exciton condensates. The interlayer excitons can also interact strongly with excess charge carriers when electron and hole densities are unequal. Here, we use optical spectroscopy to quantitatively probe the local thermodynamic properties of strongly correlated electron-hole fluids in MoSe 2 /hBN/WSe 2 heterostructures. We observe a discontinuity in the electron and hole chemical potentials at matched electron and hole densities, a definitive signature of an excitonic insulator ground state. The excitonic insulator is stable up to a Mott density of ~0.8 × 10 12 cm −2 and has a thermal ionization temperature of ~70 K. The density dependence of the electron, hole, and exciton chemical potentials reveals strong correlation effects across the phase diagram. Compared with a non-interacting uniform charge distribution, the correlation effects lead to significant attractive exciton-exciton and exciton-charge interactions in the electron-hole fluid. Our work highlights the unique quantum behavior that can emerge in strongly correlated electron-hole systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie A., Jr.↗
NASA NTRS1998

Thermal State-of-Charge in Solar Heat Receivers

A theoretical framework is developed to determine the so-called thermal state-of-charge (SOC) in solar heat receivers employing encapsulated phase change materials (PCMS) that undergo cyclic melting and freezing. The present problem is relevant to space solar dynamic power systems that would typically operate in low-Earth-orbit (LEO). The solar heat receiver is integrated into a closed-cycle Brayton engine that produces electric power during sunlight and eclipse periods of the orbit cycle. The concepts of available power and virtual source temperature, both on a finite-time basis, are used as the basis for determining the SOC. Analytic expressions for the available power crossing the aperture plane of the receiver, available power stored in the receiver, and available power delivered to the working fluid are derived, all of which are related to the SOC through measurable parameters. Lower and upper bounds on the SOC are proposed in order to delineate absolute limiting cases for a range of input parameters (orbital, geometric, etc.). SOC characterization is also performed in the subcooled, two-phase, and superheat regimes. Finally, a previously-developed physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) system is used in order to predict the SOC as a function of measurable parameters.

Hall, Carsie, A., III↗
DOE OSTI2022

MPEX High Heat Flux Plasma Dump Design, Manufacture, and Articles Test

Here, the Material Plasma Exposure eXperiment (MPEX) device is a steady-state linear plasma device currently in the final design phase at the Oak Ridge National Laboratory. This device will reach ion fluences up to 10 31 m -2 and will be used to expose neutron-irradiated materials to divertor-relevant plasmas and to study the effects of plasma-material interactions. These studies will elucidate the complex effects of plasmas with divertor candidate materials capable of withstanding high heat flux and high fluences for next-generation fusion devices. Bidirectional plasma will be generated using a high-power (200 kW) helicon source. Plasma will be confined by superconducting magnets. The last plasma-facing component on the upstream side of the MPEX device is the dump, which has been designed to intercept plasma and energetic particles. The dump will have a total heat load of 9.2 kW. A copper alloy (Glidcop AL-15) was selected for use in the water-cooled flange design because of its high thermal conductivity, its retaining strength at elevated temperatures, and its ability to be used in the high-temperature braze joints used in this application. Titanium-zirconium-molybdenum (TZM) tiles are brazed to the Glidcop AL-15 flange using high-temperature braze alloy. External water-cooling channels are used on the dump flange to prevent water leakage inside the vacuum space. This article discusses the details of the high heat flux dump design, including the computational fluid dynamics (CFD) and structural analyses performed to validate the design to meet the operational requirements of the MPEX device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Continued performance improvement and integration of MOOSE's thermal-hydraulics capabilities (M3 Milestone Report)

This work introduces performance, robustness and workflow improvements to Multiphysics Object-Oriented Simulation Environment (MOOSE)-based thermal-hydraulics solvers. It presents work related to the acceleration of segregated fluid dynamics algorithms, which show approximately a factor of 10 speedup compared to the preceding implementation. Additionally, we discuss approaches to use advanced, Schurr complement-based, field split preconditioners for monolithic solution algorithms relying on the finite volume method. The presence of the Rhie-Chow interpolation makes the utilization of this preconditioner challenging, but the results indicate that for a moderately large problem a factor of 3.4 speedup can be achieved in conjunction with a factor of 3.5 reduction in memory usage. Furthermore, we introduce several pseudo-time stepping approaches to MOOSE for the robust convergence to steady-state solutions when steady-state solves don't converge due to the initial guesses being too far from the solution in Newton's method. Every MOOSE-based application has access this algorithm and can benefit from its use. Moreover, several new avenues have been presented for importing meshes from commercial software which make meshing easier. Lastly, the Component system within the Thermal-Hydraulics Module (THM) of MOOSE is abstracted by separating geometry- and physics-related properties.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

A new approach to the study of the turbulence layer

Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.

Levenson, Ruth A.↗