• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.31 no.7 s.262
• /
• pp.620-627
• /
• 2007

In this paper, we adapt a modified internal energy non-equilibrium first-order extrapolation thermal boundary condition to the thermal lattice Boltzmann model (TLBM). This model is the double populations approach to simulate hydrodynamic and thermal fields. The bounce-back boundary condition which is a traditional boundary condition of lattice Boltzmann method has only a first order in numerical accuracy at the boundary and numerical instability. A non-equilibrium first-order extrapolation boundary condition has been verified to be of better numerical stability than the bounce-back boundary condition and this boundary condition is proved to be of second-order accuracy for the flat boundaries. The two-dimensional natural convection flow in a square cavity with Pr=0.71 and various Rayleigh numbers are simulated. The results are found to be in good agreement with those of previous studies.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.37-44
• /
• 2000

Numerical investigation is made to study the effects of natural convection on the formation of macrosegregation of a Pb-Sn alloy solidification process in a 2-D confined rectangle mold. The governing equations are calculated using previous continuum models with SIMPLE algorithm doring the solidification process. In addition. to track the solid-liquid interface with time variations. the moving boundary condition Is adopted and irregular interface shapes are treated with Boundary-Fitted Coordinate system. As the temperature reduce from the liquidus to the solidus, the liquid concentration of Sn. the lighter constituent, increases. Then the buoyancy-driven flow due to temperature and liquid composition gradients, called thermosolutal convection or double diffusion, occurs in the mushy region and forms the complicated macrosegregation maps. Related to this phnomena, effects on the macrosegregation formation depending on the cooling condition and gravity values are described.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.4 no.3
• /
• pp.217-232
• /
• 1992

A Numerical study on two-dimensional laminar natural convection with and without surface radiation in fully or partially open square cavity was performed. The cavity has one vertical heated wall facing a vertical opening and two horizontal insulated walls. The pressure boundary condition was applied to the opening instead of the velocity boundary condition. The results of this study showed that the increase of partition length decreased the convective and the radiative Nusselt numbers. It was also found that the increase of wall emissivity decreased the convective Nusselt numbers but increased the radiative Nusselt numbers.

• East Asian mathematical journal
• /
• v.40 no.1
• /
• pp.109-118
• /
• 2024

This paper proposes a numerical method based on domain decomposition to find approximate solutions for one-dimensional convection-diffusion equations with Neumann boundary conditions. First, the equations are transformed into convection-diffusion equations with Dirichlet conditions. Second, the author introduces the Prediction/Correction Domain Decomposition (PCDD) method and estimates errors for the interface prediction scheme, interior scheme, and correction scheme using known error estimations. Finally, the author compares the PCDD algorithm with the fully explicit scheme (FES) and the fully implicit scheme (FIS) using three examples. In comparison to FES and FIS, the proposed PCDD algorithm demonstrates good results.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.1 no.1
• /
• pp.64-72
• /
• 1989

Laminar natural convection heat transfer from a horizontal heat exchanger tube with one infinitely long vertical plate fin has been studied by a finite-difference numerical procedure. In predicting convective heat transfer from a circular tube, the thermal boundary condition at solid fluid interface is usually assumed to be isothermal. However, in reality, the thermal boundary condition is not isothermal, and the tube has the thickness and the conductivity. So the temperature at the interface is not known a priori to the calculation. This problem has the conjugate phenomena which occur between the tube conduction and external natural convection, and between the fin conduction and external natural convection. Numerical results are obtained to determine the effects of the conductivity of solid wall and the thickness of tube wall on heat transfer. It is found that the conduction causes significant influence on the natural convection heat transfer at low K and high ${\delta}$.

• Journal of the Korean Society of Visualization
• /
• v.20 no.1
• /
• pp.29-37
• /
• 2022

We carried out a laboratory experiment about the thermo-fluidic characteristics of natural convection boundary layer over a vertical heated plate under constant heat flux condition. Particle image velocimetry has been applied to observe the surface convection velocity close to the vertical plate submerged in the water chamber with the condition of Ra = 7 × 10⁹ and Pr = 8.1. The velocity distributions indicate that the distinct stripe-like structures appears in the upstream (earlier transition region) and the distinct negative-positive and Λ(λ)-shaped flow structures in the downstream (mid-transition region). In addition, the temporal variation of spanwise and streamwise velocity is also presented.

• Communications of the Korean Mathematical Society
• /
• v.34 no.3
• /
• pp.1015-1027
• /
• 2019

A class of systems of singularly perturbed convection diffusion type equations with integral boundary conditions is considered. A numerical method based on a finite difference scheme on a Shishkin mesh is presented. The suggested method is of almost first order convergence. An error estimate is derived in the discrete maximum norm. Numerical examples are presented to validate the theoretical estimates.

• Journal of Mechanical Science and Technology
• /
• v.18 no.6
• /
• pp.1026-1035
• /
• 2004

An immersed boundary method for solving the Navier-Stokes and thermal energy equations is developed to compute the heat transfer over or inside the complex geometries in the Cartesian or cylindrical coordinates by introducing the momentum forcing, mass source/sink, and heat source/sink. The present method is based on the finite volume approach on a staggered mesh together with a fractional step method. The method of applying the momentum forcing and mass source/sink to satisfy the no-slip condition on the body surface is explained in detail in Kim, Kim and Choi (2001, Journal of Computational Physics). In this paper, the heat source/sink is introduced on the body surface or inside the body to satisfy the iso-thermal or iso-heat-flux condition on the immersed boundary. The present method is applied to three different problems : forced convection around a circular cylinder, mixed convection around a pair of circular cylinders, and forced convection around a main cylinder with a secondary small cylinder. The results show good agreements with those obtained by previous experiments and numerical simulations, verifying the accuracy of the present method.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.10
• /
• pp.1337-1345
• /
• 2001

Numerical investigation is made to study the effects of thermosolutal convection on the formation of macrosegregation in a Pb-Sn alloy solidification process in a two dimensional confined rectangluar mold. The basic equations are sovled using the Contrinum Model theory with the SIMPE algorithm during the solidification process. In addition, to track the liquid-solid interface with time variations, the moving boundary condition was adopted and moving irregular interface shapes were treated with the time-dependent, boundary-fitted coordinate system. As the temperature reduces from the liquidus to the solidus, the liquid concentration of Sn, the lighter constituent, increases. Then the buoyancy-driven flow due to temperature and liquid composition gradients occurs in the mushy region and forms the complicated macrosegregation maps. belated to this phenomena, effects on the macrosegregation formation depending on the cooling condition and gravity values are examined.

• Journal of IKEEE
• /
• v.13 no.3
• /
• pp.55-61
• /
• 2009

The heat transfer analysis applying finite element method has been carried out. Particularly, the convection boundary condition associated with the mixed boundary condition is numerically formulated by the Galerkin method analogous to the magnetic field problem. Also, the coupled electromagnet-thermal field analysis by the proposed heat transfer coefficient computation algorithm is executed to enhance the accuracy of solutions. Finally, the validity of the proposed results is verified by comparison with the measured ones.