• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.181-186
• /
• 1999

The effects of the frequency of upstream gust on the unsteady boundary characteristics on a downstream blade was simulated by using a Navier-Stokes code. The Navier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds k-e turbulence model to close the momentum equations. The MIT flapping foil experiment set-up is used to simulate the interaction between the upstream wake and a blade. The frequency of the upstream wake is simulated by varying rate of pitching motion of the flapping airfoils. Three reduced frequencies. 3.62. 7.24. and 10.86. are simulated. As the frequency increases, the unsteady fluctuation on the surfaces of the downstream hydrofoil is shown to decrease while the upstream flapper wake has larger first harmonics of y-velocity component. The unsteady vortices are shown to interact with each other and. as a result. the upstream wake becomes undiscernible inside the inner layer. The turbulence kinetic energy shows a similar behavior. Limiting streamlines around the trailing edge of the flapper are shown to conform with the unsteady Kutta condition for a round trailing edge. while limiting streamlines around the trailing edge of the hydrofoil conforms with the unsteady Kutta condition for a sharp edge.

• Journal of the Society of Naval Architects of Korea
• /
• v.58 no.3
• /
• pp.144-157
• /
• 2021

Brake Horse Power (BHP) reduction ratios by air injection to the underside of the hull surface in an actual ship are predicted using an unstructured finite-volume CFD solver and compared with the sea trial results. In addition, air lubrication system installed on the existing vessel is investigated to find a good solution for additional drag reduction. As a results, it is found that the thickness of the air layer should be minimized within a stable range while securing the area covered by the air layer as much as possible. Furthermore, the amount of frictional drag reduced by air injection is found to be independent of surface roughness and still effective on rough surface. Based on the results of this study, it is expected that systematic and reliable air lubrication system can be designed and evaluated using the proposed method.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.86-92
• /
• 2014

In this paper, prediction of separation trajectory for Two-stage-To-Orbit space launch vehicle has been numerically simulated by using an aerodynamic database based on steady state analysis. Aerodynamic database were obtained for matrix of longitudinal and vertical positions. The steady flow simulations around the launch vehicle have been made by using a 3-D RANS flow solver based on unstructured meshes. For this purpose, a vertex-centered finite-volume method was adopted to discretize inviscid and viscous fluxes. Roe's finite difference splitting was utilized to discretize the inviscid fluxes, and the viscous fluxes were computed based on central differencing. To validate this flow solver, calculations were made for the wind-tunnel experiment model of the LGBB TSTO vehicle configuration on steady state conditions. Aerodynamic database was constructed by using flow simulations based on test matrix from the wind-tunnel experiment. ANN(Artificial Neural Network) was applied to construct interpolation function among aerodynamic variables. Separation trajectory for TSTO launch vehicle was predicted from 6-DOF equation of motion based on the interpolated function. The result of present separation trajectory calculation was compared with the trajectory using experimental database. The predicted results for the separation trajectory shows fair agreement with reference[4] solution.

• Journal of Korea Water Resources Association
• /
• v.43 no.1
• /
• pp.105-117
• /
• 2010

This paper presents a new and simple technique to perform MUSCL reconstruction for solving 2D shallow water equations. The modified MUSCL scheme uses weighted area ratio to apply non-uniform grid in stead of the previous method that equally distributed the difference of conservation variables to each interface. The suggested method can physically reconstruct conservation variables in case of uniform grid as well as non-uniform grid. In this study, Unsplit scheme applicable to unstructured grid is used and efficient slope limiter of TVD scheme is used to control numerical oscillation which can be occurred in modified MUSCL scheme. For accurate and efficient treatment of bed slope term, the modified MUSCL scheme is coupled with the surface gradient method. The finite volume model applied to suggested scheme is verified through a comparison between numerical solution and laboratory measurements data such as the simulations of isolated building test case and Bellos's dam break test case.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.8 no.1
• /
• pp.39-45
• /
• 2005

In the situation which Russia's ratification of the Kyoto protocol at February,2004, ANNEX I nations must reduce GHG(Green House Gas) discharge rate from 2008 by 2012 to the reduction level at 1990. We introduce the CO₂ ocean sequestration that is one of promising method for getting the stable CO₂ concentration in the atmosphere. There are four categories : ocean transportation technique, ocean initial dissolution technique, ocean deep current evaluation technique, and ocean biological evaluation technique. In this paper, we carried out the fundamental numerical study on the ocean initial dissolution technique, when the Liquidized CO₂ is emitted at the deep ocean, It is very important to the dissolution rate of movable CO₂ interface because it Is directly impact to the ocean organism. In order to investigate the relation of the interface movement and rate of the dissolution, we develope CR(Computational Fluid Dynamics) code that was constructed by the finite volume method based on the unstructured mesh, and a droplet's boundary surface can move and one direction dissolution from disperse phase into continuous phase adopted as its physics be. This study clarifies hydrodynamic relation between solubility and movement of the droplet through the verification of the Cm code.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.6
• /
• pp.458-467
• /
• 2019

The purpose of this study is the suggestion of optimized parameters in OI (Optimal Interpolation) by experimental study. The observation of applying optimal interpolation is ADCP (Acoustic Doppler Current Profiler) data at the southwestern sea of Korea. FVCOM (Finite Volume Coastal Ocean Model) is used for the barotropic model. OI is to the estimation of the gain matrix by a minimum value between the background error covariance and the observation error covariance using the least square method. The scaling factor and correlation radius are very important parameters for OI. It is used to calculate the weight between observation data and model data in the model domain. The optimized parameters from the experiments were found by the Taylor diagram. Constantly each observation point requires optimizing each parameter for the best assimilation. Also, a high accuracy of numerical model means background error covariance is low and then it can decrease all of the parameters in OI. In conclusion, it is expected to have prepared the foundation for research for the selection of ocean observation points and the construction of ocean prediction systems in the future.

• Journal of Advanced Marine Engineering and Technology
• /
• v.36 no.8
• /
• pp.1091-1096
• /
• 2012

The aerodynamic characteristics of circular cylinder in a channel are studied to make clear the flow feature by solving the Navier-Stokes equation based on the finite volume method with unstructured grids. Reviews are made on with the vorticity, velocity, dynamic pressure, residual and drag, where the Reynolds numbers are 50 and 100. The flows for $Re{\succeq}50$ shows the vortex shedding in the wake, and the result is the same as the case of moving cylinder. The ground effect of flat bottom results in the growth of vortex, being generated in the upper side of the cylinder and elongated in the rear. As the cylinder approaches to wall, for example 0.6, the cylinder plays as a role of blockage to obstruct the flow between the cylinder and wall. The drag coefficients are compared with others' results to confirm the validity of the present numerical simulation.

• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1303-1313
• /
• 2002

Effects of the reduced frequency of upstream wake on downstream unsteady boundary layer flow were simulated by using a Wavier-Stokes code. The Wavier-Stokes code is based on an unstructured finite volume method and uses a low Reynolds number turbulence model to close the momentum equations. The geometry used in this paper is the MIT flapping foil experimental set-up and the reduced frequency of the upstream wake is varied in the range of 0.91 to 10.86 to study its effect on the unsteady boundary layer flow. Numerical solutions show that they can be divided into two categories. One is so called the low frequency solution, and behaves quite similar to a Stokes layer. Its characteristics is found to be quite similar to those due to either a temporal or spatial wave. The low frequency solutions are observed clearly when the reduced frequency is smaller than 3.26. The other one is the high frequency solution. It is observed for the reduced frequency larger than 7.24. It shows a sudden shift of the phase angle of the unsteady velocity around the edge of the boundary layer. The shift of phase angle is about 180 degree, and leads to separation of the boundary layer flow from corresponding outer flow. The high frequency solution shows the characteristics of a temporal wave whose wave length is half of the upstream frequency. This characteristics of the high frequency solution is found to be caused by the strong interaction between unsteady vortices. This strong interaction also leads to destroy of the upstream wake strips inside the viscous sublayer as well as the buffer layer.

• The Bulletin of The Korean Astronomical Society
• /
• v.39 no.1
• /
• pp.70.2-70.2
• /
• 2014

We developed a three-dimensional (3D) magnetohydrodynamic (MHD) simulation code to reproduce the structure of a solar wind and the propagation of a coronal mass ejection (CME) through it. This code is constructed by a finite volume method based on a total variation diminishing (TVD) scheme using an unstructured grid system (Tanaka 1994). The grid system can avoid the singularity arising in the spherical coordinate system. In this study, we made an improvement of the code focused on the propagation of a CME through a solar wind, which extends a previous work done by Nakamizo et al. (2009). We first reconstructed a solar wind in a steady state from physical values obtained at 50 solar radii away from the Sun via an MHD tomography applied to interplanetary scintillation (IPS) data (Hayashi et al. 2003). We selected CR2057 and inserted a spheromak-type CME (Kataoka et al. 2009) into a reconstructed solar wind. As a result, we found that our simulation well captures the velocity, temperature and density profiles of an observed solar wind. Furthermore, we successfully reproduce the general characteristics of an interplanetary coronal mass ejection (ICME) obtained by the Helios 1/2 spacecraft (R. J. FORSYTH et al. 2006).