Multi-Temperature, Thermal & Ion Fraction Effects over Wedge and Bluff Body Shapes in Hypervelocity Flow. Ward W. Vuillemot, Uri Shumlak

Transcription

1 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

2 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

3 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

4 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

5 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

6 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

7 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

8 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

9 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

10 Wedge Flow Region 1 (Freestream) Conditions Mach #: 3 Pressure: 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 Computational d University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000 δ θ

11 Wedge Flow : Computational Results 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

12 Bluff Body Flow Nitrogen (N 2 ) gas bluff body flow simulated at 50,000 feet standard atmosphere (density = 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

13 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

14 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 Drag versus Mach # Standard Atmosphere Conditions , , , , , Normalized Drag Mach # 3.0 1, ,000 5, : 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

15 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

16 Bluff Body Flow : Subsonic rho Density Contours Velocity Profile Streamtrace Mach 0.25 Mach rho University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000

17 Bluff Body Flow : Sonic & Supersonic rho Density Contours Velocity Profile Streamtrace Mach 1.0 Mach rho University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000

18 Bluff Body Flow : Mach 3 Comparison Density Contours Velocity Profile Streamtrace rho Mach 3.0, T = 5000 K & T = 217 K Ion Fraction Contours Velocity Profile Streamtrace fi Mach 3.0, T = 5000 K 18 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000

19 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

20 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

21 Poster Talk Request You can download a PDF version of this poster talk. It can be found at: < vuillemot,_ward/icopsconf.pdf> 21 University of Washington - Computational Fluid Dynamics Laboratory ICOPS2000