Ji-Chuan Liu¹, Quan-Guo Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.93, No.3, pp. 203-220, 2013, DOI:10.3970/cmes.2013.093.203

Abstract In this paper, we propose an algorithm to solve a Cauchy problem of the Laplace equation in doubly connected domains for 2D and 3D cases in which the Cauchy data are given on the outer boundary. We want to seek a solution in the form of the single-layer potential and discrete it by parametrization to yield an ill-conditioned system of algebraic equations. Then we apply the Tikhonov regularization method to solve this ill-posed problem and obtain a stable numerical solution. Based on the regularization parameter chosen suitably by GCV criterion, the proposed method can get the approximate temperature and heat… More >

Chein-Shan Liu¹, Chung-Lun Kuo²

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.2, pp. 109-132, 2011, DOI:10.3970/cmes.2011.076.109

Abstract The Tikhonov method is a famous technique for regularizing ill-posed systems. In this theory a regularization parameter α needs to be determined. Based-on an invariant-manifold defined in the space of (x,t) and from the Tikhonov minimization functional, we can derive an optimal vector driven system of nonlinear ordinary differential equations (ODEs). In the Optimal Vector Driven Algorithm (OVDA), the optimal regularization parameter αk is presented in the iterative solution of x, which means that a dynamical Tikhonov regularization method is involved in the solution of nonlinear ill-posed problem. The OVDA is an extension of the Landweber-Scherzer iterative algorithm. Numerical examples… More >

C. Shi¹, C. Wang¹, T. Wei^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.86, No.2, pp. 71-92, 2012, DOI:10.3970/cmes.2012.086.071

Abstract In this paper, we consider a typical ill-posed inverse heat source problem, that is, we determine a space-dependent heat source term in a multi-dimensional heat equation from a pair of Cauchy data on a part of boundary. By a simple transformation, the inverse heat source problem is changed into a Cauchy problem of a homogenous heat conduction equation. We use the method of fundamental solutions (MFS) coupled with the Tikhonov regularization technique to solve the ill-conditioned linear system of equations resulted from the MFS discretization. The generalized cross-validation rule for determining the regularization parameter is used. Numerical results for four… More >

Roman Chapko¹, B. Tomas Johansson²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 105-128, 2012, DOI:10.3970/cmes.2012.085.105

Abstract We consider a Cauchy problem for the Laplace equation in a 3-dimen -sional semi-infinite domain that contains a bounded inclusion. The canonical situation is the upper half-space in I\tmspace -.1667em R3 containing a bounded smooth domain. The function value of the solution is specified throughout the plane bounding the upper half-space, and the normal derivative is given only on a finite portion of this plane. The aim is to reconstruct the solution on the surface of the bounded inclusion. This is a generalisation of the situation in Chapko and Johansson (2008) to three-dimensions and with Cauchy data only partially given.… 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 of the MCGM, namely the… More >

Leevan Ling¹, Tomoya Takeuchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.29, No.1, pp. 45-54, 2008, DOI:10.3970/cmes.2008.029.045

Abstract The method of fundamental solutions is coupled with the boundary control technique to solve the Cauchy problems of the Laplace Equations. The main idea of the proposed method is to solve a sequence of direct problems instead of solving the inverse problem directly. In particular, we use a boundary control technique to obtain an approximation of the missing Dirichlet boundary data; the Tikhonov regularization technique and the L-curve method are employed to achieve such goal stably. Once the boundary data on the whole boundary are known, the numerical solution to the Cauchy problem can be obtained by solving a direct… More >