Using H-Convergence to Calculate the Numerical Errors for 1D Unsaturated Seepage Under Steady-State Conditions

Robert P. Chapuis, Coline Taveau, François Duhaime, Simon Weber, Vahid Marefat, Lu Zhang, Daniela Blessent, Najib Bouaanani, Dominique Pelletier

Producción científica: Contribución a una revistaArtículorevisión exhaustiva

Resumen

Unsaturated zones are important for geotechnical design, geochemical reactions, and microbial reactions. The numerical analysis of unsaturated seepage is complex because it involves highly nonlinear partial differential equations. The permeability can vary by orders of magnitude over short vertical distances. This article defines and uses H-convergence tests to quantify numerical errors made by uniform meshes with element size (ES) for 1D steady-state conditions. The quantitative H-convergence should not be confused with a qualitative mesh sensitivity study. The difference between numerical and mathematical convergences is stated. A detailed affordable method for an H-convergence test is presented. The true but unknown solution is defined as the asymptote of the numerical solutions for all solution components when ES decreases to zero. The numerical errors versus ES are then assessed with respect to the true solution, and using a log–log plot, which indicates whether a code is correct or incorrect. If a code is correct, its results follow the rules of mathematical convergence in a mathematical convergence domain (MCD) which is smaller than the numerical convergence domain (NCD). If a code is incorrect, it has an NCD but no MCD. Incorrect algorithms of incorrect codes need to be modified and repaired. Existing codes are shown to converge numerically within large NCDs but generate large errors, up to 500%, in the NCDs, a dangerous situation for designers.

Idioma originalInglés
Páginas (desde-hasta)345-358
Número de páginas14
PublicaciónInternational Journal for Numerical and Analytical Methods in Geomechanics
Volumen49
N.º1
DOI
EstadoAceptada/en prensa - 2024

Tipos de productos de Minciencias

  • Artículo A1 - Q1

Huella

Profundice en los temas de investigación de 'Using H-Convergence to Calculate the Numerical Errors for 1D Unsaturated Seepage Under Steady-State Conditions'. En conjunto forman una huella única.

Citar esto