Ridha Djebali^1,2, Bernard Pateyron², Mohamed El Ganaoui³

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.4, pp. 293-308, 2012, DOI:10.3970/cmes.2012.088.293

Abstract We present in this study a numerical investigation of unsteady two-dimensional natural convection of an electrically conducting fluid in a square cavity under an externally imposed magnetic field. A temperature gradient is applied between the two opposing side walls parallel to y-direction, while the floor and ceiling parallel to x-direction are adiabatic. The flow is characterized by the Rayleigh number Ra raged in 103-106, the Prandtl number Pr ranged in 0.01-10, the Hartman number Ha determined by the strength of the imposed magnetic field ranged in 0-100 and its tilting angle from x-axis ranging from 0 to 90 . The… More >

Taoufik Naffouti^1,2 and Ridha Djebali^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.3, pp. 211-228, 2012, DOI:10.3970/cmes.2012.088.211

Abstract This paper reports numerical results of natural convection flow evolving inside confined medium defined by two-dimensional square enclosure containing isothermal hot source placed asymmetrically at bottom wall. The sides-walls are isothermally cooled at a constant temperature; however the ceiling and the rest of bottom wall are insulated. The lattice Boltzmann method is used to solve the dimensionless governing equations with the associated boundary conditions. The flow is monitored by the Grashof number and the Prandtl number taken here 0.71. Numerical simulations are performed to study the effects of Grashof number ranging from 104 to 106, hot source length from 0.1… More >

A. Arefmanesh¹, M. Najafi², M. Nikfar³

CMES-Computer Modeling in Engineering & Sciences, Vol.69, No.2, pp. 91-118, 2010, DOI:10.3970/cmes.2010.069.091

Abstract Procuring a numerical solution through an application of the meshless local Petrov-Galerkin method (MLPG) on the fluid flow and mixed convection in a complex geometry cavity filled with a nanofluid is the scope of the present study. The cavity considered is a square enclosure having a lower temperature sliding lid at the top, a differentially higher temperature wavy wall at the bottom, and two thermally insulated walls on the sides. The nanofluid medium used is a water-based nanofluid, Al₂O₃-water with various volume fractions of its solid. To carry out the numerical simulations, the developed governing equations are determined in terms… More >

Wei-Chung Tien¹, Yue-Tzu Yang¹, Cha’o-Kuang Chen^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.5, pp. 447-462, 2012, DOI:10.3970/cmes.2012.085.447

Abstract Both buoyancy effects of thermal and mass diffusion in the natural convection flow about a vertical plate are considered in this paper. The non-linear coupled differential governing equations for velocity, temperature and concentration fields are solved by using the homotopy analysis method. Without the need of iteration, the obtained solution is in the form of an infinite power series which indicates those series have high accuracy when comparing it with other-generated by the traditional method. The impact of the Prandtl number, Schmidt number and the buoyancy parameter on the flow are widely discussed in detail. More >

Ahmad Shirzadi¹

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.1, pp. 45-64, 2012, DOI:10.3970/cmes.2012.085.045

Abstract This paper presents a meshless method based on the meshless local integral equation (LIE) method for solving the two-dimensional diffusion and diffusion-convection equations subject to a non-local condition. Suitable finite difference scheme is used to eliminate the time dependence of the problem. A weak formulation on local subdomains with employing the fundamental solution of the Laplace equation as test function transforms the resultant elliptic type equations into local integral equations. Then, the Moving Least Squares (MLS) approximation is employed for discretizing spatial variables. Two illustrative examples with exact solutions being used as benchmark solutions are presented to show the efficiency… More >

C.S. Nor Azwadi¹, N.I.N. Izual²

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.6, pp. 559-574, 2012, DOI:10.3970/cmes.2012.084.559

Abstract Natural convection in a differentially heated enclosure plays vital role in engineering applications such as nuclear reactor, electronic cooling technologies, roof ventilation, etc. The developed thermal flow patterns induced by the density difference are expected to be critically dependence on the inclination angles of the cavity. Hence, thermal and fluid flow pattern inside a differentially heated side enclosure walls with various inclination angles have been investigated numerically using the mesoscale lattice Boltzmann scheme. Three different dimensionless Rayleigh numbers were used, and a dimensionless Prandtl number of 0.71 was set to simulate the circulation of air in the system. It was… More >

Ridha Djebali^1,2, Mohamed ElGanaoui³, Taoufik Naffouti¹

CMES-Computer Modeling in Engineering & Sciences, Vol.84, No.6, pp. 499-527, 2012, DOI:10.3970/cmes.2012.084.499

Abstract A thermal lattice Boltzmann model for incompressible flow is developed and extended to investigate the natural convection flow in porous media under the effect of uniform magnetic field. The study shows that the flow behaviour is various parameters dependent. The Rayleigh number (Ra), Hartmann number (Ha), Darcy number (Da) and the medium inclination angle from the horizontal (Φ), the magnetic field orientation (ψ) and the medium porosity (ε) effects are carried out in wide ranges encountered in industrial and engineering applications. It was found that the flow and temperature patterns change significantly when varying these parameters. To confirm the accuracy… More >

Yiming Chen, Mingxu Yi, Chen Chen, Chunxiao Yu

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.6, pp. 639-654, 2012, DOI:10.3970/cmes.2012.083.639

Abstract In this paper, Bernstein polynomials method is proposed for the numerical solution of a class of space-time fractional convection-diffusion equation with variable coefficients. This method combines the definition of fractional derivatives with some properties of Bernstein polynomials and are dispersed the coefficients efficaciously. The main characteristic behind this method is that the original problem is translated into a Sylvester equation. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples show that the method is effective. More >

D. Ngo-Cong^1,2, N. Mai-Duy¹, W. Karunasena², T. Tran-Cong^1,3

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.3, pp. 275-310, 2012, DOI:10.3970/cmes.2012.083.275

Abstract In this paper, natural convection flows in concentric and eccentric annuli are studied using a new numerical method, namely local moving least square - one dimensional integrated radial basis function networks (LMLS-1D-IRBFN). The partition of unity method is used to incorporate the moving least square (MLS) and one dimensional-integrated radial basis function (1D-IRBFN) techniques in an approach that leads to sparse system matrices and offers a high level of accuracy as in the case of 1D-IRBFN method. The present method possesses a Kronecker-Delta function property which helps impose the essential boundary condition in an exact manner. The method is first… More >

Alfredo Nicolás¹, Blanca Bermúdez²

CMES-Computer Modeling in Engineering & Sciences, Vol.82, No.3&4, pp. 163-178, 2011, DOI:10.32604/cmes.2011.082.163

Abstract Mixed convection viscous incompressible fluid flows, under a gravitational system, in rectangular cavities are reported using the unsteady Boussinessq approximation in velocity-vorticity variables. The results are obtained using a numerical method based on a fixed point iterative process to solve the nonlinear elliptic system that results after time discretization; the iterative process leads to the solution of uncoupled, well-conditioned, symmetric linear elliptic problems for which efficient solvers exist regardless of the space discretization. Results with different aspect ratios A up to Grashof numbers Gr = 100000 and Reynolds numbers Re = 1000 for the lid driven cavity problem are reported.… More >