• Journal of computational fluids engineering
• /
• v.2 no.1
• /
• pp.1-7
• /
• 1997

A FA-FD hybrid method, developed for solving three-dimensional incompressible Navier-Stokes equations, is applied to calculate three-dimensional laminar flows through a square duct with a 90° bend. The method discretizes the convective terms in the primary flow direction with 3rd-order upwind finite-differences and the convective and diffusive terms in the transverse directions with the two-dimensional finite analytic method. The non-staggered grid system is used and the pressure-velocity coupling is achieved by a global iteration procedure based on the PISO algorithm. Detailed comparisons between the computed solutions and the available experimental data are given mainly for the velocity distributions at cross-sections in a 90° bend of a square duct with both fully developed and developing entry flows. Although the computational result shows generally a good agreement with the experimental data, there are some significant discrepancies underlining the necessity of more accurate numerical methods as well as reliable experimental data for their validation.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.19 no.3
• /
• pp.220-227
• /
• 2007

Heat transfer and momentum transfer under conditions of both oscillating flow and oscillating pressure within pulse tubes show very different behavior from those for steady state conditions. The analytic solutions of axial velocity and temperature of the gas within pulse tubes were obtained by assuming that the variations in pressure and temperature were purely sinusoidal and small. The shear stress and the heat flux at the tube wall obtained from the solutions are expressed in terms of the cross-sectional averaged velocity, the difference between mean temperature and instantaneous cross-sectional averaged temperature and the difference between mean pressure and instantaneous pressure. It is shown that the complex shear factor, which has been applied to momentum transfer of incompressible oscillating flow, and the complex Nusselt number, which has been applied to either heat transfer with oscillating pressure only or heat transfer of incompressible oscillating flow, could also be used for momentum transfer and heat transfer subjected to both oscillating flow and oscillating pressure, respectively.

• 한국기계연구소 소보
• /
• s.19
• /
• pp.15-21
• /
• 1989

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.239-244
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved by successively in the moving least square sense. Some weighing functions were tested in order to investigate the up-winding effect for the convection term. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as FVM.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.66-72
• /
• 1997

Circular-to-rectangular transition ducts are used as exhaust components of high performance fighter aircraft with vectored thrust nozzles. Three-dimensional incompressible Navier-Stokes solver is used to analyze the transition duct. Cross sections of transition duct are defined by superelliptic equation. The grid system is generated by Non-Uniform Rational B-Spline, after generating surface grid by blending the cross sections. Good agreement between the results of the computational simulation and the experimental data is observed.

• Proceedings of the KSME Conference
• /
• 2007.05b
• /
• pp.3156-3161
• /
• 2007

In the field of yarn dyeing, the most generally employed method is a type of package dyeing which uses a package of cheeses stacked on a spindle made of a perforated tube. In order to understand the process of level dyeing, it is essential to perform a study of the porous flow through the spindle for the cheese dyeing method. In this paper, the axisymmetric, incompressible, Navier-Stokes equations are solved for several spindle configurations using a fully implicit finite volume scheme. For investigating the flow patterns through the spindle, porous diameter and porosity is varied in the present study. The computational results show that the total pressure loss depends only on the velocity of inflow regardless of porous diameter and porosity and a large percentage of the mass flow rate through the spindle is discharged at the upside of the spindle. Therefore, it is required to design a new spindle to obtain the level dyeing.

• Journal of computational fluids engineering
• /
• v.19 no.4
• /
• pp.20-28
• /
• 2014

An efficient solution algorithm for simulating free surface problem is presented. Navier-Stokes equations for variable density incompressible flow are employed as the governing equation on Cartesian meshes. In order to describe the free surface motion efficiently, VOF(Volume Of Fluid) method utilizing THINC(Tangent of Hyperbola for Interface Capturing) scheme is employed. The most time-consuming part of the current free surface flow simulations is the solution step of the linear system, derived by the pressure Poisson equation. To solve a pressure Poisson equation efficiently, the PCG(Preconditioned Conjugate Gradient) method is utilized. This study showed that the proper application of the preconditioner is the key for the efficient solution of the free surface flow when its pressure Poisson equation is solved by the CG method. To demonstrate the efficiency of the current approach, we compared the convergence histories of different algorithms for solving the pressure Poisson equation.

• 한국전산유체공학회:학술대회논문집
• /
• 1995.10a
• /
• pp.55-63
• /
• 1995

The flow field around a high-speed train including cross-wind effects has been simulated. This study solves 3-D unsteady incompressible Navier-Stokes equations in the inertial frame using the iterative time marching scheme. The governing equations are differenced with 1st-order accurate backward difference scheme for the time derivatives, 3th-order accurate QUICK scheme for the convective terms and 2nd-order accurate central difference scheme for the viscous terms. The Marker-and-Cell concept was applied to efficiently solve continuity equation, which is differenced with 2nd-order accurate central difference scheme. The 4th-order artificial damping is added to the continuity equation for numerical stability. A C-H type of elliptic grid system is generated around a high-speed train including ground. The Baldwin-Lomax turbulent model was implemented to simulate the turbulent flows. To validate the present procedure, the flow around a high speed train at constant yaw angle of $45^{\circ}\;and\;90^{\circ}$ has been simulated. The simulation shows 3-D vortex generation in the lee corner. The flow separation is also observed around the rear of the train. It has concluded that the results of present study properly agree with physical flow phenomena.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.11 no.1
• /
• pp.495-509
• /
• 2019

The objective of this study is to analyze the compressibility effects of multiphase cavitating flow during the water-entry process. For this purpose, the water-entry of a projectile at transonic speed is investigated computationally. A temperature-adjusted Tait equation is used to describe the compressibility effects in water, and air and vapor are treated as ideal gases. First, the computational methodology is validated by comparing the simulation results with the experimental measurements of drag coefficient and the theoretical results of cavity shape. Second, based on the computational methodology, the hydrodynamic characteristics of flow are investigated. After analyzing the cavitating flow in compressible and incompressible fluids, the characteristics under compressible conditions are focused upon. The results show that the compressibility effects play a significant role in the development of cavitation and the pressure inside the cavity. More specifically, the drag coefficient and cavity size tend to be larger in the compressible case than those in the incompressible case. Furthermore, the influence of entry velocities on the hydrodynamic characteristics is investigated to provide an insight into the compressibility effects on cavitating flow. The results show that the drag coefficient and the impact pressure vary with the entry velocity, and the prediction formulas for drag coefficient and impact pressure are established respectively in the present study.

• The Journal of the Acoustical Society of Korea
• /
• v.40 no.4
• /
• pp.261-269
• /
• 2021

Without any validation of the incompressible assumption, most of previous studies on cavitation flow and its noise have utilized numerical methods based on the incompressible Reynolds Average Navier-Stokes (RANS) equations because of advantage of its efficiency. In this study, to investigate the effects of the flow compressibility on the Tip Vortex Cavitation (TVC) flow and noise, both the incompressible and compressible simulations are performed to simulate the TVC flow, and the Ffowcs Williams and Hawkings (FW-H) integral equation is utilized to predict the TVC noise. The DARPA Suboff submarine body with an underwater propeller of a skew angle of 17 degree is targeted to account for the effects of upstream disturbance. The computation domain is set to be same as the test-section of the large cavitation tunnel in Korea Research Institute of Ships and Ocean Engineering to compare the prediction results with the measured ones. To predict the TVC accurately, the Delayed Detached Eddy Simulation (DDES) technique is used in combination with the adaptive grid techniques. The acoustic spectrum obtained using the compressible flow solver shows closer agreement with the measured one.