The Numerical Solution of the External Dirichlet Generalized Harmonic Problem for a Sphere by the Method of Probabilistic Solution

Mamuli Zakradze,Zaza Tabagari,Nana Koblishvili,Tinatin Davitashvili,Jose Maria Sanchez,Francisco Criado-Aldeanueva
DOI: https://doi.org/10.3390/math11030539
2024-05-30
Abstract:In the present paper, an algorithm for the numerical solution of the external Dirichlet generalized harmonic problem for a sphere by the method of probabilistic solution (MPS) is given, where generalized indicates that a boundary function has a finite number of first kind discontinuity curves. The algorithm consists of the following main stages: (1) the transition from an infinite domain to a finite domain by an inversion; (2) the consideration of a new Dirichlet generalized harmonic problem on the basis of Kelvin theorem for the obtained finite domain; (3) the numerical solution of the new problem for the finite domain by the MPS, which in turn is based on a computer simulation of the Weiner process; (4) finding the probabilistic solution of the posed generalized problem at any fixed points of the infinite domain by the solution of the new problem. For illustration, numerical examples are considered and results are presented.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the exterior Dirichlet generalized harmonic problem of solving the boundary function with the first - kind discontinuous curves on the external sphere. Specifically, the paper focuses on how to effectively solve this kind of problem numerically in the infinite domain, especially when the boundary function has the first - kind discontinuities on a finite number of curves. Traditional numerical methods often encounter difficulties in dealing with this type of problem, such as slow convergence speed or instability. Therefore, new algorithms need to be developed to improve the accuracy and stability of the solution. The paper proposes an algorithm based on the Probabilistic Solution (MPS), which is implemented through the following steps: 1. **Transformation from infinite domain to finite domain**: Transform the original problem from the infinite domain to the finite domain through inversion techniques. 2. **Construction of a new problem based on Kelvin's theorem**: Apply Kelvin's theorem on the transformed finite domain to reconstruct a new Dirichlet generalized harmonic problem. 3. **Application of the probabilistic solution**: Use the Probabilistic Solution (MPS) to numerically solve the new problem, which involves computer simulation of the Wiener process. 4. **Solving the original problem**: Use the solution of the new problem to find the probabilistic solution of the original problem at any fixed point in the infinite domain. The paper also verifies the effectiveness and accuracy of the proposed algorithm through specific numerical examples, showing the potential of this method in practical applications.