We investigated a nonlinear singularly perturbed elliptic reaction-diffusion coupled system having non-smooth data networked by a $ k $-star graph. We considered all possible boundary conditions at the free boundary located at the tail of the edge and imposed the continuity condition with Kirchhoff's junction law at the junction point of the $ k $-star graph to obtain a continuous solution for this coupled system. First, we showed the existence and uniqueness of the solution using the variational formulation approach. Then, we reformulated it into a minimization problem over a function space to conclude the uniqueness of the solution. For the approximation of the continuous problem, note that the upwind scheme for the flux condition at the free boundary leads to a parameter uniform first-order approximation. To obtain a higher-order uniform accuracy, we utilized a three-point scheme for first-order derivatives and a five-point approximation at the point of discontinuity. These approximations typically did not yield an M-matrix or strict diagonally dominant structure of the stiffness matrix. Hence, we provided a suitable transformation that could lead to a sufficient condition for preserving the strict diagonally dominant structure of the stiffness matrix. We performed a comprehensive convergence analysis to demonstrate the almost second-order uniform accuracy on each edge of the $ k $-star graph. Numerical experiments highly validate the theory on the $ k $-star graph.
Citation: Dilip Sarkar, Shridhar Kumar, Pratibhamoy Das, Higinio Ramos. Higher-order convergence analysis for interior and boundary layers in a semi-linear reaction-diffusion system networked by a $ k $-star graph with non-smooth source terms[J]. Networks and Heterogeneous Media, 2024, 19(3): 1085-1115. doi: 10.3934/nhm.2024048
We investigated a nonlinear singularly perturbed elliptic reaction-diffusion coupled system having non-smooth data networked by a $ k $-star graph. We considered all possible boundary conditions at the free boundary located at the tail of the edge and imposed the continuity condition with Kirchhoff's junction law at the junction point of the $ k $-star graph to obtain a continuous solution for this coupled system. First, we showed the existence and uniqueness of the solution using the variational formulation approach. Then, we reformulated it into a minimization problem over a function space to conclude the uniqueness of the solution. For the approximation of the continuous problem, note that the upwind scheme for the flux condition at the free boundary leads to a parameter uniform first-order approximation. To obtain a higher-order uniform accuracy, we utilized a three-point scheme for first-order derivatives and a five-point approximation at the point of discontinuity. These approximations typically did not yield an M-matrix or strict diagonally dominant structure of the stiffness matrix. Hence, we provided a suitable transformation that could lead to a sufficient condition for preserving the strict diagonally dominant structure of the stiffness matrix. We performed a comprehensive convergence analysis to demonstrate the almost second-order uniform accuracy on each edge of the $ k $-star graph. Numerical experiments highly validate the theory on the $ k $-star graph.
[1] | Y. V. Pokornyi, A. V. Borovskikh, Differential equations on networks (geometric graphs), J. Math. Sci., 119 (2004), 691–718. https://doi.org/10.1023/B:JOTH.0000012752.77290.fa doi: 10.1023/B:JOTH.0000012752.77290.fa |
[2] | B. S. Pavlov, M. D. Faddeev, Model of free electrons and the scattering problem, Theor. Math. Phys., 55 (1983), 485–492. https://doi.org/10.1007/BF01015809 doi: 10.1007/BF01015809 |
[3] | T. Nagatani, Traffic flow on star graph: Nonlinear diffusion, Physica A, 561 (2021), 125251. https://doi.org/10.1016/j.physa.2020.125251 doi: 10.1016/j.physa.2020.125251 |
[4] | D. B. West, Introduction to graph theory, Prentice Hall, Inc., Upper Saddle River, NJ, 1996. |
[5] | W. C. Connor, J. Wengong, R. Luke, F. J. Timothy, S. J. Tommi, H. G. William, et al., A graph-convolutional neural network model for the prediction of chemical reactivity, Chem. Sci., 10 (2019), 370–377. https://doi.org/10.1039/C8SC04228D doi: 10.1039/C8SC04228D |
[6] | J. D. Murray, Mathematical biology: Ⅱ: Spatial models and biomedical applications, Interdiscip. Appl. Math., 2003. |
[7] | L. O. Müller, G. Leugering, P. J. Blanco, Consistent treatment of viscoelastic effects at junctions in one-dimensional blood flow models, J. Comput. Phys., 314 (2016), 167–193. https://doi.org/10.1016/j.jcp.2016.03.012 doi: 10.1016/j.jcp.2016.03.012 |
[8] | I. Rodriguez-Iturbe, R. Muneepeerakul, E. Bertuzzo, S. A. Levin, A. Rinaldo, River networks as ecological corridors: A complex systems perspective for integrating hydrologic, geomorphologic, and ecologic dynamics, Water Resour. Res., 45 (2009), 1–22. https://doi.org/10.1029/2008WR007124 doi: 10.1029/2008WR007124 |
[9] | J. V. Below, A. J. Lubary, Instability of stationary solutions of reaction-diffusion equations on graphs, Results Math., 68 (2015), 171–201. https://doi.org/10.1007/s00025-014-0429-8 doi: 10.1007/s00025-014-0429-8 |
[10] | S. Iwasaki, S. Jimbo, Y. Morita, Standing waves of reaction-diffusion equations on an unbounded graph with two vertices, SIAM J. Appl. Math., 82 (2022), 1733–1763. https://doi.org/10.1137/21M1454572 doi: 10.1137/21M1454572 |
[11] | H. M. Srivastava, A. K. Nain, R. K. Vats, P. Das, A theoretical study of the fractional-order $p$-Laplacian nonlinear Hadamard type turbulent flow models having the Ulam–Hyers stability, Rev. Real Acad. Cienc. Exactas Fis. Nat. Ser. A-Mat., 117 (2023), 1–19. https://doi.org/10.1007/s13398-023-01488-6 doi: 10.1007/s13398-023-01488-6 |
[12] | V. Mehandiratta, M. Mehra, G. Leugering, Existence and uniqueness results for a nonlinear Caputo fractional boundary value problem on a star graph, J. Math. Anal. Appl., 477 (2019), 1243–1264. https://doi.org/10.1016/j.jmaa.2019.05.011 doi: 10.1016/j.jmaa.2019.05.011 |
[13] | D. G. Gordeziani, M. Kupreishvili, H. V. Meladze, T. D. Davitashvili, On the solution of boundary value problem for differential equations given in graphs, Appl. Math. Inform. Mech., 13 (2008), 80–91. |
[14] | G. M. Gie, M. Hamouda, C. Y. Jung, R. M. Temam, Singular Perturbations and Boundary Layers, Springer International Publishing, 2018. https://doi.org/10.1007/978-3-030-00638-9 doi: 10.1007/978-3-030-00638-9 |
[15] | P. Das, S. Rana, J. Vigo-Aguiar, Higher-order accurate approximations on equidistributed meshes for boundary layer originated mixed type reaction-diffusion systems with multiple scale nature, Appl. Numer. Math., 148 (2020), 79–97. https://doi.org/10.1016/j.apnum.2019.08.028 doi: 10.1016/j.apnum.2019.08.028 |
[16] | P. Das, Comparison of a priori and a posteriori meshes for singularly perturbed nonlinear parameterized problems, J. Comput. Appl. Math., 290 (2015), 16–25. https://doi.org/10.1016/j.cam.2015.04.034 doi: 10.1016/j.cam.2015.04.034 |
[17] | P. Das, An a posteriori based convergence analysis for a nonlinear singularly perturbed system of delay differential equations on an adaptive mesh, Numerical Algorithms, 81 (2019), 465–487. https://doi.org/10.1007/s11075-018-0557-4 doi: 10.1007/s11075-018-0557-4 |
[18] | P. Das, A higher order difference method for singularly perturbed parabolic partial differential equations, J. Differ. Equations Appl., 24 (2018), 452–477. https://www.tandfonline.com/doi/full/10.1080/10236198.2017.1420792 doi: 10.1080/10236198.2017.1420792 |
[19] | D. Shakti, J. Mohapatra, P. Das, J. Vigo-Aguiar, A moving mesh refinement based optimal accurate uniformly convergent computational method for a parabolic system of boundary layer originated reaction-diffusion problems with arbitrary small diffusion terms, J. Comput. Appl. Math., 404 (2022), 113167. https://doi.org/10.1016/j.cam.2020.113167 doi: 10.1016/j.cam.2020.113167 |
[20] | S. Kumar, P. Das, K. Kumar, Adaptive mesh-based efficient approximations for Darcy scale precipitation–dissolution models in porous media, Int. J. Numer. Methods Fluids, 96 (2024), 1415–1444. https://doi.org/10.1002/fld.5294 doi: 10.1002/fld.5294 |
[21] | S. Saini, P. Das, S. Kumar, Parameter uniform higher order numerical treatment for singularly perturbed Robin type parabolic reaction-diffusion multiple scale problems with large delay in time, Appl. Numer. Math., 196 (2024), 1–21. https://doi.org/10.1016/j.apnum.2023.10.003 doi: 10.1016/j.apnum.2023.10.003 |
[22] | S. Sain, P. Das, S. Kumar, Computational cost reduction for coupled system of multiple scale reaction-diffusion problems with mixed type boundary conditions having boundary layers, Rev. Real Acad. Cienc. Exactas Fís. Nat. Ser. A Mat., 117 (2023), 1–27. https://doi.org/10.1007/s13398-023-01397-8 doi: 10.1007/s13398-023-01397-8 |
[23] | S. Kumar, P. Das, A uniformly convergent analysis for multiple scale parabolic singularly perturbed convection-diffusion coupled systems: Optimal accuracy with less computational time, Appl. Numer. Math., 207 (2025), 534–557. https://doi.org/10.1016/j.apnum.2024.09.020 doi: 10.1016/j.apnum.2024.09.020 |
[24] | B. P. Andreianov, G. M. Coclite, C. Donadello, Well-posedness for vanishing viscosity solutions of scalar conservation laws on a network, Discrete Contin. Dyn. Syst., 37 (2017), 5913–5942. http://dx.doi.org/10.3934/dcds.2017257 doi: 10.3934/dcds.2017257 |
[25] | M. Musch, U. S. Fjordholm, N. H. Risebro, Well-posedness theory for nonlinear scalar conservation laws on networks, Networks Heterogen. Media, 17 (2022), 101–128. https://doi.org/10.3934/nhm.2021025 doi: 10.3934/nhm.2021025 |
[26] | G. M. Coclite, C. Donadello, Vanishing viscosity on a star-shaped graph under general transmission conditions at the node, Networks Heterogen. Media, 15 (2020), 197–213. https://doi.org/10.3934/nhm.2020009 doi: 10.3934/nhm.2020009 |
[27] | J. D. Towers, An explicit finite volume algorithm for vanishing viscosity solutions on a network, Networks Heterogen. Media, 17 (2022), 1–13. https://doi.org/10.3934/nhm.2021021 doi: 10.3934/nhm.2021021 |
[28] | S. F. Pellegrino, On the implementation of a finite volumes scheme with monotone transmission conditions for scalar conservation laws on a star-shaped network, Appl. Numer. Math., 155 (2020), 181–191. https://doi.org/10.1016/j.apnum.2019.09.011 doi: 10.1016/j.apnum.2019.09.011 |
[29] | F. R. Guarguaglini, R. Natalini, Vanishing viscosity approximation for linear transport equations on finite star-shaped networks, J. Evol. Equations, 21 (2021), 2413–2447. https://doi.org/10.1007/s00028-021-00688-0 doi: 10.1007/s00028-021-00688-0 |
[30] | H. Egger, N. Philippi, On the transport limit of singularly perturbed convection-diffusion problems on networks, Math. Methods Appl. Sci., 44 (2021), 5005–5020. https://doi.org/10.1002/mma.7084 doi: 10.1002/mma.7084 |
[31] | H. Egger, N. Philippi, A hybrid discontinuous Galerkin method for transport equations on networks, in Finite volumes for complex applications IX, Bergen, Norway, Springer, 323 (2020), 487–495. https://doi.org/10.1007/978-3-030-43651-3_45 |
[32] | H. Egger, N. Philippi, A hybrid-DG method for singularly perturbed convection-diffusion equations on pipe networks, ESAIM Math. Model. Numer. Anal., 57 (2023), 2077–2095. https://doi.org/10.1051/m2an/2023044 doi: 10.1051/m2an/2023044 |
[33] | V. Kumar, G. Leugering, Singularly perturbed reaction-diffusion problems on a k-star graph, Math. Methods Appl. Sci., 44 (2021), 14874–14891. https://doi.org/10.1002/mma.7749 doi: 10.1002/mma.7749 |
[34] | P. A. Farrell, J. J. H. Miller, E. O'Riordan, G. I. Shishkin, Singularly perturbed differential equations with discontinuous source terms, in Analytical and Numerical Methods for Convection-Dominated and Singularly Perturbed Problems (eds. J.J.H. Miller, G.I. Shishkin and L. Vulkov), Nova Science Publishers, New York, (2000), 23–32. |
[35] | Z. Cen, A hybrid difference scheme for a singularly perturbed convection-diffusion problem with discontinuous convection coefficient, Appl. Math. Comput., 169 (2005), 689–699. https://doi.org/10.1016/j.amc.2004.08.051 doi: 10.1016/j.amc.2004.08.051 |
[36] | S. Kumar, S. Kumar, P. Das, Second-order a priori and a posteriori error estimations for integral boundary value problems of nonlinear singularly perturbed parameterized form, Numerical Algorithms, 2024. https://doi.org/10.1007/s11075-024-01918-5 |
[37] | S. Santra, J. Mohapatra, P. Das, D. Choudhari, Higher-order approximations for fractional order integro-parabolic partial differential equations on an adaptive mesh with error analysis, Comput. Math. Appl., 150 (2023), 87–101. https://doi.org/10.1016/j.camwa.2023.09.008 doi: 10.1016/j.camwa.2023.09.008 |
[38] | V. Kumar, G. Leugering, Convection dominated singularly perturbed problems on a metric graph, J. Comput. Appl. Math. 425 (2023), 115062. https://doi.org/10.1016/j.cam.2023.115062 doi: 10.1016/j.cam.2023.115062 |
[39] | H. Zhu, Z. Li, Z. Yang, Analysis and computation for a class of semilinear elliptic boundary value problems, Comput. Math. Appl., 64 (2012), 2735–2743. https://doi.org/10.1016/j.camwa.2012.08.004 doi: 10.1016/j.camwa.2012.08.004 |
[40] | J. J. Nieto, J. M. Uzal, Nonlinear second-order impulsive differential problems with dependence on the derivative via variational structure, J. Fixed Point Theory Appl., 22 (2020), 1–19. https://doi.org/10.1007/s11784-019-0754-3 doi: 10.1007/s11784-019-0754-3 |
[41] | K. Atkinson, W. Han, Theoretical Numerical Analysis: A Functional Analysis Framework, Springer, New York, 2009. |
[42] | G. V. Alekseev, R. V. Brizitskii, Z. Y. Saritskaya, Stability estimates of solutions to extremal problems for a nonlinear convection-diffusion-reaction equation, J. Appl. Ind. Math., 10 (2016), 155–167. https://doi.org/10.17377/sibjim.2016.19.201 doi: 10.17377/sibjim.2016.19.201 |
[43] | M. Manikandan, R. Ishwariya, Robust computational technique for a class of singularly perturbed nonlinear differential equations with Robin boundary conditions, J. Math. Model., 11 (2023), 411–423. https://doi.org/10.22124/jmm.2023.23515.2100 doi: 10.22124/jmm.2023.23515.2100 |
[44] | R. Shiromani, V. Shanthi, P. Das, A higher order hybrid-numerical approximation for a class of singularly perturbed two-dimensional convection-diffusion elliptic problem with non-smooth convection and source terms, Comput. Math. Appl., 142 (2023), 9–30. https://doi.org/10.1016/j.camwa.2023.04.004 doi: 10.1016/j.camwa.2023.04.004 |
[45] | M. Chandru, T. Prabha, V. Shanthi, A hybrid difference scheme for a second-order singularly perturbed reaction-diffusion problem with non-smooth data, Int. J. Appl. Comput. Math., 1 (2015), 87–100. https://doi.org/10.1007/s40819-014-0004-8 doi: 10.1007/s40819-014-0004-8 |
[46] | J. J. H. Miller, E. O'Riordan, G. I. Shishkin, Fitted Numerical Methods for Singular Perturbation Problems: Error Estimates in the Maximum Norm for Linear Problems in One and Two Dimensions, World Scientific, 1996. |
[47] | P. A. Farrell, A. F. Hegarty, J. J. H. Miller, E. O'Riordan, G. I. Shishkin, Robust Computational Techniques for Boundary Layers, Chapman and Hall, CRC Press, Boca Raton, USA, 2000. |
[48] | N. Kopteva, M. Stynes, Numerical analysis of a singularly perturbed nonlinear reaction-diffusion problem with multiple solutions, Appl. Numer. Math., 51 (2004), 273–288. https://doi.org/10.1016/j.apnum.2004.07.001 doi: 10.1016/j.apnum.2004.07.001 |
[49] | S. Kumar, R. Ishwariya, P. Das, Impact of mixed boundary conditions and non-smooth data on layer originated non-premixed combustion problems: Higher order convergence analysis, Stud. Appl. Math., (2024), e12673. https://doi.org/10.1111/sapm.12763 doi: 10.1111/sapm.12763 |