• Journal of Ship and Ocean Technology
• /
• v.1 no.2
• /
• pp.19-30
• /
• 1997

The accurate prediction of the flow and the pressure distribution near the tip of the blade is crucial in determining the tip vortex cavitation inception which usually occurs on the blade tip or inside the core of the tip vortex just downstream of the blade tip. An improved boundary element method is applied to the prediction of the flow around propeller blades, with emphasis at the tip region. In the method, the Blow adapted grid and a higher order panel method, which combines a hyperboloidal panel geometry with a hi-quadratic dipole distribution, are used in order to accurately model the trailing wake geometry and the highly rolled-up regions in the wake. The method is applied to several propeller geometries and the results have been found to agree well to the existing experimental data. Inviscid flow methods are able to predict the pressures at the tip as well as the shape of the trailing wake. On the other hand, they are unable to determine the flow inside the viscous core of the tip vortex, where cavitation inception often occurs. Thus, a method is presented that treats the flow inside the viscous core. The inner flow is treated with a 2-D Clavier-stokes solution without making any assumptions for axisymmetric flow and conicity of the flow along the tip trajectory. The method can thus allow the treatment of general propeller blade configurations. The velocity and pressure distributions inside the core are shown and compared to those from other numerical methods.

• 한국전산유체공학회:학술대회논문집
• /
• 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).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.5
• /
• pp.526-534
• /
• 2004

A new second-moment closure model for turbulent heat fluxes is proposed on the basis of the elliptic equation. The new model satisfies the near-wall balance between viscous diffusion, viscous dissipation and temperature-pressure gradient correlation, and also has the characteristics of approaching its respective conventional high Reynolds number model far away from the wall. The predictions of turbulent heat transfer in a channel flow have been carried out with constant wall heat flux and constant wall temperature difference boundary conditions respectively. The velocity field variables are supplied from the DNS data and the differential equations only fur the mean temperature and the scalar flux are solved by the present calculations. The present model is tested by direct comparisons with the DNS to validate the performance of the model predictions. The prediction results show that the behavior of the turbulent heat fluxes in the whole region is well captured by the present model.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.2
• /
• pp.162-170
• /
• 2014

A CFD(Computational Fluid Dynamics) analysis is presented to predict hydrodynamic characteristics of a marine propeller. A commercial RANS(Reynolds Averaged Navier-Stokes equation) solver, namely FLUENT, is utilized in conjunction with fully unstructured meshes around rotating propeller. Mesh generation process is greatly accelerated by using fully unstructured meshes composed of both isotropic and anisotropic tetrahedral elements. The anisotropic tetrahedral elements were used in the flow domain near the blade and shaft, where the viscous effect is important, having complex shape yet resolving the thin boundary layers. For other regions, isotropic tetrahedral elements are utilized. Two different approaches simulating rotational effect of the propeller are employed, namely Moving reference frame technique for steady simulation, and Sliding mesh technique for unsteady simulation. Both approaches are applied to the propeller open water (POW) test simulation. The current results, which are thrust and torque coefficients, are compared with available experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.39 no.3
• /
• pp.41-47
• /
• 2002

Turbulent wake behind a ship propeller has been investigated using the adaptive hybrid 2-frame PTV(Particle Tracking Velocimetry). 400 instantaneous velocity fields were measured according to 4 different blade phases and ensemble-averaged to investigate the spatial evolution of the vortical structure of near wake within one propeller diameter downstream. The phase averaged mean velocity fields show the potential wake and the viscous wake formed by the boundary layers developed on the blade surfaces. As the tip vortex evolves downstream, the slipstream is contracted and the turbulent intensity is decreased with viscous dissipation and turbulent diffusion.

• Journal of Energy Engineering
• /
• v.26 no.1
• /
• pp.9-22
• /
• 2017

A combined thermal and flow analysis was carried out to study the behavior and performance of a small, commercial LTD (Low-Temperature-Differential) heat engine. Laminar-flow solutions for annulus and channel flows were employed to estimate the viscous drags on the piston and the displacer and the pressure difference across the displacer. Temperature correction factors were introduced to account for the departure from the ideal heat transfer processes. The analysis results indicate that the work required to overcome the viscous drags on engine moving parts and to move the displacer is much smaller than the moving-boundary work produced by the power piston for temperature differentials in the neighborhood of $20^{\circ}C$ and engine speeds below 10 RPS. A comparison with experimental data reveals large degradations from the ideal heat transfer processes. Thus, heat-transfer devices inside the displacer cylinder are recommended.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.3
• /
• pp.527-539
• /
• 1992

본 연구에서는 터빈익렬의 입구유동면에 주어지는 끝벽 경계층유동에 의하여 익렬 내의 유동에서 발생하는 여러 와류들에 의한 2차 유동과 이와 연관된 여러가지 3차원 점성유동 현상 그리고 이에 따른 유동손실을 보다 정확히 예측하기 위한 수치해 석적 연구를 수행하였으며, 이에 필요한 수치해석적 연구를 수행하였으며, 이에 필요 한 수치해석코드를 작성하였다.유동특성에 대하여 상세한 연구결과가 보고되어 있 는 UTRC(United Technologies Research Center) 평면 터빈익렬을 연구대상으로 채택하 여 익렬 내의 3차원 유동특성을 연구하고 계산한 결과를 기존의 결과와 비교 검토하였 다. 강한 2차유동이 존재하는 경우에 발생하는 수치확산을 감소시키기 위하여 대류 항에 대하여 2차 정확도(second-order accuracy)의 선형상류도식(linear upwind sche- me)을 사용하여 일반적으로 널리 사용되는 하이브리드도식(hybrid scheme)에 의한 해 석결과와 비교하였다. 터빈익렬 내의 난류 유동은 익렬의 회전과 유선의 만곡 등에 의한 영향으로 복잡한 유동현상을 나타내지만, 터빈익렬 내의 난류유동 특성에 대한 실험결과가 아직까지는 부족하고 또한 본 연구에서는 평균유동값의 정확한 해석에 중 점을 두었으므로 표준 k-.epsilon. 모델을 사용하였다.

• Journal of Korea Water Resources Association
• /
• v.45 no.4
• /
• pp.349-360
• /
• 2012

A horizontally curvilinear non-hydrostatic free surface model that was applicable to three-dimensional viscous flows was developed. The proposed model employed a top-layer equation to close kinematic free-surface boundary condition, and an isotropic k-${\varepsilon}$ model to close turbulence viscosity in the Reynolds averaged Navier-Stokes equation. The model solved the governing equations with a fractional step method, which solved intermediate velocities in the advection-diffusion step, and corrects these provisional velocities by accounting for source terms including pressure gradient and gravity acceleration. Numerical applications were implemented to the wind-driven currents in a two-dimensional closed basin, the flow in a steep-sided trench, and the flow in a strongly-curved channel accounting for secondary current by the centrifugal force. Through the numerical simulations, the model showed its capability that were in good agreement with experimental data with respect to free surface elevation, velocity, and turbulence characteristics.

• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.439-464
• /
• 2017

The hydro-elastic system consisting of a pre-stretched highly elastic plate, compressible Newtonian viscous fluid, and the rigid wall is considered and it is assumed that on the plate a lineal-located time-harmonic force acts. It is required to investigate the dynamic behavior of this system and determine how the problem parameters and especially the pre-straining of the plate acts on this behavior. The elasticity relations of the plate are described through the harmonic potential and linearized (with respect to perturbations caused by external time-harmonic force) form of these relations is used in the present investigation. The plane-strain state in the plate is considered and the motion of that is described within the scope of the three-dimensional linearized equations of elastic waves in elastic bodies with initial stresses. The motion of the fluid is described by the linearized Navier-Stokes equations and it is considered the plane-parallel flow of this fluid. The Fourier transform with respect to the space coordinate is applied for a solution to the corresponding boundary-value problem. Numerical results on the frequency response of the interface normal stress and normal velocity and the influence of the initial stretching of the plate on this response are presented and discussed. In particular, it is established that the initial stretching of the plate can decrease significantly the absolute values of the aforementioned quantities.

• Earthquakes and Structures
• /
• v.8 no.1
• /
• pp.275-304
• /
• 2015

The objective of this paper is to investigate the validity of transmitting boundaries in dynamic analysis of soil-structure interaction problems. As a case study, the proposed Baghdad metro line is considered. The information about the dimensions and the material properties of the concrete tunnel and surrounding soil were obtained from a previous study. A parametric study is carried out to investigate the effect of several parameters including the peak value of the horizontal component of earthquake displacement records and the frequency of the dynamic load. The computer program (Mod-MIXDYN) is used for the analysis. The numerical results are analyzed for three conditions; finite boundaries (traditional boundaries), infinite boundaries modelled by infinite elements (5-node mapped infinite element) presented by Selvadurai and Karpurapu, 1988), and infinite boundaries modelled by dashpot elements (viscous boundaries). It was found that the transmitting boundary absorbs most of the incident energy. The distinct reflections observed for the "fixed boundaries" disappear by using "transmitted boundaries". This is true for both cases of using viscous boundaries or mapped infinite elements. The type and location of the dynamic load represent two controlling factors in deciding the importance of using infinite boundaries. It was found that the results present significant differences when earthquake is applied as a base motion or a pressure load is applied at the surface ground. The peak value of the vertical displacement at nodes A, B, E and F (located at the tunnel's crown and side walls, and at the surface above the tunnel and at the surface 6.5 m away from tunnel's centre respectively) increases with the frequency of the surface pressure load for both cases 1 and 2 (traditional boundaries and mapped infinite elements respectively) while it decreases for case 3 (viscous boundaries). The modular ratio Ec/Es (modulus of elasticity of the concrete lining to that of the surrounding soil) has a considerable effect on the peak value of the horizontal displacement at node B (on the side wall of the tunnel lining) increase about (17.5) times, for the three cases (1, 2, and 3).