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

Three-Dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for the various contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreements.

• Proceedings of the KSME Conference
• /
• 2003.04a
• /
• pp.1919-1922
• /
• 2003

Bypass anastomosis are frequently adopted for surgical treatments. After the bypass grafting, the bypass artery is often occluded due to restenosis and/or anastomotic neointimal fibrous hyperplasia phenomena. Optimal coronary bypass anastomosis should be investigated to improve the patency for the arterial bypass techniques. The objective of this study is to investigate the influence of bypass with sequential bypass effects in the stenosed coronary artery. Numerical analyses are focused on the understanding of the flow patterns for different sequential anastomosis techniques. Blood flow field is treated as two-dimensional incompressible laminar flow. The finite volume method is adopted for discretization of the governing equations. The Carreau model is employed as the constitutive equation for blood. To find the optimal sequential bypass anastomotic configurations, the mass flow rates at the outlet of different models are compared quantitatively.

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

The three-dimensional Euler equations are solved numerically for the analysis of contraction flows in wind or water tunnels. A second-order finite difference method is used for the spatial discretization on the nonstaggered grid system and the 4-stage Runge-Kutta scheme for the numerical integration in time. In order to speed up the convergence, the local time stepping and the implicit residual-averaging schemes are introduced. The pressure field is obtained by solving the pressure-Poisson equation with the Neumann boundary condition. For the evaluation of the present Euler solver, numerical computations are carried out for three contraction geometries, one of which was adopted in the Large Cavitation Channel for the U.S. Navy. The comparison of the computational results with the available experimental data shows good agreement.

• Steel and Composite Structures
• /
• v.24 no.5
• /
• pp.523-536
• /
• 2017

A layerwise (LW) formulation based on the Galerkin method is presented to investigate the three-dimensional stress state in long sandwich plate which is subjected to tension force and pure bending moment. Based on the Galerkin method and the LW discretization approach, the equilibrium equations of elasticity for the long plate are written in the weak form and discretized through the thickness of the plate. The discretized equations are written in terms of displacement components of the numerical layers. The governing equations of the plate are solved analytically for the free edge boundary conditions. The distribution of stress state especially the 3D stress state in the vicinity of the edges of the sandwich plate which is subjected to tension and pure bending is studied. In order to increase the accuracy, the out of plane stresses are obtained by integrating the equilibrium equations of elasticity. The convergence and accuracy of the predictions are studied and various numerical results are presented for distribution of the in-plane and out of plane stresses in symmetric and un-symmetric sandwich plates.

• East Asian mathematical journal
• /
• v.33 no.1
• /
• pp.53-65
• /
• 2017

Numerical solutions of the Rosenau-Burgers equation based on the cubic B-spline finite element method are introduced. The backward Euler method is used for discretization in time, and the obtained nonlinear algebraic system is changed to a linear system by the Newton's method. We show that those methods are unconditionally stable. Two test problems are studied to demonstrate the accuracy of the proposed method. The computational results indicate that numerical solutions are in good agreement with exact solutions.

• Journal of Power System Engineering
• /
• v.13 no.5
• /
• pp.11-17
• /
• 2009

A numerical study of unsteady mixed convection in a cavity with high viscous fluid is presented. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Rayleigh number ($10^3\;{\leq}\;Ra\;{\leq}\;10^5$), the Reynolds number (0 < Re $\leq$ 1), and the aspect ratio (0.5 $\leq$ AR $\leq$ 2). The fluid used is silicon oil, a high prandtl number fluid, Pr = 909.1. The results show velocity vectors and temperature distributions. It is found that the periodic flows in a cavity are observed at very low Reynolds numbers, and the period of periodic flow decreases with increasing Reynolds and Rayleigh numbers, and increases with increasing aspect ratio. Also, the Reynolds number range of periodic flow increases with increasing Rayleigh numbers and aspect ratio.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.139-144
• /
• 1999

There have been many studies for heat transfer enhancement. Particularly, the study of flow in heat exchangers which have fin device has been main theme in heat transfer area. Practically, the circular tube which has internal fins is widely used for developing heat transfer rate. In this study, flow and heat transfer analysis of the circular tube with fins are investigated. The height and the number of fins are arbitrary. The flow field is assumed to be laminar. The conformal mapping is used for analytic solution of the laminar flow field. Discretization of governing equation, namely, FDM was used for numerical analysis. The velocity field, flow rate and shear stress are calculated for some numbers of fins in circular tube and for some heights of fin. Temperature fields are plotted along the tube length. It can be shown that the numerical solution agrees with the analytical solution.

• Journal of the Korean Mathematical Society
• /
• v.48 no.3
• /
• pp.529-550
• /
• 2011

In this work, a numerical solution of the incompressible Navier-Stokes equations is proposed. The method suggested is based on an algorithm of discretization by mixed finite elements with a posteriori error estimation of the computed solutions. In order to evaluate the performance of the method, the numerical results are compared with some previously published works or with others coming from commercial code like Adina system.

• Journal of Power System Engineering
• /
• v.11 no.4
• /
• pp.56-64
• /
• 2007

A numerical simulation is carried out mixed convection in horizontal channel with a heat source from below of rectangular cavity. Finite volume method was employed for the discretization and PISO algorithm was used for calculating pressure term. The parameters governing the problem are the Reynolds number ($10^{-2}{\leq}Re{\leq}50$), the Rayleigh number ($10^3{\leq}Ra{\leq}2.06{\times}10^5$), the Prandtl number ($0.72{\leq}Pr{\leq}909$), the aspect ratio ($0.5{\leq}AR=W/H{\leq}2$) and the angle of inclination ($0^{\circ}{\theta}60^{\circ}$). Mean Nusselt number distributions were obtained and effect of Reynolds number, Rayleigh number and Prandtl number on mixed convection in the horizontal channel with rectangular cavity were investigated.

• Nuclear Engineering and Technology
• /
• v.32 no.1
• /
• pp.46-56
• /
• 2000

Numerical results are presented for the 2-dimensional turbulent transient heat transfer of the shell/tube heat exchanger with a step change of the inlet temperature in the primary side. Heat transfer boundary conditions outside the pipe are given partially by the convection heat transfer conditions and partially by insulated conditions. Calculation results were obtained by solving the unsteady two-dimensional elliptic forms for the Reynolds-averaged governing equations for the mass, momentum and energy. Finite-difference method was used to obtain discretization equations, and the SIMPLER solution algorithm was employed for the calculation procedure. Turbulent model used is the algebraic model proposed by Cebeci-Smith. Results presented include the time variant Nusselt number distribution, average temperature distribution and outlet temperatures for the various inlet temperatures and flow rates.