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

The present work was undertaken to numerically analyze the defrosting phenomena of windshield glass. In order to analysis the phase change from frost to water on windshield glass by discharging hot air from a defroster nozzle, the flow and the temperature field of the cabin interior, the heat transfer through the windshield glass, and the phase change of frost should be solve simultaneously. In the present work, the flow field was obtained by solving 3-D incompressible Navier-Stokes equations, and the temperature field was computed from the incompressible energy equation. The phase change process was solved by the enthalpy method. For the code validation, the temperature and the phase change of the driven cavity were calculated. The calculation showed a good agreement with other numerical results. Then, the present code was applied to the defrosting problem of a real automobile, and a good agreement with the experimental data was also obtained.

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

Three-dimensional steady incompressible laminar entry flows in a square duct of $90^{\circ}$ bend are numerically simulated by a new solution code(PowerCFD) using unstructured cell-centered method. Solutions are obtained with three unstructured grid types of hexahedron, prism and hybrid at a Reynolds number, based on the hydraulic diameter and bulk velocity, of 790. Interesting features of the flow are presented in detail. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a $90^{\circ}$ bend of a square duct with fully-developed entry flows. It is found that the code is capable of producing the nature of laminar flow in curved square duct with no grid type dependency.

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

The turbulence noise generated from blunt trailing-edge is numerically predicted by using the hydrodynamic/acoustic splitting method at the Reynolds number based on thickness of flat plate, $Re_h=1000$, and the freestream Mach number $M_o=0.2$. The turbulent flow field is simulated by incompressible large-eddy simulation and the acoustic field is predicted efficiently with the linearized perturbed compressible equations (LPCE) recently proposed by the authors. The turbulent flow characteristics are validated with the results of the previous experimental study and direct numerical simulation. The acoustic properties predicted from LPCE are compared with the solutions of analytical formulations.

• Journal of computational fluids engineering
• /
• v.16 no.1
• /
• pp.73-82
• /
• 2011

In this paper, the vortex-induced vibration of circular cylinders is studied using the immersed boundary method on the Cartesian mesh. The Reynolds numbers considered is from 100 to 200. Using the configuration of tendemly arranged multiple circular cylinders, the vortex shedding behind of the cylinders and their flow-induced motion are investigated. The staggered MAC grid arrangement, which is the typical grid system for the incompressible flow on the Cartesian meshes, is utilized. Pressure correction method is applied for solving the divergence-free incompressible velocity field. The body motion is described by immersed boundary technique that has advantages for moving object on the fixed computational domain. It is also discussed for the computational noise in hydrodynamic forces when body motion is represented by the immersed boundary method. The Predictor/Corrector method is used for simulating the nonlinear response of the elastically mounted cylinder excited by vortex-shedding.

• 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).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.4
• /
• pp.242-247
• /
• 2009

AILU type preconditioners for a two-dimensional combined P2P1 finite element formulation of the interaction of rigid cylinder with incompressible fluid flow have been devised and tested by solving fluid-structure interaction (FSI) problems. The FSI code simulating the interaction of a rigid cylinder with an unsteady flow is based on P2P1 mixed finite element formulation coupled with combined formulation. Four different preconditioners were devised for the two-dimensional combined P2P1 finite element formulation extending the idea of Nam et al., which was proposed for the preconditioning of a P2P1 mixed finite element formulation of the incompressible Navier-Stokes equations. It was found that PC-III or PC-IV among them perform well with respect to computational memory and convergence rate for some bench-mark problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.9
• /
• pp.814-821
• /
• 2007

A hybrid method is presented for efficient computation of turbulent flow noise at low Mach numbers. In this method, the turbulent flow field is computed by incompressible large eddy simulation (LES), while the acoustic field is computed with the linearized perturbed compressible equations (LPCE) derived in this study. Since LPCE is computed on the rather coarse acoustic grid with the flow variables and source term obtained by the incompressible LES, the computational efficiency of calculation is greatly enhanced. Furthermore, LPCE suppress the instability of perturbed vortical mode and therefore secure consistent and stable acoustic solutions. The proposed LES/LPCE hybrid method is applied to three low Mach number turbulent flow noise problems: i) circular cylinder, ii) isolated flat plate, and iii) interaction between cylinder wake and airfoil. The computed results are closely compared with the experimental measurements.

• Journal of the Korean Mathematical Society
• /
• v.51 no.3
• /
• pp.635-654
• /
• 2014

We consider the Cauchy problem for a Keller-Segel-fluid model with degenerate diffusion for cell density, which is mathematically formulated as a porus medium type of Keller-Segel equations coupled to viscous incompressible fluid equations. We establish the global-in-time existence of weak solutions and bounded weak solutions depending on some conditions of parameters such as chemotactic sensitivity and consumption rate of oxygen for certain range of diffusive exponents of cell density in two and three dimensions.

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1737-1751
• /
• 2013

An efficient and stable point collocation scheme based on a meshfree method is studied for the stationary incompressible Navier-Stokes equations. We describe the diffuse derivatives associated with the moving least square method. Using these diffuse derivatives, we propose a point collocation method to fit in solving the Navier-Stokes equations which improves the stability of the direct point collocation scheme. The convergence of the numerical solution is investigated from numerical examples. The driven cavity ow and the backward facing step ow are implemented for the reliability of the scheme. Also, the viscous ow on complicated geometry is successfully calculated such as the ow past a circular cylinder in duct.

• Structural Engineering and Mechanics
• /
• v.6 no.6
• /
• pp.711-725
• /
• 1998

Mixed finite element formulations for incompressible materials show pressure oscillations or pressure modes in four-node quadrilateral elements. The criterion for the stability in the pressure solution is the so-called Babu$\check{s}$ka-Brezzi stability condition, and the four-node elements based on mixed variational principles do not appear to satisfy this condition. In this study, a pressure continuity residual based on the pressure discontinuity at element edges proposed by Hughes and Franca is used to study the stabilization of pressure solutions in bilinear displacement-constant pressure four-node quadrilateral elements. Also, a solid mechanics problem is presented by which the stability of mixed elements can be studied. It is shown that the pressure solutions, although stable, are shown to exhibit sensitivity to the stabilization parameters.