• Nuclear Engineering and Technology
• /
• 제42권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.

• 대한조선학회논문집
• /
• 제44권2호
• /
• pp.119-129
• /
• 2007

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2005년도 창립60주년 기념 추계 학술대회
• /
• pp.65.1-65.1
• /
• 2005

• 한국가시화정보학회:학술대회논문집
• /
• 한국가시화정보학회 2006년도 추계학술대회 논문집
• /
• pp.78-85
• /
• 2006

Flows around two oscillating cylinders in side-by-side arrangement at Re=185 are simulated using immersed boundary method. The cylinders oscillate vertically with prescribed sinusoidal function in opposite directions in uniform cross-flow. Flow patterns and drag & lift forces are described by varying distance between two cylinders and oscillating frequency. Time series of flow patterns are investigated along with corresponding drag k lift coefficients.

• 대한기계학회논문집B
• /
• 제40권11호
• /
• pp.717-726
• /
• 2016

본 연구는 밀폐계 내부에 고온의 원형 실린더가 존재할 때, 밀폐계의 종횡비 변화에 따른 밀폐계 내부의 3차원 자연대류 현상에 대해 수치해석을 수행하였다. 밀폐계 내부의 원형 실린더는 유한체적법(FVM)에 기초한 가상 경계법(IBM)을 사용하여 구현하였다. 본 연구에서 고려한 Rayleigh 수의 범위는 $10^5{\leq}Ra{\leq}10^6$이며, Prandtl 수는 0.7이다. 밀폐계의 폭을 변화하여 밀폐계의 종횡비를 증가시켰으며, 밀폐계의 종횡비는 $1{\leq}W/L{\leq}4$ 범위에서 1 간격으로 고려하였다. 본 연구에서 고려한 모든 Rayleigh 수와 밀폐계의 종횡비 범위에서 열유동장은 x=0 단면을 기준으로 좌우 대칭을 이루며 정상상태에 도달하였다. 또한 밀폐계의 종횡비가 증가할수록 원형 실린더의 표면 평균 Nusselt수는 증가하는 반면, 밀폐계 벽면의 표면 평균 Nusselt수는 감소하였다.

• 대한조선학회논문집
• /
• 제51권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.

• International Journal of Ocean System Engineering
• /
• 제1권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.

• Wind and Structures
• /
• 제26권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.

• 대한조선학회논문집
• /
• 제46권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.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 춘계 학술대회논문집
• /
• pp.105-110
• /
• 2001

An Eulerian-Lagrangian method, so called immersed boundary method, is used for analysing viscous flow around arbitrary bodies, where governing equations are discretized on a regular grid by using a finite volume method. To improve the accuracy of flow near body boundaries, a second-order accurate interpolation scheme is used and a level-set based grid deformation method is presented to construct the adaptive grids around body boundaries. The present scheme is used to simulate steady flow around a semicircular cylinder mounted on the bottom of flow domain and calculated results are validated by results of a body fitted grid method. Finally, present method is applied to a complex flow around multi body and the usefulness is checked by investigating calculated results.