Electrode design for lift augmentation and power generation of atmospheric entry vehicles during aerocapture and entry, descent, and landing maneuvers

Abstract

A magnetohydrodynamic (MHD) flow control mechanism is described which substantially improves the existing processes in that smaller magnetic fields, requiring far less mass, may be placed away from the forebody of the spacecraft to produce Lorentz forces that augment the lift and the drag forces for guidance, navigation, and control of the spacecraft. The MHD flow control mechanism may also be configured to provide additional thermal protection of the electrodes therein.

Claims

1. A magnetohydrodynamic (MHD) flow control system for use with an aeroshell of a spacecraft comprising: at least a first pair of electrodes configured to be embedded in a first predetermined portion of the aeroshell, wherein the at least a first predetermined portion of the aeroshell is a channel having multiple side walls and a bottom portion connecting the multiple side walls having an outward facing surface and an inward facing surface, and further wherein each electrode is embedded within one of the multiple side walls of the channel; and a magnet configured for placement on an inward facing surface of the bottom of the channel.

2. The system of claim 1, wherein the magnet is selected from the group consisting of a permanent magnet, an electromagnet and combinations thereof.

3. The system of claim 1, wherein the electrodes are formed of at least one of tungsten and a tungsten alloy.

4. The system of claim 1, wherein the channel includes slanted side walls.

5. The system of claim 4, wherein each electrode of the first pair of electrodes is embedded in a slanted side wall.

6. The system of claim 4, wherein the channel further includes a flat bottom portion connecting the slanted side walls.

7. The system of claim 6, wherein the magnet is located on an inward facing surface of the aeroshell opposite the flat bottom portion of the channel.

8. The system of claim 1, wherein the aeroshell is a blunt-body aeroshell.

9. An aeroshell for use with a spacecraft comprising: a blunt-body configuration having at least one channel formed therein, wherein the at least one channel has multiple sides and a bottom connecting the multiple sides having an outward facing surface and an inward facing surface; at least a first pair of electrodes embedded in at least two of the multiple sides of the at least one channel; a magnet placed on an inward facing surface of the bottom of the least one channel.

10. The aeroshell of claim 9, wherein the channel includes two slanted sides and a flat bottom connecting the two slanted sides.

11. The aeroshell of claim 10, wherein the electrodes of the at least a first electrode pair are embedded in the two slanted sides, and further wherein the electrodes are located opposite one another across the flat bottom.

12. The aeroshell of claim 9, wherein the magnet is placed opposite the flat bottom of the channel.

13. The aeroshell of claim 9, wherein the magnet is selected from the group consisting of a permanent magnet, an electromagnet, and combinations thereof.

14. The system of claim 9, wherein the electrodes are formed of at least one of tungsten and a tungsten alloy.

15. An aeroshell for use with a spacecraft comprising: multiple channels formed therein, wherein each of the multiple channels has multiple sides and a bottom connecting the multiple sides including an outward facing surface and an inward facing surface; multiple electrode pairs embedded in at least two of the multiple sides of one of the multiple channels; and at least one magnet associated with each of the multiple electrode pairs, the at least one magnet being placed on an inward facing surface of the bottom of the one of the multiple channels containing the multiple electrode pairs embedded therein.

16. The aeroshell of claim 15, wherein each of the multiple channel includes two slanted sides and a flat bottom connecting the two slanted sides.

17. The aeroshell of claim 16, wherein the electrodes of each of the multiple electrode pairs are embedded in the two slanted sides of at least one of the multiple channels, and further wherein the electrodes are located opposite one another across the flat bottom.

18. The aeroshell of claim 15, wherein the at least one magnet is placed opposite the flat bottom of the channel.

19. The aeroshell of claim 15, wherein each of the multiple electrode pairs is embedded in a different one of the at least one multiple channels.

Description

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS

(1) FIG. 1 is a prior art illustration of the desired result from implementation of an embodiment described herein;

(2) FIG. 2 is a prior art illustration summarizing the physics principles behind implementation of an embodiment described herein;

(3) FIGS. 3a, 3b, 3c, 3d, 3e illustrate various views of an MHD patch implemented in accordance with an embodiment described herein;

(4) FIG. 4 shows the Lorentz Forces (F.sub.x, F.sub.y, F.sub.z) on a blunt body aeroshell in the presence of ionized flow in accordance with an embodiment described herein;

(5) FIG. 5 shows effects of magnetic field orientation on the Lorentz forces and energy values of the MHD approach applied to a specific use case in accordance with an embodiment described herein;

(6) FIG. 6 shows the resulting Lorentz forces in accordance with the magnetic field orientations of FIG. 5;

(7) FIG. 7 illustrates the effect of distance between electrodes on the current density streamlines in accordance with an embodiment described herein;

(8) FIG. 8 shows the resulting force values in accordance with the distance between electrodes of FIG. 7;

(9) FIGS. 9a, 9b, 9c are schematics of a prior art mid L/D aeroshell used in the analysis of the MHD approach of the embodiments described herein;

(10) FIG. 10 illustrates the location of the mid L/D aeroshell which was selected for analysis of the MHD approach of the embodiments described herein as applied to a first specific use case; and

(11) FIG. 11 illustrates the location of the mid L/D aeroshell which was selected for analysis of the MHD approach of the embodiments described herein as applied to a second specific use case; and

(12) FIG. 12 shows temperature contours of the flow near the current streamlines created by the electrodes and magnetic field in accordance with the embodiments described herein.

DETAILED DESCRIPTION OF THE EMBODIMENTS

(13) The MHD approach in the present embodiments substantially improves the existing processes in that smaller magnetic fields, requiring far less mass, are placed away from the forebody of the spacecraft to produce Lorentz forces that augment the lift and the drag forces for guidance, navigation, and control of the spacecraft, as explained herein. The embodiments described below facilitate lift generation for enabling the use of blunt bodies and other aeroshell shapes to meet stringent entry parameters for several mission destinations of interest to NASA, including Neptune and Mars. More particularly, the embodiments herein provide for a (MHD) flow control system which provides for additional thermal protection of the electrodes.

(14) It has been documented and well understood how the principles of MHD could be applied to entry bodies to reduce thermal stresses on the aeroshell as well as to augment drag and to harvest atmospheric gases and to generate electrical power.

(15) MHD Theory for EDL Applications

(16) As shown in the schematic of FIG. 2, the principles of MHD stem from the theory of electrical currents flowing in the presence of a magnetic field. In this case, the electrical current flow is created when the conductivity of the planet's (or other celestial body's) atmosphere is sufficient to facilitate current flow between electrodes with opposite polarity that are placed on the outer surface of the spacecraft. Hence, the conductivity of the flow is critical to the use of this MHD approach during atmospheric entries.

(17) When electrically conducting fluid crosses magnetic field lines, charged particles experience Lorentz forces {right arrow over (F)}.sub.L=q{right arrow over (u)}{right arrow over (B)}, and since the charges of electrons and ions are of opposite signs, the electrons and ions are pulled apart, which creates a Faraday electromotive force (emf). If the circuit is closed, current flows:
j=(1K)uB.sub.z(1)
Here is the electrical conductivity, and K is the load factor (load resistance vs plasma resistance). Interaction of the induced current with the B field creates body force per unit area:
{right arrow over (F)}={right arrow over (j)}{right arrow over (B)}(2)
If the B field is inclined with respect to the surface, the force would have both decelerating(drag) and normal (lift) components:
F.sub.x=(1K)uB.sub.z.sup.2(3)
F.sub.z=(1K)uB.sub.zB.sub.x(4)

(18) The scalar electrical conductivity is proportional to the number density n.sub.e of principal charge carriers (electrons) and inversely proportional to the rate, or frequency, of their collisions with ions (v.sub.ei) and neutral molecules (v.sub.en):

(19) = e 2 n e m ( v e i + v e n ) ( 5 )

(20) Here e and m are the electron charge and mass, respectively.

(21) At low ionization fraction

(22) n e n 1 0 - 3 - 1 0 - 2
(n is the number density of gas), electrons collide mostly with neutral molecules, and the conductivity is proportional to the ionization fraction:

(23) 3 1 0 5 n e n S / m ( 6 )

(24) In this regime, the conductivity in hypersonic shock and boundary layers reaches 10-100 S/m.

(25) At high ionization fraction,

(26) n e n 1 0 - 3 - 1 0 - 2 ,
electrons collide mostly with ions, and the conductivity does not depend on the ionization fraction, instead being determined by the electron temperature:
constT.sub.e.sup.3/2(7)

(27) In this regime, the conductivity reaches 1000-3000 S/m.

(28) Hall Effect

(29) With the induced electric field E perpendicular to the B field, electrons and ions have {right arrow over (E)}{right arrow over (B)} drift motion along the flow, which effectively diverts their cross-low motion and thus reduces MHD body forces. This can be expressed as an effective conductivity lower than the scalar one.

(30) Ion Slip Effect

(31) MHD effects increase with flow velocity. However, as B field increases, the motion of electrons and ions is retarded so much that the ion-electron fluid starts to slip against the bulk neutral gas, and the lower ion-electron velocity reduces the MHD effects (including body forces). This ion slip effect can also be expressed as a reduction in the effective conductivity.

(32) eff = ~ 1 + ~ 2 ( 8 )
Here the conductivity {tilde over ()} and the Hall parameter {tilde over ()} corrected for ion slip are:

(33) ~ = 1 + e i ; ~ = e 1 + e i ( 9 )
Here is the scalar conductivity, and the electron and ion Hall parameters are:

(34) e = e B m ( v e i + v e n ) ; i = e B M v i n ( 10 )
(M and v.sub.in are the ion mass and ion-neutral collision frequency, respectively).

(35) The performance-reducing Hall and ion slip effects become dominant when the Hall parameter (the product of the electron mobility and the magnetic field) is high, which occurs when the magnetic field is high and the gas density is low. However, these effects can be exploited for modulating axisymmetric forces and moments, as will be illustrated later in one case presented for Mars.

(36) For the embodiments described herein, NASA Langley's CFD code LAURA is used to calculate properties of the flow field around an entry aeroshell. A description of LAURA may be found in Thompson, K., et al., LAURA User's Manual: 5.6, NASA TM 2020-220566, 2020. These calculated properties are then fed into MHD plasma code CFDWARP to calculate MHD properties of the flow and estimate Lorentz forces generated and electrical power available on a local region of the aeroshell during that entry condition. The CFDWARP code has been described in the prior art at, for example, B Parent, et al., Modeling Weakly-Ionized Plasmas in Magnetic Field: A New Computationally-Efficient Approach, Journal of Computational Physics, Vol. 300, Pages 779-799, 2015, which is incorporated herein by reference in its entirety.

(37) In one embodiment, an MHD patch effector device is integrated with a blunt-body aeroshell for generating Lorentz forces as shown in FIGS. 3a, 3b, 3c, 3d and 3e. The MHD patch effector may consist of: 1) a pair of electrodes 10a, 10b embedded in the surface of the aeroshell 5 at a prescribed distance apart, wherein the electrodes are angled from facing directly toward the shock layer closest to the electrode in a conformal recessed configuration (FIGS. 3d, 3e) and 2) a magnet 15 placed just below the aeroshell's outer surface, with the magnetic field (B) protruding out of the surface. The magnet 15 could be a permanent magnet(s) or electromagnet(s), e.g., coil, self-powered by the electrical energy harvested from the flow or powered by another energy source. FIG. 3a shows the front view of the aeroshell 5 with (+/) electrodes 10a, 10b embedded therein. FIG. 3b is a side, cross-sectional view of FIG. 3a, wherein the aeroshell 5 with electrode 10b is shown, as well as the magnet 15 on the inside of the aeroshell, opposite the electrodes on the outside of the aeroshell. FIG. 3b also includes direction indicators for the magnetic field (B), Flow of air, plasma around the aeroshell and electrodes (Flow) and spaceship velocity vector vs and force (F). FIGS. 3c, 3d and 3e offer alternative views of the MHD patch effector device with directional indicators for magnetic field (B), (Flow), force (F) and current density (J), wherein: F=JB.

(38) FIG. 4 shows the Lorentz Forces (F.sub.x, F.sub.y, F.sub.z) on a blunt body aeroshell in the presence of ionized flow. For an initial case, the MHD results are as follows. Force on the gas by the electromagnetic fields (JB)=(2718, 103, +103) Newtons (N). That represents a force of 2718 N in the direction opposite of the bulk flow, a force of 103 N that acts normal to the surface that is pressing on the surface, and a force of 103 N that acts tangential to the surface and perpendicular to the Bulk flow direction and perpendicular to the surface normal direction. Whether the force in the X-axis acts as drag depends on the orientation of the modeled aeroshell surface to the velocity vector of the spacecraft. Likewise, whether the force in the Y-axis acts as lift and the force in the Z-axis acts as a rolling moment depends on the orientation of the modeled aeroshell surface to the velocity vector of the spacecraft. For instance, if the modeled aeroshell is located on the underbelly of the mid L/D vehicle (see FIG. 10), then these X-, Y-, and Z-forces can be considered drag, lift, and roll, respectively. However, if these X-, Y-, and Z-forces are acting along an aeroshell representative of a 70-degree cone, then the X-axis force is acting more like a lift force than a drag force, and the Y-axis force is acting more like a drag force than a lift force, while the Z-axis force is representative of a rolling force. Further to the initial example, the rate of energy input to the gas by the electromagnetic fields (E*J)=3.01 Megawatts; the rate of work done by the electromagnetic fields (V*(JB))=29.2 Megawatts; and the rate of heat deposited to the gas by electromagnetic fields (Energy inputwork done)=26.2 Megawatts.

(39) Additional cases were run for understanding the effects of magnetic field orientation on the Lorentz Forces and Energy values of the MHD approach applied to the Neptune entry case of 29.24 km/s at an atmospheric density of 1.45e4 kg/m3. A value of 1 Tesla was used for these additional cases. Shown in FIG. 5, the orientation of the magnetic field varied in the X-Y plane for 180 degrees (=90 to +90 degrees). The magnetic field remained orthogonal to the Z-axis for these additional cases.

(40) The resulting Lorentz Forces are illustrated in FIG. 6 for the values of 0 between 90 degrees and +90 degrees. The force in the Y-axis (perpendicular to the aeroshell surface) is maximal when =+/45 degrees (F.sub.emfield,y). The force in the X-axis (bulk flow direction) is maximal when is +/90 degrees (F.sub.emfield,x).

(41) A force of 4000 Newtons in these cases translates in an electromagnetic force per surface area of 200 kN/m2. Therefore, an MHD patch effector having an area of 1 m.sup.2 would produce a force of 2.0e+5 N. That value is the same order of magnitude as the whole body drag and lift forces computed by LAURA for the Neptune entry velocity and atmospheric density. The math comes down to the sizing and orientation of the magnetic field to produce the desire drag and lift vectors for controlling the entry body. The entry body could be a number of shapes, possibly including a 70-degree cone.

(42) The distance between the two electrodes (D) (See FIG. 3a) was varied between 0.07 meters and 0.97 meters. Zero current boundary conditions were imposed everywhere except at the two electrodes. Both electrodes were given the same voltage (short circuit). Flow properties are frozen to represent no coupling. Hall effect and ion slip effects are taken into consideration.

(43) The effect of distance between electrodes on the current density streamlines is illustrated in FIG. 7. The longer streamlines (S.sub.L) represent the case for D=0.97 meters spacing while the shorter streamlines (S.sub.S) represent the case for D=0.07 meters.

(44) The effect of electrode spacing (D) (See FIG. 3a) on Lorentz forces is shown in FIG. 8. Force increases when the electrode spacing is increased. This effect is due to longer length of the current streamlines being more perpendicular to the flow and the magnetic field as D is increased. For instance, the shape of the shorter streamlines S.sub.S for the shorter spacing look more circular when compared to the much longer streamlines S.sub.L for the wider spacing. Hence, a longer portion of the S.sub.L streamlines, which represent J, are orthogonal to the B-field than the portion of the shorter streamlines S.sub.S.

(45) The force values (shown in FIG. 8) reach up to 800 kN/m2. Hence, a MHD patch that encompasses 1 m.sup.2 of the aeroshell can, theoretically, produce up to 800 kN force. That value is not only the same order of magnitude but also exceeds the value computed by LAURA. Hence, the L/D ratio exceeds 1.0 for this situation. This is a ground-breaking discovery for Neptune missions. And performance during atmospheric entries at Earth, such as return from the Moon or return from Mars, is expected to be similar to this Neptune case.

(46) The following paragraphs provide further support of the description and embodiments discussed above with respect to particular missions, i.e., Neptune and Mars. Those skilled in the art with the benefit of this disclosure will appreciate that the examples are not intended to limit the scope of the disclosure, but rather to provide illustrative examples of certain embodiments.

(47) The criteria used to select a region of flow along the spacecraft for further analysis may be derived from the principles of MHD illustrated above in FIG. 2. Conductivity may be deemed the most important criteria in this example. As explained above, conductivity is dependent on temperature and density conditions of the flow field near the surface of the vehicle. A second criteria may be the velocity of the flow in terms of magnitude and direction for establishing the location and placement of the electrodes. A third criteria may be the geometric orientation of the magnetic (B-) field with respect to the current flow to produce Lorentz forces with axis components in the orthogonal coordinate system that are in the drag direction and lift direction and rolling moment direction.

(48) The LAURA Navier-Stokes solver was applied for the baseline Mars and Neptune flowfields considered in the embodiments described herein. Two-temperature thermochemical nonequilibrium was assumed, where the chemical kinetics for the Mars case are taken from Johnston and Brandis, Modeling of Nonequilibrium CO Fourth-Positive and CN Violet Emission in CO2N2 Gases, Journal of Quantitative Spectroscopy and Radiative Transfer, Vol. 149, 2014, pp. 303-317 and for the Neptune case from a combination of Park, C. Nonequilibrium Ionization and Radiation in Hydrogen Helium Mixtures, Journal of Thermophysics and Heat Transfer, Vol. 26, No. 2, 2012, pp. 231-243, Gocken, T., N2-CH4-Ar Chemical Kinetic Model for Simulations of Atmospheric Entry to Titan, Journal of Thermophysics and Heat Transfer, Vol. 21, No. 1, 2007, pp. 9-1, Fujita K, Yamada T, and Ishii N., Impact of Ablation Gas Kinetics on Hyperbolic Entry Radiative Heating, AIAA Paper 2006-1185, 2006, and Johnston and Brandis. All references are incorporated herein by reference. The wall was assumed fully catalytic and in radiative equilibrium. The solution and grid convergence criteria applied were consistent with the state-of-the-art applied for NASA flight programs.

(49) For an example Neptune case, the team reviewed the aerothermal analysis and other mission results of the 2004 Neptune study for the mid L/D aeroshell shown in FIGS. 9a (top view), 9b (side view) and 9c (aft view). The location selected LA.sub.N for further analysis is shown in FIG. 10 based on the criteria stated above. The entry conditions used in the 2004 Neptune study were selected.

(50) For a Mars example, the team reviewed several mission cases and selected a blunt body aeroshell. The location LA.sub.M selected for further analysis is shown in FIG. 11 based on the criteria stated above. However, for Mars, a higher entry velocity representing a faster flight to Mars was selected.

(51) For the Neptune entry, a two-temperature thermochemical nonequilibrium flow field was modeled in LAURA for the geometry from Edquist, Karl T., et al., Configuration, Aerodynamics, and Stability Analysis for a Neptune Aerocapture Orbiter, AIAA 2004-4953, AIAA Atmospheric Flight Mechanics Conference and Exhibit, 16-19 Aug. 2004, Providence, Rhode Island, L/D=0.8, 40 degrees Angle of Attack, 29.24 km/s entry velocity at an atmospheric density of 1.45e4 kg/m3. Free-stream mass fractions of 0.6246 H2, 0.2909 He, and 0.0846 CH4 were used. The species included in the simulation are H2, H, H+, He, He+, e, CH4, CH3, CH2, CH, C2, C, and C+. Temperature, pressure, and electron number densities were calculated. For the 29 km/s Neptune case, the whole body drag=5.6487e+05 N and lift=4.4911e+05 N. An additional case was run in LAURA for comparison to the 2004 case. For the 35 km/s Neptune case, the whole body drag=8.3943e+05 N and lift=6.5381e+05 N.

(52) For the Mars entry, a two-temperature thermochemical nonequilibrium flow field was modeled in LAURA to represent a HIAD at 0 degrees Angle of Attack, 9 km/s entry velocity at an atmospheric density of 4.4e-5 kg/m3. Free-stream mass fractions of 0.97 CO2 and 0.03 N2 were used. The species included in the simulation are CO2, N2, CO, NO, O2, CN, C2, C, N, O, N2+, CO+, NO+, C+N+, 0+, and e. Temperature, pressure, and electron number densities were calculated. The HIAD geometry consists of a 70-degree sphere-cone with a 10-meter nose radius, 17.2-meter maximum diameter, and an 0.35-meter shoulder radius. The whole body drag value of 7.0887e+05 N was calculated by LAURA for the Mars case at 9.0 km/s, 4.4e5 kg/m3. Additional Mars cases were run in LAURA for comparison.

(53) At 7.5 km/s, 4.4e5 kg/m3: Drag=4.9227e+05 N. At 9.0 km/s, 4.4e4 kg/m3: Drag=7.0887e+06 N. At 7.5 km/s, 4.4e4 kg/m3: Drag=4.9227e+06 N.

(54) FIG. 12 illustrates the temperature contours of the flow near the current streamlines created by the electrodes and magnetic fields in certain electrode arrangements, also described in co-pending U.S. patent application Ser. No. 17/727,009, entitled System and Method for Lift Augmentation of Atmospheric Entry Vehicles During Aerocapture and Entry, Descent, and Landing Maneuvers. The present embodiments address possible thermal interruptions by provide larger physical separation between the electrodes and the shock layer by lowering the electrodes into the aeroshell surface while also angling the electrodes from facing directly toward the portion of the shock layer closest to the electrode as shown in FIGS. 3a, 3b, 3c, 3d and 3e. In addition to changes in the electrode configuration, different electrode materials may be selected in accordance with mission requirements and expectations.

(55) Lowering the electrodes into a channel configuration the protrudes into the aeroshell structure may add volume to that region of the spacecraft to accommodate their integration. This may lead to additional heat shield and aeroshell designs that accommodate multiple channels at clocked radii originating from the stagnation region (center) of the aeroshell to the outer shoulder of the aeroshell. The length of the channel can depend on the depth of the recession and size of the electrodes. Further, one skilled in the art with the benefit of this disclosure will appreciate that multiple electrode pairs and corresponding magnets may be used in a single channel or in multiple channels in order to provide additional guidance, navigation, and control of the spacecraft.

(56) One skilled in the art will appreciate that while the embodiments herein are generally directed to a single MHD patch, the embodiments are not so limited. Multiple patches having different configurations, e.g., dimensions, electrode spacing, electrode material, etc., may be utilized on an aeroshell in accordance with expected environment. Such alternatives are considered to be well within the scope of the embodiments.

(57) The performance of the MHD Patch for Earth entries such as aerocapture missions returning from the Moon or from Mars is expected to be very similar to the large lift forces shown for the Neptune cases illustrated herein.

(58) It is to be understood that the novel concepts described and illustrated herein may assume various alternative configurations, except where expressly specified to the contrary. It is also to be understood that the specific systems, devices and processes illustrated in the attached drawings, and described herein, are simply exemplary embodiments of the embodied concepts defined in the appended claims. Accordingly, specific dimensions and other physical characteristics relating to the embodiments disclosed herein are not to be considered as limiting, unless the claims expressly state otherwise.

(59) Reference in the specification to one embodiment or to an embodiment means that a particular element, feature, structure, or characteristic described in connection with the embodiments is included in at least one embodiment. The appearance of the phrases in one embodiment, in some embodiments, and in other embodiments in the specification are not necessarily all referring to the same embodiment or the same set of embodiments.

(60) As used herein, the terms comprises, comprising, includes, including, has, having or any other variation thereof, are intended to cover a non-exclusive inclusion. For example, a process, method, article, system or apparatus that comprises a list of elements is not necessarily limited to only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Further, unless expressly stated to the contrary, or refers to an inclusive or and not to an exclusive or. Additionally, use of the a or an are employed to describe elements and components of the embodiments herein. This is done merely for convenience and to give a general sense of the invention. This detailed description should be read to include one or at least one and the singular also includes the plural unless it is obviously meant otherwise.