• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.36 no.10
• /
• pp.939-946
• /
• 2008

Euler codes or Navier-Stokes codes for compressible flows suffer severe degradation in convergence as Mach number approaches zero. The convergence problem arose from the wide disparity in characteristic speeds can be solved using preconditioning methods without large modifications. In this paper, a preconditioned RANS(Reynolds Averaged Navier-Stokes) solver is developed for analysis of low Mach number flows. In order to validate the method, computational examples are chosen and the results are compared with the experimental data and the existing computed results showing a good accuracy and convergence characteristics for steady inviscid, laminar and turbulent flows at low Mach number.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.3 no.2
• /
• pp.18-24
• /
• 2000

A two-dimensional numerical method for inviscid two-fluid flows with a significant entrainment into both directions is established, and the oil leakage from a non-pressurized underwater pipe is studied. The interface between two fluids is modeled at a vortex sheet. The flow field and the subsequent interface evolution are solved by using the vortex-in-cell method. For longer flow simulation with a realistic two fluids interaction, an efficient merging scheme is introduced. In the Boussinesq limit, the speed of the external fluid intrusion into the pipe is very close to the existing mathematical models, and the lock exchange is observed in spite of a significant roll-up of the interface and entrainments. It is believed that the developed method can be utilized effectively for further detailed studies on various two-fluid flows which are encountered in many different marine oil spill problems.

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

• International Journal of Aeronautical and Space Sciences
• /
• v.17 no.2
• /
• pp.260-267
• /
• 2016

The design process of blended inlet body (BIB) for the preliminary design of a near-space high supersonic unmanned aerial vehicle (HSUAV) is presented. The mass flow rate and cowl area of inlet at a design point are obtained according to the cruise condition of the HSUAV. A mixed-compression axisymmetric supersonic inlet section with a fixed geometry reasonably matching the high supersonic cruise state is created by using the inviscid theory of aerodynamics. The inlet section is optimized and used as a baseline section for the BIB design. Three BIB concepts for the HSUAV are proposed, and their internal aerodynamic characteristics of inlet are evaluated using Euler computational fluid dynamics (Euler CFD) solver. The preferred concept is identified, in which the straight leading edge of the baseline HSUAV configuration is modified into the convex leading edge to accommodate the inlet and meet the requirements of the cowl area to capture the sufficient air flow. The total recovery of inlet for the preferred BIB concept and the aerodynamic characteristics of the modified HSUAV configuration are verified using Navier-Stokes computational fluid dynamics (NS CFD) solver. The validation indicates that the preferred BIB concept can meet both the requirements of the inlet and aerodynamic performance of the HSUAV.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.37 no.11
• /
• pp.1055-1065
• /
• 2009

In this study, stability analysis of boundary layers on airfoils is performed by using parabolized stability equations(PSE). Boundary layer edge conditions are obtained by compressible inviscid flow calculations. Mean velocity and temperature profiles of the laminar boundary layer are obtained by solving compressible boundary layer equations in generalized curvilinear coordinates with fourth order accuracy in the wall normal direction. Laminar mean flow profiles are used as input data for PSE to investigate growth rates of disturbances and stability characteristics. For the cases of boundary layer on NACA0012 and HSNLF(1)-0213 airfoils at Mach number 0.5, growth rates with respect to disturbance frequencies and profiles of disturbance amplitude are investigated. The effect of angle of attack on stability characteristics are examined at both upper and lower surfaces. The neutral stability curves, effect of Mach number and effect of airfoil section shapes are also analyzed.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.3
• /
• pp.80-86
• /
• 1994

A numerical method for simulations of inviscid incompressible flow fields around a ship advancing on the free surface is developed. A body fitted coordinate system, generated by numerically solving elliptic type partial differential equations is used to conform the ship and free surface configurations. Three dimensional Euler equations transformed to the non-staggered body fitted coordinate system are discretised by finite difference method. Time and spatial derivatives are discretised by forward and centered differencings, respectively, and artificial dissipations are added to discretised convection terms for improvements of numerical stability. At each time steps, free surface elevations are recomputed to satisfy nonlinear free surface conditions. Poisson equations for pressure field are solved iteratively and the velocity field for next time step is extrapolated. To verify the developed numerical method, flow fields around a Wigley model are simulated(Fn=0.250-0.408) and compared with experimental data to show good agreements.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.3
• /
• pp.161-167
• /
• 1992

We consider the preparation of long finite amplitude nondispersive waves over a step bottom between two regions of finite different depths. Two dimensional motion is assumed. with the wave crests parallel to the step, and irrotational flow in the inviscid fluid is considered. To describe the transformation of finite amplitude waves we use the finite-amplitude shallow-water equations, the conditions of mass flow conservation and pressure continuity at the cut above the step in Riemann's variables. The equations define four families of curves-characteristics on which the values of the Riemann's invariants remain constant and a system of two nonlinear equations that relates the amplitudes of incident reflected and transmitted waves. The system obtained is difficult to analyze in common form. Thus we consider some special cases having practical usage for tsunami waves. The results obtained are compared with the long wave theory and significant nonlinear effects are found even for quite small amplitude waves.

• Journal of the Korean earth science society
• /
• v.22 no.6
• /
• pp.512-527
• /
• 2001

Separated flows passed complex geometries are modeled by discrete vortex techniques. The flows are assumed to be rotational and inviscid, and a new techlnique is described to determine the stream functions for linear shear profiles. The geometries considered are the snow cornice and the backward-facing step, whose edges allow for the separation of the flow and reattachment downstream of the recirculation regions. A point vortex has been added to the flows in order to constrain the separation points to be located at the edges, while the conformal mappings have been modified in order to smooth the sharp edges and to let the separation points free to oscillate around the points of maximum curvature. Unsteadiness is imposed to the flow by perturbing the vortex location, either by displacing the vortex from the equilibrium, or by imposing a random perturbation with zero mean to the vortex in equilibrium. The trajectories of passive scalars continuously released upwind of the separation point and trapped by the recirculating bubble are numerically integrated, and concentration time series are calculated at fixed locations downwind of the reattachment points. This model proves to be capable of reproducing the trapping and intermittent release of scalars, in agreement with the simulation of the flow passed a snow cornice performed by a discrete multi-vortex model, as well as with direct numerical simulations of the flow passed a backward-facing step. The results of simulation indicate that for flows undergoing separation and reattachment the unsteadiness of the recirculating bubble is the main mechanism responsible for the intense large-scale concentration fluctuations downstream.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.8
• /
• pp.37-46
• /
• 2004

Several notes on ground effects drawn from Navier-Stokes analyses and their aerodynamic interpretations were addressed here; For two-dimensional ground effect, the change of surface pressure due to image vortex, the venturi effect due to thickness and the primary inviscid flow phenomena of ground effect, and for three-dimensional ground effect, strengthened wing tip vortices, increased effective span and the outward drift of trailing vortices. Irodov's criteria were evaluated to investigate the static longitudinal stability of conventional NACA 6409 and DHMTU 8-30 airfoils. The analysis results demonstrated superior static longitudinal stability of DHMTU 8-30 airfoil. The DHMTU airfoil has quite lower value of lrodov's criterion than the conventional NACA airfoil, which require much smaller tail volume to stabilize the whole WIG-craft at its design stage.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.24 no.5
• /
• pp.65-76
• /
• 2020

On the initial design process of a scramjet vehicle such as the trajectory prediction, it is inevitable to estimate the aerodynamic performance of a three-dimensional effect. Despite the necessity of intensive computing for the three-dimensional model, it is inefficient in predicting a wide range of aerodynamic performance. In this study, an engineering model for aerodynamic performance was developed based on two-dimensional computational fluid analysis and linearized supersonic inviscid flow theory. Correspondingly, the three-dimension aerodynamic performance relations are presented based on the two-dimensional results. And the additional three-dimensional computation was performed to evaluate the adequacy for the extended relations.