User menu

Accès à distance ? S'identifier sur le proxy UCLouvain

Application of modified Patankar schemes to stiff biogeochemical models for the water column

Primary tabs

Burchard, Hans Deleersnijder, Eric [UCL] Meister, Andreas

In this paper, we apply recently developed positivity preserving and conservative Modified Patankar-type solvers for ordinary differential equations to a simple stiff biogeochemical model for the water column. The performance of this scheme is compared to schemes which are not unconditionally positivity preserving (the first-order Euler and the second- and fourth-order Runge–Kutta schemes) and to schemes which are not conservative (the first- and second-order Patankar schemes). The biogeochemical model chosen as a test ground is a standard nutrient–phytoplankton–zooplankton–detritus (NPZD) model, which has been made stiff by substantially decreasing the half saturation concentration for nutrients. For evaluating the stiffness of the biogeochemical model, so-called numerical time scales are defined which are obtained empirically by applying high-resolution numerical schemes. For all ODE solvers under investigation, the temporal error is analysed for a simple exponential decay law. The performance of all schemes is compared to a high-resolution high-order reference solution. As a result, the second-order modified Patankar–Runge–Kutta scheme gives a good agreement with the reference solution even for time steps 10 times longer than the shortest numerical time scale of the problem. Other schemes do either compute negative values for non-negative state variables (fully explicit schemes), violate conservation (the Patankar schemes) or show low accuracy (all first-order schemes).

Document type Article de périodique (Journal article) – Article de recherche
Access type Accès restreint
Publication date 2005
Language Anglais
Journal information "Ocean Dynamics : theoretical, computational oceanography and monitoring" - Vol. 55, no.3-4, p. 326-337 (2005)
Peer reviewed yes
Publisher Springer (Heidelberg)
issn 1616-7341
e-issn 1616-7228
Publication status Publié
Affiliation UCL - SST/IMMC/MEMA - Applied mechanics and mathematics
Links
  1. Abdulle A (2001) Fourth order Chebyshev methods with recurrence relations. SIAM J Sci Comp 23:2042–2055
  2. Abdulle A, Medovikov AA (2001) Second order Chebyshev methods based on orthogonal polynomials. Numer Math 90:1–18
  3. Bolding K, Burchard H, Pohlmann T, Stips A (2002) Turbulent mixing in the Northern North Sea: a numerical model study. Cont Shelf Res 22:2707–2724
  4. Burchard H, Deleersnijder E, Meister A (2003) A high-order conservative Patankar-type discretisation for stiff systems of production-destruction equations. Appl Numer Math 47:1–30
  5. Burchard H, Bolding K, Kühn W, Meister A, Neumann T,Umlauf L (2005) Description of a flexible and extendable physical-biogeochemical model system for the water column. J Mar Sys (in press)
  6. Deleersnijder E, Beckers J-M, Campin J-M, El Mohajir M, Fichefet T, Luyten P (1997) Some mathematical problems associated with the development and use of marine models. In: Diaz JI (ed) The mathematics of models for climatology and environment, vol 48. of NATO ASI Series, Springer, Berlin Heidelberg New York, pp 41–86
  7. Fasham M. J. R., Ducklow H. W., McKelvie S. M., A nitrogen-based model of plankton dynamics in the oceanic mixed layer, 10.1357/002224090784984678
  8. Fennel W, Neumann T (1996) The mesoscale variability of nutrients and plankton as seen in a coupled model. Dt Hydrogr Z 48:49–71
  9. Hairer E, Wanner G (2004) Solving ordinary differential equations II, Series in computational mathematics 14, 3rd edn. Springer, Berlin Heidelberg New York
  10. Hairer E, Nørsett S, Wanner G (2000) Solving ordinary differential equations I, Series in computational mathematics 8, 2nd edn. Springer, Berlin Heidelberg New York
  11. Harrison WG, Harris L, Irwin BD (1996) The kinetics of nitrogen utilization in the oceanic mixed layer: nitrate and ammonium interactions at nanomolar concentrations. Limnol Oceanogr 41:16–32
  12. Horváth Z (1998) Positivity of Runge–Kutta and diagonally split Runge–Kutta methods. Appl Numer Math 28:309–326
  13. van der Houwen PJ (1996) The development of Runge–Kutta methods for parabolic differential equations. Appl Num Math 20:261–273
  14. van der Houwen PJ, Sommeijer BP (1980) On the internal stability of explicit m-stage Runge–Kutta methods for large m-values. ZAMM 60:479–485
  15. Hundsdorfer W, Verwer JG (2003) Numerical solution of time-dependent advection-diffusion-reaction equations0, vol 33 of Series in computational mathematics. Springer, Berlin Heidelberg New York
  16. Kondo J (1975) Air-sea bulk transfer coefficients in diabatic conditions. Bound Layer Meteor 9:91–112
  17. Kühn W, Radach G (1997) A one-dimensional physical-biological model study of the pelagic nitrogen cycling during the spring bloom in the northern North Sea (FLEX’76). J Mar Res 55:687–734
  18. Lebedev V (2000) Explicit difference schemes for solving stiff problems with a complex or separable spectrum. Comp Math Math Phys 40:1729–1740
  19. Leonard BP (1991) The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection. Comput Meth Appl Mech Eng 88:17–74
  20. Medovikov AA (1998) High order explicit methods for stiff ordinary differential equations. BIT 38:372–390
  21. Meister A (1998) Comparison of different Krylov subspace methods embedded in an implicit finite volume scheme for the computation of viscous and inviscid flow fields on unstructured grids. J Comput Phys 140:311–345
  22. Oschlies A, Kähler P (2004) Biotic contribution to air-sea fluxes of CO2 and O2 and its relation to new production, export production, and net community production. Global Biogeochemical Cycles 18. GB1015, doi:10.1029/2003GB002094
  23. Patankar SV (1980) Numerical heat transfer and fluid flow. McGraw-Hill, New York
  24. Pietrzak J (1998) The use of TVD limiters for forward-in-time upstream-biased advection schemes in ocean modeling. Mon Weather Rev 126:812–830
  25. Popova EE, Ryabchenko VA, Fasham MJR (2000) Biological pump and vertical mixing in the southern ocean: their impact on atmospheric CO2. Global Biogeochem Cycles 14:477–498
  26. Robertson HH (1966) The solution of a set of reaction rate equations. In: Walsh J (ed) Numerical analysis, an introduction. Academic, New York, pp 178–182
  27. Sandu A, Verwer J, van Loon M, Carmichael G, Potra F, Dabdub D, Seinfeld J (1997) Benchmarking stiff ODE solvers for atmospheric chemistry problems I: Implicit versus explicit. Atmos Environ 31:3151–3166
  28. Soetaert K, Herman PMJ, Middelburg JJ (1996) A model of early diagenetic processes from the shelf to abyssal depths. Geochim Cosmochim Acta 60:1019–1040
  29. Verwer JG (1996) Explicit Runge–Kutta methods for parabolic partial differential equations. Appl Num Math 22:359–379
  30. Weber L, Völker C, Schartau M, Wolf-Gladrow DA (2005) Modelling the speciation and biochemistry of iron at the Bermuda Atlantik Time-series Study site, Global Biogeochemical Cycles, 19. doi:10.1029/2004GB002340
Bibliographic reference Burchard, Hans ; Deleersnijder, Eric ; Meister, Andreas. Application of modified Patankar schemes to stiff biogeochemical models for the water column. In: Ocean Dynamics : theoretical, computational oceanography and monitoring, Vol. 55, no.3-4, p. 326-337 (2005)
Permanent URL http://hdl.handle.net/2078.1/173772