Xiao Zhao¹, Haitian Yang^1,2, Qiang Gao¹

CMES-Computer Modeling in Engineering & Sciences, Vol.74, No.2, pp. 139-160, 2011, DOI:10.3970/cmes.2011.074.139

Abstract A piecewised adaptive algorithm in the time domain is presented to solve the transient convection-diffusion heat transfer problem. By expanding all variables at a time interval, an initial and boundary value problem is decoupled into a series of recursive boundary value problems which can be solved by FEM or other well developed numerical schemes to deal with boundary value problems. A steady computing accuracy can be adaptively maintained via the power increase of the expansion, particularly when the step size varies in the whole computing process. Additionally for the nonlinear cases, there is no requirement of iteration and additional assumption… More >

D. Mirzaei¹, M. Dehghan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.72, No.3, pp. 185-210, 2011, DOI:10.3970/cmes.2011.072.185

Abstract The meshless local Petrov-Galerkin (MLPG) method with an efficient technique to deal with the time variable are used to solve the heat conduction problem in this paper. The MLPG is a meshless method which is (mostly) based on the moving least squares (MLS) scheme to approximate the trial space. In this paper the MLS is used for approximation in both time and space domains, and we avoid using the time difference discretization or Laplace transform method to overcome the time variable. The technique is applied for continuously nonhomogeneous functionally graded materials (FGM) in a finite strip and a hallow cylinder.… More >

Jong Keun Moon¹, Seung Jo Kim²

CMES-Computer Modeling in Engineering & Sciences, Vol.67, No.3, pp. 239-264, 2010, DOI:10.3970/cmes.2010.067.239

Abstract Large structural problems with high precision and complexity require a high-performance computation using the efficient parallel algorithm. The purpose of this paper is to present the parallel performance of thermal-structural coupled analysis tested on a parallel cluster system. In the coupled analysis, the heat transfer analysis is carried out, and then the structural analysis is performed based on temperature distribution. For the automatic and efficient connection of two parallel analysis modules, the several communication patterns were studied. The parallel performance was demonstrated for the sample and the real application problems, such as a laminated composite material by the DNS(Direct Numerical… More >

J.H. Tian^1,2, X. Han², S.Y. Long², G.Y. Sun², Y. Cao¹, G.Q. Xie³

CMES-Computer Modeling in Engineering & Sciences, Vol.63, No.2, pp. 101-116, 2010, DOI:10.3970/cmes.2010.063.101

Abstract A transient heat conduction analysis of the functionally graded material (FGM) plates has been investigated based on the hybrid numerical method (HNM). HNM combines the layer element method with the method of Fourier transforms and proves to be efficient and reliable. The FGM plates are infinite large and the material properties vary continuously through thickness. The transient heat source acted on the FGM plates. The temperature distribution of the FGM plates is obtained in different time and different position. Some useful results for transient heat conduction are shown in figures. Applications of HNM to transient heat conduction are firstly presented… More >

K. Han¹, Y. T. Feng¹, D. R. J. Owen¹

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.1, pp. 1-30, 2010, DOI:10.3970/cmes.2010.059.001

Abstract A computational framework is established for effective modelling of fluid-thermal-particle interactions. The numerical procedures comprise the Discrete Element Method for simulating particle dynamics; the Lattice Boltzmann Method for modelling the mass and velocity field of the fluid flow; and the Discrete Thermal Element Method and the Thermal Lattice Boltzmann Method for solving the temperature field. The coupling of the three fields is realised through hydrodynamic interaction force terms. Selected numerical examples are provided to illustrate the applicability of the proposed approach. More >

J. Sladek¹, V. Sladek¹, P.H. Wen², Y.C. Hon³

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.2, pp. 191-218, 2009, DOI:10.3970/cmes.2009.048.191

Abstract The meshless local Petrov-Galerkin (MLPG) method is used to solve the inverse heat conduction problem of predicting the distribution of the heat transfer coefficient on the boundary of 2-D and axisymmetric bodies. Using this method, nodes are randomly distributed over the numerical solution domain, and surrounding each of these nodes, a circular sub-domain is introduced. By choosing a unit step function as the test function, the local integral equations (LIE) on the boundaries of these sub-domains are derived. To eliminate the time variation in the governing equation, the Laplace transform technique is applied. The local integral equations are nonsingular and… More >

M. De Munck¹, D. Moens², W. Desmet³, D.Vandepitte³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 119-166, 2009, DOI:10.3970/cmes.2009.047.119

Abstract Deterministic simulation tools enable a very precise simulation of physical phenomena using numerical models. In many real life situations however, a deterministic analysis is not sufficient to assess the quality of a design. In a design stage, some physical properties of the model may not be determined yet. But even in a design ready for production, design tolerances and production inaccuracies introduce variability and uncertainty. In these cases, a non-deterministic analysis procedure is required, either using a probabilistic or a non-probabilistic approach. The authors developed an intelligent Kriging response surface based optimisation procedure that can be used in combination with… More >

J. Sladek¹, V. Sladek¹, P. Solek², C.L. Tan³, Ch. Zhang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 61-96, 2009, DOI:10.3970/cmes.2009.047.061

Abstract The meshless local Petrov-Galerkin (MLPG) method for transient linear thermoelastic analysis is presented. Orthotropic material properties are considered here. In uncoupled thermoelasticity, the temperature field is not influenced by displacements. Therefore, in the first step, the heat conduction equation is solved for the temperature distribution in the domain. The equations of motion are then solved with the inertial term considered. A Heaviside step function as the test functions is applied in the weak-form to derive local integral equations for solving two- and three-dimensional problems. Local integral equations are written on small sub-domains with circular or spherical shapes. They surround nodal… More >

Yueting Zhou¹, Xing Li², Dehao Yu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 191-222, 2009, DOI:10.3970/cmes.2009.043.191

Abstract The transient response of an orthotropic functionally graded strip with a partially insulated crack under convective heat transfer supply is considered. It is modeled there exists thermal resistant in the heat conduction through the crack region. The mixed boundary value problems of the temperature field and displacement field are reduced to a system of singular integral equations in Laplace domain. The expressions with high order asymptotic terms for the singular integral kernel are considered to improve the accuracy and efficiency. The numerical results present the effect of the material nonhomogeneous parameters, the orthotropic parameters and dimensionless thermal resistant on the… More >

E. J. Sellountos¹, A. Sequeira¹, D. Polyzos²

CMES-Computer Modeling in Engineering & Sciences, Vol.41, No.3, pp. 215-242, 2009, DOI:10.3970/cmes.2009.041.215

Abstract A Meshless Local Petrov-Galerkin (MLPG) method based on Local Boundary Integral Equation (LBIE) techniques is employed here for the solution of transient elastic problems with damping. The Radial Basis Functions (RBF) interpolation scheme is exploited for the meshless representation of displacements throughout the computational domain. On the intersections between the local domains and the global boundary, tractions are treated as independent variables via conventional boundary interpolation functions. The MLPG(LBIE)/RBF method is applied to both transient and steady-state Fourier transform elastodynamic domains. In both cases the LBIEs employ the simple elastostatic fundamental solution instead of the complicated time and frequency dependent… More >