Flux difference splitting and the balancing of source terms and flux gradients.

*(English)*Zbl 0972.65056Authors’ summary: Flux difference splitting methods are widely used for the numerical approximation of homogeneous conservation laws where the flux depends only on the conservative variables. However, in many practical situations this is not the case. Not only are source terms commonly part of the mathematical model, but also the flux can vary spatially even when the conservative variables do not.

It is the discretisation of the additional terms arising from these two situations which is addressed in this work, given that a specific flux difference splitting method has been used to approximate the underlying conservation law. The discretisation is constructed in a manner which retains an exact balance between the flux gradients and the source terms when this is appropriate.

The effectiveness of these new techniques, in both one and two dimensions, is illustrated using the shallow water equations, in which the additional terms arise from the modelling of bed slope and, in one dimension, breadth variation. Roe’s scheme is chosen for the approximation of the conservation laws and appropriate discrete forms are constructed for the additional terms, not only in the first-order case but also in the presence of flux- and slope-limited high-resolution corrections. The method is then extended to two-dimensional flow where it can be applied on both quadrilateral and triangular grids.

It is the discretisation of the additional terms arising from these two situations which is addressed in this work, given that a specific flux difference splitting method has been used to approximate the underlying conservation law. The discretisation is constructed in a manner which retains an exact balance between the flux gradients and the source terms when this is appropriate.

The effectiveness of these new techniques, in both one and two dimensions, is illustrated using the shallow water equations, in which the additional terms arise from the modelling of bed slope and, in one dimension, breadth variation. Roe’s scheme is chosen for the approximation of the conservation laws and appropriate discrete forms are constructed for the additional terms, not only in the first-order case but also in the presence of flux- and slope-limited high-resolution corrections. The method is then extended to two-dimensional flow where it can be applied on both quadrilateral and triangular grids.

Reviewer: Dinh Nho Hao (Frankfurt am Main)

##### MSC:

65M06 | Finite difference methods for initial value and initial-boundary value problems involving PDEs |

35L65 | Hyperbolic conservation laws |

76B15 | Water waves, gravity waves; dispersion and scattering, nonlinear interaction |

76M12 | Finite volume methods applied to problems in fluid mechanics |

##### Keywords:

source terms; finite volume schemes; upwind schemes; flux limiters; slope limiters; conservation laws; flux difference splitting methods; shallow water equations; Roe’s scheme##### Software:

MPDATA
PDF
BibTeX
XML
Cite

\textit{M. E. Hubbard} and \textit{P. Garcia-Navarro}, J. Comput. Phys. 165, No. 1, 89--125 (2000; Zbl 0972.65056)

Full Text:
DOI

##### References:

[1] | Bermúdez, A.; Vázquez, M.E., Upwind methods for hyperbolic conservation laws with source terms, Comput. fluids, 23, 1049, (1994) · Zbl 0816.76052 |

[2] | Bermúdez, A.; Dervieux, A.; Desideri, J.-A.; Vázquez, M.E., Upwind schemes for the two-dimensional shallow water equations with variable depth using unstructured meshes, Comput. methods appl. mech. eng., 155, 49, (1998) · Zbl 0961.76047 |

[3] | Cunge, J.A.; Holly, F.M.; Verwey, A., Practical aspects of computational river hydraulics, (1980) |

[4] | Goutal, N.; Maurel, F., Proceedings of the 2nd workshop on dam-break wave simulation, (1997) |

[5] | Garcia-Navarro, P.; Vazquez-Cendon, M.E., Some considerations and improvements on the performance of Roe’s scheme for 1d irregular geometries, (1997) |

[6] | Glaister, P., Difference schemes for the shallow water equations, (1987) · Zbl 0967.76068 |

[7] | Glaister, P., Prediction of supercritical flow in open channels, Comput. math. appl., 24, 69, (1992) · Zbl 0763.76049 |

[8] | Hubbard, M.E., On the accuracy of one-dimensional models of steady converging/diverging open channel flows, (1999) · Zbl 1014.76051 |

[9] | Huang, J.; Song, C.S., Stability of dynamic flood routing schemes, J. hydraul. eng. ASCE, 111, (1985) |

[10] | M. E. Hubbard, and, N. Dodd, A 2d numerical model of wave run-up, overtopping and regeneration, in preparation. |

[11] | Hudson, J., Numerical techniques for the shallow water equations, (1999) |

[12] | Jenny, P.; Muller, B., Rankine – hugoniot – riemann solver considering source terms and multidimensional effects, J. comput. phys., 145, 575, (1998) · Zbl 0926.76079 |

[13] | LeVeque, R.J., Numerical methods for conservation laws, (1992) · Zbl 0847.65053 |

[14] | LeVeque, R.J., Balancing source terms and flux gradients in high-resolution Godunov methods: the quasi-steady wave-propagation algorithm, J. comput. phys., 146, 346, (1998) · Zbl 0931.76059 |

[15] | Priestley, A., Roe-type schemes for super-critical flows in rivers, (1989) |

[16] | Roe, P.L., Approximate Riemann solvers, parameter vectors, and difference schemes, J. comput. phys., 43, 357, (1981) · Zbl 0474.65066 |

[17] | P. L. Roe, Fluctuations and signals—A framework for numerical evolution problems, in Numerical Methods for Fluid Dynamics, edited by K. W. MortonOxford Univ. Press, London, 1982, 219-257. |

[18] | Roe, P.L., Characteristic-based schemes for the Euler equations, Ann. rev. fluid mech., 18, 337, (1986) · Zbl 0624.76093 |

[19] | Smolarkiewicz, P.K.; Margolin, L.G., MPDATA: A finite-difference solver for geophysical flows, J. comput. phys., 140, 459, (1998) · Zbl 0935.76064 |

[20] | Sweby, P.K., High resolution schemes using flux limiters for hyperbolic conservation laws, SIAM J. numer. anal., 21, 995, (1984) · Zbl 0565.65048 |

[21] | Vázquez-Cendón, M.E., Improved treatment of source terms in upwind schemes for the shallow water equations in channels with irregular geometry, J. comput. phys., 148, 497, (1999) · Zbl 0931.76055 |

[22] | van Leer, B., Towards the ultimate conservative difference scheme V. A second order sequel to Godunov’s method, J. comput. phys., 32, 101, (1979) · Zbl 1364.65223 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.