• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.4
• /
• pp.241-249
• /
• 1991

A quasi-three dimensional calculation method is presented on the basis of Wu's idea using finite element methods. In B-B flow the governing equations are cast into a single equation to overcome the restriction of the type of turbomachinery, and Kutta condition is exactly assured by introducing a combination of two kinds of stream functions. In H-S flow a dissipative force which is assumed to be opposed to the relative velocity is added to the governing equation for a consistent loss model. The entropy change along each streamline is then calculated by assuming that the dissipative force may be a force coming from laminar viscous stresses with inviscid velocity distributions. Both the flow solvers are combined to build a three-dimensional flow field through a few iterations. For an effect of the distortion of H-S flow surface the body forces are computed after each B-B flow calculation is finished. Mizuki's centrifugal impellers are tested numerically. The reliability of the numerical solution compared with experimental data is guaranteed.

• Water for future
• /
• v.27 no.3
• /
• pp.95-105
• /
• 1994

The boundary layer flow under a sluice gate is numerically solved by the random vortex sheet method combined with the vortex-in-cell method in a boundary-fitted coordinate system. The numerical solution shows that the boundary layer developed along the vertical sluice gate wall is the primary cause for the discrepancy in the contraction ratio between the laboratory experiments and inviscid theory; the bottom boundary layer plays much a smaller role in the discrepancy. By dimensional analysis it is concluded that the discrepancy is inversely proportional to the 3/4th power of the gate opening, as analyzed by Benjamin(1956). The results of the numerical simulation and dimensional analysis show a good agreement with experimental results obtained by Benjamin(1956).

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.20 no.6
• /
• pp.1994-2004
• /
• 1996

Numerical calculation was applied to supersonic under-expanded jets, and compared with the results of a linear theory and other experiments. TVD difference scheme was employed to solve 2-dimensional and axisymmetric inviscid Euler equation. This paper aims to explore the effects of angle of divergence and design Mach number of nozzle on the structure of under-expanded jets. The angle of divergence was varied from 0 to 20 deg. The results show that the length of the first cell of the under-expanded jets decreases and Mach disk generates at lower nozzle pressure ratio, if the angle of divergence or design Mach number of nozzle increases. The distance from the nozzle exit to Mach disk in 2-dimensional jets becomes much larger than that of axisymmetric jets, and the widths of the jet boundary and the barrel shock wave are also larger than that of axisymmetric jets. Calculation results indicate that the configuration of the under-expanded jets is strongly dependent on the nozzle pressure ratio.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.100-105
• /
• 2000

The present work describes the prediction method for the unsteady flow field and the acoustic pressure field of a ducted axial fan. The prediction method is comprised of time-marching free-wake method, acoustic analogy. and the Helmholtz-Kirchhoff BEM. The predicted sound signal of a rotor is similar to the experiment one. We assume that the rotor rotates with a constant angular velocity and the flow field around the rotor is incompressible and inviscid. Then, a time-marching free-wake method is used to model the fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lawson's method is used to predict the acoustic source. The newly developed Helmholtz-Kirchhoff BEM for thin body is used to calculate the sound field of the ducted fan. The ducted fan with 6 blades is analysed and the sound field around the duct is calculated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.113-118
• /
• 2000

Centrifugal pump are widely used and the noise generated by these machines causes one of the most serious problems. In general, the centrifugal pump noise is often dominated by tones at BPF(blade passage frequency) and its higher harmonics. This is a consequence of the strong interaction between the flow discharged from the impeller and the cutoff in the casing. However, only a few researches have been carried out on predicting the noise because of the difficulty in obtaining detailed information about the flow field and casing effects on noise radiation. The objective of this study is to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal pump, and to calculate the effects of small vanes that are attached in original impeller-splitter impeller. We assume that the impeller rotates with a constant angular velocity and the flow field around the impeller is incompressible and inviscid. So, a discrete vortex method(DVM) is used to model the centrifugal pump and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The splitter impeller changes the acoustic characteristics as well as performance. Two-splitter type impeller is good for acoustic characteristics.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.1
• /
• pp.123-130
• /
• 2008

An analytical method for the free vibration of a flexible rectangular plate in contact with water is developed by the Rayleigh-Ritz method. The plate clamped along the edges is partially contacted with water at both sides. It is assumed that the contained water is incompressible and inviscid. The wet mode shape of the plate is assumed as a combination of the dry mode shapes of a clamped beam. The liquid motion is described by using the liquid displacement potential and determined by using the compatibility conditions along the liquid interface with the plate. Minimizing the Rayleigh quotient based on the energy conservation gives an eigenvalue problem. It is found that the theoretical results can predict excellently the fluid-coupled natural frequencies comparing with the finite element analysis result.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.10a
• /
• pp.47-52
• /
• 2004

In the design of scramjet intake for hypersonic flight, a variety of aerothermodynamics phenomena are encountered. These phenomena include blunt leading - edge effects, boundary layer development issues, transition, inviscid / viscous coupling, shock - shock interactions, shock / boundary - layer interactions, and flow profile effects. For intakes that are designed to operate within a narrow Mach number / altitude envelope, an understanding of a few of these phenomena might be required. In this work several predominant flowfield phenomena (viscous phenomena, boundary - layer separation, and combustor entrance profile) are discussed to investigate the performance of the intake at the altitude and angle of attack extremes of the HyShot flight experiment.

• Journal of the Society of Naval Architects of Korea
• /
• v.37 no.4
• /
• pp.11-18
• /
• 2000

This paper is concerned with the generation of an optimal forward hull form by a nonlinear programming method. A Rankine source panel method based on the inviscid and potential flow approximation is employed to calculate the wave-making resistance and SQP method is also used for the optimization. The hull form is represented by a spline function. The forward hull form of a minimum wave resistance with the given design constraints is generated. In addition, the forward hull form of a minimum total resistance by considering the frictional resistance together with an empirical form factor is produced and compared with the former result.

• Journal of the Society of Naval Architects of Korea
• /
• v.53 no.2
• /
• pp.92-100
• /
• 2016

A comparative method is applied to evaluate well-known formulas for estimating the size of supercavities of axisymmetric cavitators for the supercavitating underwater vehicle. Basic functional forms of these formulas are derived first for the cavity diameter from a momentum integral estimate and second for the cavity length from an asymptotic analysis of inviscid supercavity flows. The length and the diameter of axisymmetric supercavities estimated by each formula are compared, with available experimental data for a disk and a 45° conical cavitators, and also with computational results obtained by a CFD code, ‘fluent’, for conical cavitators of wide range of cone angles. Results for estimating the length and the diameter of the supercavities show in general a good agreement, which confirms the size of the supercavities for disk and conical cavitators can be estimated accurately by these simple formulas of an elementary function of cavitation number and drag coefficient of the cavitator. These formulas will be useful for from conceptual design of the cavitator to real-time control of the supercavitating underwater vehicle.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.23 no.2
• /
• pp.93-114
• /
• 2019

In this article, we introduce a finite difference method for solving the Navier-Stokes equations in rectangular domains. The method is proved to be energy stable and shown to be second-order accurate in several benchmark problems. Due to the guaranteed stability and the second order accuracy, the method can be a reliable tool in real-time simulations and physics-based animations with very dynamic fluid motion. We first discuss a simple convection equation, on which many standard explicit methods fail to be energy stable. Our method is an implicit Runge-Kutta method that preserves the energy for inviscid fluid and does not increase the energy for viscous fluid. Integration-by-parts in space is essential to achieve the energy stability, and we could achieve the integration-by-parts in discrete level by using the Marker-And-Cell configuration and central finite differences. The method, which is implicit and second-order accurate, extends our previous method [1] that was explicit and first-order accurate. It satisfies the energy stability and assumes rectangular domains. We acknowledge that the assumption on domains is restrictive, but the method is one of the few methods that are fully stable and second-order accurate.