• Journal of the Society of Naval Architects of Korea
• /
• v.41 no.4
• /
• pp.1-8
• /
• 2004

As a computer has been continuously progressed to reduce R&D time and cost, the study of the flow physics has been significantly relied on the numerical method. Recently, Large Eddy Simulation(LES) has been widely used in CFD community to accurately capture the turbulent flows. The LES code requires high accuracy in time, as well as in space. Also, it should have strong robustness to ensure the convergence in various complicated flows. The objective of the present work is to develop a base code for LES simulation, having 2$^{nd}$ order accuracy in time and 4$^{th}$ order accuracy in space. To achieve the present objective, the four-step fractional step method was enhanced by adopting compact Pade'scheme. The standard Smagorinsky model was implemented for the first stage of the present code development. The flows over a cylinder and an airfoil were successfully simulated. and an airfoil were successfully simulated.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.27 no.7
• /
• pp.344-353
• /
• 2015

In the present study the characteristics of turbulent oscillatory flows in a square-sectional $180^{\circ}$curved duct were investigated experimentally. A series of experiments for air flow were conducted to measure axial velocity profiles, secondary flow velocity profiles and pressure distributions. The measurements were made by a Laser Doppler Velocimeter (LDV) system with a data acquisition and processing system which includes Rotating Machinery Resolve (RMR) and PHASE software. The results from the experiment are summarized as follows. (1) The maximum velocity moved toward the outer wall from the region of a bend angle of $30^{\circ}$. The velocity distribution had a positive value extended over the total phase in the region of a bend angle of $150^{\circ}$. (2) Secondary flows were generally proportional to the velocity of the main flow. The intensity of the secondary flow was about 25% as much as that in the axial direction. (3) Pressure distributions were effects of the oscillatory Dean number and respective region.

• Journal of Ship and Ocean Technology
• /
• v.7 no.4
• /
• pp.16-31
• /
• 2003

The finite volume based multi-block RANS code, WAVIS developed at KRISO, is used to simulate the turbulent flows around a submarine with the realizable $\textsc{k}-\varepsilon$ turbulence model. RANS methods are verified and validated at the level of validation uncertainty 1.54% of the stagnation pressure coefficient for the solution of the turbulent flows around SUBOFF submarine model without appendages. Another SUBOFF configuration, axisymmetric body with four identical stem appendages, is also computed and validated with the experimental data of the nominal wake and hydrodynamic coefficients. The hydrodynamic forces and moments for SUBOFF model and a practical submarine are predicted at several drift and pitch angles. The computed results are in extremely good agreement with experimental data. Furthermore, it is noteworthy that all the computations at the present study were carried out in a PC and the CPU time required for 2.8 million grids was about 20 hours to get fully converged solution. The current study shows that CFD can be a very useful and cost effective tool for the prediction of the hydrodynamic performance of a submarine in the basic design stage.

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

Computational simulations of turbulent flows around a model ship have been performed to investigate an influence of TVD schemes on the accuracy of advective terms associated with ship resistances. Several TVD schemes including upwind, second-order upwind, vanLeer, and QUICK as well as a nonTVD linear scheme were studied by examining temporal and spatial characteristics of accuracy transition in adjacent cells to the hull. Even though vanLeer scheme was the most accurate among TVD schemes in both structured and unstructured grid systems, the ratio of accuracy switch from 2^nd order to 1^st order in vanLeer scheme was considerable compared with the 2^nd order linear scheme. Also, the accuracy transition was observed to be overally scattered in the unstructured grid while the accuracy transition in the structured grid appeared relatively clustered. It concluded that TVD schemes had to be carefully used in computational simulations of turbulent flows around a model ship due to the loss of accuracy despite its attraction of numerical stability.

• 한국연소학회:학술대회논문집
• /
• 2006.10a
• /
• pp.202-208
• /
• 2006

The purpose of this research is to model numerically the turbulent gas-particle flows in a rectangular chamber using Eulerian-Eulerian Method. A computer code using the ${\kappa}-{\varepsilon}-Ap$ two-phase turbulence model is developed for the numerical study. This code and the Eulerian multiphase model in FLUENT were used for the numerical simulations of the two-phase flow in a rectangular chamber. The numerical results calculated by the two different turbulent gas-particle codes have shown that the ${\kappa}-{\varepsilon}-Ap$ model results in a stronger diffusion of the flow momentum in the gas-particle turbulence interaction than the Eulerian multiphase model in FLUENT.

• Journal of Advanced Marine Engineering and Technology
• /
• v.25 no.1
• /
• pp.131-138
• /
• 2001

In this paper, an experimental investigation of characteristics of developing turbulent steady flows in a square-sectional $180^{\circ}$curved duct is presented. The experimental study using air in a square-sectional $180^{\circ}$ curved duct carryed out to measure axials direction velocity and wall shear stress distrbutions by using Laser Dopper Velocimeter(LDV) system with data acquistion and processing the system of FIND6260 softwere at 7 sections from the inlet($\phi=0^{\circ}$) to the outlet($\phi=180^{\circ}$) in $301^{\circ}$ intervals of a curved duct.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.2
• /
• pp.206-214
• /
• 2003

Mean flow and turbulence properties of developing turbulent flows in a 90 degree square bend with spanwise rotation were measured by a hot-wire anemometer. A slanted wire is rotated into 6 orientations and the voltage outputs from them are combined to obtain the mean velocity and Reynolds stress components. The combinative effects of the centrifugal and Coriolis forces due to the curvature and the rotation of bend on the mean motion and turbulence structures are investigated experimentally. Results show that the two body forces can either enhance or counteract each other depending on the flow direction in the bend.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.66-72
• /
• 2000

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.12
• /
• pp.3345-3356
• /
• 1994

A general purpose program, TURBO-3D, for the analysis of 3-D fluid flow in complex geometry has been developed, which employs a standard $k-\varepsilon$ turbulence model and a general nonorthogonal grid system. For the purpose of verification of the program and testing the applicability, turbulent flows in an S-shaped diffuser and turbulent flows over an backward facing step are solved and compared with the earlier results. Comparison with the results by the STAR-CD program has been also made for the same flow configuration and grid structure. The agreements are excellent and hence the program has been verified. Since the present program is applicable only on limited flow phenomena and lacks the pre-and post processor, further improvements toward these directions are being made.

• Journal of Korean Society for Atmospheric Environment
• /
• v.12 no.3
• /
• pp.231-241
• /
• 1996

The effect of radiation and buoyancy on the thermophoresis phenomenon owing to the presence of highly absorbing, emitting particles (such as soot or pulverized coal) suspended in a two phase flow system was investigated numerically for a turbulent mixed convection flow. The analysis of conservation equations for a gas-particle flow system was performed on the basis of a two-fluid model from a continuum Eulerian viewpoint. The modified van Driest and Cebeci mixing length turbulence model was adopted in the anaylsis of turbulent flow. In addition, the P-1 approximation was used to evaluate the radiation heat transfer. As expected from the particle concentration and drift velocity distribution, the cumulative collection efficiency E (x) becomes larger when the buoyancy effect increases (i.e. higher Grashof number), while smaller as the radiation effect increases (i.e. higher optical thickness).