• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• Journal of the Korean Society of Industry Convergence
• /
• v.5 no.1
• /
• pp.45-52
• /
• 2002

In this study, an exact solution of governing differential equation for the bending problem of partially loaded orthotropic rectangular plates is presented and also its computer program is developed. The method requires that two opposite edges be clamped or simply supported, or one edge clamped and the other simply supported. Any combination of boundary conditions could exist along the other edges. The plate could he subjected to uniform, partially uniform, and line loads. The solution for the deflection of rectangular plate is expressed as a Levy type single Fourier series and the loads arc expressed as a corresponding series. The advantage of the solution is that it overcomes the limitations of the previous Navier's and Levy's methods (limitation of boundary condition and loading conditions of plate), it is easy to program on a computer and it becomes fast to solve the bending problem with computer program. Calculations are presented for isotropic and orthotropic plates with different loading and boundary conditions. Comparisons are made for the isotropic plate with various boundary conditions between the result of this paper and the result of Navier, Levy and Szilard. The deflections were in excellent agreement.

• 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 the Korean Society of Propulsion Engineers
• /
• v.15 no.5
• /
• pp.10-18
• /
• 2011

A new film cooling method to protect the pylon injector from aerodynamic heating for a scramjet combustor is proposed and verified with numerical methods. The conditions for the Mach 8 flight at an altitude of 35km are considered. Air is considered as a coolant. Three-dimensional Navier-Stokes equations with $k-{\omega}$ SST turbulence model are used. A downward injection of coolant from the top of the pylon gives higher cooling effects with less mass flow rate of coolant than the upward coolant injection from bottom of the pylon. Also, the downward injection shows little flow separation due to the favorable pressure gradient and does not disturb the flowfields near pylon injector, which results in reduction of pressure losses.

• The KSFM Journal of Fluid Machinery
• /
• v.19 no.1
• /
• pp.37-44
• /
• 2016

In this study, a parametric study of a centrifugal pump with double volute has been performed numerically using three-dimensional Reynolds-averaged Navier-Stokes equations. The shear stress transport model was selected as turbulence closure through turbulence model test. The finite volume method and unstructured grid system were used for the numerical analysis. The optimal grid system in the computational domain was determined through a grid dependency test. The expansion coefficient, circumferential and radial starting positions and length of divider were selected as the geometric parameters to be tested. And, the hydraulic efficiency and the radial thrust coefficient were considered as performance parameters. It was found that the radial thrust and hydrualic efficiency are more sensitive to the expansion angle and circumferential starting position of the divider than the other geometrical parameters.

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

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.1
• /
• pp.49-56
• /
• 1997

Heat-transfer enhancement is seeked through modifications of fin surface. Real life plate-fin heat exchangers have complex three-dimensional geometries. Fins can have arrays of dimples and are attached to rows of penetrating tubes. To isolate the effect of surface modification, we model the real flow by a two-dimensional channel flow with a dimple on one side. The flow is analysed by solving the incompressible Navier-Stokes equation by a finite volume method on a generalized boundary-fitted coordinate. Results show a trapped vortex inside the dimple for all cases computed. Local maximum of Nusselt number occurs near the downstream end of the dimple, due to such a vortex. Location of the vortex does not change with respect to the wall temperature change, but moved downstream when Reynolds number increases. This, together with the results that in all cases vortex core is somewhat downstream of the dimple center, suggests that the mean flow above continuously feeds the kinetic energy to the recirculating flow. Heat transfer enhancement and pressure losses are studied through analysing the relevant dimensionless parameters like, Nusselt number and friction factor. In all cases computed, dimpled channel flow experiences less pressure loss than two-dimensional Poiseuille flow.

• Journal of Advanced Marine Engineering and Technology
• /
• v.35 no.6
• /
• pp.809-813
• /
• 2011

Turbulent flows are numerically simulated in the three dimensional inlet duct for heat recovery steam generator. The present study is aimed to analyze the effect of a variation in turbulent flow pattern by the change of roof angle in the transition duct. The finite volume based Navier-Stokes equations with unstructured grids are solved to make clear the flow dynamic phenomena. Reviews are made on with the data of path lines, velocity vectors, dynamic pressure, residuals for numerical convergence and so on. The k-epsilon, k-omega, Reynolds stress and RNG k-epsilon are used for generation of turbulence. Two types of roof angle are applied with and without the swirl in the duct. Turbulent flow patterns could be investigated for the optimum duct design based on the computational results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.8
• /
• pp.663-671
• /
• 2007

In this work, shape optimization of a wire-wrapped fuel assembly in a liquid metal reactor has been carried out by combining a three-dimensional Reynolds-averaged Navier-Stokes analysis with the radial basis neural network method, a well known surrogate modeling technique for optimization. Sequential Quadratic Programming is used to search the optimal point from the constructed surrogate. Two geometric design variables are selected for the optimization and design space is sampled using Latin Hypercube Sampling. The optimization problem has been defined as a maximization of the objective function, which is as a linear combination of heat transfer and friction loss related terms with a weighing factor. The objective function value is more sensitive to the ratio of the wire spacer diameter to the fuel rod diameter than to the ratio of the wire wrap pitch to the fuel rod diameter. The optimal values of the design variables are obtained by varying the weighting factor.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.4
• /
• pp.44-52
• /
• 2002

Aerodynamic characteristics of a flapping airfoil in low Reynolds number flows are numerically studied using the unsteady, incompressible Navier-Stokes flow solver with a two-equation turbulence model. For more efficient computation of unsteady flows over flapping airfoil, the flow solver is parallel-implemented by MPI programming method Unsteady computations are performed for low Reynolds number flows over a NACA four-digit series airfoils. Effects of pitching, plunging, and flapping motion with different reduced frequency, amplitude, thickness and camber on aerodynamic characteristics are investigated. Present computational results yield a better agreement in thrust at various reduced frequency with experimental data.