J.R. Zabadal¹, B. Bodmann¹, V. G. Ribeiro², A. Silveira², S. Silveira²

CMES-Computer Modeling in Engineering & Sciences, Vol.103, No.4, pp. 215-227, 2014, DOI:10.3970/cmes.2014.103.215

Abstract This work presents a new analytical method to transform exact solutions of linear diffusion equations into exact ones for nonlinear advection-diffusion models. The proposed formulation, based on Bäcklund transformations, is employed to obtain velocity fields for the unsteady two-dimensional Helmholtz equation, starting from analytical solutions of a heat conduction type model. More >

J. Lin¹, W. Chen^1,2, C. S. Chen³, X. R. Jiang⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.6, pp. 485-505, 2013, DOI:10.3970/cmes.2013.094.485

Abstract To alleviate the difficulty of dense matrices resulting from the boundary knot method, the concept of the circulant matrix has been introduced to solve axi-symmetric Helmholtz problems. By placing the collocation points in a circular form on the surface of the boundary, the resulting matrix of the BKM has the block structure of a circulant matrix, which can be decomposed into a series of smaller matrices and solved efficiently. In particular, for the Helmholtz equation with high wave number, a large number of collocation points is required to achieve desired accuracy. In this paper, we present an efficient circulant boundary… More >

Chia-Cheng Tsai^1,2, D. L. Young³

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.2, pp. 175-205, 2013, DOI:10.3970/cmes.2013.094.175

Abstract It is well known that the method of fundamental solutions (MFS) is a numerical method of exponential convergence. In this study, the exponential convergence of the MFS is demonstrated by obtaining the eigensolutions of the Helmholtz equation. In the solution procedure, the sought solution is approximated by a superposition of the Helmholtz fundamental solutions and a system matrix is resulted after imposing the boundary condition. A golden section determinant search method is applied to the matrix for finding exponentially convergent eigenfrequencies. In addition, the least-squares method of fundamental solutions is applied for solving the corresponding eigenfunctions. In the solution procedure,… More >

N Pimprikar¹, J Teresa², D Roy^1,3, R M Vasu⁴, K Rajan⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.92, No.1, pp. 33-61, 2013, DOI:10.3970/cmes.2013.092.033

Abstract Nearly pollution-free solutions of the Helmholtz equation for k-values corresponding to visible light are demonstrated and verified through experimentally measured forward scattered intensity from an optical fiber. Numerically accurate solutions are, in particular, obtained through a novel reformulation of the H1 optimal Petrov-Galerkin weak form of the Helmholtz equation. Specifically, within a globally smooth polynomial reproducing framework, the compact and smooth test functions are so designed that their normal derivatives are zero everywhere on the local boundaries of their compact supports. This circumvents the need for a priori knowledge of the true solution on the support boundary and relieves the… More >

J.Y. Zhang¹, F.Z. Wang^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 141-154, 2012, DOI:10.3970/cmes.2012.088.141

Abstract The boundary knot method (BKM) is a kind of boundary-type meshless method, only boundary points are needed in the solution process. Since the BKM is mathematically simple and easy to implement, it is superior in dealing with Helmholtz problems with high wavenumbers and high dimensional problems. In this paper, we give an overview of the traditional BKM with collocation approach and provide three novel approaches for the BKM, as far as they are relevant for the other boundary-type techniques. The promising research directions are expected from an improved BKM aspect. More >

V. Mallardo¹, M.H. Aliabadi²

CMES-Computer Modeling in Engineering & Sciences, Vol.87, No.2, pp. 77-96, 2012, DOI:10.3970/cmes.2012.087.077

Abstract The present paper intends to couple the Fast Multipole Method (FMM) with the Boundary Element Method (BEM) in the 2D scalar wave propagation. The procedure is aimed at speeding the computation of the integrals involved in the governing Boundary Integral Equations (BIEs) on the basis of the distance between source point and integration element. There are three main contributions. First, the approach is of adaptive type in order to reduce the number of floating-point operations. Second, most integrals are evaluated analytically: the diagonal and off-diagonal terms of the H and G matrices by consolidated techniques, whereas the moment Mk by… More >

Wei-Ming Lee¹, Jeng-Tzong Chen^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.37, No.3, pp. 243-273, 2008, DOI:10.3970/cmes.2008.037.243

Abstract In this paper, a semi-analytical approach is proposed to solve the scattering problem of flexural waves and to determine dynamic moment concentration factors (DMCFs) in an infinite thin plate with multiple circular holes. The null-field integral formulation is employed in conjunction with degenerate kernels, tensor transformation and Fourier series. In the proposed direct formulation, all dynamic kernels of plate are expanded into degenerate forms and further the rotated degenerate kernels have been derived for the general exterior problem. By uniformly collocating points on the real boundary, a linear algebraic system is constructed. The results of dynamic moment concentration factors for… More >

C.C. Tsai¹, Y.C. Lin², D.L. Young^2,3, S.N. Atluri⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.2, pp. 103-114, 2006, DOI:10.3970/cmes.2006.016.103

Abstract In the applications of the method of fundamental solutions, locations of sources are treated either as variables or a priori known constants. In which, the former results in a nonlinear optimization problem and the other has to face the problem of locating sources. Theoretically, farther sources results in worse conditioning and better accuracy. In this paper, a practical procedure is provided to locate the sources for various time-independent operators, including Laplacian, Helmholtz operator, modified Helmholtz operator, and biharmonic operator. Wherein, the procedure is developed through systematic numerical experiments for relations among the accuracy, condition number, and source positions in different… More >

Zai You Yan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 81-90, 2006, DOI:10.3970/cmes.2006.013.081

Abstract Boundary element method in acoustics for Neumann boundary condition problems including sharp edges & corners is investigated. In previous acoustic boundary element method, acoustic pressure and normal velocity are the two variables at sharp edges & corners. However, the normal velocity at sharp edges & corners is discontinuous due to the indefinite normal vector. To avoid the indefinite normal vector and the discontinuous normal velocity at sharp edges & corners, normal vector of elemental node is defined and applied in the numerical implementation. Then the normal velocity is transformed to velocity which is unique even at sharp edges & corners.… More >

S.Yu. Reutskiy¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 27-40, 2006, DOI:10.3970/cmes.2006.012.027

Abstract In this paper a new boundary method for eigenproblems with the Helmholtz equation in simply and multiply connected domains is presented. The solution of an eigenvalue problem is reduced to a sequence of inhomogeneous problems with the differential operator studied. The method shows a high precision in simply and multiply connected domains and does not generate spurious eigenvalues. The results of the numerical experiments justifying the method are presented. More >