J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.5, pp. 477-490, 2004, DOI:10.3970/cmes.2004.006.477

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for solution of static and elastodynamic problems in a homogeneous anisotropic medium. The Heaviside step function is used as the test functions in the local weak form. It is leading to derive local boundary integral equations (LBIEs). For transient elastodynamic problems the Laplace transfor technique is applied and the LBIEs are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry such as circles in 2-d problems. The final form of local integral equations has a pure contour character only in… More >

Xiaomin Wang^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.102, No.2, pp. 169-182, 2014, DOI:10.3970/cmes.2014.102.169

Abstract In this paper, an efficient and robust wavelet Laplace inversion method of solving the fractional differential equations is proposed. Such an inverse function can be applied to any reasonable function categories and it is not necessary to know the properties of original function in advance. As an example, we have applied the proposed method to the solution of the Bagley–Torvik equations and Numerical examples are given to demonstrate the efficiency and accuracy of the proposed. More >

P. Jose¹, R.Kanapady², K.K.Tamma³

CMES-Computer Modeling in Engineering & Sciences, Vol.5, No.5, pp. 409-422, 2004, DOI:10.3970/cmes.2004.005.409

Abstract In this article, a novel hybrid finite element and Laplace transform formulation is presented for the computations of transient electromagnetic fields. The formulation is first based on application of Laplace transform technique for the pertinent differential equations, namely the Maxwell's equation in the non-integral form with subsequently, employing the Galerkin finite element formulations on the transformed equations to maintain the modeling versatility of complex geometries and numerical features for computational analysis. In addition, in conjunction with the above, proper scaling of the field quantities is applied to improve the condition of the effective global stiffness matrix. The problem is first… More >

Zhou YH^1,2, Wang XM², Wang JZ^1,2 , Liu XJ²

CMES-Computer Modeling in Engineering & Sciences, Vol.77, No.2, pp. 137-160, 2011, DOI:10.3970/cmes.2011.077.137

Abstract In this paper, we present an efficient wavelet-based algorithm for solving a class of fractional vibration, diffusion and wave equations with strong nonlinearities. For this purpose, we first suggest a wavelet approximation for a function defined on a bounded interval, in which expansion coefficients are just the function samplings at each nodal point. As the fractional differential equations containing strong nonlinear terms and singular integral kernels, we then use Laplace transform to convert them into the second type Voltera integral equations with non-singular kernels. Certain property of the integral kernel and the ability of explicit wavelet approximation to the nonlinear… More >

J. Sladek¹, V. Sladek¹, P.H. Wen², Y.C. Hon³

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.2, pp. 191-218, 2009, DOI:10.3970/cmes.2009.048.191

Abstract The meshless local Petrov-Galerkin (MLPG) method is used to solve the inverse heat conduction problem of predicting the distribution of the heat transfer coefficient on the boundary of 2-D and axisymmetric bodies. Using this method, nodes are randomly distributed over the numerical solution domain, and surrounding each of these nodes, a circular sub-domain is introduced. By choosing a unit step function as the test function, the local integral equations (LIE) on the boundaries of these sub-domains are derived. To eliminate the time variation in the governing equation, the Laplace transform technique is applied. The local integral equations are nonsingular and… More >

Chein-Shan Liu¹, Satya N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 253-276, 2009, DOI:10.3970/cmes.2009.043.253

Abstract Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a₀ + a₁x + ... + a_nxⁿ expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a new polynomial interpolation… More >

J. Sladek¹, V. Sladek¹, P. Solek²

CMES-Computer Modeling in Engineering & Sciences, Vol.43, No.3, pp. 223-252, 2009, DOI:10.3970/cmes.2009.043.223

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for solution of static and elastodynamic problems in 3-D continuously non-homogeneous anisotropic bodies. Functionally graded materials (FGM) are multi-phase materials with the phase volume fractions varying gradually in space, in a pre-determined profile. The Heaviside step function is used as the test functions in the local weak form resulting into the derived local integral equations (LIEs). For transient elastodynamic problems either the Laplace transform or the time difference techniques are applied. Nodal points are randomly distributed in the 3D analyzed domain and each node is surrounded by a spherical… More >

J. Sladek¹, V. Sladek¹, C.L. Tan², S.N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.32, No.3, pp. 161-174, 2008, DOI:10.3970/cmes.2008.032.161

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed for the solution of steady-state and transient heat conduction problems in a continuously non-homogeneous anisotropic medium. The Laplace transform is used to treat the time dependence of the variables for transient problems. The analyzed domain is covered by small subdomains with a simple geometry. A weak formulation for the set of governing equations is transformed into local integral equations on local subdomains by using a unit test function. Nodal points are randomly distributed in the 3D analyzed domain and each node is surrounded by a spherical subdomain to which… More >

A. J. Davies¹, D. Crann¹, S. J. Kane¹, C-H. Lai²

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.2, pp. 79-86, 2007, DOI:10.3970/cmes.2007.018.079

Abstract The solution process for diffusion problems usually involves the time development separately from the space solution. A finite difference algorithm in time requires a sequential time development in which all previous values must be determined prior to the current value. The Stehfest Laplace transform algorithm, however, allows time solutions without the knowledge of prior values. It is of interest to be able to develop a time-domain decomposition suitable for implementation in a parallel environment. One such possibility is to use the Laplace transform to develop coarse-grained solutions which act as the initial values for a set of fine-grained solutions. The… More >

J. Sladek¹, V. Sladek¹, S.N. Atluri²

CMES-Computer Modeling in Engineering & Sciences, Vol.6, No.3, pp. 309-318, 2004, DOI:10.3970/cmes.2004.006.309

Abstract Meshless methods based on the local Petrov-Galerkin approach are proposed for solution of steady and transient heat conduction problem in a continuously nonhomogeneous anisotropic medium. Fundamental solution of the governing partial differential equations and the Heaviside step function are used as the test functions in the local weak form. It is leading to derive local boundary integral equations which are given in the Laplace transform domain. The analyzed domain is covered by small subdomains with a simple geometry. To eliminate the number of unknowns on artificial boundaries of subdomains the modified fundamental solution and/or the parametrix with a convenient cut-off… More >