• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.40 no.12
• /
• pp.1025-1031
• /
• 2012

The present paper deals with concept design of supersonic inlet based on compressible flow theory and flow control of bleeding in order to guarantee stability of supersonic inlet of ramjet engine in broad range of operating conditions. Shock instability, shock wave-boundary layer interaction and flow separation should be properly controlled to improve performance of the supersonic inlet. Considering shock strength, boundary layer and flow separation, the supersonic inlet is modified from the basic model which is designed under inviscid theory. Consequently, shock is stabilized, and required mass flow rate is obtained. Furthermore, bleeding is applied to the supersonic inlet to maintain performance in off-design conditions. Mass flow condition is adopted for modeling of bleeding effect, and performance of the supersonic inlet is evaluated by changing bleeding locations and numbers.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.4
• /
• pp.12-17
• /
• 2018

A natural cavitation-type high-speed underwater vehicle with solid rocket motor is tested, and its speed and running distance are measured. The outputs from pressure sensors on the surface of the vehicle reveal a pressure-time history reflecting the development of supercavitation. Underwater cameras installed on the wall of the test pool record the entire process from the onset of supercavitation to its full development. CNU-SuperCT, based on two-dimensional inviscid theoretical analysis, is used to simulate test results. Considering CNU-SuperCT does not include the control fins of the vehicle, simulation results agree with test results very well. Additionally, pictures from underwater cameras support the test results.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.3
• /
• pp.29-40
• /
• 1993

The fundamental characteristics of nonlinear free-surface waves generated by a shallowly submerged 3-dimensional hydrofoil are investigated. The fluid is assumed inviscid, incompressible and its motion irrotational. The surface tension on the free-surface is neglected. The hydrofoil is represented by a horseshoe vortex system whose shape is assumed fixed. Also the strengths of vortices are assumed given. The exact problem for the wave potential due to the horseshoe vortex system is formulated by the variational principle based on the classical Hamilton's principle. The localized finite element method is used in the numerical computations. In order to increase the numerical efficiency, an intermediate nonlinear-to-linear transition buffer subdomain for a smooth matching is introduced between the fully nonlinear computation subdomain and the truncated linear infinite subdomain. Also used is the modal analysis to reduce the computation tome drastically. The effect of inflow velocity, submergence depth of the hydrofoil and the shape of circulation distribution on the wave profiles are thoroughly examined. Especially it was possible to investigate the nonlinear influence of the free vortex on the free vortex. The nonlinear free-surface effect on the induced forces on the hydrofoil is also investigated.

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

• Ocean Systems Engineering
• /
• v.8 no.4
• /
• pp.381-407
• /
• 2018

This paper presents a hybrid numerical approach, which combines a two-phase Navier-Stokes model (NS) and the fully nonlinear potential theory (FNPT), for modelling wave-structure interaction. The former governs the computational domain near the structure, where the viscous and turbulent effects are significant, and is solved by OpenFOAM/InterDyMFoam which utilising the finite volume method (FVM) with a Volume of Fluid (VOF) for the phase identification. The latter covers the rest of the domain, where the fluid may be considered as incompressible, inviscid and irrotational, and solved by using the Quasi Arbitrary Lagrangian-Eulerian finite element method (QALE-FEM). These two models are weakly coupled using a zonal (spatially hierarchical) approach. Considering the inconsistence of the solutions at the boundaries between two different sub-domains governed by two fundamentally different models, a relaxation (transitional) zone is introduced, where the velocity, pressure and surface elevations are taken as the weighted summation of the solutions by two models. In order to tackle the challenges associated and maximise the computational efficiency, further developments of the QALE-FEM have been made. These include the derivation of an arbitrary Lagrangian-Eulerian FNPT and application of a robust gradient calculation scheme for estimating the velocity. The present hybrid model is applied to the numerical simulation of a fixed horizontal cylinder subjected to a unidirectional wave with or without following current. The convergence property, the optimisation of the relaxation zone, the accuracy and the computational efficiency are discussed. Although the idea of the weakly coupling using the zonal approach is not new, the present hybrid model is the first one to couple the QALE-FEM with OpenFOAM solver and/or to be applied to numerical simulate the wave-structure interaction with presence of current.

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

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.49 no.6
• /
• pp.449-456
• /
• 2021

In this paper, we present a strategy to extend solution capability of an existing low Mach number preconditioned compressible solver to incompressible flows with a little modification. To this end, the energy equation that is of the same form of the total energy equation of compressible flows is used. The energy equation is obtained by a linear combination of the thermal energy equation, the continuity equation and the mechanical energy equation. Subsequently, a modified artificial compressibility method in conjunction with a time marching technique is applied to these incompressible governing equations for steady flow solutions. It is found that the Roe average of the common governing equations is equally valid for both the compressible and incompressible flow conditions. The extension of an existing compressible solver to incompressible flows does not affect the original compressible flow analysis. Validity for incompressible flow analysis of the extended solver is examined for various inviscid, laminar and turbulent flows.

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

• Ocean Systems Engineering
• /
• v.6 no.3
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• Journal of the Korean Society for Marine Environment & Energy
• /
• v.7 no.4
• /
• pp.192-198
• /
• 2004

A numerical method is developed to solve a two-dimensional diffraction problem for a body located in a sediment pocket where a heavier muddy water is trapped. In the present study, the wave exciting forces acting on a submerged body in the water-sediment interface by an incident wave is investigate. It is assumed that the heavier mud is trapped locally in a sediment pocket. A mathematical formulation is made in the scope of the potential theory. The fluid is assumed to be inviscid, incompressible and its motion irrotational. The boundary conditions on the unknown free surface and interface are linearized. As a method of solution, the localized finite-element method is adopted. In the method, the computation domain is reduced by utilizing the complete set of analytic solutions known in the infinite subdomain to be truncated by introduction of an appropriate juncture conditions. The main advantage of this method is that any complex geometry of the boundaries can be easily accommodated. Computations are carried out for mono-chromatic plane progressive surface waves normally incident on the domain. Numerical results are compared with those obtained by Lassiter based on Schwingers variational method. Good Agreements are obtained in general. Another numerical computations are made for the cases with and without a body in the sediment pocket.