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

The interaction between hypersonic free stream and side jet flow at high altitudes is investigated by using the direct simulation Monte Carlo (DSMC) method. In order to alleviate the difficulty associated with the large density difference between the free stream and the side jet flow and to simulate the two flows simultaneously, a weighting factor technique is applied. For validation, the corner flow over a pair of plates perpendicularly attached is calculated with and without a side jet, and the results are compared with experiment. For a more realistic configuration, the flow past a blunted cone cylinder shape is solved. The leeward or windward jet is injected into the free stream and the effect on the aerodynamic force and moment is observed at various flow angles. The lambda shock effect and the wake structure are studied in terms of the surface pressure differential. A higher interaction between the free stream and the side jet flow is observed when the side jet is injected in the windward direction.

• Journal of computational fluids engineering
• /
• v.10 no.1
• /
• pp.56-62
• /
• 2005

Large eddy simulation of vortexing flow of molten steel in the continuous casting mold with and without DC magnetic field was conducted. The influence of the position of magnetic field to the residence time and depth of the vortex was analyzed. The mechanism of the influence of magnetic field to the vortexing flow was found. The computational results show that the vortexing flow is the result of shearing of the two un-symmetric surface flows from the mold narrow faces when they meet adjacent to the SEN; the un-symmetric flow for turbulent vortex is caused by turbulent energy of the fluid and that for biased vortex is caused by biased flow and the turbulent energy of fluid; with the moving of the magnetic field from the centerline of the outlet of the SEN to the free surface, the surface velocity is decreased gradually and the depth of the turbulent vortex and the biased vortex is decreased, the residence time is increased with the magnetic field moves from DL=120mm to DL=60mm and then decreased; the turbulent vortex and the biased vortex can be eliminated when the magnetic field is located at the free surface.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.1
• /
• pp.43-57
• /
• 1999

The fluid flow and heat transfer in a thin liquid film are investigated numerically. The flow Is assumed to be two-dimensional laminar and surface tension is considered. The most important characteristics of this flow is the existence of a hydraulic jump through which the flow undergoes very sharp and discontinuous change. Arbitrary Lagrangian-Eulerian(ALE) method is used to describe moving free boundary and a modified SIMPLE algorithm based on streamline upwind Petrov-Galerkin(SUPG) finite element method is used for time marching iterative solution. The numerical results obtained by solving unsteady full Navier-Stokes equations are presented for planar and radial flows subject to constant wall temperature or constant wall heat flux, and compared with available experimental data. It Is discussed systematically how the inlet Reynolds and Froude numbers and surface tension affect the formation of a hydraulic jump. In particular, the effect of temperature dependent fluid properties is also discussed.

• Kyungpook Mathematical Journal
• /
• v.57 no.4
• /
• pp.651-665
• /
• 2017

Some computational fluid dynamics problems concerning the thin films flow of viscous fluid with a free surface and draining or coating fluid-flow problems can be delineated by third-order ordinary differential equations. In this paper, the aim is to introduce the numerical solutions of the boundary value problems of such equations by variational iteration method. In this paper, it is shown that the third-order boundary value problems can be written as a system of integral equations, which can be solved by using the variational iteration method. These solutions are gleaned in terms of convergent series. Numerical examples are given to depict the method and their convergence.

• Proceedings of the KSME Conference
• /
• 2004.11a
• /
• pp.1466-1471
• /
• 2004

In this paper, we report the numerical and experimental solutions of the axi-symmetric flows in the axial plane driven by an impingement of fluid from the bottom wall of a circular cylinder. We managed to visualize successfully the flow pattern shown on the vertical plane through the container axis. The numerical results are not show to compare well with the experimental results for the case of the Rossby number 3. Because the numerical results calculate on the assumption that vortex flows are axi-symmetric flow on the other hand real experimental results are show asymmetric flow. The numerical solutions reveal that inertial oscillation plays an important role at small Rossby numbers, or at a larger background rotation.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.435-447
• /
• 2011

A Computational Fluid Dynamics model is presented in this study for the simulation of the complex fluid flows with free surfaces inside the Tuned Liquid Column Dampers in horizontal motion. The characteristics of the fluid model of the TLCD in horizontal motion include the free surface of the multiphase flow and the horizontal moving frame. In this study, the time depend unsteady Standard ${\kappa}-{\varepsilon}$ turbulent model based on Navier-Stokes equations is chosen. The volume of fluid (VOF) method and sliding mesh technique are adopted to track the free surface of water inside the vertical columns of TLCD and treat the moving boundary of the walls of TLCD in horizontal motion. Several model solution parameters comprising different time steps, mesh sizes, convergence criteria and discretization schemes are examined to establish model parametric independency results. The simulation results are compared with the experimental data in the dimensionless amplitude of the water column in four different configured groups of TLCDs with four different orifice areas. The predicted natural frequencies and the head loss coefficient of TLCDs from CFD model are also compared with the experimental data. The predicted numerical results agree well with the available experimental data.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.35 no.7
• /
• pp.723-733
• /
• 2011

A numerical method of the CSF(Continuum Surface Force) model is presented for the calculation of the surface tension force and implemented in an in-house solution code(PowerCFD). The present method(code) employs an unstructured cell-centered method based on a conservative pressure-based finite-volume method with volume capturing method(CICSAM) in a volume of fluid(VOF) scheme for phase interface capturing. The application of the present method to a 2-D liquid drop problem is illustrated by an equilibrium and nonequilibrium oscillating drop calculation. It is found that the present method simulates efficiently and accurately surface tension-dominant multiphase flows.

• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.1
• /
• pp.38-52
• /
• 1991

In this paper, a semi-Lagrangian method is used to solve the nonlinear hydrodynamics of a three-dimensional body beneath the free surface in the time domain. The boundary value problem is solved by using the boundary integral method. The geometries of the body and the free surface are represented by the curved panels. The surfaces are discretized into the small surface elements using a bi-cubic B-spline algorithm. The boundary values of $\phi$ and $\frac{\partial{\phi}}{\partial{n}}$ are assumed to be bilinear on the subdivided surface. The singular part proportional to $\frac{1}{R}$ are subtracted off and are integrated analytically in the calculation of the induced potential by singularities. The far field flow away from the body is represented by a dipole at the origin of the coordinate system. The Runge-Kutta 4-th order algorithm is employed in the time stepping scheme. The three-dimensional form of the integral equation and the boundary conditions for the time derivative of the potential Is derived. By using these formulas, the free surface shape and the equations of motion are calculated simultaneously. The free surface shape and fille forces acting on a body oscillating sinusoidally with large amplitude are calculated and compared with published results. Nonlinear effects on a body near the free surface are investigated.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.126-134
• /
• 2008

We developed a three dimensional Navier-Stokes code based on the level set method to simulate two phase flows with high density ratio. The Navier-Stokes equations with consideration of the surface tension effects are solved by using SIMPLE algorithm on a non-staggered grid. The present code is validated by simulating two test problems. First one is to simulate a rising bubble inside a cube. The thickness of the interface of the bubble is shown to affect the pressure distribution around the interface. As the thickness decreases, the pressure field around the interface becomes more oscillatory. As the bubble rises, a ring vortex is shown to form around the interface and the bubble eventually develops into an ellipsoidal shape. Merge of two bubbles inside a container is secondly tested to show the robustness of the present code for two phase flow simulation. Numerical results show stable and reliable behavior during the process of merging of two bubbles. The velocity and pressure fields around the interface of bubbles are shown oscillation free during the merging of two bubbles.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.6
• /
• pp.823-832
• /
• 2002

A modified $textsc{k}$-$\varepsilon$model is proposed for calculation of transitional boundary-layer flows with changing pressure gradient. In order to develop the model for this problem, the flow is divided into three regions; pre-transition region, transition region and fully turbulent region. The effect of pressure gradient is taken into account in stream-wise intermittency factor, which bridges the eddy-viscosity models in the pre-transition region and the fully turbulent region. From intermittency data in various flows, Narashima＇s intermittency function, F(${\gamma}$), has been found to be proportional to $\chi$$^{n}$ according to the extent of pressure gradient. Three empirical correlations of intermittency factor being analyzed, the best one was chosen to calculate three benchmark cases of bypass transition flows with different free-stream turbulence intensity under arbitrary pressure gradient. It was found that the variations of skin friction and shape factor as well as the profiles of mean velocity in the transition region were very satisfactorily predicted.