• Nuclear Engineering and Technology
• /
• 제41권3호
• /
• pp.279-286
• /
• 2009

Improving over a previous study [1], this paper provides a Monte Carlo method for the heat conduction analysis of problems with complicated geometry (such as a pebble with dispersed fuel particles). The method is based on the theoretical results of asymptotic analysis of neutron transport equation. The improved method uses an appropriate boundary layer correction (with extrapolation thickness) and a scaling factor, rendering the problem more diffusive and thus obtaining a heat conduction solution. Monte Carlo results are obtained for the randomly distributed fuel particles of a pebble, providing realistic temperature distributions (showing the kernel and graphite-matrix temperatures distinctly). The volumetric analytic solution commonly used in the literature is shown to predict lower temperatures than those of the Monte Carlo results provided in this paper.

• Journal of Advanced Marine Engineering and Technology
• /
• 제22권5호
• /
• pp.670-679
• /
• 1998

The solution of hyperbolic equation with wave nature has sharp discontinuties in the medium at the wave front. Difficulties encounted in the numrtical solution of such problem in clude among oth-ers numerical oscillation and the representation of sharp discontinuities with good resolution at the wave front. In this work inviscid Burgers equation and modified heat conduction equation is intro-duced as hyperboic equation. These equations are caculated by numerical methods(explicit method MacCormack method Total Variation Diminishing(TVD) method) along various Courant numbers and numerical solutions are compared with the exact analytic solution. For inviscid Burgers equa-tion TVD method remains stable and produces high resolution at sharp wave front but for modified heat Conduction equation MacCormack method is recommmanded as numerical technique.

• Structural Engineering and Mechanics
• /
• 제4권1호
• /
• pp.61-72
• /
• 1996

The time-stepping boundary element method has been so far applied by the authors to transient heat conduction in isotropic solids as well as in orthotropic solids. In this paper, attempt is made to extend the method to 2-D transient heat conduction in arbitrarily anisotropic solids. The resulting boundary integral equation is discretized by means of the boundary element with quadratic interpolation. The final system of equations thus obtained is solved by advancing the time step from the given initial state to the final state. Through numerical compuation of a few examples the potential usefulness of the proposed method is demonstrated.

• 설비공학논문집
• /
• 제1권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}$.

• 대한기계학회논문집
• /
• 제14권2호
• /
• pp.453-462
• /
• 1990

The interaction of forced convection-conduction with thermal radiation in laminar boundary layer over a plate fin is studied numerically. The analysis is based on complete solution whereby the heat conduction equation for the fin is solved simultaneously with the conservation equations for mass, momentum and energy in the fluid boundary layer adjacent to the fin. The fluid is a gray medium and diffusion(Rosseland) approximation is used to describe the radiative heat flux in the energy equation. The resulting boundary value problem are convection-conduction parameter N$_{c}$ and radiation-conduction parameter m, Prandtl number Pr. Numerical results are presented for gases with the Prandtl numbers of 0.7 & 5 with values of N$_{c}$ and M ranging from 0 to 10 respectively. The object of this study is to provide the first results on forced convection-radiation interaction in boundary layer flow over a semi-infinite flay plate which can be used for comparisons with future studies that will consider a more accurate expression for the radiative heat flux. The agreement of the results from the complete solution presented by E. M. Sparrow and those from this paper for the special case of M=0 is good. The overall rate of heat transfer from the fin considering radiative effect is higher than that from the fin neglecting radiative effect. The local heat transfer coefficient with radiative effect is higher than that without radiative effect. In the direction from tip to base, those coefficients decrease at first, attain minimum, and then increase. The larger values of N$_{c}$ M, Pr give rise to larger fin temperature variations and the fin temperature without radiative effect is always higher than that with radiative effect.

• Journal of Advanced Marine Engineering and Technology
• /
• 제31권7호
• /
• pp.833-841
• /
• 2007

The aim of this paper is to analyze the conjugate problems of heat conduction in solid walls coupled with laminar free convection flow adjacent to a vertical flat plate under boundary layer approximation. Using the similarity transformations the governing boundary layer equations for momentum and energy are reduced to a system of partial differential equations and then solved numerically using Finite Difference Method(FDM) known as the Keller-box scheme. Computed solutions to the governing equations are obtained for a wide range of non-dimensional parameters that are present in this problem, namely the coupling parameter P. the Prandtl number Pr and the heat generation parameter Q. The variations of the local heat transfer rate as well as the interface temperature and the friction along the plate and typical velocity and temperature profiles in the boundary layer are shown graphically. Numerical solutions have been consider for the Prandtl number Pr=0.70

• Structural Engineering and Mechanics
• /
• 제11권4호
• /
• pp.373-392
• /
• 2001

The main idea behind the paper is to present two alternative methods of homogenization of the heat conduction problem in composite materials, where the heat conductivity coefficients are assumed to be random variables. These two methods are the Monte-Carlo simulation (MCS) technique and the second order perturbation second probabilistic moment method, with its computational implementation known as the Stochastic Finite Element Method (SFEM). From the mathematical point of view, the deterministic homogenization method, being extended to probabilistic spaces, is based on the effective modules approach. Numerical results obtained in the paper allow to compare MCS against the SFEM and, on the other hand, to verify the sensitivity of effective heat conductivity probabilistic moments to the reinforcement ratio. These computational studies are provided in the range of up to fourth order probabilistic moments of effective conductivity coefficient and compared with probabilistic characteristics of the Voigt-Reuss bounds.

• 대한기계학회논문집B
• /
• 제33권3호
• /
• pp.170-177
• /
• 2009

In the present study, a finite element analysis of conjugate heat transfer problem inside a cavity with a heat-generating conducting body, where constant heat flux is generated, is conducted. A conduction heat transfer problem inside the solid body is automatically coupled with natural convection inside the cavity by using a finite element formulation. A finite element formulation based on SIMPLE type algorithm is adopted for the solution of the incompressible Navier-Stokes equations coupled with energy equation. The proposed algorithm is verified by solving the benchmark problem of conjugate heat transfer inside a cavity having a centered body. Then a conjugate natural heat transfer problem inside a cavity having a heat-generating conducting body with constant heat flux is solved and the effect of the Rayleigh number on the heat transfer characteristics inside a cavity is investigated.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제9권2호
• /
• pp.75-82
• /
• 2005

An inverse problem of transient heat conduction in a thin finite circular plate with the given temperature distribution on the interior surface of a thin circular plate being a function of both time and position has been solved with the help of integral transform technique and also determine the thermal deflection on the outer curved surface of a thin circular plate defined as $0\;{\leq}\;r\;{\leq}\;a,\;0\;{\leq}\;z\;{\leq}\;h$. The results, obtained in the series form in terms of Bessel's functions, are illustrated numerically.

• Journal of Mechanical Science and Technology
• /
• 제20권1호
• /
• pp.158-166
• /
• 2006

Extended Graetz problem in microtube is analyzed by using eigenfunction expansion to solve the energy equation. For the eigenvalue problem we applied the shooting method and Galerkin method. The hydrodynamically isothermal developed flow is assumed to enter the microtube with uniform temperature or uniform heat flux boundary condition. The effects of velocity and temperature jump boundary condition on the microtube wall, axial conduction and viscous dissipation are included. From the temperature field obtained, the local Nusselt number distributions on the tube wall are obtained as the dimensionless parameters (Peclet number, Knudsen number, Brinkman number) vary. The fully developed Nusselt number for each boundary condition is obtained also in terms of these parameters.