• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.24 no.4
• /
• pp.277-286
• /
• 2012

In the numerical simulation of wave fields using a multi-phase flow model that considers simultaneous flows of materials with different states such as gas, liquid and solid, there is need of an accurate representation of the interface separating the fluids. We adopted an algebraic interface capturing method called tangent of hyperbola for interface-capturing(THINC) method for the capture of the free-surface in computations of multi-phase flow simulations instead of geometrical-type methods such a volume of fluid(VOF) method. The THINC method uses a hyperbolic tangent functions to represent the surface, and compute the numerical flux for the fluid fraction functions. One of the remarkable advantages of THINC method is its easy applicability to incorporate various numerical codes based on Navier-Stokes solver because it does not require the extra geometric reconstruction needed in most of VOF-type methods. Several tests were carried out in order to investigate the advection of interfaces and to verify the applicability of the THINC method to wave fields based on the one-field model for immiscible two-phase flows (TWOPM). The numerical results revealed that the THINC method is able to track the interface between air and water separating the fluids although its algorithm is fairly simple.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.60-66
• /
• 2007

Separation control has been performed using synthetic jets on airfoil at high angle of attack. Computed results demonstrated that stall characteristics and control surface performance could be substantially improved by resizing separation vortices. It was observed that the actual flow control mechanism and flow structure is fundamentally different depending on the range of synthetic jet frequency. For low frequency range, small vortices due to synthetic jet penetrated to the large leading edge separation vortex, and as a result, the size of the leading edge vortex was remarkably reduced. For high frequency range, however, small vortex did not grow up enough to penetrate into the leading edge separation vortex. Instead, synthetic jet firmly attached the local flow and influenced the circulation of the virtual airfoil shape which is the combined shape of the main airfoil with the separation vortex. Theses results show the characteristic of unsteady flow of single synthetic jet. Beside, we researched on multi-array synthetic jet to obtain applicable synthetic jet velocity. Multi-location synthetic jet is proposed to eliminate small vortex on suction surface of airfoil. With the results, we concluded that the flow around airfoil is stable by high frequency synthetic jet with elimination of small vortex and confirmation of stable flow. Moreover, performance of multi-array/multi-location synthetic jet can be improved by changing phase angle of multi-location synthetic jet.

• Journal of Advanced Research in Ocean Engineering
• /
• v.1 no.1
• /
• pp.14-24
• /
• 2015

A three-dimensional multidirectional nonlinear numerical wave tank (NWT) based on the Navier-Stokes equations and the Finite Volume Method (FVM) is developed by using the two-phase hydrodynamic flow solver naoe-FOAM-SJTU based on the open source toolbox OpenFOAM. The free surface is capturing with the Volume Of Fluids (VOF). The directional wave including Stokes wave, solitary wave and nonlinear wave are simulated and verified. The multi-directional waves are also simulated with particular wave spectral such as JONSWAP and wave directional spreading function. The obtained numerical results show the capability of the solver to generate different type of multidirectional nonlinear waves accurately. Meanwhile, it implies that the presented NWT can easily extend to model the wave-structures interactions, which will be great help to the offshore structures design.

• International Journal of Ocean System Engineering
• /
• v.4 no.1
• /
• pp.49-56
• /
• 2014

A three-dimensional multidirectional nonlinear numerical wave tank (NWT) based on the Navier-Stokes equations and the Finite Volume Method (FVM) is developed by using the two-phase hydrodynamic flow solver naoe-FOAM-SJTU based on the open source toolbox OpenFOAM. The free surface is capturing with the Volume Of Fluids (VOF). The directional wave including Stokes wave, solitary wave and nonlinear wave are simulated and verified. The multi-directional waves are also simulated with particular wave spectral such as JONSWAP and wave directional spreading function. The obtained numerical results show the capability of the solver to generate different type of multidirectional nonlinear waves accurately. Meanwhile, it implies that the presented NWT can easily extend to model the wave-structures interactions, which will be great help to the offshore structures design.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.334-336
• /
• 2008

A novel adaptive mesh refinement (AMR) strategy based on the Moment-of-Fluid (MOF) method for volume-tracking dynamic interface computation is presented. The Moment-of-Fluid method is a new interface reconstruction and volume advection method using volume fraction as well as material centroid. The mesh refinement is performed based on the error indicator, the deviation of the actual centroid obtained by interface reconstruction from the reference centroid given by moment advection process. Using the AMR-MOF method, the accuracy of volume-tracking computation with evolving interfaces is improved significantly compared to other published results.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.334-336
• /
• 2008

A novel adaptive mesh refinement (AMR) strategy based on the Moment-of-Fluid (MOF) method for volume-tracking dynamic interface computation is presented. The Moment-of-Fluid method is a new interface reconstruction and volume advection method using volume fraction as well as material centroid. The mesh refinement is performed based on the error indicator, the deviation of the actual centroid obtained by interface reconstruction from the reference centroid given by moment advection process. Using the AMR-MOF method, the accuracy of volume-tracking computation with evolving interfaces is improved significantly compared to other published results.

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

In level set method, material properties are made to change smoothly across an interface of two materials with different properties by introducing an interpolation or smoothing scheme. So far, the weighted arithmetic mean (WAM) method has been exclusively adopted in level set method, without complete assessment for its validity. We showed here that the weighted harmonic mean (WHM) method for rate constants of various rate processes, including viscosity, thermal conductivity, electrical conductivity, and permittivity, gives much more accurate results than the WAM method. The selection of interpolation scheme is particularly important in multi-phase electrohydrodynamic problems in which driving force for fluid flow is electrical force exerted on the phase interface. Our analysis also showed that WHM method for both electrical conductivity and permittivity gives not only more accurate, but also more physically realistic distribution of electrical force at the interface. Our arguments are confirmed by numerical simulations of drop deformation under DC electric field.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.214-219
• /
• 2007

Both the bubble rising in a fully filled container and the droplet splash are simulated by a solution code(PowerCFD). This code employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with interface capturing method (CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The present results are compared with other numerical solutions found in the literature. It is found that the present code simulate complex free surface flows such as multi phase flows due to large density difference efficiently and accurately.

• Nuclear Engineering and Technology
• /
• v.44 no.8
• /
• pp.871-888
• /
• 2012

This paper aims at presenting the state of the art, the recent progress, and the perspective for the future, in the modelling of two-phase flow in the horizontal legs of a PWR. All phenomena relevant for safety analysis are listed first. The selection of the modelling approach for system codes is then discussed, including the number of fluids or fields, the space and time resolution, and the use of flow regime maps. The classical two-fluid six-equation one-pressure model as it is implemented in the CATHARE code is then presented and its properties are described. It is shown that the axial effects of gravity forces may be correctly taken into account even in the case of change of the cross section area or of the pipe orientation. It is also shown that it can predict both fluvial and torrential flow with a possible hydraulic jump. Since phase stratification plays a dominant role, the Kelvin-Helmholtz instability and the stability of bubbly flow regime are discussed. A transition criterion based on a stability analysis of shallow water waves may be used to predict the Kelvin-Helmholtz instability. Recent experimental data obtained in the METERO test facility are analysed to model the transition from a bubbly to stratified flow regime. Finally, perspectives for further improvement of the modelling are drawn including dynamic modelling of turbulence and interfacial area and multi-field models.

• Journal of computational fluids engineering
• /
• v.16 no.2
• /
• pp.17-23
• /
• 2011

The turbulent flow and vortex shedding phenomena around pantograph panhead of high speed train were investigated and compared with available experimental data and other simulations. The pantograph head was simplified to be a square-cross-section pillar and assumed to be no interference with other bodies. The Reynolds number (Re) was 22,000. The LES(large eddy simulation) of FDS code was applied to solve the momentum equations and the Wener-Wengle wall model was employed to solve the near wall turbulent flow. Smagorinsky model($C_s$=0.2) was used as SGS(subgrid scale) model. The total grid numbers were about 9 millions and the analyzed domain was divided into 12 multi blocks which were communicated with each other by MPI. The time-averaged mainstream flows were calculated and well compared with experimental data. The phased-averaged quantities had also a good agreement with experimental data. The near-wall turbulence should be carefully treated by wall function or direct resolution to get successful application of LES methods.