Patent Grant
11724836
Not Applicable
The present invention pertains generally to safely return spacecrafts and space vehicles to earth's surface. The term “returnable” refers to the ability of spacecraft to leave and later reenter planetary atmospheres and reach the planet's surface intact and safe. More specifically, the invention pertains to control systems and methods of control for returnable space vehicles re-entering planetary atmosphere.
Atmospheric re-entry refers to the process by which vehicles outside of a planet's atmosphere can enter or reenter the atmosphere and reach the surface of the planet intact. The technology for atmospheric re-entry owes its origins to the development of intercontinental ballistic missiles (ICBMs) during the Cold War. Early calculations showed that the kinetic energy of a nuclear warhead returning from orbit was sufficient to completely vaporize the warhead. In order to prevent a spacecraft from suffering a similar fate, several obstacles had to be overcome. There are difficulties that are inherent aspects of atmospheric entry caused by high velocity transit through an atmosphere. These problems include high temperatures generated by friction between the entry system and the atmosphere while the atmospheric entry system is traveling through the atmosphere at high velocities, high acceleration loads experienced by the reentry system and its payload due to rapid changes in the system's velocity, and difficulty of controlling the flight path as the system descends through an atmosphere due to several causes including poor aerodynamics, high velocities, and the requirement of fast reacting guidance methods. Each of these problems can be mitigated by allowing the atmospheric transit to occur over a longer period of time to allow for better flight control and less extreme aeroheating.
The transportation of a space object through an atmosphere is not necessarily difficult. All that must occur is that the gravitational attraction between the object in space and the celestial body overcome the other forces acting on the object. This results in a ballistic entry into and through the atmosphere of the celestial body, followed by the object reaching the surface of the celestial body. However, entering and traversing an atmosphere ballistically is relatively uncontrolled and can be physically damaging to the object entering the atmosphere, its payload, or to the celestial body. The energy put into a spacecraft by the launch system, usually a powerful rocket, must be removed from the spacecraft before it can land without damage. The spacecraft starts at zero altitude and zero velocity. It is then accelerated to a high velocity by the launch system. For example, the velocity required to keep a spacecraft in low Earth orbit is about 18,000 miles per hour while a spacecraft on a direct return from the Moon will have a reentry velocity of about 25,000 miles per hour. The kinetic energy (Ek) of the returning spacecraft is equal to half the mass (m) of the spacecraft multiplied by the square of the spacecraft's velocity (v), expressed as Ek=½mv2. The potential energy (Ep) of a returning spacecraft is its altitude (h) multiplied by its mass (m) and a gravitational constant (g), expressed as Ep=mgh. Most of a spacecraft's kinetic energy must be removed during reentry by aerodynamic drag between the spacecraft and the atmosphere for the spacecraft to land safely. This is why reentering spacecraft must have effective thermal protection systems, “heat shields,” and even with effective thermal protection the temperature on the aerodynamic surface of the spacecraft can rise to thousands of degrees because all of the kinetic energy of the spacecraft's velocity is removed in only a few minutes of flight. This heating causes the plasma “fireball” that surrounds a reentering spacecraft.
Deceleration is one obstacle in atmospheric re-entry. The Earth's rotational velocity is approximately 1000 miles per hour at the equator, and slower as a cosine function at higher latitudes. When a spacecraft leaves Earth orbit and prepares for re-entry, it can be traveling as fast as 18,000 mph. From the time of re-entry to just before landing, the spacecraft must decelerate to match the Earth's rotational velocity as closely as possible. Deceleration is generally accomplished by relying on atmospheric friction and drag (i.e. using wind resistance), a technique known as aerobraking. Aerobraking, however, presents another obstacle. The air in the Earth's atmosphere is composed mostly of nitrogen and oxygen. When an object such as a spacecraft moves through the atmosphere at supersonic speeds, it generates shockwaves as it collides with particles in the atmosphere. Since the orbital velocity for most spacecraft is about 25 times the speed of sound, air molecules in the path of the spacecraft are “shock heated”, i.e. compressed so violently that the temperature of the shockwaves increases to several thousand degrees. To provide an example, meteorites entering the Earth's atmosphere are often vaporized by their own shock waves. Because of shock heating, re-entry space vehicles require a heat shield of some form, usually comprising ablative materials and highly effective insulators. The spacecraft must also have great structural strength to bear the deceleration mechanical loads generated by these forces which have been measured to be 6 to 40 gravities or more (60-400 m/s2)
Additionally, a spacecraft's speed and the resulting collisions with molecules in the atmosphere break up neutral atoms and molecules into electrons and ions. In high speed flows, shock heating and viscous dissipation in the boundary layers will first lead to a dissociation of the participating molecules (breaking the molecules into their individual atoms). At higher temperatures the collisions are so violent that the electrons are knocked clear of the nucleus. These free electrons and ions form plasma. Re-entry vehicles all generate plasmas due to shock and boundary layer heating, Radio waves cannot penetrate the highly conductive plasma, and therefore, re-entry vehicles suffer from a temporary radio blackout during re-entries.
Resulting heat from collisions with molecules in the atmosphere makes it imperative for a spacecraft's re-entry “angle of attack,” relative to the atmosphere, to fall within a certain range. The angle of attack must be oriented so that the heatshield absorbs the bulk of the re-entry heat. If the angle of attack is too shallow, the spacecraft will skip off the atmosphere and back towards space (similar to a stone being skipped across the surface of a lake. If the spacecraft's angle of attack is too steep, the spacecraft risks burning up due to extremely high heat loads from excess friction with air molecules. The window of a successful angle of attack depends on the spacecraft's geometry, speed, and surrounding air density. Air is less dense in the upper atmosphere, and thus, a spacecraft will encounter less friction. By approaching the Earth at a shallow angle, the spacecraft can spend more time in the upper atmosphere and increase the duration of deceleration. As the spacecraft moves lower into the atmosphere, it may have to adjust its angle of attack by means of a reaction control system, a mechanism for attitude control.
A reaction control system (RCS) is a subsystem of a spacecraft that is used for reentry flight dynamics. Its purpose is for attitude control and steering. Attitude control refers to control of the angular position or rotation of the spacecraft relative to the object it is orbiting. These angles are referred to as pitch, yaw and roll. An RCS system is capable of providing small amounts of thrust in a desired direction or combination of directions. An RCS system is also capable of providing torque to allow control of rotation. This is in contrast to a spacecraft's main engine, which is only capable of providing thrust in one direction, but is much more powerful. RCS systems can be used not only for attitude control during re-entry, but also for station keeping in orbit, close maneuvering during docking procedures, control of orientation, and as a backup means of de-orbiting.
On Apr. 12, 1961, the Soviet space program launched Vostok I, the first manned space flight. Although the service module of Vostok I had a nitrogen gas RCS system, the re-entry capsule lacked an RCS system and was unable to perform attitude control in flight. Hence, the Vostok I re-entry capsule was designed as a sphere with a heat-shield covering the entire outer surface of the capsule.
U.S. Pat. No. 3,093,346, by Faget, discloses a space re-entry capsule with attitude control via means of generating torque thrusts. The '346 patent was the design for the space capsules for Project Mercury, the first United States manned flight space program. The Mercury capsules were equipped with hydrogen peroxide (H202) RCS systems providing thrust for attitude control.
Over the years, RCS systems have improved, but generally maintain the same operating concepts. Conventional RCS systems are limited by the amount of fuel the spacecraft can carry. If the fuel is completely exhausted before reentry is achieved, the spacecraft will lose its ability for attitude control.
Beginning with the United States Gemini Program, the standard for fuel on many RCS systems has been a hi-propellant hypergolic liquid combination of an oxidizer and hydrazine or a hydrazine derivative. One example is monomethyl hydrazine (CH3N2H2) with nitrogen tetroxide (N202) as an oxidizer. As a hypergolic combination, the constituents ignite on contact with each other and create the thrust force for the RCS system. The disadvantages are that hydrazine, hydrazine derivatives, and nitrogen tetraoxide are generally expensive and toxic. Orion Propulsion, Inc. has developed an oxygen and methane RCS thruster for use on spacecraft. However, the inventor is not aware of an RCS system that uses oxygen or a mixture of oxygen and other gases (i.e. nitrogen) breathable by humans having the further feature of serving as part of a secondary life support system. Because spacecraft can only carry a finite amount of fuel with little chance to refill, some alternative RCS systems have been developed so that fuel can be conserved. One such alternative RCS system used momentum wheels which spin to control rotational rates on a vehicle. U.S. Pat. Nos. 6,834,561, 6,463,365, and 5,386,738 describe a control moment gyroscope (“CMG”) for spacecraft attitude control consisting of a spinning rotor and one or more motorized gimbals that tilt the rotor's angular momentum. As the rotor tilts, the changing angular momentum causes a gyroscopic torque that rotates the spacecraft.
A tether is a long cable usually made of thin strands of high-strength fibers and/or conducting wires. It is known in the art to use tethers to decelerate and deorbit an object in orbit around a celestial body. The inventor is not aware of any publication describing the use of a tether for spacecraft attitude control.
The NASA Hypersonic Inflatable Heat Shield Prototype which was tested suborbitally on Jul. 23, 2012 is an example of an inflatable reentry system. This inflatable aeroshell flew on a suborbital trajectory to an altitude of 450 km, was inflated from a 22 inch wide 308 kg capsule into a 10 foot wide aeroshell, and reentered successfully. This inflatable aeroshell, although having the great benefit of being lightweight, experienced 20 gravities of force and 1,000 degrees Fahrenheit of heat load during the suborbital reentry.
Another example of an inflatable aeroshell reentry system is the Russian Inflatable Re-entry and Descent Technology system that was launched to Low Earth Orbit in February 2000 from the Baikonur Cosmodrome. After completing six orbits of the Earth, the system reentered, but it also experienced extreme and uncontrollable thermal and force loads.
It should be noted that both of these inflatable reentry aeroshell vehicles incorporate robust thermal protection systems adapted to withstand ballistic reentry from hypersonic and orbital velocities. Neither of these aeroshell inflatable reentry systems has any means of controlling its trajectory, angle of attack or G force load.
The approaches described in this section are approaches that could be pursued, but not necessarily approaches that have been previously conceived or pursued. Therefore, unless otherwise indicated, it should not be assumed that any of the approaches described in this section qualify as prior art merely by virtue of their inclusion in this section.
The present invention is directed at a system, method, or apparatus satisfying the need to overcome the fuel limitations and costs of conventional RCS systems.
An embodiment allows spacecraft attitude control by generating moments about the center of mass, producing an angular acceleration. An embodiment does not require motorized gimbals and spinning rotors to generate a torque force. Rather, torque is generated by producing a friction force on a tether at a distance from an axis of rotation. In accordance with an embodiment, there is a length of tether extending from the forward section of a space capsule in a manner such that the tether's force line is aligned perpendicularly through the capsule's center of mass. “Forward section” as used herein, is in reference to the forward section of a spacecraft during the launching phase. The tether is held by a reel which can vary the length of the section of tether extended from the body of a spacecraft. The tether can be kept enclosed within the capsule during flight and deployed or extended using a reel when the capsule is preparing for re-entry. The reel can be operated by crew members or remotely from a mission control facility on Earth. An excessive length of tether is wound about the reel so that additional tether can be unreeled in the event that part or all of the extended portion of tether is severed or burned away from the spacecraft during re-entry. The tether, when extended from the forward section of the space capsule, functions as a hypersonic parachute, decreasing the capsule's velocity via drag forces.
In accordance with an embodiment, the tether can be made of a heat resistant material with conductive properties making it suitable to function as an antenna. The unexposed end of the tether can be integrated with the spacecraft's radio communications system. The entire length of tether can serve as an antenna for transmitting and receiving radio communications. The reel can feasibly extend the tether to a sufficient length for radio communications to avoid the conductive plasma generated from shock and boundary heating during re-entry, allowing continuous radio communications.
An embodiment, can comprise a “tether direction control apparatus” (hereinafter referred to as “control apparatus”) for offsetting the tether's force line away from a space capsule's center of mass for the purpose of attitude control. If there is no offset, the control apparatus remains in a “zero position”, or centrally aligned with the capsule's z-axis. During re-entry flight, due to drag forces, the extended tether will hold a position parallel to the direction of air resistance relative to the capsule. By offsetting the tether while a drag force is exerted upon it, a moment equal to the drag force multiplied by the distance of the offset from the control apparatus's zero position will be produced about the capsule's center of mass. Moments created by the offset and drag forces can be used to alter pitch, yaw and roll with regards to attitude control.
Certain embodiments can comprise a propellant-based RCS system to serve as a backup system for spacecraft attitude control in the event that a primary tether-based RCS system fails to operate successfully. The propellant-based RCS system functions on the same principle as traditional propellant-based RCS systems, in that it provides attitude control via a series of torque thrusts generated by exhaust of a gas propellant through translation thrusters. The placement of the translation thrusters generally requires that the thrust vector be aligned to pass through the z-axis of the spacecraft or an unwanted roll or rotation will result when the thrusters are fired.
An embodiment's propellant-based RCS system can comprise a ventilation line connected to a backup life support system. For certain embodiments, the propellant comprises compressed oxygen gas. The propellant can comprise a combination of compressed oxygen and other gases (i.e. nitrogen) in a ratio suitable for sustaining human life. In the event that oxygen levels in the passenger cabin fall below a safe level for humans, the backup life support system can be activated by controls within the spacecraft or automatically activated when one or more sensors detect oxygen and/or carbon dioxide levels.
The present invention is illustrated by way of examples, not by way of limitation, in the figures of the accompanying drawings and in which like reference numerals refer to similar elements and in which:
In the following description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the present invention. It will be apparent, however, that the present invention may be practiced without these specific details. In other instances, well-known structures and devices are shown in block diagram form to avoid unnecessarily obscuring the present invention.
A lid 101 covering the opening at the forward section of the capsule detaches from the capsule body 100 during re-entry and prior to the release of the tree tether 105 from the capsule body 100, the detachment may be triggered via controls within the capsule or remotely from a mission control facility on Earth. At this point, the lid 101 can be considered to be a discarded expendable component. The free tether 105 is released from the capsule 100. The tether 105 remains threaded through the center hole 106 of the control apparatus 104 and held in place by the reel 102. The control apparatus 104 is initially in a zero position.
Each end of the cylindrical shaft 109 is connected to the center of a circular flange 108. During flight, the tether 103 is stored by being wound around the cylindrical shaft 109 on the reel 102, between the flanges 108. One circular flange 108 is connected to a driveshaft 110, said driveshaft 110 being further connected to a motor 111. The motor 111 is bi-directional and turns the driveshaft 110 clockwise or counter-clockwise. When driven by the motor 111, the driveshaft 110 turns the reel 102 accordingly and retracts or extends the tether 103 depending on the direction of the motor. In this embodiment, turning the driveshaft 110 counter-clockwise extends the tether, or causes free tether 105 to lengthen, and turning the driveshaft 110 clockwise retracts the tether 103 shortening the free tether 105, The motor 111 is powered by a power source 112, and can be controlled via an interface 113 from the spacecraft's controls.
The free tether 105 remains threaded through the center hole 106 of a control apparatus 104, initially in a zero position. The free tether 105 remains loosely coiled above the control apparatus 104 and remains so until the spacecraft is ready for atmospheric re-entry. On re-entry, the free tether 105 is released into the atmosphere and drags behind the spacecraft, serving as a hypersonic parachute. The unexposed end 114 of the tether 103 is interfaced with the spacecraft's radio communications system 115, allowing the tether 103 to function as a radio antenna during atmospheric re-entry. Highly conductive plasma forms as a result of shock and boundary layer heating during atmospheric re-entry, blocking radio waves. The free tether 105 extends far behind the spacecraft, avoiding the conductive plasma and allowing the spacecraft to send and receive radio transmissions during re-entry. The tether 103 is made of a heat-resistant conductive material such as aluminum or steel for sending and receiving radio communications and withstanding the high temperature of the conductive plasma.
While the spacecraft has a primary life-support system, the backup RCS system further functions as part of a backup life support system, supplying the breathable gas propellant to the human occupants of the spacecraft should the oxygen or carbon dioxide levels reach a predetermined danger threshold and/or the primary life-support system fails. The gas propellant tank 122 is further connected to a backup life-support system 133 by a ventilation line 134. The backup life-support system comprises carbon dioxide canisters, fans, and filters. The carbon dioxide canisters remove carbon dioxide by reacting it with another chemical (i.e. lithium hydroxide, calcium hydroxide, sodium hydroxide), and the fans and filters remove dust and trace odors from within the spacecraft. The backup life-support system can be activated manually by controls within the spacecraft, or automatically when one or more sensors detect that oxygen or carbon dioxide concentrations have reached an unsafe level.
The backup RCS also functions as part of a backup life support system 133 comprising fans and filters for removing dust, odors, and carbon dioxide. The gas propellant tank 122 is connected to the backup life support system by a ventilation line 134 and can be activated manually by controls within the spacecraft, or automatically when one or more sensors detect that oxygen or carbon dioxide concentrations have reached an unsafe level.
In
The springs 121 attaching the positive-x and negative-x arms 116, 117 to the control apparatus 104 are also pulled in the negative-y direction by the ends of the springs 121 attached to said control-apparatus 104. The effect is that the center hole 106 is offset from the z-axis 107 in the negative-y direction by the measurement equal to the distance 135 from the center hole's 106 current position and its original zero position. A tether threaded through the center hole 106 would similarly be offset from the z-axis 107.
Relative to the x-y axes as depicted in
In
A moment is a rotational effect produced by a force at some distance from an axis of rotation. The moment (M) is equal to the product of the force (F) and the distance (d) from the axis of rotation about which it is applied. In
The line that passes through the capsule's center of gravity and is perpendicular to both the capsule's z-axis 136 and said moment arm 137 is the axis of rotation and in this case happens to be the y-axis. Moment 139 is the product of the friction force produced by tether 138 and the offset 137, a distance equal to the difference between the center hole's 106 current position and z-axis. Here, the moment 139 acts about the capsule's center of mass in the clockwise direction, as viewed in
A friction force 138 is produced by the entire length of tether 105 colliding with molecules in the atmosphere 140. Because the friction force 138 is produced at a distance 135 from the z-axis 107, a moment 139 is produced, causing a rotational effect about the capsule's 100 center of mass 141. The z axis 107 moves in a counterclockwise direction, further from the atmosphere 140 by using the center of mass 141 as a pivot point, increasing θA. The end result is that by offsetting the tether 105, the capsule 100 has adjusted pitch by increasing θA for atmospheric reentry using the rotational effect generated by a moment 139.
A friction force 138 is produced by the entire length of tether 105 colliding with molecules in the atmosphere 140. Because the friction force 138 is produced at a distance 135 from the z-axis 107, a moment 139 is produced, causing a rotational effect about the capsule's 100 center of mass 141. The z-axis 107 is shifted in a clockwise direction, relative to the capsule's 100 orientation in this illustration, using the center of mass 141 as a pivot point. The z-axis's new position 143 acts as the capsule's new approach vector 144. The end result is that by offsetting the tether 105, the capsule 100 has adjusted yaw and steered its direction to a new approach vector 144 which is to the right of the original approach vector 142, relative to the capsule's 100 orientation in this illustration.
A friction force 138 is produced by the entire length of tether 105 colliding with molecules in the atmosphere 140. Because the friction force 138 is produced at a distance 135 from the z-axis 107, a moment 139 is produced, causing a rotational effect about the capsule's center of mass 141. The z-axis 107 is shifted in a counter-clockwise direction relative to the capsule's orientation in this illustration, using the center of mass 141 as a pivot point. The z-axis's new position 143 acts as the capsule's new approach vector 144. The end result is that by offsetting the tether 105, the capsule 100 has adjusted yaw and steered its direction to a new approach vector 144 which is to the left of the original approach vector 142, relative to the capsule's orientation in this illustration. An example device comprises an aerodynamic decelerator for producing a tension force from aerodynamic drag, a vehicle capable of producing lift force that changes as the angle of attack varies, and a controller connected to the decelerator and also connected to the aerodynamic vehicle, for applying the tension force produced by the decelerator to the vehicle in a controlled manner to change the vehicle's angle of attack, thus varying the lift force produced by the vehicle. An example method is to adjust the load vector applied to a reentry vehicle by an aerodynamic decelerator to accomplish controlled skip entry through the upper atmosphere thereby reducing the vehicle's forward velocity vector and minimizing dynamic and thermal loads.
An illustrative embodiment comprises an aerodynamic decelerator for producing a tension force from aerodynamic drag, a vehicle capable of producing lift force that varies with angle of attack, and a control device connecting the aerodynamic decelerator to the aerodynamic vehicle, for applying the tension force produced by the decelerator to the vehicle in a controlled manner to change the vehicle's angle of attack, thus varying the lift force produced by the vehicle. An example method comprises using the devices described herein to accomplish controlled skip entry through the upper atmosphere in order to more gradually reduce the vehicle's forward velocity vector and thereby reduce dynamic and thermal loads.
A pull point is a point where drag generated by a tether is exerted on a space craft. For inflatable structures, rigidity is the degree of pressurization sufficient to retain structural integrity despite forces encountered during atmospheric entry. Skip entry is a method of atmospheric entry comprising one or more “skips” off of the atmosphere where energy is lost.
An example atmospheric entry device comprising a controllable tether and an aerodynamic body can use a skip entry technique.
Frictional drag force applied to a tether induces tension in the tether. An example device can comprise a tether attached to an aerodynamic body and a tether controller that shifts the pull point's location.
Atmospheric re-entry, even when initiated from a circular low-Earth orbit, typically requires a thermal protection system comprising a heat shield, ablative material, or radiative dissipation techniques. Semi-analytical and numerical simulations of the atmospheric re-entry from low-Earth orbits of a capsule with a 20-km heat resistant tether attached have shown that the thermal input flux on the capsule is reduced by more than one order of magnitude with respect to a comparable re-entry without a tether.
Long tethers have low ballistic coefficients and a large surface for heat dissipation. Moreover, a long tether is stabilized by gravity gradient and consequently tends to maintain a high angle of attack with respect to the wind velocity. The exposed surface of a 20 kilometer long 1 millimeter diameter tether is 20 square meters, which is larger than the typical cross section of a re-entry capsule. For example, the Apollo command module's cross section is under 12 square meters. The resulting strong drag decelerates the capsule during re-entry. Where an example embodiment allows variance of the application of drag force to the re-entry vehicle, so that the force vector is offset from the vehicle's center of mass, to change the vehicle's angle of attack, allowing control of the vehicle's reentry flight path. By using this method to allow the reentry vehicle to skip in and out of the atmosphere, especially during the portion of the reentry process where the greatest thermal and mechanical loads are produced, extends the time of reentry flight in a controlled manner and reduces the peak mechanical and thermal forces acting on the reentry vehicle.
A device, and the method of its use, for controlled atmospheric entry allows for an atmospheric entry system (“system”) that may be controlled in order to produce a flight path having reduced deceleration loads over an extended period of time in order to decrease mechanical and thermal loads and stresses that occur during atmospheric entry. This will allow for a gentler more precisely controlled transit through an atmosphere.
An embodiment can be used with any spacecraft design that produces lifting force that can be varied by varying the vehicle's angle of attack. A winged shuttle, a lifting body, and a space capsule are all examples of such spacecraft. An embodiment can work with an inflatable reentry system that uses a pressurized flexible toroid to provide rigid support to a flexible conical aeroshell payload section. Synergistic benefits may be obtained by using a lightweight inflatable reentry vehicle in conjunction with an embodiment.
Aerodynamic Vehicle (100):
An inflatable aeroshell is an example of the type of aerodynamic vehicle discussed in the embodiment below. Examples of inflatable aeroshells comprise the NASA inflatable aeroshell and the Russian inflatable aeroshell. These inflatable aeroshells have the benefit of being lightweight. They also have the disadvantage of only operating in an uncontrolled ballistic reentry trajectory, which produces large thermal and mechanical loads.
An aerodynamic reentry vehicle capable of producing variable lift force when its angle of attack is varied may also be used. Examples of such vehicles include the USAF X37B robotic reusable reentry spacecraft as shown in
The inflatable aeroshell skin 209 is fixedly attached to the exterior surface of toroid 207 at one or more locations about the circumference of the exterior surface of the toroid (possibly continuously attached such that there is no fluid communication between the inflatable toroid and the skin), and extends below the center of the toroid, forming an inverted payload volume 225. The skin being positioned in such a way that it forms a conical shape extending below the toroid. Skin 209 may be made of any material having mechanical strength sufficient to support the payload during reentry. It may be adapted to protect payload 225 against heating and may optionally comprise a thermal protection material 211.
The aerodynamic body 100 comprises the toroid 207 and the skin 209 that extends from the toroid, forming a conical nose 211 that sits underneath the toroid. This aerodynamic body automatically orients the entry system due to the action of natural aerodynamic forces in such a way as to have the aerodynamic body sit between the payload and the source of gravity with the nose-cone 211 pointing toward the gravitational source. Although an embodiment can reduce the thermal and mechanical loads acting on the reentry system and payload, it does not eliminate these loads entirely. The aerodynamic body therefore also shields the payload of the system from the majority of the frictional and thermal forces generated by the device's controlled transit through the atmosphere. The inflatable aerodynamic body acts as a lift generating body, aerodynamically orients and stabilizes the system, protects the rest of the system from friction with the atmosphere, can act as a backup life support system if it is pressurized with breathable gas, and may be provided with air bags to cushion the impact of the system when it impacts the surface of the celestial body.
The tether 105 is attached to the vehicle 100 by control system 203. The tether can be deployed by any number of different tether deployment means. The tether 105 is attached by any suitable mechanical means to the tether controller 203. The tether extends from the tether controller 203 to the tether end point 223. The tether 105 is deployed by a tether deployment means 158 and extends out behind the aerodynamic body 100. When the tether is extended it begins colliding with the molecules that constitute the atmosphere, causing friction. The friction generated by the tether exerts a drag force which is communicated to the rest of the atmospheric entry device at the Pull Point. If the tether is 20-km-long and 1-mm diameter tether, then its surface area is 20 square meters, which is much larger than the cross section of a re-entry capsule. For example the aerodynamic surface area of the 10 foot diameter NASA inflatable reentry system is about seven square meters.
The tether controller 203 is attached to tether 105 and through controllable lines 213, 215 and 217 to the exterior surface of the toroid. In
The tether control lines 213, 215 and 217 and the tether 105 itself can be extended or retracted by electric motors, by hydraulic or pneumatic actuators or even by manually pulling or winding the lines by mechanical means. This opens up the rather interesting possibility of extreme sports enthusiasts surfing the upper atmosphere from orbit.
In the lower atmosphere and at low, subsonic, speeds tether 105 may not generate a great deal of drag. Certain embodiments could incorporate a parachute into tether controller 203. This parachute could be deployed prior to landing. Alternatively, a deployable lighter than air balloon could be incorporated into the system to inflate to let the payload float in the atmosphere. Finally, the payload, with a balloon or parachute system, could be ejected after reentry, but prior to landing. An embodiment could also incorporate inflatable air bags or foam cushions to reduce the effect of the landing impact on the payload.
Inflation Means:
In one or more embodiments, the aerodynamic body comprises toroid, a skin, a payload volume, and one or more tether attachment points (possibly comprising a controller). There could be a way of going from a compacted deflated state to a rigid inflated state). The inflation of the toroid should be completed by the transfer of some material, gas or combination of gasses by an inflation means from a volume outside of the system to the un-inflated toroid of the aerodynamic body (preferably delivered from a pressure vessel). These details are not shown as inflating a flexible toroid is within the skill of those expert in the art. If an embodiment is used for human reentry, all or part of the toroid could be pressurized with oxygen or a breathing gas mixture to serve as a source of life support breathing gas.
Skip Entry:
Skip entry is a technique for entering an atmosphere. It is beneficial for entry systems that have a relatively low lift-to-drag ratio since these sorts of entry systems have difficulty extending their landing range and deceleration period due to their aerodynamic flight characteristics. When engaging in skip entry a space object makes one or more successive “skips” off of (or through) the atmosphere. Each successive “skip” reduces the energy of the space object relative to the celestial body whose atmosphere is being entered. The skip entry provides a space object entering an atmosphere a longer period of time and course of transit through the atmosphere. The increased period of transit increases the duration of time during which the entering object can shed energy relative to the celestial body. By increasing the duration of the atmospheric transit the energy of the space object can be released more gradually. This gradual reduction of the space object's energy is advantageous because it reduces both heating and rapid deceleration due to frictional forces that result from the space object's physical interaction with the molecules of gas and other particulates that comprise the atmosphere. Methods of achieving skip entry require precise guidance and control of the atmospheric entry system. Without precise guidance and control the atmospheric entry system attempting to achieve skip entry could fail to sustain its intended trajectory, which could result in one of a number of problems. If the atmospheric entry system takes too shallow of an entry trajectory, or achieves too much lift upon entry, the atmospheric entry system could skip entirely out of the atmosphere, and possibly out of the celestial body's gravity well. This could result in the complete loss of the atmospheric entry system and its payload. If the atmospheric entry system takes too steep of a trajectory, or has too small a velocity, the aerodynamics of the system may not generate enough lift for the atmospheric entry system to perform skip entry. This could cause the atmospheric entry system to engage in ballistic entry which could potentially destroy the atmospheric entry system and its payload due to excessive heating, high acceleration loads, or a high velocity impact with the surface of the celestial body. A third problematic possibility is that the atmospheric entry system achieves skip entry, but does so in such a way as to have the system and payload move off of its intended trajectory and land in an unintended location. The increased transit duration is effectively an increased flight path which allows for the atmosphere entry system to select a landing location from a larger potential landing area.
The method of skip entry into an atmosphere is achieved by calculating an appropriate trajectory, then initially descending into the outermost region of the atmosphere. After the initial descent, the aerodynamic profile of the atmospheric entry system (with or without help from some other forces including thrust or drag) generates lift which causes the entry system to ascend. As the object gradually ascends the gravitational force overrides the lift force and the object begins another descent into the atmosphere. This process may be repeated more than once before the atmospheric entry system loses the velocity (or other flight characteristics) required to generate sufficient lift to make another “skip.” When the object cannot, or does not wish to, make another “skip,” the atmospheric entry system travels along a ballistic trajectory through the remainder of the atmosphere. These “skips” increase the duration of the atmospheric entry system's transit in the upper, less dense, atmosphere. The increased duration of the atmospheric transit and the lower instantaneous deceleration gives this method of atmospheric entry many advantages as compared to fully ballistic atmospheric entry.
Increased duration of flight in the upper atmosphere is desirable because this is where most of the energy of reentry is dissipated. If the total energy release is made over more time, the effect is a much gentler and less stressful reentry.
The deceleration can be made more gradually resulting in lower acceleration loads being put on the payload and atmospheric entry system as a whole. Slower deceleration results in less intense aeroheating of the atmospheric entry system. The lower velocities and increased transit duration also reduce heat buildup on the entry system and acceleration loading on the system, which in turn allows for less mass of the system being dedicated to shielding. The increased atmospheric transit time coupled with the smaller velocities that the entry system achieves during the atmospheric entry allows the entry system more time to and ease of maneuvering.
An example method could facilitate atmospheric entry into the Earth's atmosphere from low Earth orbit. The devices and methods described can be used to provide a controlled reentry into any planetary atmosphere from any trajectory. Of course the control rules for each entry would be unique and would have to be calculated according to means well known to those skilled in the art of trajectory planning and atmospheric reentry.
The International Space Station (“ISS”) orbits the earth in low Earth orbit (“LEO”) 617. The outer bound of LEO is approximately 2,000 kilometers above the surface of the earth. In this orbit the ISS is traveling at approximately 18,000 miles per hour. In the event that some cargo needs to be safely, gently, and precisely transported from the ISS to a location on the surface of the Earth an embodiment could be used to accomplish the task.
Referring now to
As the system continues on a decaying orbit it begins to interact with the top layer of the atmosphere at an altitude of about 100 km 623. As the system begins to descend into the outer limits of the Earth's atmosphere the aerodynamic drag characteristics of the aerodynamic body will cause the system to orient so that said aerodynamic body shields the payload from aeroheating. At the same time the molecules of the outer atmosphere colliding with the tether generates a frictional drag force that induces a tension in the tether. The force of this tension is exerted in the opposite direction from the system's path through atmosphere. This drag force causes the system to further decelerate. If no control is exerted on the reentry system, it will follow a pure ballistic trajectory 625 and experience large thermal and structural loads before impacting on the surface at point 627.
As the tether decelerates the system the tether controller exerts forces onto the tether in such a way as to cause the tether drag force to impart a torque on the system. This torque force alters the system's aerodynamic flight characteristics, including angle of attack, generating a greater lift/drag ratio. As the lift force generated increases it will eventually offset the gravitational force between the system and the planetary surface at point 605, causing the system to gain altitude for a period of time and extend its flight distance and time through the atmosphere before the lift force is no longer greater than the gravitational force at which point the system begins to descend again at point 607. Skip entry may be used as many times as required to minimize the deceleration loads/rates on the payload, or as many times as are desired for the flight profile as is illustrated by points 609, 611, 613. With each skip the system loses energy. The control can be as simple as measuring the deceleration force and changing the angle of attack of the reentry vehicle when the deceleration exceeds some predefined limit. In theory this could yield the result of a gentle return to the surface.
The increased period of time that skip entry allows for the system to be traversing through the atmosphere allows for an increased period of time in which the tether's drag force can gently decelerate the system. Once the system's velocity, altitude, or other flight characteristic preclude the system from engaging in any further “skips” off of the atmosphere as shown at point 615 the system would begin a ballistic trajectory, although the tether and tether controller could still be used for some course modification. At this point in the system's descent should be sufficiently slow as to allow parachutes, or some other deceleration means to gently lower the system to a point on the surface of the Earth, or to simply let the system fall to the surface of the celestial body at point 617.
The flight path and controls exerted on the system can be optimized to provide for the slowest deceleration possible or controlled, constant, sustained deceleration rates. This allows for the atmospheric entries to be completed entirely within predetermined parameters (time of descent, acceleration load requirements, landing location).
In an embodiment the tether controller comprises one or more control cords attached to a point on the tether and a point on the toroid. The means of control would be by varying the tension on the one or more control cords so that the tether's drag force is deflected through the one or more control cords and exerted on the control cord's point(s) of attachment to the toroid. This action will cause the sum of the tensions on the tether and the one or more control cords to be functionally exerted on a point different than the device's center of gravity. This deflection would cause there to be a torque force impart to the device which would alter its aerodynamic flight characteristics. By changing the aerodynamic body's aerodynamic flight characteristics, including but not limited to its angle of attack, the lift to drag ratio of the system can be modified to correspond to a desired flight plan. In certain embodiments, the tether controller comprises a control device that exerts a force on the tether such that there is a controllable change in the angle between the tether and the aerodynamic body. When the force exerted by the tether does not pass through the spacecraft's center of mass a torque alters' the spacecraft's angle of attack to change the lift forces acting on the system.
In an example embodiment the tether controller comprises a means by which the pull point of the tether may be altered in at least a two dimensions to shift the tether's force vector away from the vehicle's center of gravity to impart a moment on the vehicle to adjust pitch (angle of attack) or yaw. Furthermore, the tether's length may be extended or retracted through the use of a winch to vary the magnitude of the drag force imparted by the tether. When combined this method of controlling the length of the extended tether and the location of the pull point would allow for control of the amount of drag force that the tether would generate and where that drag force is functionally imputed on the reentry vehicle. This will provide a reactive guidance means for the atmospheric entry system as it descends through an atmosphere. If the tether is electrodynamic it may also be possible to generate drag forces outside of the atmosphere.
In certain embodiments, the tether controller comprises a means of moving the tether's physical attachment point to the rest of the system (most likely to the toroid). By moving the tether's attachment point the tether controller also moves the point at which the tether's drag force is functionally imparted to the rest of the system. If the vector along which the tether's drag force is functionally imparted to the rest of the system is moved away from the system's center of mass the drag force will impart a moment on the aerodynamic body. This moment would generally result in a change in the system's angle of attack, resulting in a change in the system's aerodynamic flight characteristics.
Inflation:
In an example embodiment, the gas that is used to inflate the toroid to rigidity is Oxygen (O2). In this embodiment, there would be a hose equipped with a regulator allowing controllable fluid communication between the interior of the toroid and the payload volume. In this embodiment, the O2 used to pressurize the toroid to rigidity could be tapped into as a life support system. In another example embodiment, the gas that is used to inflate the toroid to rigidity is Helium (He). Helium, being an inert gas, is unlikely to react with any other chemicals that it may be exposed to, and thus is unlikely to be dangerous. Helium gas is light weight which could increase the spacecraft's buoyancy. The increased buoyancy of the system could result in a smoother, gentler deceleration and stop than a less buoyant system would allow.
In an example embodiment, the toroid has a plurality of discrete compartments in fluid communication with one another through valves. The valves could be configured so that discrete compartments could be inflated by the release of pressurization material into discrete compartments, but would not allow the decompression of one compartment to deflate any other of the discrete compartments. This configuration would help the toroid retain some structural rigidity in the event of a damage that would cause depressurization. In certain embodiments, the skin of the aerodynamic body comprises one or more compartments in fluid communication with the toroid. These compartments inflate with the same pressurization fluid as the toroid. This would aid in the prevention of communication of heat from the outer surface of the skin of the aerodynamic body to the inner surface of the skin of the aerodynamic body. Inflatable components are light weight and can be deflated and stored in a relatively small volume.
Certain embodiments comprise an aerodynamic body, a hypersonic decelerator, and a controller that is connected to both the hypersonic decelerator and the aerodynamic body such that the forces acting on the hypersonic decelerator are transmitted to the aerodynamic body in a controllable manner so the hypersonic decelerator and control means can be used as a RCS (reaction control system) for the aerodynamic body.
In this embodiment the aerodynamic body may comprise a capsule, a shuttle, a heat shield, or an inflatable aeroshell. The aerodynamic body will preferably have flight characteristic that can be influenced by the controlled use of the hypersonic decelerator in such a way as to result in an alteration of the aerodynamic body's flight characteristics (angle of attack, lift/drag ratio, aerodynamic profile, shape, etc.).
The aerodynamic body allows the forces generated by the hypersonic decelerator to be translated from pure decelerating drag into a means for altering the flight characteristics of the system in order for the system as a whole to be controllably maneuverable.
The hypersonic decelerator refers not to a specific decelerator system but instead could comprise one, more than one, or a combination of hypersonic decelerators that rely on drag to decelerate an object traveling through an atmosphere at hypersonic speeds. This group comprises: single line tethers, multi-line tethers, tapes, ribbons, ballutes, parachutes, wings, and sails. Different hypersonic decelerators could be used simultaneously or in sequence. The hypersonic decelerator uses the drag generated by friction between the hypersonic decelerator and the atmosphere to both decelerate the entire system and to generate forces that can be controllably transferred to the aerodynamic body in order to alter the spacecraft's flight dynamics in a controllable manner.
Multiple hypersonic decelerators may be used in conjunction with one another, either having multiple distinct hypersonic decelerators used at the same time, or with different hypersonic decelerators used at different periods during the system flight in order to correspond to different deceleration requirements. This may include but is not limited to the use of different hypersonic decelerators depending on the velocity range that the system is traveling within, or the use of different hypersonic decelerators depending on the density of the atmosphere through which the system is traveling.
The control means can be any device, system, or method of use that allows for the drag forces generated by the hypersonic decelerator to be used to alter any one or more flight characteristics of the system in such a way as to allow the flight path of the system to be controlled by alteration of the drag forces acting on the hypersonic decelerator.
In a illustrative embodiment, the aerodynamic body comprises an inflatable aeroshell, the hypersonic decelerator comprises a tether, and the controller allows the spacecraft to engage in skip (single or multiple) entry of an atmosphere.
Mathematical Model of the Invention
The present invention has not yet been flight tested, so it is not possible to provide experimental examples to inform those skilled in the art. In order to provide information about the present invention to and those skilled in the art, the inventor has provided below a detailed mathematical model of the present invention. The general conclusion from this mathematical model is that de-orbit forces and heating rates are much less for the controlled reentry system made possible by the present invention.
Introduction
Of greater importance to the use of the tether as an entry trajectory controller is the observation that in the past attempts to use tethers on reentering objects have been unsuccessful because the tension distribution along the length of the tether becomes chaotic. Many slack regions develop and the tension at any single point is only intermittently positive. Moreover, the tension at the attachment point, which actually provides the control force, varies widely and is often zero. Trajectory control is impractical or impossible in such a system unless the system can be made to produce a non-zero, smoothly varying and predictable tension at the attachment point. One way to ameliorate this problem is to choose materials that offer diffuse reflection of atmospheric molecules in free molecular flow. In this case, the magnitude of the aerodynamic force is still proportional to the normal component of stream velocity, but it is also proportional to stream velocity itself, and thus can act in both the normal and tangential directions. As will be discussed in more detail below, this causes the tether to be smoother and the developments occur slower, with fewer slack regions. Thus the use of non-metallic tether material that diffusely reflects atmospheric molecules in free molecular flow is recommended for the present invention. In the following analysis only diffuse reflection is assumed. In practice this means the use of non-metallic tether material that diffusely reflects atmospheric molecules in free molecular flow is recommended.
As will be discussed in detail below, the tether tension at the attachment point to the spacecraft is initially small and it is sensitive to initial conditions. It is helpful to add an aft-body attached to the free end of the tether. The addition of an aft-body as part of the present invention is useful in order to reduce sensitivity to initial conditions and to generate non-negligible tether tension in the early part of the entry trajectory.
It may take as long as a half-hour for the maximum departure from the direction of spacecraft flight, i.e. the x-axis, to fall below 1 km in altitude, i.e. the y-axis, even with a large aft-body. Because of this it is desirable to gradually deploy the tether behind the spacecraft at a rate that maintains positive tension in the tether. When the towing speed is nearly the orbital speed, the deployed tether will be aligned with the x-axis. This method of operation is recommended. At the same time, if the gradual deployment method should fail and the tether deploys all at once, the action of the aft-body will eventually align it with the stream direction anyway.
In summary:
a. To promote faster stabilization in the initial portion of atmospheric entry, the tether material should be capable of diffuse reflection of atmospheric molecules.
b. In order to prevent large uncertainties in initial entry, the tether should be deployed gradually a so that the initial shape is parallel to the stream direction.
c. To provide insensitivity to initial conditions, guarantee some initial tension and devise some mechanism whereby the tether settles into to a steady shape and tension distribution for all flight regimes, it is necessary to add an aft-body attached to the free end of the tether.
d. The aft-body design parameters must be chosen to prevent flag flutter instability.
In the stability analysis to follow, it is assumed that the atmospheric density is constant along the tether and equal to its value at the attachment point, despite the very large tether length of 20 km in the nominal design concept. This is justified on a number of counts. First, for the flight regimes, (H,V) wherein the attachment point tension is largest, the tether is swept back toward the x axis, so that the variation of altitude along the tether is a small fraction of the nominal length. Secondly, models of heating and breakage conditions described in detail below show that a 20 km tether will experience extreme heating rate gradients and frequent breakage near the free end—thus necessitating the use of a shorter tether. A shorter tether will also be less sensitive to varying atmospheric conditions with altitude since it will “sample” variations over a smaller altitude range. Finally, a shorter tether, by reducing the spatial dependence of the density will render the dynamics more predictable.
The coefficients of the lateral modes behave as uncoupled, damped oscillators and are obviously exponentially stable. For sufficiently long time, the coupling becomes negligible so that the longitudinal mode coefficient also become uncoupled damped oscillators. Thus local stability can be concluded for the tether with aft-body as a matter of general principle.
This conclusion can be verified by assuming a 0.5 m radius aft-body and spot checking the plane instability regions determined previously with no aft-body. Starting each case with a several kilometer deviation from the equilibrium shape, it is found that all cases are stable. An example is provided by the case H=100 km, V/V0=0.85 where the tether was started in a severe 4 km oscillation of the first lateral mode, resulting in significant variation in tension. After 730 seconds, the oscillation amplitudes in both the shape and tension distribution are significantly reduced. Finally, after approximately 1500 seconds, the shape and tension settle into equilibrium. The foregoing results show that a modestly sized, 0.5 meter radius aft-body, and assuming initial departures of 3-4 km or less in the initial deviation from the equilibrium, the tether/aft-body system will converge to a stabile equilibrium. This does not mean that the equilibrium state will be stable for larger initial deviations. However, the results discussed in detail in the mathematical model below do show that if initial deployment is slow and smooth and if altitude does not vary very rapidly, the system can stay close to a stable equilibrium shape and tension during the course of entry. Because the equilibrium state is relatively simple and predictable, this is a great advantage for the development of a robust, dependable entry trajectory controller.
In Summary:
1. A small aft-body avoids chaotic motion initial deviations from the equilibrium state
2. For flight regimes in which the variations of atmospheric properties over the tether length are small, the equilibrium shape is simple and relatively insensitive to temporal and spatial variations in atmospheric density. To secure this advantage over all altitude and velocity ranges, the tether length should be just long enough to provide the control force required. Moreover, a long tether will experience extreme heating rate gradients and frequent breakage near the free end. In other words, if we do not shorten the tether, Mother Nature will do it for us! Since the key parameter determining tension is the tether length x radius (for a round tether) it is possible to reduce the length while increasing the radius. The tether cross section need not be round, but may be a rectangular tape or an oval. A circular cross section tether is assumed in the mathematical model of the present invention set forth below for ease of calculation of this embodiment.
The entry profile for the present invention should be planned so density variations in the atmosphere are slow enough to allow the tether to follow its equilibrium configuration fairly closely. The mathematical model of the present invention shown in detail below teaches that the tether is efficient in redirecting the angle-of-attack of the spacecraft to match the tether tangent direction near the attachment point. If the tether is maintained in a sagged equilibrium shape, without oscillation, then the main effect is to impose a positive angle-of-attack on the fore-body. Moreover the moment arm of the tether tension about the spacecraft center-of-mass will be much larger than can be attained by only adjusting the attachment location of the tether to the back of the spacecraft. Thus changing the tether length allows more control force than adjusting the attachment point.
In Summary:
1. The present invention's entry profile must be planned so atmospheric density variations are slow enough to allow the tether to follow its equilibrium configuration fairly closely.
2. The tether appears to be quite efficient in redirecting the angle-of-attack of the fore-body to match the tether tangent direction near the attachment point. It is probable that changing the tether length might be a more effective control strategy than adjusting the attachment point. So how can the present invention be used to aid atmospheric entry? The overall goal of the present invention is to prolong the atmospheric entry process so heating rates on the fore-body are significantly reduced. To control entry the present invention can modulate the tether length, which can be increased or decreased and can vary the position of the attachment point of the tether relative to the spacecraft's center of mass. If the present invention includes sensors that can measure altitude and both the horizontal and vertical components of velocity and measure the tether's sag angle, then one possibility for using the tether to reduce entry speed is to periodically deploy and retract the tether so as to stimulate skipping in order to periodically “cool off” the vehicle by repeatedly moving to higher altitude. The model shown in detail below compares two cases, uncontrolled and controlled by the present invention, that feature entry from low Earth orbit. For both cases:
The tether is 500 m long and 2 mm in diameter,
In fore-body model, we set A=10 so that the lift-to-drag ratio is approximately unity
The system starts at orbital speed at 1701 cm altitude.
The high altitude behavior is essentially the same in both cases because of the very low atmospheric density and the small frontal area of the tether. The significant differences between the two cases occur once they reach approximately 110 km altitude. Once the sensible atmosphere is reached, the pitch angle aligns with the velocity direction.
The uncontrolled reentry system sustains a shallow glide path and then very abruptly dives starting at about 95 km. The velocity versus time plot shows that the vehicle sustains nearly orbital velocity until it very suddenly comes to rest. On the whole, the entire entry event occupies only about a dozen minutes. This causes the “reentry fireball” whose thermal and mechanical loads are a major challenge to reusable spacecraft design designers.
In contrast, the controlled system develops larger amplitude skips before reaching 110 km. Well before this the controller makes the angle-of-attack converge to 45 degrees, thereby producing maximum lift. At about 100 km altitude, near the “Kármán line”, the controlled system skips decrease in amplitude and increase in frequency. For a period of nearly 25 minutes, the controlled system maintains nearly horizontal flight (taking into account the skipping motion) that is centered on the Kármán line. It is during this period, the sequence of skips produces nearly step-changes in velocity. Once the skip amplitude dies out, the system completes entry to the surface with relatively gentle declines in altitude and velocity. Overall, the final entry, starting at 100 km, takes roughly 20 minutes nearly twice as long as in the uncontrolled case. The consequences for entry heating of the spacecraft depend on the details of the entry body. However, the skip amplitudes decrease and then cease after about 15 minutes. The maximum heating values for the controlled system are about two orders-of-magnitude below the uncontrolled case. Although the significant heating of the uncontrolled system lasts longer for the controlled system, the times spent near maximum heating values is only about three times that of the uncontrolled system. Thus the integrated heating rate, which is a measure of maximum temperature, is about an order-of-magnitude less for the controlled system than for the uncontrolled system.
The g-loads on the forebody structure, i.e. the inertial acceleration minus gravity are very different for the uncontrolled and controlled systems. The uncontrolled system shows a single large peak at about 5.5 gs. The controlled system experiences a lengthy sequence of maxima, corresponding to the skips, at values between 0.6 to 1.1 g, with very brief pulse at 1.9 g. The present invention can also be used by spacecraft for reentry from lunar or deep space trajectories, high elliptical orbits and generally any trajectory or orbit that is more energetic than low Earth orbit. First the tether is used as an aerobrake at the orbits perigee, which lowers the apogee. Then, when the orbit is circular low Earth orbit, the tether is used to control the angle of attack of the spacecraft as it reenters the atmosphere. The ‘deep space’ tether used for aerobraking the spacecraft into low Earth orbit could be a long tether and the ‘low Earth orbit’ tether used for controlling the spacecraft angle of attack could be a short tether. They can be two separate systems or they may be one system of controllable lengths.
If we inspect the velocity time histories for the uncontrolled deep space system there is a single, abrupt decrease in speed down to 2 km/s. The undulations seen in the later half of the entry process are due to uncontrolled oscillation of the angle-of-attack. The velocity of the controlled system drops suddenly on penetrating 100 km, but only drops to about 5.5 km/s. There follows a lengthy plateau corresponding to the first skip, then the second skip reduces the velocity to 2 km/s. Subsequently the vehicle glides smoothly to the ground with very small final velocity. Obviously, the controlled entry process is much more gentle and gradual than for the uncontrolled system.
A similar pattern is consistently repeated for the heating rate estimates. The uncontrolled system has a lone spike at 2.6×109, whereas the controlled system shows a very brief spike with a maximum at 0.8×109. Then there is a very small spike due to the second skip. Overall, the integrated heating rate of the uncontrolled system is about 2.5 times that of the controlled system. This is largely due to the fact that both systems must inevitably hit the Karman line as the first obstacle in their descent. This reduction in heating can be improved by using a tether aerobrake to make the high energy orbit into a low Earth orbit.
Finally, the uncontrolled system exhibits a maximum of over 300 gs. This is the deceleration profile for a meteorite. The vehicle is a fireball until it is nearly in free fall until it hits the ground at high velocity. The controlled system experiences a maximum of 30 g s, and this can be reduced to a short 1.8 g load by lowering the orbit from elliptical to low earth circular by using the present invention first as an aerobreak at orbital perigee and then, when the orbit is circular, using the present invention to control the skipping of the payload to lose energy near the 100 km line.
Prior to experimental testing of the present invention the inventors plan to include breaking conditions in the simulation code; examine thermal performance of the tether and investigate combined control of tether length as well as attachment point location.
Skipper Atmospheric Entry System Dynamics
1. Introduction
The present invention, which the inventors have named the “Skipper” is an atmospheric entry system comprising a fore-body, which may be a spacecraft, an ore body produced by space mining, or any other physical object to be reentered though an atmosphere, trailing a long, thin aerodynamic decelerator tether, which may be of any length and shape. The fore-body may be shaped to provide various desirable aerodynamic properties, e.g. a particular lift-to-drag ratio, etc. The tether is placed in tension due to aerodynamic forces applied during entry. By varying the point of attachment on the fore-body and/or the length of the tether, the tension force on the fore-body can be changed and made to alter the angle of attack of the fore-body. This change of angle of attack allows the entry trajectory of the fore-body and tether to be controlled so mechanical stress on and aero-thermal heating of the fore-body significantly reduced. Flight tests have not yet been conducted on this invention. The inventors therefore seek to substantiate that the postulated benefits of the present invention can be achieved by conducting a comprehensive mathematical modeling analysis of the dynamics, mechanics and heat transfer of the invention, the skipper system. As the invention is further developed, these mathematical models will be verified or refuted and refined by real world experience; experimental results from flight tests of the present invention.
This analysis starts by investigating the dynamics of the tether during high altitude flight through the atmosphere. This is done to examine two issues. First, past towed tether systems have exhibited a combination of a Bernard-von Karman vortex street and Kelvin-Helmholtz instability. In these cases the tether exhibited an unstable fluttering motion rather than assuming a simple, steady shape behind the towing body. The resulting instability built up large, erratic motions that rendered the past towed tether system useless. This past experience involved low altitude, subsonic flight at continuum flow conditions, and thus was significantly different from the system taught by the present invention. It is important to show that the present invention either will not exhibit such instabilities or can be modified to avoid them. The second issue explored is whether or not the aerodynamics of the tether will produce a steady, predictable tension force at the attachment point to the fore-body that is insensitive to reasonably anticipated initial conditions. This steady, predictable tension force is essential to successful design of the trajectory controller since the attachment point tension is the control actuation force. This condition is uncertain, may or may not hold, even when there are no large scale instabilities. To examine the relatively fast time-scale dynamics of the towed tether, we begin by making some simplifying assumptions about the fore-body motion. We assume the fore-body is much more massive than the tether and sustains horizontal flight at a constant speed that is equal to or less than the circular orbit speed at the given altitude. The attachment point is assumed to be fixed at the system center-of-mass. With these assumptions, we derive a system of equations for the tether motion, starting from various initial conditions. A simulation code based on the equations of motion is used to examine the tether motion. Later in this mathematical model, we will remove the above assumptions and generalize the formulation to include coupled fore-body/tether dynamics, aerodynamic heating and tether breaking conditions, leading ultimately to analysis and simulation of the trajectory control system.
2. Coordinate System and Reference Frame Definition
We refer the system motion to a rotating frame, the origin of which is centered at the system center-of-mass, which is identical to the fore-body center-of-mass. For this initial investigation, we assume that the fore-body exerts whatever lift forces are required to maintain horizontal flight at constant speed, V, over a non-rotating Earth. The rotating frame and coordinate axes are illustrated in
The reference frame rotates relative to inertial space at angular rate, Ω{circumflex over (z)}, where:
Ω=V/(RE+H)
RE=Earth radius
H=Altitude above the surface (2-1.a-c)
Note that if V is less that the orbital speed for a circular orbit at altitude H, the system is subject to a gravitational acceleration, {tilde over (g)}, in the negative y-direction, given by:
At this point, we need to relate the time derivatives of various quantities relative to the inertial frame to their apparent rates of change as seen in the rotating frame. Denote by
the inertial derivative of a quantity, and by a dot over the quantity for the time derivative relative to the rotating frame. If
Likewise, the inertial acceleration of the point is given by:
Of principal importance is the difference between the centripetal and Coriolis accelerations and the gravitational acceleration, which we denote by
Or using (2-2.b) to set μE=(RE+H)3Ω02, we have:
Thus, we obtain:
Note that if
Since the spatial extent of the system is likely to be much smaller than (RE+H), we can approximate the expression (2.7) by expanding it out to linear terms in x, y, or z and neglecting higher order terms:
If
3. Lumped Mass Model of the Tether
The initial skipper system evaluation by this model presumes the use of a long and very thin (of the order of millimeters) filament. The length-to-radius ratio is such that the tether is unlikely to have significant bending stiffness. Thus the appropriate model is a string; understood to be a filament that can support only axial tension. Moreover, it is also reasonable to assume that the filament offers no resistance to compression.
The inventor's experience with the dynamics of cable and tether systems indicates that the simplest and most effective model is a lumped mass model. This is illustrated in
As illustrated in
1=L/N (3-1)
Assume that the tether is of circular cross-section and allow for the possibility that the radius, r(s), can vary with the arc-lengths, where s=0 at the attachment point. The N segments are labeled with indices starting at unity near the attachment point, and likewise for the N+1 point masses. The position vectors of the point masses in the rotating frame are denoted by
ρr=Volumetric mass density of the tether material
For later use in the aerodynamic model, we also define the area, Ak, projected by each segment onto a plane perpendicular to the segment direction:
Ak=l(r(kl)+r((k−1)l));k=1, . . . ,N (3-3)
Each segment offers resistance only to stretching, with a stiffness, Gk, given by:
E=Elastic modulus of the tether material
The tension, Tk(t), in segment k is this stiffness multiplied by the segment (positive) strain. We also model a small extensional damping in the segment. This adds a term to the tension equal to Gkη multiplied by the strain rate, where η<<1 is a viscous damping coefficient. These quantities are now calculated explicitly in terms of the position vectors of the N+1 point masses.
The axial strain in segment k and the strain rate are given by:
δk=(|
{dot over (δ)}=(
Where
Then the tension, Tk(t), in segment k is given by:
Tk(t)=max{0,Gk(δk+η{dot over (δ)}k)} (3-7)
The elastic force per unit mass on the kth point mass is:
{circumflex over (l)}k=(
Note that {circumflex over (l)}k is the unit vector parallel to segment k.
We can now state the equations of motion for the tether system. We equate the inertial acceleration of mass k minus the gravitational acceleration; namely
to the sum of the elastic force per unit mass, (3-8.a), and the aerodynamic force per unit mass. Making use of (2-8) and (3-6), and setting
{dot over (v)}kx=2Ωvky−(Ω02−Ω2)xk+{circumflex over (x)}□
{dot over (v)}ky=−2Ωvkx−(
{dot over (v)}kz=−Ω02zk+{circumflex over (z)}□
where
4. Aerodynamic Forces on the Tether
To formulate the correct model for the aerodynamic forces acting on the tether during entry, an important parameter is the Knudsen number, Kn. This is defined as):
Kn=λ/D
λ=molecular mean free path
D=characteristic length scale of the tether (4-1.a-c)
[See further: J. A. Storch, “Aerodynamic Disturbances on Spacecraft in Free Molecular flow”. Aerospace Corporation Report No. TR-2003 (3397)-1, 17 Oct. 2002.]
In this case, the characteristic length scale is that which determines the largest gradients of the flow, namely the tether diameter. The nominal concept is to assume a uniform diameter of 1 mm. The Knudsen number classifies the flow across the spectrum from continuum to highly rarified flow. By convention, Kn<0.01 indicates the fully continuum regime, while flows with Kn>10 are free molecular. In free molecular flow, momentum transfer occurs when individual molecules arrive at a body surface and are subsequently reflected or re-emitted, such that during this process the molecules do not interact among themselves. The skipper system starts operation in Earth orbit and even for extremely low orbits, λ is almost 1 km. Thus the flow is certain to be in the free molecular regime initially.
To assess the limitations of this formulation, it is necessary to consider the variation of mean free path with altitude. λ is determined by:
Using temperature and pressure profiles according to the U.S. Standard Atmosphere, one obtains the altitude as a function of mean free path shown in
We can expect that the shape assumed by the tether centerline will have radii of curvature on the order of total length, and this is many orders of magnitude larger than r(s). Therefore, to calculate the aerodynamic forces per unit length on the tether in the vicinity of some arc length coordinate, s, one can approximate the tether as a right circular cylinder, as illustrated in
In order to calculate the aerodynamic forces within the free molecular flow theory, we need to introduce two surface interaction parameters:
where τi and τr are the incident and re-emitted tangential momentum fluxes, and pi, pr, pw are respectively the normal fluxes of momentum that are incident, re-emitted, and that which would be re-emitted if all molecules left in Maxwellian equilibrium with the surface. In this model, we shall deal mainly with the two extreme cases: perfectly diffuse reflection (σ=σ′=1), and perfectly specular reflection (σ=σ′=0).
With the above definitions, the normal and tangential force components, denoted by N and T, respectively, can be given in terms of the normal and tangential force coefficients:
N=ρar(s)|
Where ρa is the atmospheric mass density at the given altitude. This is calculated by interpolation of the U.S. Standard Atmosphere density versus altitude tables.
For free molecular flow around the right circular cylinder illustrated in
[See further: L. Talbot. “Free Molecule Flow Forces and Heat Transfer for an Infinite Circular Cylinder at an Angle of Attack”. University of California Institute of Engineering Research, Berkeley, Calif., Technical report No. HE-150-136, 26 Mar., 1956.] where I0(ξ2/2) and I1(ξ2/2) are modified Bessel functions of the first kind, Tw is the temperature of the tether surface, and:
ξ=S sin β
S=the speed ratio=|
and where Vav is the average speed of atmospheric molecules:
Vav=√{square root over (2RTa/M)}
Ta=Atmospheric temperature
M=Average molecular weight (4-7.a-c)
The above expressions are complex to implement. Moreover, through (4-5.a) the force equations are coupled to the thermal dynamics of the system. Fortunately, we can invoke the large speed ratio approximation. Even at temperature ˜27K (with M=28 g/mole) the average molecular speed is less than 500 m/s. At higher altitudes as can be judged from the U.S. Standard Atmosphere, the molecular speed is considerably less. On the other hand, during the most critical portion of atmospheric entry, the entry vehicle speed is of the order of several kilometers per second. Thus, except for very small angles of attack, both S and ξ are ˜10 or greater. Series expansions of the Bessel functions for large arguments yield:
Substituting these expressions into (4-5) and neglecting terms of order 1/S2 or smaller, we obtain:
‘Now, returning to (4-4), we express the force per unit length, denoted here by
And define the unit normal vector:
{circumflex over (n)}(s)=
Then combining (4-4) and (4-9) and making use of (4-10,11), we have:
Where the minus sign arises because
|
we obtain:
We conclude this section of the analysis by relating this model to the discretized tether model introduced in the last section. Using the definition {circumflex over (l)}k=(
This takes the place of {circumflex over (t)}(s). (4-10) are replaced by:
for k=2, . . . , N+1. Using (2-3) and (3-6), we find that
The total aerodynamic force on a segment of length l is applied to mass k. We replace r(s) in (4-14) by the average ½(r(kl)+r((k−1)l)). Thus the total force on mass k, using (4-15)-(4-17) is:
where in the last line, we have used definition (3-3). The force per unit mass,
In summary, equations (3-9), with (3-8) and (4-19) and associated definitions constitute the equations of motion of the system. These are incorporated into the “skippertow” simulation code, wherein the dynamical equations are integrated using a modified third-order Adams method. The integration method necessitates use of a smaller time step than the popular fourth-order Runge-Kutta algorithm, but requires only two derivative evaluations per time step. In the following sections of this model, we use the simulation code to examine various dynamics issues. This part of the model analysis closes with a summary of what the remainder of this analysis model must consider relating to the addition of fully coupled fore-body dynamic, formulation of the thermal dynamics of the system, and complete simulations of the entry trajectory in order to refine our choice of system parameters.
5. Effects of Initial Conditions and Surface Interaction Parameters
In this section the analytical model, simulation code is used to study the dynamical motion resulting from various initial tether deployment states in the high-altitude, nearly orbital portion of flight. This is done for the two extreme cases of the surface interaction parameters, namely: perfectly diffuse reflection (σ=σ′=1), and perfectly specular reflection (σ=σ′=0). The tether parameters assumed are:
L=20 km;r(s)=constant=0.5 mm
E=248×109N−m−2;ρT=2100 kg−m−3
η=1.0×10−3 (5-1.a,c)
The values for L and r are consistent with [1]. The elasticity modulus and mass density are approximately those of Kevlar. For future reference, we also stipulate the breaking tenacity, Ebreak, which when multiplied by the strain yields the stress at which failure occurs. For the stipulated material:
Ebreak=7.0×109N−m−2 (5-2)
In the following discussion, the tether is modeled with 40 segments:
N=40 (5-3)
For the initial skipper flight condition, we assume a 100 km altitude and that the towing speed is the orbit speed:
H=100 km
Ω=Ω0 (5-4.a,b)
We first study the tether behavior when the initial configuration is at a large angle to the stream direction (the x-axis). In particular, suppose that the tether, due to gravity gradient is aligned along the negative y-axis, as shown in
At the outset, (
The following table summarizes the maximum excursions along the y axis. It is evident that the departures from the x axis slowly damp out, but at a rate that is much less than exponential. This can be expected from the fact that in the specular reflection case, the magnitude of the aerodynamic forces depends only on the normal component of velocity, and the force acts only in the normal direction. There is negligible force when the tether is tangent to the stream because (4.5.b) shows that CT=0 when σ=0. Thus as the tether assumes a shape that is closer to the x axis, the stabilizing aerodynamic force is strongly reduced in magnitude. Also, for similar reasons, the whip-like curling motions are relatively weakly opposed by the aerodynamics since those portions of the curl that are parallel to the stream experience almost no opposing force.
Of greater importance to the use of the tether as an entry trajectory controller is the observation that the tension distribution becomes chaotic. Many slack regions develop and the tension at any single point is only intermittently positive. For example,
One way to ameliorate this problem is to choose materials that offer diffuse reflection of atmospheric molecules in free molecular flow. In this case, the magnitude of the aerodynamic force is still proportional to the normal component of stream velocity, but it is also proportional to stream velocity itself, and thus can act in both the normal and tangential directions.
On the other hand, chaotic motion near the free end still occurs and convergence to a straight shape along the stream direction is still very slow. Moreover, even when such convergence occurs, the angle of attack approaches zero, and the aerodynamic drag is very small (CT=0 in the large speed ratio approximation, but equal to √{square root over (π)}/S according to the precise formula (4.5.b)). Thus the tension at the attachment point is very small and sensitive to initial conditions. To remedy this difficulty, it is necessary to add an aft-body attached to the free end of the tether. To model this we increase the mass of the N+1 point mass and treat it as a sphere of given radius to obtain the aerodynamic model. More specifically, for point mass N+1, we set:
Where Aaft is the effective frontal area of the aft-body which is assumed to be a sphere, but it may be another more optimal shape as may be determined by further detailed models or by flight testing. The sphere is used in this model as an example and to simplify the simulation.
Assuming a two meter radius,
The convergence situation is summarized by
In this example, it still takes a half-hour for the maximum departure from the x-axis to fall below 1 km, even with a rather large aft-body. Because of this it is desirable to gradually deploy the tether behind the fore-body at a slow rate such that there is always some tension in the tether. In this case, when the towing speed is nearly the orbital speed, we can expect the deployed tether to be aligned with the x-axis. This method of operation is recommended. At the same time, if the gradual deployment method should fail and the tether deploy all at once, the action of the aft-body will still align it with the stream direction.
6. Tether Dynamics at Sub-Orbital Velocities
In this section of the mathematical model of the present invention, we examine tether motion when the fore-body travels at a speed that is below orbital speed. It is assumed that the foregoing recommendations are in place, i.e. specularly reflective materials, an aft-boy, and gradual deployment from the forebody. Based on the last assumption, we assume in the following calculations that the tether/aft-body are initially along the x axis.
In this case, Ω<Ω0, and the local acceleration of gravity, {tilde over (g)}, is nonzero. Thus we expect the tether to sag under this partial gravity. The tether will therefore present a nonzero angle-of-attack along its entire length, thereby increasing the tension at the attachment point. There are two questions to be examined. First, will the tether assume a static shape, without any occurrence of a Helmholtz or “flag fluttering” instability? Secondly, assuming it does settle into a static shape, how is the attachment point tension affected?
Rather than assume a large aft-body, let us consider the result if the size of the aft body is reduced to 1 meter diameter and suppose it is a mylar balloon of 1 mm thickness. Assuming a mass density of 1390 kg/m3, the aft body mass is 4.367 kg. Thus, using (3.2),
To provide an initial example, the fore-body speed is assumed to be 80% of the orbit speed, i.e., Ω=0.8Ω0. As
To see if this chaotic motion can be stabilized, let us revert to the case of a 2 m radius aft-body.
The foregoing results show there is a combination (or combinations) of altitude and towing speed that results in what we shall call flag flutter instability. While the amplitude of flutter does not grow beyond a certain point, there results rapidly changing tension distributions, which are not desirable for the entry trajectory control application. We also see that it is possible to select aft-body parameters that will stabilize the system. However the stabilizing design of the aft-body remains in question. Simulation studies of every possible flight condition are impractical. Clearly, analysis is needed to determine the general conditions under which the flag flutter instability occurs in order to attain a stable predictable design.
At this point in the analysis we have developed a fairly complete model of the dynamics of a tether being towed by a massive fore-body in steady horizontal flight, including aerodynamic forces due to free molecular flow. The dynamic simulations performed according to this model indicate:
a. To promote faster stabilization in the initial portion of atmospheric entry, the tether material should be capable of diffuse reflection of atmospheric molecules.
b. In order to prevent large uncertainties in initial entry, the tether should be deployed gradually a so that the initial shape is parallel to the stream direction.
c. To provide insensitivity to initial conditions, guarantee some initial tension and devise some mechanism whereby the tether settles into to a steady shape and tension distribution for all flight regimes, it is necessary to add an aft-body attached to the free end of the tether.
d. The aft-body design parameters must be chosen to prevent flag flutter instability.
This leads to the next part of the mathematical model analysis of the present invention Skipper deorbit system:
A. Analytical investigation of the conditions leading to flag-flutter instability, and determination of corresponding design rules.
B. Development of a fore-body model and coupling it to the tether dynamics in the simulation code.
C. Development of an aerodynamic heating model, and breaking conditions for the simulation code.
D. Exercise the complete simulation code in order to design the entry trajectory control law.
7. Equations of Motion for Stability Analysis
As mentioned in Section 6, above, this analysis of the present invention must determine the general conditions under which “flag flutter instability” occurs in order to attain a stable predictable design. To accomplish this, we first formulate the continuum dynamics problem by allowing N to approach infinity in the lumped mass model previously developed for the simulation code.
This process amounts to assembling Equations (3-19) with elastic forces given by (3-2), (3-5), (3-7), (3-8), and aerodynamic forces defined by (4-16), (4-19) and (5-5-), and taking the limit N→∞ while maintaining L=Nl fixed. In performing this limiting process, we identify a material point on the filament centerline by its arclength, s, from the attachment point when the string is un-stretched. The resulting equations of motion in the continuum formulation are:
with boundary conditions:
Here
of both time and arclength that are the continuum counterparts of the lumped mass position, link segment tension and stream velocity at the mass locations. [Compare this with D. Yong, “Strings, Chains and Ropes”. SLAM Review, Vol. 48, No. 4, pp 771-781, 2006.]
This analysis is closely related to the problem of a periodically driven chain hanging under gravity. This apparently simple system displays complex dynamics; including parametric resonance, instability, chaotic motion, and even self-knotting [See further: A. Belmonte, M. J. Shelley, S. T. Eldakar, and C. H. Wiggens, “Dynamic Patterns and Self-Knotting of a Driven Hanging Chain”. Physical Review Letters. Vol. 87, No. 11, 10 Sep. 2001, which displays some rather spectacular examples. This reference describes a motion similar to one observed in connection with a towed cable system:
“We also observe a state, occurring only at large L (small s), in which the chain simultaneously rotates and swings three dimensionally around its suspension point, which we term stellate as the chain traces a multipointed pattern during its motion (Reference
The hanging chain results indicate that one must carefully analyze the skipper system to ensure that unstable or chaotic motion is avoided.
Before embarking on this task, one can introduce reasonable approximations for the gravity and inertial terms in (7-1), (7-2), and (7-4.b). While {tilde over (g)} is typically a substantial fraction of one “g”, even given a tether length of 20 km, the terms −3Ω02y, and (Ω02−Ω2)
Making these simplifications, we get:
with boundary conditions:
In our stability analysis to follow, it is assumed that the atmospheric density is constant along the tether and equal to its value at the attachment point, despite the very large tether length of 20 km in the nominal design concept. This is justified on a number of counts. First, for the flight regimes, (H,V) wherein the attachment point tension is largest, the tether is swept back toward the x axis, so that the variation of altitude along the tether is a small fraction of the nominal length. Secondly, models of heating and breakage conditions in the sequel will show that a 20 km tether will experience extreme heating rate gradients and frequent breakage near the free end—thus necessitating the use of a shorter tether. A shorter tether will also be less sensitive to varying atmospheric conditions with altitude since it will “sample” variations over a smaller altitude range. Finally, a shorter tether, by reducing the spatial dependence of the density will render the dynamics more predictable, e.g. the equilibrium solution for given (H,V) takes a particularly simple form, as shown below.
In the following three sections of this mathematical analysis of the present invention, the problem without the aft body is considered. First, in Section 8, the equilibrium solution is obtained. Then, in Section 9, the dynamics of small departures from equilibrium are formulated, and stability conditions for this static state are obtained in Section 10. Subsequently the full problem, comprising Equations (7-5)-(7-8), is considered, viewing the aft body terms in (7-8.b) as a perturbation.
8. Equilibrium Solution
Here, all the time derivatives in (7-5)-(7-8) are set to zero, and we consider the case maf→∞ in (7-8.b) to obtain:
s=L:T=0 (8-1)
which replaces (7-8.b).
Under these conditions, Equations (7-5)-(7-8) become:
Where the subscript “s” indicates the quantity pertains to its value for the static equilibrium
configuration, and {circumflex over (t)}s denotes the unit tangent vector,
In this case, one can see that the tether equilibrium shape is a straight line, as suggested by the simulation results in Section 6. Let us proceed on the assumption that {circumflex over (t)}s is a constant vector. Let {circumflex over (n)}s ({circumflex over (n)}s□{circumflex over (t)}s=0, {circumflex over (n)}s□ŷ>0) be the unit normal vector to the filament curve. If θ is the angle between the x axis and the unit tangent vector, then one can write:
{circumflex over (t)}s=cos θ{circumflex over (x)}−sin θŷ,{circumflex over (n)}s=sin θ{circumflex over (x)}+cos θŷ (8-3.a,b)
Forming the dot product of (8-2.a) with {circumflex over (t)}s assuming {circumflex over (t)}s is a constant vector, and using (8-2.d) after integrating the expression for Ts(s) one obtains:
Ts(s)=rρaV2(L−s)[ cos θ+2γ] sin θ (8-4)
where we define:
and where the last line follows by substitution of the previous expressions for {tilde over (g)}, V, and V0. Next, taking the dot product of (8-2.a) with {circumflex over (n)}s and assuming {circumflex over (t)}s and hence {circumflex over (n)}s are constant gives:
sin2θ=2γ cos θ0 (8-6)
From (8-6) we obtain:
sin θ=√{square root over (2[−γ2+γ√{square root over (1+γ2)}])} (8-7)
Finally, substituting all the above results back into (8-2.a) gives d{circumflex over (t)}s/ds=0 identically, showing that the only unique solution for the equilibrium shape is a straight line.
Collecting results, we can summarize as follows:
9. Dynamics of Small Departures from Equilibrium
To investigate the stability of the equilibrium state found in the last section, we specialize the equations of motion, (7-5)-(7-8), to the case of small departures from equilibrium. Let:
Where u and w are of order ϵL, ϵ<<1. Unless otherwise noted, we evaluate the various terms in (7-5)-(7-8) up to contributions of order unity and order ϵ, but neglect higher order terms. This prescription leads to:
where we also treat the term
as or order ϵ since me material damping coefficient may be considered quite small. Substituting these expressions into (7-5) and noting that the time-independent terms cancel out by virtue of the static equilibrium solution, we obtain:
Note that the term
might appear to be of order ϵ2, but it will be seen that the coefficient multiplying it is quite large, and moreover, this term is the most critical to the stability issue since it represents a parametric excitation similar to that occurring in the vertically excited hanging chain problem. In addition, since Ts(L)=0 this term dominates in the region near the free end. Thus, ignoring the order ϵ2 terms except as noted, and taking the dot products of (9-3) with {circumflex over (t)}s and {circumflex over (n)}s gives the equations of motion for small departures from the static equilibrium:
with boundary conditions:
To clarify the parameter dependencies, we conclude the section with the dimensionless forms of (9-4) and (9-5). Define the dimensionless arclength and time as:
Also, introduce the dimensionless displacements from Equilibrium:
U=u/L,W=w/L (9-7.a,b)
And define the following parameters:
With these definitions, (9-4) and (9-5) can be written:
We see that aside from the small (but not unimportant) damping factor, λ, the system is driven by just two parameters, θ (alternately γ) and κ.
10. Stability Analysis Without an Aft Body
To gain insight, we perform eigen-analysis of various operators in (9-9) to reduce these relations to ordinary differential equations for the various mode coefficients. First consider the eigen-problem associated with the terms
in (9-9.a). Setting W=Ŵ(ξ)exp(iωτ), this takes the form:
From which we deduce:
Where
is the Bessel motion of the first kind, of order cos2 θ and C1 and C2 are constants. Since, the second term approaches infinity as ξ approaches zero, C2=0. We choose C1 below so as to simplify the orthonormality conditions. Further we must satisfy the boundary condition Ŵ(ξ=1)=0. This yields:
In summary, the eigen-solution for (10-1) is:
With the choice of the coefficient in (10-4.a), the Ŵks obey the following orthonormality condition:
∫01Ŵk(ξ)ξcos
Next, in connection with (9-9.b), consider the eigen-problem
κ2Û″(ξ)=−
W(ξ,τ) and U(ξ,τ) may be expressed as:
Now we substitute (10-7.a) and (10-7.b) into (9-9.a), and (9-9.b), respectively; multiply by Ŵk(ξ) and Ûk(ξ), integrate over ξϵ[0,1], and make use of the orthonormality conditions (10-5) and (10-6.d) to obtain:
where:
The above equations give some insight regarding the possible modes of instability. As in the case of the vertically forced hanging chain, the small amplitude lateral motion is expressed as a Bessel series, where, in the absence of forcing Bn=0 in (10-8.a), the mode coefficients in the expansion behave as damped oscillators. In the hanging chain problem, the forcing term, of the formed ϵω2 Ak cos ωτ, replaces the last term on the right in (10-8.a), making it a Mathieu equation with viscous damping. For the Mathieu equation, and for each mode, there exist an infinite number of unstable driving frequencies which branch into tongues of parametric instability in the (ϵ,ω) plane [7]. Note that because it couples the lateral motion to the longitudinal dynamics, the term
on the right of (10-8.a) can play the roll (at least for a time) of a harmonic forcing term. If there are small initial deflections such that AmBnΓkmn≠0 and the combination frequencies ωm±{circumflex over (ω)}n≈ωk, then this term initially appears as a harmonic disturbance near the resonant frequency of the kth mode and we may expect initial growth in amplitude as in the hanging chain. Subsequently, the growing amplitude, Ak, further serves to excite the longitudinal modes via the coupling term on the right side of (10-9.b). This feedback effect raises the possibility of a true instability with initial exponential growth. Note, however, in distinction from the hanging chain case, the present model does not permit negative tension. Thus amplitudes cannot grow indefinitely, but must ultimately produce regions of slack, whereupon the system breaks into chaotic motion with limited amplitudes but highly variable tension distributions.
To obtain specific results, recourse must be had to the simulation code. Here we sample a grid of points in the (H,V/V0) plane, in increments of 0.1 along the V/V0 axis and 10 km along the H axis. The tether is initially in its equilibrium position plus a random, uniformly distributed deviation at each lumped mass with maximum variation of one km. For completeness, altitudes down to 20 km were investigated, although fully continuum flow can be anticipated at below 30 km. Above 170 km, the aerodynamic forces are so weak that the equilibrium shape is near vertical and the tension is very small.
By 646 s (
11. Stability Analysis Aft Body Design
Now consider the addition of an aft body to the present invention Skipper tether deorbit system. The objective here is to determine if a reasonably small aft body can prevent the instabilities and chaotic motions noted above. As in Section 5, we assume the aft-body to be a thin-walled sphere with approximately the same mass density as the tether. To model this we increase the mass of the N+1 point mass by maft, where:
maft=4πraft2haftρT
raft=Aft-body radius
haft=Aft-body sphere thickness≈r (11-1.a-c)
Also, we employ the aerodynamic force model given by (5-5), i.e.:
With these assumptions, we revisit (7-5) through (7-8), the only new feature being (7-8.b):
The first task is to see how the revised end condition affects the equilibrium solution. Setting all time derivatives to zeros, the revised version of (8-2) is:
Note that the normal components of (11-4.a) and (11-4.d) express the equilibration of the lifting aerodynamic forces and the downward local gravity force. For a reasonably sized aft-body, the contribution of the aft-body to the equilibrium in the normal direction, which is expressed by the normal component of (11-4.d) is negligible. Therefore this component may be discarded and only the tangential component retained:
s=L:T=ρaV2Aaf{circumflex over (x)}□{circumflex over (t)}s−{tilde over (g)}mafŷ□{circumflex over (t)}s (11-5)
This ultimately implies that the unit tangent vector, {circumflex over (t)}s, is constant as in the previous case. This is also consistent with the results shown in Section 5 assuming a rather large aft-body.
Following Section 8, because of (11-5), equation (8.4) must be modified by the addition of a nonzero tension at the aft end:
Taking the dot product of (11-4.a) we recover (8.7) once again. In summary, the equilibrium solution is:
The partial differential equations, (9-4) and boundary conditions (9-5) describing small departures from equilibrium remain unchanged, except that Ts(s) is now given by (11-7.d). Denote by
Then, adopting the dimensionless variables (9-6) and (9-7), and defining parameters as in (9-8), we find that the dimensionless equations of motion can be written:
The presence of
such that the term
is completely dominated by
Thus, for analysis of local stability involving sufficiently small initial departures from equilibrium, the term
may indeed be considered of order ϵ2 and dropped from (11-9.a). Therefore for local stability, we consider:
Adopting the same procedure as in Section 10, we address the eigen-problem:
By setting χ=ξ tan θ+
where again, vk is the kth zero of Jcos
With the above results, we introduce the mode expansions, (10-7):
and derive equations analogous to (10-8):
The coefficients of the lateral modes behave as uncoupled, damped oscillators and are obviously exponentially stable. For sufficiently long time, the coupling term in (11-14.b) becomes negligible so that the longitudinal mode coefficient also become uncoupled, damped oscillators. Thus local stability can be concluded for the tether with aft-body as a matter of general principle. We verify the above prediction by assuming a 0.5 m radius aft-body and spot checking the (V/V0, H) plane instability regions determined previously with no aft-body. Starting each case with a several kilometer deviation from the equilibrium shape, it is found that all cases are stable. An example is provided by the case H=100 km, V/V0=0.85 (marked by the red cross in
The foregoing results show that with a modestly sized, 0.5 meter radius aft-body (in contrast to the 2 m radius aft-body postulated in Section 6), and assuming initial departures of 3-4 km or less in the initial deviation from the equilibrium, the tether/aft-body system will converge to a stabile equilibrium such is not the case. This does not mean that the equilibrium state will be stable for larger initial deviations. However, the results do show that if initial deployment is slow and smooth and if altitude does not vary very rapidly, the system can stay close to a stable equilibrium shape and tension during the course of entry. Because the equilibrium state is relatively simple and predictable, this is a great advantage for the development of a robust, dependable entry trajectory controller.
We determined the equilibrium condition of the tether during steady flight under free molecular flow aerodynamics. The stability of this state was investigated and insight was achieved concerning the conditions for instability and chaotic motion. The following recommendations are offered:
1. A spherical aft-body having radius on the order of 0.5 meter appears to avoid chaotic motion at least for modest initial deviations from the equilibrium state
2. For flight regimes in which the variations of atmospheric properties over the tether length are small, the equilibrium shape is extremely simple and relatively insensitive to temporal and spatial variations in atmospheric density. To secure this advantage over all altitude and velocity ranges, the tether length should be reduced. Moreover, in Section 7, we argued that models of heating and breakage conditions to be explored in the sequel will show that a 20 km tether will experience extreme heating rate gradients and frequent breakage near the free end. In other words, if we do not shorten the tether, Mother Nature will do it for us! Since the key parameter determining tension is rL there is plenty of latitude for reducing the length while increasing the radius.
The next section of this model of the present invention Skipper tether deorbit system considers the following points:
1. Development of a fore-body model and coupling it to the tether dynamics in the simulation code.
2. Development of an aerodynamic heating model, and breaking conditions for the simulation code.
3. Exercise the complete simulation code in order to design the entry trajectory control law.
12. Kinematics of Coupled Tether/Fore-body Motion
So far, this model of the present invention has considered the motion of the tether when it is towed by an infinitely massive fore-body moving at constant speed and altitude. The model now extends this earlier formulation to include coupled fore-body/tether dynamics. This development begins with consideration of the kinematics of the fore-body motion.
As
If
Likewise, the inertial acceleration of the point is given by:
As in section 2, we determine the difference between the inertial velocity and the gravitational acceleration,
Subtracting this from (12-2), we obtain:
Note that if
Since the spatial extent of the system is much smaller than (RE+H), we can approximate the expression (12.4) by expanding it out to linear terms in x, y, or z and neglecting higher order terms:
If
{dot over (v)}kx={dot over (Ω)}(yk+RE+H)−{dot over (Λ)}zk+2Ω(vky+{dot over (H)})−2Λvkz+Λ2xk−(Ω02−Ω2)xk+{circumflex over (x)}□
{dot over (v)}ky=−{umlaut over (H)}−{dot over (Ω)}xk−2Ωvkx−ΩΛzk−{tilde over (g)}+3Ω02yk−(Ω02−Ω2)yk+ŷ□
{dot over (v)}kz={dot over (Λ)}xk+2Λvkx−ΛΩ(yk+RE+H)+Λ2zk−Ω02zk+{circumflex over (z)}□
The elastic and aerodynamic forces are the same as previously, except that (12-1) is now used to determine particle velocities with respect to the atmosphere. In addition, treatment of the aft-body remains the same. Recall that to model this we increase the mass of the N+1 point mass and treat it as a sphere of given radius to obtain the aerodynamic model given by (5-5), where (12.1) now determines the velocity,
In order to complete the treatment of coupled motion kinematics, and because the fore-body is treated as a rigid body, a system of body-fixed axes must be established. Let (xb,yb,zb) be orthogonal axes with origin at the center-of-mass of the fore-body (which is also the origin of the (x,y,z) system). Also, let ({circumflex over (x)}b,ŷb,{circumflex over (z)}b) be the triad of unit basis vectors parallel to the body-fixed axes. The triad ({circumflex over (x)}b,ŷb,{circumflex over (z)}b) is obtained from by (1) a rotation of ({circumflex over (x)},ŷ,{circumflex over (z)}) about the x axis by angle γ to produce ({circumflex over (x)}′,ŷ′,{circumflex over (z)}′); followed by (2) a rotation about the y′ axis by angle α to produce ({circumflex over (x)}″,ŷ″,{circumflex over (z)}″); concluded by (3) a rotation about the z″ axis by angle β to finally yield ({circumflex over (x)}b,ŷb,{circumflex over (z)}b). This rotation sequence is illustrated in
With the above definition, the unit basis vectors, ({circumflex over (x)}b,ŷb,{circumflex over (z)}b, can be related to ({circumflex over (x)},ŷ,{circumflex over (z)}) by:
Next, let us denote the angular velocity of the body-fixed frame relative to the (x,y,z) frame by
Referring to
Or, effecting the transformations to the ({circumflex over (x)}b,ŷb,{circumflex over (z)}b) basis:
ωx={dot over (γ)} cos α cos β+{circumflex over (α)} sin β
ωy=−{dot over (γ)} cos α sin β+{dot over (α)} cos β
ωz={dot over (γ)} sin α+{dot over (β)} (12-10.a-c)
We now invert these relations to obtain expressions for the derivatives of the yaw, roll and pitch angles in terms of the angular velocity components in the ({circumflex over (x)}b,ŷb,{circumflex over (z)}b) basis.
{dot over (α)}=ωx sin β+ωy cos β
{dot over (γ)}=(ωx cos β−ωy sin β)/cos α
{dot over (β)}=ωz−{dot over (γ)} sin α (12-11.a-c)
To conclude this section, we note that when in the following, a vector, Ā, is accompanied by a superscript, as in Ā(b), it means that its components are given in the ({circumflex over (x)}b,ŷb,{circumflex over (z)}b) basis.
13. Fore-Body Equations of Translational and Rotational Motion
To obtain the equations governing the translational motion of the fore-body center-of-mass, we
equate the quantity in (12-5) evaluated for
where Mf is the mass of the fore-body and
{circumflex over (l)}1=(
where,
The aerodynamic forces will be given in the next section.
Next we consider the rotational equations of motion for the fore-body. We assume that the (xb,yb,zb) axes are the principal axes of inertia, and write Euler's equations of rotational motion for the components of angular velocity in the (xb,yb,zb) frame. Note that the inertial angular velocity consists of
Let (Ψx,Ψy,Ψz) denote the components of this angular velocity in the (xb,yb,zb) frame. Then:
Ψx=Λ(sin β cos γ+cos β sin α sin γ)+Ω(sin β sin γ−cos β sin α cos γ)
Ψy=Λ(cos β cos γ−sin β sin α sin γ)+Ω(cos β sin γ+sin β sin α cos γ)
Ψz=−Λ(cos α sin γ)+Ω cos α cos γ (13-4.a-c)
And we can state Euler's equations as:
Ix{dot over (ω)}x=−Ix{dot over (Ψ)}x−(Iz−Iy)(ωy+Ψy)(ωz+Ψz)+Txtether+Txaero
Iy{dot over (ω)}y=−Iy{dot over (Ψ)}y−(Ix−Iz)(ωx+Ψx)(ωz+Ψz)+Tytether+Tyaero
Iz{dot over (ω)}z=−Iz{dot over (Ψ)}z−(Iy−Ix)(ωx+Ψx)(ωy+Ψy)+Tztether+Tzaero (13-5.a-c)
where Ix, Iy, and Iz are the moments of inertia, and
where
14. Model of the Fore-body and Its Aerodynamics
A variety of different types of fore-body may be employed in the skipper system, this general model of the present invention postulates a generic fore-body characterized by the fewest possible parameters. In essence, as
cpx>0 (14-1.a,b)
where, in the nominal model, cpy=cpz=0.
For derivation of inertial, and aerodynamic properties we ignore the more detailed features in the sketch of
In addition, define the characteristic area, S, as the area of the forward face normal to the xb axis, and let A denote a characteristic aspect ratio. These quantities are defined by:
S=lylz
A=lx/ly>1 (14-3.a,b)
where we assume ly<lz<lx.
Since the dimensions of the fore-body will of the order of several meters, the flow is not in the free molecular flow regime. This would appear to necessitate a complicated continuum flow model. However, for the blunt body postulated here, one may assume complete boundary layer or shock separation at the edges of the various flat surfaces. Moreover, the pressure on a given surface is the result of pressure drag (or form drag) which is generated by the component of stream velocity normal to the surface.
The first step in computing the aerodynamic model is to evaluate the fore-body velocity relative to the atmosphere,
Or, in the body-fixed frame:
Under the pressure drag assumption the aerodynamic force on the forward rectangular face, for example, is given by:
In other words, the force on a given face is directed toward the inward-pointing normal to that face (note that the stream velocity relative to the fore-body is −
The aerodynamic force appearing in the equations of translational motion, (13-1.a-c) is simply the above expression transformed into the ({circumflex over (x)},ŷ,{circumflex over (z)}) basis. The torque,
Note that (14-7) and (14-8) constitute a very simple model that at least crudely approximates flow around a blunt body and has no dependence on Mach number.
Before leaving this topic, we comment on the qualitative aspects of the above aerodynamic model in terms of conventional parameters, such as lift and drag coefficients, etc. These aspects can best be illustrated by considering planar motion in which only the pitch angle differs from zero, and in which the vehicle velocity is along the x axis. In that case, if a is the angle-of-attack, then (14-7) gives the following expressions for the lift and drag components of force (denoted by L and D, respectively):
From these formulae it is clear that the conventional lift and drag coefficients are:
CL=(A−1)sin α cos α
CD=(cos2α+A sin2α) (14-10.a,b)
The drag polar is illustrated (for A=3) in
at 45 degrees angle-of-attack, at which point the drag coefficient is
These are qualitatively reasonable and generic features. Note also that the drag polar has the form of a circle with center at CL=0,
and radius
15. Preliminary Coupled Tether/Fore-body Simulation Results
Here we consider some preliminary results for the fully coupled system. In the cases shown, the length of the tether is reduced to a somewhat more realistic length of 2.0 km, while we increase the radius to 0.5 cm. This preserves the maximum frontal area presented to the airstream as in the previously considered 20 km length and 0.5 mm case. All material properties are kept the same. The fore-body parameters are assumed as follows:
lx=0.5 m,ly=2.0 m,lz=4.5 m
ρf=9000 kg/m3 (15-1.a-d)
The center-of-pressure position is given by (14-1) with:
cpx=2.0 m,cpy=cpz=0 (15-2.a-c)
Also, the tether attachment point location is assumed to be:
Regarding initial conditions, it is assumed that the tether starts in its equilibrium position in the x-y plane. Likewise, the fore-body initial position and velocity are within the x-y plane, i.e.:
α=0,β=β0,γ=0
ωx=0,ωy=0,ωz=ωz0
Λ=0 (15-4.a-g)
Under these conditions, motion is confined to the x-y plane. This is reasonable since our first order of business is to design the entry strategy in the initial vertical plane (the plane defined by the initial velocity and the center of the Earth). The relatively smaller deviations in heading angle can be investigated later.
First, consider the motion of the fore-body alone.
Next,
After 82.79 s, the fore-body/tether combination is still mostly vertical but the tether begins to bow out in the leeward direction (
The observed oscillations occur because the entry is steep and the resulting rate of change in atmospheric density is too rapid for the tether to assume its equilibrium shape. This is unfavorable because during half of each cycle, the fore-body is forced into a positive angle-of-attack with greater lift resulting. But during the next half cycle, this is reversed; the angle-of-attack is negative. The end result is that the two effects nearly cancel out the presence of the tether in altering the altitude time history. Thus, in our design considerations we must take care to ensure that the entry profile involves density variations that are slow enough to allow the tether to follow its equilibrium configuration fairly closely.
Another observation is that the tether appears to be quite efficient in redirecting the angle-of-attack of the fore-body to match the tether tangent direction near the attachment point. Now if the tether is maintained in its sagged equilibrium shape, without oscillation, the main effect is to impose a positive angle-of-attack on the fore-body. Moreover the moment arm of the tether tension about the fore-body center-of-mass is likely much larger than can be attained by adjusting the attachment location in the back of the fore-body. Thus, it is probable that changing the tether length might be a more effective control strategy than adjusting the attachment point.
1. Care must be taken to ensure that the entry profile involves density variations that are slow enough to allow the tether to follow its equilibrium configuration fairly closely.
2. The tether appears to be quite efficient in redirecting the angle-of-attack of the fore-body to match the tether tangent direction near the attachment point. It is probable that changing the tether length might be a more effective control strategy than adjusting the attachment point. The capability for controlled length change will be installed into the simulation code.
The following section of this model of the present invention examines:
Development of an aerodynamic heating model, and breaking conditions for the simulation code. Exercises the complete simulation code in order to design the entry trajectory control law.
16. Aerothermal Model of the Tether
By far the most fragile component of the present invention Skipper tether reentry system is likely to be tether itself. If the tether survives the aerodynamic heating during atmospheric entry, so will the spacecraft fore-body. Thus we concentrate attention in this part of the model of the present invention on modeling the aerothermal dynamics of the tether. As in Section 4 (and referring to
Under these conditions, the heat input, to a differential element of area, dA, is given by [4]:
where S is the speed ratio, Ta is the atmospheric temperature, M the average molecular weight, ρa the atmospheric density, y the ratio of specific heats (˜1.4), and Tw denotes the tether wall temperature. θ is the angle between
where dEi is the total energy input, dEr is the total energy flux of the molecules re-emitted from the surface, and dEw is the energy flux which would be re-emitted if all molecules departed in Maxwellian equilibrium at the suface temperature, Tw. In accordance with our earlier recommendation that the tether surface material approximate perfectly diffuse reflection, we assume α to be unity.
With the assumption S>>1 the exponential terms in (16-1) are negligible and the error function approaches unity. Therefore we obtain the approximation:
Hence the heat transfer per unit length is approximately:
Since |
Next, we note that the thermal inertia of the thin tether is likely to be very small so that the interior of the tether is nearly in thermal equilibrium with the wall temperature. In other words, the temperature is uniformly equal to Tw over a cross-section and is a function only of the arc length coordinate. Further, as a worst case we neglect the thermal conductivity of the tether so that there is no heat flow along the tether. This means that (16-5) is the only heat input and the temperature Tw evolves according to:
Where CpT is the constant pressure specific heat of the tether material. (In the following, we assume the tether material to have high specific heat typical of such materials as Kevlar, where CpT≈1400 J/kg−K. See [8].)
In terms of the discretized tether model, we replace {circumflex over (t)}(s) by {circumflex over (t)}k=(
17. Consolidation of the AeroThermal Dynamic Equations
With (16-6) and (16-7), we have completed formulation of the dynamical model and its corresponding simulation code. At this point it is necessary to collect together the dynamical equations and consolidate them into forms suitable for analysis and for the formulation of the skipper system control algorithms. Since the principal phenomena addressed by the system pertain to the motion within the initial orbit plane, we specialize the equations developed in Sections 12 through 14 to the case of planar motion. Furthermore, we have seen that if rates of change of altitude and velocity are relatively slow, then at any instant, the tether approximately conforms to its equilibrium configuration. Therefore, the tether parameters as they affect the interaction with the for-body are modeled as in Section 8 and equations (8-8.a-f).
We begin by specializing the fore-body translational equations of motion. For planar motion, Λ=0, so that (13-1.a-d) become simply:
{dot over (Ω)}=−2Ω{dot over (H)}/(RE+H)−(Fxtether+Fxaero)/(Mf(RE+H))
{umlaut over (H)}=−{tilde over (g)}+(Fytether+Fyaero)/Mf (17-1.a,b)
where the force components are given in the (x,y) coordinate frame. These are now explicitly stated. The tether force is;
Fxtether=Ts(0){circumflex over (t)}s□{circumflex over (x)}=Ts(0)cos(θ)
Fytether=Ts(0){circumflex over (t)}s□ŷ=−Ts(0)sin(θ) (17-2.a,b)
where {circumflex over (t)}s is the unit tangent vector to the tether at the attachment point, given by (8-8.b) and Ts(0) is the tension at the attachment point given by (8-8.d) at s=0. Note that θ is the angle of sag due to gravity, i.e. at the attachment the tether points downward as s increases. θ is the angle between the horizontal and the tangent vector.
In the body-fixed frame the aerodynamic forces greatly simplify to:
where (17-3.b) repeats (14-4). We need to express
Using this relation to transform
Substituting these relations into (17-3.a) and transforming the force into the (x, y) frame we obtain:
Next consider the rotational dynamics. For planar motion, the pitch angle β, is the only Euler angle, ωz is the only rotational component and (12-11.c) gives {dot over (β)}=ωz, while (13-4) gives Ψz=Ω. Hence (13-5) yields:
Recall from (13-6) that
Tztether=−xattachTs(0)sin(θ+β)−yattachTs(0)cos(θ+β) (17-8)
According to (14-8),
Now consider the tether equilibrium equations. In the context of general planar motion of the fore-body reentering spacecraft, the ratio
appearing in (8-8.a-f) is actually the ratio of the total initial value of the gravitational energy and kinetic energy per unit mass to its current value. Notice that this energy is given by
Assuming the system is initially in orbit at altitude H(0)=H0, E(0)=Ω02(RE+H0)2. Hence the ratio
must be replaced by:
At this point we assemble the explicit forms of all equations before introducing certain conventions and alternate variables that greatly simplify their form. Collecting relations (17-1), (17-2), (17-6)-(17-9) and (8-8) with replacement (17-10) we have:
Translational Equations
Rotational Equation
Tether Attachment Point Conditions
Armed with the above results, we introduce various conventions that simplify their form. First, let us describe the attachment point coordinates as shown in
Ω(RE+H)=V cos ψ
{dot over (H)}=V sin ψ (17-14.a,b)
Using 17-11 and the above conventions; and approximating ratios such as (RE+H)/(RE+H0) by unity, we get the following equations of translational motion in terms of V and ψ:
Likewise, the rotational equation and the attachment conditions become:
Lastly the heat flow equation becomes:
18. Analysis of the Dynamical Equations as Pertains to the Control Problem
With Equations (17-15) through (17-18), together with various simulation results, we can begin to gather some insight on how the tether can be used to aid atmospheric entry. The overall goal is to prolong the entry process so that heating rates on the fore-body are significantly reduced. To control entry we have at our disposal the tether length L which can be increased or decreased and the position of the attachment point. Let us also suppose that the system is equipped with sensors that can measure altitude and both horizontal and vertical components of velocity (so that both Ω and {dot over (H)} can be known). Further, it is assumed that the sag angle, θ, can be measured. First the “skipping” motion that occurs when starting from sufficiently high orbit, or entering with sufficiently small inclination needs to be considered. At high altitude, with very low atmospheric density, both the tether and aerodynamic forces are dominated by the first terms on the right sides of (17-15a,b). Neglecting these terms, and recognizing that both ψ and V−V0 are small, the linearized versions of these equations are:
This system has eigenvalues ±iV0/(RE+H0)=±iΩ0 with periodic solutions. However, this does not reflect orbit eccentricity since {dot over (H)}=V sin ψ, and hence H are periodic with period half that of the orbit; so this motion is an atmospheric skipping phenomenon. This has been observed in the simulation results shown in
One possibility for using the tether to reduce entry speed is to periodically deploy and retract the tether so as to stimulate skipping in order to periodically “cool off” the vehicle by repeatedly moving to higher altitude. To examine this consider the simulation of the system initially at orbital velocity for 230 km initial altitude. In this case we postulate a tether having 500 m maximum length and 2.5 cm radius.
Here, the tether is fully deployed when the rate of descent increases and is retracted to essentially zero length when altitude begins to increase. This certainly encourages skipping and the variable length tether prolongs the entry process by more than a factor of two. However, the entry of the fore-body without the tether takes slightly longer. Thus this strategy seems to have merely reduced the drag of the system to the level of the fore-body alone. This does not mean that varying the length cannot be beneficial, rather it indicates that the tether dimensions make the tether-induced drag dominate the drag of the system.
Indeed, if we take the ratio of the tether force to the aerodynamic force (both of which are aerodynamic in origin) in (17-15.a), and use (17-17) to evaluate T0, we get:
Thus, aside from the angle terms, the ratio is of the order of the maximum total frontal area of the tether, 2rL to that of the fore-body. Even in the above example, 2rL=25 m2, whereas the area of the fore-body is just a few meters squared. We notice that in (17-15.a) the tether force is essentially always negative so it produces only enhanced drag which tends to shorten the entry process. The ratio of forces is the same for the {dot over (ψ)} equation, (17-15.b). But ψ is negative so that the term
is also negative and the tether produces only drag and no lift—and that drag dominates for the tether dimensions so far assumed.
Clearly using the tether just to increase drag is no way to prolong entry. Instead, 2rL must be chosen so that the tether forces in the translational equations are dominated by the aerodynamic forces on the fore-body. Thus, for example, a length of 500 m and diameter of a millimeter might be more appropriate. In such a case, the system does not have appreciably more drag than the fore-body alone.
To see what strategy remains to us, consider the aerodynamic force in the {dot over (ψ)} equation. If the angle ψ+β can be made negative, the aero force is positive—a lift force that decreases the angle of descent. The key to making this happen is to use the tether to control the pitch angle and thus enable the system to generate appreciable lift.
With this view in mind, examine (17-16), the pitch rotational equation. Note that the term
tends to make the pitch axis follow the stream direction and this will tend to oppose the effect of the tether force. It is recommended that we ensure that ra(rL)>>cpxAS that the first term on the right in (17-16) dominates.
In accordance with the above recommendations and the fact that {dot over (Ω)} is extremely small, (17-16) become approximately:
Now suppose we use the knowledge of the pitch angle rate, and the sag angle θ to adjust the attachment point angle such that:
θa=−θ−β*−η{dot over (β)} (18-4.a,b)
Here β* is the desired or targeted pitch angle that we shall specify in terms of measurable quantities below, and η is a positive constant of order unity. Substituting into (18-3), we obtain:
This system has dynamics with a time scale of order
and this is generally a few tens of seconds and much faster than the evolution of the altitude or inclination angle (upon which β* will depend). This means that β* varies much more slowly in (18-5) than (β+Γ) and that β* may be treated as a constant in (18-5). In that case, β closely approaches the current value of β*(t), that is:
In other words, the control law (18-4) creates a stable, invariant manifold that is rapidly approached by the system such that β(t)≈β*(t).
The above discussion means that we may approximately replace β(t) in the translational equations of motion with β*(t). Hence (17-15.a,b) may be replaced by:
Now it is evident that one way to prolong the entry process is to choose β* so as to maximize the vertical force term in (18-7.b), for given V; namely the expression
This obviously is maximum for β*=−ψ−π/4, attaining the value
Thus our proposed tether attachment control law is as follows:
19. Preliminary Simulation Results for the Attachment Controller
To assess the suitability of (18-8), we compare two cases that feature entry from low Earth orbit. For both cases, (1) the tether is 500 m long and 2 mm in diameter, (2) in our fore-body model, we set A=10 so that the lift-to-drag ratio is approximately unity, (3) the system starts at orbital speed at 170 km altitude. In the uncontrolled system, θa is fixed at zero, while the controlled system employs (18-8).
The high altitude behavior is essentially the same in both cases because of the very low atmospheric density and the small frontal area of the tether.
The significant differences between the two cases occur once they reach approximately 110 km altitude.
In contrast, the controlled system (see
The consequences for entry heating are difficult to determine precisely since they depend on the details of the entry body. However (examine Equation (17-18)) the combination ρaSV3 at least establishes the order of magnitude of overall heating rate. Since the frontal area, S, is common to both cases, we take ρaV3 as a crude measure of comparative heating rate.
Finally, we consider the g-loads on the forebody structure, i.e. the inertial acceleration minus gravity.
We round out this section by considering initial entry conditions typical of entry from interplanetary space. Denoting the entry angle and speed by αentry and Ventry, respectively, we set initial conditions for Ω and {dot over (H)} as follows:
{dot over (H)}(0)=−Ventry sin αentry
Ω(0)=Ventry cos αentry (19-1.a,b)
In the results shown below, we take αentry=25°, and Ventry=12 km/s. The initial altitude is set far beyond the sensible atmosphere, at 350 km. Note that if αentry<12.5°, the vehicle skips and is sent on an escape trajectory.
As
If we inspect the velocity time histories, not surprisingly, the uncontrolled system (see
The above noted patterns are consistently repeated for the heating rate estimates. The uncontrolled system (
Finally, we compare the g-loads for the uncontrolled and controlled systems. The uncontrolled system (
This model developed an aerodynamic heating model, and breaking conditions for the simulation code model of the present invention. This analysis developed an aerodynamic heating model for the tether. Although inclusion of breaking conditions in the simulation code is relatively straightforward, we deferred this task. Instead, the equations of motion were consolidated and placed in forms suitable for analysis of the entry control problem. A preliminary control algorithm for moving the tether attachment point location was proposed. Simulation results show that de-orbit and heating rates are much gentler for the controlled system.
Prior to experimental testing of the present invention the inventors plan to include breaking conditions in the simulation code; examine thermal performance of the tether and investigate combined control of tether length as well as attachment point location.
While only certain features of the selected embodiments have been illustrated and described herein, many modifications and changes will occur to those skilled in the art. It is, therefore, to be understood that the appended claims are intended to cover all such modifications and changes as fall within the true spirit of the embodiments.
This application is a continuation in part of and claims the benefit of U.S. application Ser. No. 14/025,822 filed on Sep. 13, 2013, which application claims the benefit of Patent Cooperation Treaty application number US2012/028841 filed Mar. 12, 2012 and U.S. provisional application No. 61/801,092 filed Mar. 15, 2013, which applications are hereby incorporated by reference.
| Number | Name | Date | Kind |
|---|---|---|---|
| 3241142 | Raabe | Mar 1966 | A |
| 3582016 | Sherman | Jun 1971 | A |
| 4097010 | Colombo | Jun 1978 | A |
| 4580747 | Pearson | Apr 1986 | A |
| 4824051 | Engelking | Apr 1989 | A |
| 5310144 | Salvatore et al. | May 1994 | A |
| 8052092 | Atmur et al. | Nov 2011 | B2 |
| Entry |
|---|
| Penzo et al, Tethers in Space Handbook Second Edition, May 1989, NASA (Year: 1989). |
| Carroll et al, Tether for Small Space Application, 1995, AIAA/ASU small satellite conference (Year: 1995). |
| Gavira Izquierdo et al, The tether System Experiment Preparing for ESA's First Tether Mission, May 2000, ESA bulletin (Year: 2000). |
| Number | Date | Country | |
|---|---|---|---|
| 20150076287 A1 | Mar 2015 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 14025822 | Sep 2013 | US |
| Child | 14469941 | US |
