• Proceedings of the KSME Conference
• /
• 2000.11b
• /
• pp.745-750
• /
• 2000

A numerical code which can simulate unsteady, incompressible and 3-dimensional flows in an engine cylinder has been developed. The governing equations based on the cylindrical coordinate are discretized by the finite volume method with staggered variable arrangements. A geometric conservation rule is also incorporated into the simulation code in order to deal with a moving boundary problem. For the unsteady simulation, a fractional step method is adopted. The law of wall is applied to the wall boundaries and standard $k-\;{\varepsilon}$ model is used to describe the in-cylinder turbulent flow. The model cylinder has one eccentric port, flat piston and flat cylinder-head. The comparisons with experimental data show fairly well qualitative agreement.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.27-30
• /
• 2002

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. In this study, 3-dimensional unsteady incompressible Wavier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade' scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.08a
• /
• pp.152-157
• /
• 2003

As computer capacity has been progressed continuously, the studies of the flow characteristics have been performing by the numerical methods actively. Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In this study, 3-dimensional unsteady Incompressible Navier-Stokes equation was solved by numerical method using the fractional step method with the fourth order compact pade scheme to achieve high accuracy To validate the present code and algorithm, 3D flow-field around a cylinder was simulated. The drag coefficient and lift coefficient were computed and, then, compared with experiment. The present code will be tailored to LES simulation for more accurate turbulent flow analysis.

• Journal of computational fluids engineering
• /
• v.8 no.3
• /
• pp.50-55
• /
• 2003

Recent numerical simulation has a tendency to require the higher-order accuracy in time, as well as in space. This tendency is more true in LES and acoustic noise simulation. In the present work, the accuracy of a Fractional step method, which is widely used in LES simulation, has been increased to the fourth-order accurate compact Pade discretization. To validate the present code, the flow-field past a cylinder was simulated and compared with experiment. A good agreement with experiment was achieved.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.30-37
• /
• 2005

To evaluate LBM we performed the simulation of the unsteady two dimensional flow over a square cylinder in a channel in moderate Reynolds number range, $100\~500$ by using LBM and Fractional-Step method. Frist of all we compared LBM solution of Poiseuille flow applied Farout and periodic boundary conditions with the analytical solution to verify the applicability of the boundary conditions. For LBM simulation the calculation domain was formed by structured 500x100 grids. Prescribed maximum velocity and density inlet and Farout boundary conditions were imposed on the in-out boundaries. Bounceback boundary condition was applied to the channel and the cylinder waifs. The flow patterns and vortex shedding strouhal numbers were compared with previous research results. The flow patterns by LBM were in agreement with the flow pattern by fractional step method. Furthermore the strouhal number computed by LBM simulation result was more accurate than that of fractional step method through the comparison of the previous research results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.4
• /
• pp.458-465
• /
• 2003

The main objective of the research is to develop a research code solving transient incompressible Navier-Stokes equation. In this research code, Adams-Bashforth method was applied to the convective terms of the navier stokes equation and the splitted equations were discretized spatially by finite element methods to solve the complex geometry problems easily. To reduce the divergence on the boundaries of pressure poisson equation due to the unsuitable pressure boundary conditions, multi step approximation pressure boundary conditions derived from the boundary linear momentum equations were used. Simulations of Lid Driven Flow and Flow over Cylinder were conducted to prove the accuracy by means of the comparison with results of the previous workers.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.194-199
• /
• 2007

IB (immersed boundary) method is one of the prominent tool in computational fluid dynamics for the analysis of flows over complex geometries. The IB technique simplyfies the solution procedure by eliminating the requirement of complex body fitted grids and it is also superior in terms of memory requirement. In this study we have developed numerical code (FOTRAN) for the analysis of 2D flow over a cylinder using IB technique. The code is validated by comparing the wake lengths and separation angles given by Guo et. al. We employed fractional-step procedure for solving the Navier-Stokes equations governing the flow and discrete forcing IB technique for imposing boundary conditions. Also we have developed a 3D code for the backward-facing-step flow and flow over a sphere. The reattachment length in backward-facing-step flow was compared with the one given by Nie and Armaly, which has proven the validity of our code.

• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.1-8
• /
• 2004

As a computer has been continuously progressed to reduce R&D time and cost, the study of the flow physics has been significantly relied on the numerical method. Recently, Large Eddy Simulation(LES) has been widely used in CFD community to accurately capture the turbulent flows. The LES code requires high accuracy in time, as well as in space. Also, it should have strong robustness to ensure the convergence in various complicated flows. The objective of the present work is to develop a base code for LES simulation, having 2$^{nd}$ order accuracy in time and 4$^{th}$ order accuracy in space. To achieve the present objective, the four-step fractional step method was enhanced by adopting compact Pade'scheme. The standard Smagorinsky model was implemented for the first stage of the present code development. The flows over a cylinder and an airfoil were successfully simulated. and an airfoil were successfully simulated.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.71-76
• /
• 2004

An extension of our recently developed locally linear reconstruction scheme to 2 dimensional incompressible flow solver is presented. The solver is based on a semi-implicit fractional step method in which the convective term is discretized by Adams-Bashforth method and the diffusion term by Crank-Nicolson method. Several numerical examples are tested to demonstrate the mesh type independent accuracy of the solver, which include decaying vortex flow, square cavity flow, and flow around a circular cylinder. The above examples are solved on quadrilateral or hybrid meshes. For all numerical examples, we obtained reasonable results.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.483-487
• /
• 2011

In the present study, a CFD program based on a finite volume method was developed by using an unstructured polyhedral grid system for the accurate simulation with the complex geometry of computational domain. To simulate the transient flow induced by the moving solid object, the program used a fractional step method and a ALE (Algebric Lagrangian-Eulerian) method. The grid deformation for the moving of solid object were performed with a spring analogy based on the center coordinate of each computational grid. To verify the present program with these methodologies, the numerical results of the flow around the fixed and oscillating circular cylinder were compared with the previous numerical results.