Engineering topics

Li, C. P.

Publications and source records attributed to Li, C. P..

At least 19 records

NASA NTRS1993

A comparison of high resolution upwind solvers on 3-D inviscid hypersonic flows

A detailed comparison of numerical results obtained by solving the Euler equations for the inviscid flow over delta wings in reentry configurations is presented. The investigation involves a side-by-side comparison between independently developed upwind Euler solvers at VKI (M3D) and NASA Johnson Space Center (E3D) using identical grids. In both solvers the governing equations are integrated by means of time marching finite volume shock capturing methods, based on a cell centered upwind evaluation of the cell face fluxes and nonlinear limiters. High resolution schemes are obtained via MUSCL characteristic variable extrapolation to ensure total variation diminishing (TVD) properties and therefore monotonic discontinuity capturing. Further comparisons are made with several published results by other authors based on both upwind and central discretizations. The present results contribute in asserting the high resolution upwind TVD schemes as the most reliable numerical technique to handle the strong discontinuities typical of high speed flows. The performances of the two upwind solvers are satisfactory and the numerical results in good agreement. However, the important issue of reaching a grid converged solution for the present complex three-dimensional inviscid flow problems was not achieved.

Manna, M.↗

NASA NTRS1993

A multigrid nonoscillatory method for computing high speed flows

A multigrid method using different smoothers has been developed to solve the Euler equations discretized by a nonoscillatory scheme up to fourth order accuracy. The best smoothing property is provided by a five-stage Runge-Kutta technique with optimized coefficients, yet the most efficient smoother is a backward Euler technique in factored and diagonalized form. The singlegrid solution for a hypersonic, viscous conic flow is in excellent agreement with the solution obtained by the third order MUSCL and Roe's method. Mach 8 inviscid flow computations for a complete entry probe have shown that the accuracy is at least as good as the symmetric TVD scheme of Yee and Harten. The implicit multigrid method is four times more efficient than the explicit multigrid technique and 3.5 times faster than the single-grid implicit technique. For a Mach 8.7 inviscid flow over a blunt delta wing at 30 deg incidence, the CPU reduction factor from the three-level multigrid computation is 2.2 on a grid of 37 x 41 x 73 nodes.

Li, C. P.↗

NASA NTRS1993

Computational techniques for high-speed flows with viscous and chemical effects

Algorithms for solving the Euler and the Navier-Stokes equations in conjunction with chemical kinetic equations are presented. The convective flux is estimated from a quasi one dimensional interpolation procedure. Shock, contact, and expansion waves and thermochemical nonequilibrium phenomena are captured by the Lax-Friedrichs technique. Relaxation techniques were developed to enhance their effectiveness in dealing with spatial and temporal stiffness associated with the physical problems. Both explicit and implicit smoothers were implemented into the standard multigrid time stepping method. Unsteady and steady scalar problems are discussed. A perfect gas and equilibrium air shock tube problem is investigated. Numerical schemes and techniques are compared for the problems of shock and boundary layer interaction and three dimensional viscous, nonequilibrium flow encompassing an aerobrake. The results are comparable in accuracy against other high order non-oscillatory techniques. The multigrid methods are assessed using a Mach 8 flow over a complete planar body, a sphere, and a blunt delta wing at 20 deg incidence. Applying an implicit multigrid method on a nested grid of 128 by 64 nodes the reduction factor is 0.25. The central processing unit reduction factor is 2.2 after both the single and multigrid Runge-Kutta solutions converged to machine zero on a grid of 37 by 41 by 73 nodes.

Li, C. P.↗

NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗

NASA NTRS1991

Flowfield computations and comparison with Shuttle aerodynamic data

An in-house developed flow solver, E3D, has been applied to investigate the flow field around the Shuttle Orbiter over an angle-of-attack range of 0 deg to 60 deg at Mach 3.5 and 10.0. The 3D Euler equations are integrated by means of a time-marching finite-volume shock-capturing method, based on cell-centered and upwind evaluation of the cell face fluxes with van Leer or Roe's flux-splitting techniques. A comparison of CFD, wind tunnel, and aeronautical data base is presented. Solutions are in very good agreement with available experimental data. Through the comparisons, it is demonstrated that the E3D code is capable of evaluating total aerodynamic performance.

Ma, E. C.↗

NASA NTRS1991

Numerical simulation of Shuttle ascent transonic flow using an unstructured-grid approach

An unstructured-grid, finite-volume method has been developed for simulating the inviscid flow over spacecrafts of realistic configuration. The grid generation is accomplished by a new technique on the basis of the advancing-front concept. This simple technique is shown to be equally as powerful for a complex multibody as for a single vehicle. Second- or third-order accuracy is obtained via an innovative interpolation procedure similar to the conventional MUSCL approach. This method has been applied to the Shuttle orbiter and a representative Shuttle launch vehicle consisting of the orbiter, the external tank, and the solid rocket boosters. A comparison is discussed between the present results and other results obtained from structured- and unstructured-grid methods.

Wey, T. C.↗

NASA NTRS1991

Computational aspects of chemically reacting flows

An implicit, time-dependent approach is taken to solve the inviscid conservative equations in a vertex-based finite-volume formulation. The convective fluxes are estimated by a modified Riemann solver, and higher-order accuracy is obtained by means of the MUSCL. Implicit techniques are adopted to deal with the equation stiffness as a result of the clustered grid and chemical source term. The real-gas Roe matrix and the ADI and LU schemes have been tested on the sample problems: hypersonic flow over blunt configurations, such as the forebody of a double ellipse, and the near wake of a conic aerobrake. Temperature distribution and standoff distance are validated against the shock-fitting, central-differencing results. Equilibrium temperature may be recovered by artificially increasing the chemical rates at the wall. The present method satisfactorily resolves the flow structure, including the strong expansion behind the flat base. The relative merits and performance of the real-gas TVD and implicit schemes are discussed briefly.

Li, C. P.↗

NASA NTRS1989

Three-dimensional thermochemical nonequilibrium flow modeling for hypersonic flows

A three-dimensional thermochemical nonequilibrium model has been developed and applied to the study of entry flows surrounding space vehicles. The model accounts for both chemical and vibrational nonequilibrium phenomena behind the bow shock. The thermodynamic state of a real gas is modeled with a translational-rotational temperature and a electron-vibrational temperature. Their internal energies are averaged to determine the temperature used in the reaction rates calculation. In order to establish the validity of the selected models, both one- and two-temperature models with seven and/or eleven species were investigated. Several numerical experiments that include a sphere, the RAMC vehicle and 3D AFE forebody flows were performed. Preliminary results were compared with RAMC-II experimental data. Good agreement was obtained after a two-temperature model with eleven species and thirty reactions was incorporated into the study.

Tam, L. T.↗

NASA NTRS1989

Calculation of convective heat transfer on highly blunt bodies at flow incidence

An implicit finite-difference code is used to study three-dimensional viscous heat-conducting flows over the forebody of hypersonic vehicles. In the method, adaptive grids are generated to the shock and body contour, and local flow gradients and total enthalpy are used to control numerical dissipation. Typical axisymmetric configurations of a sphere, ellipsoid, and flat-face disk are considered, along with the cases of a blunt 70-deg cone at 0 and 20 deg flow incidences and an asymmetric 60-deg cone raked off at a 73-deg angle.

Li, C. P.↗

NASA NTRS1989

Computing viscous reactive flow over blunt winged vehicles

A hypersonic flow simulation methodology is developed for analyzing thermochemical nonequilibrium and viscous phenomena associated with space-flight vehicles at realistic conditions. The computer code consists of grid generation, Navier-Stokes and species models, and equation solvers for a complete configuration at arbitrary angles of attack. Provisions are implemented in the code to consider chemically frozen or equilibrium chemistry, inviscid model, conic configurations and simplifications of physical parameters for verifying solution consistency and accuracy. Implications of physical modeling and input parameters on a blunt-cone flow are investigated in details. Applications to USA Orbiter-like vehicles have emphasized the leeside canopy flow as well as the bow-wing shock interaction problem at Mach 22 and 40 deg angle of attack. Comparison of windward centerline pressures and heating rates with available wind-tunnel and flight data is in progress.

Li, C. P.↗

NASA NTRS1988

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

NASA NTRS1988

Comparison of flux-vector and flux-difference splitting techniques for hypersonic flow

Some numerical aspects of flux-vector splitting (FVS) and flux-difference splitting (FDS) schemes are investigated to determine the accuracy for the shock and expansion waves in hypersonic blunt body flow fields. The analysis includes the implementation of Steger-Warming's three-component FVS (Reklis and Thomas, 1981), van Leer's (1986) FVS, and Yang's (1985) FDS to the right-hand side of a time marching scheme where the left-hand side is a variant of the incremental line Gauss-Seidel scheme. Numerous computations are performed by blending two simple flux limiters to test the various schemes on a generic aerobrake at M(infinity) = 10 and zero angle of attack. It is found that the flux limiter plays a very important role in the accuracy of shock capturing. The steady state results are compared to those of a shock-fitting scheme and advantages and disadvantages of each scheme are discussed briefly.

Wey, T. C.↗

NASA NTRS1988

Chemical nonequilibrium and viscous flow computation for conic aerobrake bodies

Three-dimensional analyses are presented for the viscous, reactive flow over a complete entry-body configuration with a wide-angle conic surface. The predictive method uses a split approach that solves iteratively the Navier-Stokes and the continuity equations of chemical species. The finite-difference formulation and the computational grid are adapted to the bow shock and the conformally mapped body such that the velocity components are in the computational spherical-polar space. Combinations of several conic forebody and afterbody configurations have been studied using wind-tunnel, Space Shuttle, and aerobraking orbital transfer vehicle (AOTV) entry conditions. The effects of the borebody bluntness and of finite-rate chemical reactions on the shock layer, the wall catalycity on the boundary layer, the shear-layer impingement on the afterbody, and the base-flow environment are discussed.

Li, C. P.↗

NASA NTRS1987

Chemistry-split techniques for viscous reactive blunt body flow computations

The weak-coupling structure between the fluid and species equations has been exploited and resulted in three, closely related, time-iterative implicit techniques. While the primitive variables are solved in two separated groups and each by an Alternating Direction Implicit (ADI) factorization scheme, the rate-species Jacobian can be treated in either full or diagonal matrix form, or simply ignored. The latter two versions render the split technique to solving for species as scalar rather than vector variables. The solution is completed at the end of each iteration after determining temperature and pressure from the flow density, energy and species concentrations. Numerical experimentation has shown that the split scalar technique, using partial rate Jacobian, yields the best overall stability and consistency. Satisfactory viscous solutions were obtained for an ellipsoidal body of axis ratio 3:1 at Mach 35 and an angle of attack of 20 degrees.

Li, C. P.↗

NASA NTRS1986

Implicit methods for computing chemically reacting flow

The backward Euler scheme was used to solve a large system of inviscid flow and chemical rate equations in three spatial coordinates. The flow equations were integrated simultaneously in time by a conventional ADI factorization technique, then the species equations were solved by either simultaneous or successive techniques. The methods were evaluated in their efficiency and robustness for a hypersonic flow problem involving an aerobrake configuration. It was found that both implicit methods can effectively reduce the stiffness associated with the chemical production term and that the successive solution for the species was as stable as the simultaneous solution. The latter method is more economical because the computation time varies linearly with the number of species.

Li, C. P.↗

NASA NTRS1986

Computation of three-dimensional flow about aerobrake configurations

Ellipsoid, cone and cylinder aerobrake configurations are analyzed to provide comparison data between experimental and model predictions. An analytical model was devised to account for the shock layer ahead of the body and in the near-wake region in terms of the Navier-Stokes equations expressed in conformal polar and azimuthal-angle coordinates. Using polar coordinates simplified the equations by mapping the body onto a sphere, a procedure which also reduced the magnitude of the discretization errors. The equations are then solved using an alternating direction implicit (ADI) factorization technique. Computations were carried out for Mach 3-10 at various grid resolutions and compared with available wind tunnel data. The model generated pressure distributions, heat transfer coefficients and velocity profile data that agreed relatively well with experimental data at a reduced computational cost. Further work is necessary to identify the location of shocks and to model flows about asymmetric configurations.

Li, C. P.↗

NASA NTRS1985

Numerical procedure for three-dimensional hypersonic viscous flow over aerobrake configuration

A numerical method, which is simpler and more efficient than others currently in use, is proposed for the computation of the full viscous flow over an aerobrake body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt forebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined by the body wall, outflow surface and the shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained with a 28 x 36 x 7 grid and 400 time-marching iterations. Excellent agreement of shock location is found between the present result and the schlieren photograph. Details of the base flow and shear layer impingement on the cylindrical aft body are presented for an adiabatic wall case.

Li, C. P.↗

NASA NTRS1985

Computational methods for hypersonic viscous flow over finite ellipsoid-cones at incidence

A numerical method, which is simpler than others currently in use, is proposed for determining the full viscous flow over a finite body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt foebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined to the body wall, outflow surface and the bow shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained on a 30 by 36 by 6 grid with 400 time-marching iterations. Results for body shapes belonging to the ellipsoid-cone family are compared with the experimental data for the Apollo command module and the Viking aeroshell. Validation of the method based on self-consistency is also discussed.

Li, C. P.↗