Abstract
We present a new algorithm of the finite element method (FEM) implemented as KANTBP 5M code in MAPLE for solving boundary-value problems (BVPs) for systems of second-order ordinary differential equations with continuous or piecewise continuous real or complex-valued coefficients. The desired solution in a finite interval of the real-valued independent variable is subject to mixed homogeneous boundary conditions (BCs). To reduce a BVP or a scattering problem with different numbers of asymptotically coupled or entangled open channels in the two asymptotic regions to a BVP on a finite interval, the asymptotic BCs for large absolute values of the independent variable are approximated by homogeneous Robin BCs. The BVP is discretized by means of the FEM using the Hermite interpolation polynomials with arbitrary multiplicity of the nodes, which preserves the continuity of derivatives of the desired solutions. The relevant algebraic problems are solved using the built-in linear algebra procedures. To calculate metastable states with complex eigenvalues of energy or to find bound states with the BCs depending on a spectral parameter, the Newton iteration scheme is implemented. Benchmark examples of the code application to BVPs and scattering problems of quantum mechanics are given.
V. Gerdt—It is painful to think Professor V. Gerdt is no longer among us, and this paper is his last contribution to development of solving BVPs on the base of the FEM which owes remarkable results to him. We are deeply grateful to him for his intuition, insight and support, which were invaluable during our long-standing collaboration.
This is a preview of subscription content, log in via an institution to check access.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Similar content being viewed by others
References
Alder, B., Fernbach, S., Rotenberg, M. (eds.): Atomic and Molecular Scattering: Methods In Computational Physics. Academic Press, New York, London (1971)
Katsenelenbaum, B.Z., Mercader del Rio, L., Pereyaslavets, M., Sorolla Ayza, M., Thumm, M.: Theory of Nonuniform Waveguides the Cross-Section Method. The Institution of Electrical Engineers (1998)
Brekhovskikh, L.M., Lysanov, Y.P.: Fundamentals of Ocean Acoustics. Springer-Verlag, Berlin (2003)
Ramdas Ram-Mohan, L.: Finite Element and Boundary Element Applications in Quantum Mechanics. Oxford University Press, New York (2002)
Kantorovich, L.V., Krylov, V.I.: Approximate Methods of Higher Analysis. Wiley, New York (1964)
Berezin, I.S., Zhidkov, N.P.: Computing Methods. Pergamon Press, Oxford (1965)
Derbov, V.L., et al.: Spectrum of beryllium dimer in ground \(X ^1\varSigma ^{+}_{g}\) state. J. Quant. Spectrosc. Radiat. Transf. 262, 107529–1-10 (2021). https://doi.org/10.1016/j.jqsrt.2021.107529
Wen, P.W., et al.: Near-barrier heavy-ion fusion: role of boundary conditions in coupling of channels. Phys. Rev. C 101, 014618-1-10 (2020). https://doi.org/10.1103/PhysRevC.101.014618
Divakov, D.V., Tiutiunnik, A.A., Sevastianov, A.L.: Symbolic-numeric study of geometric properties of adiabatic waveguide modes. In: Boulier, F., England, M., Sadykov, T.M., Vorozhtsov, E.V. (eds.) CASC 2020. LNCS, vol. 12291, pp. 228–244. Springer, Cham (2020). https://doi.org/10.1007/978-3-030-60026-6_13
Gusev, A.A., Gerdt, V.P., Hai, L.L., Derbov, V.L., Vinitsky, S.I., Chuluunbaatar, O.: Symbolic-numeric algorithms for solving BVPs for a system of ODEs of the second order: multichannel scattering and eigenvalue problems. In: Gerdt, V.P., Koepf, W., Seiler, W.M., Vorozhtsov, E.V. (eds.) CASC 2016. LNCS, vol. 9890, pp. 212–227. Springer, Cham (2016). https://doi.org/10.1007/978-3-319-45641-6_14
Chuluunbaatar, G., Gusev, A.A., Chuluunbaatar, O., Vinitsky, S.I., Hai, L.L.: KANTBP 4M: Program for solving the scattering problem for a system of ordinary second-order differential equations. EPJ Web Conf. 226, 02008 (2020). https://doi.org/10.1051/epjconf/202022602008
Gusev, A.A., et al.: Symbolic-numeric solution of boundary-value problems for the Schrödinger equation using the finite element method: scattering problem and resonance states. In: Gerdt, V.P., Koepf, W., Seiler, W.M., Vorozhtsov, E.V. (eds.) CASC 2015. LNCS, vol. 9301, pp. 182–197. Springer, Cham (2015). https://doi.org/10.1007/978-3-319-24021-3_14
Gusev, A.A., et al.: Symbolic-numerical solution of boundary-value problems with self-adjoint second-order differential equation using the finite element method with interpolation hermite polynomials. In: Gerdt, V.P., Koepf, W., Seiler, W.M., Vorozhtsov, E.V. (eds.) CASC 2014. LNCS, vol. 8660, pp. 138–154. Springer, Cham (2014). https://doi.org/10.1007/978-3-319-10515-4_11
Chuluunbaatar, O., et al.: KANTBP: A program for computing energy levels, reaction matrix and radial wave functions in the coupled-channel hyperspherical adiabatic approach. Comput. Phys. Commun. 177, 649–675 (2007). https://doi.org/10.1016/j.cpc.2007.05.016
Chuluunbaatar, O., Gusev, A.A., Abrashkevich, A.G., Vinitsky, S.I.: KANTBP 2.0: new version of a program for computing energy levels, reaction matrix and radial wave functions in the coupled-channel hyperspherical adiabatic approach. Comput. Phys. Commun. 179, 685–693 (2008). https://doi.org/10.1016/j.cpc.2008.06.005
Chuluunbaatar, O., Gusev, A.A., Abrashkevich, A.G., Vinitsky, S.I.: KANTBP 3.0: new version of a program for computing energy levels, reflection and transmission matrices, and corresponding wave functions in the coupled-channel adiabatic approach. Comput. Phys. Commun. 185, 3341–3343 (2014). https://doi.org/10.1016/j.cpc.2014.08.002
Gusev A.A., Hai L.L., Chuluunbaatar O., Vinitsky S.I.: KANTBP 4M: program for solving boundary problems of the system of ordinary second order differential equations. JINRLIB (2015). http://wwwinfo.jinr.ru/programs/jinrlib/kantbp4m/indexe.html
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Appendices
A Generation of IHPs on the Standard Interval
This appendix presents an algorithm for constructing IHPs according to their characteristics: p is the number of partitions of a finite element, \(z_r\) are the IHP nodes with multiplicities \(\kappa _r^{\max }\). They are applied to construct IHPs in the FEM scheme, then the conditions \(z_r\in [0,1]\), \(z_0=0\), \(z_p=1\), \(\kappa _0^{\max }=\kappa _p^{\max }\) are to be satisfied. For further implementation it is convenient to number IHPs with \(n''\).
B Calculation of the FEM Scheme Characteristics
Note that when calculating the matrices (9) of the algebraic problem (8), we do it without explicitly calculating \( N_s (z) \) from (7) by introducing global numbering \(\varphi _r^\kappa \) on each of the finite elements \(\varDelta _j \), i.e. \( \varphi _S \equiv \varphi _ {n''} (z \in \varDelta _j) \equiv \varphi _r ^ \kappa (z \in \varDelta _j) \). In our implementation, the FEM IHP schemes are numbered so that S increases with an increase in j, or with a constant j and an increase in \( n ''\), or with constant j and \( n'' \) and an increase in i. For convenience, arrays of length \( n \times 3 \) are introduced: E(j, 1) is the minimum S at which \( \varphi _S \) is defined on \( \varDelta _j \), E(j, 2) is the minimum S for which \( r = p \) and \( \varphi _S \) is defined on \( \varDelta _j \), E(j, 3) is the maximum S at which \( \varphi _S \) is defined on \( \varDelta _j \) and a two-dimensional array C with dimension \( S^{\max } \times 3 \), where depending on S, C(S, 1) , C(S, 2) , C(S, 3) correspond to \(\mu \) (the number of element of eigenvector \( \boldsymbol{\varPhi }^ h \)), \( n''\) (the number of IHP) and i (the number of equation in the system of ODEs from Eq. (1)).
C FEM generation of Algebraic Eigenvalue Problem
Rights and permissions
Copyright information
© 2021 Springer Nature Switzerland AG
About this paper
Cite this paper
Chuluunbaatar, G. et al. (2021). A Maple Implementation of the Finite Element Method for Solving Boundary-Value Problems for Systems of Second-Order Ordinary Differential Equations. In: Corless, R.M., Gerhard, J., Kotsireas, I.S. (eds) Maple in Mathematics Education and Research. MC 2020. Communications in Computer and Information Science, vol 1414. Springer, Cham. https://doi.org/10.1007/978-3-030-81698-8_11
Download citation
DOI: https://doi.org/10.1007/978-3-030-81698-8_11
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-030-81697-1
Online ISBN: 978-3-030-81698-8
eBook Packages: Computer ScienceComputer Science (R0)