R. Avila¹, F.J. Solorio¹

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 177-202, 2009, DOI:10.3970/cmes.2009.044.177

Abstract Heat transfer processes involving phase change either, solidification or melting, appear frequently in nature and in industrial applications. In this paper the convective patterns that arise from a 2D shear driven annular flow (without and with melting), are presented. The convective annular flow with radial gravity can be considered as a simplified model of the atmospheric flow in the terrestrial equatorial plane (bounded by the warm surface of the Earth and the cold tropopause). The governing equations have been numerically solved by the Spectral Element Method. The numerical results reported in this paper, for the cases without melting (at two… More >

Delfim Soares Jr.¹

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 157-176, 2009, DOI:10.3970/cmes.2009.044.157

Abstract The numerical modelling of electromagnetic waves by finite element - boundary element coupling procedures is discussed here, taking into account time-domain approaches. In this study, the global model is divided into different sub-domains and each sub-domain is analysed independently and explicitly at each time-step of the analysis: the interaction between the different sub-domains of the global model is accomplished by interface procedures. A multi-level time-step algorithm is considered in order to improve the flexibility, accuracy and stability (especially when conditionally stable time-marching procedures are employed) of the coupled analysis. At the end of the paper, numerical examples are presented, illustrating… More >

Chih-Fung Ho¹, Cheng Chang¹, Kuen-Hau Lin¹, Chao-An Lin^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 137-156, 2009, DOI:10.3970/cmes.2009.044.137

Abstract Consistent formulations of 2D and 3D pressure and velocity boundary conditions along both the stationary and non-stationary plane wall and corner for lattice Boltzmann simulations are proposed. The unknown distribution functions are made function of local known distribution functions and correctors, where the correctors at the boundary nodes are obtained directly from the definitions of density and momentum. This boundary condition can be easily implemented on the wall and corner boundary using the same formulation. Discrete macroscopic equation is also derived for steady fully developed channel flow to assess the effect of the boundary condition on the solutions, where the… More >

Sirajul Haq^1,2, Siraj-Ul-Islam³, Marjan Uddin^1,4

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.2, pp. 115-136, 2009, DOI:10.3970/cmes.2009.044.115

Abstract A mesh free method for the numerical solution of the nonlinear Schrodinger (NLS) and coupled nonlinear Schrodinger (CNLS) equation is implemented. The presented method uses a set of scattered nodes within the problem domain as well as on the boundaries of the domain along with approximating functions known as radial basis functions (RBFs). The set of scattered nodes do not form a mesh, means that no information of relationship between the nodes is needed. Error norms L₂, L_∞ are used to estimate accuracy of the method. Stability analysis of the method is given to demonstrate its practical applicability. More >

Rui Wang¹, Jianxun Zhang¹, Hisashi Serizawa², Hidekazu Murakawa²

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 97-114, 2009, DOI:10.3970/cmes.2009.044.097

Abstract In this paper, the inherent deformation of fillet welded thin plate T-joints is studied. The prediction procedure of inherent deformation consists of three parts: part one, a three dimensional (3D) thermo-elastic-plastic analysis using an in house finite element (FE) code of interactive substructure method (ISM) is utilized to obtain the welding distortions; part two, corresponding experiments are carried out to verify the computational results of ISM; part three, using the verified computational results, the inverse analysis is utilized to evaluate the welding inherent deformation. Based on the results in this study, an inherent deformations database of fillet welded thin plate… More >

R. Vertnik¹, B. Šarler²

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 65-96, 2009, DOI:10.3970/cmes.2009.044.065

Abstract The application of the mesh-free Local Radial Basis Function Collocation Method (LRBFCM) in solution of incompressible turbulent flow is explored in this paper. The turbulent flow equations are described by the low - Re number k-emodel with Jones and Launder [Jones and Launder (1971)] closure coefficients. The involved velocity, pressure, turbulent kinetic energy and dissipation fields are represented on overlapping 5-noded sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBF's. The velocity, turbulent kinetic energy and dissipation equations are solved through… More >

S Reaz Ahmed¹, S K Deb Nath¹

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 35-64, 2009, DOI:10.3970/cmes.2009.044.035

Abstract The paper presents a simplified analysis of stresses and deformations at critical sections of a tire-tread. Displacement potential formulation is used in conjunction with the finite-difference method to model the present contact problem. The solution of the problem is obtained for two limiting cases of the contact boundary - one allows the lateral slippage and the other conforms to the no-slip condition along the lateral direction. The influential effects of tire material and tread aspect-ratio are discussed. The reliability and accuracy of the solution is also discussed in light of comparison made with the usual computational approach. More >

L.M.J.S. Dinis¹, R.M. Natal Jorge², J. Belinha³

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 1-34, 2009, DOI:10.3970/cmes.2009.044.001

Abstract The Natural Neighbour Radial Point Interpolation Method (NNRPIM) is extended to the large deformation analysis of non-linear elastic structures. The nodal connectivity in the NNRPIM is enforced using the Natural Neighbour concept. After the Voronoï diagram construction of the unstructured nodal mesh, which discretize the problem domain, small cells are created, the "influence-cells". These cells are in fact influence-domains entirely nodal dependent. The Delaunay triangles are used to create a node-depending background mesh used in the numerical integration of the NNRPIM interpolation functions. The NNRPIM interpolation functions, used in the Galerkin weak form, are constructed with the Radial Point Interpolators.… More >

Q.W. Ma^1,2, J.T. Zhou¹

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 277-304, 2009, DOI:10.3970/cmes.2009.043.277

Abstract Following our previous work, the Meshless Local Petrov-Galerin me -thod based on Rankine source solution (MLPG_R) will be extended in this paper to deal with breaking waves. For this purpose, the governing equation for pressure is improved and a new technique called Mixed Particle Number Density and Auxiliary Function Method (MPAM) is suggested for identifying the free surface particles. Due to complexity of the problem, two dimensional (2D) breaking waves are only concerned here. Various cases are investigated and some numerical results are compared with experimental data available in literature to show the newly developed method is robust. More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 253-276, 2009, DOI:10.3970/cmes.2009.043.253

Abstract Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a₀ + a₁x + ... + a_nxⁿ expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a new polynomial interpolation… More >