Multi-Temperature, Thermal & Ion Fraction Effects over Wedge and Bluff Body Shapes in ypervelocity Flow Ward W. Vuillemot, Uri Shumlak 1 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Abstract We numerically investigate the effect of ionization on hypersonic flows by using an approximate Riemann fully three-dimensional MD solver, WARP3, which will include multiple temperature in the near future. The code calculates the ionization fraction as a function of temperature. We treat the flow as a single fluid with three (3) constituents, or namely ions, electrons, and neutrals. It is believed that the inclusion of multi-temperatures effects may explain the experimentally measured increase of shock stand-off distance encountered when an ionized hypersonic flow stagnates bluff body. 2 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Motivation Current MD code assumes a fully, singlely ionized fluid, where ion temperature is equal to electron temperature. While this simplifies the system of equations, it does not provide an adequate means to describe flow regimes where full ionization cannot be assumed. These said flow regimes include, but are not limited to: hypersonic vehicles similar to National Aerospace Plane ballistic re-entry vehicles advanced plasma thrusters for plasma thrusters portable, pulsed power systems 3 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Generalized MD Equations (0 <= f i <= 1) In order to modify WARP3 to include single fluid, multiple constituent flow, the single fluid equations have been re-derived to include ion fraction, f i, effects. We begin by assuming that ionization level, Z, is always of unity. Consequently, the number densities relations between ions, electrons, and neutrals are readily obtained. n ions = n electrons = f ions n n neutrals = (1 - f ions )n 4 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Non-dimensional, Conserved MD Equations Continuity Equation Conservation of Momentum Magnetic Pressure Energy Conservation ρ + ρv = 0 t ρv + (ρv V f t i B B + (p + f ib B ) I J = ρ ˆ 2 B t + f i( V B B V )= B ˆ res + B ˆ all + ˆ e t + e + p + B B V ( B V ) B = 2 ˆ B res + ˆ B all + ˆ B diamag B diamag ν visc 5 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Generalized Ohm s Law Following the prototypical methodology utilized to derive the generalized Ohm's law from two fluid equations, it can shown that the ion fraction effects are: E + U B = ne 1 LS E, electric field U, velocity field B, magnetic field D f i R ( 1 f m ions f i )N i RS, all term D, Diamagnetic term R, Resistive term (electron ion collision term) N, Neutral Collision term 6 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Generalized Ohm s Law (details) Resistive Neutral Collision all 1 Diamagnetic 1 η ( B ) B J Cn 2 m e u e J B P e m e m i 1 2 u i, C is a collision frequency constant, [L 3 t -1 ] 1: Not Supported in WARP3 due to dependence on small time scales in comparison to other terms computationally costly. Typically they are assumed to be too small to affect our solution 7 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Thermal Effects Explanation of Various Effects We can expect these effects to diffuse out, forcing transport of energy from energy-rich regions to energy-poor regions. We should not observe significant effects for hyperbolic systems (e.g. fluid speed greater than sonic speed.) Ion Fraction Effects An ion fraction less than one will adversely effect the effective magnetic field strength by the square root of the ion fraction. We can also expect the resistive diffusion to increase with a decrease in the ion fraction. We assume that the resistive diffusion will go to the vacuum resistive value asymptotically. 8 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Explanation of Various Effects (continued) Multi-Temperature Effects A model for a single fluid composed of three constituents (neutrals, ions, and electrons), will naturally have multiple methods of energy transfer for different heating mechanisms. eating Mechanisms & Associated Constituents Viscosity Neutrals Ions Resistivity & Radiation Electrons Thermal Conduction Electrons Neutrals ions 9 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Wedge Flow Region 1 (Freestream) Conditions Mach #: 3 Pressure: 0.714 Pressure: 1.0 γ: 1.4 Shock Flow 1 2 Wedge Region 2 Comparison of Analytical Results to Computational Results for Supersonic Flow over Wedge (δ = 25 ) δ θ Density X-velocity Y-velocity Pressure Analytical 25 44.5 2.796 2.279 0 3.517 Computational 25 44.3 2.790 2.276 1.0d-4 3.517 10 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000 δ θ
Wedge Flow : Computational Results 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 Mach 3.0, 25 degrees Density contours Velocity profiles streamtraces Flow enters at 25º angle from left and top boundaries, equivalent to a 25º wedge grid, but finite volume grid is simpler to generate 11 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow Nitrogen (N 2 ) gas bluff body flow simulated at 50,000 feet standard atmosphere (density = 0.186 kg m -3, temperature = 217 K) for the following Mach numbers: 0.25; 0.5; 1.0; 1.5; 3.0; and 3.0. At this altitude, the sonic velocity is approximately 300 m s -1. A case of increased enthalpy at Mach 3.0 for temperature of 5000 K. At this Mach number and increased temperature, the sonic velocity is approximately 460 m s -1. Pressure is varied so that density matches standard atmosphere at 50,000 feet and temperature matches artificially high enthalpy The exact mechanism to induce this increased enthalpy is not of concern, however something similar to K. Chadwick et al tuned frequency waveguide might be applicable 12 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Generalized Grid Use of generalized, algebraic grid for our initial study helps validate results. Bow shock formation is not a consequence of the grid, but of physics. Future studies will use elliptic grid to minimize computational costs associated with the freestream flow region 13 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Drag Calculated for Standard Atmosphere at 50,000 feet unless otherwise noted Mach [ # ] Sonic Speed [ m s -1 ] Density [ kg m -3 ] Pressure [ Pa ] Thermal Energy [ K ] Drag 1 2.5 Drag versus Mach # Standard Atmosphere Conditions 0.25 75 0.186 11,600 217 2.10 0.5 150 0.186 11,600 217 1.94 1.0 300 0.186 11,600 217 1.00 1.5 450 0.186 11,600 217 1.22 3.0 900 0.186 11,600 217 2.01 Normalized Drag 2 1.5 1 0.5 0 0.25 0.5 1 1.5 3 Mach # 3.0 1,441 0.186 276,000 5,000 48 1 : Normalized to sonic drag value (Mach # = 1) Pressure Drag along the hemisphere is lowest at sonic velocity, where the bow shock is strongest 14 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Blast Theory Density contours red-green interface represent approximate stagnation region Yellow line represents Blast theory prediction Blast Theory Prediction at Mach 3.0 WARP3 results coincide with first-order approximation Blast Theory prediction Shock stand-off distance is chosen to match WARP3 results, as there is no theory that accurately predicts distance 15 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Subsonic rho 0.19 0.19 0.19 0.18 0.18 0.18 Density Contours Velocity Profile Streamtrace Mach 0.25 Mach 0.5 0.18 0.18 rho 0.20 0.20 0.19 0.18 0.18 0.17 0.17 0.16 16 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Sonic & Supersonic rho 0.26 0.24 0.22 0.20 0.18 0.15 Density Contours Velocity Profile Streamtrace Mach 1.0 Mach 1.5 0.13 0.11 rho 0.39 0.35 0.31 0.28 0.24 0.20 0.17 0.13 17 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Bluff Body Flow : Mach 3 Comparison Density Contours Velocity Profile Streamtrace rho 0.66 0.58 0.51 0.44 0.36 0.29 0.22 0.14 Mach 3.0, T = 5000 K & T = 217 K Ion Fraction Contours Velocity Profile Streamtrace fi 0.04 0.03 0.03 0.02 0.02 0.01 0.01 0.00 Mach 3.0, T = 5000 K 18 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Conclusions Present state of WARP3 has been benchmarked to analytic results for wedge flow and first-order approximation blast theory for bluff body flow WARP3 is a robust CFD tool that includes such effects as: resistive, viscous, and thermal diffusion, and ion fraction effects Relatively small increase in enthalpy results in significant gains in ion fraction at stagnation region In the future, we will combine codes so that WARP3 will be an explicit/implicit MD solver with three-constituent single-fluid flow 19 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Conclusions (continued) There is evidence from experiments both in the United States and former Soviet Union that an ionized medium can reduce the drag on hypersonic vehicles. These effects may be explained by differences in the three constituents' temperatures, and subsequent coupled, energy transport evolution. Multiple temperature and ion fraction effects coupled with electrically conducting fluid may results in many benefits in such diverse areas as: reduced wave drag reduced skin friction reduced control-effector cross-section and/or divergement angle vectorable thrust 20 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000
Poster Talk Request You can download a PDF version of this poster talk. It can be found at: <http://www.aa.washington.edu/cfdlab/publications/ vuillemot,_ward/icopsconf.pdf> 21 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000