Cheng-Yu Ku^1,2,3,Weichung Yeih^1,2, Chein-Shan Liu⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.76, No.2, pp. 83-108, 2011, DOI:10.3970/cmes.2011.076.083

Abstract In this paper, a general dynamical method based on the construction of a scalar homotopy function to transform a vector function of Non-Linear Algebraic Equations (NAEs) into a time-dependent scalar function by introducing a fictitious time-like variable is proposed. With the introduction of a transformation matrix, the proposed general dynamical method can be transformed into several dynamical Newton-like methods including the Dynamical Newton Method (DNM), the Dynamical Jacobian-Inverse Free Method (DJIFM), and the Manifold-Based Exponentially Convergent Algorithm (MBECA). From the general dynamical method, we can also derive the conventional Newton method using a certain fictitious time-like function. The formulation presented… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.3, pp. 179-196, 2006, DOI:10.3970/cmes.2006.015.179

Abstract This paper studies the numerical computations of the second-order singularly perturbed boundary value problems (SPBVPs). In order to depress the singularity we consider a coordinate transformation from the x-domain to the t-domain. The relation between singularity and stiffness is demonstrated, of which the coordinate transformation parameter λ plays a key role to balance these two tendencies. Then we construct a very effective Lie-group shooting method to search the missing initial condition through a weighting factor r ∈ (0,1) in the t-domain formulation. For stabilizing the new method we also introduce two new systems by a translation of the dependent variable.… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.2, pp. 69-86, 2006, DOI:10.3970/cmes.2006.015.069

Abstract In this paper we will study the numerical integrations of second order boundary value problems under the imposed conditions at t=0 and t=T in a general setting. We can construct a compact space shooting method for finding the unknown initial conditions. The key point is based on the construction of a one-step Lie group element G(u₀,u_T) and the establishment of a mid-point Lie group element G(r). Then, by imposing G(u₀,u_T) = G(r) we can search the missing initial conditions through an iterative solution of the weighting factor r ∈ (0,1). Numerical examples were examined to convince that the new approach… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 149-164, 2006, DOI:10.3970/cmes.2006.013.149

Abstract The present paper provides a Lie-group shooting method for the numerical solutions of second order nonlinear boundary value problems exhibiting multiple solutions. It aims to find all solutions as easy as possible. The boundary conditions considered are classified into four types, namely the Dirichlet, the first Robin, the second Robin and the Neumann. The two Robin type problems are transformed into a canonical one by using the technique of symmetric extension of the governing equations. The Lie-group shooting method is very effective to search unknown initial condition through a weighting factor r ∈ (0,1) Furthermore, the closed-form solutions are derived… More >

Shuncai Li^1,2,3, Zhengzhu Dong², Dan Ma²

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.2, pp. 93-110, 2011, DOI:10.3970/cmes.2011.071.093

Abstract Weak local support is a very common phenomenon in roadway support engineering. It is a problem that needs to be studied thoroughly at the theoretical level. So far, the literature on stress field theory of rock surrounding a roadway is largely restricted to analytical solutions of stress for roadways with a uniform support or no support at all. The corresponding stress solution under conditions of local or weak local support has not been provided. Based on a mechanical model of weak local support at the boundary of a circular roadway and the boundary element method on boundary value problems of… More >

Sunilkumar N¹, D Roy^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.65, No.2, pp. 107-154, 2010, DOI:10.3970/cmes.2010.065.107

Abstract The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries. Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials. There is thus a case for combining these advantages in a so-called hybrid scheme or a 'smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C^p(p ≥ 1) continuity. One such recent attempt, a… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.1, pp. 1-22, 2009, DOI:10.3970/cmes.2009.047.001

Abstract In this paper an analytical method for approximating the solution of backward heat conduction problem is presented. The Fourier series expansion technique is used to formulate a first-kind Fredholm integral equation for the temperature field u(x,t) at any time t < T, when the data are specified at a final time T. Then we consider a direct regularization, instead of the Tikhonov regularization, by adding the term αu(x,t) to obtain a second-kind Fredholm integral equation. The termwise separable property of kernel function allows us by transforming it to a two-point boundary value problem, and thus a closed-form solution is derived.… More >

Cheng-Yu Ku^1,2, Weichung Yeih^1,2, Chein-Shan Liu³, Chih-Chang Chi²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.2, pp. 173-190, 2009, DOI:10.3970/cmes.2009.043.173

Abstract In this paper, a new time-like function with the incorporation of the fictitious time integration method (FTIM) is proposed. The new time-like function is modified from the original time-like function in the FTIM by adding a control parameter m, which dramatically improves the performance of the FTIM for solving highly nonlinear boundary value problems (BVPs) and plays as an important controller to assure the convergence of the solution during the time integration process. The requirements and the characteristics of the new time-like function are presented by examining the FTIM through the perspective of the new time-like function in which the… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.39, No.2, pp. 125-154, 2009, DOI:10.3970/cmes.2009.039.125

Abstract We propose a new numerical method for solving the boundary value problems of ordinary differential equations (ODEs) under multipoint boundary conditions specified at t = T_i, i = 1,...,m, where T₁ < ... < T_m. The finite difference scheme is used to approximate the ODEs, which together with the m-point boundary conditions constitute a system of nonlinear algebraic equations (NAEs). Then a Fictitious Time Integration Method (FTIM) is used to solve these NAEs. Numerical examples confirm that the new approach is highly accurate and efficient with a fast convergence. The FTIM can also be used to find the periods of… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.33, No.2, pp. 179-198, 2008, DOI:10.3970/cmes.2008.033.179

Abstract Dirichlet boundary value problem of quasilinear elliptic equation is numerically solved by using a new concept of fictitious time integration method (FTIM). We introduce a fictitious time coordinate t by transforming the dependent variable u(x,y) into a new one by (1+t)u(x,y) =: v(x,y,t), such that the original equation is naturally and mathematically equivalently written as a quasilinear parabolic equation, including a viscous damping coefficient to enhance stability in the numerical integration of spatially semi-discretized equation as an ordinary differential equations set on grid points. Six examples of Laplace, Poisson, reaction diffusion, Helmholtz, the minimal surface, as well as the explosion… More >