• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.355-362
• /
• 2009

Most material of engineering interest undergoes solidification process from liquid to solid state. Identifying the underlying mechanism during solidification process is essential to determine the microstructure of material which governs the physical properties of final product. In this paper, we expand our previous two-dimensional numerical technique to three-dimensional simulation for computing dendritic solidification process with fluid convection. We used Level Contour Reconstruction Method to track the moving liquid-solid interface and Sharp Interface Technique to correctly implement phase changing boundary condition. Three-dimensional results showed clear difference compared to two-dimensional simulation on tip growth rate and velocity.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.40 no.11
• /
• pp.717-726
• /
• 2016

This study evaluated the effect of aspect ratio of an enclosure with a heated inner circular cylinder on three-dimensional natural convection. The immersed boundary method was used to model the inner circular cylinder based on the finite volume method. The Rayleigh number was varied between $10^5$ and $10^6$, and the Prandtl number was maintained at 0.7. The aspect ratio of the three-dimensional enclosure was changed in steps of 1 within a range of 1-4 by increasing the width of the enclosure. In this study, the flow and thermal fields in the enclosure reached the steady state, and showed a mirror-symmetric pattern with respect to the center plane (x=0). In addition, the surface-averaged Nusselt number of the inner circular cylinder increased, while the total surface-averaged Nusselt number of the enclosure walls decreased with increase in the aspect ratio of the enclosure.

• Wind and Structures
• /
• v.26 no.5
• /
• pp.279-292
• /
• 2018

Inspired by the energy harvesting eel, a flexible flag behind a D-shape cylinder in a uniform viscous flow was simulated by using the immersed boundary method (IBM) along with low-speed wind tunnel experimentation. The flag in the wake of the cylinder was strongly influenced by the vortices shed from the upstream cylinder under the vortex-vortex and vortex-body interactions. Geometric and flow parameters were optimized for the flexible flag subjected to passive flapping. The influence of length and bending coefficient of the flexible flag, the diameters (D) of the cylinder and the streamwise spacing between the cylinder and the flag, on the energy generation was examined. Constructive and destructive vortex interaction modes, unidirectional and bidirectional bending and the different flapping frequency were found which explained the variations in the energy of the downstream flag. Voltage output and flapping behavior of the flag were also observed experimentally to find a more direct relationship between the bending of the flag and its power generation.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.1
• /
• pp.1-9
• /
• 2009

Three-dimensional characteristics of fluid flow and heat transfer around a wavy circular cylinder having sinusoidal variation in cross sectional area along the spanwise direction are numerically investigated using the immersed boundary method. The three different wavelengths of ${\pi}4$, ${\pi}3$ and ${\pi}2$ at the fixed wavy amplitude of 0.1 have been considered to investigate the effects of waviness especially on the forced convection heat transfer around a wavy cylinder when the Reynolds and Prandtl numbers are 300 and 0.71, respectively. The present computational results for a wavy cylinder are compared with those for a smooth cylinder. The time- and total surface-averaged Nusselt number for a wavy cylinder with ${\lambda}={\pi}/2$ is larger than that for a smooth cylinder, whereas that with ${\lambda}={\pi}/4$ and ${\pi}/3$ is smaller than that for a smooth cylinder. However, because the surface area exposed to heat transfer for a wavy cylinder is larger than that for a smooth cylinder, the total heat transfer rate for a wavy cylinder with different wavelengths of ${\lambda}={\pi}/4$, ${\pi}/3$ and ${\pi}/2$ is larger than that for a smooth cylinder.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.3
• /
• pp.248-257
• /
• 2011

In the present study, we proposed a new numerical wave tank model to analyze the vertical tension-leg circular floating bodies, using a 2-D Navier-Stokes solver. An IBM(Immersed Boundary Method) capable of handling interactions between waves and moving structures with complex geometry on a standard regular Cartesian grid system is coupled to the VOF(Volume of Fluid) method for tracking the free surface. Present numerical results for the motions of the floating body were compared with existing experimental data as well as numerical results based on FAVOR(Fractional Area Volume Obstacle Representation) algorithm. For detailed examinations of the present model, the additional hydraulic experiments for floating motions and free surface transformations were conducted. Further, the versatility of the proposed numerical model was verified via the numerical and physical experiments for the general rectangular floating bodies. Numerical results were compared with experiments and good agreement was archived.

• Nuclear Engineering and Technology
• /
• v.42 no.6
• /
• pp.620-635
• /
• 2010

In this paper we describe current activities on the project Multi-Scale Modeling and Analysis of convective boiling (MSMA), conducted jointly by the Paul Scherrer Institute (PSI) and the Swiss Nuclear Utilities (Swissnuclear). The long-term aim of the MSMA project is to formulate improved closure laws for Computational Fluid Dynamics (CFD) simulations for prediction of convective boiling and eventually of the Critical Heat Flux (CHF). As boiling is controlled by the competition of numerous phenomena at various length and time scales, a multi-scale approach is employed to tackle the problem at different scales. In the MSMA project, the scales on which we focus range from the CFD scale (macro-scale), bubble size scale (meso-scale), liquid micro-layer and triple interline scale (micro-scale), and molecular scale (nano-scale). The current focus of the project is on micro- and meso-scales modeling. The numerical framework comprises a highly efficient, parallel DNS solver, the PSI-BOIL code. The code has incorporated an Immersed Boundary Method (IBM) to tackle complex geometries. For simulation of meso-scales (bubbles), we use the Constrained Interpolation Profile method: Conservative Semi-Lagrangian $2^{nd}$ order (CIP-CSL2). The phase change is described either by applying conventional jump conditions at the interface, or by using the Phase Field (PF) approach. In this work, we present selected results for flows in complex geometry using the IBM, selected bubbly flow simulations using the CIP-CSL2 method and results for phase change using the PF approach. In the subsequent stage of the project, the importance of effects of nano-scale processes on the global boiling heat transfer will be evaluated. To validate the models, more experimental information will be needed in the future, so it is expected that the MSMA project will become the seed for a long-term, combined theoretical and experimental program.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.12 no.1
• /
• pp.343-353
• /
• 2020

Based on the Precise Transfer Matrix Method (PTMM), the dynamic model is constructed to observe the vibration behaviour of cylindrical shells with variable thickness by solving a set of first-order differential equations. The free vibration of stiffened cylindrical shells with variable thickness can be obtained to compare with the exact solution and FEM results. The reliability of the present method of free vibration is well proved. Furthermore, the effect of thickness on the vibration responses of the cylindrical shell is also discussed. The acoustic response of immersed cylindrical shells is analyzed by a Pluralized Wave Superposition Method (PWSM). The sound pressure coefficient can be gained by collocating points along the meridian line to satisfy the Neumann boundary condition. The mode convergence analysis of the cylindrical shell is carried out to guarantee calculation precision. Also, the reliability of the present method on sound radiation is verified by comparing with experimental results and numerical results.

• International Journal of Ocean System Engineering
• /
• v.1 no.3
• /
• pp.136-143
• /
• 2011

In this paper, three-dimensional free-surface flows are simulated by using two different numerical methods, the constrained interpolation profile (CIP)-based and finite volume (FV)-based methods. In the CIP-based method, the governing equations are solved on stationary staggered Cartesian grids by a finite difference method, and an immersed boundary technique is applied to deal with wave-body interactions. In the FV-based method, the governing equations are solved by applying collocated finite volume discretization, and body-fitted meshes are used. A free-surface boundary is considered as the interface of the multi-phase flow with air and water, and a volumeof-fluid (VOF) approach is applied to trace the free surface. Among many variations of the VOF-type method, the tangent of hyperbola for interface capturing (THINC) and the compressive interface capturing scheme for arbitrary meshes (CICSAM) techniques are used in the CIP-based method and FV-based method, respectively. Numerical simulations have been carried out for dam-breaking and wave-body interaction problems. The computational results of the two methods are compared with experimental data and their differences are observed.

• International Journal of High-Rise Buildings
• /
• v.8 no.4
• /
• pp.255-264
• /
• 2019

This paper describes large eddy simulation of wind pressures on a square cylinder in a uniform flow and a high-rise building immersed in an atmospheric turbulent boundary layer. For the atmospheric boundary layer case, the inflow turbulence is generated by a numerical wind tunnel. In the numerical simulation, particular attention is devoted to the performance of an auto hexahedral non-structural mesh. Both simulations are performed for three grid systems: an auto hexahedral non-structured grid, a structured Cartesian grid and a non-structured triangular prism grid, and for three grid numbers. The present study shows that the auto hexahedral unstructured mesh achieves the best simulation results for wind pressures on the square cylinder and the high-rise building. When the grid number is sufficiently large, the differences among the results obtained from the three investigated grid systems are not significant. However, the advantage of the auto hexahedral unstructured mesh becomes clear when the grid number decreases, because it enables a balanced distribution of orthogonal grids. The results described in this paper demonstrate that the auto hexahedral non-structured mesh has good potential applicability to simulation of urban flows.