Numerical solution of the modified Bessel equation

DOI10.1093/IMANUM/DRS031zbMATH Open1271.65144arXiv1110.2089OpenAlexW2161110508MaRDI QIDQ2841065FDOQ2841065

Authors: Michael Carley 

Publication date: 24 July 2013

Published in: IMA Journal of Numerical Analysis (Search for Journal in Brave)

Abstract: A Green's function based solver for the modified Bessel equation has been developed with the primary motivation of solving the Poisson equation in cylindrical geometries. The method is implemented using a Discrete Hankel Transform and a Green's function based on the modified Bessel functions of the first and second kind. The computation of these Bessel functions has been implemented to avoid scaling problems due to their exponential and singular behavior, allowing the method to be used for large order problems, as would arise in solving the Poisson equation with a dense azimuthal grid. The method has been tested on monotonically decaying and oscillatory inputs, checking for errors due to interpolation and/or aliasing. The error has been found to reach machine precision and to have computational time linearly proportional to the number of nodes.

Full work available at URL: https://arxiv.org/abs/1110.2089

Recommendations

• Computing solutions of the modified Bessel differential equation for imaginary orders and positive arguments
• On some integrals from modified Bessel functions
• Mathematical software for modified Bessel functions
• Publication:4867821
• A code evaluate modified Bessel functions based on the continued fraction method

zbMATH Keywords

error boundsnumerical examplesPoisson equationbiharmonic equationmodified Bessel functiondiscrete Hankel transformcylindrical geometriesGreen's-function-based solver

Mathematics Subject Classification ID

Computation of special functions and constants, construction of tables (65D20) Laplace operator, Helmholtz equation (reduced wave equation), Poisson equation (35J05) Bessel and Airy functions, cylinder functions, ({}_0F_1) (33C10) Boundary value problems for higher-order elliptic equations (35J40) Error bounds for boundary value problems involving PDEs (65N15) Boundary element methods for boundary value problems involving PDEs (65N38) Numerical approximation and evaluation of special functions (33F05)

Cited In (1)

• Computation of potentials from current electrodes in cylindrically stratified media: a stable, rescaled semi-analytical formulation

This page was built for publication: Numerical solution of the modified Bessel equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2841065)