• 한국전산유체공학회지
• /
• 제11권4호
• /
• 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).

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2006년도 추계 학술대회논문집
• /
• pp.70-73
• /
• 2006

Mesh generation for the region of interest is prerequisite for numerical analysis of governing partial differential equations describing phenomena with proper physic. Mesh generation is, however, usually considered as a major obstacle for a routine application of numerical approaches in Engineering applications. Therefore automatic mesh generation is highly pursued. In this paper automated quadrilateral surface mesh generation is proposed. According to the present method, Cartesian cells of proper resolution for a region bounding the whole region of interest are first generated and the interior cells are identified. Then projecting their surface meshes onto the boundary surfaces gives surface mesh consisting of quadrilateral cells. This method has been implemented as an application program, and example cases are given.

• 천문학회보
• /
• 제41권1호
• /
• pp.50.1-50.1
• /
• 2016

The incorporation of radiation from massive stars is essential for modeling the dynamics and chemistry of star-forming clouds, yet it is a computationally demanding task for three-dimensional problems. We describe the implementation and tests of radiative transfer module due to point sources on a three-dimensional Cartesian grid in the Eulerian MHD code Athena. To solve the integral form of the radiation transfer equation, we adopt a widely-used long characteristics method with spatially adaptive ray tracing in which rays are split when sampling of cells becomes coarse. We use a completely asynchronous communication pattern between processors to accelerate transport of rays through a computational domain, a major source of performance bottleneck. The results of strong and weak scaling tests show that our code performs well with a large number of processors. We apply our radiation hydrodynamics code to some test problems involving dynamical expansion of HII regions.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2008년도 추계학술대회B
• /
• pp.2477-2482
• /
• 2008

Effect of computational domain size in simulating of periodic obstacle flow has been investigated for the flow past tube banks. Reynolds number, defined by freestream velocity (U) and cylinder diameter (d), was fixed as 200, and center-to-center distance (P) as 1.5d. In-line square array was considered. Drag coefficient, lift coefficient and Strouhal number were calculated depending on domain size. Circular cylinders were implemented on a Cartesian grid system by using an immersed boundary method. Boundary condition is periodic in both streamwise and lateral directions. Previous studies in literature often use a square domain with a side length of P, which contains only one cylinder. However, this study reveals that size is improper. Especially, RMS values of flow-induced forces are most sensitive to the domain size.

• 대한조선학회논문집
• /
• 제54권6호
• /
• pp.492-500
• /
• 2017

In the study, numerical simulation on the greenwater impact load of free surface offshore structure in the regular waves using fixed cartesian grid system and Modified Marker-Density (MMD) method were carried out and the results were reviewed. In order to compare numerical simulation and experimental results, the FPSO with the scale ratio of 1/100 model ship with fixed rectangular deck was selected and turbulence characteristic of the flow was considered by applying the Sub-Grid Scale (SGS) in laminar flow. As a result, it is reviewed how the greenwater impact load inflowed from bow in regular headsea wave influence the flow on the deck and the flow characteristic by numerical simulation and the experiment results were compared and reviewed. Based on this study, it would be useful to numerically study the effect of greenwater on offshore structure.

• 한국전산유체공학회지
• /
• 제16권2호
• /
• pp.1-10
• /
• 2011

A new computational program, which is based on the CIP/CCUP(Constraint Interpolation Profile/CIP Combined Unified Procedure) method, has been developed to numerically analyse sloshing phenomena dealt as multiphase-flow problems. For the convection terms of Navier-Stokes equations, the RCIP(Rational function CIP) method was adopted and the THINC-WLIC(Tangent of Hyperbola for Interface Capturing-Weighted Line Interface Calculation) method was used to capture the air/water interface. To validate the present numerical method, two-dimensional dam-breaking and sloshing problems in a rectangular tank were solved by the developed method in a stationary Cartesian grid system. In the case of sloshing problems, simulations by using a improved MPS(Moving Particle Simulation) method, which is named as PNU-MPS(Pusan National University-MPS), were also carried out. The computational results are compared with those of experiments and most of the comparisons are reasonably good.

• 대한기계학회논문집
• /
• 제15권2호
• /
• pp.668-676
• /
• 1991

Turbulent flows around tube banks and in the diffuser were studied using a non-orthogonal boundary fitted coordinate system and the modified K-.epsilon. turbulence model. In these cases, many problems emerge which stem from the geometrical complexity of the flow domain and the physical complexity of turbulent flow itself. To treat the complex geometry, governing equations were reformulated in a non-orthogonal coordinate system with Cartesian velocity components and discretised by the finite volume method with a non-staggered variable arrangement. The modified K-.epsilon. model of Hanjalic and Launer was applied to solve above two cases under the condition of strong and mild pressure gradient. The results using the modified K-.epsilon. model results in both test cases.

• 대한조선학회논문집
• /
• 제53권1호
• /
• pp.1-9
• /
• 2016

It is important to research and understand the physical phenomenon around a semi-submersible offshore structure on waves and currents because the wave run-up and load occurs owing to the waves and currents. In this study, the numerical simulations are performed about flow around a fixed semi-submersible offshore structure. The Modified Marker-density method is adopted in the present computation procedure, this method is one of the various methods to define the free-surface. The present computation results are compared with existing experimental and numerical simulation(VOF method) results. And, the computation results are relatively coincident with the existing results of model test and numerical simulation by VOF method.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• 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.

• 한국전산유체공학회지
• /
• 제17권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.