Jianxin Zhu, Zengsi Chen, Zheqi Shen

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.5, pp. 547-560, 2012, DOI:10.3970/cmes.2012.083.547

Abstract An acoustic waveguide with continuously varying sound speed is discussed in this paper. When the waveguide is open along the depth, the perfectly matched layer (PML) is used to terminate the infinite domain. Since the sound speed is gradually varied, the density is assumed as constant in each fluid layer. For this waveguide, it is shown that the mode relation is derived by using the differential transfer matrix method (DTMM). To solve leaky and PML modes, Newton's iteration is applied, and Chebyshev pseudospectral method is used for obtaining initial guesses. The solutions are with high More >

L. Dong¹, S. N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.2, pp. 183-220, 2012, DOI:10.32604/cmes.2012.083.183

Abstract In this paper, we develop T-Trefftz Voronoi Cell Finite Elements (VCF -EM-TTs) for micromechanical modeling of composite and porous materials. In addition to a homogenous matrix in each polygon-shaped element, three types of arbitrarily-shaped heterogeneities are considered in each element: an elastic inclusion, a rigid inclusion, or a void. In all of these three cases, an inter-element compatible displacement field is assumed along the element outer-boundary, and interior displacement fields in the matrix as well as in the inclusion are independently assumed as T-Trefftz trial functions. Characteristic lengths are used for each element to scale… More >

Minghui Wang¹, Musheng Wei², Shanrui Hu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.3&4, pp. 173-182, 2011, DOI:10.3970/cmes.2011.077.173

Abstract The mapping from the symmetric solution set to its independent parameter space is studied and an iterative method is proposed for the least-squares minimum-norm symmetric solution of AXB = E. Numerical results are reported that show the efficiency of the proposed methods. More >

G.A. Gravvanis¹, K.M. Giannoutakis²

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.4, pp. 305-330, 2011, DOI:10.3970/cmes.2011.071.305

Abstract In this paper we present parallel explicit approximate inverse matrix techniques for solving sparse linear systems on shared memory systems, which are derived using the finite element method for biharmonic equations in three space variables. Our approach for solving such equations is by considering the biharmonic equation as a coupled equation approach (pair of Poisson equation), using a FE approximation scheme, yielding an inner-outer iteration method. Additionally, parallel approximate inverse matrix algorithms are introduced for the efficient solution of sparse linear systems, based on an anti-diagonal computational approach that eliminates the data dependencies. Parallel explicit More >

Ying-Hsiu Shen¹, Chein-Shan Liu^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.2, pp. 157-178, 2011, DOI:10.3970/cmes.2011.071.157

Abstract When the local differential quadrature (LDQ) has been successfully applied to solve two-dimensional problems, the global method of DQ still has a problem by requiring to solve the inversions of ill-posed matrices. Previously, when one uses (n-1)th order polynomial test functions to determine the weighting coefficients with n grid points, the resultant n ×n Vandermonde matrix is highly ill-conditioned and its inversion is hard to solve. Now we use (m-1)th order polynomial test functions by n grid points that the size of Vandermonde matrix is m×n, of which m is much less than n. We More >

J.A. Ferreira^2,3, P. de Oliveira², P. Silva⁴, J.N. Murta⁵

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.1, pp. 1-14, 2011, DOI:10.3970/cmes.2011.071.001

Abstract Mathematical models to describe drug concentration profiles of topically administered drug in the anterior chamber aqueous humor have been proposed by several authors. The aim of this paper is to present a mathematical model to predict the drug concentration in the anterior chamber when a therapeutical contact lens with the drug entrapped in nanoparticles is used. More >

Chein-Shan Liu¹, Hong-Ki Hong¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.3, pp. 279-308, 2010, DOI:10.3970/cmes.2010.060.279

Abstract We propose novel algorithms to calculate the inverses of ill-conditioned matrices, which have broad engineering applications. The vector-form of the conjugate gradient method (CGM) is recast into a matrix-form, which is named as the matrix conjugate gradient method (MCGM). The MCGM is better than the CGM for finding the inverses of matrices. To treat the problems of inverting ill-conditioned matrices, we add a vector equation into the given matrix equation for obtaining the left-inversion of matrix (and a similar vector equation for the right-inversion) and thus we obtain an over-determined system. The resulting two modifications… More >

Chein-Shan Liu¹, Weichung Yeih², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.3, pp. 281-312, 2009, DOI:10.3970/cmes.2009.044.281

Abstract Here we develop a general purpose pre/post conditionerT, to solve an ill-posed system of linear equations,Ax=b. The conditionerTis obtained in the course of the solution of the Laplace equation, through a boundary-collocation Trefftz method, leading to:Ty=x, whereyis the vector of coefficients in the Trefftz expansion, andxis the boundary data at the discrete points on a unit circle. We show that the quality of the conditionerTis greatly enhanced by using multiple characteristic lengths (Multiple Length Scales) in the Trefftz expansion. We further show thatTcan be multiplicatively decomposed into a dilationTDand a rotationTR. For an odd-orderedA, we More >

Dong Ju Woo¹, Jin Oh Yang¹, Beom-Soo Kim¹, Seungsoo Lee¹, Jin Yeon Cho²

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.3, pp. 149-162, 2008, DOI:10.3970/cmes.2008.030.149

Abstract A new orphan-cell-free overset method is proposed to carry out the coupled analysis of overlapping finite element substructures. In the proposed overset method, the meshless MLS (Moving Least Squares) approximation is used to obtain the boundary data for the overlapped interface, whereas the Lagrange interpolation scheme has been commonly used in the conventional overset methods. The meshless character of MLS approximation makes it possible to overcome the problem of orphan-cell, which is often encountered in the conventional overset methods. Further, a new connectivity matrix solution procedure is developed to reduce the computational time in the 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 More >