• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.13 no.3
• /
• pp.119-127
• /
• 1993

Recently, ADI scheme has been a most common tool for solving shallow water equation numerically. But ADI models of tidal flow is likely to cause so called ADI effect in such a region of the Yellow Sea which shows complex topography and has submarine canyons especially. To overcome this, a finite difference algorithm is developed which adopts fully implicit method and preconditioned conjugate gradient squared method. Applying the algorithm including simulation of intertidal zone to Sae-Man-Keum. velocity fields and flooding/drying phenomena are simulated well in spite of complex topography.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.4 no.1
• /
• pp.74-92
• /
• 2000

A numerical analysis has been conducted to study the transient flowfield during the transition from the booster to sustainer phase in an integrated rocket ramjet (IRR) propulsion system. Emphasis is placed on the unsteady inlet aerodynamics, fuel/air mixing in an entire ramjet engine during the flow transient phase. The computational geometry consists of the entire IRR engine, including the inlet, the combustion chamber, and the exhaust nozzle. Turbulence closure is achieved using a low-Reynolds-number two-equation model. The governing equations are solved numerically by means of a finite-volume, preconditioned flux-differencing scheme over a wide range of Mach umber. Various important physical processes are investigated systemically, including terminal shock train.

• The KSFM Journal of Fluid Machinery
• /
• v.11 no.3
• /
• pp.7-13
• /
• 2008

The cavitating flow simulation is of practical importance for many engineering systems, such as pump, turbine, nozzle, Infector, etc. In the present work, a solver for two-phase flows has been developed and applied to simulate the cavitating flows past hydrofoils. The governing equation is the two-phase Navier-Stokes equation, comprised of the continuity equation of liquid and vapor phase. The momentum and energy equation is in the mixture phase. The solver employs an implicit, dual time, preconditioned algorithm using finite difference scheme in curvilinear coordinates. An experimental data and other numerical data were compared with the present results to validate the present solver. It is concluded that the present numerical code has successfully accounted for two-phase Navier-Stokes model of cavitation flow.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2008.03a
• /
• pp.260-265
• /
• 2008

A numerical scheme for solid propellant rocket has been studied using preconditioning method to research unsteady combustion processes for the double-base propellant with a converging-diverging nozzle. The Navier-Stokes equation is solved by dualtime stepping method with finite volume method. The turbulence model uses a shear stress transport modeling. The species equation follows up the method of Xinping WI, Mridul Kumar and Kenneth K. Kuo. A preconditioned algorithm is applied to solve incompressible regime inside the combustor and compressible flow at nozzle. Mass flux was evaluated using modified advective upwind splitting method. The simulated result the comparison a fully coupled implicit method and a semi implicit method in terms of accuracy and efficiency. This report shows the result of solid rocket propellant combustion.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.5 s.176
• /
• pp.1215-1230
• /
• 2000

Parallel Approaches using Cray T3E which is NIPP (Massively Parallel Processors) machine are presented for the efficient computation of the finite element analysis of 3-D shape rolling processes. D omain decomposition method coupled with parallel linear equation solver is used. Domain decomposition is applied for obtaining element tangent stifffiess matrices and residual vectors. Direct and iterative parallel algorithms are used for solving the linear equations. Direct algorithm is_parallel version of direct banded matrix solver. For iterative algorithms, the well-known preconditioned conjugate gradient solver with Jacobi preconditioner is also employed. Moreover a new effective iterative scheme with block inverse matrix preconditioner, which is named by present authors, is presented and its results are compared with the one using Jacobi preconditioner. PVM and MPI are used for message passing and synchronization between processors. The performance and efficiency of each algorithm is discussed and comparisons are made among different algorithms.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.3
• /
• pp.322-329
• /
• 2001

A study of the shock wave turbulent boundary layer interaction is presented. The focus of the study is the interactions of the shock waves with the turbulent boundary layer on the falt plate. Three examples are investigated. The computations are performed, using mixed explicit-implicit generalized Galerkin finite element method. The linear equations at each time step are solved by a preconditioned GMRES algorithm. Numerical results indicate that the implicit scheme converges to the asymptotic steady state much faster than the explicit counterpart. The computed surface pressures and skin friction coefficients display good agreement with experimental data. The flowfield manifests a complex shock wave system and a pair of counter-rotating vortices.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.6 no.3
• /
• pp.216-225
• /
• 1994

By use of laboratory experiment data set for an elliptic shoal by Berkhoff et al. (1982), both accuracy and Performance tests of numerical results between PCGM (Preconditioned Conjugate Gradient Method) and PA(Parabolic Approximation) are compared. Although both results show good agreement with the experimental data the PA model gives better reproduction of the relatively high amplitudes in the section 4-5 downwave of the shoal, in comparison with the PCGM. The PA model has been proved to be a useful tool for predicting wave transformationsin large shallow water region, but it can be applied only to the case of negligible reflection. On the other hand, there is a need to improve the computational efficiency of the PCGM model which is a finite difference scheme directly derived from the mild slope equation and can handle reflection. By taking the results of th PA model as an input data of the PCGM, the CPU time can be reduced by about 40%.