• Journal of Mechanical Science and Technology
• /
• v.20 no.10
• /
• pp.1741-1752
• /
• 2006

A combined Streamline Upwind Petrov-Galerkin method (SUPG) and segregated finite element algorithm for solving conjugate heat transfer problems where heat conduction in a solid is coupled with heat convection in viscous fluid flow is presented. The Streamline Upwind Petrov-Galerkin method is used for the analysis of viscous thermal flow in the fluid region, while the analysis of heat conduction in solid region is performed by the Galerkin method. The method uses the three-node triangular element with equal-order interpolation functions for all the variables of the velocity components, the pressure and the temperature. The main advantage of the presented method is to consistently couple heat transfer along the fluid-solid interface. Four test cases, which are the conjugate Couette flow problem in parallel plate channel, the counter-flow in heat exchanger, the conjugate natural convection in a square cavity with a conducting wall, and the conjugate natural convection and conduction from heated cylinder in square cavity, are selected to evaluate efficiency of the presented method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.34 no.8
• /
• pp.739-747
• /
• 2010

Numerical calculations are carried out for evaluating the natural convection induced by the temperature difference between a hot inner circular cylinder and a cold outer square enclosure. A two-dimensional solution for unsteady natural convection is obtained by using the finite volume method to model an inner circular cylinder that was designed by using the immersed boundary method (IBM) for a Rayleigh number of $10^7$. In this study, we investigate the effect of the location ($\delta$) of the inner cylinder, which is located along the vertical central axis of the outer enclosure, on the heat transfer and fluid flow. The natural convection changes from unsteady to steady state depending on the $\delta$. The two critical lower bound and upper bound positions are ${\delta}_{C,L}$ = 0.05 and ${\delta}_{C,U}$ = 0.18, respectively. Within these defined bounds, the thermal and flow fields are in steady state.

• Journal of Ocean Engineering and Technology
• /
• v.20 no.1 s.68
• /
• pp.48-54
• /
• 2006

In this paper, flowfield and acoustic-field around moving bodies are simulated by the Arbitrary Lagrangian Eulerian (ALE) formulation in the finite difference lattice Boltzmann method. Some effects are checked by comparing flaw about a square cylinder in ALE formulation and that in the fixed coordinates, and both agree very well. Matching procedure between the moving grid and fixed grid is also considered. The applied method in which the both grids are connected through buffer region is shown to be superior to moving overlapped grid. Dipole-like emissions of sound wave from harmonically vibrating bodies in two- and three-dimensional cases are simulated.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.04a
• /
• pp.71-78
• /
• 1997

The high precision analysis by the p-version of the finite element method are fairly well established as highly efficient method for linear elastic problems, especially in the presence of stress singularity. It has been noted that the merits of p-version are accuracy, modeling simplicity, robustness, and savings in user's and CPU time. However, little has been done to exploit their benefits in elasto-plastic analysis. In this paper, the p-version finite element model is proposed for the materially nonlinear analysis that is based on the incremental theory of plasticity, the associated flow rule, and von-Mises yield criteria. To obtain the solution of nonlinear equation, the Newton-Raphson method and initial stiffness method, etc are used. Several numerical examples are tested with the help of the square plates with cutout, the thick-walled cylinder under internal pressure, and the center cracked plate under tensile loading. Those results are compared with the there cal solutions and the numerical solutions of ADINA software.

• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• pp.71-76
• /
• 2004

An extension of our recently developed locally linear reconstruction scheme to 2 dimensional incompressible flow solver is presented. The solver is based on a semi-implicit fractional step method in which the convective term is discretized by Adams-Bashforth method and the diffusion term by Crank-Nicolson method. Several numerical examples are tested to demonstrate the mesh type independent accuracy of the solver, which include decaying vortex flow, square cavity flow, and flow around a circular cylinder. The above examples are solved on quadrilateral or hybrid meshes. For all numerical examples, we obtained reasonable results.

• Journal of Advanced Marine Engineering and Technology
• /
• v.30 no.5
• /
• pp.580-588
• /
• 2006

Two-dimensional turbulent flows past a square cylinder and cavity noise are simulated by the finite difference lattice Boltzmann method with subgrid turbulence model. The method, based on the standard Smagorinsky subgrid model and a single-time relaxation lattice Boltzmann method, incorporates the advantages of FDLBM for handling arbitrary boundaries. The results are compared with those by the experiments carried out by Noda & Nakayama and Lyn et al. Numerical results agree with the experimental ones. Besides, 2D computation of the cavity noise generated by flow over a cavity at a Mach number of 0.1 and a Reynolds number based on cavity depth of 5000 is calculated. The computation result is well presented a understanding of the physical phenomenon of tonal noise occurred primarily by well-jet shear layer and vortex shedding and an aeroacoustic feedback loop.

• Journal of the Korean Society for Heat Treatment
• /
• v.5 no.4
• /
• pp.224-232
• /
• 1992

Coolingability of coolants is important factor in cooling processor heat treatment of steel. Using standard apparatus and method defined in the Korean Industrial Standard three different shapes of probe were designed, ie, cylinderical, spherical and square on shape with same volume of standard probe. Distilled water and sodium aquious solutions with different concentration of NaOH, NaCl and $Na_2CO_3$ were examined. Estimation of coolingability of each quenchants for the probes of cylinderical, spherical or square shape, the cooling rate is greater square, cylinder and sphere in order. Coolingability of sodium aquious solution of NaCl, $Na_2CO_3$ and NaOH is found generally greater then that of distilled water. Effectiveness of ingredients is in the order of $Na_2CO_3$, NaOH and NaCl. In both solutions coolingability increases in 20%, 5%, and 10%in order. Analytical results obtained from Finite Element Method were compared with experimental ones and found as practically satisfactional.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.24 no.3
• /
• pp.247-255
• /
• 2012

Numerical calculations are carried out for the natural convection in a square enclosure with two hot cylinders induced by temperature difference between a cold outer rectangular cylinder and two hot circular cylinders. A two-dimensional solution for unsteady natural convection is obtained, using the immersed boundary method (IBM) to model two inner circular cylinders based on finite volume method, for different Rayleigh numbers varying over the range of $10^3$ to $10^5$. The study goes further to investigate the effect of the location of two cylinders on the heat transfer and fluid flow. The location of inner circular cylinders is changed vertically along the center-line of square enclosure. The changes of heat transfer quantities have been presented.

• Wind and Structures
• /
• v.14 no.5
• /
• pp.465-480
• /
• 2011

In this paper, a stabilized large eddy simulation technique is developed to predict turbulent flow with high Reynolds number. Streamline Upwind Petrov-Galerkin (SUPG) stabilized method and three-step technique are both implemented for the finite element formulation of Smagorinsky sub-grid scale (SGS) model. Temporal discretization is performed using three-step technique with viscous term treated implicitly. And the pressure is computed from Poisson equation derived from the incompressible condition. Then two numerical examples of turbulent flow with high Reynolds number are discussed. One is lid driven flow at Re = $10^5$ in a triangular cavity, the other is turbulent flow past a square cylinder at Re = 22000. Results show that the present technique can effectively suppress the instabilities of turbulent flow caused by traditional FEM and well predict the unsteady flow even with coarse mesh.

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

As an alternative for solving the incompressible Navier-Stokes equations, we present a vorticity-based integro-differential formulation for vorticity, velocity and pressure variables. One of the most difficult problems encountered in the vorticity-based methods is the introduction of the proper value-value of vorticity or vorticity flux at the solid surface. A practical computational technique toward solving this problem is presented in connection with the coupling between the vorticity and the pressure boundary conditions. Numerical schemes based on an iterative procedure are employed to solve the governing equations with the boundary conditions for the three variables. A finite volume method is implemented to integrate the vorticity transport equation with the dynamic vorticity boundary condition . The velocity field is obtained by using the Biot-Savart integral derived from the mathematical vector identity. Green's scalar identity is used to solve the total pressure in an integral approach similar to the surface panel methods which have been well-established for potential flow analysis. The calculated results with the present mettled for two test problems are compared with data from the literature in order for its validation. The first test problem is one for the two-dimensional square cavity flow driven by shear on the top lid. Two cases are considered here: (i) one driven both by the specified non-uniform shear on the top lid and by the specified body forces acting through the cavity region, for which we find the exact solution, and (ii) one of the classical type (i.e., driven only by uniform shear). Secondly, the present mettled is applied to deal with the early development of the flow around an impulsively started circular cylinder.