• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.56-59
• /
• 2009

In the present work, LES with new variational multiscale method is conducted on the fully developed channel flow with Reynolds number is 180 based on the friction velocity and the channel half width. Incompressible Navier-Stokes equations are integrated using finite element method with the basis function of NURBS. To solve space-time equations, Newton's method with two stage predictor multicorretor algorithm is employed. The code is parallelized using MPI. The computational domain is a rectangular box of size $2{\pi}{\times}2{\times}4/3{\pi}$ in the streamwise, wall normal and spanwise direction. Mean velocity profiles and velocity fluctuations are compared with the data of DNS. The results agree well with those of DNS and other traditional LES.

• Journal of Advanced Marine Engineering and Technology
• /
• v.32 no.8
• /
• pp.1215-1220
• /
• 2008

To obtain the shape of the free surface more accurately, computations are carried out by a finite volume method using unstructured meshes and an interface capturing method. Free-surface flow, which is very important in the fields of ship and marine engineering, is numerically simulated for flows of both water and air. Control volumes are used with an arbitrary number of faces and allows a local mesh refinement. The integration is of second order, with a midpoint rule integration and linear interpolation. The method is fully implicit and uses quadratic interpolation. The solution method of pressure-correction type solves sequentially equations of momentum, continuity, conservation, and two-equations turbulence model. Comparison are quantitatively made between the computation and experiment in order to confirm the solution method.

• The KSFM Journal of Fluid Machinery
• /
• v.1 no.1 s.1
• /
• pp.72-80
• /
• 1998

Computational analysis of incompressible flows by numerically solving Navier-Stokes equations using multi-block finite volume method is conducted on a parallel computing system. Numerical algorithms adopted in this study $include^{(1)}$ QUICK upwinding scheme for convective $terms,^{(2)}$ central differencing for other terms $and^{(3)}$ the second-order Euler differencing for time-marching procedure. Structured grids are used on the body-fitted coordinate with multi-block concept which uses overlaid grids on the block-interfacing boundaries. Computational code is parallelized on the MPI environment. Numerical accuracy of the computational method is verified by solving a benchmark test case of the flow inside two-dimensional rectangular cavity. Computation in the axial compressor cascade is conducted by using 4 PE's md, as results, no numerical instabilities are observed and it is expected that the present computational method can be applied to the turbomachinery flow problems without major difficulties.

• Journal of Ocean Engineering and Technology
• /
• v.15 no.4
• /
• pp.46-52
• /
• 2001

SMAC method, one of the numerical simulation techniques, is modified from the original MAC method for the time-dependent variation of fluid flows. The Navier-Stokes equations for incompressible time-dependent viscous flow are applied, and Also marker particles which present the visualization of fluid flows are used. In this study, two-dimensional numerical simulations of the overtopping are carried out by SMAC method, and the simulation results are visualized, In addition to, motion pictures are made for efficient visualization of the simulation results. This numerical simulation could also be applied to the design of coastal structures as dike and revetment.

• 한국전산유체공학회:학술대회논문집
• /
• 2005.10a
• /
• pp.213-217
• /
• 2005

In this study the water including sediment is assumed that the density is different from fresh water. And the phenomenon inhaled by low pressure around the pipe is numerically simulated in two dimension. The simulation is done using finite difference method in rectangular staggered mesh system and Navier-Stokes equations and continuty equation are employed as governing equations. The method of Irregular leg lengths and stars are adopted to satisfy boundary condition of body boundary. Marker-Density method is used to calculate the density of mixed flow. Also SGS turbulence model is applied to consider vortex smaller then grid at high Renolds number. This study is to analyze inhalation phenomenon of mixed flow with sediment and to verify the numerical method for mixed flow. To verify the numerical results are compared with experimental results

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.10
• /
• pp.1317-1325
• /
• 2000

A new efficient numerical method for computing three-dimensional, unsteady, incompressible flows is presented. To eliminate the restriction of CFL condition, a fully-implicit time advancement in which the Crank-Nicolson method is used for both the diffusion and convection terms, is adopted. Based on an approximate block LU decomposition method, the velocity -pressure decoupling is achieved. The additional decoupling of the intermediate velocity components in the convection term is made for the fully -implicit time advancement scheme. Since the iterative procedures for the momentum equations are not required, the velocity components decouplings bring forth the reduction of computational cost. The second-order accuracy in time of the present numerical algorithm is ascertained by computing decaying vortices. The present decoupling method is applied to minimal channel flow unit with DNS (Direct Numerical Simulation).

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

We present an immersed boundary (IB) method for 3D simulation of flappingflags in a uniform flow. The proposed formulation is manipulated on the basis of an efficient Navier-Stokes solver adopting the fractional step method and a staggered Cartesian grid system. A direct numerical method is developed to calculate the flag motion, with the elastic force treated implicitly. The fluid motion defined on an Eulerian grid and the flag motion defined on a Lagrangian grid are independently solved and the mass of flag is handled in a natural way. An additional momentum forcing is formulated from the flag motion equation in a way similar with the direct-forcing IB formulation and acts as the interaction force between the flag and ambient fluid. A series of numerical tests are performed and the present results are compared qualitatively and quantitatively with previous studies. The instantaneous flag motion is analyzed under different conditions and surrounding vortical structures are identified. The effects of physical parameters on the flapping frequency are studied.

• Journal of the Korean Society of Visualization
• /
• v.14 no.2
• /
• pp.46-51
• /
• 2016

In this study, a Masked Omni-Directional Integration(MODI) method for pressure calculation is proposed using the Particle Image Velocimetry (PIV) data. To obtain the velocity field, the Affine PIV method was adopted. Synthetic images were generated for a solid body rotation. Calculation on the pressure was based on the Navier-Stokes equation. The results obtained by the MODI were compared with those obtained by theoretical pressure and by the Omni-Directional Integration(ODI) method. It was shown that the minimum error by the proposed MODI method was attained when the mask size was 1.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

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

A level set method is proposed to simulate the incompressible two-phase flow considering the effect of surface tension. For reinitialization of level set junction, a direct approach method is employed, instead of solving hyperbolic type equation. A mixed element is adopted, so that the continuity mid Navier-Stokes equations are solved by using the quadratic elements (six-node triangular element mid nine-node quadrilateral element), mid the level set function is solved by using the linear elements (three-node triangular element mid four-node quadrilateral element). In order to verify the accuracy mid robustness of the codes, the present methods are applied to a few benchmark problems. It is confirmed that the present results are in good qualitative mid quantitative agreements with the existing studies.