Chein-Shan Liu^1,2, Hong-Hua Dai¹, Satya N. Atluri¹

CMES-Computer Modeling in Engineering & Sciences, Vol.81, No.2, pp. 195-228, 2011, DOI:10.3970/cmes.2011.081.195

Abstract In this continuation of a series of our earlier papers, we define a hyper-surface h(x,t) = 0 in terms of the unknown vector x, and a monotonically increasing function Q(t) of a time-like variable t, to solve a system of nonlinear algebraic equations F(x) = 0. If R is a vector related to ∂h / ∂x, , we consider the evolution equation x^· = λ[αR + βP], where P = F − R(F·R) / ||R||² such that P·R = 0; or x^· = λ[αF + βP^∗], where P^∗ = R − F(F·R) / ||F||² such that P^*·F = 0. From these evolution More >

Sushil Mohan Ratnaker¹, Sanjay Mittal^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.80, No.3&4, pp. 179-200, 2011, DOI:10.3970/cmes.2011.080.179

Abstract Flow past a circular cylinder looses stability at Re ~ 47 via Hopf bifurcation. The eigenmode responsible for the instability leads to the von Kármán vortex shedding. In this work the linear stability of the flow to other modes, near the critical Re, is investigated. In particular, the study explores the possibility of modes other than the Kármán mode having the largest growth rate for Re < Re_cr. To this extent, global linear stability analysis (LSA) of the steady flow past a circular cylinder is carried out for Re = 45 and 48. In addition to the Kármán… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.3, pp. 279-304, 2011, DOI:10.3970/cmes.2011.071.279

Abstract For solving a system of nonlinear algebraic equations (NAEs) of the type: F(x)=0, or F_i(x_j) = 0, i,j = 1,...,n, a Newton-like algorithm has several drawbacks such as local convergence, being sensitive to the initial guess of solution, and the time-penalty involved in finding the inversion of the Jacobian matrix ∂F_i/∂x_j. Based-on an invariant manifold defined in the space of (x,t) in terms of the residual-norm of the vector F(x), we can derive a gradient-flow system of nonlinear ordinary differential equations (ODEs) governing the evolution of x with a fictitious time-like variable t as an independent variable. … More >

Ying-Hsiu Shen¹, Chein-Shan Liu^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.2, pp. 157-178, 2011, DOI:10.3970/cmes.2011.071.157

Abstract When the local differential quadrature (LDQ) has been successfully applied to solve two-dimensional problems, the global method of DQ still has a problem by requiring to solve the inversions of ill-posed matrices. Previously, when one uses (n-1)th order polynomial test functions to determine the weighting coefficients with n grid points, the resultant n ×n Vandermonde matrix is highly ill-conditioned and its inversion is hard to solve. Now we use (m-1)th order polynomial test functions by n grid points that the size of Vandermonde matrix is m×n, of which m is much less than n. We More >

A. Idesman ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.71, No.2, pp. 111-148, 2011, DOI:10.3970/cmes.2011.071.111

Abstract The paper deals with the following issues of existing time-integration methods for a semi-discrete system of elastodynamics equations: a) the quantification and the suppression of spurious high frequencies; b) the selection of the amount of numerical dissipation for a time-integration method; and c) accurate time integration of low modes. The finite element method used in the paper or other methods can be applied for the space discretization. A new two-stage time-integration procedure consisting of basic computations and the filtering stage is developed. For accurate integration of all frequencies, a time-integration method with zero (or small)… More >

Chih-Wen Chang¹

CMC-Computers, Materials & Continua, Vol.21, No.2, pp. 87-106, 2011, DOI:10.3970/cmc.2011.021.087

Abstract We address a new numerical approach to deal with these multi-dimensional backward wave problems (BWPs) in this study. A fictitious time τ is utilized to transform the dependent variable u(x, y, z, t) into a new one by (1+τ)u(x, y, z, t)=: v(x, y, z, t, τ), such that the original wave equation is written as a new hyperbolic type partial differential equation in the space of (x, y, z, t, τ). Besides, a fictitious viscous damping coefficient can be employed to strengthen the stability of numerical integration of the discretized equations by using a group preserving scheme. More >

Chein-Shan Liu¹, Weichung Yeih², Satya N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.3, pp. 301-324, 2010, DOI:10.3970/cmes.2010.059.301

Abstract When problems in engineering and science are discretized, algebraic equations appear naturally. In a recent paper by Liu and Atluri, non-linear algebraic equations (NAEs) were transformed into a nonlinear system of ODEs, which were then integrated by a method labelled as the Fictitious Time Integration Method (FTIM). In this paper, the FTIM is enhanced, by using the concept of arepellorin the theory ofnonlinear dynamical systems, to situations where multiple-solutions exist. We label this enhanced method as MSFTIM. MSFTIM is applied and illustrated in this paper through solving boundary-layer problems, boundary-value problems, and eigenvalue problems with More >

Chia-Cheng Tsai¹, Chein-Shan Liu², Wei-Chung Yeih³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 131-152, 2010, DOI:10.3970/cmes.2010.056.131

Abstract The fictitious time integration method (FTIM) previously developed by Liu and Atluri (2008a) is combined with the method of fundamental solutions and the Chebyshev polynomials to solve Poisson-type nonlinear PDEs. The method of fundamental solutions with Chebyshev polynomials (MFS-CP) is an exponentially-convergent meshless numerical method which is able to solving nonhomogeneous partial differential equations if the fundamental solution and the analytical particular solutions of the considered operator are known. In this study, the MFS-CP is extended to solve Poisson-type nonlinear PDEs by using the FTIM. In the solution procedure, the FTIM is introduced to convert 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 More >

Chein-Shan Liu¹

CMC-Computers, Materials & Continua, Vol.15, No.3, pp. 221-240, 2010, DOI:10.3970/cmc.2010.015.221

Abstract We propose a simple numerical scheme for solving the space- and time-fractional derivative Burgers equations: D_t^αu + εuu_x = vu_xx + ηD_x^βu, 0 < α, β ≤ 1, and u_t + D_*^β(D_*^1-βu)²/2 = vu_xx, 0 < β ≤ 1. The time-fractional derivative D_t^αu and space-fractional derivative D_x^βu are defined in the Caputo sense, while D_*^βu is the Riemann-Liouville space-fractional derivative. A fictitious time τ is used to transform the dependent variable u(x,t) into a new one by (1+τ)^γu(x,t) =: v(x,t,τ), where 0 < γ ≤ 1 is a parameter, such that the original equation is written as a new functional-differential type partial differential equation More >