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 >

Nicolas Ali Libre^1,2, Arezoo Emdadi², Edward J. Kansa^3,4, Mohammad Rahimian², Mohammad Shekarchi²

CMES-Computer Modeling in Engineering & Sciences, Vol.24, No.1, pp. 61-80, 2008, DOI:10.3970/cmes.2008.024.061

Abstract The numerical solution of partial differential equations (PDEs) with Neumann boundary conditions (BCs) resulted from strong form collocation scheme are typically much poorer in accuracy compared to those with pure Dirichlet BCs. In this paper, we show numerically that the reason of the reduced accuracy is that Neumann BC requires the approximation of the spatial derivatives at Neumann boundaries which are significantly less accurate than approximation of main function. Therefore, we utilize boundary treatment schemes that based upon increasing the accuracy of spatial derivatives at boundaries. Increased accuracy of the spatial derivative approximation can be achieved by h-refinement reducing the… More >

Chein-Shan Liu ¹

CMES-Computer Modeling in Engineering & Sciences, Vol.20, No.2, pp. 111-122, 2007, DOI:10.3970/cmes.2007.020.111

Abstract A highly accurate new solver is developed to deal with interior and exterior mixed-boundary value problems for two-dimensional Laplace equation, including the singular ones. To promote the present study, we introduce a circular artificial boundary which is uniquely determined by the physical problem domain, and derive a Dirichlet to Robin mapping on that artificial circle, which is an exact boundary condition described by the first kind Fredholm integral equation. As a consequence, we obtain a modified Trefftz method equipped with a characteristic length factor, ensuring that the new solver is stable because the condition number can be greatly reduced. Then,… More >

Tarun Kant¹, Sandeep S. Pendhari², Yogesh M. Desai³

CMES-Computer Modeling in Engineering & Sciences, Vol.17, No.2, pp. 135-162, 2007, DOI:10.3970/cmes.2007.017.135

Abstract A two-point boundary value problem (BVP) is formed in the present work governed by a set of first-order coupled ordinary differential equations (ODEs) in terms of displacements and the transverse stresses through the thickness of laminate (in domain -h/2 < z < h/2) by introducing partial discretization methodology only in the plan area of the three dimensional (3D) laminate. The primary dependent variables in the ODEs are those which occur naturally on a plane z=a constant. An effective numerical integration (NI) technique is utilized for tackling the two-point BVP in an efficient manner. Numerical studies on cross-ply and angle-ply composite… More >

A.Kurenkov¹, A.Thess², H.Babovsky³

FDMP-Fluid Dynamics & Materials Processing, Vol.2, No.1, pp. 47-58, 2006, DOI:10.3970/fdmp.2006.002.047

Abstract A possibility for the determination of the interface between two electrically conducting fluids in cylindrical geometry is presented. The fluids with different conductivities are situated in an infinite cylinder. Along the axis of the cylinder a homogeneous electrical current is applied. The perturbation of the interface leads to an inhomogeneous electrical current and, therefore, results in an electrical potential change in the fluids and a magnetic field modification outside the fluids. The dependence of the electrical potential on the interface shape is obtained analytically. The interface profile is then recovered from data of the electrical potential measurements which have been… More >

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

CMC-Computers, Materials & Continua, Vol.9, No.2, pp. 111-136, 2009, DOI:10.3970/cmc.2009.009.111

Abstract In this paper, we take the advantage of an analytical method to solve the advection-dispersion equation (ADE) for identifying the contamination problems. First, the Fourier series expansion technique is employed to calculate the concentration field C(x, t) at any time t< T. Then, we consider a direct regularization by adding an extra term αC(x,0) on the final condition to carry off a second kind Fredholm integral equation. The termwise separable property of the kernel function permits us to transform itinto a two-point boundary value problem. The uniform convergence and error estimate of the regularized solution C^α(x,t) are provided and a… 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¹

CMC-Computers, Materials & Continua, Vol.11, No.1, pp. 15-32, 2009, DOI:10.3970/cmc.2009.011.015

Abstract Motivated by the evolutionary and dissipative properties of parabolic type partial differential equation (PDE), Liu (2008a) has proposed a natural and mathematically equivalent approach by transforming the quasilinear elliptic PDE into a parabolic one. However, the above paper only considered a rectangular domain in the plane, and did not treat the difficulty arisen from the quasilinear PDE defined in an arbitrary plane domain. In this paper we propose a new technique of internal and boundary residuals in a fictitious rectangular domain, which are driving forces for the ordinary differential equations based on the Fictitious Time Integration Method (FTIM). Several numerical… More >