• Journal of Computational Design and Engineering
• /
• v.3 no.4
• /
• pp.349-362
• /
• 2016

The numerical solutions of unsteady hydromagnetic natural convection Couette flow of a viscous, incompressible and electrically conducting fluid between the two vertical parallel plates in the presence of thermal radiation, thermal diffusion and diffusion thermo are obtained here. The fundamental dimensionless governing coupled linear partial differential equations for impulsive movement and uniformly accelerated movement of the plate were solved by an efficient Finite Element Method. Computations were performed for a wide range of the governing flow parameters, viz., Thermal diffusion (Soret) and Diffusion thermo (Dufour) parameters, Magnetic field parameter, Prandtl number, Thermal radiation and Schmidt number. The effects of these flow parameters on the velocity (u), temperature (${\theta}$) and Concentration (${\phi}$) are shown graphically. Also the effects of these pertinent parameters on the skin-friction, the rate of heat and mass transfer are obtained and discussed numerically through tabular forms. These are in good agreement with earlier reported studies. Analysis indicates that the fluid velocity is an increasing function of Grashof numbers for heat and mass transfer, Soret and Dufour numbers whereas the Magnetic parameter, Thermal radiation parameter, Prandtl number and Schmidt number lead to reduction of the velocity profiles. Also, it is noticed that the rate of heat transfer coefficient and temperature profiles increase with decrease in the thermal radiation parameter and Prandtl number, whereas the reverse effect is observed with increase of Dufour number. Further, the concentration profiles increase with increase in the Soret number whereas reverse effect is seen by increasing the values of the Schmidt number.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.12
• /
• pp.1-8
• /
• 2005

A parametric study has been accomplished to figure out the effects of the elliptic cylinder thickness, angle of attack, and Reynolds number on the lift and drag forces exerted on the elliptic cylinder. A two-dimensional incompressible Navier-Stokes flow solver is developed using SIMPLER method to analyze the unsteady viscous flow over elliptic cylinder. Thickness-to-chord ratios of 0.2, 0.4, and 0.6 elliptic cylinders are simulated at different Reynolds numbers of 400 and 600, and angles of attack of $10^{\circ}$, $20^{\circ}$, and $30^{\circ}$. Through this study, it is observed that the elliptic cylinder thickness, angle of attack, and Reynolds number affect significantly not only the time-mean values and the amplitudes of the drag and lift forces but also the frequencies of the force oscillations.

• Journal of the Society of Naval Architects of Korea
• /
• v.54 no.1
• /
• pp.18-25
• /
• 2017

This study deals with numerically comparing performance according to rudder shape called 'Twisted rudder and Wavy twisted rudder'. In comparison with conventional rudder, rudder with wavy shape has showed a better performance at high angles of attack($30^{\circ}{\sim}40^{\circ}$) due to delaying stall. But most of study concerned with wavy shape had been performed in uniform flow condition. In order to identify the characteristics behind a rotating propeller, the present study numerically carries out an analysis of resistance and self-propulsion for KCS with twisted rudder and wavy twisted rudder. The turbulence closure model, Realizable $k-{\epsilon}$, is employed to simulate three-dimensional unsteady incompressible viscous turbulent and separation flow around the rudder. The simulation of self-propulsion analysis is performed in two step, because of finding optimization case of wavy shape. The first step presents there are little difference between twisted rudder and case of H_0.65 wavy twisted rudder in delivered power. So two kind of rudders are employed from first step to compare lift-to-drag ratio and torque at high angles of attack. Consequently, the wavy twisted rudder is presented as a possible way of delaying stall, allowing a rudder to have a better performance containing superior lift-to-drag ratio and torque than twisted rudder at high angles of attack. Also, as we indicate the flow visualization, check the quantity of separation flow around the rudder.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• pp.164-169
• /
• 2009

Study on swimming of microorganisms like, sperm motility, cilia beating, bacterial flagellar propulsion has found immense significance in the field of biological fluiddynamics. Because of the complexity involved, it is challenging for the researchers to model such problems. Immersed boundary method has proved its efficacy in the field of biological fluiddynamics, The present work aims at performing a numerical study on the microorganism locomotion using the immersed boundary method proposed by Peskin[1]. A two-dimensional model of the microorganism is modeled as thin elastic filament described as a sine wave. The neutrally buoyant organism undergoing deformations is immersed in a viscous and incompressible fluid. The fluid quantities are described using Eulerian coordinates and the immersed body is represented by Lagrangian coordinates. The Eulerian and Lagrangian variables are connected by the Dirac delta function. The Navier-Stokes equations governing the fluid flow are solved using the fractional step method on a staggered Cartesian grid system. The developed numerical code in FORTRAN will be validated by comparing the numerical results with the available results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.10 s.241
• /
• pp.1163-1171
• /
• 2005

The shape optimization of an airfoil has been performed for an incompressible viscous flow. In this study, Pareto frontier sets, which are global and non-dominated solutions, can be obtained without various weighting factors by using the multi-objective genetic algorithm An NACA0012 airfoil is considered as a baseline model, and the profile of the airfoil is parameterized and rebuilt with four Bezier curves. Two curves, front leading to maximum thickness, are composed of five control points and the rest, from maximum thickness to tailing edge, are composed of four control points. There are eighteen design variables and two objective functions such as the lift and drag coefficients. A generation is made up of forty-five individuals. After fifteenth evolutions, the Pareto individuals of twenty can be achieved. One Pareto, which is the best of the . reduction of the drag furce, improves its drag to $13\%$ and lift-drag ratio to $2\%$. Another Pareto, however, which is focused on increasing the lift force, can improve its lift force to $61\%$, while sustaining its drag force, compared to those of the baseline model.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.4
• /
• pp.328-333
• /
• 2014

A hybrid particle-mesh method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations is a combination of the Vortex-In-Cell(VIC) method for convection and the penalization method for diffusion. The key feature of the numerical methods is to determine velocity and vorticity fields around a solid body on a temporary grid, and then the time evolution of the flow is computed by tracing the convection of each vortex element using the Lagrangian approach. Assuming that the vorticity and velocity fields are to be computed in time domain analysis, pressure fields are estimated through a complete set of solutions at present time step. It is possible to obtain vorticity and velocity fields prior to any pressure calculation since the pressure term is eliminated in the vorticity-velocity formulation. Therefore, pressure field is explicitly treated by solving a suitable Poisson equation. In this paper, we propose a simple way to numerically implement the vorticity-velocity-pressure formulation including a penalty term. For validation of the proposed numerical scheme, we illustrate the early development of viscous flows around an impulsive started circular cylinder for Reynolds number of 9500.

• Journal of Korea Water Resources Association
• /
• v.34 no.6
• /
• pp.605-615
• /
• 2001

SMAC method, one of the numerical simulation techniques, is modified from the original MAC method for the time-dependent variation of fluid flows. The Navier-Stokes equations for incompressible time-dependent viscous flow is applied, and marker particles which present the visualization of fluid flaws are used. In this study, two-dimensional numerical simulations of the water column collapse are carried out by SMAC method, and the simulation results are compared with Martin and Moyce's experimental data and result of the MPS method. A good results are obtained. This numerical simulation could also be applied to the breaking phenomenon of hydraulic structures such as dam break.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.38-46
• /
• 2010

This study aims at investigating the effect of the baffle height on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finite-volume method which has been well verified by comparing with the results of the relevant previous researches. The ratio of the baffle height ($h_B$) to filling level (h) has been changed in the range of $0{\leq}h_B/h{\leq}1.2$ to observe the effect on the impact loads on the side wall and free surface behavior. Generally, as baffle height increases, the impact pressure on the wall decreases and the deformation of free surface becomes weaker. However it seemed that a critical ratio of the baffle height existed to reveal the lowest impact pressure on the wall. Consequently, $h_B/h=0.8$ among $h_B/hs$ considered in the study showed the lowest impact pressure.

• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.4
• /
• pp.508-516
• /
• 2010

This study aims at investigating the effect of the chamfer on the liquid sloshing in the three-dimensional (3D) rectangular tank. In order to simulate the 3D incompressible viscous two-phase flow in the 3D tank with partially filled liquid, the present study has adopted the volume of fluid (VOF) method based on the finitevolume method which has been well verified by comparing with the results of the relevant previous researches. The effects of the chamfering top corners of the tank on the liquid sloshing characteristics have been investigated. The angle of the chamfering top corners (${\theta}$) has been changed in the range of $0^{\circ}{\leq}{\theta}{\leq}60^{\circ}$(${\Delta}{\theta}=15^{\circ}$) to observe the free surface behavior, and the effect on wall impact load. Generally, as the angle of the chamfering top corners increases, the impact pressure on the upper knuckle point decreases. However it seemed that a critical angle of the chamfering top corners exists to reveal the lowest impact pressure on the wall.