• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1601-1606
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.7
• /
• pp.809-817
• /
• 2004

In this paper, the unsteady behavior of the viscous flow field past an impulsively started elliptic cylinder is studied numerically. In order to analyze flow field, we introduce vortex particle method. The vorticity transport equation is solved by fractional step algorithm which splits into convection term and diffusion term. The convection term is calculated with Biot-Savart law, the no-through boundary condition is employed on solid boundaries. The diffusion term is modified based on the scheme of particle strength exchange. The particle redistributed scheme for general geometry is adapted. The flows around an elliptic cylinder are investigated for various attack angles at Reynolds number 200. The comparison between numerical results of present study and experimental data shows good agreements.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.3 no.3
• /
• pp.216-224
• /
• 2011

In This study develops a quasi-three dimensional numerical model of wave driven coastal currents with accounting the effects of the wave-current interaction and the surface rollers. In the wave model, the current effects on wave breaking and energy dissipation are taken into account as well as the wave diffraction effect. The surface roller associated with wave breaking was modeled based on a modification of the equations by Dally and Brown (1995) and Larson and Kraus (2002). Furthermore, the quasi-three dimensional model, which based on Navier-Stokes equations, was modified in association with the surface roller effect, and solved using frictional step method. The model was validated by data sets obtained during experiments on the Large Scale Sediment Transport Facility (LSTF) basin and the Hazaki Oceanographical Research Station (HORS). Then, a model test against detached breakwater was carried out to investigate the performance of the model around coastal structures. Finally, the model was applied to Akasaki port to verify the hydrodynamics around coastal structures. Good agreements between computations and measurements were obtained with regard to the cross-shore variation in waves and currents in nearshore and surf zone.

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.1026-1035
• /
• 2004

An immersed boundary method for solving the Navier-Stokes and thermal energy equations is developed to compute the heat transfer over or inside the complex geometries in the Cartesian or cylindrical coordinates by introducing the momentum forcing, mass source/sink, and heat source/sink. The present method is based on the finite volume approach on a staggered mesh together with a fractional step method. The method of applying the momentum forcing and mass source/sink to satisfy the no-slip condition on the body surface is explained in detail in Kim, Kim and Choi (2001, Journal of Computational Physics). In this paper, the heat source/sink is introduced on the body surface or inside the body to satisfy the iso-thermal or iso-heat-flux condition on the immersed boundary. The present method is applied to three different problems : forced convection around a circular cylinder, mixed convection around a pair of circular cylinders, and forced convection around a main cylinder with a secondary small cylinder. The results show good agreements with those obtained by previous experiments and numerical simulations, verifying the accuracy of the present method.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.17 no.1
• /
• pp.21-31
• /
• 2005

In this study, derivation is made of a one-grid source function for the extended Boussinesq equations of Nwogu (1993) in order to generate nonlinear waves internally. The energy velocity approach used in the line source method is verified analytically by the fractional step splitting method. The source function method is verified by generating accurately nonlinear waves as well as linear waves for horizontally one-dimensional cases. It is found that numerical solutions by the source function method are the same as those by the line source method.

• International Journal of Reliability and Applications
• /
• v.7 no.2
• /
• pp.167-175
• /
• 2006

Effective and quick methods that are easy to carry out even by hand, or easy to be programmed by hand-held calculators are needed for assessing the sizes of contrasts of unreplicated $2^P$ factorial designs. Moreover, they have the advantage to use the original numerical measurements which makes the analysis easier to explain. Basically, Lenth (1989) is one of the most familiar of such quick and powerful methods. Later on, Aboukalam (2001) proposes under constant effects an alternative sophisticated method to Lenth's method. The proposed method is the supreme from two considerable powers. The first utmost indicates less inflation deficiency while the other utmost indicates less shrinkage deficiency. Also under constant effects, Al-Shiha (2006) introduces an alternative quick method which is less shrinkage deficiency while the inflation deficiency is the same. If effects are random, Aboukalam (2005) introduces an alternative quick method in which the first power is favored as long as the second power is within a small margin. In the spirit of quickness and fixed effects, this article adds another method which is supreme from the two considerable powers. The method is based on a one step of the scale-part of a suggested M-estimate for location. Explicitly, we suggest adapting the skipped median (ASKM) estimate. Critical values of ASKM-method, for several sample sizes often used, are empirically computed.

• Wind and Structures
• /
• v.15 no.3
• /
• pp.189-208
• /
• 2012

Lock-in and drag amplification phenomena are studied for a flexible cantilever using a simplified fluid-structure interaction approach. Instead of solving the 3D domain, a simplified setup is devised, in which 2D flow problems are solved on a number of planes parallel to the wind direction and transversal to the structure. On such planes, the incompressible Navier-Stokes equations are solved to estimate the fluid action at different positions of the line-like structure. The fluid flow on each plane is coupled with the structural deformation at the corresponding position, affecting the dynamic behaviour of the system. An Arbitrary Lagrangian-Eulerian (ALE) approach is used to take in account the deformation of the domain, and a fractional-step scheme is used to solve the fluid field. The stabilization of incompressibility and convection is achieved through orthogonal quasi-static subscales, an approach that is believed to provide a first step towards turbulence modelling. In order to model the structural problem, a special one-dimensional element for thin walled cross-section beam is implemented. The standard second-order Bossak method is used for the time integration of the structural problem.

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

Study on swimming of microorganisms like, sperm motility, cilia beating, bacterial flagellar propulsion has found immense significance in the field of biological fluiddynamics. Because of the complexity involved, it is challenging for the researchers to model such problems. Immersed boundary method has proved its efficacy in the field of biological fluiddynamics, The present work aims at performing a numerical study on the microorganism locomotion using the immersed boundary method proposed by Peskin[1]. A two-dimensional model of the microorganism is modeled as thin elastic filament described as a sine wave. The neutrally buoyant organism undergoing deformations is immersed in a viscous and incompressible fluid. The fluid quantities are described using Eulerian coordinates and the immersed body is represented by Lagrangian coordinates. The Eulerian and Lagrangian variables are connected by the Dirac delta function. The Navier-Stokes equations governing the fluid flow are solved using the fractional step method on a staggered Cartesian grid system. The developed numerical code in FORTRAN will be validated by comparing the numerical results with the available results.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.44-53
• /
• 2012

Immersed boundary method(IBM) is a numerical scheme proposed to simulate flow field around complex objectives using simple Cartesian grid system. In the previous studies, the IBM has mostly been implemented to fractional step method based Navier-Stokes solvers. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. The weight coefficients of the bi-linear and quadratic interpolation equations were formulated by using only geometric information of boundary to reconstruct velocities near IB. Flow around 2D circular cylinder at Re=40 and 100 was solved by using these formulations. It was found that the pressure buildup was not observed even when the bi-linear interpolation was adopted. The use of quadratic interpolation made the predicted aerodynamic forces in good agreement with those of previous studies. For an analysis of moving boundary, we smulated an oscillating circular cylinder with Re=100 and KC(Keulegan-Carpenter) number of 5. The predicted flow fields were compared with experimental data and they also showed good agreements.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.6
• /
• pp.753-765
• /
• 2003

A finite element code for the numerical solution of the Navier-Stokes equation is parallelized by vertex-oriented domain decomposition. To accelerate the convergence of iterative solvers like conjugate gradient method, parallel block ILU, iterative block ILU, and distributed ILU methods are tested as parallel preconditioners. The effectiveness of the algorithms has been investigated when P1P1 finite element discretization is used for the parallel solution of the Navier-Stokes equation. Two-dimensional and three-dimensional Laplace equations are calculated to estimate the speedup of the preconditioners. Calculation domain is partitioned by one- and multi-dimensional partitioning methods in structured grid and by METIS library in unstructured grid. For the domain-decomposed parallel computation of the Navier-Stokes equation, we have solved three-dimensional lid-driven cavity and natural convection problems in a cube as benchmark problems using a parallelized fractional 4-step finite element method. The speedup for each parallel preconditioning method is to be compared using upto 64 processors.