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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Performance Comparison of Helium and sCO 2 as Coolants for Modular Finger-Type Divertors

Several studies at the Georgia Institute of Technology have evaluated the thermal and fluids performance of the helium-cooled modular divertor with multiple jets (HEMJ) over the past decade. This finger-type divertor was studied both experimentally at nearly prototypical conditions and numerically at fully prototypical operating conditions using experimentally validated simulations. Recently, supercritical carbon dioxide (sCO 2 ) has been studied as the primary coolant in power cycles and other applications in various systems, in part because CO 2 achieves the high densities typical of supercritical fluids at relatively low temperatures and pressures, with a critical point of (7.38 MPa, 31°C). This density makes it possible to realize very compact and efficient sCO 2 power cycles. The feasibility of sCO 2 as a coolant for plasma-facing components, specifically the divertor, was therefore evaluated as part of the Fusion Energy System Studies design study activities. This work compares the thermal-fluid performance of helium and sCO 2 in the HEMJ divertor geometry using numerical simulations at prototypical conditions: inlet temperatures T i = 600°C to 700°C, pressures p ≈ 10 MPa, and steady-state incident heat fluxes on the tile q″ < 17 MW/m 2 . The performance is quantified here as the maximum heat flux that can be accommodated by the plasma-facing tile, the pumping power fraction, defined as the ratio of the coolant pumping power to the incident thermal power, and the operating stress limits based on ASME pressure vessel criteria. As expected, helium requires lower mass flow rates and pumping power fractions within imposed maximum temperature limits for the HEMJ pressure boundary. However, it also appears that neither helium nor sCO 2 can remove 10 MW/m 2 of incident heat flux while meeting ASME pressure vessel criteria. Finally, the numerical modeling reveals that sCO 2 may remove slightly higher incident heat fluxes than helium due to the imposed stress limits due to the sCO 2 coolant resulting in smaller local temperature gradients, albeit at a considerably higher pumping power fraction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Validation of Loci-Stream for Autogenous Pressurization of Cryogenic Propellant Tank

Autogenous pressurization of cryogenic propellant tanks eliminates the need to have an additional pressurant tank on the space vehicle, which is highly advantageous due to reduced vehicle mass and design complexity. Autogenous pressurization therefore is one of the key technologies for deep space exploration and long-term space missions. The complex interaction of thermal gradients, turbulence and phase change near the interface make the problem a challenging one to model. Nodal analysis tools and reduced order models are unable to capture the necessary physics. 3-D CFD analyses are necessary to fully characterize autogenous pressurization. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. NASA MSFC's Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. We demonstrate our modeling approach and validation of Loci-Stream for predicting autogenous pressurization of a flight scale propellant tank in order for the solver to serve as a reliable design and analysis tool for NASA's CFM application needs. Liquid hydrogen tank pressurization tests carried out at the MSFC test stand 300 facilities provide reliable validation data for this purpose. These tests were modeled using the Loci-Stream solver with a newly implemented two-phase sharp interface treatment. We show that our modeling approach and CFD solver predict the autogenous pressurization phenomena satisfactorily, and document challenging aspects of modeling this problem.

cryogenic fluid management↗
NASA NTRSDate not supplied

Validation of Loci-Stream for Autogenous Pressurization of Cryogenic Propellant Tank

Autogenous pressurization of cryogenic propellant tanks eliminates the need to have an additional pressurant tank on the space vehicle, which is highly advantageous due to reduced vehicle mass and design complexity. Autogenous pressurization therefore is one of the key technologies for deep space exploration and long-term space missions. The complex interaction of thermal gradients, turbulence and phase change near the interface make the problem a challenging one to model. Nodal analysis tools and reduced order models are unable to capture the necessary physics. 3-D CFD analyses are necessary to fully characterize autogenous pressurization. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. NASA MSFC's Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. We demonstrate our modeling approach and validation of Loci-Stream for predicting autogenous pressurization of a flight scale propellant tank in order for the solver to serve as a reliable design and analysis tool for NASA's CFM application needs. Liquid hydrogen tank pressurization tests carried out at the MSFC test stand 300 facilities provide reliable validation data for this purpose. These tests were modeled using the Loci-Stream solver with a newly implemented two-phase sharp interface treatment. We show that our modeling approach and CFD solver predict the autogenous pressurization phenomena satisfactorily, and document challenging aspects of modeling this problem.

CFM↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2010

Orion Entry Performance-Based Center-of-Gravity Box

The Orion capsule is designed both for Low Earth Orbit missions to the ISS and for missions to the moon. For ISS class missions, the capsule will use an Apollo-style direct entry. For lunar return missions, depending on the timing of the mission, the capsule could perform a direct entry or a skip entry of up to 4800 n.mi. in order to land in the coastal waters of California. The physics of atmospheric re-entry determine the capability of the Orion vehicle. For a given vehicle mass and shape, physics tells us that the driving parameters for an entry vehicle are the hypersonic lift-to-drag ratio (L/D) and the flight path angle at entry interface (gamma(sub EI)). The design of the Orion atmospheric re-entry must meet constraints during both nominal and dispersed flight conditions on landing accuracy, heating rate, total heat load, sensed acceleration, and proper disposal of the Service Module. These constraints define an entry corridor in the space of L/D-gamma(sub EI); if the vehicle falls within this corridor, then all constraints are met. The gamma(sub EI) dimension of the corridor can be further constrained by the gloads experienced during emergency entries. Thus, the entry performance for the Orion vehicle can be described completely by the L/D. Bounds on the hypersonic L/D necessary to achieve all the mission requirements can be defined for the given entry corridor. Landing accuracy performance drives the lower limit on L/D. In order to achieve the desired landing accuracy, a minimum L/D must be ensured. The design of the Thermal Protection System (TPS) drives the upper limit on L/D. A higher L/D can drive mass into the design of the TPS. Conversely, once the TPS is designed, the L/D must be ensured to stay below a certain limit in order for the TPS to stay within its design envelop. The L/D must stay within its upper and lower bounds during dispersed flight conditions. L/D is a function of both the aerodynamics and the center-of-gravity (CG) of the vehicle. The aerodynamics of the vehicle are determined by Computational Fluid Mechanics (CFD) and wind tunnel tests. However, the aerodynamics are not known precisely. Instead, an aerodynamic database has been developed where the aerodynamic coefficients are known to fall within a probabilistic band defined by upper and lower bounds. It is expected that the probabilistic band will shrink after the first missions are flown and real-world data is collected. Until that time, the Orion must be designed to the current aerodynamic database. Thus, for a given aerodynamic database with given uncertainties, the allowable range in L/D can be mapped to an allowable box for the CG location. The CG box is used to set requirements on the dispersions allowed for vehicle packaging and cargo storage. As the aerodynamic uncertainties decrease, the size of the CG box can increase. This paper discusses the technique used to map the minimum and maximum L/D bounds set by the entry performance requirements to the allowable dispersions in CG while accounting for aerodynamic uncertainties. The L/D is defined as the ratio of the lift force to the drag force. It is equivalent to the ratio of lift coefficient (C(sub L)) over drag coefficient (C(sub D)). C(sub L) and C(sub D) are functions of Mach number (M) and angle of attack (alpha). A Mach number of 25 is used as a measuring point of the hypersonic L/D. Variations in C(sub L), C(sub D) and alpha cause variations in L/D. Equation (1) shows the three contributions to the variation in L/D.

Rea, Jeremy R.↗
NASA NTRS2011

Z-Pinch Magneto-Inertial Fusion Propulsion Engine Design Concept

Fusion-based nuclear propulsion has the potential to enable fast interplanetary transportation. Due to the great distances between the planets of our solar system and the harmful radiation environment of interplanetary space, high specific impulse (Isp) propulsion in vehicles with high payload mass fractions must be developed to provide practical and safe vehicles for human spaceflight missions. Magneto-Inertial Fusion (MIF) is an approach which has been shown to potentially lead to a low cost, small fusion reactor/engine assembly (1). The Z-Pinch dense plasma focus method is an MIF concept in which a column of gas is compressed to thermonuclear conditions by an estimated axial current of approximately 100 MA. Recent advancements in experiments and the theoretical understanding of this concept suggest favorable scaling of fusion power output yield as I(sup 4) (2). The magnetic field resulting from the large current compresses the plasma to fusion conditions, and this is repeated over short timescales (10(exp -6) sec). This plasma formation is widely used in the field of Nuclear Weapons Effects (NWE) testing in the defense industry, as well as in fusion energy research. There is a wealth of literature characterizing Z-Pinch physics and existing models (3-5). In order to be useful in engineering analysis, a simplified Z-Pinch fusion thermodynamic model was developed to determine the quantity of plasma, plasma temperature, rate of expansion, energy production, etc. to calculate the parameters that characterize a propulsion system. The amount of nuclear fuel per pulse, mixture ratio of the D-T and nozzle liner propellant, and assumptions about the efficiency of the engine, enabled the sizing of the propulsion system and resulted in an estimate of the thrust and Isp of a Z-Pinch fusion propulsion system for the concept vehicle. MIF requires a magnetic nozzle to contain and direct the nuclear pulses, as well as a robust structure and radiation shielding. The structure, configuration, and materials of the nozzle must meet many severe requirements. The configuration would focus, in a conical manner, the Deuterium-Tritium (D-T) fuel and Lithium-6/7 liner fluid to meet at a specific point that acts as a cathode so the Li-6 can serve as a current return path to complete the circuit. In addition to serving as a current return path, the Li liner also serves as a radiation shield. The advantage to this configuration is the reaction between neutrons and Li-6 results in the production of additional Tritium, thus adding further fuel to the fusion reaction and boosting the energy output. To understand the applicability of Z-Pinch propulsion to interplanetary travel, it is necessary to design a concept vehicle that uses it. The propulsion system significantly impacts the design of the electrical, thermal control, avionics, radiation shielding, and structural subsystems of a vehicle. The design reference mission is the transport of crew and cargo to Mars and back, with the intention that the vehicle be reused for other missions. Several aspects of this vehicle are based on a previous crewed fusion vehicle study called Human Outer Planet Exploration (HOPE), which employed a Magnetized Target Fusion (MTF) propulsion concept. Analysis of this propulsion system concludes that a 40-fold increase of Isp over chemical propulsion is predicted. This along with a greater than 30% predicted payload mass fraction certainly warrants further development of enabling technologies. The vehicle is designed for multiple interplanetary missions and conceivably may be suited for an automated one-way interstellar voyage.

Miernik, Janie H.↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inferring thermal ion temperature and residual kinetic energy from nuclear measurements in inertial confinement fusion implosions

In inertial con finement fusion implosion experiments, the presence of residual anisotropic fluid motion within the stagnating hot spot leads to signifi cant variations in ion-temperature measurements using neutron time- of-flight detectors along different lines of sight. The minimum ion-temperature measurement is typically used as representative of the thermal temperature. In the presence of isotropic flows, however, even the minimum DT neutron-inferred ion temperature can be well above the plasma thermal temperature. Using both DD and DT neutron-inferred ion-temperature measurements, we show that it is possible to determine the contribution of isotropic flows and infer the DT burn-averaged thermal ion temperature. The contribution of large isotropic flows on driving the ratio of DD to DT neutron-inferred ion temperatures well below unity and approaching the lower bound of 0.8 is demonstrated in multimode simulations. The minimum DD neutron- inferred ion temperature is determined from the velocity variance analysis, accounting for the presence of isotropic flows. Being close to the DT burn-averaged thermal ion temperature, the inferred DD minimum ion temperatures demonstrate a strong correlation with the experimental yields in the OMEGA implosion database. An analytical expression is also derived to explain the effect of mode ℓ = 1 ion-temperature measurement asymmetry on yield degradations caused by the anisotropic flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗
NASA NTRS2000

The Life Cycle Application of Intelligent Software Modeling for the First Materials Science Research Rack

Marshall Space Flight Center (MSFC) has been funding development of intelligent software models to benefit payload ground operations for nearly a decade. Experience gained from simulator development and real-time monitoring and control is being applied to engineering design, testing, and operation of the First Material Science Research Rack (MSRR-1). MSRR-1 is the first rack in a suite of three racks comprising the Materials Science Research Facility (MSRF) which will operate on the International Space Station (ISS). The MSRF will accommodate advanced microgravity investigations in areas such as the fields of solidification of metals and alloys, thermo-physical properties of polymers, crystal growth studies of semiconductor materials, and research in ceramics and glasses. The MSRR-1 is a joint venture between NASA and the European Space Agency (ESA) to study the behavior of different materials during high temperature processing in a low gravity environment. The planned MSRR-1 mission duration is five (5) years on-orbit and the total design life is ten (IO) years. The MSRR-1 launch is scheduled on the third Utilization Flight (UF-3) to ISS, currently in February of 2003). The objective of MSRR-1 is to provide an early capability on the ISS to conduct material science, materials technology, and space product research investigations in microgravity. It will provide a modular, multi-user facility for microgravity research in materials crystal growth and solidification. An intelligent software model of MSRR-1 is under development and will serve multiple purposes to support the engineering analysis, testing, training, and operational phases of the MSRR-1 life cycle development. The G2 real-time expert system software environment developed by Gensym Corporation was selected as the intelligent system shell for this development work based on past experience gained and the effectiveness of the programming environment. Our approach of multi- uses of the simulation model and its intuitive graphics capabilities is providing a concurrent engineering environment for rapid prototyping and development. Operational schematics of the MSRR-1 electrical, thermal control, vacuum access, and gas supply systems, and furnace inserts are represented graphically in the environment. Logic to represent first order engineering calculations is coded into the knowledge base to simulate the operational behavior of the MSRR-1 systems. An example of engineering data provided includes electrical currents, voltages, operational power, temperatures, thermal fluid flow rates. pressures, and component status indications. These type of data are calculated and displayed at appropriate instrumentation points, and the schematics are animated to reflect the simulated operational status of the MSRR-1. The software control functions are also simulated to represent appropriate operational behavior based on automated control and response to commands received by the crew or ground controllers. The first benefit of this simulation environment is being realized in the high fidelity engineering analysis results from the electrical power system G2 model. Secondly, the MSRR-1 simulation model will be embedded with a hardware mock-up of the MSRR-1 to provide crew training on MSRR-1 integrated payload operations. G2 gateway code will output the simulated instrumentation values, termed as telemetry, in a flight-like data stream so that the crew has realistic and accurate simulated MSRR-1 data on the flight displays which will be designed for crew use. The simulation will also respond appropriately to crew or ground initiated commands, which will be part of normal facility operations. A third use of the G2 model is being planned; the MSRR-1 simulation will be integrated with additional software code as part of the test configuration of the primary onboard computer, or Master Controller, for MSRR-1. We will take advantage of the G2 capability to simulate the flight like data stream to test flight software responses and behavior. A fourth use of the G2 model will be to train the Ground Support Personnel that will monitor the MSRR-1 systems and payloads while they are operating aboard the ISS. The intuitive, schematic based environment will provide an excellent foundation for personnel to understand the integrated configuration and operation of the MSRR-1, and the anticipated telemetry feedback based on operational modes of the equipment. Expert monitoring features will be enhanced to provide a smart monitoring environment for the operators. These features include: (1) Animated, intuitive schematic-based displays which reflect telemetry values, (1) Real-time plotting of simulated or incoming sensor values, (3) High/Low exception monitoring for analog data, (4) Expected state monitoring for discrete data, (5) Data trending, (6) Automated malfunction procedure execution to diagnose problems, (7) Look ahead capability to planned MSRR-1 activities in the onboard timeline. And finally, the logic to calculate telemetry values will be deactivated, and the same environment will interface to the incoming data for the real-time telemetry stream to schematically represent the onboard hardware configuration. G2 will be the foundation for the real-time monitoring and control environment. In summary, our MSRR-1 simulation model spans many elements of the life cycle development of this project: Engineering Analysis, Test and Checkout, Training of Crew and Ground Personnel, and Real-time monitoring and control. By utilizing the unique features afforded by an expert system development environment, we have been able to synergize a powerful tool capable of addressing our project needs at every phase of project development.

Rice, Amanda↗
DOE OSTI2025

An intercomparison of wall fluxes in a turbulent thermal convection chamber: Direct numerical simulations and wall-modeled large-eddy simulations enhanced by machine learning

Thermal convection in a closed chamber is driven by a warm bottom, a cold top, and side walls at various temperatures. Although wall fluxes are the source of convection energy, accurately modeling these fluxes (i.e., the wall model) is challenging. In large-eddy simulations (LESs), many wall models are traditionally derived from the canonical boundary layer, which may be unsuitable for thermal convection bounded by both horizontal and vertical walls. This study conducts a model intercomparison of dry convection in a cubic-meter chamber using three direct numerical simulations (DNSs) and four LESs with different wall models. The LESs employ traditional wall models, a new wall model employing physics-aware neural networks, and a refined grid near the walls. The experiment involves four cases with varying sidewall temperatures. Our results show that LESs capture the main flow features and the trends of mean fluxes. The physics-aware neural networks and refined wall grids can improve the temporally averaged local fluxes when the large-scale circulation has a preferred direction. Even without the local improvement of wall fluxes, the LES flow quantities (temperature and velocities) can still largely match those in DNSs, provided the mean flux largely matches the DNSs. Additionally, DNSs reveal that a variation in corner treatments has minimal impacts on the flow quantities away from corners. Finally, LESs underestimate the mean fluxes of the entire wall due to their inability to resolve corner regions, but their mean flux away from the corner can better match DNS.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Hamiltonian kinetic-Hall magnetohydrodynamics with fluid and kinetic ions in the current and pressure coupling schemes

We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Qualification and Commissioning of Helium Flow Loop Experiment for Blanket Design Measurements

Sufficient cooling of plasma-facing materials remains an outstanding challenge in the design of fusion reactor blankets in commercial power demonstration plants. Due to its chemical inertness and low neutron interaction cross section, pressurized helium is a candidate coolant fluid for such systems; however, helium has a small thermal mass compared to liquid coolants, potentially reducing heat removal performance. To address this need, a number of heat transfer enhancements have been proposed to improve the cooling efficiency of such components, thereby decreasing pumping power needs and improving overall plant efficiency. Toward this end, a helium flow loop experiment (HFLE) has been designed and commissioned to test advanced passive heat transfer enhancements in unit-cell test sections, providing necessary data for model validation and subsequent system design. The HFLE is designed to provide flow of pressurized (up to 4 MPa) helium at flow rates up to 80 g/s, enabling heat transfer and pressure drop measurements in test pieces at Reynolds numbers in excess of 180 000. To explore the effects of novel and complex heat transfer enhancements, test sections are produced via additive manufacturing, providing geometries not typically obtainable by conventional machining. Here in this work, we present results from HFLE commissioning and the initial thermal-hydraulic tests of an additively manufactured rifled-rib test section. Results are compared to smooth pipe correlations, and plans are described for future HFLE measurements. These preliminary experiments indicate the utility of the HFLE for heat transfer enhancement testing and simulation validation activities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

Performance of Simple Gas Foil Thrust Bearings in Air

Foil bearings are self-acting hydrodynamics devices used to support high speed rotating machinery. The advantages that they offer to process fluid lubricated machines include: high rotational speed capability, no auxiliary lubrication system, non-contacting high speed operation, and improved damping as compared to rigid hydrodynamic bearings. NASA has had a sporadic research program in this technology for almost 6 decades. Advances in the technology and understanding of foil journal bearings have enabled several new commercial products in recent years. These products include oil-free turbochargers for both heavy trucks and automobiles, high speed electric motors, microturbines for distributed power generation, and turbojet engines. However, the foil thrust bearing has not received a complimentary level of research and therefore has become the weak link of oil-free turbomachinery. In an effort to both provide machine designers with basic performance parameters and to elucidate the underlying physics of foil thrust bearings, NASA Glenn Research Center has completed an effort to experimentally measure the performance of simple gas foil thrust bearing in air. The database includes simple bump foil supported thrust bearings with full geometry and manufacturing techniques available to the user. Test conditions consist of air at ambient pressure and temperatures up to 500 C and rotational speeds to 55,000 rpm. A complete set of axial load, frictional torque, and rotational speed is presented for two different compliant sub-structures and inter-pad gaps. Data obtained from commercially available foil thrust bearings both with and without active cooling is presented for comparison. A significant observation made possible by this data set is the speed-load capacity characteristic of foil thrust bearings. Whereas for the foil journal bearing the load capacity increases linearly with rotational speed, the foil thrust bearing operates in the hydrodynamic high speed limit. In this case, the load capacity is constant and in fact often decreases with speed if other factors such as thermal conditions and runner distortions are permitted to dominate the bearing performance.

Bruckner, Robert J.↗