• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.

• Journal of computational fluids engineering
• /
• v.8 no.2
• /
• pp.24-32
• /
• 2003

In this work, we propose a new idea of flapping airfoil design for optimal aerodynamic performance from detailed computational investigations of flow physics. Generally, flapping motion which is combined with pitching and plunging motion of airfoil, leads to complex flow features such as leading edge separation and vortex street. As it is well known, the mechanism of thrust generation of flapping airfoil is based on inverse Karman-vortex street. This vortex street induces jet-like flow field at the rear region of trailing edge and then generates thrust. The leading edge separation vortex can also play an important role with its aerodynamic performances. The flapping airfoil introduces an alternative propulsive way instead of the current inefficient propulsive system such as a propeller in the low Reynolds number flow. Thrust coefficient and propulsive efficiency are the two major parameters in the design of flapping airfoil as propulsive system. Through numerous computations, we found the specific physical flow phenomenon which governed the aerodynamic characteristics in flapping airfoil. Based on this physical insight, we could come up with a new kind of airfoil of tadpole-shaped and more enhanced aerodynamic performance.

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

This paper presents the numerical results of fluid flow and mixing in a microfluidic device for electro-osmotic flow (EOF) with an trapezoidal electrode array on the bottom wall (ETZEA). Differently from previous EOF in a channel which only transports fluid in colloidal system. ETZEA can also be utilized to mix a target liquid with a reagent. In this study we propose a method of controlling fluid flow and mixing enhancement. To obtain the flow and mixing characteristics, numerical computations are performed by using a commercial code, CFX-10, and a self-made code LBM-D. It was found that the flow near the trapezoidal electrode in the ETZEA is of 3-D complex flows due to the zeta potential difference between the trapezoidal electrode and channel walls, and as a consequence the hetrogeneous zeta potential on the electrodes plays an important role in mixing the liquid.

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

Flow around a guided missile in high maneuver, i.e. at a high angle of attack, shows complex phenomena. It is well known that even in geometrically symmetric conditions the flow around a missile at high angles of attack can generate unexpected large side forces and yaw moments due to asymmetric vortices. In this paper, a CFD code (FLUENT) based on the Navier-Stokes equations was used for the numerical analysis to find a suitable numerical mechanism for generation of asymmetric vortices. It is shown that a numerical technique of applying different surface roughness to a specific area of the missile nose surface gives the best fit in comparison with the experimental results. In addition, a numerical investigation of variations of side forces and pressure distributions with angle of attack and roll angle was conducted for the purpose of identifying the source of vortex asymmetries.

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

NSBT(No Slip Boundary Treatment) is a newly developed scheme for the treatment of a no slip condition on the solid wall of obstacle in a flow field. In our research, NSBT was used to perform LBM simulation of a flow over a circular cylinder to determine the flow feature and aerodynamics characteristic of the cylinder. To ascertain the applicability of NSBT on the complex shape of the obstacle, it was first simulated for the case of the flow over a circular and square cylinder in a channel and the results were compared against the solution of Navier-Stokes equation. The simulations were performed in a moderate range of Reynolds number at each cylinder position to identify the flow feature and aerodynamic characteristics of circular cylinder in a channel. The drag coefficients of the cylinder were calculated from the simulation results. We have numerically confirmed that the critical reynolds number for vortex shedding is in the range of 200$\sim$250. For the gap parameter $\gamma$ = 2 cases at Re > 240, the vortex shedding were symmetric and it resembled the Karmann vortex. As the cylinder approached to one wall, the vorticity significantly reduced in length while the vorticity on the other side elongated and the vorticity combined with the wall boundary-layer vorticity. The resultant $C_d$ by LBM concurred with the results of DNS simulation performed by previous researchers.

• Journal of computational fluids engineering
• /
• v.13 no.1
• /
• pp.35-40
• /
• 2008

Generally flight vehicles have many cavities such as wheel wells, bomb bays and windows on their external surfaces and the flow around these cavities makes separation, vortex, shock and expansion waves, reattachment and other complex flow phenomenon. The flow around the cavity makes abnormal and three-dimensional noise and vibration even thought the aspect ratio (L/D) is small. The cavity giving large effects to the flow might make large noise, cause structural damage or breakage, harm the aerodynamic performance and stability, or damage the sensitive devices. In this study, numerical analysis was performed for cavity flows by the unsteady compressible three dimensional Reynolds-Averaged Navier-Stokes (RANS) equations with Wilcox's $\kappa-\omega$ turbulence model. The MPI(Message Passing Interface) parallelized code was used for calculations by PC-cluster. The cavity has the aspect ratios of 2.5, 3.5 and 4.5 with the W/D ratio of 2 for three-dimensional cavities. The Sound Pressure Level (SPL) analysis was done with FFT to check the dominant frequency of the cavity flow. The dominant frequencies were analyzed and compared with the results of Rossiter's formula and Ahuja& Mendoza's experimental datum.

• Journal of computational fluids engineering
• /
• v.15 no.1
• /
• pp.81-87
• /
• 2010

In this study, we conducted resistance and propulsion performance test of ship composed of the Resistance Test, Propeller Open Water Test and Self Propulsion Test using the CFD(Computational Fluid Dynamics). We used commercial RANS(Reynolds Averaged Navier Stokes equation) solver, as a calculating tool. The unstructured grids were used in a bow and stern of ship, having complex shape, for a convenience of generating grids, and the structured grids were adopted in a central hull and rest of hull having a relatively simple shape which is called hybrid grid method. In addition, The sliding mesh method was adopted to rotate a propeller directly in the Propeller Open Water and Self Propulsion Test. The Resistance Test and Self Propulsion Test were calculated using Volume of Fluid (VOF) model and considering a free surface. And all The three cases were applied realizable k-epsilon model as the turbulence model. The results of calculations were verified for the suitability of calculations by comparing MOERI's EFD results.

• Journal of computational fluids engineering
• /
• v.17 no.2
• /
• pp.11-20
• /
• 2012

In this paper, the authors test various open codes and commercial codes based on CFD technology on the wind field around the complex terrain, which is a very important transport physics in the event of forrest fire. To study the physical mechanism inside the transition from surface fire to crown fire, the wake flow behind a parallel array of trees is studied numerically to show the flow separation in the turbulent boundary layer. Two sites near to Kunsan National University are chosen for the measurement of real wind field, and obtained data are compared with those from various computational codes such as Wind-Ninja, NIST-FDS, ANSYS-CFX, and ANSYS-FlUENT, etc. Through this research, feasibility and accuracy of the present CFD codes are investigated quantitatively, compared with the measured data with AWS.

• Journal of computational fluids engineering
• /
• v.6 no.3
• /
• pp.32-40
• /
• 2001

Viscous solutions of supersonic jet impinging on a flat plate normal to the flow are simulated using three-dimensional Navier-Stokes solver. The jet impinging flow structure exhibits such complex nature as shock shell, plate shock and Mach disk depending on the flow parameters. Among others, the dominant parameters are the ratio of the nozzle exit pressure to the ambient pressure and the distance between the nozzle exit plane and the impinging plane. In the present study, the nozzle contour and the pressure ratio are held fixed, while the jet impinging distance is varied to illuminate the characteristics of the jet plume with the distance. As the plate is placed close to the nozzle at 3D high, the computed wall pressure at or near the jet center oscillates with large amplitude with respect to the mean value. Here D is the nozzle exit diameter. The amplitude of wall pressure fluctuations subsides as the distance increases, but the maximum mean pressure level at the plate is achieved when the distance is about 4D high. The frequency of the wall pressure is estimated at 6.0 kHz, 9.3 kHz, and 10.0 kHz as the impinging distance varies from 3D, 4D, to 6D, respectively.

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.55-63
• /
• 2000

In this paper, the numerical study has been done for the improvement of the superdetonative ram accelerator performance and for the design optimization of the system. The objective function to optimize the premixture composition is the ram tube length, required to accelerate projectile from initial velocity V/sub 0/ to target velocity V/sub e/. The premixture is composed of H₂, O₂, N₂ and the mole numbers of these species are selected as design variables. RSM(Response Surface Methodology) which is widely used for the complex optimization problems is selected as the optimization technique. In particular, to improve the non-linearity of the response and to consider the accuracy and the efficiency of the solution, design space stretching technique has been applied. Separate sub-optimization routine is introduced to determine the stretching position and clustering parameters which construct the optimum regression model. Two step optimization technique has been applied to obtain the optimal system. With the application of stretching technique, we can perform system optimization with a small number of experimental points, and construct precise regression model for highly non-linear domain. The error compared with analysis result is only 0.01% and it is demonstrated that present method can be applied to more practical design optimization problems with many design variables.