• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.407-413
• /
• 2011

VIC (Vortex-In-Cell) method for viscous incompressible flow is presented to simulate the wake behind a modified NACA16 foil. With uniform rectangular grid, the velocity in field is calculated using streamfunction from vorticity field by solving the Poisson equation in which FFT(Fast Fourier Transform) is combined with 2nd order finite difference scheme. Here, LES(Large Eddy Simulation) with Smagorinsky model is applied for turbulence calculation. Effective viscosity is formulated using magnitude of strain tensor(or vorticity). Then the turbulent diffusion as well as viscous diffusion becomes particle strength exchange(PSE) with averaged eddy viscosity. The well-established panel method is combined to obtain the irrotational velocity and to apply the no-penetration boundary condition on the body panel. And wall diffusion is used for no-slip condition numerical results of turbulent stresses are compared with experimental results (Bourgoyne, 2003). Before comparing process, LES(Large Eddy Simulation) SGS(Subgrid scale) stress is transformed Reynolds averaged stress (Winckelmans, 2001).

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.2
• /
• pp.275-283
• /
• 2006

Experiment and three dimensional numerical investigations of incompressible turbulent flow through square channels with one- and two-sided ribbed walls are performed to determine pressure drop and heat transfer. The CFX(version 5.7) software package is used for the computation. The ribbed walls have a $45^{\circ}$ inclined square rib. Uniform heat flux is maintained on whole inner heat transfer channel area. The numerical results coincide with experimental data that obtained for $7,600{\le}Re{\le}24.900$, the pitch-to-rib height ratio (p/e) of 8.0. and the rib height-to-channel hydraulic diameter ratio ($e/D_h$) of 0.0667. The results show that values of local heat transfer coefficient and friction factor in the channel with two-sided ribbed wall are higher than those in the channel with one-sided ribbed walls.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.825-830
• /
• 2004

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

• Journal of computational fluids engineering
• /
• v.11 no.3 s.34
• /
• pp.59-63
• /
• 2006

Accurate prediction of a supersonic missile base drag continues to defy even well-rounded CFD codes. In an effort to address the accuracy and predictability of the base drags, the influence of grid system and competitive turbulence models on the base drag is analyzed. Characteristics of some turbulence models is reviewed through incompressible turbulent flow over a flat plate, and performance for the base drag prediction of several turbulence models such as Baldwin-Loman(B-L), Spalart-Allmaras(S-A), k-$\varepsilon$, k-$\omega$ model is assessed. When compressibility correction is injected into the S-A model, prediction accuracy of the base drag is enhanced. The NSWC wind tunnel test data are utilized for comparison of CFD and semi-empirical codes on the accuracy of base drag predictability: they are about equal, but CFD tends to perform better. It is also found that, as angle of attack of a missile with control fins increases, even the best CFD analysis tool we have lacks the accuracy needed for the base drag prediction.

• Journal of Advanced Research in Ocean Engineering
• /
• v.2 no.2
• /
• pp.61-66
• /
• 2016

This study presents the wave run-up height and depression depth around offshore cylindrical structures according to the wave period. The present study employs the volume of fluid method with the realizable turbulence model based on a commercial computational fluid dynamics software called the "STAR-CCM+" to simulate a 3D incompressible viscous two-phase turbulent flow. The present results for the wave run-up height and depression depth with regard to the wave period are compared with those of the relevant previous experimental and numerical studies.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.84-85
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for a large scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady computation of the incompressible Navier-Stokes equation, and Smagorinsky LES(Large Eddy Simulation) model is chosen for turbulent flow computation. Both METIS and MPI library are used for domain partitioning and data communication between processors respectively. Tiburon of Hyundai-motor is chosen as the computational model at $Re=7.5{\times}10^{5}$, which is based on the car height. It is confirmed that the drag under road condition is smaller than that of wind tunnel condition.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2005.11a
• /
• pp.178-179
• /
• 2005

Numerical and experimental investigation of incompressible turbulent flow and heat transfer through square channels with varying number of ribbed walls were conducted to determined pressure drop and heat transfer. The CFX solver used for the computation. The rough walls have a $45^0$ inclined square rib. Uniform heat flux is maintained on whole inner heat transfer channel area. The numerical results agreed well with experimental data that obtained for 7600$D_h$) of 0.0667. The results show that values of local heat transfer coefficient and friction factor increase with an increasing number of ribbed walls.

• Journal of Hydrospace Technology
• /
• v.2 no.2
• /
• pp.1-13
• /
• 1996

This paper presents the result of a computational study on the wake characteristics of two tanker models, i.e. HSVA and Mystery hull forms. The focus of the study is on the distributions of axial, radial and tangential velocities of the two hull forms in way of the propeller, especially over the propeller disk. The effect of bilge vortices on the velocity distribution is also concerned. For the computation of stern and wake flows of the two hull farms, the incompressible Reynolds-Averaged Wavier-Stokes(RANS) equations are numerically solved by the second order finite difference method, which employs a four stage Runge-Kutta scheme with a residual averaging technique and the Baldwin-Lomax model. The calculated pressure distributions on the hull surface and the axial, radial and tangential velocity distributions over the propeller disk are presented for the two hull forms. Finally, the result of wake analysis for the computed wake distribution over the propeller disk is given in comparison with those for the experimental wake distribution fur the both hull forms.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.229-240
• /
• 1991

The large camber angle theory of turbomachine blade of compressor has been developed recently for the two-dimensional flow by Hawthorn, et al. However, in the above theory it was assumed that the fluid was incompressible and inviscid, and the blades had no thickness. In this study, the flow in a blade cascade being mounted in parallel fashion with blade of arbitrary thickness is studied in order to determine the effects of the camber angle on the performance characteristic of the blade section under the consideration of compressibility and viscosity of fluid. The panel method is used for potential flow analysis. The flow in the boundary-layer is obtained by solving the integral boundary-layer structure through the laminar, transitional , and turbulent flow using the pressure field determined from the potential flow. And then the viscous-inviscid interaction scheme is used for interaction of these two flows. For the determination of the variation in the outlet fluid angle influenced by deviation in cascade flow, the superposition method which is used for single foil is introduced in this analysis. By the introduction of this method, the effects of the deviation on outlet fluid angle and the resulting fluid angle are made to adjust for oneself through the calculation. As the result of this study, the blade of large camber angle, large incidence angle, large pitch-chord ratio has large viscous and compressible effect than those of small camber angle. Lift force increase as camber angle increases, but above 60.deg. of camber angle, lift force decrease as camber angle increases. But drag force increases linearly with camber angle increases in the entire region.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.