Research Article | | Peer-Reviewed

Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems

Received: 8 June 2025     Accepted: 30 June 2025     Published: 15 August 2025
Views:       Downloads:
Abstract

Several recently published studies regarding flow problems propose schemes of high order of accuracy designed as evolution of traditional methods. A drawback common to these new schemes is the necessity to adopt uniform mesh refinement for solving sharp problems, by increasing the computational cost. Even the so called essentially non-oscillatory and weight essentially non-oscillatory methods suffer of the same drawback and are not suitable to cope with h-adaptive methods due to their definition on finite volumes necessarily of equal diameter. Therefore, in order to overcome the above drawback, the formulation of dynamically locally self h-adaptive processes is designed to achieve the dual purpose to increase the accuracy and to keep as small as possible the number of finite volumes. To define a locally h-adaptive finite volume (FV) scheme need two simple but important tools, namely a particular FV named Bridge FV positioned between two adjacent subdomains and the definition of suitable profiles approximating the fluxes on the FV faces. In this article a new FV method for the numerical solution of convective-diffusive 1D problems is developed. It is conservative, second order in time and space for equal FV, and allows the partitioning of the domain by equal or unequal finite volumes, thus dynamically locally self h-adaptive. The definition of the monotonic profiles is accomplished by means of cubic weighted ν-splines and Taylor expansions. The profile analysis respect to the numerical properties is conducted in the normalized plane with the velocity varying in time and space and gives the flux value on the FV faces. Moreover the flux is assigned by Upwind or by second order back-ward Characteristics if the estimated flux is outside of the unit square or the transformation into the normalized plane is not possible, respectively. The initial-boundary stability and convergence properties of the new method are examined in detail, also in presence of h-adaptivity. In addition, a generalization of the new scheme to 2D and 3D problems is presented. Finally, some numerical test are carried out to verify the properties of the new method, including two CFD problems.

Published in Pure and Applied Mathematics Journal (Volume 14, Issue 4)
DOI 10.11648/j.pamj.20251404.11
Page(s) 69-92
Creative Commons

This is an Open Access article, distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium or format, provided the original work is properly cited.

Copyright

Copyright © The Author(s), 2025. Published by Science Publishing Group

Keywords

FVM for Convective-diffusive Problems, High Order Monotonic Schemes, Convergence for Initial-boundary Values Problems, Dynamically Local Self h-Adaptivity

References
[1] Franca, L. P., Frey, S. L., Hughes, T. J. R.: Stabilized finite element methods: I. Application to the advective-diffusive model. Comput. Methods Appl. Mech. Engrg. 95, 253-276 (1992).
[2] Hughes, T. J. R., Franca, L. P., Hulbert, G. M.: A new finite element formulation for computational fluid dynamics: VIII. The Galerkin/least-squares method for advective-diffusive equations. Comput. Methods Appl. Mech. Engrg. 73, 173-189 (1989).
[3] Kengni Jotsa, A. C., Pennati, V. A.: Solution of the 2D Navier-Stokes equations by a new FE fractional step method. IMACS/ISGG WORKSHOP, 141-150, Roma, Italy, 2011.
[4] Kengni Jotsa, A. C.: Solution of the 2D Navier-Stokes equations by a new FE fractional step method. Doctoral dissertation, Università degli Studi dell’Insubria, Como, Italy. (2012). Available:
[5] Kengni Jotsa, A. C., Pennati, V. A.: A cost-effective FE method for 2D Navier–Stokes equations. Engineering Applications of Computational Fluid Mechanics. 9, 66-83 (2015).
[6] Khelifa, A., Robert, J.-L., Ouellet, Y.: A Douglas-Wang finite element approach for transient advection-diffusion problems. Comput. Methods Appl. Mech. Engrg. 10, 113-129 (1993).
[7] Siegel, P., Mosé, R., Ackerer, P., Jaffre, J.: Solution of the Advection-Diffusion Equation Using a Combination of Discontinuous and Mixed Finite Elements. Int. J. Numer. Meth. Fluids. 24, 595-613 (1997).
[8] Harten, A., Engquist, B., Osher, S., Chakravarthy, S. R.: Uniformly high order accurate essentially non-oscillatory schemes, III. J. Comput. Phys. 71, 231-303 (1987).
[9] Quarteroni, A.: Modellistica numerica per problemi differenziali. Springer-Verlag, Milano (2008).
[10] Leonard, B. P.: A stable and accurate convective modelling procedure based on quadratic upstream interpolation. Comput. Methods Appl. Mech. Engrg. 19, 59-98 (1979).
[11] Leonard, B. P.: Simple high-accuracy resolution program for convective modelling of discontinuities. Int. J. Numer. Meth. Fluids. 8, 1291-1318 (1988).
[12] De Biase, L., Feraudi, F., Pennati, V. A.: A finite volume method for the solution of convection-diffusion 2D problems by a quadratic profile with smoothing. Int. J. Num. Meth. Heat Fluid Flow. 6, 3-24 (1996).
[13] Deponti, A., Pennati, V. A., De Biase, L.: A fully 3D finite volume method for incompressible Navier–Stokes equations. Int. J. Num. Meth. Fluids. 52, 617-638 (2006).
[14] Deponti, A., Pennati, V. A., De Biase, L., Maggi, V., Berta, F.: A new fully three-dimensional numerical model for ice dynamics. Journal of Glaciology. 52, 365-376 (2006).
[15] Feraudi, F., Pennati, V. A.: Trasporto del Calore in Fluidi Incomprimibili: un Nuovo Approccio Numerico per Problemi Bidimensionali Non Stazionari. L’Energia Elettrica. 74, 238-251 (1993).
[16] Leonard, B. P.: The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection. Comput. Methods Appl. Mech. Engrg. 88, 17-74 (1991).
[17] Leonard, B. P., Mokhtari, S.: Beyond first-order upwinding: The ultra-sharp alternative for non-oscillatory steady-state simulation of convection. Int. J. Numer. Meth. Eng., 30, 729-766 (1990).
[18] Babu, G., Bansal, K.: A high order robust numerical scheme for singularly perturbed delay parabolic convection diffusion problems. Journal of Applied Mathematics and Computing. 68, 363-389 (2022).
[19] Du, Y., Wang, Y., Yuan, L.: A high-order modified finite-volume method on Cartesian grids for nonlinear convection–diffusion problems. Comp. Appl. Math. 39, 214 (2020).
[20] Xu, M.: A high-order finite volume scheme for unsteady convection-dominated convection–diffusion equations. Numerical Heat Transfer, Part B: Fundamentals. 76, 253-272 (2019).
[21] Zhu, X., Rui, H.: High-order compact difference scheme of 1D nonlinear degenerate convection–reaction–diffusion equation with adaptive algorithm. Numerical Heat Transfer, Part B: Fundamentals. 75, 43-66 (2019).
[22] Huang, L., Jiang, Z., Lou, S., Zhang, X., Yan, C.: Simple and robust h-adaptive shock-capturing method for flux reconstruction framework. Chinese Journal of Aeronautics. 36, 348-365 (2023).
[23] Nishikawa, H.: The QUICK Scheme is a Third-Order Finite-Volume Scheme with Point-Valued Numerical Solutions. arXiv preprint arXiv:2006.15143 [math. NA]. (2020).
[24] Oh, J.: Adaptive Mesh-free Methods for Partial Differential Equations. Doctoral dissertation. Department of Mathematics of the University of Southern Mississippi, USA. (2018). Available:
[25] Sabbagh-Yazdi, S.-R., Amiri, T., Asil Gharebaghi, S.: A Proposed Damping Coefficient of Quick Adaptive Galerkin Finite Volume Solver for Elasticity Problems. Journal of the Serbian Society for Computational Mechanics. 13, 56-79 (2019).
[26] Zheng, Q., Liu, Z.: Uniform convergence analysis of a new adaptive upwind finite difference method for singularly perturbed convection reaction - diffusion boundary value problem. Journal of Applied and Computing. 70, 601-618 (2024).
[27] Kumar, S., Das, P.: A uniformly convergent analysis for multiple scale parabolic singularly perturbed convection diffusion coupled systems: Optimal accuracy with less computational time. Appl. Numer. Math. 207, 534-557 (2025).
[28] Ayalew, M., Aychluh, M., Suthar, D. L., Purohit, S. D.: Quadratic upwind differencing scheme in the finite volume method for solving the convection-diffusion equation. Mathematical and Computer Modelling of Dynamical Systems. 29, 265-285 (2023).
[29] Peng, G.: A positivity-preserving finite volume scheme for convection–diffusion equation on general meshes. International Journal of Computer Mathematics. 99, 355-369 (2022).
[30] Kumar, N., Thije Boonkkamp, J. H. M., Koren, B.: Approximation of the convective flux in the incompressible Navier–Stokes equations using local boundary-value problems. Journal of Computational and Applied Mathematics. 340, 523-536 (2018).
[31] Sukhinov, A., Chistyakov, A., Kuznetsova, I., Belova, Y., Rahimbaeva, E.: Development and Research of a Modified Upwind Leapfrog Scheme for Solving Transport Problems. Mathematics. 10, 3564 (2022).
[32] Kurganov, A.: Finite-volume schemes for shallow-water equations. Cambridge University Press. Acta Numerica. 27, 289-351 (2018).
[33] Neelan, A. A. G., Nair, M. T., Bürger, R.: Three-level order-adaptive weighted essentially non-oscillatory schemes. Results in Applied Mathematics. 12, 100217 (2021).
[34] Kossaczká, T., Ehrhardt, M., Günther, M.: Enhanced fifth order WENO shock-capturing schemes with deep learning. Results in Applied Mathematics. 12, 100201 (2021).
[35] Jiang, G.-S., Shu, C.-W.: Efficient Implementation of Weighted ENO Schemes. J. Comput. Phys. 126, 202-228 (1996).
[36] Xu, M.: A finite volume scheme for unsteady linear and nonlinear convection-diffusion-reaction problems. International Communications in Heat and Mass Transfer. 139, 106417 (2022).
[37] Pennati, V. A., Marelli, M., De Biase, L. M.: FAST - A generalized FV solution of convection-diffusion problems by means of monotonic v-splines profiles. Int. J. Num. Meth. Heat Fluid Flow. 6, 3-30 (1996).
[38] Roache, P. J.: Computational fluid dynamics. Hermosa Publishers, Albuquerque, NM (1972).
[39] Godunov, S. K.: A difference scheme for numerical solution of discontinuous solution of hydrodynamic equations. Mathematics of the USSR-Sbornik. 47, 271-306 (1959). Translated US Joint Publ. Ross. Service, JPRS7226, 1969.
[40] Dullemond, C. P.: Lecture Numerical Fluid Dynamics. Chapter 4. University of Heidelberg: Heidelberg University Press (2011).
[41] Strikwerda, J. C.: Finite difference schemes and partial differential equations. Wadsworth & Brooks/Cole Advanced Books & Software, Pacific Grove, California, USA (1989).
[42] Jury, E. I.: On the roots of a real polynomial inside the unit circle and a stability criterion for linear discrete systems. IFAC Proceedings Volumes. 1, 142-153 (1963).
[43] Marden, M.: The Geometry of the Zeros of a Polynomial in a Complex Variable. American Mathematical Society, New York (1949).
[44] Reali, M., Dassie, G., Pennati, V. A.: Direct general finite difference techniques for elliptic problems defined on bounded or unbounded 2-D domains. Lewis, R. W. et al. (ed.). Numerical Methods for Transient and Coupled Problems, pp. 43-58. Wiley. (1987).
[45] Corti, S., De Biase, L., Pennati, V. A.: A multistep three dimensional generalized finite difference technique. Sixth International Conference on Numerical methods in Thermal problems. Pineridge Press, Swansea, U. K. (1989).
[46] Kelly, D. W., Mills, R. J., Reizes, J. A., Miller, A. D.: A posteriori error estimates in finite difference techniques. J. Comput. Phys. 74, 214-232 (1988).
[47] Pennati, V., Corti, S.: Generalized finite-differences solution of 3D elliptical problems involving Neumann boundary conditions. Commun. Numer. Meth. Engng. 10, 43-58 (1994).
[48] Thompson, J. F., Soni, B. K., Weatherill, N. P.: Handbook of grid generation. CRC press (1998).
[49] Deponti, A.: Mass and thermal flows in Alpine glaciers. Application to Lys Glacier (Monte Rosa, Italian Alps). Doctoral dissertation, University of Milano–Bicocca. Italy (2003).
[50] Cordero, E., De Biase, L., Pennati, V. A.: A new finite volume method for the solution of convection–diffusion equations: analysis of stability and convergence. Commun. Numer. Meth. Engng. 13, 923-940 (1997).
[51] Cunge, J. A., Holly, F. M., Verwey, A.: Practical Aspects of Computational River Hydraulics. Pitman, London (1980).
[52] Ambrosi, D., Corti, S., Pennati, V., Saleri, F.: Numerical Simulation of Unsteady Flow at Po River Delta. Journal of Hydraulic Engineering. 122, 735-743 (1996).
[53] Kengni Jotsa, A. C., Pennati, V. A., Di Guardo, A., Morselli, M.: Shallow Water 1D Model for Pollution River Study. Pure and Applied Mathematics Journal. 6, 76-88 (2017).
[54] Vignoli, G., Titarev, V. A., Toro, E. F.: ADER schemes for the shallow water equations in channel with irregular bottom elevation. J. Comput. Phys. 227, 2463-2480 (2008).
Cite This Article
  • APA Style

    Pennati, V. A., Jotsa, A. C. K., Tagoudjeu, J. (2025). Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems. Pure and Applied Mathematics Journal, 14(4), 69-92. https://doi.org/10.11648/j.pamj.20251404.11

    Copy | Download

    ACS Style

    Pennati, V. A.; Jotsa, A. C. K.; Tagoudjeu, J. Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems. Pure Appl. Math. J. 2025, 14(4), 69-92. doi: 10.11648/j.pamj.20251404.11

    Copy | Download

    AMA Style

    Pennati VA, Jotsa ACK, Tagoudjeu J. Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems. Pure Appl Math J. 2025;14(4):69-92. doi: 10.11648/j.pamj.20251404.11

    Copy | Download

  • @article{10.11648/j.pamj.20251404.11,
      author = {Vincenzo Angelo Pennati and Antoine Celestin Kengni Jotsa and Jacques Tagoudjeu},
      title = {Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems
    },
      journal = {Pure and Applied Mathematics Journal},
      volume = {14},
      number = {4},
      pages = {69-92},
      doi = {10.11648/j.pamj.20251404.11},
      url = {https://doi.org/10.11648/j.pamj.20251404.11},
      eprint = {https://article.sciencepublishinggroup.com/pdf/10.11648.j.pamj.20251404.11},
      abstract = {Several recently published studies regarding flow problems propose schemes of high order of accuracy designed as evolution of traditional methods. A drawback common to these new schemes is the necessity to adopt uniform mesh refinement for solving sharp problems, by increasing the computational cost. Even the so called essentially non-oscillatory and weight essentially non-oscillatory methods suffer of the same drawback and are not suitable to cope with h-adaptive methods due to their definition on finite volumes necessarily of equal diameter. Therefore, in order to overcome the above drawback, the formulation of dynamically locally self h-adaptive processes is designed to achieve the dual purpose to increase the accuracy and to keep as small as possible the number of finite volumes. To define a locally h-adaptive finite volume (FV) scheme need two simple but important tools, namely a particular FV named Bridge FV positioned between two adjacent subdomains and the definition of suitable profiles approximating the fluxes on the FV faces. In this article a new FV method for the numerical solution of convective-diffusive 1D problems is developed. It is conservative, second order in time and space for equal FV, and allows the partitioning of the domain by equal or unequal finite volumes, thus dynamically locally self h-adaptive. The definition of the monotonic profiles is accomplished by means of cubic weighted ν-splines and Taylor expansions. The profile analysis respect to the numerical properties is conducted in the normalized plane with the velocity varying in time and space and gives the flux value on the FV faces. Moreover the flux is assigned by Upwind or by second order back-ward Characteristics if the estimated flux is outside of the unit square or the transformation into the normalized plane is not possible, respectively. The initial-boundary stability and convergence properties of the new method are examined in detail, also in presence of h-adaptivity. In addition, a generalization of the new scheme to 2D and 3D problems is presented. Finally, some numerical test are carried out to verify the properties of the new method, including two CFD problems.},
     year = {2025}
    }
    

    Copy | Download

  • TY  - JOUR
    T1  - Fastest - A New High Order FV Method Dynamically Locally Self h-Adaptive for Convective-diffusive Problems
    
    AU  - Vincenzo Angelo Pennati
    AU  - Antoine Celestin Kengni Jotsa
    AU  - Jacques Tagoudjeu
    Y1  - 2025/08/15
    PY  - 2025
    N1  - https://doi.org/10.11648/j.pamj.20251404.11
    DO  - 10.11648/j.pamj.20251404.11
    T2  - Pure and Applied Mathematics Journal
    JF  - Pure and Applied Mathematics Journal
    JO  - Pure and Applied Mathematics Journal
    SP  - 69
    EP  - 92
    PB  - Science Publishing Group
    SN  - 2326-9812
    UR  - https://doi.org/10.11648/j.pamj.20251404.11
    AB  - Several recently published studies regarding flow problems propose schemes of high order of accuracy designed as evolution of traditional methods. A drawback common to these new schemes is the necessity to adopt uniform mesh refinement for solving sharp problems, by increasing the computational cost. Even the so called essentially non-oscillatory and weight essentially non-oscillatory methods suffer of the same drawback and are not suitable to cope with h-adaptive methods due to their definition on finite volumes necessarily of equal diameter. Therefore, in order to overcome the above drawback, the formulation of dynamically locally self h-adaptive processes is designed to achieve the dual purpose to increase the accuracy and to keep as small as possible the number of finite volumes. To define a locally h-adaptive finite volume (FV) scheme need two simple but important tools, namely a particular FV named Bridge FV positioned between two adjacent subdomains and the definition of suitable profiles approximating the fluxes on the FV faces. In this article a new FV method for the numerical solution of convective-diffusive 1D problems is developed. It is conservative, second order in time and space for equal FV, and allows the partitioning of the domain by equal or unequal finite volumes, thus dynamically locally self h-adaptive. The definition of the monotonic profiles is accomplished by means of cubic weighted ν-splines and Taylor expansions. The profile analysis respect to the numerical properties is conducted in the normalized plane with the velocity varying in time and space and gives the flux value on the FV faces. Moreover the flux is assigned by Upwind or by second order back-ward Characteristics if the estimated flux is outside of the unit square or the transformation into the normalized plane is not possible, respectively. The initial-boundary stability and convergence properties of the new method are examined in detail, also in presence of h-adaptivity. In addition, a generalization of the new scheme to 2D and 3D problems is presented. Finally, some numerical test are carried out to verify the properties of the new method, including two CFD problems.
    VL  - 14
    IS  - 4
    ER  - 

    Copy | Download

Author Information
  • Sections