• Journal of computational fluids engineering
• /
• v.5 no.1
• /
• pp.1-7
• /
• 2000

In general, FDM(finite difference method) and FVM(finite volume method) are used for analyzing the fluid flow numerically. However it is difficult to apply them to problems involving complex geometries, multi-connected domains, and complex boundary conditions. On the contrary, FEM(finite element method) with coordinates transformation for the unstructured grid is effective for the complex geometries. Most of previous studies have used commercial codes such as KIVA or STAR-CD for the flow analyses in the engine cylinder, and these codes are mostly based on the FVM. In the present study, using the FEM for three-dimensional, unsteady, and incompressible Navier-Stokes equation, the velocity and pressure fields in the engine cylinder have been numerically analyzed. As a numerical algorithm, 4-step time-splitting method is used and ALE(arbitrary Lagrangian Eulerian) method is adopted for moving grids. In the Piston-Cylinder, the calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• Coupled systems mechanics
• /
• v.2 no.2
• /
• pp.191-214
• /
• 2013

The unsteady flow past a circular cylinder which starts rotating or rotary oscillating impulsively from rest in a viscous fluid is investigated for Reynolds numbers Re=200 and 1000, rectilinear speed ratios ${\alpha}$ between 0.5 and 5.0, and forced oscillating frequencies $f_s$ between 0.1 and 2.0. Numerical solutions of the Navier-Stokes equations are obtained by using a finite volume method on an unstructured colocated grid. The objective of the study is to examine the effect of the rotating and rotary oscillating circular cylinder on the flow patterns and dynamics loads. The numerical results reveal that the $K\acute{a}rm\acute{a}n$ vortex street vanishes entirely behind the rotating cylinder when the ratio ${\alpha}$ exceeds the critical value, and the vortex shedding behind the rotary oscillating cylinder undergoes mainly three modes named 'synchronization', 'competition' and 'natural shedding' with the increase of $f_s$. Based on the amplitude spectra analysis of the lift coefficients, the regions of the classification of flow structure modes are presented, which provide important references for the flow control in the ocean engineering.

• 한국연소학회:학술대회논문집
• /
• 2003.12a
• /
• pp.3-9
• /
• 2003

By utilizing a semi-empirical soot model, the applicability of the laminar flamelet concept for simulating the formation and oxidation of soot in the laminar diffusion flame has been studied. The source terms for two transport equations of the soot formation and oxidation are calculated in the mixture fraction/scalar dissipation rate space for laminar flamelets and stored in a library. In this study, emphasis is given to the interaction associated with radiation and soot formation. The radiative heat loss is obtained by solving the radiative transfer equation using the unstructured grid finite volume method with the WSGGM. The calculated temperatures and soot volume fractions agree relatively well with the experimental data and the previous numerical results of Kaplan et al. using the detailed chemistry.

• Journal of the Society of Naval Architects of Korea
• /
• v.44 no.1 s.151
• /
• pp.8-15
• /
• 2007

In this paper, numerical calculations are performed to analyze the unsteady flow of NACA airfoil sections. In order to ease the flow computation for the fluid region changing in time, improve the quality of solution and simplify the grid generation for the oscillating foil flow, the computational method adopts a moving and deforming mesh with the multi-block grid topology. The multi-block, structured-unstructured hybrid grid is generated using the commercial meshing software Gridgen V15. The MDM (Moving & Deforming Mesh) and the UDF (User Define function) function of FLUENT 6 are adopted for computing turbulent flows of the foil in pitching motion. Computed unsteady lift and drag forces are compared with experimental data. in general, the characteristics of unsteady lift and drag of the experiments are reproduced well in the numerical analysis.

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

An adaptive Cartesian grid method having the best elements of structured, unstructured, and Cartesian grids is developed to solve the steady two-dimensional Euler equations. The solver is based on a cell-centered finite-volume method with Roe's flux-difference splitting and implicit point Gauss-seidel time integration method. Calculations of several compressible flows are carried out to show the efficiency of the developed computer code. The results were generally in good agreements with existing data in the literature and the developed code has the good ability to capture important feature of the flows.

• Journal of Korea Water Resources Association
• /
• v.43 no.1
• /
• pp.105-117
• /
• 2010

This paper presents a new and simple technique to perform MUSCL reconstruction for solving 2D shallow water equations. The modified MUSCL scheme uses weighted area ratio to apply non-uniform grid in stead of the previous method that equally distributed the difference of conservation variables to each interface. The suggested method can physically reconstruct conservation variables in case of uniform grid as well as non-uniform grid. In this study, Unsplit scheme applicable to unstructured grid is used and efficient slope limiter of TVD scheme is used to control numerical oscillation which can be occurred in modified MUSCL scheme. For accurate and efficient treatment of bed slope term, the modified MUSCL scheme is coupled with the surface gradient method. The finite volume model applied to suggested scheme is verified through a comparison between numerical solution and laboratory measurements data such as the simulations of isolated building test case and Bellos's dam break test case.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.7 s.250
• /
• pp.603-611
• /
• 2006

This paper presents the response surface optimization method using three-dimensional Wavier-Stokes analysis to optimize the blade shape of an axial flow fan. Reynolds-averaged Wavier-Stokes equations with $k-{\epsilon}$ turbulence model are discretized with finite volume approximations using the unstructured grid. Regression analysis is used for generating response surface, and it is validated by ANOVA and t-statistics. Four geometric variables, i.e., sweep and lean angles at mean and tip respectively were employed to improve the efficiency. The computational results are compared with experimental data and the comparisons show generally good agreements. As a main result of the optimization, the total efficiency was successfully improved. Also, detailed effects of sweep and lean on the axial flow fan are discussed.

• Journal of computational fluids engineering
• /
• v.9 no.4
• /
• pp.48-54
• /
• 2004

Vortex-shedding flows past a circular cylinder for 200≤ Re ≤ 5000 are numerically simulated with the PowerCFD code, using a finite volume method and an unstructured grid system, developed by the author. The simulation is peformed by solving the unsteady 2-D Wavier-Stokes equations with both no model and turbulence model. The resulting Reynolds number dependence of the Strouhal number and of the drag and lift coefficients is compared with both experiments and previous numerical results. It is found that, in the range of 200≤ Re ≤ 5000 the calculation method with a turbulence model is capable of producing reasonably more accurate results than that with no model for the main practically relevant parameters such as Strouhal number, drag and lift coefficients.

• Journal of computational fluids engineering
• /
• v.9 no.4
• /
• pp.41-47
• /
• 2004

The paper presents numerical simulations of laminar vortex-shedding flows past a circular cylinder for Re ≤ 500. The simulations are performed by solving the unsteady 2-D Navier-Stokes equations with a finite volume method using unstructured grid system. The resulting Reynolds number dependence of the Strouhal number and of the drag and lift coefficients is compared with experiments and with previous numerical results, showing good agreement. It is found that, for the truly laminar Reynolds number range the present calculation method described is capable of producing reasonably accurate results for the main practically relevant parameters such as Strouhal number, drag and lift coefficients.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).