• 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.

• Coupled systems mechanics
• /
• v.2 no.2
• /
• pp.191-214
• /
• 2013

The unsteady flow past a circular cylinder which starts rotating or rotary oscillating impulsively from rest in a viscous fluid is investigated for Reynolds numbers Re=200 and 1000, rectilinear speed ratios ${\alpha}$ between 0.5 and 5.0, and forced oscillating frequencies $f_s$ between 0.1 and 2.0. Numerical solutions of the Navier-Stokes equations are obtained by using a finite volume method on an unstructured colocated grid. The objective of the study is to examine the effect of the rotating and rotary oscillating circular cylinder on the flow patterns and dynamics loads. The numerical results reveal that the $K\acute{a}rm\acute{a}n$ vortex street vanishes entirely behind the rotating cylinder when the ratio ${\alpha}$ exceeds the critical value, and the vortex shedding behind the rotary oscillating cylinder undergoes mainly three modes named 'synchronization', 'competition' and 'natural shedding' with the increase of $f_s$. Based on the amplitude spectra analysis of the lift coefficients, the regions of the classification of flow structure modes are presented, which provide important references for the flow control in the ocean engineering.

• Journal of computational fluids engineering
• /
• v.4 no.1
• /
• pp.53-61
• /
• 1999

Two-dimensional, unsteady, incompressible and compressible Navier-Stokes codes are developed for the computation of the viscous turbulent flow over high-lift airfoils. The compressible code involves a conventional upwind-differenced scheme for the convective terms and LU-SGS scheme for temporal integration. The incompressible code with pseudo-compressibility method also adopts the same schemes as the compressible code. Three two-equation turbulence models are evaluated by computing the flow over single and multi-element airfoils. The compressible and incompressible codes are validated by predicting the flow around the RAE 2822 transonic airfoil and the NACA 4412 airfoil, respectively. In addition, both the incompressible and compressible code are used to compute the flow over the NLR 7301 airfoil with flap to study the compressible effect near the high-loaded leading edge. The grid systems are efficiently generated using Chimera overlapping grid scheme. Overall, the κ-ω SST model shows closer agreement with experiment results, especially in the prediction of adverse pressure gradient region on the suction surfaces of high-lift airfoils.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.2
• /
• pp.22-28
• /
• 1999

Small amplitude in-phase oscillations of multi-cylinders are considered both numerically and experimentally. Flow field is separated into inner and outer regions. The basic unsteady solution is obtained analytically and the secondary flow, termed as steady streaming flow, can be obtained numerically by using Finite Volume Code with Panel Method. The Particle Induced Velocimetry, one of whole field measurements, is introduced for comparison with numerical flow visualization quantitatively. Among the algorithms for PIV, the Three Step Vector Searching Technique is applied to reduce CPU time. Small but non-zero damping coefficient, that is important in lightly damped system can be obtained with varying number of bodies and distances.

• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.46-53
• /
• 1997

The flow field around a three dimensional minivan-like body has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations on a non-orthogonal curvilinear coordinate system using second-order accurate schemes for the time derivatives, and third/second-order scheme for the spatial derivatives. The Marker-and-Cell concept is applied to efficiently solve continuity equation. A H-H type of multi-block grid system is generated around a three dimensional minivan-like body. Turbulent flows have been modeled by the Baldwin-Lomax turbulent model. To validate present procedure, the flows around the Ahmed body with 12.5° of slant angle are simulated. A good agreement with other numerical results is achived. After code validation, the flows around a mimivan-like body are simulated. The simulation shows three dimensional vortex-pair just behind body. The flow separation is also observed on the rear of the body. It has concluded that the results of present study properly agreed with physical flow phenomena.

• Journal of computational fluids engineering
• /
• v.17 no.1
• /
• pp.94-100
• /
• 2012

In this paper, the cavitating flows around a hydrofoil have been numerically investigated by using a 2-d multi-phase RANS flow solver based on pseudo-compressibility and a homogeneous mixture model on unstructured meshes. For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras one equation model was employed for the closure of turbulence. A dual-time stepping method and the Gauss-Seidel iteration were used for unsteady time integration. The phase change rate between the liquid and vapor phases was determined by Merkle's cavitation model based on the difference between local and vapor pressure. Steady state calculations were made for the modified NACA66 hydrofoil at several flow conditions. Good agreements were obtained between the present results and the experiment for the pressure coefficient on a hydrofoil surface. Additional calculation was made for cloud cavitation around the hydrofoil. The observation of the vapor structure, such as cavity size and shape, was made, and the flow characteristics around the cavity were analyzed. Good agreements were obtained between the present results and the experiment for the frequency and the Strouhal number of cavity oscillation.

• Proceedings of the Korean Society of Marine Engineers Conference
• /
• 2001.05a
• /
• pp.40-45
• /
• 2001

Recently, the use of numerical simulation has been increased rapidly because of the development of high performance computer systems. The present study is aimed to investigate flow characteristics of a two-dimensional sharp plane. Unsteady calculation by FDM(Finite Difference Method) based upon SOLA scheme which was performed at $Re=2{\times}10^4$in viscous incompressible flow within a finite domain on the irregular grid formation. Total numbers of irregular grids are $8{\times}10^4$. The minimum grid size is 1/100 of the plane length L which is the representative length. The inclined angles of every objects are $15^{\circ}, \;30^{\circ}\;and\; 45^{\circ}.$ And, the edge angle of the plane is $30^{\circ}.$ This study discussed the flow characteristics in term of the turbulent intensity, vorticity and frequency analysis. Developed flows show that the periodic Karman vortices occur at the back of the plane.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.38 no.9
• /
• pp.731-737
• /
• 2014

A vortex chamber is a simple device that separates compressed gas into a high-temperature stream and a low-temperature stream. It is increasing in popularity as a next-generation heat exchanger, but the flow physics associated with it is not yet well understood. In the present study, both experimental and numerical analyses were performed to investigate the temperature separation phenomenon inside the vortex chamber. Static pressures and temperatures were measured using high-sensitivity pressure transducers and thermocouples, respectively. Computational fluid dynamics was applied to simulate 3D unsteady compressible flows. The simulation results showed that the temperature separation is strongly dependent on the diameter of the vortex chamber and the supply pressure at the inlet ports, where the latter is closely related to the viscous work. The previous concept of a pressure gradient wave may not be a reasoning for temperature separation phenomenon inside the vortex chamber.

• Journal of the Society of Naval Architects of Korea
• /
• v.40 no.5
• /
• pp.36-42
• /
• 2003

A CFD simulation technique has been developed to handle the unsteady body motion with large amplitude by use of overlapping multi-block grid system. The three-dimensional, viscous and incompressible flow around body is investigated by solving the Navier-Stokes equations, and the motion of body is represented by moving effect of the grid system. Composite grid system is employed in order to deal with both the body motion with large amplitude and the condition of numerical wave maker in convenience at the same time. The governing equations, Navier-Stokes (N-S) and continuity equations, are discretized by a finite volume method, in the framework of an O-H type boundary-fitted grid system (inner grid system including test model) and a rectangular grid system (outer grid system including simulation equipments for generation of wave environments). If this study, several flow configurations, such as an oscillating cylinder with large KC number, are studied in order to predict and evaluate the hydrodynamic forces. Furthermore, the motion simulation of a Series 60 model advancing in a uniform flow under the condition of enforced roll motion of angle 20$^{\circ}$ is performed in the developed numerical wave tank.

• Journal of computational fluids engineering
• /
• v.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.