Results for “Electron Fluid Model” — Page 16

DOE OSTI2024

Predicting thermophysical properties of molten salts in the MgCl 2 -NaCl-KCl-LiCl system with a shell-model potential

Ternary eutectic salts composed of MgCl 2 , NaCl, and KCl, referred to as MNK salts, have recently emerged as promising candidates as high-temperature heat transfer fluids and thermal energy storage media. Here we performed classical molecular dynamics (MD) simulations to predict the densities, specific heat capacities, viscosities, and ionic self-diffusivities for MNK salts over a wide temperature range. The impact of LiCl additive on their thermophysical properties was also investigated. To capture the electronic polarization of Cl anions by neighboring cations, we developed a novel shell-model potential using the force-matching method and a dataset of ab initio calculated interatomic forces. Our extensive MD simulations predict structure and properties for pure salts and binary/ternary salt mixtures in the MgCl 2 -NaCl-KCl-LiCl system in overall good agreement with available experimental and theoretical data, which corroborates the accuracy and reliability of our developed potential.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Design and Analysis of an Integrated Additively Manufactured Test Article for Plasma-Facing Components

Current plasma-facing components (PFCs) used in helium-cooled divertor modules are complex structures with tungsten tile, steel sleeve components, and cartridges, all assembled in a helium-cooled multiple jet (HEMJ) structure. The goal of this project is to simplify the complex PFC design using additive manufacturing techniques to create a single integrated tungsten test article. Apart from the flexibility this opens up in exploring a wide array of geometries for the article, having a single integrated article significantly reduces the number of joints and parts in the article, thus reducing chances of leaks. A process called electron beam melting has shown to produce very high-density samples and unique geometries, enabling HEMJ or similar designs. To validate and optimize this novel design, the model underwent a series of computational fluid dynamics and finite element analysis simulations to replicate steady-state heat flux in the divertors. The simulations presented in this study consider a steady-state base heat flux of 5 MW/m 2 , with water serving as the coolant. Future research will explore the use of helium as a coolant, simulate edge-localized-mode conditions, and include experimental validation. Since 3D-printed tungsten is anisotropic, the build direction versus build plane of the article are taken into consideration for the test article strength. Because of the high operating temperatures and low ductility of tungsten, thermal creep and brittle fracture are important failure mechanisms to consider. In conclusion, the cap is evaluated with various flow velocities and nozzle diameters, and an optimal design choice is made for which this cap will survive the divertor conditions with a conservative safety margin.

Divertor↗

DOE OSTI2023

CLASSY. VI. The Density, Structure, and Size of Absorption-line Outflows in Starburst Galaxies

Galaxy formation and evolution are regulated by the feedback from galactic winds. Absorption lines provide the most widely available probe of winds. However, since most data only provide information integrated along the line of sight, they do not directly constrain the radial structure of the outflows. In this paper, we present a method to directly measure the gas electron density in outflows (n e ), which in turn yields estimates of outflow cloud properties (e.g., density, volume filling factor, and sizes/masses). We also estimate the distance (r n ) from the starburst at which the observed densities are found. We focus on 22 local star-forming galaxies primarily from the COS Legacy Archive Spectroscopic SurveY (CLASSY). In half of them, we detect absorption lines from fine-structure excited transitions of Si $\tiny{II}$ (i.e., Si $\tiny{II}$ * ). We determine n e from relative column densities of Si $\tiny{II}$ and Si $\tiny{II}$ * , given Si $\tiny{II}$ * originates from collisional excitation by free electrons. We find that the derived n e correlates well with the galaxy's star formation rate per unit area. From photoionization models or assuming the outflow is in pressure equilibrium with the wind fluid, we get r n ~ 1–2r * or ~5r * , respectively, where r * is the starburst radius. Based on comparisons to theoretical models of multiphase outflows, nearly all of the outflows have cloud sizes large enough for the clouds to survive their interaction with the hot wind fluid. Most of these measurements are the first ever for galactic winds detected in absorption lines and, thus, will provide important constraints for future models of galactic winds.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Dielectric Fluids for the Direct Forced Convection Cooling of Power Electronics

The future of electrification of vehicles and other systems will require the creation of high-power density power electronics with low junction-to-fluid thermal resistance cooling solutions. One way to create this solution is to move high-heat-flux liquid cooling (single- or two-phase) as close to the power electronics components as possible. One novel approach involves submersion in dielectric fluids as the cooling solution. We first provide the range of fluid properties and develop a figure of merit (FOM) to aid in dielectric fluid selection. Next, we perform computational fluid dynamics/heat transfer (CFD/HT) modeling using single-phase cooling (submerged jet impingement on an enhanced surface) to validate the dielectric fluid FOM. Results of the study show that the developed FOM is a good representation of the performance of the fluids when compared to the results of the CFD/HT analysis. Both FOM and the CFD/HT analysis show that based on pure thermohydraulic considerations, several commercially available fluids present higher performance, on the order of 5% of water. Finally, the FOM can be used to quickly assess the thermohydraulic performance of a dielectric fluid, as well as the secondary application-specific properties such as boiling point, saturation pressure, flash point, and global warming potential, thereby allowing for fluid candidates to be readily compared.

computational fluid dynamics↗

NASA NTRSDate not supplied

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe↗

NASA NTRSDate not supplied

A Simplified Model of VIPER Thermal Management System. Part I: Loop Heat Pipe

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) uses a combination of loop heat pipes (LHPs) and heaters as the primary hardware to modulate the temperature of the electronic boxes and payloads instruments. LHPs have been designed and instrumented to be capable of passive and active shutdown. The goal of the passive shutdown, using a thermal control valve (TCV) to limit the flow of the working fluid, is to preserve thermal energy and avoid unnecessary heat leaks during lunar night. The goal of the active shutdown on the other hand, is to forcefully increase the vehicle internal temperature to bakeout the payload instruments. Several numerical modeling tools are available to analyze LHP performance under the complex set of operational requirements and environments levied on the VIPER Thermal Management System (TMS). Yet, while those numerical modeling tools are capable of detailed performance analysis, they tend to take large computational resources and long computational time. In this paper, a model of a single LHP developed based on control volume approach is presented. By using the control volume approach, the complex geometry of the LHP, and the components attached to it, are reduced to thermal lump capacitances, thus reducing the computation effort. This paper also describes the effort to correlate the proposed LHP model using thermal vacuum (TVAC) test data.

Loop Heat Pipe↗

DOE OSTI2023

Two-fluid theory of composite bosons and fermions and the quantum Hall proximity effect

Here we propose a two-fluid description of fractional quantum Hall systems, in which one component is a condensate of composite bosons and the other a Fermi liquid formed by composite fermions (or simply electrons). We employ the theory to model the interface between a fractional quantum Hall liquid and a (composite) Fermi liquid metal, where we find a penetration of quantum Hall condensate into the metallic region reminiscent of the proximity effect in superconductor-metal interfaces. We also find a novel and physically reasonable set of gapped quasielectron and neutral modes in fractional quantum Hall liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Design Optimization of A generic Fissile Solution for Mo99 production using Electron Beam-based Neutron Generator using MCNP+CFD

Molybdenum-99 ( 99 Mo) is a medical isotope used in 80% of all medical imaging procedures today. However, 99 Mo is not directly administered to patients for imaging; its decay product, technetium-99m ( 99 mTc), is a metastable isotope with a half-life of ~6 hours. Such a short half-life does not allow for production, separation, and application to a patient before decaying further. Therefore, the solution is to produce the precursor to 99 mTc, 99 Mo. 99 Mo has a half-life of 66 hours, long enough to produce the material and send it to the desired location before it fully decays. Under National Nuclear Security Administration’s (NNSA) highly enriched uranium (HEU) minimization mission, the Material Management and Minimization (M3) office leads the Molybdenum-99 (Mo-99) program. A major goal of the Mo-99 program is to develop methods of producing 99 Mo using low enriched uranium (LEU) opposed to the traditionally used HEU domestically in the U.S. Internationally, LEU is the new standard for use in research reactors and isotope production facility. In this study, we introduce a generic fissile solution system for 99 Mo production facility meeting the design requirement of 1) LEU usage as fissile material and 2) subcritical fissile system coupled with a horizontal electron beam (E-beam) accelerator. The focus of this research is to optimize a fissile solution configuration for maximum 99 Mo production yield using MCNP model, and ultimately implement the optimized design with computational fluid dynamics (CFD) model to understand thermal hydraulic behavior and solution convection characteristic. The report consists of two parts: 1) isotope yield calculation and criticality analysis with a wide range of parametric test matrix, and 2) thermal hydraulic analysis for solution mixing and cooling assessment.

60 APPLIED LIFE SCIENCES↗

NASA NTRSDate not supplied

Fluid & Thermal Analysis of a Manifold Microchannel Heat Sink

A three-dimensional (3-D) numerical model has been developed to study the fluid flow and heat transfer through a manifold-microchannel (MMC) heat sink, a microfluidic cooling system designed as a heat dissipation solution for compact, high-powered electronic systems. An MMC consists of a 3-D manifold structure that distributes working fluid through alternating inlet and outlet manifold channels that guide the fluid to and from the microchannels, shortening the flow length and thereby, reducing the pressure drop across the system. Wide bandgap semiconductors such as silicon carbide (SiC) are increasingly utilized due to their ability to enable smaller electronic devices and components to operate more efficiently at higher frequencies. These semiconductors can greatly benefit from the MMC, which can extract all the heat produced by compact, high-powered electronics, enabling them to achieve optimal performance. The MMC also utilizes an embedded cooling approach which allows for fabrication directly within SiC-based electronics, resulting in significant mass and weight savings and lower thermal resistance than conventional heat sink designs. This paper outlines the methods used to conduct fluid and thermal analysis of different MMC configurations for flow parameters such as flow rate and inlet flow conditions. This study involves both steady state and transient thermal analysis as a gradually increasing heat flux is applied to the system over time, which is representative of the heat dissipation from the electronics over a given period. Results from the analysis show that the MMC is capable of effectively dissipating heat flux of ~100 W/cm2 while maintaining considerably low pressure drop (< 1 kPa) and thermal resistance (< 1 K/W) for certain single-phase working fluids.

manifold microchannel↗

NASA NTRSDate not supplied

Fluid & Thermal Analysis of a Manifold Microchannel Heat Sink

A three-dimensional (3-D) numerical model has been developed to study the fluid flow and heat transfer through a manifold-microchannel (MMC) heat sink, a microfluidic cooling system designed as a heat dissipation solution for compact, high-powered electronic systems. An MMC consists of a 3-D manifold structure that distributes working fluid through alternating inlet and outlet manifold channels that guide the fluid to and from the microchannels, shortening the flow length and thereby, reducing the pressure drop across the system. Wide bandgap semiconductors such as silicon carbide (SiC) are increasingly utilized due to their ability to enable smaller electronic devices and components to operate more efficiently at higher frequencies. These semiconductors can greatly benefit from the MMC, which can extract all the heat produced by compact, high-powered electronics, enabling them to achieve optimal performance. The MMC also utilizes an embedded cooling approach which allows for fabrication directly within SiC-based electronics, resulting in significant mass and weight savings and lower thermal resistance than conventional heat sink designs. This paper outlines the methods used to conduct fluid and thermal analysis of different MMC configurations for flow parameters such as flow rate and inlet flow conditions. This study involves both steady state and transient thermal analysis as a gradually increasing heat flux is applied to the system over time, which is representative of the heat dissipation from the electronics over a given period. Results from the analysis show that the MMC is capable of effectively dissipating heat flux of ~100 W/cm2 while maintaining considerably low pressure drop (< 1 kPa) and thermal resistance (< 1 K/W) for certain single-phase working fluids.

manifold microchannel↗

DOE OSTI2022

Transport of an intense proton beam from a cone-structured target through plastic foam with unique proton source modeling

Laser-accelerated proton beams are applicable to several research areas within high-energy density science, including warm dense matter generation, proton radiography, and inertial confinement fusion, which all involve transport of the beam through matter. Here, we report on experimental measurements of intense proton beam transport through plastic foam blocks. The intense proton beam was accelerated by the 10 ps, 700 $\textit{J}$ OMEGA EP laser irradiating a curved foil target, and focused by an attached hollow cone. The protons then entered the foam block of density 0.38 g/cm 3 and thickness 0.55 or 1.00 mm. At the rear of the foam block, a Cu layer revealed the cross section of the intense beam via proton- and hot electron-induced Cu-K α emission. Images of x-ray emission show a bright spot on the rear Cu film indicative of a forward-directed beam without major breakup. 2D fluid-PIC simulations of the transport were conducted using a unique multi-injection source model incorporating energy-dependent beam divergence. Along with postprocessed calculations of the Cu – $K_α$ emission profile, simulations showed that protons retain their ballistic transport through the foam and are able to heat the foam up to several keV in temperature. The total experimental emission profile for the 1.0 mm foam agrees qualitatively with the simulated profile, suggesting that the protons indeed retain their beamlike qualities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗

DOE OSTI2024

Multiphysics Co-Optimization Design and Analysis of Double-Side Cooled Silicon Carbide-Based Power Module: Preprint

With the rapid growth of Electric Vehicles (EVs) and Hybrid Electric Vehicles (HEVs), much more rigorous design targets have been set for automotive power electronics, including high power density, high reliability, and low cost. Novel power module and inverter technologies based on wide bandgap (WEG) semiconductors have been developed to meet these design targets, while providing optimal power semiconductor operating temperature and promising thermomechanical performance. Compared with conventional cooling techniques which are normally applied only on one side of power module, double-side cooling approach is now believed to be the solution to enable high power density and low thermal resistance of WEG semiconductor-based power electronics. In this work, we develop a three-phase power module that is double-sided cooled using dielectric fluid jet impingement. In each phase, four silicon carbide (SiC) power semiconductors are bonded to copper busbars without electrical insulation layers. A finite element analysis (FEA) model is created for thermal and thermomechanical analysis. Based on FEA modeling results, we select particular dimensions for a parametric study to optimize thermal and mechanical performance. Using a multi-objective genetic algorithm (MOGA)-based optimization method, we have minimized the maximum junction temperature and thermal stresses within the power module. The multiphysics co-optimization approach has enabled an efficient design process of power modules with greatly reduced computational cost, as compared to conventional processes that rely on exhaustive numerical simulations and iterations.

ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2022

3D Exact Analytical Solutions of Two-fluid Plasma, Magnetohydrodynamics, and Neutral Fluid Equations for the Creation of Ordered Structures as well as Jet-like Flows

The 3D exact analytical solutions of ideal two-fluid plasma, single-fluid plasma, and neutral fluid equations have been found using physically justifiable assumptions. Surprisingly these solutions satisfy all nonlinearities in the systems. It is pointed out that these solutions explain the fundamental mechanism behind the creation of a vast variety of ordered structures in plasmas and fluids. In the limiting case of 2D dependence of fields, the theoretical model for plasma is applied to explain the formation of spicules in the solar chromosphere. It is pointed out that the main contribution of electron (ion) baroclinic vectors is to produce vorticity in the plasma, and that magnetic field generation is coupled with the flow of both electrons and ions.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1977

Strong convective and shock wave behavior in solar flares

A model has been developed to study the gasdynamics of a flare region heated by a stream of energetic electrons. It is shown that the energy deposition can introduce strong chromospheric dynamical effects. As a result of fluid motion into rarefied regions, there is considerable redistribution of mass causing a profound influence on the emitted line radiation.

Bloomberg, H. W.↗

NASA NTRS1977

A reexamination of two-fluid solar wind models

The two-fluid solar-wind equations have been solved by a method which is approximately 50 times faster than any previously developed, through the use of asymptotic expansions which are self-consistently iterated upon to find a solution that passes through the critical point. The energy assumptions in two-fluid solar-wind models are reexamined, and the conclusions are as follows: (1) proton thermal conduction may not be neglected, (2) the Coulomb logarithm must be calculated as a function of radius, and (3) the electron and proton temperatures at the base need not be equal, even when the time scale for energy exchange between the species is an order of magnitude smaller than the expansion time at the base. It is possible to reproduce reasonable quiet-time solar-wind parameters at 1 AU, but only if the proton temperature is approximately twice the electron temperature at 1 solar radius. This may indicate that extended proton heating is important in the outer solar corona. Winds with velocities at 1 AU of 450 km/s are generated without nonthermal energy deposition but require high proton temperatures as well as very low densities at the base. Higher-velocity solutions are not possible in a spherically symmetric geometry for reasonable particle fluxes at 1 AU, and it is suggested that these higher-velocity states probably require additional heating, acceleration mechanisms, or nonradial flow.

Nerney, S.↗

NASA NTRS1983

A theory of Jovian decameter radiation

A theory of the Jovian decameter radiation is presented based on the assumed existence of beams of energetic electrons in the inner Jovian magnetosphere. Beam-like electron distributions are shown to be unstable to the growth of both upper hybrid and lower hybrid electrostatic waves. The upconversion of these waves to fast extraordinary mode electromagnetic radiation is calculated by using a fluid model. Two possibilities are considered. First, a random phase approximation is made which leads to a very conservative estimate of intensity that can be expected in decameter radiation. The alternative possibility is also considered, viz, that the upconversion process is coherent. A comparison of both processes suggests that an incoherent interaction may be adequate to account for the observed intensity of decametric radiation, except perhaps near the peak of the spectrum (8 MHz). The coherent process is intrinsically more efficient and can easily produce the observed intensity near 8 MHz if only 0.01% of the energy in the beam is converted to electrostatic energy.

Goldstein, M. L.↗

DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗