Engineering topics
Publications and source records attributed to Venkateswaran, S..
An implementation of both finite rate and equilibrium chemistry have been completed for the OVERFLOW code, a chimera capable, complex geometry flow code widely used to predict transonic flow fields. The implementation builds on the computational efficiency and geometric generality of the solver.
Computational fluid dynamic analyses of the Navier-Stokes equations coupled with solid-phase pyrolysis, gas-phase combustion, turbulence and radiation are performed to study hybrid rocket combustor flowfields. The computational study is closely co-ordinated with a companion experimental program using a planar slab burner configuration with HTPB as fuel and gaseous oxygen. Computational predictions agree reasonably well with measurement data of fuel regression rates and surface temperatures. Additionally, most of the parametric trends predicted by the model are in general agreement with experimental trends. The computational model is applied to extend the results from the lab-scale to a full-scale axisymmetric configuration. The numerical predictions indicate that the full-scale configuration burns at a slower rate than the lab-scale combustor under identical specific flow rate conditions. The results demonstrate that detailed CFD analyses can play a useful role in the design of hybrid combustors.
The effects of high-aspect-ratio grids on algorithm convergence are considered by means of vector stability theory and computational experiments. The results indicate that approximately factored implicit schemes experience convergence deterioration because of nonoptimum local time-stepping procedures and increased need for viscous preconditioning. Based on this insight, an enhanced algorithm is devised using improved selection of the local time step, appropriate definition of the viscous preconditioning matrix, and proper implementation of the boundary conditions. The new algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.
The objective of the study is to assess the capability of two-equation turbulence models for simulating propulsion-related flowfields. The standard kappa-epsilon model with Chien's low Reynolds number formulation for near-wall effects is used as the baseline turbulence model. Several experimental test cases, representative of rocket combustor internal flowfields, are used to catalog the performance of the baseline model. Specific flowfields considered here include recirculating flow behind a backstep, mixing between coaxial jets and planar shear layers. Since turbulence solutions are notoriously dependent on grid and numerical methodology, the effects of grid refinement and artificial dissipation on numerical accuracy are studied. In the latter instance, computational results obtained with several central-differenced and upwind-based formulations are compared. Based on these results, improved turbulence modes such as enhanced kappa-epsilon models as well as other two-equation formulations (e.g., kappa-omega) are being studied. In addition, validation of swirling and reacting flowfields are also currently underway.
Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 seconds are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies-done by scaling the frequency from 2.45 GHz to 0.91 GHz-indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on axis.
The performance of Navier-Stokes codes are influenced by several phenomena. For example, the robustness of the code may be compromised by the lack of grid resolution, by a need for more precise initial conditions or because all or part of the flowfield lies outside the flow regime in which the algorithm converges efficiently. A primary example of the latter effect is the presence of extended low Mach number and/or low Reynolds number regions which cause convergence deterioration of time marching algorithms. Recent research into this problem by several workers including the present authors has largely negated this difficulty through the introduction of time-derivative preconditioning. In the present paper, we employ the preconditioned algorithm to address convergence difficulties arising from sensitivity to grid stretching and high aspect ratio grid cells. Strong grid stretching is particularly characteristic of turbulent flow calculations where the grid must be refined very tightly in the dimension normal to the wall, without a similar refinement in the tangential direction. High aspect ratio grid cells also arise in problems that involve high aspect ratio domains such as combustor coolant channels. In both situations, the high aspect ratio cells can lead to extreme deterioration in convergence. It is the purpose of the present paper to address the reasons for this adverse response to grid stretching and to suggest methods for enhancing convergence under such circumstances. Numerical algorithms typically possess a maximum allowable or optimum value for the time step size, expressed in non-dimensional terms as a CFL number or vonNeumann number (VNN). In the presence of high aspect ratio cells, the smallest dimension of the grid cell controls the time step size causing it to be extremely small, which in turn results in the deterioration of convergence behavior. For explicit schemes, this time step limitation cannot be exceeded without violating stability restrictions of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios.
Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 sec are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies (done by scaling the frequency from 2.45 GHz to 0.91 GHz) indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on-axis.
Vector stability theory and computational experiments are used to examine the effects of high aspect ratio grids on algorithm convergence. Emphasis is placed on identifying the degree of convergence deterioration that ADI schemes experience as grid aspect ratio increases and on deriving methods for circumventing this problem. Approximately factored implicit schemes are found to experience convergence deterioration because of nonoptimum local time stepping procedures and the increased need for viscous preconditioning. An enhanced algorithm is devised on the basis of improved selection of the local time step, appropriate definition of viscous preconditioning matrix, and proper implementation of boundary conditions. The novel algorithm provides uniformly efficient convergence at all aspect ratios for both Euler and Navier-Stokes computations for a variety of test problems.
The Langley 8 foot high temperature tunnel (8 ft HTT) is used to test components of hypersonic vehicles for aerothermal loads definition and structural component verification. The test medium of the 8 ft HTT is obtained by burning a mixture of methane and air under high pressure; the combustion products are expanded through an axisymmetric conical contoured nozzle to simulate atmospheric flight at Mach 7. This facility was modified to raise the oxygen content of the test medium to match that of air and to include Mach 4 and Mach 5 capabilities. These modifications will facilitate the testing of hypersonic air breathing propulsion systems for a wide range of flight conditions. A computational method to predict the thermodynamic, transport, and flow properties of the equilibrium chemically reacting oxygen enriched methane-air combustion products was implemented in a computer code. This code calculates the fuel, air, and oxygen mass flow rates and test section flow properties for Mach 7, 5, and 4 nozzle configurations for given combustor and mixer conditions. Salient features of the 8 ft HTT are described, and some of the predicted tunnel operational characteristics are presented in the carpet plots to assist users in preparing test plans.
We have implemented preconditioning for multi-species reacting flows in two independent codes, an implicit (ADI) code developed in-house and the RPLUS code (developed at LeRC). The RPLUS code was modified to work on a four-stage Runge-Kutta scheme. The performance of both the codes was tested, and it was shown that preconditioning can improve convergence by a factor of two to a hundred depending on the problem. Our efforts are currently focused on evaluating the effect of chemical sources and on assessing how preconditioning may be applied to improve convergence and robustness in the calculation of reacting flows.
Two widely used family of algorithms, pressure-based and density-based methods, have been developed for computational fluid dynamics (CFD) problems over the years. Pressure-based methods (such as SIMPLE and PISO) use a Poisson-like equation for updating pressure instead of the continuity equation, while density-based methods use the continuity equation to update density (an equation of state is used to provide density in pressure based schemes and pressure in density based schemes). Pressure-based methods were developed originally for incompressible flows at low Reynolds numbers and were then extended to high Reynolds numbers and compressible applications. On the other hand, density based methods were originally developed for transonic flows and have been extended down to low Mach numbers through the use of preconditioning techniques. We compare these two very different approaches to solving the Navier-Stokes equations in order to gain an understanding of their similarities and differences. Specifically, we consider the PISO scheme as a representative pressure-based method and contrast it with a recently developed preconditioning scheme. We also compare the relative performance of the PISO algorithm with a Euler implicit algorithm that is employed to solve the preconditioned equations by means of a vector stability analysis.
A previous time-derivative preconditioning procedure for solving the Navier-Stokes is extended to the chemical species equations. The scheme is implemented using both the implicit ADI and the explicit Runge-Kutta algorithms. A new definition for time-step is proposed to enable grid-independent convergence. Several examples of both reacting and non-reacting propulsion-related flowfields are considered. In all cases, convergence that is superior to conventional methods is demonstrated. Accuracy is verified using the example of a backward facing step. These results demonstrate that preconditioning can enhance the capability of density-based methods over a wide range of Mach and Reynolds numbers.
Explore the source record for details and available documents.
Combustion instability in liquid rocket engines is simulated computationally by using a simple two-parameter model for the combustion response function. The objectives of the study are to assess the capabilities of CFD algorithms for instability studies and to investigate the response to parametric effects such as bombs and distributed combustion. Results indicate that numerical solutions of high accuracy can be obtained if a sufficient number of grid points are used per wavelength of the disturbance. The short-term response to bombs or pulses triggers a large number of modes in the combustor whose faithful resolution requires highly dense grids, although there is evidence that correct long-term solutions can be obtained even if all the short-term frequencies are not resolved. Long-term responses to pulses are shown to decay to the most unstable mode in small amplitude cases, and to exhibit limit cycles in large amplitude cases. Comparison of distributed with concentrated heat release indicates the former is more stable for given values of the combustion response parameters, and that the distributed heat release gives rise to higher frequency disturbances. Wave steepening is observed in the solutions, but its effect is less pronounced in multidimensional waves than in one-dimensional waves.
An experimental study of an axial-compression interior-corner region with shock impingement was conducted in the NASA Langley 20-inch Mach 6 Tunnel to obtain aerodynamic heating rate distributions. The model, a generic hypersonic engine inlet, had a sharp leading edge forebody, cowl, and interchangeable struts with sharp and blunt leading edges. Experiments were conducted at free-stream unit Reynolds numbers of 1.68 x 10 to the 6th and 3.35 x 10 to the 6th per foot for the blunt and sharp-strut configurations, respectively. For the sharp-strut configuration, with and without trips, the peak heating occurred on the cowl, near the compression corner, downstream of the forebody oblique shock impingement region. The peak heating levels, for both the tripped and untripped boundary-layer flows were about 45 and 50 times the reference laminar heating levels, respectively. The peak heating for the blunt leading edge strut model, which was 30 times the reference value, occurred on the cowl, away from the corner, where the bow shock from the strut leading edge intersected the cowl surface.