Numerical modeling for primary atomization of liquid jets
Explore the source record for details and available documents.
Chuech, S. G.↗
Engineering topics
Przekwas, A. J.
Publications and source records attributed to Przekwas, A. J..
Direct simulation for the instability and breakup of laminar liquid jets
A direct numerical simulation method is described for predicting the deformation of laminar liquid jets. In the present nonlinear direct simulation, the convective term, which has been discarded in past linear analyses by Rayleigh and others, is included in the hydrodynamic equations. It is shown that only by maintaining full complexity of the nonlinear surface tension term accurate drop formation can be predicted. The continuity and momentum equations in the transient form are integrated on an adaptive grid, conforming the jet and surface wave shape. The equations, which are parabolic in time and elliptic in space, are solved by a TVD scheme with characteristic flux splitting. The results of the present work are discussed and compared with available measurements and other analyses. The comparison shows that among the predictions, the current 1-D direct simulation results agree best with the experimental data. Furthermore, the computer time requirements are much (an order of magnitude) smaller than those of previously reported multidimensional analyses.
Chuech, S. G.↗
Assessment of shock capturing schemes for resonant flows in nonlinear instability analysis
The paper presents computational assessment of advanced numerical schemes for nonlinear acoustic problems related to combustion instabilities in liquid rocket engines. Several time-accurate, shock capturing schemes have been evaluated on a benchmark, closed-end resonant pipe flow problem. It involves the numerical solution of inviscid, compressible gas dynamics equations to predict acoustic wave propagation, wave steepening, formation of shocks, acoustic energy dissipation and wave-wall reflection for several hundred wave cycles. It was demonstrated that high accuracy TVD type schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time. The selected scheme was then applied to analyze the acoustic responses of resonant pipe-resonator, radial acoustic modes and hub-baffle configurations. Interesting observations of wave shape and damping characteristics have been drawn from presented computational studies.
Przekwas, A. J.↗
A comparative study of advanced shock-capturing schemes applied to Burgers' equation
Several variations of the TVD scheme, ENO scheme, FCT scheme, and geometrical schemes, such as MUSCL and PPM, are considered. A comparative study of these schemes as applied to the Burgers' equation is presented. The objective is to assess their performance for problems involving formation and propagation of shocks, shock collisions, and expansion of discontinuities.
Yang, H. Q.↗
Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems
The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.
Przekwas, A. J.↗
Numerical modeling for primary atomization of liquid jets
In the proposed numerical model for primary atomization, surface-wave dispersion equations are solved in conjunction with the jet-embedding technique of solving mean flow equations of a liquid jet. Linear and approximate nonlinear models have been considered. In each case, the dispersion equation is solved over the whole wavelength spectrum to predict drop sizes, frequency, and liquid-mass breakup rates without using any empirical constants. The present model has been applied to several low-speed and high-speed jets. For the high-speed case (the LOX/H2 coaxial injector of the Space Shuttle Main Engine Preburner), predicted drop sizes and liquid breakup rates are in good agreement with the results of the CICM code, which have been calibrated against measured data.
Przekwas, A. J.↗
Numerical grid generation and flow simulation in SSME thrust chamber
The development of liquid and solid rocket engines for future space projects demands a detailed optimization process for highly efficient performance and cost reasons. Also, testing of full size engines may not be feasible when the large size requires test facilities which are cost prohibitive or if vacuum operation cannot be acquired. For such situations only scaling from small test scale measurements or accurate analytical predictions will provide the performance prior to actually flying the mission. A rigorous approach for simulating the combustion processes in liquid rocket engines by employing a direct solution of Navier-Stokes equations within the entire volume of the thrust chambers is presented. This method is illustrated in the solution of reactive flow in the Space Shuttle Main Engine (SSME) thrust chamber. The objective is to review recent improvements in the mathematical model and to present the grid generation methodology suitable for rocket thrust chamber geometries.
Gross, K. W.↗
Numerical and analytical study of fluid dynamic forces in seals and bearings
A numerical model based on a transformed, conservative form of the three dimensional Navier-Stokes equation and an analytical model based on lumped fluid parameters are presented and compared with studies of modeled rotor bearing seal systems. The rotor destabilizing factors are related to the rotative character of the flow field. It is shown that these destabilizing factors can be reduced through a decrease in the fluid average circumferential velocity. However, the rotative character of the flow field is a complex three dimensional system with bifurcated secondary flow patterns that significantly alter the fluid circumferential velocity. By transforming the Navier-Stokes equations to those for a rotating observer and using the numerical code PHOENICS-84 with a nonorthogonal body fitted grid, several numerical experiments were carried out to demonstrate the character of this complex flow field. In general, fluid injection and/or preswirl of the flow field opposing the shaft rotation significantly intensified these secondary recirculation zones and thus reduced the average circumferential velocity; injection or preswirl in the direction of rotation significantly weakened these zones.
Tam, L. T.↗
A multidomain global modeling technique for analysis of Space Shuttle main engine
A multidomain global modeling technique is presented. The purpose of this multidomain methodology is to simultaneously analyze the three-dimensional compressible flow field in Space Shuttle main engine components. This technique provides flexibility and modularity for individual domain and also simplifies the grid generation task for each domain. A two-way coupling procedure is adopted for exchanging interface boundary conditions. This algorithm which requires that the adjacent domains share an overlapping grid slab, is fully conservative and implicit. For non-matching grid systems and for non-overlapping interfaces, an interpolation and/or extrapolation procedure is used when transferring the interface boundary conditions at the adjacent domains. To maintain overall conservation, a flux correction method is subsequently used. The coupling procedure was verified by modeling several geometries with features similar to SSME components. Satisfactory results for all test cases were obtained. Finally, this algorithm has been applied to the SSME flow simulation.
Tam, L. T.↗
Liquid oxygen sloshing in Space Shuttle External Tank
This paper describes a numerical simulation of the hydrodynamics within the liquid oxygen tank of the Space Shuttle External Tank during liftoff. Before liftoff, the tank is filled with liquid oxygen (LOX) to approximately 97 percent with the other 3 percent containing gaseous oxygen (GOX) and helium. During liftoff, LOX is drained from the bottom of the tank, and GOX is pumped into the tank's ullage volume. There is a delay of several seconds before the GOX reaches the tank which causes the ullage pressure to decrease for several seconds after liftoff; this pressure 'slump' is a common phenomenon in rocket propulsion. When four slosh baffles were removed from the tank, the ullage gas pressure dropped more rapidly than in all previous flights. The purpose of this analysis was to determine whether the removal of the baffles could have caused the increased pressure 'slump' by changing the LOX surface dynamics. The results show that the LOX surface undergoes very high vertical accelerations (up to 5 g) and, therefore, splashing almost certainly occurs. The number of baffles does not affect the surface if the structural motion is assumed; but, the number of baffles may affect the structural motion of the tank.
Kannapel, M. D.↗
Numerical and analytical study of fluid dynamic forces in seals and bearings
Study results for rotor/bearing/seal system models are compared with both a numerical model based on a transformed, conservative version of the three-dimensional Navier-Stokes equation, and an analytical model based on 'lumped' fluid parameters. Rotor-destabilizing factors are related to flow field rotation; it is shown that these factors can be reduced through a decrease in the fluid average circumferential velocity. A transformation of the Navier-Stokes equations to those for a rotating observer and a numerical code with a nonorthogonal body-fitted grid are used to carry out numerical experiments demonstrating the character of the complex flow field.
Tam, L. T.↗
Three-dimensional analysis of liquid oxygen sloshing in Space Shuttle external tank
A three-dimensional numerical simulation has been performed in order to investigate the hydrodynamics within the liquid oxygen tank of the Space Shuttle external tank after liftoff. The results indicate that the LOX surface undergoes very high vertical accelerations (up to 5 g) and that splashing almost certainly occurs. Although the number of slosh baffles is not found to affect the surface, it is noted that the number of baffles may affect the structural motion of the tank. Surface accelerations are similar to those observed in previous two-dimensional simulations.
Kannapel, M. D.↗
Computational simulation of liquid rocket injector anomalies
A computer model has been developed to analyze the three-dimensional two-phase reactive flows in liquid fueled rocket combustors. The model is designed to study the influence of liquid propellant injection nonuniformities on the flow pattern, combustion and heat transfer within the combustor. The Eulerian-Lagrangian approach for simulating polidisperse spray flow, evaporation and combustion has been used. Full coupling between the phases is accounted for. A nonorthogonal, body fitted coordinate system along with a conservative control volume formulation is employed. The physical models built into the model include a kappa-epsilon turbulence model, a two-step chemical reaction, and the six-flux radiation model. Semiempirical models are used to describe all interphase coupling terms as well as chemical reaction rates. The purpose of this study was to demonstrate an analytical capability to predict the effects of reactant injection nonuniformities (injection anomalies) on combustion and heat transfer within the rocket combustion chamber. The results show promising application of the model to comprehensive modeling of liquid propellant rocket engines.
Przekwas, A. J.↗
Numerical analysis of flow in the hot gas manifold of the Space Shuttle main engine
This paper presents a numerical model and results of analyses carried out to characterize the flow through the two duct configuration of the Hot Gas Manifold of the Space Shuttle Main Engine. Three dimensional computations have been performed for a half-scale air test model using a nonorthogonal body-fitted coordinate system. The calculation domain is extended from the inlet of the turnaround duct to the exit of the transfer duct. Three test cases, one for laminar flow and two for turbulent flow, have been considered. For turbulent flows, constant eddy viscosity and the k-epsilon model of turbulence were employed. As expected, laminar flow calculation predicts much larger pressure drop than turbulent flow cases. The turbulent flow results are in good agreement with the available flow-visualization data. This study and experimental data indicate that the two-transfer duct design will significantly improve the flow distribution in the Hot Gas Manifold and thereby enhance the overall performance of the SSME.
Owens, S. F.↗
SSME thrust chamber modeling with Navier Stokes equations
The capability of predicting two-dimensional, compressible and reacting flow in the combustion chamber and nozzle of the Space Shuttle Main Engine (SSME) is demonstrated. A nonorthogonal body fitted coordinate system has been used to represent the combustor and nozzle geometry. The Navier-Stokes equations are solved for the entire thrust chamber with the k-epsilon turbulence model accounting for compressibility and large pressure gradients effects. Results of the computational test cases reveal all expected features of the transonic nozzle flows including location of sonic line, internal shock and boundary layer build-up. Calculated performance parameters such as thrust, flow rate, and specific impulse are also in reasonble agreement with available data. The results show promising potential of solving full Navier-Stokes equations with heat transfer and two-phase combustion in truly comprehensive modeling of rocket engines.
Przekwas, A. J.↗
Analysis of physical-chemical processes governing SSME internal fluid flows
The basic issues concerning the physical chemical processes of the Space Shuttle Main Engine are discussed. The objectives being to supply the general purpose CFD code PHOENICS and the associated interactive graphics package - GRAFFIC; to demonstrate code usage on SSME related problems; to perform computations and analyses of problems relevant to current and future SSME's; and to participate in the development of new physical models of various processes present in SSME components. These objectives are discussed in detail.
Singhal, A. K.↗
Analysis of physical-chemical processes governing SSME internal fluid flows
The efforts to adapt CHAM's computational fluid dynamics code, PHOENICS, to the analysis of flow within the high pressure fuel turbopump (HPFTP) aft-platform seal cavity of the SSME are summarized. In particular, the special purpose PHOENICS satellite and ground station specifically formulated for this application are listed and described, and the preliminary results of the first part two-dimensional analyses are presented and discussed. Planned three-dimensional analyses are also briefly outlined. To further understand the mixing and combustion processes in the SSME fuelside preburners, a single oxygen-hydrogen jet element was investigated.
Singhal, A. K.↗
SSME thrust chamber simulation using Navier-Stokes equations
The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.
Przekwas, A. J.↗