E. Majchrzak¹, B. Mochnacki², A. L. Greer³, J. S. Suchy⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.2, pp. 131-146, 2009, DOI:10.3970/cmes.2009.041.131

Abstract Multi-layered thin metal film subjected to a short-pulse laser heating is considered. Mathematical description of the process discussed bases on the equation in which there appear the relaxation time and the thermalization time (dual-phase-lag-model). In this study we develop a three level implicit finite difference scheme for numerical modelling of heat transfer in non-homogeneous metal film. At the interfaces an ideal contact between successive layers is assumed. At the stage of computations a solution of only one three-diagonal linear system corresponds to transition from time t to t + Δt. The mathematical model, numerical algorithm and examples More >

Donghuan Liu¹, Xiaoping Zheng^1,2, Yinghua Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.3, pp. 263-300, 2009, DOI:10.3970/cmes.2009.039.263

Abstract A discontinuous Galerkin (DG) finite element method for the heat conduction problems with local high gradient and thermal contact resistance is presented. The DG formulation is constructed by employing the stabilization term and the Bassi-Rebay numerical flux term. The stabilization term is defined by a penalization of the temperature jump at the interface. By eliminating the penalization term of the temperature jump in the region of local high gradient and imperfect contact interfaces, the present DG method is applied to solve problems involving local high gradient and thermal contact resistance where the numerical flux is… More >

Ahmed Mezrhab¹, Hassan Naji²

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.3, pp. 283-296, 2009, DOI:10.3970/fdmp.2009.005.283

Abstract The objective of this paper is to assess the effectiveness of the coupled Lattice Boltzmann Equation (LBE) and finite volume method strategy for the simulation of the interaction between thermal radiation and laminar natural convection in a differentially heated square cavity. The vertical walls of the cavity are adiabatic, while its top and bottom walls are cold and hot, respectively. The air velocity is determined by the lattice Boltzmann equation and the energy equation is discretized by using a finite volume method. The resulting systems of discretized equations have been solved by an iterative procedure More >

R. Djebali^1,2, M. El Ganaoui^2,3, H. Sammouda¹, R. Bennacer⁴

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.3, pp. 261-282, 2009, DOI:10.3970/fdmp.2009.005.261

Abstract The simplified thermal lattice Boltzmann model (STLBM) developed by Peng, Shu and Chew (2003) is used in this work to simulate low-Rayleigh-number natural convection in a heated rectangular cavity on a uniform grid. It is shown how by resorting to the double populations approach both hydrodynamic and thermal fields can be effectively simulated. Furthermore, a general benchmark is carried out to account for the effect of different parameters in relatively wide ranges. Results are compared with previous works available in the literature. More >

M.R. Arab^1,2, E. Semma³, B. Pateyron¹, M. El Ganaoui¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.2, pp. 161-176, 2009, DOI:10.3970/fdmp.2009.005.161

Abstract In this work, flows in porous media are simulated by using a Lattice Boltzmann Method (LBM). A model D2Q9 with a single collision operator is proposed. This method is applied on 2D digital images obtained by a Scanning Electron Microscope technique (SEM), and followed by a special treatment in order to obtain an image of synthesis that is finally read by the numerical code. The first results tested on two-dimensional configurations show the reliability of this strategy in simulating with a good accuracy phenomena of heat and mass transport. The numerical study is extended to More >

Aouachria Z¹

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.2, pp. 137-148, 2009, DOI:10.3970/fdmp.2009.005.137

Abstract This work treats heat and mass transfer by natural convection along a vertical wall in porous media imbibed by fluid, using an integral method. The problem governing parameters are the buoyancy ratio,N, and the Lewis number, Le The results for the local Nusselt and Sherwood numbers are presented for a large range of these parameters. The concentration and thermal boundary layer thickness are also determined. We observe that our results are in good agreement with those obtained by Bejan and Khair (1985). More >

M. El Alami¹, E. A. Semma^2,3, M. Najam¹, R. Boutarfa²

FDMP-Fluid Dynamics & Materials Processing, Vol.5, No.1, pp. 23-36, 2009, DOI:10.3970/fdmp.2009.005.023

Abstract This study is devoted to the investigation of natural convection in a two dimensional horizontal channel with rectangular heated blocks at the bottom. The aspect ratio of the channel is A = L'/H' = 5. The blocks are heated with a constant temperature while the upper plane of the channel is cold. The governing equations are solved using a finite volumes method and the SIMPLEC algorithm is used for the treatment of the pressure-velocity coupling. Special emphasis is given to detail the effect of the Rayleigh number and blocks height on the heat transfer and the More >

W.F.Yuan¹, K.H.Tan ¹

CMC-Computers, Materials & Continua, Vol.13, No.1, pp. 49-62, 2009, DOI:10.3970/cmc.2009.013.049

Abstract A simple cellular automaton (CA) scheme is proposed to simulate heat conduction in anisotropic domains. The CA is built on random nodes rather than an irregular grid. The local rule used in the CA is defined by physical concepts instead of differential equations. The accuracy of the proposed approach is verified by classical examples. As an application of the proposed method, the CA approach is incorporated into fibre model which is widely used in finite element analysis to calculate the temperature distribution on the cross-section of composite beams. Numerical examples demonstrate that the proposed scheme More >

Chih-Ping Wu^1,2, Shao-En Huang²

CMC-Computers, Materials & Continua, Vol.12, No.3, pp. 251-282, 2009, DOI:10.3970/cmc.2009.012.251

Abstract A modified Pagano method is developed for the three-dimensional (3D) coupled analysis of simply-supported, doubly curved functionally graded (FG) piezo-thermo-elastic shells under thermal loads. Four different loading conditions, applied on the lateral surfaces of the shells, are considered. The material properties of FG shells are regarded as heterogeneous through the thickness coordinate, and then specified to obey an exponent-law dependent on this. The Pagano method, conventionally used for the analysis of multilayered composite elastic plates/shells, is modified to be feasible for the present analysis of FG piezo-thermo-elastic plates/shells. The modifications include that a displacement-based formulation More >

Silieti, M.¹, Divo, E.², Kassab, A.J.¹

CMC-Computers, Materials & Continua, Vol.12, No.2, pp. 121-144, 2009, DOI:10.3970/cmc.2009.012.121

Abstract A hybrid singularity superposition/boundary element-based inverse problem method for the reconstruction of multi-dimensional heat flux distributions is developed. Cauchy conditions are imposed at exposed surfaces that are readily reached for measurements while convective boundary conditions are unknown at surfaces that are not amenable to measurements such as the walls of the cooling holes. The purpose of the inverse analysis is to determine the heat flux distribution along cooling hole surfaces. This is accomplished in an iterative process by distributing a set of singularities (sinks) inside the physical boundaries of the cooling hole (usually along cooling… More >