Nicolas Ali Libre^1,2, Arezoo Emdadi², Edward J. Kansa^3,4, Mohammad Shekarchi², Mohammad Rahimian²

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.3, pp. 263-284, 2008, DOI:10.3970/cmes.2008.038.263

Abstract We present a wavelet based adaptive scheme and investigate the efficiency of this scheme for solving nearly singular potential PDEs. Multiresolution wavelet analysis (MRWA) provides a firm mathematical foundation by projecting the solution of PDE onto a nested sequence of approximation spaces. The wavelet coefficients then were used as an estimation of the sensible regions for node adaptation. The proposed adaptation scheme requires negligible calculation time due to the existence of the fast Discrete Wavelet Transform (DWT). Certain aspects of the proposed adaptive scheme are discussed through numerical examples. It has been shown that the proposed adaptive scheme can detect… More >

Zhi-Zhong Yan^1,2, Yue-Sheng Wang^1,3, Chuanzeng Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.1, pp. 59-88, 2008, DOI:10.3970/cmes.2008.038.059

Abstract In this paper, a numerical method based on the wavelet theory is developed for calculating band structures of 2D phononic crystals consisting of general anisotropic materials. After systematical consideration of the appropriate choice of wavelets, two types of wavelets, the Haar wavelet and Biorthogonal wavelet, are selected. Combined with the supercell technique, the developed method can be then applied to compute the band structures of phononic crystals with point or line defects. We illustrate the advantages of the method both mathematically and numerically. Particularly some representative numerical examples are presented for various material combinations (solid-solid, solid-fluid and fluid-fluid) with complex… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.1, pp. 25-38, 2008, DOI:10.3970/cmes.2008.038.025

Abstract In this paper a complex matrix operator and a complex field are used to express the Maxwell equations, of which the complex field embraces all field variables and the matrix operator embraces the time and space differential operators. By left applying the operator on the complex field one can get all the four Maxwell equations, which are usually expressed by the vector form. The new formulation matches the Lorenz gauge condition, and its mathematical advantage is that it can incorporate the Maxwell equations into a single equation. The introduction of four-potential is possible only under the Lorenz gauge. In terms… More >

Sirajul Haq¹, Siraj-ul-Islam², Arshed Ali³

CMES-Computer Modeling in Engineering & Sciences, Vol.38, No.1, pp. 1-24, 2008, DOI:10.3970/cmes.2008.038.001

Abstract In this paper we propose a meshfree technique for the numerical solution of the modified equal width wave (MEW) equation. Combination of collocation method using the radial basis functions (RBFs) with first order accurate forward difference approximation is employed for obtaining meshfree solution of the problem. Different types of RBFs are used for this purpose. Performance of the proposed method is successfully tested in terms of various error norms. In the case of non-availability of exact solution, performance of the new method is compared with the results obtained from the existing methods. Propagation of a solitary wave, interaction of two… More >

Shueei-Muh Lin^1,3, Sen-Yung Lee², Cheng-Chuan Tsai², Chien-Wi Chen²,Wen-Rong Wang³, Jenn-Fa Lee³

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.1, pp. 33-54, 2008, DOI:10.3970/cmes.2008.034.033

Abstract The conventional theories for circulation of arteries are emphasized on fluid behavior or some simplified models for experimental utility. In this study, a new mathematical theory is proposed to describe the wave propagation through the elastic tube filled with viscous and incompressible fluid. The radial, longitudinal and flexural vibrations of a tube wall are introduced simultaneously. Meanwhile, the linearlized momentum and continuity equations of tube flow field are expressed in the integral form. Based on these considerations, three wave modes are obtained simultaneously. These wave modes are the flexural, Young and Lamb modes, respectively. The characteristics of these modes are… More >

Delfim Soares Jr. ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.2, pp. 57-68, 2008, DOI:10.3970/cmes.2008.032.057

Abstract In this work, a time-domain finite element -- boundary element iterative coupling technique is presented in order to analyse electromagnetic scattering from two-dimensional inhomogeneous objects. In the iterative coupling algorithm, the domain of the original problem is subdivided into sub-domains and each sub-domain is analysed independently (as an uncoupled model) taking into account successive renewals of variables at common interfaces. In order to improve the effectiveness of the iterative coupling approach, the evaluation of an optimised relaxation parameter is discussed, taking into account the minimisation of a square error functional. The algorithm that arises is efficient, accurate and flexible. At… More >

A.I. Abreu¹, W.J. Mansur¹, D. Soares Jr^1,2, J.A.M. Carrer³

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 123-132, 2008, DOI:10.3970/cmes.2008.030.123

Abstract This paper is concerned with the numerical computation of space derivatives of a time-domain (TD-) Boundary Element Method (BEM) formulation for the analysis of scalar wave propagation problems. In the present formulation, the Convolution Quadrature Method (CQM) is adopted, i.e., the basic integral equation of the TD-BEM is numerically substituted by a quadrature formula, whose weights are computed using the Laplace transform of the fundamental solution and a linear multi-step method. In order to numerically compute space derivatives, the present work properly transforms the quadrature weights of the CQM-BEM, adopting the so-called Complex-Variable-Differentiation Method (CVDM). Numerical examples are presented at… More >

D. Kourounis¹, L.N. Gergidis¹, A. Charalambopoulos¹

CMES-Computer Modeling in Engineering & Sciences, Vol.28, No.3, pp. 185-202, 2008, DOI:10.3970/cmes.2008.028.185

Abstract The sensitivity of analytical solutions of the direct acoustic scattering problem in prolate spheroidal geometry on the wavenumber and shape, is extensively investigated in this work. Using the well known Vekua transformation and the complete set of radiating "outwards'' eigensolutions of the Helmholtz equation, introduced in our previous work ([Charalambopoulos and Dassios(2002)], [Gergidis, Kourounis, Mavratzas, and Charalambopoulos (2007)]), the scattered field is expanded in terms of it, detouring so the standard spheroidal wave functions along with their inherent numerical deficiencies. An approach is employed for the determination of the expansion coefficients, which is optimal in the sense, that minimizes the… More >

Mira Mitra¹, S. Gopalakrishnan²

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 125-136, 2008, DOI:10.3970/cmes.2008.027.125

Abstract In this paper, the wave characteristics, namely, the spectrum and dispersion relations of multi-wall carbon nanotubes (MWNTs) are studied. The MWNTs are modeled as multiple thin shells coupled through van der Waals force. Each wall of the MWNT has three displacements, i.e, axial, circumferential and radial with variation along the axial and circumferential directions. The wave characteristics are obtained by transforming the governing differential wave equations to frequency domain via Fourier transform. This transformation is first done in time using fast Fourier transform (FFT) and then in one spatial dimension using Fourier series. These transformed equations are solved by posing… More >

T. Kar, M.L. Munjal¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.1&2, pp. 37-48, 2008, DOI:10.3970/cmes.2008.027.037

Abstract In the present work, a wedge structure made of absorbing panels has been analyzed by making use of the matrizant analysis with the help of the Boundary-Condition-Transfer (BCT) algorithm. The rectangular panel wedge, as it is called in this manuscript, is simple in geometry. The theoretical model, based on the plane wave acoustical coupling between multiple interacting ducts of variable cross sectional area, is applied to predict the pressure reflection coefficient of the present wedge configuration. Bulk reaction and hence wave propagation in the wedge material has been assumed in the proposed model. An asymptotic solution using the Peano-Baker series… More >