• Journal of Industrial Technology
• /
• v.22 no.A
• /
• pp.9-17
• /
• 2002

This analysis is aimed to find out how the conditions of secondary air injection affects the residence time and the turbulence energy of flue gas and flow field in a small incinerator. A commercial code, PHOENICS, is used to simulate the flow field of an Incinerator. The computational grid system is constructed in a cartesian coordinate system In this numerical experiment, an independent numerical variable is the conditions of secondary air injection and dependants are the residence time of flue gas and the mean value of turbulence energy in a primary combustion chamber. The flow field and the distribution of turbulence energy are analysed to evaluate the residence time of flue gas and the turbulence energy The computational results say that the tangential injection of secondary air make the residence time much longer than the radial injection and that the radial injection of secondary make turbulence much stronger than the tangential injection.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.49-57
• /
• 2013

Analysis of fluid-structure interaction for two nearby underwater vehicles immersed in the sea is quite challenging because simulation of flow around them is very difficult due to the complexity of underwater vehicle shapes. The conventional approach using body-fitted or unstructured grids demands much time in dynamic grid generation, and yields slow convergence of solution. Since an analysis of fluid-structure interaction must be based on accurate simulation results, a more efficient way of simulating flow around underwater vehicles, without sacrificing accuracy, is desirable. An immersed boundary method facilitates implementation of complicated underwater-vehicle shapes on a Cartesian grid system. An LES modeling is also incorporated to resolve turbulent eddies. In this paper, we will demonstrate the effectiveness of the immersed boundary method we adopted, by presenting the simulation results on the flow around a modeled high-speed underwater vehicle interacting with a modeled low-speed one.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.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.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.2737-2742
• /
• 2007

Strong wind flow around a building complex was numerically studied by LES. The original motivation of this work stemmed from the efforts to develop a risk assessment technique for windstorm hazards. Lagrangian-averaged scale-invariant dynamic subgrid-scale model was used for turbulence modeling, and a log-law-based wall model was employed on all the solid surfaces including the ground and the surface of buildings to replace the no-slip condition. The shape of buildings was implemented on the Cartesian grid system by an immersed boundary method. Key flow quantities for the risk assessment such as mean and RMS values of pressure on the surface of the selected buildings are presented. In addition, characteristics of the velocity field at some selected locations vital to safety of human beings is also reported.

• 한국전산유체공학회:학술대회논문집
• /
• 2010.05a
• /
• pp.397-403
• /
• 2010

The Immersed boundary method(IBM) is one of CFD techniques which can 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 these cases, pressure buildup near IB was found to occur when linear interpolation and stadard mass conservation is used and the interpolation scheme became complicated when higher order of interpolation is adopted. In this study, we implement the IBM to an incompressible Navier-Stokes solver which uses SIMPLE algorithm. 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.

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

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.1-10
• /
• 2010

In the present paper, a numerical method based on the constraint interpolation profile (CIP) method is applied for simulating two-dimensional violent sloshing problems. The free surface boundary value problem is considered as a multiphase problem which includes water and air. A stationary Cartesian grid system is adopted, and an interface capturing method is used to trace the shape of free surface profile. The CIP combined unified procedure (CCUP) scheme is applied for flow solver, and the tangent of hyperbola for interface capturing (THINC) scheme is used for interface capturing. Numerical simulations have been carried out for partially-filled 2D tanks under forced sway and roll motions at various filling depths and frequencies. The computational results are compared with experiments and/or the other numerical results to validate the present numerical method.

• Journal of computational fluids engineering
• /
• v.20 no.4
• /
• pp.51-57
• /
• 2015

In the present study, two-dimensional numerical investigation of flow past a square cylinder beneath the free surface has been performed to identify the effects of presence of the free surface. An immersed boundary method was adopted for implementation of the cylinder cross-section in a Cartesian grid system. Also, a level-set method was used to capture the interface of two fluids. To prevent transition to three-dimensional flow, Reynolds number chosen for this simulation was 150. The cases for Froude number 0.2 and gap ratio(h/D) between 0.25 and 5.00 were examined. At the specific Reynolds number, we study the effects of gap ratio on flow characteristics around a square cylinder by computing flow fields, force coefficients and Strouhal number.

• Proceedings of the KSME Conference
• /
• 2008.11b
• /
• 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.

• Nuclear Engineering and Technology
• /
• v.41 no.10
• /
• pp.1347-1360
• /
• 2009

A multi-dimensional component for the thermal-hydraulic system analysis code, MARS, was developed for a more realistic three-dimensional analysis of nuclear systems. A three-dimensional and two-fluid model for a two-phase flow in Cartesian and cylindrical coordinates was employed. The governing equations and physical constitutive relationships were extended from those of a one-dimensional version. The numerical solution method adopted a semi-implicit and finite-difference method based on a staggered-grid mesh and a donor-cell scheme. The relevant length scale was very coarse compared to commercial computational fluid dynamics tools. Thus a simple Prandtl's mixing length turbulence model was applied to interpret the turbulent induced momentum and energy diffusivity. Non drag interfacial forces were not considered as in the general nuclear system codes. Several conceptual cases with analytic solutions were chosen and analyzed to assess the fundamental terms. RPI air-water and UPTF 7 tests were simulated and compared to the experimental data. The simulation results for the RPI air-water two-phase flow experiment showed good agreement with the measured void fraction. The simulation results for the UPTF downcomer test 7 were compared to the experiment data and the results from other multi-dimensional system codes for the ECC delivery flow.