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

The three-dimensional laminar mixed convection flow between the conductive printed circuit boards. on which the heat generating rectangular blocks are uniformly distributed, has been examined in the present study. The flow and heat-transfer characteristics are assumed to be pseudo periodic in the streamwise direction and symmetric in the cross-stream direction. Using an algorithm of SIMPLER, the continuity equation. the Navier-Stokes equations and the energy equation are solved numerically in the three-dimensional domain Inside the channel. The convective derivative terms are discretized by the QUICK scheme to accurately capture the flow field. The flow and the heat transfer characteristics are thoroughly examined for various Re and Gr.

• Journal of computational fluids engineering
• /
• v.11 no.4 s.35
• /
• pp.39-47
• /
• 2006

In this paper, the authors developed a multidimensional interpolation method inside a finite volume cell in the computation of high-order accurate numerical flux such as the fifth order WEND (weighted essentially non-oscillatory) scheme. This numerical method starts from a simple Taylor series expansion in a proper spatial order of accuracy, and the WEND filter is used for the reconstruction of sharp nonlinear waves like shocks in the compressible flow. Two kinds of interpolations are developed: one is for the cell-averaged values of conservative variables divided in one mother cell (Type 1), and the other is for the vertex values in the individual cells (Type 2). The result of the present study can be directly used to the cell refinement as well as the convective flux between finer and coarser cells in the Cartesian adaptive grid system (Type 1) and to the post-processing as well as the viscous flux in the Navier-Stokes equations on any types of structured and unstructured grids (Type 2).

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.39-46
• /
• 1999

The FVM(Finite Volume Method) have been used mainly for the flow analyses in the piston-cylinder. The objective of the present study is to analyze numerically the piston-driven intake flows using the FEM(Finite Element Method). The FEM algorithm used in this study is 4-step time-splitting method which requires much less execution time and computer storage than the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the momentum equations to prevent checkerboard pressure oscillations. Also, the ALE(arbitrary Lagrangian Eulerian) method is adopted for the moving grids. The calculated results show good agreement in comparison with those by the FVM and the experimental results by the LDA.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.98-105
• /
• 1999

The Pressure-based Unstructured-grid Navier-Stokes Solver PUNS-2/3D for incompressible steady and unsteady viscous flows has been developed. It is based on nonstaggered cell-centered finite volume method. Second-order upwind scheme with least-square reconstruction is used for convective fluxes. The SIMPLE method is implemented to couple the pressure and velocity fields. And the time derivatives in the momentum equations are discretised using a second-order Euler backward-differencing scheme. The discretised linear equations are solved by the preconditioned Biconjugate Gradient Stabilized method(Bi-CGSTAB). The developed solver is applied to validation problems using hybrid meshes.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.11
• /
• pp.3031-3040
• /
• 1995

Two-dimensional incompressible Navier-Stokes equations have been solved by the node-centered finite volume method with the unstructured triangular meshes. The pressure-velocity coupling is handled by the artificial compressibility algorithm due to its computational efficiency associated with the hyperbolic nature of the resulting equations. The convective fluxes are obtained by the Roe's flux difference splitting scheme using edge-based connectivities and higher-order differences are achieved by a reconstruction procedure. The time integration is based on an explicit four-stage Runge-Kutta scheme. Numerical procedures with local time stepping and implicit residual smoothing have been implemented to accelerate the convergence for the steady-state solutions. Comparisons with experimental data and other numerical results have proven accuracy and efficiency of the present unstructured approach.

• 한국가시화정보학회:학술대회논문집
• /
• 2006.12a
• /
• pp.149-152
• /
• 2006

This paper describes the flow visualization around the model of a commercial passenger airplane, Boeing 747-400, which cruises in the transonic speed. The geometry was realized through the reverse engineering based on the photogrammetry. The three-dimensional inviscid steady compressible governing equations are solved in the unstructured grid system under the cruise condition and in a finite volume method. The convective term is processed by the Crank-Nicholson scheme and first order upwind scheme is applied. The lift and drag forces in the wing with engines increase by 1.49% End 3.9%, respectively compared with the wing without engines.

• Structural Engineering and Mechanics
• /
• v.47 no.5
• /
• pp.599-622
• /
• 2013

Seismic response of two dimensional liquid tanks is numerically simulated using fully nonlinear velocity potential theory. Galerkin-weighted-residual based finite element method is used for solving the governing Laplace equation with fully nonlinear free surface boundary conditions and also for velocity recovery. Based on mixed Eulerian-Lagrangian (MEL) method, fourth order explicit Runge-Kutta scheme is used for time integration of free surface boundary conditions. A cubic-spline fitted regridding technique is used at every time step to eliminate possible numerical instabilities on account of Lagrangian node induced mesh distortion. An artificial surface damping term is used which mimics the viscosity induced damping and brings in numerical stability. Four earthquake motions have been suitably selected to study the effect of frequency content on the dynamic response of tank-liquid system. The nonlinear seismic response vis-a-vis linear response of rectangular liquid tank has been studied. The impulsive and convective components of hydrodynamic forces, e.g., base shear, overturning base moment and pressure distribution on tank-wall are quantified. It is observed that the convective response of tank-liquid system is very much sensitive to the frequency content of the ground motion. Such sensitivity is more pronounced in shallow tanks.

• 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.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.13 no.2
• /
• pp.299-306
• /
• 1989

The combined convective heat transfer from vertical parallel plates with constant temperature has been studied by numerical method. The governing equations for the system are solved by the finite difference method and successive over relaxation scheme for Re$_{L}$ = 50 - 500, Gr = 10$^{4}$, Pr = 0.7. Results for various plate spacings and plate lengths are as follows ; For various plate spacings, the mean Nusselt number increases and then decreases as the dimensionless plate spacing increases. The optimum plate spacing for maximum mean Nusselt number decreases with increasing Reynolds number and can be expressed as a function of Reynolds number. For various plate lengths, the mean plate Nudest number increases as the dimensionless plate length decreases and Reynolds number increases.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.4
• /
• pp.71-78
• /
• 2005

The vaporization characteristics of a liquid heptane droplet in a supercritical nitrogen flow is numerically analyzed. The present model can account for real gas effects, liquid-phase internal circulation, variable thermodynamic properties and high-pressure effects. Time marching method with preconditioning scheme is employed to handle the low Mach number flows in dense heptane droplet region. Computations are made for the wide range of convective velocity and ambient pressure. Numerical results indicate that the droplet deformation becomes stronger by increasing the Reynolds number and it becomes relatively weak by increasing the pressure.