Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Impulse Approximation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

Comments on SEE: Comparative Advantages and and Experimental Consequences

The Satellite Energy Exchange experiment measures the periodic, near-miss encounters between a sheppard satellite and a small test body (satellite) in approximately the same orbit about a primary. Several important experimental requirements have been chosen to enhance capabilities: (a) The satellite be flown in a sun-synchronous orbit at an altitude of about 1350 Km, (b) Passive temperature system stabilized by spacecraft axial rotation with sunshade baffles at the end of the spacecraft, (c) Test bodies with different material composition be available for experiments, (d) The containment spacecraft fly about the sheppard mass in a zero-g environment whereas the test bodies, experience average zero-g environment over an orbital period, (e) Primary attitude and station-keeping uses magnetic field alignment plus micro-Newton thrusters such as Field Emission Electric Propulsion, and (f) Very low power (nW) laser tracking systems minimize impulse delivered to test bodies. With the above conditions, SEE has the capabilities: (1) Long duration (several years life-time) flight experiment (2) Long-term, active (with historical time record), self-calibration of satellite mass distribution (capsule geodesy) over lifetime of the spacecraft. (3) Novel passive thermal stabilization systems designed to attain cryogenic temperatures around 78K. (4) Novel spacecraft stabilization systems. (5) Ability to measure G to 1 part in 10(exp 6-7) depending on ultimate duration of experiment. (6) Ability to place limits on both temporal and spacial variations on G. (7) Ability to set experimental limits on the Post Newtonian parameters (PPN) alpha(2) and zeta(2). (8) Ability to measure (or place limits on) the non Einsteinian eccentricity of the Earth-Sun system (and the parameter alpha(1)) for long duration flight. (9) Ability to measure Delta((dot)-G)/G to 1 part in 10(exp 12-13). The MiniSTEP, competes in a limited way with Project SEE. It is designed to improve the measurement of the equivalence principle by seven orders of magnitude using active, low temperature (1.8 K) cooling for SQUID based, differential superconducting circuits. The experiment consists of a small cylinder concentrically located within a larger cylinder at its null gravitational point. The satellite is operated in zero-g mode using four differential accelerometers consisting to two test bodies of different material composition. The SQUIDS are needed to measure test body motion to precisions of 10(exp -18) over a four orbit period. The entire satellite moves in a very precise zero-g mode since the accelerometers are rigidly attached to the satellite. This limits the experiment to an approximately six month due to limitations on helium storage used in cryogenic cooling and thrust control to maintain the zero-g operation.

Smalley, Larry L.↗

Impact of space shuttle orbiter reentry on mesospheric NOx.

Shuttle orbiter reentry will produce large amounts of nitric oxide (NO) through shock heating at mesospheric altitudes. The effects of this reentry are modeled by considering an impulsive source of NO in time-dependent equations including photochemistry and transport. The model uses an odd-nitrogen approach similar to Strobel's. Parameterized flux boundary conditions are imposed at the mesopause and stratopause and concentrations of the odd nitrogen species (NO, NO2, NO+, N(4S) and N(2D) are found numerically by a finite difference scheme. A localized disturbance can last many hours, depending on the effectiveness of eddy transport. Although the model is one-dimensional, effects of horizontal transport are also included approximately.

Stolarski, R. S.↗

Impulsive acceleration and scatter-free transport of about 1 MeV per nucleon ions in (He-3)-rich solar particle events

Impulsive solar energetic particle (SEP) events are studied to: (1) describe a distinct class of SEP ion events observed in interplanetary space, and (2) test models of focused transport through detailed comparisons of numerical model prediction with the data. An attempt will also be made to describe the transport and scattering properties of the interplanetary medium during the times these events are observed and to derive source injection profiles in these events. ISEE 3 and Helios 1 magnetic field and plasma data are used to locate the approximate coronal connection points of the spacecraft to organize the particle anisotropy data and to constrain some free parameters in the modeling of flare events.

Mason, G. M.↗

Thermal elastic shock and its effect on TOPEX spacecraft attitude control

Thermal elastic shock (TES) is a twice per orbit impulsive disturbance torque experienced by low-Earth orbiting spacecraft. The fundamental equations used to model the TES disturbance torque for typical spacecraft appendages (e.g., solar arrays and antenna booms) are derived in detail. In particular, the attitude-pointing performance of the TOPEX spacecraft, when subjected to the TES disturbance, is analyzed using a three-axis nonlinear time-domain simulation. Results indicate that the TOPEX spacecraft could exceed its roll-axis attitude-control requirement during penumbral transitions, and remain in violation for approximately 150 sec each orbit until the umbra collapses. A localized active-control system is proposed as a solution to minimize and/or eliminate the degrading effects of the TES disturbance.

Zimbelman, Darrell F.↗

Computational Optimization of a Rotary Valved Pulse Combustor Concept

A resonant pulse combustor valve concept is introduced that utilizes two slotted, coaxial counterrotating discs for mechanical actuation at the combustor inlet. The intended test article for prototype demonstration is a small, 22 in. long, propane fueled laboratory combustor, flowing approximately 0.006 lbm/s of air. The objective is to develop an externally actuated (i.e., active) valve that yields better performance and longer life than the traditional internally actuated (i.e., passive), reed-type valve found on most pulse combustors. The rotary valve motion is optimized using an axisymmetric, two-dimensional computational fluid dynamic simulation with a domain that includes the valve as a moveable interior wall. Parameters such as slew rate, dwell period in the open position, and total closed period are varied using fuel specific impulse as the figure of merit. Variations in fuel injector location and air fuel ratio are also examined. Additionally, the performance impact of leakage from the rotary valve is quantified since leakage is endemic to the design. The optimized simulation results indicate that the rotary valve concept can deliver the desired performance attributes using disc rotational speeds and stresses that are well within the realm of modern materials. A preliminary mechanical valve design is included in the report.

pressure gain combustion↗

Computational Optimization of a Rotary Valved Pulse Combustor Concept

A resonant pulse combustor valve concept is introduced that utilizes two slotted, coaxial counterrotating discs for mechanical actuation at the combustor inlet. The intended test article for prototype demonstration is a small, 22 in. long, propane fueled laboratory combustor, flowing approximately 0.006 lbm/s of air. The objective is to develop an externally actuated (i.e., active) valve that yields better performance and longer life than the traditional internally actuated (i.e., passive), reed-type valve found on most pulse combustors. The rotary valve motion is optimized using an axisymmetric, two-dimensional computational fluid dynamic simulation with a domain that includes the valve as a moveable interior wall. Parameters such as slew rate, dwell period in the open position, and total closed period are varied using fuel specific impulse as the figure of merit. Variations in fuel injector location and air fuel ratio are also examined. Additionally, the performance impact of leakage from the rotary valve is quantified since leakage is endemic to the design. The optimized simulation results indicate that the rotary valve concept can deliver the desired performance attributes using disc rotational speeds and stresses that are well within the realm of modern materials. A preliminary mechanical valve design is included in the report.

pressure gain combustion↗

Pulse Detonation Engine Modeled

Pulse Detonation Engine Technology is currently being investigated at Glenn for both airbreathing and rocket propulsion applications. The potential for both mechanical simplicity and high efficiency due to the inherent near-constant-volume combustion process, may make Pulse Detonation Engines (PDE's) well suited for a number of mission profiles. Assessment of PDE cycles requires a simulation capability that is both fast and accurate. It should capture the essential physics of the system, yet run at speeds that allow parametric analysis. A quasi-one-dimensional, computational-fluid-dynamics-based simulation has been developed that may meet these requirements. The Euler equations of mass, momentum, and energy have been used along with a single reactive species transport equation, and submodels to account for dominant loss mechanisms (e.g., viscous losses, heat transfer, and valving) to successfully simulate PDE cycles. A high-resolution numerical integration scheme was chosen to capture the discontinuities associated with detonation, and robust boundary condition procedures were incorporated to accommodate flow reversals that may arise during a given cycle. The accompanying graphs compare experimentally measured and computed performance over a range of operating conditions for a particular PDE. Experimental data were supplied by Fred Schauer and Jeff Stutrud from the Air Force Research Laboratory at Wright-Patterson AFB and by Royce Bradley from Innovative Scientific Solutions, Inc. The left graph shows thrust and specific impulse, Isp, as functions of equivalence ratio for a PDE cycle in which the tube is completely filled with a detonable hydrogen/air mixture. The right graph shows thrust and specific impulse as functions of the fraction of the tube that is filled with a stoichiometric mixture of hydrogen and air. For both figures, the operating frequency was 16 Hz. The agreement between measured and computed values is quite good, both in terms of trend and magnitude. The error is under 10 percent everywhere except for the thrust value at an equivalence ratio of 0.8 in the left figure, where it is 14 percent. The simulation results shown were made using 200 numerical cells. Each cycle of the engine, approximately 0.06 sec, required 2.0 min of CPU time on a Sun Ultra2. The simulation is currently being used to analyze existing experiments, design new experiments, and predict performance in propulsion concepts where the PDE is a component (e.g., hybrid engines and combined cycles).

Paxson, Daniel E.↗

Deflection and Fragmentation of Near-Earth Asteroids

Collisions by near-earth asteroids or the nuclei of comets pose varying levels of threat to man. A relatively small object, approximately 100 m diameter, which might be found on an impact trajectory with a populated region of the Earth, could potentially be diverted with a velocity of approximately1 cm/sec from an Earth impacting trajectory by impact (at 12 km/sec) by a rocket launched, 10^2 to 10^3 Kg impactor. For larger bodies, the use of kinetic energy impactors appear impractical because of the larger mass requirement. For any size object, nuclear explosions appear to be more efficient, using either the prompt blow-off from neutron radiation, the impulse from ejecta of a near-surface explosion for deflection, or, least efficiently, as a fragmenting charge.

Harris, Alan W.↗

Evidence for impulsive ion acceleration during the 0312 UT flare of 1980 June 7

One of the basic problems concerning the physics of solar flares is related to the process which is responsible for the acceleration of both electrons and ions. It has been proposed that the acceleration process proceeds in two different phases. The first phase results in the acceleration of electrons to energies in the range from 10 to 100 KeV. Electrons and ions with energies exceeding 30 MeV are finally produced during the second phase. Attention is given to the observational evidence which shows that these two phases, if in fact they are separate, must operate within seconds of each other, and that this process must be able to repeat itself producing pulses of electrons and ions approximately every 10 seconds. The observations were made with the Gamma Ray Spectrometer (GRS) on the SMM satellite during the 0312 UT flare on June 7, 1980.

Forrest, D. J.↗

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane-parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of nonthermal electrons and X-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick H I, Ca II and Mg II transitions using the non-local thermodynamic equilibrium (non-LTE) prescription in the program MULTI. Hydrogen ionization rates from X-ray photoionization and collisional ionization by nonthermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in the pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring. We have computed the chromospheric radiation that results from a range of coronal heating rates, with particular emphasis on the widely observed diagnostic H(alpha). Our conclusion is that the H(alpha) fluxes and profiles actually observed in flares can only be produced under conditions of a low-pressure corona with strong beam heating. Therefore we suggest that H(alpha) in flares is produced primarily at the footprints of newly heated loops where significant evaporation has not yet occurred. As a single loop evolves in time, no matter how strong the heating rate may become, the H(alpha) flux will diminish as the corona becomes denser and hence more effective at stopping the beam. This prediction leads to several observable consequences regarding the spatial and temporal signatures of the X-ray and H(alpha) radiation during flares.

Hawley, Suzanne L.↗

Beam-induced pressure gradients in the early phase of proton-heated solar flares

The pressure gradient induced in a coronal loop by proton beam momentum deposition is calculated and compared with the thermal pressure gradient arising from nonuniform deposition of beam energy; it is assumed that the transfer of momentum and energy from beam to target occurs via the Coulomb interaciton. Results are presented for both a low mean energy and a high mean energy proton beam injected at the loop apex and characterized by a power-law energy spectrum. The present treatment takes account of the breakdown of the cold target approximation for the low-energy proton beam in the corona, where the thermal speed of target electrons exceeds the beam speed. It is found that proton beam momentum deposition plays a potentially significant role in flare dynamics only in the low mean energy case and only in the corona, where it may dominate the acceleration of target material for as long as several tens of seconds. This conclusion suggest that the presence of low-energy nonthermal protons may be inferred from velocity-sensitive coronal observations in the early impulsive phase.

Tamres, David H.↗

Stability of an oscillated fluid with a uniform density gradient

Instabilities in a fluid with a constant density gradient that is subject to arbitrarily oriented oscillatory accelerations are considered. With the Boussinesq approximation and for the case of an unbounded fluid, transformation to Lagrangian coordinates allows the reduction of the problem to an ordinary differential equation for each three-dimensional wavenumber. The problem has three parameters: the nondimensional amplitude R of the base-state oscillation, the nondimensional level of background steady acceleration, which for some cases can be represented in terms of a local (in time) Richardson number Ri, and the Prandtl number Pr. Some general bounds on stability are derived. For Pr = 1 closed-form solutions are found for impulse (delta function) accelerations and a general asymptotic solution is constructed for large R and general imposed accelerations. The asymptotic solution takes advantage of the fact that at large R wave growth is concentrated at 'zero points'. These are times when the effective vertical wavenumber passes through zero. Kelvin-Helmholtz instabilities are found to dominate at low R, while Rayleigh-Taylor instabilities dominate at high R. At high R, the uniform shear of the Kelvin-Helmholtz case tends to distort and weaken instability waves. With unsteady flows, Ri = 1/4 is no longer an instability limit. Significant instabilities have been found for sinusoidal forcing for Ri up to 0.6.

Jacqmin, David↗

Performance of Small (100-lb Thrust) Rocket Motors Using Coaxial Injection of Hydrazine and Nitrogen Tetroxide

An investigation was conducted on a small (approximately 100-lb thrust) rocket using coaxial injection of hydrazine and nitrogen tetroxide. Characteristic-velocity efficiencies of 94 percent of the theoretical shifting equilibrium value were obtained at a chamber pressure of about 300 pounds per square inch using a 21-tube injector and a combustion chamber characteristic length of 10 inches. Performance at lower chamber pressures could be improved by reducing contraction ratio and thereby increasing the combustion chamber length and injector pressure drop, which would tend to promote better mixing. Calculations based on experimental data showed a vacuum specific impulse of 305 seconds with a nozzle area ratio of 50.

Wasserbauer, Joseph F.↗

Characteristics of Whipple Shield Performance in the Shatter Regime

Between the onset of projectile fragmentation and the assumption of rear wall failure due to an impulsive load, multi-wall ballistic limit equations are linearly interpolated to provide reasonable yet conservative predictions of perforation thresholds with conveniently simple mathematics. Although low velocity and hypervelocity regime predictions are based on analytical expressions, there is no such scientific foundation for predictions in the intermediate (or shatter) regime. As the debris flux in low earth orbit (LEO) becomes increasingly dominated by manmade pollution, the profile of micrometeoroid and orbital debris (MMOD) risk shifts continually towards lower velocities. For the International Space Station (ISS), encounter velocities below 7 km/s now constitute approximately 50% of the penetration risk. Considering that the transition velocity from shatter to hypervelocity impact regimes described by common ballistic limit equations (e.g. new non-optimum Whipple shield equation [1]) occurs at 7 km/s, 50% of station risk is now calculated based on failure limit equations with little analytical foundation. To investigate projectile and shield behavior for impact conditions leading to projectile fragmentation and melt, a series of hypervelocity impact tests have been performed on aluminum Whipple shields. In the experiments projectile diameter, bumper thickness, and shield spacing were kept constant, while rear wall thickness was adjusted to determine spallation and perforation limits at various impact velocities and angles. The results, shown in Figure 1 for normal and 45 impacts, demonstrated behavior that was not sufficiently described by the simplified linear interpolation of the NNO equation (also shown in Figure 1). Hopkins et al. [2] investigated the performance of a nominally-identical aluminum Whipple shield, identifying the effects of phase change in the shatter regime. The results (conceptually represented in Figure 2) were found to agree well with those obtained in this study at normal incidence, suggesting that shielding performance in the shatter regime could be well described by considering more complex phase conditions than currently implemented in most BLEs. Furthermore, evidence of these phase effects were found in the oblique test results, providing the basis for an empirical description of these effects that can be applied in MMOD risk assessment software. In this paper, results of the impact experiments are presented, and characteristics of target damage are evaluated. A comparison of intermediate velocity impact failure mechanisms in current BLEs are discussed and compared to the findings of the experimental study. Risk assessment calculations have been made on a simplified structure using currently implemented penetration equations and predicted limits from the experimental program, and the variation in perceived mission risk is discussed. It was found that ballistic limit curves that explicitly incorporated phase change effects within the intermediate regime lead to a decrease in predicted MMOD risk for ISS-representative orbits. When considered for all Whipple-based shielding configurations onboard the ISS, intermediate phase change effects could lead to significant variations in predicted mission risk.

Ryan, Shannon↗

A circuit model for filament eruptions and two-ribbon flares

A circuit model is derived for solar filament eruptions and two-ribbon flares. In the model the filament is approximated as a line current and the current sheet as infinitely thin. The model reproduces the slow energy buildup and eruption of the filament and the energy dissipation in a current sheet at the top of postflare loops during the two-ribbon flare. The two circuits considered are that of the filament and its return current and that of the current sheet and its return current. These circuits are inductively coupled, and free energy stored in the filament in the pre-flare phase is found to be transferred to the sheet during the impulsive phase and rapidly dissipated there. In the solutions for the evolution of the filament current sheet system four phases are distinguished: (1) a slow energy buildup, (2) a 'metastable' state, (3) the eruptive phase, and (4) a postflare phase. These phases are described in detail.

Martens, P. C. H.↗

On optimal infinite impulse response edge detection filters

The authors outline the design of an optimal, computationally efficient, infinite impulse response edge detection filter. The optimal filter is computed based on Canny's high signal to noise ratio, good localization criteria, and a criterion on the spurious response of the filter to noise. An expression for the width of the filter, which is appropriate for infinite-length filters, is incorporated directly in the expression for spurious responses. The three criteria are maximized using the variational method and nonlinear constrained optimization. The optimal filter parameters are tabulated for various values of the filter performance criteria. A complete methodology for implementing the optimal filter using approximating recursive digital filtering is presented. The approximating recursive digital filter is separable into two linear filters operating in two orthogonal directions. The implementation is very simple and computationally efficient, has a constant time of execution for different sizes of the operator, and is readily amenable to real-time hardware implementation.

Sarkar, Sudeep↗

COMPASS Final Report: Radioisotope Electric Propulsion (REP) Centaur Orbiter New Frontiers Mission

Radioisotope Electric Propulsion (REP) has been shown in past studies to enable missions to outer planetary bodies including the orbiting of Centaur asteroids. Key to the feasibility for REP missions are long life, low power electric propulsion (EP) devices, low mass Radioisotope Power System (RPS) and light spacecraft (S/C) components. In order to determine the key parameters for EP devices to perform these REP missions a design study was completed to design an REP S/C to orbit a Centaur in a New Frontiers (NF) cost cap. The design shows that an orbiter using several long lived (approx.200 kg xenon (Xe) throughput), low power (approx.700 W) Hall thrusters teamed with six (150 W each) Advanced Stirling Radioisotope Generators (ASRG) can deliver 60 kg of science instruments to a Centaur in 10 yr within the NF cost cap. Optimal specific impulses (Isp) for the Hall thrusters were found to be around 2000 s with thruster efficiencies over 40 percent. Not only can the REP S/C enable orbiting a Centaur (when compared to an all chemical mission only capable of flybys) but the additional power from the REP system can be used to enhance science and simplify communications. The mission design detailed in this report is a Radioisotope Power System (RPS) powered EP science orbiter to the Centaur Thereus with arrival 10 yr after launch, ending in a 1 yr science mapping mission. Along the trajectory, approximately 1.5 yr into the mission, the REP S/C does a flyby of the Trojan asteroid Tlepolemus. The total (Delta)V of the trajectory is 8.9 km/s. The REP S/C is delivered to orbit on an Atlas 551 class launch vehicle with a Star 48 B solid rocket stage

Oleson, Steven R.↗

An Investigation Into Transecting Satellites in Future Space Traffic Management Scenarios

The number of satellites expected to populate the near-Earth space environment is set to dramatically increase in the coming decade as new large constellations are approved and deployed. Current strategies for placing new batches of these satellites on orbit often involve launching into an initial orbit, and then performing apogee raising maneuvers to reach a target altitude. Similarly, end-of-life planning for these constellation satellites can consist of de-orbit burns that lower perigee to permit disposal via re-entry. Both the raising and de-orbiting maneuvers can result in the individual satellites traveling in high-eccentricity orbits that have the potential to transect other spacecraft trajectories. While individual large constellations may be able to coexist in separate altitude and inclination bands, having thousands of satellites moving between these bands as new satellites are replaced and old ones are removed could pose additional collision risks. Similar concerns have been raised regarding the impact that large numbers of university-class CubeSats might have in terms of their overall collision risk, especially as these satellites typically do not have propulsion systems for active maneuvering. To assess the impact that transecting satellites might have to future space traffic management strategies, this study explored a variety of future realistic scenarios using a high-fidelity simulation tool. The model can simulate the orbit of tens of thousands of resident space objects (RSOs) simultaneously, to include active satellites, debris, rocket bodies, or even future hypothetical scenarios, using a realistic force model that incorporates non-spherical gravity, atmospheric drag, solar radiation pressure, and more. As the model is run forward in time, various statistics and meta-data are gathered on any predicted conjunction event, providing insight into the nature and frequency of potential collisions, e.g., what size are the two satellites, who operates the satellites, are they active or passive objects, etc. Additional customization is available in terms of how probability of collision is computed, and how the probability ellipsoids and screening volumes are determined. The simulation tool also allows for rule-based maneuvers for active satellites, e.g, given an advance conjunction “warning,” one or both of the satellites can maneuver to a safe distance. A wide range of maneuvers can be implemented using impulsive or low-thrust methods, and the latencies can also be varied, e.g., using maneuver lead times of 48 hours, 24 hours, or 12 hours. Validation of the simulation results is performed against current and historical datasets available, to include comparisons to prior conjunction data messages (CDMs), object properties (mass, volume, etc.), and two-line-element records from both public and internal sources. Using the simulation environment, an assessment on the general risks that transecting satellites might pose for hypothetical future space object environments will be presented. This will involve the simulation of approximately 50,000 new large constellation satellites, in addition to the existing catalog of approximately 20,000 known resident space objects (RSOs), over propagation periods of one month to one year. A description of the simulation methodologies, scenarios evaluated, and validation methods will also be presented, as well as a preliminary assessment of the effectiveness of several candidate maneuver strategies that have the potential to reduce collision risk between active satellites.

Brian Gunter↗