• 한국전산유체공학회:학술대회논문집
• /
• 1998.11a
• /
• pp.90-95
• /
• 1998

The two-dimensional incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. Incompressible code using pseudo-compressibility and dual-time stepping method involves a conventional upwind differencing scheme for the convective terms and LU-SGS scheme for time integration. Compressible code also adopts an FDS scheme and LU-SGS scheme. Several two-equation turbulence models (the standard $k-{\varepsilon}$ model, the $k-{\omega}$ model. and $k-{\omega}$ SST model) are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by computing the flow around the transonic RAE2822 airfoil and the NACA4412 airfoil, respectively. Both the results show a good agreement with experimental surface pressure coefficients and velocity profiles in the boundary layers. Also, the GA(W)-1 single airfoil and the NLR7301 airfoil with a flap are computed using the two-equation turbulence models. The grid systems around two- and three-element airfoil are efficiently generated using Chimera grid scheme, one of the overlapping grid generation methods.

• Journal of the Korea Institute of Building Construction
• /
• v.19 no.1
• /
• pp.39-46
• /
• 2019

A three-dimensional flow simulation model for high flow concrete was developed using Incompressible Smoothed Particle Hydrodynamics (ISPH), which can solved Navier-Stokes equation with the assumption of a fluid to be incompressible. For the simulation, a computer program code for ISPH was implemented with MATALB programming code. A piecewise cubic spline function was used for the kernel function of ISPH. Projetion method was used to calculate the velocity and pressure of particles as a function of time. Fixed ghost particle was used for wall boundary condition. Free surface boundaries were determined by using virtual density of particles. In order to validate the model and the code, the simulation results of slump flow test, $T_{500}$ test and L-box test were compared with experimental ones. The simulation results were well matched with the experimental results. The simulation described successfully the characteristics of the flow phenomenon according to the change of the viscosity and yield stress of high flow concrete.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.299-300
• /
• 2003

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.16 no.1
• /
• pp.65-84
• /
• 2012

We present an application of the variational multiscale methodology to the computation of concentric annular pipe flow. Isogeometric analysis is utilized for higher order approximation of the solution using Non-Uniform Rational B-Splines (NURBS) functions. The ability of NURBS to exactly represent curved geometries makes NURBS-based isogeometric analysis attractive for the application to the flow through the curved channel.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.3
• /
• pp.341-349
• /
• 1997

The characteristics of defect correction method are discussed in a sample heat conduction problem showing the numerical solution of the error correction equation can predict the error of the numerical solution of the original governing equation. A way of using defect correction method combined with the existing algorithm for the incompressible fluid flow, is proposed and subsequently tested for the driven square cavity problem. The error correction equations for the continuity equation and the momentum equations are considered to estimate the errors of the numerical solutions of the original governing equations. With this new approach, better velocity and pressure fields can be obtained by correcting the original numerical solutions using the estimated errors. These calculated errors also can be used to estimate the orders of magnitude of the errors of the original numerical solutions.

• Journal of applied mathematics & informatics
• /
• v.22 no.1_2
• /
• pp.67-82
• /
• 2006

In this paper, we present a finite difference method for the implementation of the rotation of a circular cylinder in the incompressible flow field by solving the two-dimensional unsteady Navier-Stokes equations. The approach is to use staggered grid method so that the accuracy and order of convergence of the associated algorithms can be maintained. The proposed method is easy to be implemented and is effective. A set of simulations for the flow dynamics is provided to show the computational results.

• Korean Journal of Computational Design and Engineering
• /
• v.1 no.1
• /
• pp.56-64
• /
• 1996

A quadrilateral-triangular mixed grid method for the solution of incompressible viscous flow is presented. The solution domain near the body surface is meshed using elliptic grid geneator to acculately simulate the viscous flow. On the other hand, we used unstructured triangular grid system generated by advancing front technique of a simple automatic grid generation algorithm in the rest of the computational domain. The present method thus is capable of not only handling complex geometries but providing accurate solutions near body surface. The numerical technique adopted here is PISO type finite element method which was developed by the present author. Investigations have been made of two-dimensional unsteady flow of Re=550 past a circular cylinder. In the case of use of the unstructured grid only, there exists a considerable amount of difference with the existing results in drag coefficient and vorticity at the cylinder surface; this may be because of the lack of the grid clustering to the surface that is a inevitable requirement to resolve the viscous flow. However, numerical results on the mixed grid show good agreements with the earlier computations and experimental data.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.21 no.5
• /
• pp.71-79
• /
• 2017

A numerical analysis of the liquid-gas two-phase flows has been conducted. The incompressible equations of the two-phase flows were solved by the artificial compressibility method with the CLSVOF interface capturing method. To analyze the grid dependency of CLSVOF, a numerical analysis of Zalesak's disk and three-dimensional liquid deformation problem were carried out, and the reconstruction of deformation was investigated. The Rayleigh-Taylor instability was numerically analyzed by applying the equations of incompressible two-phase flow, and the surface instability was observed.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.389-396
• /
• 2010

In this paper, a new computational code was developed using Chorin's artificial compressibility method to solve the two-dimensional incompressible Navier-Stokes equations. In spatial derivatives, Roe's flux difference splitting was used for the inviscid flux, while central differencing was used for the viscous flux. Furthermore, AF-ADI with dual time stepping method was implemented for accurate unsteady computations. Two-equation turbulence models, Menter's $k-{\omega}$ SST model and Coakley's $q-{\omega}$ model, hae been adopted to solve high-Reynolds number flows. A number of numerical simulations were carried out for steady laminar and turbulent flow problems as well as unsteady flow problem. The code was verified and validated by comparing the results with other computational results and experimental results. The results of numerical simulations showed that the present developed code with the artificial compressibility method can be applied to slve steady and unsteady incompressible flows.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.