On the Analytical Inversion of Solver Matrices Used in Numerical Approximations for the Diffusion Equation

Abstract:

The goal of this paper is to introduce an analytical approach for the inversion of nxn solver matrices, which are typically used in Finite Difference Method approximations. In the present case, they are used to solve the Diffusion Equation numerically, since in many physics and engineering fields, partial differential equations cannot be solved analytically. The method presented in this work is primarily formulated for cylindrical coordinates, which are often used in Gas Release Experiments as those described in [8]. However, it is possible to introduce a generalized method, which also allows solutions for Cartesian solvers. The advantage of having the explicit inverse is considerable, since the computational effort is reduced. In this paper we also carry out an investigation on the eigenvalues of the backward and forward solver matrix in order to determine an optimal range for the discretization parameters.

• Periodical: Defect and Diffusion Forum (Volume 413)
• Pages: 3-18
• DOI: https://doi.org/10.4028/www.scientific.net/DDF.413.3
• Online since: December 2021
• Authors: Till Glage*, Axel von der Weth, Frederik Arbeiter, Daniela Piccioni Koch
• Keywords: Fusion Research, Gas Release Experiment, Hydrogen Diffusion, Matrix Solver, Partial Differential Equation (PDE)

References

[1] Bickley, W. G.;Thompson, R. S. H. G., Matrices Their meaning and manipulation, The English Universities Press Ltd., London, (1964).

[2] Fischer, Gerd, Lineare Algebra, 7.Auflage, Friedr. Vieweg & Sohn Braunschweig/Wiesbaden, (1981).

[3] Kulkarni, Devadatta; Schmidt, Darrell; Tsui, Sze-Kai. Eigenvalues of pseudo-Toeplitz matrices, Department of Mathematical Science, Oakland University, Rochester.

[4] Schulz, Marvin R. et al., Analytical Solution of a Gas Release problem considering permeation with time-dependent boundary conditions, Journal of Computational and Theoretical Transport, (2021).

DOI: 10.1080/23324309.2020.1828469

[5] Schwarz, H. R., Numerische Mathematik, 2. Auflage, B. G. Teubner, Stuttgart, (1988).

[6] Smith, G. D. Numerical Solution of Partial Differential Equations, Oxford University Press, (1969).

[7] Bronstein, I.N.; Semendjajew K. A.; Taschenbuch der Mathematik, 23. Auflage; Grosche, G.; Ziegler, V.; Ziegler, D. (publishers); Verlag Harri Deutsch Thun und Frankfurt/Main (1987).

DOI: 10.1002/piuz.19800110415

[8] von der Weth, A., DSL 2021 Paper: Numerical Solution Strategies in Permeation Processes.

[9] Klimenko, D., Arbeiter, F., Pasler, V., Schlindwein, G., von der Weth, A. and Zinn, K., Definition of the Q-PETE experiment for investigation of hydrogen isotopes permeation through the metal structures of a DEMO HCPB breeder zone, Fusion engineering and design, vol. 136A, 563-568, (2018).

DOI: 10.1016/j.fusengdes.2018.03.024

[10] Axel von der Weth, A., Arbeiter, F., Klimenko, D., Pasler, V., Schlindwein, G., Permeation data analysis including a nonzero hydrogen concentration on the low pressure detector side for a purged permeation experiment, Presentation, DSL Conference 2018, Amsterdam, 27th June (2018).

DOI: 10.4028/www.scientific.net/ddf.391.18