Boundary value problems and medical imaging

2nd International Conference on Mathematical Modeling in Physical Sciences 2013

物理学;数学

Fokas, Athanasios S.^1,2 ; Kastis, George A.^2

Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB30WA, United Kingdom^1

Research Center of Mathematics, Academy of Athens, Soranou Efessiou 4, Athens 11527, Greece^2

关键词: Complex planes;    Infinite series;    Mellin transform;    Non-linear PDEs;    Non-self-adjoint;    Numerical computations;    Numerical techniques;    Uniformly convergent;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/490/1/012017/pdf
DOI  :  10.1088/1742-6596/490/1/012017

来源: IOP

【 摘 要 】

The application of appropriate transform pairs, such as the Fourier, the Laplace, the sine, the cosine and the Mellin transforms, provides the most well known method for constructing analytical solutions to a large class of physically significant boundary value problems. However, this method has several limitations. In particular, it requires the given PDE, domain and boundary conditions to be separable, and also may not be applicable if the given boundary value problem is non-self-adjoint. Furthermore, it expresses the solution as either an integral or an infinite series, neither of which are uniformly convergent on the boundary of the domain (for nonvanishing boundary conditions), which renders such expressions unsuitable for numerical computations. Here, we review a method recently introduced by the first author which can be applied to certain nonseparable and non-self-adjoint problems. Furthermore, this method expresses the solution as an integral in the complex plane which is uniformly convergent on the boundary of the domain. This method, which also suggests new numerical techniques, is illustrated for both evolution and elliptic PDEs. Athough this method was first applied to certain nonlinear PDEs called integrable and was originally formulated in terms of the so-called Lax pairs, it can actually be applied to linear PDEs without the need to analyse the associated Lax pair. The existence of Lax pairs is used here in order to motivate a related development, namely the emergence of a novel formalism for analysing certain inverse problems arising in medical imaging. Examples include PET and SPECT.