Weichung Yeih^1,2, Chein-Shan Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.46, No.2, pp. 107-128, 2009, DOI:10.3970/cmes.2009.046.107

Abstract We consider an inverse problem for estimating an unknown time dependent heat source H(t) in a heat conduction equation u_t(x,t) = u_xx(x,t) + H(t). First this inverse problem is formulated as a three-point boundary value problem (BVP) for ODEs discretized from the transformed homogeneous governing equation. To treat this three-point BVP we develop a two-stage Lie-group shooting method (TSLGSM). The novel approach is examined through numerical examples to convince that it is rather accurate and efficient; the estimation error is small even for identifying discontinuous and oscillatory heat sources under noise. More >

Chein-Shan Liu¹, Chih-Wen Chang², Jiang-Ren Chang^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.1, pp. 1-16, 2008, DOI:10.3970/cmes.2008.032.001

Abstract In this paper, we propose a new method to tackle of two famous boundary layer equations in fluid mechanics, namely, the Falkner-Skan and the Blasius equations. We can employ this method to find unknown initial conditions. The pivotal point is based on the erection of a one-step Lie group element$\mathbf {G}(T)$ and the formation of a generalized mid-point Lie group element$\mathbf {G}(r)$. Then, by imposing$\mathbf {G}(T) = \mathbf {G}(r)$ we can seek the missing initial conditions through a minimum discrepancy from the target in terms of a weighting factor$r \in (0, 1)$. Numerical examples are worked out to persuade that… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.30, No.1, pp. 1-16, 2008, DOI:10.3970/cmes.2008.030.001

Abstract In this paper we propose a new numerical integration method of second-order boundary value problems (BVPs) resulting from the elastica of slender rods under different loading conditions and boundary conditions. We construct a compact space shooting method for finding unknown initial conditions. The key point is based on the construction of a one-step Lie group element G(T) and the establishment of a generalized mid-point Lie group element G(r) by using the mean value theorem. Then, by imposing G(T) = G(r) we can search the missing initial condition through a closed-form solution in terms of the weighting factor r ∈ (0,1).… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.27, No.3, pp. 137-150, 2008, DOI:10.3970/cmes.2008.027.137

Abstract For the inverse vibration problem, a Lie-group shooting method is proposed to simultaneously estimate the time-dependent damping and stiffness functions by using two sets of displacement as inputs. First, we transform these two ODEs into two parabolic type PDEs. Second, we formulate the inverse vibration problem as a multi-dimensional two-point boundary value problem with unknown coefficients, allowing us to develop the Lie-group shooting method. For the semi-discretizations of PDEs we thus obtain two coupled sets of linear algebraic equations, from which the estimation of damping and stiffness coefficients can be written out explicitly. The present approach is very interesting, which… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.26, No.3, pp. 157-168, 2008, DOI:10.3970/cmes.2008.026.157

Abstract For the Sturm-Liouville eigenvalues problem we construct a very effective Lie-group shooting method (LGSM) to search the eigenvalues, and when eigenvalue is determined we can also search a missing left-boundary condition of the slope through a weighting factor r ∈ (0,1). Hence, the eigenvalues and eigenfunctions can be calculated with a better accuracy. Because a closed-form formula is derived to calculate unknown slope in terms of λ for the estimation of eigenvalues, the present method is easy to implement and has a low computational cost. Similarly by applying the LGSM to find a corresponding eigenfunction in terms of λ is… 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 space are numerically integrated towards… More >

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

CMC-Computers, Materials & Continua, Vol.7, No.3, pp. 139-154, 2008, DOI:10.3970/cmc.2008.007.139

Abstract In this paper, we propose a Lie-group shooting method to deal with the classical Blasius flat-plate problem and to find unknown initial conditions. The pivotal point is based on the erection of a one-step Lie group element$\mathbf G(T) and the formation of a generalized mid-point Lie group element$\mathbf G(r). Then, by imposing G(T) = G(r) we can derive some algebraic equations to recover the missing initial conditions. It is the first time that we can apply the Lie-group shooting method to solve the classical Blasius flat-plate problem. Numerical examples are worked out to persuade that the novel approach has better… More >

Chein-Shan Liu^1,2, Chung-Lun Kuo³, Dongjie Liu⁴

CMC-Computers, Materials & Continua, Vol.24, No.2, pp. 105-124, 2011, DOI:10.3970/cmc.2011.024.105

Abstract The inverse Cauchy problems for elliptic equations, such as the Laplace equation, the Poisson equation, the Helmholtz equation and the modified Helmholtz equation, defined in annular domains are investigated. The outer boundary of the annulus is imposed by overspecified boundary data, and we seek unknown data on the inner boundary through the numerical solution by a spring-damping regularization method and its Lie-group shooting method (LGSM). Several numerical examples are examined to show that the LGSM can overcome the ill-posed behavior of inverse Cauchy problem against the disturbance from random noise, and the computational cost is very cheap. More >