J. Sladek¹, V. Sladek², P. Solek¹, S.N. Atluri³

CMES-Computer Modeling in Engineering & Sciences, Vol.34, No.3, pp. 273-300, 2008, DOI:10.3970/cmes.2008.034.273

Abstract A meshless method based on the local Petrov-Galerkin approach is proposed, to solve boundary and initial value problems of piezoelectric and magneto-electric-elastic solids with continuously varying material properties. Stationary and transient dynamic 2-D problems are considered in this paper. The mechanical fields are described by the equations of motion with an inertial term. To eliminate the time-dependence in the governing partial differential equations the Laplace-transform technique is applied to the governing equations, which are satisfied in the Laplace-transformed domain in a weak-form on small subdomains. Nodal points are spread on the analyzed domain, and each node is surrounded by a… More >

J. Sladek¹, V. Sladek¹, Ch. Zhang², C.L. Tan³

CMES-Computer Modeling in Engineering & Sciences, Vol.16, No.1, pp. 57-68, 2006, DOI:10.3970/cmes.2006.016.057

Abstract The Meshless Local Petrov-Galerkin (MLPG) method for linear transient coupled thermoelastic analysis is presented. Orthotropic material properties are considered here. A Heaviside step function as the test functions is applied in the weak-form to derive local integral equations for solving two-dimensional (2-D) problems. In transient coupled thermoelasticity an inertial term appears in the equations of motion. The second governing equation derived from the energy balance in coupled thermoelasticity has a diffusive character. To eliminate the time-dependence in these equations, the Laplace-transform technique is applied to both of them. Local integral equations are written on small sub-domains with a circular shape.… More >

Paola Gardano¹, Peter Dabnichki^1,2

CMES-Computer Modeling in Engineering & Sciences, Vol.15, No.2, pp. 87-98, 2006, DOI:10.3970/cmes.2006.015.087

Abstract The work presents a numerical simulation of hydrodynamic forces generated in front crawl swimming. The three dimensional Laplace's equation is used for the analysis of the flow around a moving body in an infinite domain and considers the effect of the added mass and the acceleration on the hydrodynamic forces (Drag and Lift) generated by the interaction between the flow and the body at different geometric configurations of the arm -- variable elbow angle. Boundary Element Method (BEM) was used to obtain the solution of the three dimensional equation numerically. The aim of the work was two-fold:
1) to… More >

J. Sladek¹, V. Sladek¹, P. H. Wen², M.H. Aliabadi³

CMES-Computer Modeling in Engineering & Sciences, Vol.13, No.2, pp. 103-118, 2006, DOI:10.3970/cmes.2006.013.103

Abstract A meshless local Petrov-Galerkin (MLPG) method is applied to solve bending problems of shear deformable shallow shells described by the Reissner theory. Both static and dynamic loads are considered. For transient elastodynamic case the Laplace-transform is used to eliminate the time dependence of the field variables. A weak formulation with a unit test function transforms the set of governing equations into local integral equations on local subdomains in the mean surface of the shell. Nodal points are randomly spread on that surface and each node is surrounded by a circular subdomain to which local integral equations are applied. The meshless… 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 >

Roman Chapko¹, B. Tomas Johansson²

CMES-Computer Modeling in Engineering & Sciences, Vol.85, No.2, pp. 105-128, 2012, DOI:10.3970/cmes.2012.085.105

Abstract We consider a Cauchy problem for the Laplace equation in a 3-dimen -sional semi-infinite domain that contains a bounded inclusion. The canonical situation is the upper half-space in I\tmspace -.1667em R3 containing a bounded smooth domain. The function value of the solution is specified throughout the plane bounding the upper half-space, and the normal derivative is given only on a finite portion of this plane. The aim is to reconstruct the solution on the surface of the bounded inclusion. This is a generalisation of the situation in Chapko and Johansson (2008) to three-dimensions and with Cauchy data only partially given.… 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 >

Liviu Marin¹

CMES-Computer Modeling in Engineering & Sciences, Vol.48, No.2, pp. 121-154, 2009, DOI:10.3970/cmes.2009.048.121

Abstract We investigate the numerical implementation of the alternating iterative algorithm originally proposed by ` 12 ` 12 `$12 `&12 `#12 `^12 `_12 `%12 `~12 *Kozlov91 in the case of the Cauchy problem for the two-dimensional Laplace equation using a meshless method. The two mixed, well-posed and direct problems corresponding to every iteration of the numerical procedure are solved using the method of fundamental solutions (MFS), in conjunction with the Tikhonov regularization method. For each direct problem considered, the optimal value of the regularization parameter is chosen according to the generalized cross-validation (GCV) criterion. An efficient regularizing stopping criterion which ceases… More >

Chih-Chang Chi¹, Weichung Yeih^1,2, Chein-Shan Liu³

CMES-Computer Modeling in Engineering & Sciences, Vol.47, No.2, pp. 167-190, 2009, DOI:10.3970/cmes.2009.047.167

Abstract In this study, a novel method for solving the Cauchy problem of Laplace equation is developed. Through the fictitious time integration method (FTIM), the finding of the root of the resulting linear equations can be transformed into for finding the fixed point of a system of first order ordinary differential equations, in which a fictitious time variable is introduced. In such a sense, the inverse of ill-posed leading matrix is not necessary for the FTIM. This method uses the residual of each equation to control the evolution of unknowns in the fictitious time, and it is different from the conventional… More >