Engineering topics

Liou, Meng-Sing

Publications and source records attributed to Liou, Meng-Sing.

105 records · Page 6

NASA NTRS1989

Flux splitting algorithms for two-dimensional viscous flows with finite-rate chemistry

The Roe flux-difference splitting method has been extended to treat two-dimensional viscous flows with nonequilibrium chemistry. The derivations have avoided unnecessary assumptions or approximations. For spatial discretization, the second-order Roe upwind differencing is used for the convective terms and central differencing for the viscous terms. An upwind-based TVD scheme is applied to eliminate oscillations and obtain a sharp representation of discontinuities. A two-stage Runge-Kutta method is used to time integrate the discretized Navier-Stokes and species transport equations for the asymptotic steady solutions. The present method is then applied to two types of flows: the shock wave/boundary layer interaction problems and the jet in cross flows.

Shuen, Jian-Shun↗

NASA NTRS1989

Conservative treatment of boundary interfaces for overlaid grids and multi-level grid adaptations

Conservative algorithms for boundaray interfaces of overlaid grids are presented. The basic method is zeroth order, and is extended to a higher order method using interpolation and subcell decomposition. The present method, strictly based on a conservative constraint, is tested with overlaid grids for various applications of unsteady and steady supersonic inviscid flows with strong shock waves. The algorithm is also applied to a multi-level grid adaptation in which the next level finer grid is overlaid on the coarse base grid with an arbitrary orientation.

Moon, Young J.↗

NASA NTRS1989

A time accurate finite volume high resolution scheme for three dimensional Navier-Stokes equations

A time accurate, three-dimensional, finite volume, high resolution scheme for solving the compressible full Navier-Stokes equations is presented. The present derivation is based on the upwind split formulas, specifically with the application of Roe's (1981) flux difference splitting. A high-order accurate (up to the third order) upwind interpolation formula for the inviscid terms is derived to account for nonuniform meshes. For the viscous terms, discretizations consistent with the finite volume concept are described. A variant of second-order time accurate method is proposed that utilizes identical procedures in both the predictor and corrector steps. Avoiding the definition of midpoint gives a consistent and easy procedure, in the framework of finite volume discretization, for treating viscous transport terms in the curvilinear coordinates. For the boundary cells, a new treatment is introduced that not only avoids the use of 'ghost cells' and the associated problems, but also satisfies the tangency conditions exactly and allows easy definition of viscous transport terms at the first interface next to the boundary cells. Numerical tests of steady and unsteady high speed flows show that the present scheme gives accurate solutions.

Liou, Meng-Sing↗

NASA NTRS1989

A detailed analysis of inviscid flux splitting algorithms for real gases with equilibrium or finite-rate chemistry

The extension of the known flux-vector and flux-difference splittings to real gases via rigorous mathematical procedures is demonstrated. Formulations of both equilibrium and finite-rate chemistry for real-gas flows are described, with emphasis on derivations of finite-rate chemistry. Split-flux formulas from other authors are examined. A second-order upwind-based TVD scheme is adopted to eliminate oscillations and to obtain a sharp representation of discontinuities.

Shuen, Jian-Shun↗

NASA NTRS1988

A numerical study of ENO and TVD schemes for shock capturing

The numerical performance of a second-order upwind-based total variation diminishing (TVD) scheme and that of a uniform second-order essentially non-oscillatory (ENO) scheme for shock capturing are compared. The TVD scheme used is a modified version of Liou, using the flux-difference splitting (FDS) of Roe and his superbee function as the limiter. The construction of the basic ENO scheme is based on Harten, Engquist, Osher, and Chakravarthy, and the 2-D extensions are obtained by using a Strang-type of fractional-step time-splitting method. Numerical results presented include both steady and unsteady, 1-D and 2-D calculations. All the chosen test problems have exact solutions so that numerical performance can be measured by comparing the computer results to them. For 1-D calculations, the standard shock-tube problems of Sod and Lax are chosen. A very strong shock-tube problem, with the initial density ratio of 400 to 1 and pressure ratio of 500 to 1, is also used to study the behavior of the two schemes. For 2-D calculations, the shock wave reflection problems are adopted for testing. The cases presented in this report include flows with Mach numbers of 2.9, 5.0, and 10.0.

Chang, Shih-Hung↗

NASA NTRS1988

Solutions of one-dimensional steady nozzle flow revisited

A formulation is presented by which any iteration process for obtaining the entropy increase in the flow of a one-dimensional steady nozzle is eliminated, and the simple solution of a quadratic equation is obtained. The proper parameters are then explicitly seen in the equation, and their effects on the solution are easily determined. Since only one root of the equation is physically admissible, entropy production, and therefore the shock wave, are uniquely determined by this set of parameters.

Liou, Meng-Sing↗

NASA NTRS1988

Splitting of inviscid fluxes for real gases

Flux-vector and flux-difference splittings for the inviscid terms of the compressible flow equations are derived under the assumption of a general equation of state for a real gas in equilibrium. No necessary assumptions, approximations or auxiliary quantities are introduced. The formulas derived include several particular cases known for ideal gases and readily apply to curvilinear coordinates. Applications of the formulas in a TVD algorithm to one-dimensional shock-tube and nozzle problems show their quality and robustness.

Liou, Meng-Sing↗

NASA NTRS1988

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, van Leer and Roe are tested in all possible combinations on the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust; it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗

NASA NTRS1988

A computational analysis of under-expanded jets in the hypersonic regime

Underexpanded axisymmetric jets are studied numerically using a full Navier-Stokes solver. Emphasis has been given to supersonic and hypersonic jets in supersonic and hypersonic ambient flows, a phenomenon previously overlooked. It is demonstrated that the shear layers and shock patterns in a jet plume can be captured without complicated viscous/inviscid and subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief effect has been identified for underexpanded jets in supersonic ambient flows. While it is well known that an underexpanded jet in a quiescent ambience (or subsonic ambience) contains multiple shock cells, the present study shows that because of the supersonic pressure relief effect, an underexpanded jet in a supersonic or hypersonic ambience contains only one major shock cell.

Hsu, Andrew T.↗

NASA NTRS1988

Choice of implicit and explicit operators for the upwind differencing method

The flux-vector and flux-difference splittings of Steger-Warming, Van Leer and Roe are tested in all possible combinations in the implicit and explicit operators that can be distinguished in implicit relaxation methods for the steady Euler and Navier-Stokes equations. The tests include one-dimensional inviscid nozzle flow, and two-dimensional inviscid and viscous shock reflection. Roe's splitting, as anticipated, is found to uniformly yield the most accurate results. On the other hand, an approximate Roe splitting of the implicit operator (the complete Roe splitting is too complicated for practical use) proves to be the least robust with regard to convergence to the steady state. In this respect, the Steger-Warming splitting is the most robust: it leads to convergence when combined with any of the splittings in the explicit operator, although not necessarily in the most efficient way.

Liou, Meng-Sing↗

NASA NTRS1988

Splitting of inviscid fluxes for real gases

Liou, Meng-Sing↗

NASA NTRS1988

A comparison of ENO and TVD schemes

The numerical performance of a second-order upwind-based TVD scheme is compared with that of a uniform second-order ENO scheme on shock capturing. The cases considered include flows with Mach numbers of 2.9, 5.0, and 10.0. For cases with Mach numbers of 5.0 and 10.0, the computed ENO results are inferior to the corresponding TVD results.

Chang, Shih-Hung↗

NASA NTRS1988

A computational analysis of under-expanded jets in the hypersonic regime

Underexpanded axisymmetric jets are studied numerically using a full Navier-Stokes solver. Emphasis has been given to supersonic and hypersonic jets in supersonic and hypersonic ambient flows, a phehomenon previously being overlooked for the most part. The present work demonstrates that the shear layers and shock patterns in a jet plume can be captured without complicated viscous/inviscid and subsonic/supersonic coupling schemes. In addition, a supersonic pressure relief effect has been identified for underexpanded jets in supersonic ambient flows. While it is well known that an underexpanded jet in a quiescent ambience (or subsonic ambience) contains multiple shock cells, the present study shows that because of the supersonic pressure relief effect, an underexpanded jet in a supersonic or hypersonic ambience contains only one major shock cell.

Hsu, Andrew T.↗

NASA NTRS1987

A generalized procedure for constructing an upwind based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some linear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to scalar nonlinear equations and a constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance in numerical experiments. The formulation is extended to the nonlinear system. It is demonstrated that the present procedure can easily convert existing central or upwind, and second- or higher-order differencing schemes to preserve monotonicity and yield physically admissible solutions. The formulation is simple mathematically as well as numerically; both matrix-vector multiplication and Riemann solver are avoided. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made.

Liou, Meng-Sing↗

NASA NTRS1987

A generalized procedure for constructing an upwind-based TVD scheme

A generalized formulation for constructing second- and higher-order accurate TVD (total variation diminishing) schemes is presented. A given scheme is made TVD by limiting antidiffusive flux differences with some nonlinear functions, so-called limiters. The general idea of the formulation and its mathematical proof of Harten's TVD conditions is shown by applying the Lax-Wendroff method to a scalar nonlinear equation and constant-coefficient system of conservation laws. For the system of equations, several definitions are derived for the argument used in the limiter function and present their performance to numerical experiments. Then the formulation is formally extended to the nonlinear system of equations. It is demonstrated that use of the present procedure allows easy conversion of existing central or upwind, and second- or higher-order differencing schemes so as to preserve monotonicity and to yield physically admissible solutions. The formulation is simple mathematically as well as numerically; neither matrix-vector multiplication nor Riemann solver is required. Roughly twice as much computational effort is needed as compared to conventional scheme. Although the notion of TVD is based on the initial value problem, application to the steady Euler equations of the formulation is also made. Numerical examples including various ranges of problems show both time- and spatial-accuracy in comparison with exact solutions.

Liou, Meng-Sing↗