S. D. Akbarov^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.92, No.4, pp. 387-421, 2013, DOI:10.3970/cmes.2013.092.387

Abstract Flexural wave dispersion in finitely pre-stretched (or pre-compressed) solid and hollow, circular cylinders is investigated with the use of the threedimensional linearized theory of elastic waves in initially stressed bodies. It is assumed that the initial strains in the cylinders are homogeneous and correspond to the uniaxial tension, or compression, along their central axes. The elasticity relations of the cylinders’ materials are described by the harmonic potential. The analytical solution of the corresponding field equations is presented and, using these solutions, the dispersion equations for the cases under consideration are obtained. The dispersion equations are solved numerically and based on… More >

Chih-Wen Chang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.88, No.2, pp. 65-92, 2012, DOI:10.3970/cmes.2012.088.065

Abstract We utilize a regularized integral equation scheme to resolve the three-dimensional backward advection-dispersion equation (BADE) for identifying the groundwater pollution source identification problems in this research. First, the Fourier series expansion method is employed to estimate the concentration field C(x, y, z, t) at any time t < T. Second, we contemplate a direct regularization by adding an extra term a(x, y, z, 0) to transform a second-kind Fredholm integral equation for C(x, y, z, 0). The termwise separable property of the kernel function permits us to acquire a closed-form regularized solution. In addition, a tactic to determine the regularization… More >

Chih-Wen Chang¹, Chein-Shan Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.3, pp. 261-276, 2009, DOI:10.3970/cmes.2009.051.261

Abstract The backward advection-dispersion equation (ADE) for identifying the groundwater pollution source identification problems (GPSIPs) is numerically solved by employing a fictitious time integration method (FTIM). The backward ADE is renowned as ill-posed because the solution does not continuously count on the data. We transform the original parabolic equation into another parabolic type evolution equation by introducing a fictitious time coordinate, and adding a viscous damping coefficient to enhance the stability of numerical integration of the discretized equations by employing a group preserving scheme. When several numerical examples are amenable, we find that the FTIM is applicable to retrieve all past… More >

C.-S. Liu¹, C.-W. Chang², J.-R. Chang²

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 67-82, 2006, DOI:10.3970/cmes.2006.012.067

Abstract In this paper we are concerned with the backward problems governed by differential equations. It is a first time that we can construct a backward time dynamics on the past cone, such that an augmented dynamical system of the Lie type X^˙ = B(X,t)X with t ∈ R⁻, X ∈ Mⁿ⁺¹ lying on the past cone and B ∈ so(n,1), was derived for the backward differential equations system x^· =f(x,t), t ∈ R⁻, x ∈ Rⁿ. These two differential equations systems are mathematically equivalent. Then we apply the backward group preserving scheme (BGPS), which is an explicit single-step… More >

Chein-Shan Liu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.12, No.1, pp. 55-66, 2006, DOI:10.3970/cmes.2006.012.055

Abstract In this paper we are concerned with the numerical integration of Burgers equation backward in time. We construct a one-step backward group preserving scheme (BGPS) for the semi-discretization of Burgers equation. The one-step BGPS is very effectively to calculate the solution at an initial time t = 0 from a given final data at t = T, which with a time stepsize equal to T and with a suitable grid length produces a highly accurate solution never seen before. Under noisy final data the BGPS is also robust to against the disturbance. When the solution appears steep gradient, several steps… More >

C.H. Tsai¹, D.L. Young², J. Kolibal³

CMES-Computer Modeling in Engineering & Sciences, Vol.66, No.1, pp. 53-72, 2010, DOI:10.3970/cmes.2010.066.053

Abstract The time evolution method of fundamental solutions (MFS) is proposed to solve backward heat conduction problems (BHCPs). The time evolution MFS belongs to one of the mesh-free numerical methods and is essentially composed of a sequence of diffusion fundamental solutions which exactly satisfy the heat conduction equations. Through correct treatment of temporal evolution, the resulting system of the time evolution MFS is smaller, and effectively decreases the possibility of ill-conditioning induced by such strongly ill-posed problems. Both one-dimensional and two-dimensional BHCPs are examined in this study, and the numerical results demonstrate the accuracy and stability of the MFS, especially for… More >

Chih-Wen Chang¹, Chein-Shan Liu²

CMES-Computer Modeling in Engineering & Sciences, Vol.59, No.3, pp. 239-274, 2010, DOI:10.3970/cmes.2010.059.239

Abstract In this article, we propose a backward group preserving scheme (BGPS) to tackle the multi-dimensional backward heat conduction problem (BHCP). The BHCP is well-known as severely ill-posed because the solution does not continuously depend on the given data. When eight numerical examples (including nonlinear and nonhomogeneous BHCP, and Neumann and Robin conditions of homogeneous BHCP) are examined, we find that the BGPS is applicable to the multi-dimensional BHCP. Even with noisy final data, the BGPS is also robust against disturbance. The one-step BGPS effectively reconstructs the initial data from the given final data, which with a suitable grid length produces… More >

J. Sladek¹, V. Sladek¹, M. Wünsche², Ch. Zhang²

CMES-Computer Modeling in Engineering & Sciences, Vol.54, No.2, pp. 223-252, 2009, DOI:10.3970/cmes.2009.054.223

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed, to solve the interface crack problem between two dissimilar anisotropic elastic solids. Both stationary and transient mechanical and thermal loads are considered for two-dimensional (2-D) problems in this paper. A Heaviside step function as the test functions is applied in the weak-form to derive local integral equations. Nodal points are spread on the analyzed domain, and each node is surrounded by a small circle for simplicity. The spatial variations of the displacements and temperature are approximated by the Moving Least-Squares (MLS) scheme. After performing the spatial integrations, one obtains… 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 >

R. Vertnik¹, B. Šarler²

CMES-Computer Modeling in Engineering & Sciences, Vol.44, No.1, pp. 65-96, 2009, DOI:10.3970/cmes.2009.044.065

Abstract The application of the mesh-free Local Radial Basis Function Collocation Method (LRBFCM) in solution of incompressible turbulent flow is explored in this paper. The turbulent flow equations are described by the low - Re number k-emodel with Jones and Launder [Jones and Launder (1971)] closure coefficients. The involved velocity, pressure, turbulent kinetic energy and dissipation fields are represented on overlapping 5-noded sub-domains through collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second derivatives of the fields are calculated from the respective derivatives of the RBF's. The velocity, turbulent kinetic energy and dissipation equations are solved through… More >