LiviuMarin ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.37, No.3, pp. 203-242, 2008, DOI:10.3970/cmes.2008.037.203

Abstract In this paper, a meshless method for the stable solution of direct and inverse problems associated with the two-dimensional Laplace equation in the presence of boundary singularities and noisy boundary data is proposed. The governing equation and boundary conditions are discretized by the method of fundamental solutions (MFS), whilst the existence of the boundary singularity is taken into account by subtracting from the original MFS solution the corresponding singular solutions, as given by the asymptotic expansion of the solution near the singular point. However, even in the case when the boundary singularity is accounted for, More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.36, No.3, pp. 261-286, 2008, DOI:10.3970/cmes.2008.036.261

Abstract The inverse Sturm-Liouville problem finds its applications in the identification of mechanical properties and/or geometrical configurations of a vibrating continuous medium; however, this problem is hard to solve, either theoretically or numerically. Previously, Liu (2008a) has constructed a Lie-group shooting method to determine the eigenvalues, and the corresponding eigenfunctions, for the direct Sturm-Liouville problem. In this study, we are concerned with solving the inverse Sturm-Liouville problem, by developing a Lie-group of SL(2,R) to construct nonlinear algebraic equations (NAEs), when discrete eigenvalues are specified. Our purpose here is to use these NAEs to solve the unknown function More >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 99-122, 2008, DOI:10.3970/cmes.2008.030.099

Abstract The application of the method of fundamental solutions (MFS) to inverse boundary value problems associated with the steady-state heat conduction in isotropic media in the presence of sources, i.e. the Poisson equation, is investigated in this paper. Based on the approach of Alves and Chen (2005), these problems are solved in two steps, namely by finding first an approximate particular solution of the Poisson equation and then the numerical solution of the resulting inverse boundary value problem for the Laplace equation. The resulting MFS discretised system of equations is ill-conditioned and hence it is solved More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.2, pp. 65-76, 2008, DOI:10.3970/cmes.2008.030.065

Abstract Trefftz method (TM) is one of widely used meshless numerical methods in elliptic type boundary value problems, of which the approximate solution is expressed as a linear combination of T-complete bases, and the unknown coefficients are determined from boundary conditions by solving a linear equations system. However, the accuracy of TM is severely limited by its ill-conditioning. This paper is a continuation of the work of Liu (2007a). The collocation TM is modified and applied to the direct and inverse problems of biharmonic equation in a simply connected plane domain. Due to its well-conditioning of More >

S. D. Harris¹, R. Mustata², L. Elliott², D. B. Ingham², D. Lesnic²

CMES-Computer Modeling in Engineering & Sciences, Vol.25, No.2, pp. 69-80, 2008, DOI:10.3970/cmes.2008.025.069

Abstract Two homogeneous anisotropic materials are butted together to form a contact surface within a single composite material (the specimen). An inverse boundary element method (BEM) is developed to determine the components of the hydraulic conductivity tensor of each material and the position of the contact surface. A steady state flow is forced through the specimen by the application of a constant pressure differential on its opposite faces. Experimental measurements (simulated) of pressure and average hydraulic flux at exposed boundaries are then used in a modified least squares functional. This functional minimises the gap between the More >

Liviu Marin¹, Henry Power¹, Richard W. Bowtell², Clemente Cobos Sanchez², Adib A. Becker¹, Paul Glover²，Arthur Jones¹

CMES-Computer Modeling in Engineering & Sciences, Vol.23, No.3, pp. 149-174, 2008, DOI:10.3970/cmes.2008.023.149

Abstract We investigate the reconstruction of a divergence-free surface current distribution from knowledge of the magnetic flux density in a prescribed region of interest in the framework of static electromagnetism. This inverse problem is motivated by the design of gradient coils for use in magnetic resonance imaging (MRI) and is formulated using its corresponding integral representation according to potential theory. A novel boundary element method (BEM) which employs linear interpolation on quadratic surfaces and also satisfies the continuity equation for the current density, i.e. a divergence-free BEM, is presented. Since the discretised BEM system is ill-posed More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.8, No.1, pp. 1-16, 2008, DOI:10.3970/cmc.2008.008.001

Abstract In the present work we study numerical computations of inverse thermal stress problems. The unknown boundary conditions of an elastically deformable heat conducting rod are not given a priori and are not allowed to measure directly, because the boundary may be not accessible to measure. However, an internal measurement of temperature is available. We treat this inverse problem by using a semi-discretization technique, of which the time domain is divided into many sub-intervals and the physical quantities are discretized at these node points of discrete times. Then the resulting ordinary differential equations in the discretized More >

Chein-Shan Liu^1,2

CMC-Computers, Materials & Continua, Vol.7, No.2, pp. 81-96, 2008, DOI:10.3970/cmc.2008.007.081

Abstract We consider an inverse problem for estimating an unknown heat conductivity parameter α(t) in a heat conduction equation T_t(x,t) = α(t)T_xx(x,t) with the aid of an extra measurement of temperature gradient on boundary. Basing on an establishment of the one-step Lie-group elements G(r) and G(l) for the semi-discretization of heat conduction equation in time domain, we can derive algebraic equations from G(r) = G(l). The new method, namely the Lie-group estimation method (LGEM), is examined through numerical examples to convince that it is highly accurate and efficient; the maximum estimation error is smaller than 10^-5 for smooth parameter and More >

Cheng-Hung Huang^1,2, Chih-Chun Shih¹

CMES-Computer Modeling in Engineering & Sciences, Vol.21, No.3, pp. 239-254, 2007, DOI:10.3970/cmes.2007.021.239

Abstract An inverse forced vibration problem, based on the Conjugate Gradient Method (CGM), (or the iterative regularization method), is examined in this study to estimate simultaneously the unknown time-dependent applied force and moment for an Euler-Bernoulli beam by utilizing the simulated beam displacement measurements. The accuracy of this inverse problem is examined by using the simulated exact and inexact displacement measurements. The numerical experiments are performed to test the validity of the present algorithm by using different types of applied force and moment, sensor locations and measurement errors. Results show that excellent estimations on the applied More >

Leevan Ling¹, Tomoya Takeuchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.2, pp. 99-110, 2007, DOI:10.3970/cmes.2007.020.099

Abstract We aim to identify the unknown source locations in a two-dimensional heat equation from scattered measurements. In [Inverse Problems, 22(4):1289--1305, 2006], we proposed a numerical procedure that identifies the unknown source locations of 2D heat equation solely based on three measurement points. Due to the nonlinearity and complexity of the problem, the quality of the resulting estimations is often poor especially when the number of unknown is large. In this paper, we purpose a linear refinement scheme that takes the outputs of the existing nonlinear algorithm as initial guesses and iteratively improves on the accuracy More >