DERIVATION OF PARTICULAR SOLUTIONS USING CHEBYSHEV POLYNOMIAL BASED FUNCTIONS

2007 ◽  
Vol 04 (01) ◽  
pp. 15-32 ◽  
Author(s):  
C. S. CHEN ◽  
SUNGWOOK LEE ◽  
C.-S. HUANG

In this paper, we propose a simple and direct numerical procedure to obtain particular solutions for various types of differential equations. This procedure employs the power series expansion of a differential operator. Chebyshev polynomials are selected as basis functions for the approximation of the inhomogeneous terms of the given partial differential equation. This numerical scheme provides a highly efficient and accurate approximation for the evaluation of a particular solution for a variety of classes of partial differential equations. To demonstrate the effectiveness of the proposed scheme, we couple the method of fundamental solutions to solve a modified Helmholtz equation with irregular boundary configuration. The solutions were observed to have high accuracy.

2016 ◽  
Vol 20 (2) ◽  
pp. 512-533 ◽  
Author(s):  
Ji Lin ◽  
C. S. Chen ◽  
Chein-Shan Liu

AbstractThis paper describes an application of the recently developed sparse scheme of the method of fundamental solutions (MFS) for the simulation of three-dimensional modified Helmholtz problems. The solution to the given problems is approximated by a two-step strategy which consists of evaluating the particular solution and the homogeneous solution. The homogeneous solution is approximated by the traditional MFS. The original dense system of the MFS formulation is condensed into a sparse system based on the exponential decay of the fundamental solutions. Hence, the homogeneous solution can be efficiently obtained. The method of particular solutions with polyharmonic spline radial basis functions and the localized method of approximate particular solutions in combination with the Gaussian radial basis function are employed to approximate the particular solution. Three numerical examples including a near singular problem are presented to show the simplicity and effectiveness of this approach.


2011 ◽  
Vol 08 (03) ◽  
pp. 545-559 ◽  
Author(s):  
C. S. CHEN ◽  
C. M. FAN ◽  
P. H. WEN

A new version of the method of approximate particular solutions (MAPSs) using radial basis functions (RBFs) has been proposed for solving a general class of elliptic partial differential equations. In the solution process, the Laplacian is kept on the left-hand side as a main differential operator. The other terms are moved to the right-hand side and treated as part of the forcing term. In this way, the close-form particular solution is easy to obtain using various RBFs. The numerical scheme of the new MAPSs is simple to implement and yet very accurate. Three numerical examples are given and the results are compared to Kansa's method and the method of fundamental solutions.


2011 ◽  
Vol 27 (4) ◽  
pp. 493-501 ◽  
Author(s):  
C. C. Tsai ◽  
M. E. Quadir ◽  
H. H. Hwung ◽  
T. W. Hsu

ABSTRACTIn this paper, analytical particular solutions of the augmented polyharmonic spline (APS) associated with Reissner plate model are explicitly derived in order to apply the dual reciprocity method. In the derivations of the particular solutions, a coupled system of three second-ordered partial differential equations (PDEs), which governs problems of Reissner plates, is initially transformed into a single six-ordered PDE by the Hörmander operator decomposition technique. Then the particular solutions of the coupled system can be found by using the particular solution of the six-ordered PDE derived in the first author's previous study. These formulas are further implemented for solving problems of Reissner plates under arbitrary loadings. In the solution procedure, an arbitrary loading measured at some scattered points is first interpolated by the APS and a corresponding particular solution can then be approximated by using the prescribed formulas. After that the complementary homogeneous problem is formally solved by the method of fundamental solutions (MFS). Numerical experiments are carried out to validate these particular solutions.


2020 ◽  
Author(s):  
Bianca Kretz ◽  
Willi Freeden ◽  
Volker Michel

<p>The aspect of poroelasticity is anywhere interesting where a solid material and a fluid come into play and have an effect on each other. This is the case in many applications and we want to focus on geothermics. It is useful to consider this aspect since the replacement of the water in the reservoir below the Earth's surface has an effect on the sorrounding material and vice versa. The underlying physical processes can be described by partial differential equations, called the quasistatic equations of poroelasticity (QEP). From a mathematical point of view, we have a set of three (for two space and one time dimension) partial differential equations with the unknowns u (displacement) and p (pore pressure) depending on the space and the time.</p><p>Our aim is to do a decomposition of the data given for u and p in order that we can see underlying structures in the different decomposition scales that cannot be seen in the whole data.<br>For this process, we need the fundamental solution tensor of the QEP (cf. [1],[5]).<br>That means we assume that we have given data for u and p (they can be obtained for example by a method of fundamental solutions, cf. [1]) and want to investigate a post-processing method to these data. Here we follow the basic approaches for the Laplace-, Helmholtz- and d'Alembert equation (cf. [2],[4],[6]) and the  Cauchy-Navier equation as a tensor-valued ansatz (cf. [3]). That means we want to modify our elements of the fundamental solution tensor in such a way that we smooth the singularity concerning a parameter set τ=(τ<sub>x</sub>,τ<sub>t</sub>). <br>With the help of these modified functions, we construct scaling functions which have to fulfil the properties of an approximate identity.<br>They are convolved with the given data to extract more details of u and p.</p><p><strong>References</strong></p><p>[1] M. Augustin: A method of fundamental solutions in poroelasticity to model the stress field in geothermal reservoirs, PhD Thesis, University of Kaiserslautern, 2015, Birkhäuser, New York, 2015.<br>[2] C. Blick, Multiscale potential methods in geothermal research: decorrelation reflected post-processing and locally based inversion, PhD Thesis, Geomathematics Group, Department of Mathematics, University of Kaiserslautern, 2015.<br>[3] C. Blick, S. Eberle, Multiscale density decorrelation by Cauchy-Navier wavelets, Int. J. Geomath. 10, 2019, article 24.<br>[4] C. Blick, W. Freeden, H. Nutz: Feature extraction of geological signatures by multiscale gravimetry. Int. J. Geomath. 8: 57-83, 2017.<br>[5] A.H.D. Cheng and E. Detournay: On singular integral equations and fundamental solutions of poroelasticity. Int. J. Solid. Struct. 35, 4521-4555, 1998.<br>[6] W. Freeden, C. Blick: Signal decorrelation by means of multiscale methods, World of Mining, 65(5):304-317, 2013.<br><br></p>


Mathematics ◽  
2019 ◽  
Vol 7 (1) ◽  
pp. 40 ◽  
Author(s):  
Shumaila Javeed ◽  
Dumitru Baleanu ◽  
Asif Waheed ◽  
Mansoor Shaukat Khan ◽  
Hira Affan

The analysis of Homotopy Perturbation Method (HPM) for the solution of fractional partial differential equations (FPDEs) is presented. A unified convergence theorem is given. In order to validate the theory, the solution of fractional-order Burger-Poisson (FBP) equation is obtained. Furthermore, this work presents the method to find the solution of FPDEs, while the same partial differential equation (PDE) with ordinary derivative i.e., for α = 1 , is not defined in the given domain. Moreover, HPM is applied to a complicated obstacle boundary value problem (BVP) of fractional order.


2020 ◽  
Vol 2020 ◽  
pp. 1-19
Author(s):  
Yassin Belkourchia ◽  
Lahcen Azrar

The dynamic behavior of structures with piezoelectric patches is governed by partial differential equations with strong singularities. To directly deal with these equations, well adapted numerical procedures are required. In this work, the differential quadrature method (DQM) combined with a regularization procedure for space and implicit scheme for time discretization is used. The DQM is a simple method that can be implemented with few grid points and can give results with a good accuracy. However, the DQM presents some difficulties when applied to partial differential equations involving strong singularities. This is due to the fact that the subsidiaries of the singular functions cannot be straightforwardly discretized by the DQM. A methodological approach based on the regularization procedure is used here to overcome this difficulty and the derivatives of the Dirac-delta function are replaced by regularized smooth functions. Thanks to this regularization, the resulting differential equations can be directly discretized using the DQM. The efficiency and applicability of the proposed approach are demonstrated in the computation of the dynamic behavior of beams for various boundary conditions and excited by impulse and Multiharmonics piezoelectric actuators. The obtained numerical results are well compared to the developed analytical solution.


Sign in / Sign up

Export Citation Format

Share Document