Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.21, No.1, pp. 17-40, 2011, DOI:10.3970/cmc.2011.021.017

Abstract Only the left-boundary data of temperature and heat flux are used to estimate an unknown parameter function α(x) in T_t(x,t) = ∂(α(x)T_x)/∂x + h(x,t), as well as to recover the right-boundary data. When α(x) is given the above problem is a well-known inverse heat conduction problem (IHCP). This paper solves a mixed-type inverse problem as a combination of the IHCP and the problem of parameter identification, without needing to assume a function form of α(x) a priori, and without measuring extra data as those used by other methods. We use the one-step Lie-Group Adaptive Method (LGAM) for the semi-discretizations of… More >

Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.19, No.3, pp. 285-314, 2010, DOI:10.3970/cmc.2010.019.285

Abstract In this article, we propose a new numerical approach for solving these multi-dimensional nonlinear and nonhomogeneous backward heat conduction problems (BHCPs). A fictitious time t is employed to transform the dependent variable u(x, y, z, t) into a new one by (1+t)u(x, y, z, t)=: v(x, y, z, t, t), such that the original nonlinear and nonhomogeneous heat conduction equation is written as a new parabolic type partial differential equation in the space of (x, y, z, t, t). In addition, a fictitious viscous damping coefficient can be used to strengthen the stability of numerical integration of the discretized equations… More >

Chih-Wen Chang^1,2, Chein-Shan Liu³

CMC-Computers, Materials & Continua, Vol.19, No.1, pp. 17-36, 2010, DOI:10.3970/cmc.2010.019.017

Abstract The present study shows a backward group preserving scheme (BGPS) to deal with the multi-dimensional backward wave problem (BWP). The BWP is well-known as seriously ill-posed because the solution does not continuously count on the given data. When three numerical experiments are tested, we reveal that the BGPS is applicable to the multi-dimensional BWP. Even with noisy final data, the BGPS is also robust against perturbation. The numerical results are very pivotal in the computations of multi-dimensional BWP. More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.18, No.1, pp. 1-20, 2010, DOI:10.3970/cmc.2010.018.001

Abstract We are concerned with the reconstruction of an unknown space-dependent rigidity coefficient in a wave equation. This problem is known as one of the inverse scattering problems. Based on a two-point Lie-group equation we develop a Lie-group adaptive method (LGAM) to solve this inverse scattering problem through iterations, which possesses a special character that by using onlytwo boundary conditions and two initial conditions, as those used in the direct problem, we can effectively reconstruct the unknown rigidity function by aself-adaption between the local in time differential governing equation and the global in time algebraic Lie-group equation. The accuracy and efficiency… More >

Weichung Yeih¹, Chein-Shan Liu², Chung-Lun Kuo³, Satya N. Atluri⁴

CMC-Computers, Materials & Continua, Vol.17, No.3, pp. 275-302, 2010, DOI:10.3970/cmc.2010.017.275

Abstract In this paper, a multiple-source-point boundary-collocation Trefftz method, with characteristic lengths being introduced in the basis functions, is proposed to solve the direct, as well as inverse Cauchy problems of the Laplace equation for a multiply connected domain. When a multiply connected domain with genus p (p>1) is considered, the conventional Trefftz method (T-Trefftz method) will fail since it allows only one source point, but the representation of solution using only one source point is impossible. We propose to relax this constraint by allowing many source points in the formulation. To set up a complete set of basis functions, we… More >

Liviu Marin¹

CMC-Computers, Materials & Continua, Vol.17, No.3, pp. 233-274, 2010, DOI:10.3970/cmc.2010.017.233

Abstract We investigate two algorithms involving the relaxation of either the given boundary temperatures (Dirichlet data) or the prescribed normal heat fluxes (Neumann data) on the over-specified boundary in the case of the iterative algorithm of Kozlov91 applied to Cauchy problems for two-dimensional steady-state anisotropic heat conduction (the Laplace-Beltrami equation). The two mixed, well-posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV)… More >

Chih-Wen Chang¹, Chein-Shan Liu²

CMC-Computers, Materials & Continua, Vol.17, No.1, pp. 19-40, 2010, DOI:10.3970/cmc.2010.017.019

Abstract In this study, we employ a semi-analytical approach to solve a two-dimensional advection-dispersion equation (ADE) for identifying the contamination problems. First, the Fourier series expansion technique is used to calculate the concentration field C(x, y, t) at any time t < T. Then, we ponder a direct regularization by adding an extra termaC(x, y, 0) on the final time data C(x, y, T), to reach a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us obtaining a closed-form solution of the Fourier coefficients. A strategy to choose the regularization parameter is offered. The solver utilized in… More >

Hrvoje Gotovac¹, Vedrana Kozulić², Blaž Gotovac¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 173-198, 2010, DOI:10.3970/cmc.2010.015.173

Abstract The space-time Adaptive Fup Collocation Method (AFCM) for solving boundary-initial value problems is presented. To solve the one-dimensional initial boundary value problem, we convert the problem into a two-dimensional boundary value problem. This quasi-boundary value problem is then solved simultaneously in the space-time domain with a collocation technique and by using atomic Fup basis functions. The proposed method is a generally meshless methodology because it requires only the addition of collocation points and basis functions over the domain, instead of the classical domain discretization and numerical integration. The grid is adapted progressively by setting the threshold as a direct measure… More >

Chia-Ming Fan^1,2, Chein-Shan Liu³, Weichung Yeih¹, Hsin-Fang Chan¹

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 67-86, 2010, DOI:10.3970/cmc.2010.015.067

Abstract In this study, the nonlinear obstacle problems, which are also known as the nonlinear free boundary problems, are analyzed by the scalar homotopy method (SHM) and the finite difference method. The one- and two-dimensional nonlinear obstacle problems, formulated as the nonlinear complementarity problems (NCPs), are discretized by the finite difference method and form a system of nonlinear algebraic equations (NAEs) with the aid of Fischer-Burmeister NCP-function. Additionally, the system of NAEs is solved by the SHM, which is globally convergent and can get rid of calculating the inverse of Jacobian matrix. In SHM, by introducing a scalar homotopy function and… More >

Chih-Wen Chang¹, Chein-Shan Liu², Jiang-Ren Chang³

CMC-Computers, Materials & Continua, Vol.15, No.1, pp. 45-66, 2010, DOI:10.3970/cmc.2010.015.045

Abstract In this article, we propose a semi-analytical method to tackle the two-dimensional backward heat conduction problem (BHCP) by using a quasi-boundary idea. First, the Fourier series expansion technique is employed to calculate the temperature field u(x, y, t) at any time t < T. Second, we consider a direct regularization by adding an extra termau(x, y, 0) to reach a second-kind Fredholm integral equation for u(x, y, 0). The termwise separable property of the kernel function permits us to obtain a closed-form regularized solution. Besides, a strategy to choose the regularization parameter is suggested. When several numerical examples were tested,… More >