Engineering topics

Merkle, Charles L.

Publications and source records attributed to Merkle, Charles L..

At least 37 records · Page 2

NASA NTRS1993

CFD analyses of coolant channel flowfields

The flowfield characteristics in rocket engine coolant channels are analyzed by means of a numerical model. The channels are characterized by large length to diameter ratios, high Reynolds numbers, and asymmetrical heating. At representative flow conditions, the channel length is approximately twice the hydraulic entrance length so that fully developed conditions would be reached for a constant property fluid. For the supercritical hydrogen that is used as the coolant, the strong property variations create significant secondary flows in the cross-plane which have a major influence on the flow and the resulting heat transfer. Comparison of constant and variable property solutions show substantial differences. In addition, the property variations prevent fully developed flow. The density variation accelerates the fluid in the channels increasing the pressure drop without an accompanying increase in heat flux. Analyses of the inlet configuration suggest that side entry from a manifold can affect the development of the velocity profile because of vortices generated as the flow enters the channel. Current work is focused on studying the effects of channel bifurcation on the flow field and the heat transfer characteristics.

Yagley, Jennifer A.↗

NASA NTRS1993

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids.

Buelow, Philip E. O.↗

NASA NTRS1993

Computational fluid dynamic modelling of cavitation

Models in sheet cavitation in cryogenic fluids are developed for use in Euler and Navier-Stokes codes. The models are based upon earlier potential-flow models but enable the cavity inception point, length, and shape to be determined as part of the computation. In the present paper, numerical solutions are compared with experimental measurements for both pressure distribution and cavity length. Comparisons between models are also presented. The CFD model provides a relatively simple modification to an existing code to enable cavitation performance predictions to be included. The analysis also has the added ability of incorporating thermodynamic effects of cryogenic fluids into the analysis. Extensions of the current two-dimensional steady state analysis to three-dimensions and/or time-dependent flows are, in principle, straightforward although geometrical issues become more complicated. Linearized models, however offer promise of providing effective cavitation modeling in three-dimensions. This analysis presents good potential for improved understanding of many phenomena associated with cavity flows.

Deshpande, Manish↗

NASA NTRS1993

CFD modeling of microwave electrothermal thrusters

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.

Schwer, Douglas A.↗

NASA NTRS1993

Cavitation Modeling in Euler and Navier-Stokes Codes

Many previous researchers have modeled sheet cavitation by means of a constant pressure solution in the cavity region coupled with a velocity potential formulation for the outer flow. The present paper discusses the issues involved in extending these cavitation models to Euler or Navier-Stokes codes. The approach taken is to start from a velocity potential model to ensure our results are compatible with those of previous researchers and available experimental data, and then to implement this model in both Euler and Navier-Stokes codes. The model is then augmented in the Navier-Stokes code by the inclusion of the energy equation which allows the effect of subcooling in the vicinity of the cavity interface to be modeled to take into account the experimentally observed reduction in cavity pressures that occurs in cryogenic fluids such as liquid hydrogen. Although our goal is to assess the practicality of implementing these cavitation models in existing three-dimensional, turbomachinery codes, the emphasis in the present paper will center on two-dimensional computations, most specifically isolated airfoils and cascades. Comparisons between velocity potential, Euler and Navier-Stokes implementations indicate they all produce consistent predictions. Comparisons with experimental results also indicate that the predictions are qualitatively correct and give a reasonable first estimate of sheet cavitation effects in both cryogenic and non-cryogenic fluids. The impact on CPU time and the code modifications required suggests that these models are appropriate for incorporation in current generation turbomachinery codes.

Deshpande, Manish↗

NASA NTRS1993

CFD analyses of coolant channel flowfields

The flowfield characteristics in a rocket engine coolant channels are analyzed by means of a numerical model. The channels are characterized by large length to diameter ratios, high Reynolds numbers, and asymmetrical heating. At representative flow conditions, the channel length is approximately twice the hydraulic entrance length so fully developed conditions are reached. The supercritical hydrogen coolant introduces strong property variations that have a major influence on the developing flow and the resulting heat transfer. Comparisons of constant and variable property solutions show substantial differences. The density variation accelerates the fluid in the channels increasing the pressure drop without an accompanying increase in heat flux. Analyses of the inlet configuration suggest that side entry from a manifold can affect the development of the velocity profile because of vortices generated as the flow enters the channel.

Yagley, J. A.↗

NASA NTRS1993

Analysis of microwave-heated rocket engines for space propulsion

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.

Schwer, D. A.↗

NASA NTRS1993

Prediction of engine and near-field plume reacting flows in low-thrust chemical rockets

A computational model is employed to study the reacting flow within the engine and near-field plumes of several small gaseous hydrogen-oxygen thrusters. The model solves the full Navier-Stokes equations coupled with species diffusion equations for a hydrogen-oxygen reaction kinetics system and includes a two-equation q-omega model for turbulence. Predictions of global performance parameters and localized flowfield variables are compared with experimental data in order to assess the accuracy with which these flowfields are modeled and to identify aspects of the model which require improvement. Predicted axial and radial velocities 3 mm downstream of the exit plane show reasonable agreement with the measurements. The predicted peak in axial velocity in the hydrogen film coolant along the nozzle wall shows the best agreement; however, predictions within the core region are roughly 15 percent below measured values, indicating an underprediction of the extent to which the hydrogen diffuses and mixes with the core flow. There is evidence that this is due to three-dimensional mixing processes which are not included in the axisymmetric model.

Weiss, Jonathan M.↗

NASA NTRS1992

Comparison between the PISO algorithm and preconditioning methods for compressible flow

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.

Merkle, Charles L.↗

NASA NTRS1992

Multi-dimensional analysis of combustion instabilities in liquid rocket motors

A three-dimensional analysis of combustion instabilities in liquid rocket engines is presented based on a mixed finite difference/spectral solution methodology for the gas phase and a discrete droplet tracking formulation for the liquid phase. Vaporization is treated by a simplified model based on an infinite thermal conductivitiy assumption for spherical liquid droplets of fuel in a convective environment undergoing transient heating. A simple two parameter phenomenological combustion response model is employed for validation of the results in the small amplitude regime. The computational procedure is demonstrated to capture the phenomena of wave propagation within the combustion chamber accurately. Results demonstrate excellent amplitude and phase agreement with analytical solutions for properly selected grid resolutions under both stable and unstable operating conditions. Computations utilizing the simplified droplet model demonstrate stable response to arbitrary pulsing. This is possibly due to the assumption of uniform droplet temperature which removes the thermal inertia time-lag response of the vaporization process. The mixed-character scheme is sufficiently efficient to allow solutions on workstations at a modest increase in computational time over that required for two-dimensional solutions.

Grenda, Jeffrey M.↗

NASA NTRS1992

The relationship between pressure-based and density-based algorithms

The PISO, pressure-based algorithm, is compared with implicit time-marching systems to ascertain their similarities and differences. Both methods are expressed in vector form for comparison purposes. The vector form of the PISO method is triangular, allowing an uncoupled solution procedure, while the Euler implicit method requires the simultaneous solution of all equations. Upwind differencing is performed according to the direction of the eigenvalues in both systems, but this is the particle velocity in the PISO method and the acoustic velocity in the Euler method. Vector stability calculations show that the PISO method is conditionally stable depending on the time step, but that it provides adequate damping at small time steps to enable good convergence. Unconditional stability can be provided by retaining the energy coupling terms in the time derivative. Transforming the Euler implicit equations to the PISO variables along with a modification of the time derivatives gives the density-based method the same amplification factors at low speeds as the pressure-based method without affecting their behavior at high speeds.

Merkle, Charles L.↗

NASA NTRS1991

Time-derivative preconditioning for viscous flows

A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.

Choi, Yunho↗

NASA NTRS1991

Time-derivative preconditioning for viscous flows

Choi, Yunho↗

NASA NTRS1991

A review of liquid rocket propulsion programs in Japan

An assessment of Japan's current capabilities in the areas of space and transatmospheric propulsion is presented. The primary focus is upon Japan's programs in liquid rocket propulsion and in space plane and related transatmospheric areas. Brief reference is also made to their solid rocket programs, as well as to their supersonic air breathing propulsion efforts that are just getting underway.

Merkle, Charles L.↗

NASA NTRS1991

An overview of the Penn State Propulsion Engineering Research Center

An overview of the Penn State Propulsion Engineering Research Center is presented. The following subject areas are covered: research objectives and long term perspective of the Center; current status and operational philosophy; and brief description of Center projects (combustion, fluid mechanics and heat transfer, materials compatibility, turbomachinery, and advanced propulsion concepts).

Merkle, Charles L.↗

NASA NTRS1991

CFD applications in chemical propulsion engines

The present research is aimed at developing analytical procedures for predicting the performance and stability characteristics of chemical propulsion engines. Specific emphasis is being placed on understanding the physical and chemical processes in the small engines that are used for applications such as spacecraft attitude control and drag make-up. The small thrust sizes of these engines lead to low nozzle Reynolds numbers with thick boundary layers which may even meet at the nozzle centerline. For this reason, the classical high Reynolds number procedures that are commonly used in the industry are inaccurate and of questionable utility for design. A complete analysis capability for the combined viscous and inviscid regions as well as for the subsonic, transonic, and supersonic portions of the flowfield is necessary to estimate performance levels and to enable tradeoff studies during design procedures.

Merkle, Charles L.↗

NASA NTRS1991

CFD Applications in Propulsion

An overview of various applications of Computational Fluid Dynamics (CFD) algorithms to propulsion problems is given. Topics of interest include incompressible, low speed compressible, transonic, and supersonic problems. A common family of algorithms is used for all applications and emphasis is placed on maintaining accuracy and convergence efficiency for all problems. Specific problems include pump hydrodynamics, simultaneous combustion and mixing in rocket engines, viscous nozzle flow, and CFD applications to combustion stability.

Merkle, Charles L.↗

NASA NTRS1991

Solar rocket plume/mirror interactions

The extent to which the plume from a solar thermal rocket will impinge on the solar collector is studied by flow field analysis. Such interaction can adversely affect collector performance through fouling, excessive heat loading, or pressure loads that deform the delicate structures. The geometrical shape of the collector is such that only the flow from the nozzle boundary layer can reach it, but the thrust levels of interest lead to very viscous nozzle flows with thick boundary layers. Reasonable accuracy in solving these flows requires a fully coupled viscous-inviscid procedure. Results show that the fraction of the plume that hits the collector can be well estimated by continuum theory, but that transitional and rarefied phenomena will have some impact on how it is distributed over the surface. Initial results for one representative condition show that approx. 4 percent of the total flow in the jet makes its way to the collector. The pressures on the collector, however, remain quite low because of its distance from the engine. Additional work is needed to document the effect of thrust scaling and wall cooling on impingement.

Yu, Sheng-Tao↗