Seokkwan Yoon

Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations

On developpe un schema de relaxation multigrille. Application a des ecoulements transsoniques

An LU-SSOR scheme for the Euler and Navier-Stokes equations

A new multigrid relaxation scheme, lower-upper symmetric successive overrelaxation (LU-SSOR) is developed for the steady-state solution of the Euler and Navier-Stokes equations. The scheme, which is based on central differences, does not require flux splitting for approximate Newton iteration. Application to transonic flow shows that the new method is efficient and robust. The vectorizable LU-SSOR scheme needs only scalar diagonal inversions.

Numerical study of chemically reacting flows using a lower-upper symmetric successive overrelaxation scheme

A new computational fluid dynamics code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite-volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step, finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model

Three-dimensional incompressible Navier-Stokes solver using lower-upper symmetric-Gauss-Seidel algorithm

A numerical method based on the pseudocompressibility concept is developed for solving the three-dimensional incompressible Navier-Stokes equations using the lower-upper symmetric-Gauss-Seidel implicit scheme. Very high efficiency is achieved in a new flow solver, INS3D-LU code, by accomplishing the complete vectorizability of the algorithm on obIique planes of sweep in three dimensions.

Fully coupled implicit method for thermochemical nonequilibrium air at suborbital flight speeds

A CFD technique is described in which the finite-rate chemistry in thermal and chemical nonequilibrium air is fully and implicitly coupled with the fluid motion. Developed for use in the suborbital hypersonic flight speed range, the method accounts for nonequilibrium vibrational and electronic excitation and dissociation, but not ionization. The steady-state solution to the resulting system of equations is obtained by using a lower-upper factorization and symmetric Gauss-Seidel sweeping technique through Newton iteration. Inversion of the left-hand-side matrices is replaced by scalar multiplications through the use of the diagonal dominance algorithm. The code, named CENS2H (Compressible-Euler-Navier-Stokes Two-Dimensional Hypersonic), is fully vectorized and requires about 8.8 x 10 to the -5th sec per node point per iteration using a Cray X-MP computer. Converged solutions are obtained after about 2400 iterations. Sample calculations are made for a circular cylinder and a 10 percent airfoil at 5 deg angle of attack. The calculated cylinder flow field agrees with that obtained experimentally. The code predicts a 10 percent change in lift, drag, and pitching moment for the airfoil due to the thermochemical phenomena.

Calculation of real-gas effects on blunt-body trim angles

The effect of vibrational excitation and dissociation at high temperatures on the trim angle of attack of a blunt lifting body is calculated for a nonequilibrium flow regime in air using a CFD technique. The vibrational-electronic temperature and the species densities are calculated assuming the flow to be in a nonequilibrium state. The forebody flow of a two-dimensional blunt body of the shape of the Apollo Command Module at a finite angle of attack is calculated. The results show that the pitching moment around a reference point is larger and the trim angle of attack is smaller for a reacting gas than for a perfect gas. The calculated shift in the trim angle due to the real-gas effect is of the same order as that seen during the Apollo flights.

Artificial dissipation models for hypersonic external flow

Four artificial dissipation models which augment central difference schemes were examined for hypersonic external flows. The models were a first and third order dissipation model, a directionally scaled first and third order dissipation model, a flux limited dissipation model, and a flux difference split dissipation model. Each model was implemented in the lower-upper symmetric-Gauss-Seidel (LU-SGS) algorithm to solve the full Navier-Stokes equations. The latter two models can be regarded as total variation diminishing (TVD) schemes. Test results for model problems showed that the flux limited dissipation model was robust enough to predict a high speed blunt body flow with strong shock and expansion waves. The flux difference split dissipation model was capable of shock capturing with higher resolution, but was less robust. First and third order dissipation models turned out to be neither accurate nor robust enough for high Mach number flow computations.

Numerical study of chemically reacting flows using an LU scheme

A new computational fluid dynamic code has been developed for the study of mixing and chemical reactions in the flow fields of ramjets and scramjets. The code employs an implicit finite volume, lower-upper symmetric successive overrelaxation scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a fully-coupled and very efficient manner. The combustion processes are modeled by an 8-species, 14-step finite rate chemistry model whereas turbulence is simulated by a Baldwin-Lomax algebraic model. The validity of the code is demonstrated by comparing the numerical calculations with both experimental data and previous calculations of a cold flow helium injection into a straight channel and premixed hydrogen-air reacting flows in a ramped duct. The code is then used to calculate the mixing and chemical reactions of a hydrogen jet transversely injected into a supersonic airstream. Results are presented describing the flow field, the recirculation regions in front and behind the injector, and the chemical reactions.

LU-SGS implicit algorithm for three-dimensional incompressible Navier-Stokes equations with source term

A numerical method is developed for solving the incompressible Navier-Stokes equations using the concept of pseudocompressibility. A lower-upper symmetric-Gauss-Seidel implicit scheme is developed for three-dimensional incompressible viscous flow computations. The present algorithm offers additional advantages when solving the flow equations with source terms. Complete vectorizability of the algorithm on oblique planes of sweep in three-dimensions is accomplished in a new flow solver, INS3D-LU code. Spatial differencing is a second-order accurate semi-discrete finite-volume method augmented by a third-order accurate numerical dissipation model which is based on spectral-radii. Comparison of numerical solutions for a curved duct with experimental data shows good agreement. The method is applied to calculate the inducer flow of the Space Shuttle Main Engine turbopump.

Multigrid convergence of an implicit symmetric relaxation scheme

The multigrid method has been applied to an existing three-dimensional compressible Euler solver to accelerate the convergence of the implicit symmetric relaxation scheme. This lower-upper symmetric Gauss-Seidel implicit scheme is shown to be an effective multigrid driver in three dimensions. A grid refinement study is performed including the effects of large cell aspect ratio meshes. Performance figures of the present multigrid code on Cray computers including the new C90 are presented. A reduction of three orders of magnitude in the residual for a three-dimensional transonic inviscid flow using 920 k grid points is obtained in less than 4 min on a Cray C90. 23 refs.