We construct finite volume schemes of very high order of accuracy in space and time for solving the nonlinear Richards equation (RE). The general scheme is based on a three-stage predictor-corrector procedure. First, a high-order weighted essentially non-oscillatory (WENO) reconstruction procedure is applied to the cell averages at the current time level to guarantee monotonicity in the presence of steep gradients. Second, the temporal evolution of the WENO reconstruction polynomials is computed in a predictor stage by using a global weak form of the governing equations. A global space-time DG FEM is used to obtain a scheme without the parabolic time-step restriction caused by the presence of the diffusion term in the RE. The resulting nonlinear algebraic system is solved by a Newton-Krylov method, where the generalized minimal residual method algorithm of Saad and Schulz is used to solve the linear subsystems. Finally, as a third step, the cell averages of the finite volume method are updated using a one-step scheme, on the basis of the solution calculated previously in the space-time predictor stage. Our scheme is validated against analytical, experimental, and other numerical reference solutions in four test cases. A numerical convergence study performed allows us to show that the proposed novel scheme is high order accurate in space and time. © 2011 John Wiley & Sons, Ltd.

A Novel Numerical Method of High-Order Accuracy for Flow in Unsaturated Porous Media / C., Zambra; Dumbser, Michael; Toro, Eleuterio Francisco; N., Moraga. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - STAMPA. - 89:2(2012), pp. 227-240. [10.1002/nme.3241]

A Novel Numerical Method of High-Order Accuracy for Flow in Unsaturated Porous Media

Dumbser, Michael;Toro, Eleuterio Francisco;
2012-01-01

Abstract

We construct finite volume schemes of very high order of accuracy in space and time for solving the nonlinear Richards equation (RE). The general scheme is based on a three-stage predictor-corrector procedure. First, a high-order weighted essentially non-oscillatory (WENO) reconstruction procedure is applied to the cell averages at the current time level to guarantee monotonicity in the presence of steep gradients. Second, the temporal evolution of the WENO reconstruction polynomials is computed in a predictor stage by using a global weak form of the governing equations. A global space-time DG FEM is used to obtain a scheme without the parabolic time-step restriction caused by the presence of the diffusion term in the RE. The resulting nonlinear algebraic system is solved by a Newton-Krylov method, where the generalized minimal residual method algorithm of Saad and Schulz is used to solve the linear subsystems. Finally, as a third step, the cell averages of the finite volume method are updated using a one-step scheme, on the basis of the solution calculated previously in the space-time predictor stage. Our scheme is validated against analytical, experimental, and other numerical reference solutions in four test cases. A numerical convergence study performed allows us to show that the proposed novel scheme is high order accurate in space and time. © 2011 John Wiley & Sons, Ltd.
2012
2
C., Zambra; Dumbser, Michael; Toro, Eleuterio Francisco; N., Moraga
A Novel Numerical Method of High-Order Accuracy for Flow in Unsaturated Porous Media / C., Zambra; Dumbser, Michael; Toro, Eleuterio Francisco; N., Moraga. - In: INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING. - ISSN 0029-5981. - STAMPA. - 89:2(2012), pp. 227-240. [10.1002/nme.3241]
File in questo prodotto:
File Dimensione Formato  
NME-Richards1D.pdf

accesso aperto

Tipologia: Versione editoriale (Publisher’s layout)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 250.4 kB
Formato Adobe PDF
250.4 kB Adobe PDF Visualizza/Apri

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11572/88334
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 26
  • ???jsp.display-item.citation.isi??? 25
social impact