THEORY OF INITIATION OF DETONATION IN SOLID AND LIQUID EXPLOSIVES
Differences in growth of shock initiated detonation in liquid and solid explosives - energy release rate model for solid granular explosives
SEARCH · Engineering Papers
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.
Differences in growth of shock initiated detonation in liquid and solid explosives - energy release rate model for solid granular explosives
Mathematical models for solid bodies of planetary dimensions - models of Uranus and Neptune
This paper discusses the modelling by mechanical subelements, the general plasticity, the viscoplasticity, and the creep behavior of solids under multiaxial loading conditions. The formulation of a time-independent elastic-plastic analysis is based on the viscoplasticity theory and the assumed stress finite element method. An example of an in-plane stress problem is included.
Known models of Shuttle Solid Rocket Motor (SRM) performance have failed to produce pressure-time traces which accurately matched actual motor performance, especially during the first 5 seconds after ignition and during the last quarter of web burn time. Efforts to compensate for these differences in model reconstruction and actual performance resulted in resorting to the use of a Burning Anomaly Rate Function (BARF). It was suspected that propellant erosive burning was primarily responsible for the variation of model from actual results. The three dimensional Hercules Grain Design and Internal Ballistics Evaluation Program was made operational and slightly modified and an extensive trial and error effort was begun to test the hypothesis of erosive burning as an explanation of the burning anomaly. It was found that introduction of erosive burning (using Green's erosive burning equation) over portions of the aft segment grain and above a threshold gas Mach number did, in fact, give excellent agreement with the actual motor trace.
Four erosive burning models, equations (11) to (14). are developed in this work by using a power law relationship to correlate (1) the erosive burning ratio and the local velocity gradient at propellant surfaces; (2) the erosive burning ratio and the velocity gradient divided by centerline velocity; (3) the erosive burning difference and the local velocity gradient at propellant surfaces; and (4) the erosive burning difference and the velocity gradient divided by centerline velocity. These models depend on the local velocity gradient at the propellant surface (or the velocity gradient divided by centerline velocity) only and, unlike other empirical models, are independent of the motor size. It was argued that, since the erosive burning is a local phenomenon occurring near the surface of the solid propellant, the erosive burning ratio should be independent of the bore diameter if it is correlated with some local flow parameters such as the velocity gradient at the propellant surface. This seems to be true considering the good results obtained by applying these models, which are developed from the small size 5 inch CP tandem motor testing, to CFD simulations of much bigger motors.
This users manual for the solid rocket booster performance evaluation model (SRB-II) contains descriptions of the model, the program options, the required program inputs, the program output format and the program error messages. SRB-II is written in FORTRAN and is operational on both the IBM 370/155 and the MSFC UNIVAC 1108 computers.
A numerical model is developed to examine laminar flame spread and extinction over a thin solid fuel in lowspeed concurrent flows. The model provides a more precise fluid-mechanical description of the flame by incorporating an elliptic treatment of the upstream flame stabilization zone near the fuel burnout point. Parabolic equations are used to treat the downstream flame, which has a higher flow Reynolds number. The parabolic and elliptic regions are coupled smoothly by an appropriate matching of boundary conditions. The solid phase consists of an energy equation with surface radiative loss and a surface pyrolysis relation. Steady spread with constant flame and pyrolysis lengths is found possible for thin fuels and this facilitates the adoption of a moving coordinate system attached to the flame with the flame spread rate being an eigen value. Calculations are performed in purely forced flow in a range of velocities which are lower than those induced in a normal gravity buoyant environment. Both quenching and blowoff extinction are observed. The results show that as flow velocity or oxygen percentage is reduced, the flame spread rate, the pyrolysis length, and the flame length all decrease, as expected. The flame standoff distance from the solid and the reaction zone thickness, however, first increase with decreasing flow velocity, but eventually decrease very near the quenching extinction limit. The short, diffuse flames observed at low flow velocities and oxygen levels are consistent with available experimental data. The maximum flame temperature decreases slowly at first as flow velocity is reduced, then falls more steeply close to the quenching extinction limit. Low velocity quenching occurs as a result of heat loss. At low velocities, surface radiative loss becomes a significant fraction of the total combustion heat release. In addition, the shorter flame length causes an increase in the fraction of conduction downstream compared to conduction to the fuel. These heat losses lead to lower flame temperatures, and ultimately, extinction. This extinction mechanism differs from that of blowoff, where the flame is unable to be stabilized due to the high flow velocity.
A pattern search technique is presented, which is utilized in a computer program that minimizes the sum of the squares of the differences, at various times, between a desired thrust-time trace and that calculated with a special mathematical internal ballistics model of a solid propellant rocket motor. The program is demonstrated by matching the thrust-time trace obtained from static tests of the first Space Shuttle SRM starting with input values of 10 variables which are, in general, 10% different from the as-built SRM. It is concluded that an excellent match is obtained.
The experimental deflections for a transport type solid wing model were measured for several single point load conditions. These deflections were compared with those obtained by structural modeling of the wing by using plate and solid elements of Structural Performance Analysis and Redesign (SPAR) program. The solid element representation of the wing showed better agreement with the experimental deflections than the plate representation. The difference between the measured and calculated deflections is about 5 percent.
The James Webb Space Telescope Instrument Support Integration Module (ISIM) is being designed and developed at the Goddard Space Flight Center. The ISM Thermal Distortion Testing (ITDT) program was started with the primary objective to validate the ISM mechanical design process. The ITDT effort seeks to establish confidence and demonstrate the ability to predict thermal distortion in composite structures at cryogenic temperatures using solid element models. This-program's goal is to better ensure that ISIM meets all the mechanical and structural requirements by using test results to verify or improve structural modeling techniques. The first step to accomplish the ITDT objectives was to design, and then construct solid element models of a series 2-D test assemblies that represent critical building blocks of the ISIM structure. Second, the actual test assemblies consisting of composite tubes and invar end fittings were fabricated and tested for thermal distortion. This paper presents the development of the GSFC Cryo Distortion Measurement Facility (CDMF) to meet the requirements of the ISIM 2-D test. assemblies, and other future ISIM testing needs. The CDMF provides efficient cooling with both a single, and two-stage cryo-cooler. Temperature uniformity of the test assemblies during thermal transients and at steady state is accomplished by using sapphire windows for all of the optical ports on the radiation shields and by using .thermal straps to cool the test assemblies. Numerical thermal models of the test assemblies were used to predict the temperature uniformity of the parts during cooldown and at steady state. Results of these models are compared to actual temperature data from the tests. Temperature sensors with a 0.25K precision were used to insure that test assembly gradients did not exceed 2K lateral, and 4K axially. The thermal distortions of two assemblies were measured during six thermal cycles from 320K to 35K using laser interferometers. The standard deviation for all of the distortion measurements is less than 0.5 microns, which falls within the ISIM requirement of 3 microns.
Flow characteristics about submerged rocket nozzle model immersed in simulated exhaust gas
The results of a pressure test of a .00548 scale 146 inch Space Shuttle Solid Rocket Booster (SRB) with and without protuberances, conducted in a 14 x 14 inch trisonic wind tunnel are presented. Static pressure distributions for the SRB at reentry attitudes and flight conditions were obtained. Local longitudinal and ring pressure distributions are presented in tabulated form. Integration of the pressure data was performed. The test was conducted at Mach numbers of 0.40 to 4.45 over an angle of attack range from 60 to 185 degrees. Roll angles of 0, 45, 90 and 315 degrees were investigated. Reynolds numbers per foot varied for selected Mach numbers.
Significant error has been observed in the long term prediction of the Mean Local Time of the Ascending Node on the Aqua spacecraft. This error of approximately 90 seconds over a two year prediction is a complication in planning and timing of maneuvers for all members of the Earth Observing System Afternoon Constellation, which use Aqua's MLTAN as the reference for their inclination maneuvers. It was determined that the source of the prediction error was the lack of a solid Earth tide model in the operational force models. The Love Model of the solid Earth tide potential was used to derive analytic corrections to the inclination and right ascension of the ascending node of Aqua's Sun-synchronous orbit. Additionally, it was determined that the resonance between the Sun and orbit plane of the Sun-synchronous orbit is the primary driver of this error. The analytic corrections have been added to the operational force models for the Aqua spacecraft reducing the two-year 90-second error to less than 7 seconds.
All subprograms or routines associated with the solid rocket booster performance evaluation model are indexed in this computer listing. An alphanumeric list of each routine in the index is provided in a table of contents.
The relationship between solid deposit formation and liquid degradation product concentration was studied for the high temperature (400 C) stressing of three hydrocarbon model fuels. A Jet Fuel Thermal Oxidation Tester was used to simulate actual engine fuel system conditions. The effects of fuel type, dissolved oxygen concentration, and hot surface contact time (reaction time) were studied. Effects of reaction time and removal of dissolved oxygen on deposit formation were found to be different for n-dodecane and for 2-ethylnaphthalene. When ten percent tetralin is added to n-dodecane to give a simpler model of an actual jet fuel, the tetralin inhibits both the deposit formation and the degradation of n-dodecane. For 2-ethylnaphthalene primary product analyses indicate a possible self-inhibition at long reaction times of the secondary reactions which form the deposit precursors. The mechanism of the primary breakdown of these fuels is suggested and the primary products which participate in these precursor-forming reactions are identified. Some implications of the results to the thermal degradation of real jet fuels are given.
A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a plate or shell finite element model. Multi-point constraints provide a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with plate or shell finite elements. For simple double cantilever beam (DCB), end notched flexure (ENF), and single leg bending (SLB) specimens, mixed mode energy release rate distributions were computed across the width from nonlinear finite element analyses using the virtual crack closure technique. The analyses served to test the accuracy of the shell/3D technique for the pure mode I case (DCB), mode II case (ENF) and a mixed mode I/II case (SLB). Specimens with a unidirectional layup where the delamination is located between two 0 plies, as well as a multidirectional layup where the delamination is located between two non-zero degree plies, were simulated. For a local 3D model extending to a minimum of about three specimen thicknesses in front of and behind the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures modeled with plate elements, the shell/3D modeling technique offers a great potential, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.
Due to the severity of the internal solid rocket motor (SRM) environment, very few direct measurements of that environment exist; therefore, the appearance of such data provides a unique opportunity to assess current thermal/fluid modeling capabilities. As part of a previous study of SRM internal insulation performance, the internal thermal environment of a laboratory-scale SRM featuring aluminized propellant was characterized with two types of custom heat-flux calorimeters: one that measured the total heat flux to a graphite slab within the SRM chamber and another that measured the thermal radiation flux. Therefore, in the current study, a thermal/fluid model of this lab-scale SRM was constructed using ANSYS Fluent to predict not only the flow field structure within the SRM and the convective heat transfer to the interior walls, but also the resulting dispersion of alumina droplets and the radiative heat transfer to the interior walls. The dispersion of alumina droplets within the SRM chamber was determined by employing the Lagrangian discrete phase model that was fully coupled to the Eulerian gas-phase flow. The P1-approximation was engaged to model the radiative heat transfer through the SRM chamber where the radiative contributions of the gas phase were ignored and the aggregate radiative properties of the alumina dispersion were computed from the radiative properties of its individual constituent droplets, which were sourced from literature. The convective and radiative heat fluxes computed from the thermal/fluid model were then compared with those measured in the lab-scale SRM test firings and the modeling approach evaluated.
A shell/3D modeling technique was developed for which a local solid finite element model is used only in the immediate vicinity of the delamination front. The goal was to combine the accuracy of the full three-dimensional solution with the computational efficiency of a shell finite element model. Multi-point constraints provided a kinematically compatible interface between the local 3D model and the global structural model which has been meshed with shell finite elements. Double Cantilever Beam, End Notched Flexure, and Single Leg Bending specimens were analyzed first using full 3D finite element models to obtain reference solutions. Mixed mode strain energy release rate distributions were computed using the virtual crack closure technique. The analyses were repeated using the shell/3D technique to study the feasibility for pure mode I, mode II and mixed mode I/II cases. Specimens with a unidirectional layup and with a multidirectional layup were simulated. For a local 3D model, extending to a minimum of about three specimen thicknesses on either side of the delamination front, the results were in good agreement with mixed mode strain energy release rates obtained from computations where the entire specimen had been modeled with solid elements. For large built-up composite structures the shell/3D modeling technique offers a great potential for reducing the model size, since only a relatively small section in the vicinity of the delamination front needs to be modeled with solid elements.