• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.32 no.3
• /
• pp.1-9
• /
• 2004

This study concentrates on the upwind schemes for convergence acceleration of the multigrid method for the Navier-Stokes equations. Comparative study of the upwind schemes in the Fourier space has been performed to identify why the second-order upwind scheme with enlarged stencil can be preconditioned better than the classical second-order upwind scheme. The full-coarsening multigrid method with implicit preconditioned multistage scheme has been implemented for verification of analysis. Numerical simulations on the inviscid and turbulent flows with the Spalart-Allmaras turbulent model have been performed. The results showed consistent trend with the analysis.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.68-74
• /
• 2014

This computational study features the transient compressible and inviscid flow analysis on a metallic plasma discharge from the opposing composite electrodes which is subjected to pulsed electric power. The computations have been performed using the flux corrected transport algorithm on the axisymmetric two-dimensional domain of electrode gap and outer space along with the calculation of plasma compositions and thermophysical properties such as plasma electrical conductivity. The mass ablation from aluminum electrode surfaces are modeled with radiative flux from plasma column experiencing intense Joule heating. The computational results shows the highly ionized and highly under-expanded supersonic plasma discharge with strong shock structure of Mach disk and blast wave propagation, which is very similar to muzzle blast or axial plasma jet flows. Also, the geometrical effects of composite electrodes are investigated to compare the amount of mass ablation and penetration depth of plasma discharge.

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

Aerodynamic shape optimization of two-dimensional airfoils in inviscid compressible flows is performed using a continuous adjoint formulation on unstructured meshes. Accurate evaluation of the gradient is achieved by using a reconstruction scheme based on the Laplacian averaging. A least-square method with extended stencil is used for flow gradient calculations. Proper convergence criterion is studied on Euler and adjoint equations for efficient design. The present method has been applied to RAE2822 and NACA0012 airfoils such that wave drag can be minimized by removing the shock wave. An inverse design is also performed to recover the shock wave on the designed RAE2822 airfoil.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.5
• /
• pp.9-17
• /
• 2005

In the early stage of helicopter design, an optimal configuration is usually determined after a numerous parametric study about the aerodynamic performance due to geometric variation. In order to improve the aerodynamic performance of a shrouded tail rotor, optimization of the tip clearance gap between blade and shroud, the blade planform shape, and the arrangement of blade spacing is required. In the present study, the aerodynamic performance characteristics of a shrouded tail rotor due to geometric variation was investigated by using an inviscid compressible unstructured mesh flow solver for rotary wings.

• Journal of ILASS-Korea
• /
• v.11 no.2
• /
• pp.89-97
• /
• 2006

An inviscid axisymmetric model capable of predicting droplet bouncing and the detailed pre-impact motion, influenced by the ambient pressure, has been developed using boundary element method (BEM). Because most droplet impact simulations of previous studies assumed that a droplet was already in contact with the impacting substrate at the simulation start, the previous simulations could not accurately describe the effect of the gas compressed between a failing droplet and the impacting substrate. To properly account for the surrounding gas effect, an effect is made to release a droplet from a certain height. High gas pressures are computationally observed in the region between the droplet and the impact surface at instances just prior to impact. The current simulation shows that the droplet retains its spherical shape when the surface tension energy is dominant over the dissipative energy. When increasing the Weber number, the droplet surface structure is highly deformed due to the appearance of the capillary waves and, consequently, a pyramidal surface structure is formed; this phenomenon was verified with our experiment. Parametric studies using our model include the pre-impact behavior which varies as a function of the Weber number and the surrounding gas pressure.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.3
• /
• pp.590-598
• /
• 1988

Viscous flows through a screen normal to an uniform flow are numerically simulated. A .kappa.-.epsilon. model is adopted for evaluation of the Reynolds stresses. The existence of a screen is regarded as extra sources in the momentum equations. The amount of extra sources is related to the resistance coefficient and the refraction coefficient of the screen. Flows are numerically simulated for various resistance coefficients and heights of the screen and Reynolds numbers. The present method has been verified to reasonably simulate viscous wakes and shear layers of the screen, for which the inviscid theory is quite limitted. As the fluids approach the screen, the velocity is reduced and the pressure is raised to satisfy the Bernoulli equation. After passing the screen, the velocity shows its minimum value at the down-stream, but static pressure is slowly recovered. A detached separation-bubble from the screen appears as the resistance coefficient is increased to a certain level. Such results are qualitatively in agreement with limitted experimental data available. The turbulent kinetic energy shows its maximum value at further down stream and decrease thereafter.

• Journal of applied mathematics & informatics
• /
• v.13 no.1_2
• /
• pp.53-65
• /
• 2003

Here we present a study of small-amplitude, shallow water waves on sloping beds. The beds considered in this analysis are linear and quadratic in nature. First we start with stating the relevant governing equations and boundary conditions for the theory of water waves. Once the complete prescription of the water-wave problem is available based on some assumptions (like inviscid, irrotational flow), we normalize it by introducing a suitable set of non-dimensional variables and then we scale the variables with respect to the amplitude parameter. This helps us to characterize the various types of approximation. In the process, a summary of equations that represent different approximations of the water-wave problem is stated. All the relevant equations are presented in rectangular Cartesian coordinates. Then we derive the equations and boundary conditions for small-amplitude and shallow water waves. Two specific types of bed are considered for our calculations. One is a bed with constant slope and the other bed has a quadratic form of surface. These are solved by using separation of variables method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.21 no.11
• /
• pp.1538-1551
• /
• 1997

Kinetic energy conservation for fixed and moving grids is examined in time-accurate finite element computation of fully unsteady inviscid flows. As numerical algorithms, fractional step method (FSM) and modified SIMPLE are used. To simulate the flow in moving grid system, arbitrary Lagrangian-Eulerian (ALE) method is adopted. In the present study, the energy conserving time integration rule for finite element algorithm is proposed and discussed schematically. It is shown that the discretization by Crank-Nicolson in time and Galerkin (central difference) in space must be used to ensure energy conservation. The developed code has been tested for a standing vortex in fixed or moving grid system, sloshing in a tank and propagation of a solitary wave, and has been shown to be a completely energy conserving algorithm.

• Journal of Mechanical Science and Technology
• /
• v.18 no.2
• /
• pp.313-324
• /
• 2004

A rich phenomenon in the dynamics of azimuthal vortices in a circular cylinder caused by the inertial oscillation is investigated numerically at high Reynolds numbers and moderate Rossby numbers. In the actual spin-up flow where both the Ekman circulation and the bottom friction effects are included, the first appearance of a seed vortex is generated by the Ekman boundary-layer on the bottom wall and the subsequent roll-up near the corner bounded by the side wall. The existence of the small vortex then rapidly propagates toward the inviscid region and induces a complicated pattern in the distribution of azimuthal vorticity, i.e. inertial oscillation. The inertial oscillation however does not deteriorate the classical Ekman-pumping model in the time scale larger than that of the oscillatory motion. Motions of single vortex and a pair of vortices are further investigated under a slip boundary-condition on the solid walls. For the case of single vortex, repeated change of the vorticity sign is observed together with typical propagation of inertial waves. For the case of a pair of vortices with a two-step profile in the initial azimuthal velocity, the vortices' movement toward the outer region is resisted by the crescent-shape vortices surrounding the pair. After touching the border between the core and outer regions, the pair vortices weaken very fast.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.167-173
• /
• 1999

An Euler-based CFD tool has been developed for the performance evaluation of a thrust reverser mounted on a high bypass ratio turbofan engine. The computational domain surrounded by the ground and non-reflection boundary includes the whole nacelle configuration with a deployed thrust reverser. The numerical algorithm is based on the modified Godunovs scheme to allow the second order accuracy in both space and time. The grid system is generated by using eleven multi-blocks, of which the total cell number is 148,400. The thrust reverser is modeled as if it locates at the nacelle simply in all circumferential direction. The existence of a fan and an OGV(Outlet Guide Vane) is simulated by adopting the actuator disk concept, in which predetermined radial distributions of stagnation pressure ratio and adiabatic efficiency coefficient are used for the rotor type disk, and stagnation pressure losses and flow outlet angles for the stator type disk. All boundary conditions including the fan and OGV simulation are treated by Riemann solver. The developed solver is applied to a turbofan engine with a bypass ratio of about 5.7 and the diameter of the fan cowl of 83 inch. The computational results show that the Euler-based inviscid method is very useful and economical to evaluate the performance of thrust reversers.