We developed a Jacobi–spectral framework for the heat equation in a spherical domain under axial symmetry and Dirichlet boundary conditions. The angular part of the Laplacian was realized as a Jacobi Sturm–Liouville operator on a weighted $ L^{2} $ space, enabling the Jacobi transform to diagonalize the angular component and project the partial diferential equation (PDE) onto a sequence of decoupled radial problems. Each projected equation reduced to a Euler-type radial ordinary diferential equation (ODE) driven by the corresponding Jacobi coefficient of the source term. These modal equations were solved in terms of spherical-Bessel eigenfunctions and radial Green kernels, yielding explicit Duhamel-type formulas for the time-dependent coefficients and establishing convergence in the weighted $ L^{2} $ space. The Legendre case $ (\alpha, \beta) = (0, 0) $ recovered the classical axisymmetric model, while general Jacobi parameters provided a unified extension of this setting. A central result was the demonstration of a rigorous equivalence between the Jacobi–spectral representation and the classical separation-of-variables solution written in spherical harmonics and spherical-Bessel modes. The proposed framework clarified the angular–radial coupling in spherical geometries and connected naturally with modern Jacobi and ultraspherical spectral methods.
Citation: Juan Toribio Milane, José A. Gómez Hernández, Juan R. Holguín, Pedro N. Tifa de Jesús. A Jacobi–spectral framework for the heat equation with Dirichlet boundary conditions[J]. AIMS Mathematics, 2026, 11(3): 5776-5797. doi: 10.3934/math.2026238
We developed a Jacobi–spectral framework for the heat equation in a spherical domain under axial symmetry and Dirichlet boundary conditions. The angular part of the Laplacian was realized as a Jacobi Sturm–Liouville operator on a weighted $ L^{2} $ space, enabling the Jacobi transform to diagonalize the angular component and project the partial diferential equation (PDE) onto a sequence of decoupled radial problems. Each projected equation reduced to a Euler-type radial ordinary diferential equation (ODE) driven by the corresponding Jacobi coefficient of the source term. These modal equations were solved in terms of spherical-Bessel eigenfunctions and radial Green kernels, yielding explicit Duhamel-type formulas for the time-dependent coefficients and establishing convergence in the weighted $ L^{2} $ space. The Legendre case $ (\alpha, \beta) = (0, 0) $ recovered the classical axisymmetric model, while general Jacobi parameters provided a unified extension of this setting. A central result was the demonstration of a rigorous equivalence between the Jacobi–spectral representation and the classical separation-of-variables solution written in spherical harmonics and spherical-Bessel modes. The proposed framework clarified the angular–radial coupling in spherical geometries and connected naturally with modern Jacobi and ultraspherical spectral methods.
| [1] | G. Szegő, Orthogonal polynomials, American Mathematical Society, 23 1939. |
| [2] | A. P. Agudelo, J. G. López, Ecuación diferencial asociada a los polinomios ortogonales clásicos, Sci. Tech., 10 (2004), 179–184. Available from: http://www.redalyc.org/articulo.oa?id=84911640031. |
| [3] | N. Asmar, Partial differential equations, 2 Eds., Pearson Prentice Hall, 2004. |
| [4] | M. Tyn, L. Debnath, Linear partial differential equations for scientists and engineers, 4 Eds., Boston: Birkhäuser, 2007. |
| [5] | J. T. Milane, P. N. T. de-Jesús, Solución equivalente del problema de Dirichlet usando transformada de Legendre y funciones de Green, representada por la fórmula integral de Poisson, Lat.-Am. J. Phys. Educ., 15 (2021). |
| [6] | L. R. S. M. Jiménez, El problema $L^p$ Dirichlet para la ecuación de Laplace, México, Universidad de Sonora, 2011. |
| [7] | M. Lewin, Spectral theory and quantum mechanics, Springer Cham, 2024. https://doi.org/10.1007/978-3-031-66878-4 |
| [8] |
Z. Zhao, Z. H. Lu, Y. G. Zhao, Simulating non-stationary and non-Gaussian cross-correlated fields using multivariate Karhunen–Loève expansion and L-moments-based Hermite polynomial model, Mech. Syst. Signal Pr., 216 (2024), 111480. https://doi.org/10.1016/j.ymssp.2024.111480 doi: 10.1016/j.ymssp.2024.111480
|
| [9] |
F. Rösler, A. Stepanenko, Computing eigenvalues of the Laplacian on rough domains, Math. Comput., 93 (2024), 111–161. https://doi.org/10.1090/mcom/3827 doi: 10.1090/mcom/3827
|
| [10] |
J. Kileel, N. F. Marshall, O. Mickelin, A. Singer, Fast expansion into harmonics on the ball, SIAM J. Sci. Comput., 47 (2025), A1117–A1144. https://doi.org/10.1137/24M1668159 doi: 10.1137/24M1668159
|
| [11] | F. W. J. Olver, D. W. Lozier, R. F. Boisvert, C. W. Clark, NIST handbook of mathematical functions, New York: Cambridge University Press, 2010. |
| [12] | R. Álvarez-Nodarse, Polinomios ortogonales: Historia y aplicaciones, Bol. Soc. Esp. Mat. Apl., 18 (2001), 19–46. |
| [13] | J. C. Mason, D. C. Handscomb, Chebyshev polynomials, New York: Chapman & Hall/CRC, 2002. https://doi.org/10.1201/9781420036114 |
| [14] | L. Debnath, D. Bhatta, Integral transforms and their applications, 3 Eds., Chapman and Hall/CRC, 2014. https://doi.org/10.1201/b17670 |
| [15] |
J. Bremer, Fast algorithms for Jacobi expansions via nonoscillatory phase functions, IMA J. Numer. Anal., 40 (2020), 2019–2045. https://doi.org/10.1093/imanum/drz016 doi: 10.1093/imanum/drz016
|
| [16] |
Q. Yan, S. Jiang, J. Harlim, Spectral methods for solving elliptic PDEs on unknown manifolds, J. Comput. Phys., 486 (2023), 112132. https://doi.org/10.1016/j.jcp.2023.112132 doi: 10.1016/j.jcp.2023.112132
|
| [17] | G. B. Folland, Real analysis: Modern techniques and their applications, 2 Eds, New York: Wiley, 1999. |
| [18] |
J. Bremer, Q. Pang, H. Yang, Fast algorithms for the multi-dimensional Jacobi polynomial transform, Appl. Comput. Harmon. A., 51 (2021), 231–250. https://doi.org/10.1016/j.acha.2020.01.004 doi: 10.1016/j.acha.2020.01.004
|
| [19] |
J. Shen, Y. W. Wang, J. L. Xia, Fast structured Jacobi–Jacobi transforms, Math. Comput., 88 (2019), 1743–1772. https://doi.org/10.1090/mcom/3377 doi: 10.1090/mcom/3377
|
| [20] |
J. L. Aurentz, R. M. Slevinsky, On symmetrizing the ultraspherical spectral method, J. Comput. Phys., 420 (2020), 109383. https://doi.org/10.1016/j.jcp.2020.109383 doi: 10.1016/j.jcp.2020.109383
|
| [21] |
D. Fortunato, A. Townsend, H. Wilber, The ultraspherical spectral element method, J. Comput. Phys., 423 (2021), 110087. https://doi.org/10.1016/j.jcp.2020.110087 doi: 10.1016/j.jcp.2020.110087
|
| [22] | W. Cheng, Z. Xu, Understanding the ultraspherical spectral method: An error analysis, arXiv Preprint, 2024. https://doi.org/10.48550/arXiv.2405.11847 |
| [23] | R. Koekoek, P. A. Lesky, R. F. Swarttouw, Hypergeometric orthogonal polynomials and their $q$-analogues, Berlin: Springer, 2010. https://doi.org/10.1007/978-3-642-05014-5 |
| [24] |
H. Pijeira-Cabrera, J. Quintero-Roba, J. Toribio-Milane, Differential properties of Jacobi-Sobolev polynomials and electrostatic interpretation, Mathematics, 11 (2023), 3420. https://doi.org/10.3390/math11153420 doi: 10.3390/math11153420
|
| [25] | G. T. Núñez, Propiedades y aplicaciones de los polinomios clásicos en la resolución de problemas mecanocuánticos (átomo de hidrógeno y potencial de Rosen–Morse), M.Sc. thesis, Universidad Autónoma de Santo Domingo, 2017. |
| [26] | J. M. Peña-Reyes, Integral espectral del operador de Hamilton de la ecuación de Schrödinger y aplicaciones en mecánica cuántica, M.Sc. thesis, Universidad Autónoma de Santo Domingo, 2026. https://doi.org/10.13140/RG.2.2.13580.58240 |