Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD simulation”

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 289 records · Page 16

Operator learning for energy-efficient building ventilation control with computational fluid dynamics simulation of a real-world classroom

Energy-efficient ventilation control plays an important role in reducing building energy consumption while ensuring occupant health and comfort. While Computational Fluid Dynamics (CFD) simulations provide detailed and physically accurate representations of indoor airflow, their high computational cost limits their use in real-time building control. In this work, we present a neural operator learning framework that combines the physical accuracy of CFD with the computational efficiency of machine learning to enable building ventilation control with the high-fidelity fluid dynamics models. Our method jointly optimizes the airflow supply rates and vent angles to reduce energy use and adhere to air quality constraints. We train an ensemble of neural operator transformer models to learn the mapping from building control actions to airflow fields using high-resolution CFD data. This learned neural operator is then embedded in an optimization-based control framework for building ventilation control. Experimental results show that our approach achieves significant energy savings compared to maximum airflow rate control, rule-based control, as well as data-driven control methods using spatially averaged CO 2 prediction and deep learning–based reduced-order models, while consistently maintaining safe indoor air quality. These results highlight the practicality and scalability of our method in maintaining energy efficiency and indoor air quality in real-world buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Tracking Debris Shed by a Space-Shuttle Launch Vehicle

The DEBRIS software predicts the trajectories of debris particles shed by a space-shuttle launch vehicle during ascent, to aid in assessing potential harm to the space-shuttle orbiter and crew. The user specifies the location of release and other initial conditions for a debris particle. DEBRIS tracks the particle within an overset grid system by means of a computational fluid dynamics (CFD) simulation of the local flow field and a ballistic simulation that takes account of the mass of the particle and its aerodynamic properties in the flow field. The computed particle trajectory is stored in a file to be post-processed by other software for viewing and analyzing the trajectory. DEBRIS supplants a prior debris tracking code that took .15 minutes to calculate a single particle trajectory: DEBRIS can calculate 1,000 trajectories in .20 seconds on a desktop computer. Other improvements over the prior code include adaptive time-stepping to ensure accuracy, forcing at least one step per grid cell to ensure resolution of all CFD-resolved flow features, ability to simulate rebound of debris from surfaces, extensive error checking, a builtin suite of test cases, and dynamic allocation of memory.

Stuart, Phillip C.↗

Cryogenic Extension of NASA Species Polynomials Using Hydrogen and Oxygen at Stoichiometry

The Helmholtz free energy for a mixture has been used in many chemistry and combustion studies for mixtures. An approach is described where polynomial fits intended for a thermally perfect ideal gases can be modified and implemented into a Helmholtz free energy model, such as those used by NIST in their REFPROP program. However, the limitation of the polynomial fits must be respected and should only used in a valid temperature interval where the flow field can be treated as a thermally perfect ideal gas when these species are present. This approach can allow a CFD simulation to have cryogenic regions of a flow field where high fidelity cryogenic species models are needed and allows for the use of species in other areas at higher temperatures where combustion species do not have such a model available or is not needed. Presented in the paper are CFD 1D detonation wave and species mass fraction simulation results which have cryogenic upstream conditions of Hydrogen and Oxygen at stoichiometry.

Combustion↗

Computational flowfields for static testing of powered hypersonic aftbody models

Two-dimensional CFD analyses are presented related to the ground testing of hypersonic, air-breathing models which feature scramjet exhaust flow simulation. CFD analysis indicates that it is possible to test aftbody powered hypersonic airbreather configurations in a static, pumped-down environment to obtain aftbody aerodynamic performance data.

Huebner, Lawrence D.↗

Analysis of Launch Vehicle Liftoff Debris: Historical Perspective from Space Shuttle and Application to Artemis I

Human exploration-class launch vehicles are inherently prone to debris due to the extreme environments generated during pre-launch operations, liftoff, and flight. The use of cryogenic propellants often requires thermal protection system (TPS) coatings, typically foam, to maintain the propellant conditions in the tank and prevent an accumulation ice on the external surface of the vehicle. Some ice growth is to be expected at umbilical interfaces, vents, flanges, or brackets where it is difficult to apply TPS. This ice may come loose at any time due to wind on the launch pad, structural vibration and acoustics after rocket ignition, or aerodynamic forces during flight. This phenomena is especially apparent on vehicles with no TPS, such as the Saturn V rockets used in the Apollo Program, see Figure 1. During propellant tanking, the thermal contraction of the underlying substrate may generate cracks in the TPS (Figure 1). Chunks of TPS can release due to the expansion of ingested gas from cryopumping or from aerodynamic forces if the crack creates an offset surface. Most foams will also have a certain amount of “popcorning” where small pieces of foam will pop off during flight because of the differential between the static surface pressure and the pressure of the gas trapped in the foam cell structure. There are a number of other coating or closeout materials that may be shed from the vehicle and become debris. During pre-launch operations and liftoff, the vehicle may also be exposed to debris originating from the launch pad or ground support equipment. This debris is separate from foreign object debris, or FOD, which is not intended to be present and is strictly controlled through operations and maintenance procedures. In this case, debris is generated from hardware and materials that are necessary for launch and are subject to the intense vibration, acoustics, and direct plume impingement of the launch environment. Examples include ice from umbilicals, tape and tie wraps that protect cables, and rust or corrosion from the launch platform. While NASA has historically been aware of debris as a potential issue that could cause a failure resulting in loss of mission, loss of vehicle, or loss of crew, the likelihood and severity of that risk was not always well understood or given sufficient weight in program and flight decisions. After the Space Shuttle Columbia accident (STS-107), the investigation found that foam TPS debris shed from the external tank was the proximate cause of the damage to the orbiter wing. Six previous observations of debris released from the foam ramp that covered the bipod connecting the forward end of the orbiter to the external tank resulted in minor changes or were determined to be accepted flight risks. Two occurrences of bipod ramp foam loss were not identified until the STS-107 investigation. Despite the damage inflicted by these debris strikes, the Shuttle Program Requirements Control Board deemed the vehicle safe to fly. During the Return to Flight effort following the Columbia disaster, NASA Engineering developed a process for the assessment of debris transport, impact, and damage tolerance to support independent assessments of risk by NASA Safety and Mission Assurance (S&MA). Under this system, each element (vehicle or ground system) defines a catalog of all expected debris based on launch history, component testing, or analysis. Debris transport analysis (DTA) is conducted using the debris catalog characteristics and potential flow transport mechanisms (e.g., vehicle aerodynamics, gravity, wind, plume-driven). The predicted debris impact locations and velocities are provided to the hardware owners, who use available test data and analysis to determine whether each component can withstand the impacts. In cases where the element hardware may be severely damaged or fail, the options are to mitigate the debris source through some change in design or operation, or to work with S&MA to try to characterize the probability of the impact and damage for program risk acceptance. Because of the differences in debris characteristics and transport, the DTA has been divided between the Liftoff and Ascent regimes. The development and application of Liftoff DTA methodology from the Shuttle Program to the current Artemis Program is the subject of this paper. Liftoff DTA covers the time from the start of pre-launch operations at the launch pad, up until the vehicle clears the launch tower and there is no longer any interaction with ground systems. Debris transport during this period is broadly classified as either gravity, wind, and plume-entrained (GWPE) or plume driven (PD). GWPE debris is generally lower speed, travelling in a forward-to-aft direction. PD transport includes flow features from the rocket ignition transient, as well as plume impingement and recirculation that occur as the vehicle lifts off the launch platform. In these cases, the debris typically moves in an aft-to-forward direction at higher speeds. The applicable transport mechanisms must be considered for each piece of debris depending on the material, and release location and time. For example, rust or metallic debris from the tower could fall (GWPE) and impact the vehicle before landing on the launch platform deck where it could be also be transported by plume impingement (PD). However, falling ice (GWPE) from an umbilical is unlikely to survive impact with the vehicle or launch platform and be available for PD transport. Modeling of debris transport is accomplished using a set of DTA tools which simulate debris trajectories subject to a reference frame acceleration (i.e., gravity) and aerodynamic drag. Where the trajectory encounters a solid surface, the debris is allowed to rebound with a specified coefficient of restitution. The drag is calculated by interpolating the fluid state at each point in the debris trajectory from high-fidelity computational fluid dynamics (CFD) simulations of the launch vehicle and pad. The CFD data may either be static (steady state or time averaged), typically for GWPE transport, or dynamic (time-accurate) for PD flow features like the ignition transient. Examples of the CFD flow field solutions for the Space Launch System (SLS) rocket and launch pad are shown in Figure 2. Typical SLS debris trajectory predictions from DTA are illustrated in Figure 3. The final version of this paper will include a more detailed examination of the Liftoff DTA process developed during the Shuttle Program, and how it has been augmented and applied to the SLS rocket under the Artemis Program. Comparisons with debris observations from the Artemis I launch will demonstrate validation of the tools and methodology.

Debris↗

The impact of CFD on development test facilities - A National Research Council projection

The results of a National Research Council study on the effect that advances in computational fluid dynamics (CFD) will have on conventional aeronautical ground testing are reported. Current CFD capabilities include the depiction of linearized inviscid flows and a boundary layer, initial use of Euler coordinates using supercomputers to automatically generate a grid, research and development on Reynolds-averaged Navier-Stokes (N-S) equations, and preliminary research on solutions to the full N-S equations. Improvements in the range of CFD usage is dependent on the development of more powerful supercomputers, exceeding even the projected abilities of the NASA Numerical Aerodynamic Simulator (1 BFLOP/sec). Full representation of the Re-averaged N-S equations will require over one million grid points, a computing level predicted to be available in 15 yr. Present capabilities allow identification of data anomalies, confirmation of data accuracy, and adequateness of model design in wind tunnel trials. Account can be taken of the wall effects and the Re in any flight regime during simulation. CFD can actually be more accurate than instrumented tests, since all points in a flow can be modeled with CFD, while they cannot all be monitored with instrumentation in a wind tunnel.

Korkegi, R. H.↗

Uncertainty Assessments of 2D and Axisymmetric Hypersonic Shock Wave - Turbulent Boundary Layer Interaction Simulations at Compression Corners

This paper is one of a series of five papers in a special session organized by the NASA Fundamental Aeronautics Program that addresses uncertainty assessments for CFD simulations in hypersonic flow. Simulations of a shock emanating from a compression corner and interacting with a fully developed turbulent boundary layer are evaluated herein. Mission relevant conditions at Mach 7 and Mach 14 are defined for a pre-compression ramp of a scramjet powered vehicle. Three compression angles are defined, the smallest to avoid separation losses and the largest to force a separated flow engaging more complicated flow physics. The Baldwin-Lomax and the Cebeci-Smith algebraic models, the one-equation Spalart-Allmaras model with the Catrix-Aupoix compressibility modification and two-equation models including Menter SST, Wilcox k-omega 98, and Wilcox k-omega 06 turbulence models are evaluated. Each model is fully defined herein to preclude any ambiguity regarding model implementation. Comparisons are made to existing experimental data and Van Driest theory to provide preliminary assessment of model form uncertainty. A set of coarse grained uncertainty metrics are defined to capture essential differences among turbulence models. Except for the inability of algebraic models to converge for some separated flows there is no clearly superior model as judged by these metrics. A preliminary metric for the numerical component of uncertainty in shock-turbulent-boundary-layer interactions at compression corners sufficiently steep to cause separation is defined as 55%. This value is a median of differences with experimental data averaged for peak pressure and heating and for extent of separation captured in new, grid-converged solutions presented here. This value is consistent with existing results in a literature review of hypersonic shock-turbulent-boundary-layer interactions by Roy and Blottner and with more recent computations of MacLean.

Gnoffo, Peter A.↗

Development of a machine learning model for polyethylene pyrolysis using a detailed reaction mechanism

Waste plastics have recently received significant attention as the issue of waste generation continues to increase. Thermal conversion processes, such as pyrolysis and gasification, are attractive potential technologies for utilizing waste plastics and reducing overall waste generation. Efficient utilization of plastics requires a detailed understanding of the conversion process such as pyrolysis and gasification. However, a mechanistic understanding of these processes lead to large and complex kinetic schemes that are not suited for large-scale and long-time simulation methods. Currently, most modeling approaches for pyrolysis and gasification rely on globally lumped, simplified kinetic schemes that provide results that are classified by their product type and not individual species, which limit the level of fidelity achieved via modeling. A machine learning (ML) model has been developed for the primary reactions of high-density polyethylene (HDPE) in an attempt to increase computational efficiency while still maintaining a high level of detail and accuracy. The ML model is trained on a detailed reaction mechanism containing 42 total species and 737 chemical reactions. A DeepONet branch and trunk architecture was adopted to train the model using time-steps relevant to computational fluid dynamics simulations. The ML used physics-informed loss functions to ensure mass conservation. The surrogate model has been deployed in simple MFiX CFD simulations, single particle and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗

On the Numerical Prediction of Flow Transition

Knowledge from dynamical systems theory is used to study numerical uncertainties in direct numerical simulation of transition from laminar to turbulent flows. Spurious behavior due to underresolved grids and temporal discretizations is illustrated with 2-D CFD simulations.

Yee, H. C.↗

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

Experimental and Computational Study of Trapped Vortex Combustor Sector Rig with Tri-pass Diffuser

The Trapped Vortex Combustor (TVC) potentially offers numerous operational advantages over current production gas turbine engine combustors. These include lower weight, lower pollutant emissions, effective flame stabilization, high combustion efficiency, excellent high altitude relight capability, and operation in the lean burn or RQL (Rich burn/Quick mix/Lean burn) modes of combustion. The present work describes the operational principles of the TVC, and provides detailed performance data on a configuration featuring a tri-pass diffusion system. Performance data include EINOx (NO(sub x) emission index) results for various fuel-air ratios and combustor residence times, combustion efficiency as a function of combustor residence time, and combustor lean blow-out (LBO) performance. Computational fluid dynamics (CFD) simulations using liquid spray droplet evaporation and combustion modeling are performed and related to flow structures observed in photographs of the combustor. The CFD results are used to understand the aerodynamics and combustion features under different fueling conditions. Performance data acquired to date are favorable in comparison to conventional gas turbine combustors. Further testing over a wider range of fuel-air ratios, fuel flow splits, and pressure ratios is in progress to explore the TVC performance. In addition, alternate configurations for the upstream pressure feed, including bi-pass diffusion schemes, as well as variations on the fuel injection patterns, are currently in test and evaluation phases.

Hendricks, Robert C.↗

Experimental and Computational Study of Trapped Vortex Combustor Sector Rig With Tri-Pass Diffuser

The Trapped Vortex Combustor (TVC) potentially offers numerous operational advantages over current production gas turbine engine combustors. These include lower weight, lower pollutant emissions, effective flame stabilization, high combustion efficiency, excellent high altitude relight capability, and operation in the lean burn or RQL modes of combustion. The present work describes the operational principles of the TVC, and extends diffuser velocities toward choked flow and provides system performance data. Performance data include EINOx results for various fuel-air ratios and combustor residence times, combustion efficiency as a function of combustor residence time, and combustor lean blow-out (LBO) performance. Computational fluid dynamics (CFD) simulations using liquid spray droplet evaporation and combustion modeling are performed and related to flow structures observed in photographs of the combustor. The CFD results are used to understand the aerodynamics and combustion features under different fueling conditions. Performance data acquired to date are favorable compared to conventional gas turbine combustors. Further testing over a wider range of fuel-air ratios, fuel flow splits, and pressure ratios is in progress to explore the TVC performance. In addition, alternate configurations for the upstream pressure feed, including bi-pass diffusion schemes, as well as variations on the fuel injection patterns, are currently in test and evaluation phases.

Hendricks, R. C.↗

Evaluation of CFD to Determine Two-Dimensional Airfoil Characteristics for Rotorcraft Applications

The efficient prediction of helicopter rotor performance, vibratory loads, and aeroelastic properties still relies heavily on the use of comprehensive analysis codes by the rotorcraft industry. These comprehensive codes utilize look-up tables to provide two-dimensional aerodynamic characteristics. Typically these tables are comprised of a combination of wind tunnel data, empirical data and numerical analyses. The potential to rely more heavily on numerical computations based on Computational Fluid Dynamics (CFD) simulations has become more of a reality with the advent of faster computers and more sophisticated physical models. The ability of five different CFD codes applied independently to predict the lift, drag and pitching moments of rotor airfoils is examined for the SC1095 airfoil, which is utilized in the UH-60A main rotor. Extensive comparisons with the results of ten wind tunnel tests are performed. These CFD computations are found to be as good as experimental data in predicting many of the aerodynamic performance characteristics. Four turbulence models were examined (Baldwin-Lomax, Spalart-Allmaras, Menter SST, and k-omega).

Smith, Marilyn J.↗

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Drozda, Tomasz G.↗

Analysis of the Two-Level NO PLIF Model for Low-Temperature High-Speed Flow Applications

The current work compares experimentally obtained nitric oxide (NO) laser-induced fluorescence (LIF) spectra with the equivalent spectra obtained analytically. The experimental spectra are computed from captured images of fluorescence in a gas cell and from a laser sheet passing through the fuel-air mixing flowfield produced by a high-speed fuel injector. The fuel injector is a slender strut that is currently being studied as a part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center. This injector is placed downstream of a Mach 6 facility nozzle, which simulates the high Mach number airflow at the entrance of a scramjet combustor, and injects helium, which is used as an inert substitute for hydrogen fuel. Experimental planar (P) LIF is obtained by using a UV laser to excite fluorescence from the NO molecules that are present in either a gas cell or the facility air used for the EIMP experiments. The experimental data are obtained for several segments of the NO fluorescence spectrum. The selected segments encompass LIF lines with rotational quantum numbers appropriate for low-to-moderate temperature flows similar to those corresponding to the nominal experimental flow conditions. The experimental LIF spectra are then evaluated from the data and compared with those obtained from the theoretical models. The theoretical spectra are obtained from LIFBASE and LINUS software, and from a simplified version of the two-level fluorescence model. The equivalent analytic PLIF images are also obtained by applying only the simplified model to the results of the Reynolds-averaged simulations (RAS) of the mixing flowfield. Good agreement between the experimental and theoretical results provides increased confidence in both the simplified LIF modeling and CFD simulations for further investigations of high-speed injector performance using this approach.

Tomasz G Drozda↗

Combustion-Powered Actuation for Dynamic Stall Suppression - Simulations and Low-Mach Experiments

An investigation on dynamic-stall suppression capabilities of combustion-powered actuation (COMPACT) applied to a tabbed VR-12 airfoil is presented. In the first section, results from computational fluid dynamics (CFD) simulations carried out at Mach numbers from 0.3 to 0.5 are presented. Several geometric parameters are varied including the slot chordwise location and angle. Actuation pulse amplitude, frequency, and timing are also varied. The simulations suggest that cycle-averaged lift increases of approximately 4% and 8% with respect to the baseline airfoil are possible at Mach numbers of 0.4 and 0.3 for deep and near-deep dynamic-stall conditions. In the second section, static-stall results from low-speed wind-tunnel experiments are presented. Low-speed experiments and high-speed CFD suggest that slots oriented tangential to the airfoil surface produce stronger benefits than slots oriented normal to the chordline. Low-speed experiments confirm that chordwise slot locations suitable for Mach 0.3-0.4 stall suppression (based on CFD) will also be effective at lower Mach numbers.

Matalanis, Claude G.↗

Ways to Increase Launch Velocities of 2-Stage Gas Guns

The amount of space debris is rapidly increasing and the debris is distributed over a wide variety of orbits. Satellites, manned space vehicles and space stations will have to pay increasing attention to the dangers of impacts with space debris. Various armoring techniques (i.e., double or triple layer armor) will have to tested extensively to determine the most effective armor per unit weight. Intersecting near-earth orbits can lead to impact velocities up to 15 km/sec. Conventional two-stage light gas guns can launch intact, controlled-shape projectiles with a density of 1.2 gm/cc and length- to-diameter ratios of 0.5-1.0 at velocities up to 8-9 km/sec. Higher velocities (10-11 km/sec) can be obtained' for very light projectiles. The higher launch velocities tend to be very severe on the high pressure coupling and barrel of the gun and lead to short component lifetimes. Clearly, the ability to raise the launch velocity of a gun (for reasonably massive projectile shapes) from 8-9 km/sec to 11-13 km/sec (or higher), without reduction of component lifetimes, would have significant benefits. This would allow much better simulation of the higher velocity debris impacts as well as better simulation of high speed re-entry into planetary atmospheres. Several techniques for increasing the launcher muzzle velocity above 8-9 km/sec have been studied using CFD simulations and appear to offer the potential for significant gains. The first technique is to use multiple compressions, instead of a single compression, in the pump tube of the light gas gun. In a sense, this is a kind of pre-heating of the gas in the pump tube; other types of pre-heating have yielded disappointing results in the past. The dynamics of the multiple compression pump tube is very different, however, from the earlier techniques, where the pump tube was typically heated ohmically before the gun 2 cycle was started. In this paper, we present CFD calculations that show that significant increases in muzzle velocity can be obtained with multiple compressions in the pump tube. With a conventional two-stage gun, an important limitation to obtaining higher velocities is friction and heat transfer to the barrel, which typically has a length- to- diameter ratio of 200-400. These viscous losses greatly reduce the effectiveness of the regions of the barrel far removed from the second stage breech. We have studied computationally the effect of adding an additional breech (or breeches) along the barrel to reduce these viscous losses. Velocity increases from 6.5 to 7.2 km/sec have been obtained using the main breech and one additional breech. In these results, both breeches were operated with hydrogen, heated electrothermally. We have also studied a gun geometry where the main breech is operated in the conventional manner, using piston compression. The additional breech is operated either with electrothermal heating or heating by using a high explosive charge in a novel geometry. The latter option provides very effective compression, heating and acceleration of the hydrogen working gas and is fully reusable. Calculations are presented which show that very substantial increases in muzzle velocity can be obtained this way, without overstressing the projectile or the 'gun. The third technique studied is to add a section of ram accelerator tube after the barrel to further accelerate the projectile. The ram accelerator used here is not the conventional premixed gas ram accelerator, but a new technique using high explosive as the energy source and pure hydrogen as the working gas in a geometry which can be made fully reusable. Preliminary results with this new rain accelerator geometry were presented and showed that stable ram accelerator drive can be established. Herein, detailed calculations axe presented which show that substantial velocity increases can be obtained using this ram accelerator technique in tandem with a conventional light gas gun.

Bogdanoff, David W.↗