An Improved Scheme for the Finite Difference Approximation of the Advective Term in the Heat or Solute Transport Equations
Transport equations are widely used to describe the evolution of scalar quantities subject to advection, dispersion and, possibly, reactions. Numerical methods are required to solve these equations in applications, adopting either the advective or conservative formulations. Conservative formulations...
| Autores: | , |
|---|---|
| Tipo de recurso: | artículo |
| Estado: | Versión publicada |
| Fecha de publicación: | 2024 |
| País: | España |
| Institución: | Consejo Superior de Investigaciones Científicas (CSIC) |
| Repositorio: | DIGITAL.CSIC. Repositorio Institucional del CSIC |
| OAI Identifier: | oai:digital.csic.es:10261/371324 |
| Acceso en línea: | http://hdl.handle.net/10261/371324 https://api.elsevier.com/content/abstract/scopus_id/85206993906 |
| Access Level: | acceso abierto |
| Palabra clave: | Solute transport Advective transport Finite difference scheme Heat transport Mass conservation http://metadata.un.org/sdg/11 http://metadata.un.org/sdg/13 Make cities and human settlements inclusive, safe, resilient and sustainable Take urgent action to combat climate change and its impacts |
| Sumario: | Transport equations are widely used to describe the evolution of scalar quantities subject to advection, dispersion and, possibly, reactions. Numerical methods are required to solve these equations in applications, adopting either the advective or conservative formulations. Conservative formulations are usually preferred in practice because they conserve mass. Advective formulations do not, but have received more mathematical attention and are required for Lagrangian solution methods. To obtain an advective formulation that conserves mass, we subtract the discretized fluid flow equation, multiplied by concentration, from the conservative form of the transport equation. The resulting scheme not only conserves mass, but is also elegant in that it can be interpreted as averaging the advective term at cell interfaces, instead of approximating it at cell centers as in traditional centered schemes. The two schemes are identical when fluid velocity is constant, and both have second-order convergence, but the truncation errors are slightly different. We argue that the error terms appearing in the proposed scheme actually imply an improved representation of subgrid spreading/contraction and acceleration/deceleration caused by variable velocity. We compare the proposed and traditional schemes on several problems with variable velocity caused by recharge, discharge or evaporation, including two newly developed analytical solutions. The proposed method yields results that are slightly, but consistently, better than the traditional scheme, while always conserving mass (i.e., mass at the end equals mass at the beginning plus inputs minus outputs), which the traditional centered finite differences scheme does not. We conclude that this scheme should be preferred in finite difference solutions of transport. |
|---|