Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Modeling Snowflake Divertors in MAST-U Tokamak

We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗
DOE OSTI2023

The transport dynamics of tens of micrometer-sized water droplets in RF atmospheric pressure glow discharges

Abstract Charging of particles having diameters of tens of microns has been extensively studied at atmospheric pressure in the context of, for example, electrostatic precipitators where the focus was on unipolar charging. The ambipolar charging of particles in atmospheric pressure plasmas, and of droplets in particular, has received less attention. The plasma activation of droplets is of interest for water purification, fertilizer production and materials synthesis, all of which depend on the transport of the droplets through the plasma, which in turn depends on their charging. In this paper, we report on the transport dynamics of water droplets, tens of microns in diameter, carried by the gas flow through an atmospheric pressure radiofrequency glow discharge sustained in helium. The droplets pass through the plasma with minimal evaporation and without reaching the Rayleigh limit. The droplet trajectory in the presence and absence of the plasma provides insights on the forces acting on the droplet. The measurements were analyzed using results from a three-dimensional fluid model and a two-dimensional plasma hydrodynamics model. We found that the transport dynamics as the droplet enters and leaves the plasma are due to differential charging of the droplet in the plasma gradients of the bounding sheaths to the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Acoustic responses of coupled fluid-structure system by acoustic-structural analogy

The use of an analogy between structural mechanics and acoustics makes it possible to solve fluid-structural interaction (FSI) problems using an existing structural analysis computer program. This method was implemented in MSC/NASTRAN program and the FSI analysis was performed using two dimensional coupled fluid beam model to assess and evaluate the adequacy of this approach. The coupled modal analysis of 3-D model is also briefly discussed. The normal mode, modal frequency response and transient response analysis of 2-D coupled fluid beam system is presented. The significant reduction of the acoustic pressure response at the fluid structure interface is observed as a result of fluid structure interaction.

Shin, Y. S.↗
NASA NTRS1990

Investigations of a turbulence-driven solar wind model

This work presents an investigation of the properties of the one-dimensional, two-fluid turbulence-driven solar wind model introduced by Hollweg and Johnson (1988). It is found that the model has serious difficulties in reproducing the observed high-speed wind at 1 AU. In particular, the model proton temperatures are lower than those observed by a factor of 2 or more, and the highest temperature models yield excessive wave intensities at 1 AU. It appears that the problem stems from the specific spatial distribution of heat deposition in the model. Thus this study does not rule out a turbulence-driven fast solar wind, since other forms of the turbulent evolution could probably achieve better results. A three-fluid version of the model is also presented to show that the addition of alpha particles does not significantly reduce the extreme proton temperatures displayed near the sun by the Hollweg and Johnson work. Finally, it is suggested that the additional heating needs to be located well beyond the critical point, implying that the heating mechanism for the fast solar wind is likely not the same as that heating the solar corona.

Isenberg, Philip A.↗
NASA NTRS1994

Standing shocks in a two-fluid solar wind

We present a numerical study of the formation of standing shocks in the solar wind using a two-fluid time-dependent model in the presence of Alfven waves. Included in this model is the adiabatic cooling and thermal conduction of both electrons and protons. In this study, standing shocks develop in the flow when additional critical points form as a result of either localized momentum addition or rapid expansion of the flow tube below the existing sonic point. While the flow speed and density exhibit the same characteristics as found in earlier studies of the formation of standing shocks, the inclusion of electron and proton heat conduction produces different signatures in the electron and proton temperature profiles across the shock layer. Owing to the strong heat conduction, the electron temperature is nearly continuous across the shock, but its gradient has a negative jump across it, thus producing a net heat flux out of the shock layer. The proton temperature exhibits the same characteristics for shocks produced by momentum addition but behaves differently when the shock is formed by the rapid divergence of the flow tube. The adiabatic cooling in a rapidly diverging flow tube reduces the proton temperature so substantially that the proton heat conduction becomes negligible in the vicinity of the shock. As a result, protons experience a positive jump in temperature across the shock. While Alfven waves do not affect the formation of standing shocks, they contribute to the change of the mmomentum and energy balance across them. We also find that for this solar wind model the inclusion of thermal conduction and adiabatic cooling for the elctrons and protons increases significantly the range of parameters characterizing the formation of standing shocks over those previously found for isothermal and polytropic models.

Habbal, Shadia R.↗
NASA NTRS2012

Measurements and Simulations of Surface Dielectric Barrier Discharges Used as Plasma Actuators

This report is a Ph.D. dissertation performed under NRA cooperative agreement and submitted as part of the final report. Asymmetric surface dielectric barrier discharges (DBDs) have shown promise for use as aerodynamic actuators for active flow control. In this project we studied DBD actuators experimentally and numerically. Our DBDs used a symmetric triangular high voltage waveform to generate plasma in atmospheric pressure air. Time-averaged measurements indicated that the induced force of a single barrier actuator design (one electrode insulated from the plasma) can be increased exponentially above the results of previous studies by decreasing both the length and thickness of the electrode exposed to the plasma. This increased force may allow these devices to control flow separation in a wider range of flow environments. Experiments using an intensified digital camera to examine the plasma on time scales of a few nanoseconds showed that, in addition to the previously-observed filamentary and jet-like plasma structures, discharges with very thin exposed electrodes exhibited a weak but constant plasma immediately adjacent to those electrodes. In double-barrier actuators (both electrodes insulated), decreasing the diameter of the narrower electrode lead to increasing forces, and recorded images showed the simultaneous existence of both filamentary and jet-like plasma structures. The development and application of a time-dependent, two-dimensional computational fluid plasma model has aided in understanding the detailed physics of surface DBDs at all-time scales. For simulated single-barrier discharges, the model qualitatively reproduced the filamentary and jet-like micro-discharge structures. The model was somewhat successful in reproducing the observed characteristics of double-barrier actuators. For both actuator geometries, the model indicated that the majority of the forces induced on the neutral gas occur in between micro-discharges as the plasmas decay.

Hoskinson, Alan R.↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

Filling a tank with cryogenic fluid is more challenging than filling a tank with water or any other fluid that is in liquid state at atmospheric condition. Filling a tank with cryogenic fluid is a two-step process. First the tank and the transfer line mustbe chilled. Liquid cryogens start flowing into the tank only after the tank and transfer lines are chilled to the fluid saturation temperature. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor, caused by heat transfer from the warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. Filling a cryogenic tank in the absence of gravity is more challenging because in a non-stratified environment liquid propellantmay not settle at the tank bottom as it does on earth. There is a strong possibility that liquid propellant may exit through the vent valve, which is typically located at the top of the tank to vent propellant vapor. There are several methodsof filling atank in space: Charge-Hold-Vent method [1], Vented-Chill / No-Vent-Fill (VCNVF)[2],and No Vent Fill with TVS assisted injector [3].This paper addressesnumerical modelingofthe No Vent Fill (NVF) processwhere a TVS (Thermodynamic Vent System) augmented injector was used with the vent valve closed during the entire filling process. NVF tests were conducted in the CRYOgenic Orbital Testbed (CRYOTE) tankat NASA/Marshall Space Flight Centerin 2018.The Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code [4] was used to simulate the chilland fill process.

Cryogenics↗
NASA NTRS1994

Acceleration and heating of two-fluid solar wind by Alfven waves

Earlier model studies of solar wind driven by thermal pressure and Alfven waves have shown that wave amplitudes of 20-30 km/s at the coronal base are sufficient to accelerate the flow to the high speeds observed in quasi-steady streams emanating from large coronal holes. We focus on the energy balance in the proton gas and show that heat conduction from the region where the waves are dissipated may play an important role in determining the proton temperature at the orbit of Earth. In models with 'classical' heat conduction we find a correlation between high flow speed, high proton temperature, and low electron temperature at 1 AU. The effect of wave heating on the development of anisotropies in the solar wind proton gas pressure is also investigated in this study.

Sandbaek, Ornulf↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
NASA NTRS2018

kL-based Linear and Nonlinear Two-Equation Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) two-equation turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modification to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important differences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady flows in URANS mode. The present report documents the formulation of two model levels of turbulence models as implemented in the computational fluid dynamics FUN3D code. The levels are the two-equation linear k-kL and the two-equation algebraic Reynolds stress model (ARSM). Free shear, separated and corner flow cases are documented and compared with experimental, and other turbulence model data. The results show generally very good comparisons with experimental data. The results from this formulation are similar or better than results using the SST two-equation turbulence model. ARSM shows great promise with similar level of computational resources as basic two-equation turbulence models.

Abdol-Hamid, K. S.↗
NASA NTRS1998

Analysis and Modeling of a Two-Phase Jet Pump of a Thermal Management System for Aerospace Applications

Jet pumps are devices capable of pumping fluids to a higher pressure by inducing the motion of a secondary fluid employing a high speed primary fluid. The main components of a jet pump are a primary nozzle, secondary fluid injectors, a mixing chamber, a throat, and a diffuser. The work described in this paper models the flow of a two-phase primary fluid inducing a secondary liquid (saturated or subcooled) injected into the jet pump mixing chamber. The model is capable of accounting for phase transformations due to compression, expansion, and mixing. The model is also capable of incorporating the effects of the temperature and pressure dependency in the analysis. The approach adopted utilizes an isentropic constant pressure mixing in the mixing chamber and at times employs iterative techniques to determine the flow conditions in the different parts of the jet pump.

Sherif, S.A.↗
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 OSTI2020

Numerical simulation of streamer evolution in surface dielectric barrier discharge with electrode-array

Atmospheric pressure surface dielectric barrier discharges (SDBDs) may be composed of streamers fast propagating along a dielectric surface in ambient air, producing reactive oxygen and nitrogen species, and inducing a force on the neutral gas, which can find applications, respectively, in plasma medicine and aerodynamics. In this work, a two-dimensional self-consistent fluid model was developed to study SDBDs with an electrode-array. Emphasis was placed on the interaction of counter-propagating streamers and discharge uniformity for different applied voltages and geometric configurations of the electrode-array. When two counter-propagating streamers collide, the streamers come to a stop within a certain (ultimate) distance between the streamer heads. Optimizing the applied voltages is a convenient way to improve uniformity, making the streamer heads reach a minimum distance between each other. Shortening the electrode spacing can simultaneously shorten the streamer length and the ultimate distance between streamers. Under certain discharge parameters (such as applied voltage), there exists an optimum electrode spacing maximizing uniformity.

42 ENGINEERING↗
NASA NTRS2004

Real-Time Measurements of Aft Dome Insulation Erosion on Space Shuttle Reusable Solid Rocket Motor

Real-time erosion of aft dome internal insulation was measured with internal instrumentation on a static test of a lengthened version of the Space Shuffle Reusable Solid Rocket Motor (RSRM). This effort marks the first time that real-time aft dome insulation erosion (Le., erosion due to the combined effects of thermochemical ablation and mechanical abrasion) was measured in this kind of large motor static test [designated as Engineering Test Motor number 3 (ETM3)I. This paper presents data plots of the erosion depth versus time. The data indicates general erosion versus time behavior that is in contrast to what would be expected from earlier analyses. Engineers have long known that the thermal environment in the aft dome is severe and that the resulting aft dome insulation erosion is significant. Models of aft dome erosion involve a two-step process of computational fluid dynamics (CFD) modeling and material ablation modeling. This modeling effort is complex. The time- dependent effects are difficult to verify with only prefire and postfire insulation measurements. Nozzle vectoring, slag accumulation, and changing boundary conditions will affect the time dependence of aft dome erosion. Further study of this data and continued measurements on future motors will increase our understanding of the aft dome flow and erosion environment.

McWhorter, Bruce↗
DOE OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗