The ''delta p'' approach from Nike combustion instability
Applicability of delta p approach to study of combustion instability in liquid and solid propellant rocket engines
SEARCH · Search NASA
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.
Applicability of delta p approach to study of combustion instability in liquid and solid propellant rocket engines
Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.
There is considerable evidence based on both full scale firings and cold flow simulations that hydrodynamically unstable shear flows in solid propellant rocket motors can lead to acoustic pressure fluctuations of significant amplitude. Although a comprehensive theoretical understanding of this problem does not yet exist, procedures were explored for generating useful analytical models describing the vortex shedding phenomenon and the mechanisms of coupling to the acoustic field in a rocket combustion chamber. Since combustion stability prediction procedures cannot be successful without incorporation of all acoustic gains and losses, it is clear that a vortex driving model comparable in quality to the analytical models currently employed to represent linear combustion instability must be formulated.
Cold flow tests of a one-tenth scale model based on the geometry of the Space Shuttle Solid Rocket Motor, SRM, head-end section and its single-port igniter, along with two models of a four port igniter, have been conducted. The tests were done to establish quantitative as well as qualitative data on the behavior of the flow inside the star slot during the early part of the ignition transient. This paper presents the data obtained from these tests and discusses the implications of the data with respect to ignition transients of solid rocket motors which have head-end star grains.
One major drawback of the eddy viscosity subgrid-scale stress models used in large-eddy simulations is their inability to represent correctly with a single universal constant different turbulent field in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work, a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model is based on an algebraic identity (Germano 1990) between the subgrid-scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid-scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The results of large-eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.
Preliminary results are presented from a test conducted to determine how well microphone phased array processing software could pull an acoustic signal out of background noise. The array consisted of 24 microphones in an aerodynamic fairing designed to be mounted in-flow. The processing was conducted using Functional Beam forming software developed by Optinav combined with cross spectral matrix subtraction. The test was conducted in the free-jet of the Nozzle Acoustic Test Rig at NASA GRC. The background noise was produced by the interaction of the free-jet flow with the solid surfaces in the flow. The acoustic signals were produced by acoustic drivers. The results show that the phased array processing was able to pull the acoustic signal out of the background noise provided the signal was no more than 20 dB below the background noise level measured using a conventional single microphone equipped with an aerodynamic forebody.
Ballistic and aerodynamic models are proposed to explain the dynamics of the volcanic plumes on Io. A ballistic model, which assumes that eruptive gas and solid-particle flows become completely and permanently decoupled over a small distance as compared with the maximum particle height, is shown to be unable to account for the observed shape of Plume 3. An aerodynamic model which assumes the entrainment and transportation of solid particles by the gas flow into the central portion of the top of the plume, with ballistic particle trajectories only after passage through the first part of their descent, is then developed. The model accounts for plume shape by the formation of a weak decompression shock near the edge of the plume stalk, which constrains horizontal expansion until a sufficient altitude is reached. Calculated plume dimensions are then compared with observations of Ionian Plumes 1 and 3.
The results are presented of an investigation of the unsteady pressures generated on a stator due to its interaction with the wakes shed by an upstream rotor. The influence of stator solidity, incidence flow angle and rotor/stator spacing are discussed. The results show that the major influence is due to stator solidity, particularly at large values of incidence angle. Comparisons of the measured data with an existing unsteady cascade analysis show similar trends in the chordwise variation of the predicted and measured unsteady pressure difference across the blades. Comparisons with an isolated airfoil analysis indicate the influence of solidity and unsteady blade-to-blade interaction. All comparisons were conducted for an incompressible flow with a reduced frequency of approximately 5.0.
Microgravity tests varying oxygen concentration and forced flow velocity have examined the importance of transport processes on flame spread over very thin solid fuels. Flame spread rates, solid phase temperature profiles and flame appearance for these tests are measured. A flame spread map is presented which indicates three distinct regions where different mechanisms control the flame spread process. In the near-quenching region (very low characteristic relative velocities) a new controlling mechanism for flame spread - oxidizer transport-limited chemical reaction - is proposed. In the near-limit, blowoff region, high opposed flow velocities impose residence time limitations on the flame spread process. A critical characteristic relative velocity line between the two near-limit regions defines conditions which result in maximum flammability both in terms of a peak flame spread rate and minimum oxygen concentration for steady burning. In the third region, away from both near-limit regions, the flame spread behavior, which can accurately be described by a thermal theory, is controlled by gas-phase conduction.
The Project objectives of this work were focused on the need to develop, investigate, and validate for field and commercial deployment robust pipeline surface treatments to prevent / minimize hydrate and other flow assurance solids deposition in subsea oil flowlines. A robust pipeline treatment that can prevent hydrate and wax/asphaltene adhesion/deposition will be a major fundamental breakthrough in hydrate and flow assurance science and engineering, and to critical Deepwater field operations. The project addressed the current major outstanding issue of mitigating gas hydrate and other solids deposition in different pipeline conditions, a major outstanding issue that has not been solved to-date, despite previous attempts/studies. The major goal of this project was to develop and investigate for deployment robust pipeline coatings to mitigate gas hydrate deposition in Deepwater oil flowlines, which will be critical in offshore leak and spill prevention. As oil and gas production wells mature, the water content and oil composition can change and increase the risk of gas hydrate, wax, and asphaltene deposition and blockages in flowlines, with the cost of complete inhibition becoming prohibitive. The ability to mitigate flow assurance risks in maturing flowlines can extend the life of the field, significantly reduce operational costs from inhibitor injection and prevent potential safety hazards to personnel and equipment and can prevent potential environmental hazards leading to catastrophic hydrocarbon emissions. This research allows for a novel, cost-effective method of mitigating gas hydrate and other pipeline solids deposition/blockages in flowlines, extending the life of the field while being minimally disruptive to normal flowline operation. The omniphobic surface treatment developed in this work can significantly improve the economics of energy transport by providing flow assurance and limiting catastrophic blowouts.
Contemporary discussions of planetary defense involve analyzing the risks posed by smaller sized, 20 to 200 m diameter, asteroids which are capable of breaking up in the atmosphere and generating a blast wave. Consequence assessments for this size class of asteroids are performed through fast-running analytic or semi-analytic models which are informed by high-fidelity hydrocode simulations of asteroid entry and breakup. However, insufficient historical data necessitates validating the independent physical processes which dominate airburst events. Here, the Fluid Solid Interface Smoothed Particle Hydrodynamics solver was previously used by Pearl et al. in 2023 to model the Chelyabinsk airburst and is used here to perform a series of validation simulations. The first effort involves modeling a cylinder in a hypersonic flow and comparing the bow shock geometry to that predicted by analytic theory. The second effort involves modeling the separation of two spherical bodies in supersonic flow and validating against experimental footage. Combined, these exercises demonstrate the ability of the code to model the flight-path of interacting solid bodies in a hypersonic flow.
Explore the source record for details and available documents.
A program was conducted to determine if a state of the art micro-machined silicon solid state flow sensor could be used to replace the existing space shuttle orbiter flow sensors. The rather aggressive goal was to obtain a new sensor which would also be a multi-gas sensor and operate over a much wider flow range and with a higher degree of accuracy than the existing sensors. Two types of sensors were tested. The first type was a venturi throat design and the second was a bypass design. The accuracy of venturi design was found to be marginally acceptable. The bypass sensor was much better although it still did not fully reach the accuracy goal. Two main problems were identified which would require further work.
The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.
A computational model of aluminum melting is proposed which captures both the thermal fluid-solid phase transition and the mechanical effects of oxidation. The model hybridizes ideas from smoothed particle hydrodynamics and bonded particle models to simulate both hydrodynamic flows and solid elasticity. Oxidation is represented by dynamically adding and deleting spring-like bonds between surface fluid particles to represent the formation and rupture of the oxide skin. Various complex systems are simulated to demonstrate the adaptability of the method and to illustrate the significant impact of skin properties on material flow. As a result, initial comparison to experiments of a melting aluminum cantilever highlights that the computational model can reproduce key qualitative features of aluminum relocation.
We report that the eddy current flow meter (ECFM) has been used to measure velocities and temperatures of conductive flows, such as liquid metal flows in a nuclear fission reactor. The goal of this paper is to develop a finite element electromagnetic model that can characterize the ECFM sensor performance and validate this finite element model with detailed velocity measurements of a controlled, well-characterized moving conductive solid rod. Both measurements and modeling were performed for various parameters that are important for ECFM performance such as rod velocity, rod material, ECFM sensor coil length, number of sensor coils, applied alternating current (AC) current amplitude, and applied ECFM AC frequency. For all parametric scans, the measurement and modeling agree well in both magnitude and trend. The normalized root-mean-square error between measurement and modeling is less than 10% for all cases. These results suggest that electromagnetic modeling could eventually be used to cost-effectively design future ECFM sensors in arbitrary geometry for more challenging applications such as liquid metal nuclear fission reactors.
For presentation at the AIChE Annual Meeting, Orlando, FL, November 5-10, 2023.
Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.