• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.3
• /
• pp.211-219
• /
• 2008

A three-dimensional inviscid flow solver has been developed based on unstructured meshes for the investigation of the effect of the external stores on the tail surfaces of a generic fighter aircraft. The numerical method is based on a vertex-centered finite-volume discretization and an implicit point Gauss-Seidel time integration. The calculations were made for a steady flow and the computed results were compared with experimental data to validate the flow solver. An unsteady time-accurate computation of the generic fighter aircraft with external stores at transonic flight conditions showed that the external stores cause unsteady loading on the horizontal tail surface due to the mutual interference between their wake and the horizontal tail surface. It was shown that downward deflection of the trailing edge flap significantly reduces the undesirable interference effect.

• Ocean Systems Engineering
• /
• v.3 no.3
• /
• pp.219-236
• /
• 2013

The accurate prediction of motion in waves of a marine vehicle is essential to assess the maximum sea state vs. operational requirements. This is particularly true for small crafts, such as Autonomous Surface Vessels (ASV). Two different numerical methods to predict motions of a SWATH-ASV are considered: an inviscid strip theory initially developed at MIT for catamarans and then adapted for SWATHs and new a hybrid strip theory, based on the numerical solution of the radiation forces by an unsteady viscous, non-linear free surface flow solver. Motion predictions obtained by the viscous flow method are critically discussed against those obtained by potential flow strip theory. Effects of viscosity are analyzed by comparison of sectional added mass and damping calculated at different frequencies and for different sections, RAOs and motions response in irregular waves at zero speed. Some relevant conclusions can be drawn from this study: influence of viscosity is definitely non negligible for SWATH vessels like the one presented: amplitude of the pitch and heave motions predicted at the resonance frequency differ of 20% respectively and 50%; in this respect, the hybrid method with fully non-linear, viscous free surface calculation of the radiation forces turns out to be a very valuable tool to improve the accuracy of traditional strip theories, without the burden of long computational times requested by fully viscous time domain three dimensional simulations.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.17 no.11
• /
• pp.363-367
• /
• 2016

A test section of a supersonic wind tunnel was designed for the analysis of flow characteristics over a backward-facing step with Mach 1.0 slot injection in a supersonic flow of Mach 2.5. The cavity flow of a high-speed vehicle is very complex at supersonic speed, so it is necessary to do experiments using supersonic wind tunnels to verify numerical analysis methods. The previous 2D symmetrical nozzle was replaced with an asymmetrical nozzle. The inviscid nozzle contour was designed using Method of Characteristics (MOC), and the boundary layer thickness correction was reflected by experimental data from the wind tunnel. The results were compared with a CFD analysis. The PID control system was changed to be based on the change of tank pressure. This improved the control efficiency, and the run times of supersonic flow increased by about 1 second. The flow characteristics over a backward facing step with slot injection were visualized by a Schlieren device. This equipment will be used for an experimental study of the film cooling effectiveness over a cavity with various velocities, mass flows, and temperatures.

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

Wave-body interaction is simulated using a developed code based on the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. A free surface is captured as a moving contact discontinuity within a fluid domain and an approximated Riemann solver is used to estimate the inviscid flux across the discontinuity. Immersed boundary nodes are identified inside an instantaneous fluid domain near a moving body, then dependent variables are reconstructed at those immersed boundary nodes based on interpolation along local normal lines to the boundary. Free surface flows around an oscillating cylinder are simulated and the computed wave elevations are compared with other reported results. The generation of a solitary wave by a moving wave-maker is simulated and the time histories of wave elevations at two different points are compared with other results. The developed code is applied to simulate body motion of an elastically mounted circular cylinder as a solitary wave passes the body. The force acting on an elastically mounted cylinder is compared with the force acting on a fixed cylinder. Grid independency of the computed body motion is established based on a comparison of results using three different-size grids.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.4 no.2
• /
• pp.21-30
• /
• 2000

Numerical design optimization techniques are implemented for the improvement of the ram accelerator performance. The design object is to find the minimum ram tube length required to accelerate projectile from initial velocity $V_o$ to target velocity $V_e$. The premixture is composed of $H_2$, $O_2$, $N_2$ and the mole numbers of these species are selected as design variables. The objective function and the constraints are linearized during the optimization process and gradient-based Simplex method and SLP(Sequential Linear Programming) have been employed. With the assumption of two dimensional inviscid flow for internal flow field, the analyses of the nonequilibrium chemical reactions for 8 steps 7 species have been performed. To determined the tube length, ram tube internal flow field is assumed to be in a quasi-steady state and the flow velocity is divided into several subregions with equal interval. Hence the thrust coefficients and accelerations for corresponding subregions are obtained and integrated for the whole velocity region. With the proposed design optimization techniques, the total ram tube length had been reduced 19% within 7 design iterations. This optimization procedure can be directly applied to the multi-stage, multi-premixture ram accelerator design optimization problems.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.5
• /
• pp.1033-1042
• /
• 1994

Dynamic behavior of flexible rectangular liquid containers is analyzed by a two-dimensional coupled boundary element-finite element method. The irrotational motion of inviscid and incompressible ideal fluid is modeled by boundary elements and the motion of structure by finite elements. A singularity free integral formulation is employed for the implementation of boundary element method. Coupling is performed by using compatibility and equilibrium conditions along the interface between the fluid and structure. The fluid-structure interaction effects are reflected into the coupled equation of motion as added fluid mass matrix and sloshing stiffness matrix. By solving the eigen-problem for the coupled equation of motion, natural frequencies and mode shapes of coupled system are obtained. The free surface sloshing motion and hydrodynamic pressure developed in a flexible rectangular container due to horizontal and vertical ground motions are computed in time domain.

• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). 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 model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.1
• /
• pp.1-8
• /
• 2001

The computations of the turbulent flow around the ship models with the free-surface effects were carried out. Incompressible Reynolds-Averaged Navier-Stokes equations were solved by using an explicit finite-difference method with the nonstaggered grid system. The method employed second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration. For the turbulence closure, a modified Baldwin-Lomax model was exploited. The location of the free surface was determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and a free-surface conforming grid was generated at each time step so that one of the grid boundary surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition was applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method developed in the present study, the computations were carried out for beth Wigley and Series 60 $C_B=0.6$ ship model and the computational results showed good agreements with the experimental data.

• 한국해양학회지
• /
• v.28 no.3
• /
• pp.177-185
• /
• 1993

The theory for long wave scattering (Proudman, 1914: Longuet-Higgins, 1970) is applied to a tidal-frequency Kelvin wave propagating around a small cylindrical island in a shelf sea of uniform depth. The theory includes the effects of bottom friction on wave propagation. The theoretical analysis of the Kelvin wave around the island. this amplitude change results in a uniform amplitude of the total wave along the circumference of the island in an inviscid fluid, and the dynamic cause of this is explained in terms of Coriolis effects. Bottom friction attenuates the amplitude of the total wave from the frontal side of the island to the leeward side, but the amplitude variation along the coast becomes symmetric to the line connecting both idea. The phase of the scattered wave contributes to more rapid travel of the total wave in the front and leeward side than farther offshore. The effects of bottom friction on the wave phase around the island are negligible.

• Journal of computational fluids engineering
• /
• v.18 no.1
• /
• pp.83-90
• /
• 2013

OpenMP is becoming more and more useful as a simple parallel processing paradigm on SMP (Shared Memory Multi-Processors) computing environment with the development of multi-core processors. However, very few data is available publically regarding the OpenMP performance in CFD (Computational Fluid Dynamics). In the present study a CFD test suite is prepared for the performance evaluation of OpenMP on various multi-core systems. The test suite is composed of two-dimensional numerical simulations for inviscid/viscous and reacting/non-reacting flows using three different levels of grid systems. One to five test runs were carried out on various systems from dual-core dual threads to 16-core 32-threads systems by changing the number of threads engaged for each test up to 80. The results exhibit some interesting results and the lessons learned from the tests would be quite helpful for the further use of OpenMP for CFD studies using multi-core processor systems.