• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.145-155
• /
• 1998

A computer code has been developed for the computation of the viscous flow around a ship model with the free surface. In this code, the incompressible Reynolds-averaged Navier-Stokes equations are solved numerically by a finite difference method which employes second-order finite differences for the spatial discretization and a four-stage Runge-Kutta scheme for the temporal integration of the governing equations. For the turbulence closure, a modified version of the Baldwin-Lomax model is exploited. The location of the free surface is determined by solving the equation of the kinematic free-surface condition using the Lax-Wendroff scheme and the boundary-fitted grid is generated at each time step so that one of the grid surfaces always coincides with the free surface. An inviscid approximation of the dynamic free-surface boundary condition is applied as the boundary conditions for the velocity and pressure on the free surface. To validate the computational method and the computer code developed in the present study, the numerical computations are carried out for both Wigley parabolic hull and Series 60 $C_B=0.6$ ship model and the computational results are compared with the experimental data.

• Journal of computational fluids engineering
• /
• v.14 no.1
• /
• pp.70-77
• /
• 2009

In this study, time-dependent numerical analysis was carried out to investigate the plasma jet impingement on a flat plate, and a compressible form of two-dimensional inviscid gas dynamics equations were solved using the flux corrected transport algorithm. The mathematical modeling of Joule heating in the polycarbonate capillary bore and the mass ablation from the bore wall was incorporated in the numerical analysis and the series of computation was performed for three cases depending on the distance of the opposing plate from the capillary exit. The computational results reveal that the presence of the opposing plate does not affect the flow conditions inside the capillary when compared to the case of open-air plasma discharge. In the exterior region, the flow structure shows the typical supersonic underexpanded jet which consists of the strong Mach disk in front of the opposing plate and the barrel shock at the side of the jet. It is found that the shock evolution becomes more quasi-steady when the plate distance decreases. Also, the effects of the distance between the capillary bore exit and the opposing plate on the flow conditions along the opposing plate are investigated and the pressure variation on the plate shows more complicated interaction between the plasma discharge and the opposing plate when the location of plate becomes closer to the capillary exit.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.30 no.6 s.249
• /
• pp.523-530
• /
• 2006

One of the fascinating prospects is the possibility of new hydrodynamics technology on micro-scale system since oscillations of micro-droplets are of practical and scientific importance. It has been widely conceived that the lowest oscillation mode of a pendant droplet is the longitudinal vibration, i.e. periodic elongation and contraction along the longitudinal direction. Nonlinear and forced oscillations of supported viscous droplet were focused in the present study. The droplet has a free contact line with solid plate and inviscid fluid. Natural frequencies of a pendant droplet have been investigated experimentally by imposing the acoustic wave while the frequency is being increased at a fixed amplitude. It is found that a pendant droplet shows the resonant behaviors at each mode similar to the theoretical analysis. The rotation of the droplet about the longitudinal axis is the oscillation mode of the lowest resonance frequency. This rotational mode can be invoked by periodic acoustic forcing and is analogous to the pendulum rotation. It is also found that the natural frequency of a pendant droplet is independent of the drop density and surface tension but inversely proportional to the square root of the droplet size.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1384-1391
• /
• 2001

Three-dimensional flow analysis is implemented to investigate the flow through transonic axial-flow compressor rotor(NASA R67) and to evaluate the performances of Abid's low-Reynolds-number k-$\omega$ and Baldwin-Lomax turbulence models. A finite volume method is used fur spatial discretization. The equations are solved implicitly in time by the use of approximate factorization. The upwind difference scheme is used for inviscid terms and viscous terms are approximated with central difference. The flux-difference-splitting method of Roe is used to obtain fluxes at the cell faces. Numerical analysis is performed near peak efficiency and near stall. The results are compared with the experimental data for NASA R67 rotor. Blade-to-Blade Mach number distributions are compared to confirm the accuracy of the code. From the results, it is concluded that Abid'k-$\omega$ model is better for the calculation of flow rate and efficiency than Baldwin-Lomax model. But, the predictions for Mach number and shock structure are almost the same.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.10
• /
• pp.2617-2629
• /
• 1995

The flow through a centrifugal compressor rotor was calculated using the quasi-3-dimensional and fully 3-dimensional Navier-Stokes solution methods. The calculated results, obtained during the development of the computer codes for both methods are discussed. In the inviscid quasi 3-dimensional analysis, stream function formulation was used for the blade to blade (B-B) plane calculations, and the streamline curvature method was used for the meridional (H-S) plane calculations. In the viscous 3-dimensional flow analysis, a control volume method based on a general rotating curvilinear coordinate system was used to solve the time-averaged Navier-Stokes equations, and a standard k-.epsilon. model was used to obtain eddy viscosity. The quasi-3-dimensional analysis reasonably predicts the pressure distributions and requires much less computation time in the region where viscous effects are not strong; however, it fails to predict velocity field and loss mechanism through the impeller passage. The viscous 3-dimensional flow analysis shows reasonable pressure distributions and typical jet-wake flow field through the impeller passage. Secondary flow and total pressure distributions on cross-sectional planes explain the loss mechanisms through the impeller.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.12 no.1
• /
• pp.23-28
• /
• 2008

Flow analysis in a turbine cascade by Euler or Navier-Stokes equation gives relatively accurate solution, however, those method require large computer memory or computing time. On contrast, the panel method, which is applied to incompressible and inviscid flow, provides fast and reasonal solution but the compressibility correction is required for a high air velocity case. In this paper, the compressibility corrected panel method was applied in order to find velocity distribution on turbine blades. Results showed that the calculated velocity in a turbine cascade by the compressibility corrected panel method gave good agreement with the solution by finite volume method for compressible flow.

• Journal of KSNVE
• /
• v.7 no.5
• /
• pp.861-869
• /
• 1997

An approximate analytical method has been developed for estimating hydrodynamic forces acting on oscillating inner cylinder in concentric annulus. When the rigid inner cylinder executes translational oscillation, fluid inertia and damping forces on the oscillating cylinder are generated by unsteady pressure and viscous skin friction. Considering the dynamic-characteristics of unsteady viscous flow and the added mass coefficient of inviscid fluid, these hydrodynamic forces including viscous effect are dramatically simplified and expressed in terms of oscillatory Reynolds number and the geometry of annular configuration. Thus, the viscous effect on the forces can be estimated very easily compared to an existing theory. The forces are calculated by two models developed for relatively high and low oscillatory Reynolds numbers. The model for low oscillatory Reynolds number is suitable for relatively high ratio of the penetration depth to annular space while the model for high oscillatory Reynolds number is applicable to the case of relatively low ratio. It is found that the transient ratio between two models is approximately 0.2~0.25 and the forcea are expressed in terms of oscillatory Reynolds number, explicity. The present results show good agreements with an existing numerical results, especially for high and low penetration ratios to annular gap.

• Journal of KSNVE
• /
• v.9 no.5
• /
• pp.955-965
• /
• 1999

Centrifugal fans are widely used and the noise generated by the these machines causes one of the most serious problems. In general, the centrifugal fan 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 form 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 understand the generation mechanism of sound and to develop a prediction method for the unsteady flow field and the acoustic pressure field of a centrifugal fan. We assume that the impeller rotates with a constant angular velocity and the flow field of the impeller is incompressible and inviscid. So, a discrete vortex method (DVM) is used to model the centrifugal by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. A centrifugal impeller and wedge introduced by Weidemann are used in the numerical calculation and the results are compared with the experimental data. Reasonable results are obtained not only for the peak frequencies but also for the amplitudes of the tonal sound.

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

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.3 no.3
• /
• pp.176-183
• /
• 1991

In this paper. boundary element method is applied to the analysis of nonlinear free surface wave. A particular concern is given to the treatment of the open boundaries at the in-flow boundary and out-flow boundary, which uses the mass-flux and energy-flux considering the continuity of fluid. By assuming the fluid to be inviscid and incompressible and the flow to be irrotational. the problem is formulated mathematically as a two-dimentional nonlinear problem in terms of a velocity potential. The equation(Laplace equation) and the boundary conditions are transformed into two boundary integral equations. Due to the nonlinearity of the problem. the incremental method is used for the numerical analysis. Numerical results obtained by the present boundary element method are compared with those obtained by the finite element method and also with experimental values.