| (5 intermediate revisions by the same user not shown) | |||
| Line 43: | Line 43: | ||
The main issue in MFS is the source positioning in order obtain reasonable responses. The source points can be placed outside the problem domain, avoiding singularities over the space coordinates, and the solution can be approximated over time by Laplace transform [8,9,10] or by finite difference scheme [11], an interesting approach since MFS is well-posed problem in spatial domain. Otherwise, the source points can be placed at different time levels [7], and, by using transient fundamental solutions, they will not generate singularities in time and spatial domains, however generating an ill-posed problem. | The main issue in MFS is the source positioning in order obtain reasonable responses. The source points can be placed outside the problem domain, avoiding singularities over the space coordinates, and the solution can be approximated over time by Laplace transform [8,9,10] or by finite difference scheme [11], an interesting approach since MFS is well-posed problem in spatial domain. Otherwise, the source points can be placed at different time levels [7], and, by using transient fundamental solutions, they will not generate singularities in time and spatial domains, however generating an ill-posed problem. | ||
| − | For homogeneous diffusion equations, in Golberg et al. [12] for non-zero initial conditions, rectangular shaped domain and Dirichlet’s boundary conditions prescribed, the Laplace transform method and the finite difference method were studied. For both cases it leads to a non-homogeneous Helmholtz equation. Particular solutions, for different types of problems, are generally found by the dual reciprocity method (DRM), introduced by [13,14], that transform the domain integral into a boundary-type linear combination of radial basis functions (RBFs), and Golberg et al. [15] proposed particular solutions for Helmholtz-type operators. For a similar geometry, in | + | For homogeneous diffusion equations, in Golberg et al. [12] for non-zero initial conditions, rectangular shaped domain and Dirichlet’s boundary conditions prescribed, the Laplace transform method and the finite difference method were studied. For both cases it leads to a non-homogeneous Helmholtz equation. Particular solutions, for different types of problems, are generally found by the dual reciprocity method (DRM), introduced by [13,14], that transform the domain integral into a boundary-type linear combination of radial basis functions (RBFs), and Golberg et al. [15] proposed particular solutions for Helmholtz-type operators. For a similar geometry, in Young et al. [7] transient fundamental solutions were applied for homogeneous Dirichlet’s boundary conditions with more complex initial conditions. |
In Cao et al. [16], DRM-MFS based on the modified Helmholtz operator and DRM-Trefftz were applied for homogeneous diffusion equation with mixed boundary conditions and constant initial conditions, Chantasiriwan [17] compared MFS based on transient fundamental solutions with MFS based on the modified Helmholtz operator for a square with different types of boundary conditions. | In Cao et al. [16], DRM-MFS based on the modified Helmholtz operator and DRM-Trefftz were applied for homogeneous diffusion equation with mixed boundary conditions and constant initial conditions, Chantasiriwan [17] compared MFS based on transient fundamental solutions with MFS based on the modified Helmholtz operator for a square with different types of boundary conditions. | ||
| Line 270: | Line 270: | ||
| <math>\psi \left( r\right) =\, \sum _{j=1}^{N}{\alpha }_{j}\left[ -\frac{4}{{\lambda }^{4}}-\frac{4\, ln{r}_{ij}}{{\lambda }^{4}}-\frac{{r}_{ij}^{2}\, ln{r}_{ij}}{{\lambda }^{2}}-\frac{4{K}_{0}\left( \lambda {r}_{ij}\right) }{{\lambda }^{2}}\right] -</math><math>{\alpha }_{n+1}\frac{1}{{\lambda }^{2}}-{\alpha }_{n+2}\frac{x}{{\lambda }^{2}}-</math><math>{\alpha }_{n+3}\frac{y}{{\lambda }^{2}}\, ,\quad </math> for 2D and <math display="inline">r\not =0</math> | | <math>\psi \left( r\right) =\, \sum _{j=1}^{N}{\alpha }_{j}\left[ -\frac{4}{{\lambda }^{4}}-\frac{4\, ln{r}_{ij}}{{\lambda }^{4}}-\frac{{r}_{ij}^{2}\, ln{r}_{ij}}{{\lambda }^{2}}-\frac{4{K}_{0}\left( \lambda {r}_{ij}\right) }{{\lambda }^{2}}\right] -</math><math>{\alpha }_{n+1}\frac{1}{{\lambda }^{2}}-{\alpha }_{n+2}\frac{x}{{\lambda }^{2}}-</math><math>{\alpha }_{n+3}\frac{y}{{\lambda }^{2}}\, ,\quad </math> for 2D and <math display="inline">r\not =0</math> | ||
|} | |} | ||
| − | | style="text-align: right;vertical-align: center;"| ( | + | | style="text-align: right;vertical-align: center;"| (14a) |
|- | |- | ||
| | | | ||
| Line 277: | Line 277: | ||
|<math>\psi (r)=\, \sum _{j=1}^{N}{\alpha }_{j}\left[ -\frac{4}{{\lambda }^{4}}-\frac{4\, \gamma }{{\lambda }^{4}}+\, \frac{4\, }{{\lambda }^{4}}\mathrm{ln}\,(\frac{\lambda }{2})\right] -</math><math>{\alpha }_{n+1}\frac{1}{{\lambda }^{2}}-{\alpha }_{n+2}\frac{x}{{\lambda }^{2}}-</math><math>{\alpha }_{n+3}\frac{y}{{\lambda }^{2}},\quad </math> for 2D and <math display="inline">r=</math><math>0</math> | |<math>\psi (r)=\, \sum _{j=1}^{N}{\alpha }_{j}\left[ -\frac{4}{{\lambda }^{4}}-\frac{4\, \gamma }{{\lambda }^{4}}+\, \frac{4\, }{{\lambda }^{4}}\mathrm{ln}\,(\frac{\lambda }{2})\right] -</math><math>{\alpha }_{n+1}\frac{1}{{\lambda }^{2}}-{\alpha }_{n+2}\frac{x}{{\lambda }^{2}}-</math><math>{\alpha }_{n+3}\frac{y}{{\lambda }^{2}},\quad </math> for 2D and <math display="inline">r=</math><math>0</math> | ||
|} | |} | ||
| − | | style="text-align: right;vertical-align: center;"| ( | + | | style="text-align: right;vertical-align: center;"| (14b) |
|} | |} | ||
| Line 302: | Line 302: | ||
where <math display="inline">\Delta t={t}^{n+1}-{t}^{n}</math> is time step size and <math display="inline">\theta =</math><math>1</math> defining an implicit time marching method. | where <math display="inline">\Delta t={t}^{n+1}-{t}^{n}</math> is time step size and <math display="inline">\theta =</math><math>1</math> defining an implicit time marching method. | ||
| − | Substituting and rearranging terms of Eq. (16) in Eq. (1) and in conditions | + | Substituting and rearranging terms of Eq. (16) in Eq. (1) and in conditions Eqs. (3)-(4) it yields to the modified Helmholtz equation: |
{| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" | {| class="formulaSCP" style="width: 100%;border-collapse: collapse;width: 100%;text-align: center;" | ||
| Line 338: | Line 338: | ||
|} | |} | ||
| − | where <math display="inline">{\lambda }^{2}=\sqrt{\frac{1}{k\Delta t}}</math> and <math display="inline">A\left( \mathit{\boldsymbol{X}}\right) =</math><math>-\frac{{U}^{n}\left( \mathit{\boldsymbol{X}}\right) }{k\Delta t}-\frac{F(\mathit{\boldsymbol{X}})}{k}</math> . It is worth noticing that the term <math display="inline">A\left( \mathit{\boldsymbol{X}}\right)</math> is well defined and depends on the solution obtained from the previous time and the nonhomogenous condition, which is time-independent and is calculated by DRM proposed in | + | where <math display="inline">{\lambda }^{2}=\sqrt{\frac{1}{k\Delta t}}</math> and <math display="inline">A\left( \mathit{\boldsymbol{X}}\right) =</math><math>-\frac{{U}^{n}\left( \mathit{\boldsymbol{X}}\right) }{k\Delta t}-\frac{F(\mathit{\boldsymbol{X}})}{k}</math> . It is worth noticing that the term <math display="inline">A\left( \mathit{\boldsymbol{X}}\right)</math> is well defined and depends on the solution obtained from the previous time and the nonhomogenous condition, which is time-independent and is calculated by DRM proposed in Section 3.1. The initial condition should be applied in order to start the procedure. |
By using MFS for every time step, the homogenous solution is approximated by: | By using MFS for every time step, the homogenous solution is approximated by: | ||
| Line 419: | Line 419: | ||
{| style="text-align: center; margin:auto;width: 100%;" | {| style="text-align: center; margin:auto;width: 100%;" | ||
|- | |- | ||
| − | | <math>{U}^{\ast }\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }};{t}_{F},t\right) =</math><math>\, \frac{1}{{\left( 4\pi k\tau \right) }^{\frac{d}{2}}}exp\left[ \frac{-{r}^{2}\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{4k\tau }\right]</math> | + | | <math>{U}^{\ast }\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }};{t}_{F},t\right) =</math><math>\, \frac{1}{{\left( 4\pi k\tau \right) }^{\frac{d}{2}}}\exp\left[ \frac{-{r}^{2}\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{4k\tau }\right]</math> |
|- | |- | ||
| − | | <math>{Q}^{\ast }\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }};{t}_{F},t\right) =</math><math>\, \frac{-r\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) \frac{\partial r\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{\partial n}}{8\pi k{\tau }^{2}}exp\left[ \frac{-{r}^{2}\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{4k\tau }\right] \,</math> | + | | <math>{Q}^{\ast }\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }};{t}_{F},t\right) =</math><math>\, \frac{-r\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) \frac{\partial r\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{\partial n}}{8\pi k{\tau }^{2}}\exp\left[ \frac{-{r}^{2}\left( \mathit{\boldsymbol{X}},\mathit{\boldsymbol{\xi }}\right) }{4k\tau }\right] \,</math> |
|} | |} | ||
| style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(23) | | style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(23) | ||
| Line 608: | Line 608: | ||
|} | |} | ||
| − | In this example, for MSF1 and MSF2, the coefficient <math display="inline">b=0.5</math>, and, for the modified Helmholtz operator, virtual sources are 0.5 from the physical boundary. In order to explore the methodologies, two different types of RBF were considered, referenced as RBF1 the TPS described on | + | In this example, for MSF1 and MSF2, the coefficient <math display="inline">b=0.5</math>, and, for the modified Helmholtz operator, virtual sources are 0.5 from the physical boundary. In order to explore the methodologies, two different types of RBF were considered, referenced as RBF1 the TPS described on Section 3.1 and RBF2 the following functions [23]: |
{| class="formulaSCP" style="width: 100%; text-align: center;" | {| class="formulaSCP" style="width: 100%; text-align: center;" | ||
| Line 692: | Line 692: | ||
| − | [[#img-12|Figure 12]] depicts absolute errors at | + | [[#img-12|Figure 12]] depicts absolute errors at nodes A and B, for different time steps in all methods. When compared, MSF1 and MSF based on the modified Helmholtz operator present better responses as time step decreases, and MSF2, a method that is known for its ill posed matrix, presents better results for a bigger time step. All methods showed to be sensitive concerning the field node geometrical position. |
<div id='img-12'></div> | <div id='img-12'></div> | ||
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | {| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;" | ||
| Line 738: | Line 738: | ||
==References== | ==References== | ||
| + | <div class="auto" style="text-align: left;width: auto; margin-left: auto; margin-right: auto;font-size: 85%;"> | ||
| − | [1] Gingold | + | [1] Gingold R.A., Monaghan J.J. Smoothed particle hydrodynamics: theory and application to non-spherical stars. Monthly notices of the royal astronomical society, 181(3):375-389, 1977. |
| − | [2] Oñate | + | [2] Oñate E., Idelsohn S., Zienkiewicz O.C., Taylor R.L. A finite point method in computational mechanics. Applications to convective transport and fluid flow. International Journal for Numerical Methods in Engineering, 39(22):3839-3866, 1996. |
| − | [3] Belytschko | + | [3] Belytschko T., Lu Y.Y., Gu L. Element‐free Galerkin methods. International Journal for Numerical Methods in Engineering, 37(2):229-256, 1994. |
| − | [4] Atluri | + | [4] Atluri S.N., Zhu T. A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics. Computational Mechanics, 22(2):117-127, 1998. |
| − | [5] Chen | + | [5] Chen W. Meshfree boundary particle method applied to Helmholtz problems. Engineering Analysis with Boundary Elements, 26(7):577-581, 2002. |
| − | [6] Wang | + | [6] Wang H., Qin Q.H. Methods of fundamental solutions in solid mechanics. Elsevier, 2019. |
| − | [7] Young | + | [7] Young D.L., Tsai C.C., Murugesan K., Fan C.M., Chen C.W. Time-dependent fundamental solutions for homogeneous diffusion problems. Engineering Analysis with Boundary Elements, 28(12):1463-1473, 2004. |
| − | [8] Zhu | + | [8] Zhu S.P. Solving transient diffusion problems: time-dependent fundamental solution approaches versus LTDRM approaches. Engineering Analysis with Boundary Elements, 21(1):87-90, 1998. |
| − | [9] Zhu | + | [9] Zhu S.P., Liu H.W., Lu X.P. A combination of LTDRM and ATPS in solving diffusion problems. Engineering Analysis with Boundary Elements, 21(3):285-289, 1998. |
| − | [10] Sutradhar | + | [10] Sutradhar A., Paulino G.H., Gray L.J. Transient heat conduction in homogeneous and non-homogeneous materials by the Laplace transform Galerkin boundary element method. Engineering Analysis with Boundary Elements, 26(2):119-132, 2002. |
| − | [11] Bulgakov | + | [11] Bulgakov V., Šarler B., Kuhn G. Iterative solution of systems of equations in the dual reciprocity boundary element method for the diffusion equation. International Journal for Numerical Methods in Engineering, 43(4):713-732, 1998. |
| − | [12] Golberg | + | [12] Golberg M.A., Chen C.S., Muleshkov A.S. The method of fundamental solutions for time-dependent problems. WIT Transactions on Modelling and Simulation, 23, 1970. |
| − | [13] Nardini | + | [13] Nardini D., Brebbia C.A. A new approach to free vibration analysis using boundary elements. Applied Mathematical Modelling, 7(3):157-162, 1983. |
| − | [14] Partridge | + | [14] Partridge P.W., Brebbia C.A. Dual reciprocity boundary element method. Springer Science & Business Media, 2012. |
| − | [15] Golberg | + | [15] Golberg M.A., Chen C.S., Rashed Y.F. The annihilator method for computing particular solutions to partial differential equations. Engineering Analysis with Boundary Elements, 23(3):275-279, 1999. |
| − | [16] Cao | + | [16] Cao L., Qin Q.H., Zhao N. Application of DRM-Trefftz and DRM-MFS to transient heat conduction analysis. Recent Patents on Space Technology, 2:41-50, 2010. |
| − | [17] Chantasiriwan | + | [17] Chantasiriwan S. Methods of fundamental solutions for time‐dependent heat conduction problems. International Journal for Numerical Methods in Engineering, 66(1):147-165, 2006. |
| − | [18] Young | + | [18] Young D.L., Tsai C.C., Fan C.M. Direct approach to solve nonhomogeneous diffusion problems using fundamental solutions and dual reciprocity methods. Journal of the Chinese Institute of Engineers, 27(4):597-609, 2004. |
| − | [19] Young | + | [19] Young D.L., Chen C.W., Fan C.M., Tsai C.C. The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation. Numerical Methods for Partial Differential Equations: An International Journal, 22(5):1173-1196, 2006. |
| − | [20] Golberg | + | [20] Golberg M.A., Chen C.S. The theory of radial basis functions applied to the BEM for inhomogeneous partial differential equations. Boundary Elements Communications, 5(2):57-61, 1994. |
| − | [21] Chen | + | [21] Chen C.S., Rashed Y.F. Evaluation of thin plate spline based particular solutions for Helmholtz-type operators for the DRM. Mechanics Research Communications, 25(2):195-201, 1998. |
| − | [22] Golberg | + | [22] Golberg M.A., Chen C.S. The method of fundamental solutions for potential, Helmholtz and diffusion problems. Boundary Integral Methods: Numerical and Mathematical Aspects, 1:103-176, 1998. |
| − | [23] Balakrishnan | + | [23] Balakrishnan K., Sureshkumar R., Ramachandran P. A. An operator splitting-radial basis function method for the solution of transient nonlinear Poisson problems. Computers & Mathematics with Applications, 43(3-5):289-304, 2002. |
| + | </div> | ||
Latest revision as of 14:36, 9 July 2020
Abstract
This paper presents tests and applications related to the Method of the Fundamental Solutions (MFS) for nonhomogeneous diffusion equations. The dual reciprocity method was applied to approximate the particular solution and different MSF were applied to solve the homogeneous equation. The computational accuracy for three different versions submitted to diverse geometries, initial, boundary and nonhomogeneous conditions is compared. The first version is the MFS, based on the fundamental solutions to the modified Helmholtz operator; the second and third versions are based on transient fundamental solutions proposed directly for the diffusion governing equation. This paper evaluates particularities for all methods proposed as the solution advances in time until steady-state solutions are reached. The developed computer codes are applied to several examples with the intention of demonstrating the methods’ effectiveness and applicability.
Keywords: Method of fundamental solutions, meshless methods, diffusion equations, transient diffusion
1. Introduction
The transient conduction problem, governed by the diffusion equations, appears in many areas of science and technology and simulates phenomena such as heat transfer, pollutants concentration, fluid flow, molecular movements, among others. These problems have been extensively treated in literature but, although analytical solutions can be found for some cases, when complex geometries and/or boundary conditions are prescribed, it is necessary to use numerical methods.
Generating a mesh for complicated geometries is time-consuming, requires great memory storage and takes a significant computational effort for classical methods, especially those based on a domain discretization. To overcome such problems, in the last decades there has been a significant increase in the use of meshless numerical methods for solving partial differential equations.
As for the classical methods, the meshless methods can be divided in two main categories: first, the domain-type methods such as the smoothed particle dynamics method (SPH) [1], the finite point method (FPM) [2], the element free Galerkin method (EFG) [3], the meshless local Petrov-Galerkin method (MLPG) [4], etc.; second, the boundary-type methods, as the boundary particle method (BPM) [5] and the method of fundamental solutions (MFS) [6], among others.
Due to the simplified formulation and implementation process, MSF has been largely used for different types of problems. The general concept of MFS is to approximate the solution by a linear combination of the fundamental solutions generated by source points that satisfies the governing equation in the interior domain, leaving only the boundary conditions to be satisfied. Unlike the boundary element method (BEM) in MFS, the solution in the interior domain does not require any additional integral equation. If a fundamental solution exists for the given governing equation, then the method of fundamental solutions can be utilized to obtain numerical solutions of the variable at a definite number of points in the domain [7].
The main issue in MFS is the source positioning in order obtain reasonable responses. The source points can be placed outside the problem domain, avoiding singularities over the space coordinates, and the solution can be approximated over time by Laplace transform [8,9,10] or by finite difference scheme [11], an interesting approach since MFS is well-posed problem in spatial domain. Otherwise, the source points can be placed at different time levels [7], and, by using transient fundamental solutions, they will not generate singularities in time and spatial domains, however generating an ill-posed problem.
For homogeneous diffusion equations, in Golberg et al. [12] for non-zero initial conditions, rectangular shaped domain and Dirichlet’s boundary conditions prescribed, the Laplace transform method and the finite difference method were studied. For both cases it leads to a non-homogeneous Helmholtz equation. Particular solutions, for different types of problems, are generally found by the dual reciprocity method (DRM), introduced by [13,14], that transform the domain integral into a boundary-type linear combination of radial basis functions (RBFs), and Golberg et al. [15] proposed particular solutions for Helmholtz-type operators. For a similar geometry, in Young et al. [7] transient fundamental solutions were applied for homogeneous Dirichlet’s boundary conditions with more complex initial conditions.
In Cao et al. [16], DRM-MFS based on the modified Helmholtz operator and DRM-Trefftz were applied for homogeneous diffusion equation with mixed boundary conditions and constant initial conditions, Chantasiriwan [17] compared MFS based on transient fundamental solutions with MFS based on the modified Helmholtz operator for a square with different types of boundary conditions.
To solve a linear diffusion equation with time-independent source terms, Young et al. [18] applied time-dependent fundamental solutions with DRM for a rectangular domain with different internal sources equations and in Wang and Qin [6] it was used MFS with the method of particular solutions (MPS) and with the eigenfunction expansion method (EEM) to solve several geometries.
In this article, MFS-DRM based on the modified Helmholtz operator with different types of RBFs is applied for different types of geometries, as well as for a set of different internal source terms and boundary conditions. Numerical solutions based on MFS-DRM with time-dependent fundamental solutions, that is, not using any integration techniques or requiring any kind of special geometric position for the source points, are obtained for comparison purposes. Nevertheless, in Young et al. [18] it was proposed a time moving collocation of source points, and in this paper a different approach for the sources is applied. All the described methods will be compared with each other and with the analytical solution in order to evaluate their accuracy and applicability.
This paper is organized as follows: governing equations and boundary conditions are described in Section 2; different formulations for the Method of Fundamental Solutions applied to the diffusion problem is in Section 3; numerical examples and their evaluations are in Section 4; and Section 5 consists of conclusions drawn based on the present study.
2. Governing equations
The transient diffusion equation in a domain Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Omega \left( X\right)}
and a boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Gamma \left( X\right) =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\Gamma }_{u}\left( X\right) +{\Gamma }_{q}\left( X\right)
is given by:
|
(1) |
in which Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathit{\boldsymbol{X\, }}} represents the spatial coordinates, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, t}
is time, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k}
the thermal diffusivity, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle U}
is the potential under analysis, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle F}
is a spatial internal source function and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\nabla }^{2}}
is the Laplacian operator of partial derivatives. The initial condition of the problem is represented by:
|
(2) |
Two boundary conditions can be prescribed. Dirichlet or necessary condition for the potential Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle U}
|
(3) |
and Neumann or nature condition for the potential flow:
|
(4) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{u}\left( X\right) \,}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{q}\left( X\right)}
are the boundary parts, where the Dirichlet and Neumann conditions respectively are prescribed, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathit{\boldsymbol{n}}}
is the normal related to the boundary Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\Gamma }_{q}\left( X\right)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {U}_{0}}
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \overline{U}\left( \mathit{\boldsymbol{X}},t\right)} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \overline{Q}\left( \mathit{\boldsymbol{X}},t\right) \,}
are prescribed functions.
3. Numerical method
Based on the superposition principle of the linear systems, the solution of a nonhomogeneous diffusion problem with time-independent source terms can be represented as:
|
(5) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {U}_{h}}
is the homogeneous solution and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {U}_{p}}
is the particular solution. Depending on the time marching scheme, the particular solution should satisfy:
For the Laplacian operator
|
(6a) |
For the modified Helmholtz operator
|
(6b) |
and the homogeneous solution satisfy:
|
(7) |
In this paper, the particular solution is obtained by DRM with different radial basis functions (RBFs) and the homogeneous solutions are obtained by different MFS formulations.
3.1. Dual reciprocity method for particular solution
In order to obtain the particular solution Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {U}_{p}\left( \mathit{\boldsymbol{X}}\right) \, \,}
with DRM, the term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle F\left( \mathit{\boldsymbol{X}}\right)}
from the right-hand side of equation (6) is approximated by a RBF:
|
(8) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
is the number of collocation points, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\alpha }_{i}}
are the interpolation coefficients and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\varphi }_{j}\left( \mathit{\boldsymbol{X}}\right) =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \, \varphi \left( r\right) =\varphi (\left| X-{X}_{j}\right| )
are functions based on the radial distance between collocation points with reference in Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {X}_{j}}
. Depending on the method, the collocation points will be placed only in the domain [16] or both the domain and boundary [18,19].
After applying Eq. (8) at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle N}
collocation points, the interpolation coefficients Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\alpha }_{i}}
are calculated. The particular solution is expressed:
|
(9) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\psi }_{i}(\mathit{\boldsymbol{X}})}
is the approximate particular solutions in terms of RBF Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\varphi }_{i}(\mathit{\boldsymbol{X}})}
. In order to obtain these solutions the following equations should be satisfied:
|
(10a) |
|
(10b) |
The typical approach is that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varphi}
is selected first, and then using Eq. (10) the corresponding particular solutions Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi}
are determined. For Eq. (10a) the particular solution in obtained by inverting the Laplace operator, only by repeated integration. For Eq. (10b), the particular solutions for the modified Helmholtz operator is not easy to obtain directly [20]. In this paper the RBF used was the thin plate splines (TPS) [1,20-22]:
|
(11) |
and with solving the following linear system:
|
(12) |
all unknown coefficients Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \alpha}
can be calculated, the corresponding Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \psi}
for the Laplacian operator:
|
(13) |
and for the modified Helmholtz operator:
|
(14a) | |
|
(14b) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \gamma =0.5772156649015328} is Euler’s constant and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {K}_{0}}
is the modified Bessel function of the second kind order zero.
The next step is to calculate the homogenous solution with the modified conditions presented in Eq. (7). Here we consider two time-marching typical methods: the finite difference method and transient fundamental solutions, and also, for this last one, it is proposed a different approach considering the time collocation for the source points.
3.2. Modified Helmholtz operator for homogeneous solution
Using the finite difference method in Eq. (1) to approximate Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{\partial U(X,t)}{\partial t}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle U\left( X,t\right)}
with respect to time variable Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t}
within the interval Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \left[ {t}^{n},{t}^{n+1}\right] \subset [0,T]}
|
(15) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t={t}^{n+1}-{t}^{n}}
is time step size and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \theta =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1
defining an implicit time marching method.
Substituting and rearranging terms of Eq. (16) in Eq. (1) and in conditions Eqs. (3)-(4) it yields to the modified Helmholtz equation:
|
(16) |
|
(17) |
Reducing Eq. (16) for a single-step formula:
|
(18) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\lambda }^{2}=\sqrt{\frac{1}{k\Delta t}}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A\left( \mathit{\boldsymbol{X}}\right) =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): -\frac{{U}^{n}\left( \mathit{\boldsymbol{X}}\right) }{k\Delta t}-\frac{F(\mathit{\boldsymbol{X}})}{k}
. It is worth noticing that the term Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A\left( \mathit{\boldsymbol{X}}\right)}
is well defined and depends on the solution obtained from the previous time and the nonhomogenous condition, which is time-independent and is calculated by DRM proposed in Section 3.1. The initial condition should be applied in order to start the procedure.
By using MFS for every time step, the homogenous solution is approximated by:
|
(19) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\beta }_{j}}
are unknown coefficients, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, {N}_{b}}
are the source points placed only in the boundary outside the problem domain, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {U}_{j}^{\ast }}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{{\partial U}_{j}^{\ast }\left( \mathit{\boldsymbol{X}}\right) }{\partial n}}
are the fundamental solutions proposed for the modified Helmholtz operator:
|
(20) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {K}_{0}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {K}_{1}}
are the modified Bessel function of the second kind order zero and one, respectively.
By properly placing the source points (Figure 1) in order to avoid singularities, a well posed linear system with dimension Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle ({N}_{b}x{N}_{b})}
is formed:
|
(21) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {A}_{ij}}
is the matrix containing the fundamental solution (Eq. (20)) on an Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle i}
th field point generated by a Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, j} th source and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {d}_{i}}
are the prescribed boundary conditions (Eq. (7)).
By computing Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\beta }_{j}} , the field points inside the domain can be obtained by Eq. (19). It is noteworthy that, for this method, sources and field points are at the same time level but in different geometrical positions, and it is well known that the source positioning outside the domain affects accuracy and stability. Typically, the response accuracy increase as the distance between the virtual and physical boundaries increase, but there is no definitive approach to this problem. In this article the sources geometrical position will be briefly evaluated.
|
| Figure 1. Source positions for MSF based on modified Helmholtz operator |
3.3. Transient fundamental solutions for homogeneous solution
Applying the same concept described by Eq. (19), the general solution for Eq. (1) is now described as:
|
(22) |
it is known Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {N}_{i}}
is the number of initial sources (existing only in the domain), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {N}_{b}\,}
is the number of sources in the boundary and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\beta }_{j}}
are the coefficients to be determined that represent the sources intensity. Inside the time step, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t=}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {t}^{n+1}-{t}^{n} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t}
is the impulse application time at the source point and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {t}_{F}}
is the field time, which are computed at intermediated levels for domain points and boundary points.
The fundamental solutions now are Green functions, which describe the potential field generated by a unit source applied to time Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t}
and point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \boldsymbol{\xi }}
|
(23) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \mathit{\boldsymbol{X}}} is the field point coordinate, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, \boldsymbol{\xi }}
is the source point coordinate, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, \tau =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {t}_{F}-t , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle d}
is the model dimension, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, r\left( X,\xi \right)}
is the radius measured between the source point and the field point in question, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \frac{\partial r\left( X,\xi \right) }{\partial n}}
is the projection of the radius projection over the normal-boundary.
In order to determine the sources intensities, the linear system must be solved:
|
(24) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {A}_{ij}}
is a matrix with dimension (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {N}_{i}+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {N}_{b})\, x\, ({N}_{i}+{N}_{b})
containing the fundamental solutions, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {d}_{i}}
is the vector containing the initial conditions and prescribed boundary conditions (Eq. (7)).
When using fundamental solutions that are already transient, source and field points are at same geometrical positions but at different time levels even though inside the same time step. As the matrix is ill-conditioned, the difference between the time to which the source is applied in comparison with the time the field is calculated generates significant effect on the final response and will be investigated in this article.
The uniqueness of the matrix was approached through the technical use of Tikhonov’s Method:
|
(25) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \lambda >0}
is a regularization parameter whose values generally vary between Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {10}^{-14}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {10}^{-1}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle I\,}
is the identity matrix.
3.3.1. Time sources positioning
Through Eq. (24) it is noted that both the matrix size and the vector size depend on the number of field points and source points chosen at each time step. However, it should be known that each time step is divided into intermediate steps where, firstly, the coefficients are determined by the points which the initial conditions are prescribed and then by the points where the boundary conditions are known.
In relation to the temporal coordinate of these sources, two forms of distribution would be tested: fixed source points (MFS1) and time moving source points (MFS2).
- 1) MFS1
In this paper, the scheme MFS1 is referred for fixed time source applications. On this proposal, even though the field points are obtained in a time box that advances through Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle [0,T]} , the sources impulses are applied on the problem beginning, always at a same time level for domain sources and at another level for boundary sources. In order to avoid singularities, intermediate time levels are proposed as: the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (0-} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): b)\Delta t
for domain sources, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (0+b)\Delta t}
for boundary sources, while field points will be placed at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n-}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1)\Delta t
in interior domain and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n+1)\Delta t}
at the boundary. Figure 2 shows source placement during time steps.
|
|
|
| Figure 2. MSF1 evolutionary process | ||
As the diffusion problem is an evolutionary problem, the impulse generated by maintaining the sources fixed at the beginning of the problem will be weakened as time progresses. Analogue to the geometrical distance, approximation increases accuracy as time evolves.
The matrix Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
exists when sources are applied before field nodes are calculated. For this approach, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
will be fully completed, once domain source points generate fundamental solutions at domain and boundary field points, as well as, after the first time step, boundary source points also create responses on domain and boundary field points. After computing Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \, A}
and solving the linear system as in Eq. (24), the numerical response can be obtained using Eq. (22).
- 2) MFS2
The scheme here referred as MSF2 makes use of [20,21] time-moving sources proposal. During the diffusion process, the field and source points are inside a time box that advances through Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle [0,T]} , both field and sources are at different time levels on every time step. The time levels proposed are: the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n-} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): b)\Delta t
for domain sources, and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n+b)\Delta t}
for boundary sources while field points will be placed at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n-}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1)\Delta t
in interior domain and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle (n+1)\Delta t}
at the boundary as showed in Figure 3.
|
|
|
| Figure 3. MSF2 evolutionary process | ||
Nevertheless, the Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
matrix will assume a different approach. Once boundary sources are applied after domain field points, it can only influence boundary field points. Inside a time step, the coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b}
that translates the distance between sources and field will now have an important role on the response accuracy. As in MSF1 the linear system Eq. (24) will be solved, and the potential Eq. (22) can be calculated at every time step until either the terminal time or the steady state solution is achieved.
4. Numerical results
In order to explore the methodologies potential and compare them with each other and with analytical solutions, three numerical examples were tested. As a way to evaluate the factors that influence the accuracy, several time increments, geometries, time coefficients and source positions were analyzed. Three types of internal sources were considered: constant, polynomial and senoidal for all proposed methods.
4.1. Example 1
Consider a unit circular disk, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k=1}
with the following conditions:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): I.C.:\, U\left( x,y,t=0\right) =0 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B.C.: \, U\left( x,y,t\right) =0 |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): F\left( x,y\right) =\, \frac{1}{m}=10 |
The analytical solution is given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): U\left( r,t\right) =\frac{1-{r}^{2}}{4m}-\frac{2}{m}\sum _{n=1}^{\infty }\frac{{J}_{0}({\lambda }_{n}r)}{{\lambda }_{n}^{3\, }{J}_{1}({\lambda }_{n})}{e}^{-k{\lambda }_{n}^{2}t} |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {J}_{0}}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {J}_{1}}
are the Bessel functions first kind zero and one order, respectively. The parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {\lambda }_{n}}
are the positive roots of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle {J}_{0}\left( {\lambda }_{n}\right) =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0
and in this paper was used the first 100 roots.
In order to explore all methodologies nuances, three points in different geometrical positions (0,0), (0.4,0) and (0.8,0) were studied. In order to study the influence of source positions for the modified Helmholtz operator, different distances between virtual and physical boundary were proposed in Figure 4, and for MSF1 and MSF2 three different coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b}
were analyzed.
|
| Figure 4. Physical and virtual boundary for example 1 |
Figure 5 depicts the absolute error histograms for the modified Helmholtz operator for all the distances shown in Figure 4. As expected, it can be observed for all analyzed nodes that error and virtual boundary have an inverse relation, meaning that, until a distance is achieved, the error decrease as the distance increase. The method accuracy is also slightly better for nodes geometrically positioned closer to the physical boundary.
|
| (a) |
|
| (b) |
|
| (c) |
| Figure 5. Time history of absolute error for example 1. Nodes coordinates are: (a) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.8
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): y=0 . (b) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.4 , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): y=0 . (c) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0,\, y=0 |
Figure 6 shows the absolute error histograms for MSF1 and MSF2 contemplating different time distances between sources and field nodes. In this example, for both methods, error decreases as the distance increases for all nodes. This relation shall not be generalized for different problems. Different parameters Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b}
should be tested for every problem. Joined with the time increment size, it affects Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle A}
directly. It is worth mentioning that, for all nodes, MSF1 is less sensitive to the time distance between source and field and presented better performance when compared to MSF2.
|
| (a) |
|
| (b) |
|
| (c) |
| Figure 6. Time history of absolute error for MSF1 and MSF2 for example 1. Nodes coordinates are: (a) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.8
, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): y=0 . (b) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.4 , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): y=0 . (c) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0,\, y=0 |
Figure 7 shows the time evolution history at (0,0), (0.4,0) and (0.8,0) for MSF1, MSF2 and MFS based on the modified Helmholtz operator, as well as the exact solution. For a circular domain, constant initial, boundary and nonhomogeneous conditions, accuracy decreases between the following methods: MSF based on the modified Helmholtz operator, MSF1 and MSF2. Although, in general, the results show good agreement with the analytical solution in different time stages.
|
| Figure 7. Example 1 numerical results compared with analytical response |
4.2. Example 2
Let a square plate of dimension Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 1 \times 1 , with 256 nodes, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \Delta t=} Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): 0.2s , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k=1}
be submitted to the conditions:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): I.C.: \,U\left( x,y,t=0\right) =2\left[ \cos\left( \pi x\right) + \sin\left( \pi y\right) \right] +
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{x}^{3}+{y}^{3}+{x}^{2}+y}{12} |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B.C.:\, U\left( x,y,t\right) =2\left[ \cos\left( \pi x\right) + \sin\left( \pi y\right) \right] {e}^{-k{\pi }^{2}t}+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{x}^{3}+{y}^{3}+{x}^{2}+y}{12} |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): F\left( x,y\right) =\, \frac{1}{m}=-\frac{6x+6y+2}{12} |
The analytical solution is given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): U\left( x,y,t\right) =2\left[ \cos\left( \pi x\right) + \sin\left( \pi y\right) \right] {e}^{-k{\pi }^{2}t}+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{x}^{3}+{y}^{3}+{x}^{2}+y}{12} |
In this example, for MSF1 and MSF2, the coefficient Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle b=0.5} , and, for the modified Helmholtz operator, virtual sources are 0.5 from the physical boundary. In order to explore the methodologies, two different types of RBF were considered, referenced as RBF1 the TPS described on Section 3.1 and RBF2 the following functions [23]:
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varphi \left( r\right) =\, 1+r\quad}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \quad\psi \left( r\right) =\frac{{r}^{2}}{4}+}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{r}^{3}}{9}\qquad for Laplacian operator |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \varphi \left( r\right) =\, 1+r-{\lambda }^{2}(\frac{{r}^{2}}{4}+\frac{{r}^{3}}{9})\quad}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \quad\psi \left( r\right) =}
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): \frac{{r}^{2}}{4}+\frac{{r}^{3}}{9}\qquad for modified Helmholtz operator |
Figure 8 depicts the absolute error at nodes with distinct geometric positions for all methods proposed with RBF1 and RBF2. For all nodes, at initial stages, MSF1 and MSF2 presented better results. Moreover, between the radial basis functions, RBF1 showed greater accuracy. For senoidal and exponential initial and boundary conditions, polynomial non-homogeneous condition, MFS2-RBF1 revealed less error the closer the node is to the physical boundary, and MSF1 and MSF based on the modified Helmholtz operator revealed to be better the closer the node is to center domain.
|
| (a) |
|
| (b) |
|
| (c) |
| Figure 8. Time history of absolute error for example 2. Nodes coordinates are: (a) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.1,\, y=0.1
. (b) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.5,\, y=0.5 . (c) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): x=0.9, \, y=0.9 |
Figure 9 presents the results obtained by the three tested methods and the analytical solution for a line of nodes at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle y=0.5}
on initial stage and final stage. For initial stages MFS based on modified Helmholtz operator present less accuracy than MFS1 and MFS2, while at final times all methods present good precision. Figure 10 shows potential at all nodes at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle t=3s}
.
|
|
| (a) | (b) |
| Figure 9. Example 2 numerical results compared with analytical response. (a) Initial stage. (b) Final stage | |
|
| Figure 10. Example 2 potential for all nodes |
4.3. Example 3
Consider a L-shaped irregular meshless MFS nodes, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle k=1}
with the following conditions:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): I.C.:\, U\left( x,y,t=0\right) =\mathrm{cos}\,\left( 3x\right) \mathrm{sinh}\,\left( 2y\right) xsin\left( y\right) |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): B.C.: \, U\left( x,y,t\right) =\mathrm{cos}\,\left( 3x\right) \mathrm{sinh}\,\left( 2y\right) \mathrm{exp}\,(-5t)+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): xsin(y) |
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): F\left( x,y\right) =\, xsin(y) |
The analytical solution is given by:
| Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): U\left( x,y,t\right) =\mathrm{cos}\,\left( 3x\right) \mathrm{sinh}\,\left( 2y\right) \mathrm{exp}\,(-5t)+
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): xsin(y) |
In this example, two nodes –point A (1.03,0.56) and point B (0.35,1.36) depicted in Figure 11– and different time increments were tested to verify its influence on the final response.
|
| Figure 11. Meshless MFS nodes for example 3 |
Figure 12 depicts absolute errors at nodes A and B, for different time steps in all methods. When compared, MSF1 and MSF based on the modified Helmholtz operator present better responses as time step decreases, and MSF2, a method that is known for its ill posed matrix, presents better results for a bigger time step. All methods showed to be sensitive concerning the field node geometrical position.
|
|
| (a) | (b) |
|
|
| (c) | (d) |
| Figure 12. Time history of absolute error for example 3. (a) MSF1 and MSF2 for point A. (b) MSF1 and MSF2 for point B. (c) Helmholtz for point A. (d) Helmholtz for point B | |
Figure 13 presents the results obtained by the three tested methods and the analytical solution at nodes A and B. For initial stages MFS1 presents less accuracy than the other methods. However, at final times MFS2 is worst. In general, for a senoidal non-homogeneous condition and trigonometrical initial and boundary conditions, all methods present reasonable results.
|
|
| (a) | (b) |
| Figure 13. Example 3 numerical results compared with analytical response. (a) Node A. (b) Node B | |
5. Conclusions
Transient nonhomogeneous diffusion problems with different MFS methods and several types of initial, boundary and time-independent source functions were solved. The DRM technique was used to solve the nonhomogeneous equation, and different RBFs were analyzed. For the homogeneous solution, circular, squared and hard shaped meshes were tested, as well as different time increments and source positioning. It was described MFS based on the modified Helmholtz operator, MFS based on transient fundamental solutions, and it was proposed a variation of this last method.
The numerical procedures presented were validated by comparing results between methods and also with analytical solutions. The main issue for all methods is the sources positioning. For MFS modified Helmholtz operator, sources are geometrically positioned at a virtual boundary, and for MFS modified transient fundamental solutions, sources are time positioned at intermediate stages. In the proposed method sources are fixed in time, and its intensities, translated at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://mathoid.scipedia.com/localhost/v1/":): {\textstyle \beta}
coefficients, are calculated for every time step.
Regarding the geometrical source position, the Helmholtz-based method showed to be sensitive until a certain distance between virtual and physical boundary is achieved, pattern that repeats for all mesh types. The so-called MSF1 showed to be less sensitive then MSF2 to temporal source positioning, although shape complexity interfered with the accuracy in both methods. Concerning time increments, MSF2 presented great discrepancies in accuracy between the times studied, mainly due to the disturbances generated when solving an ill-conditioned system, and for MFS1 and Helmholtz method, smaller time increments increase accuracy. The developed methods have a formulation that proved to be efficient in solving the various potential problems. Despite their limitations, it was realized that the methods present a controllable level of precision, versatility, computational gain and easy implementation.
Acknowledgements
The authors thankfully acknowledge the support of Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq). This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001.
References
[1] Gingold R.A., Monaghan J.J. Smoothed particle hydrodynamics: theory and application to non-spherical stars. Monthly notices of the royal astronomical society, 181(3):375-389, 1977.
[2] Oñate E., Idelsohn S., Zienkiewicz O.C., Taylor R.L. A finite point method in computational mechanics. Applications to convective transport and fluid flow. International Journal for Numerical Methods in Engineering, 39(22):3839-3866, 1996.
[3] Belytschko T., Lu Y.Y., Gu L. Element‐free Galerkin methods. International Journal for Numerical Methods in Engineering, 37(2):229-256, 1994.
[4] Atluri S.N., Zhu T. A new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics. Computational Mechanics, 22(2):117-127, 1998.
[5] Chen W. Meshfree boundary particle method applied to Helmholtz problems. Engineering Analysis with Boundary Elements, 26(7):577-581, 2002.
[6] Wang H., Qin Q.H. Methods of fundamental solutions in solid mechanics. Elsevier, 2019.
[7] Young D.L., Tsai C.C., Murugesan K., Fan C.M., Chen C.W. Time-dependent fundamental solutions for homogeneous diffusion problems. Engineering Analysis with Boundary Elements, 28(12):1463-1473, 2004.
[8] Zhu S.P. Solving transient diffusion problems: time-dependent fundamental solution approaches versus LTDRM approaches. Engineering Analysis with Boundary Elements, 21(1):87-90, 1998.
[9] Zhu S.P., Liu H.W., Lu X.P. A combination of LTDRM and ATPS in solving diffusion problems. Engineering Analysis with Boundary Elements, 21(3):285-289, 1998.
[10] Sutradhar A., Paulino G.H., Gray L.J. Transient heat conduction in homogeneous and non-homogeneous materials by the Laplace transform Galerkin boundary element method. Engineering Analysis with Boundary Elements, 26(2):119-132, 2002.
[11] Bulgakov V., Šarler B., Kuhn G. Iterative solution of systems of equations in the dual reciprocity boundary element method for the diffusion equation. International Journal for Numerical Methods in Engineering, 43(4):713-732, 1998.
[12] Golberg M.A., Chen C.S., Muleshkov A.S. The method of fundamental solutions for time-dependent problems. WIT Transactions on Modelling and Simulation, 23, 1970.
[13] Nardini D., Brebbia C.A. A new approach to free vibration analysis using boundary elements. Applied Mathematical Modelling, 7(3):157-162, 1983.
[14] Partridge P.W., Brebbia C.A. Dual reciprocity boundary element method. Springer Science & Business Media, 2012.
[15] Golberg M.A., Chen C.S., Rashed Y.F. The annihilator method for computing particular solutions to partial differential equations. Engineering Analysis with Boundary Elements, 23(3):275-279, 1999.
[16] Cao L., Qin Q.H., Zhao N. Application of DRM-Trefftz and DRM-MFS to transient heat conduction analysis. Recent Patents on Space Technology, 2:41-50, 2010.
[17] Chantasiriwan S. Methods of fundamental solutions for time‐dependent heat conduction problems. International Journal for Numerical Methods in Engineering, 66(1):147-165, 2006.
[18] Young D.L., Tsai C.C., Fan C.M. Direct approach to solve nonhomogeneous diffusion problems using fundamental solutions and dual reciprocity methods. Journal of the Chinese Institute of Engineers, 27(4):597-609, 2004.
[19] Young D.L., Chen C.W., Fan C.M., Tsai C.C. The method of fundamental solutions with eigenfunction expansion method for nonhomogeneous diffusion equation. Numerical Methods for Partial Differential Equations: An International Journal, 22(5):1173-1196, 2006.
[20] Golberg M.A., Chen C.S. The theory of radial basis functions applied to the BEM for inhomogeneous partial differential equations. Boundary Elements Communications, 5(2):57-61, 1994.
[21] Chen C.S., Rashed Y.F. Evaluation of thin plate spline based particular solutions for Helmholtz-type operators for the DRM. Mechanics Research Communications, 25(2):195-201, 1998.
[22] Golberg M.A., Chen C.S. The method of fundamental solutions for potential, Helmholtz and diffusion problems. Boundary Integral Methods: Numerical and Mathematical Aspects, 1:103-176, 1998.
[23] Balakrishnan K., Sureshkumar R., Ramachandran P. A. An operator splitting-radial basis function method for the solution of transient nonlinear Poisson problems. Computers & Mathematics with Applications, 43(3-5):289-304, 2002.
Document information
Published on 09/07/20
Accepted on 27/06/20
Submitted on 21/01/20
Volume 36, Issue 3, 2020
DOI: 10.23967/j.rimni.2020.07.003
Licence: CC BY-NC-SA license
Share this document
Keywords
claim authorship
Are you one of the authors of this document?