In this paper, the unsteady isothermal flow of a gas through a semi-infinite micro nano porous medium described by a non-linear two-point boundary value problem on a semi-infinite interval has been considered. We solve this problem by a nonstandard finite difference method defined on quasi-uniform grids in order to derive a new numerical approximation. By introducing a stencil that is constructed in such a way that the boundary conditions at infinity are exactly assigned, the proposed method is effectively used to determine the numerical solution. In addition, a mesh refinement and the Richardson's extrapolation allow to improve the accuracy of the numerical solution and to define a posteriori estimator for the global error of the proposed numerical scheme. We determine the accurate initial slope du/dx(0) = 1.1917906497194208 calculated for a = 0.5 with good capturing the essential behavior of u(x). This clearly demonstrates that the numerical solutions presented in this paper result highly accurate and in excellent agreement with the existing solutions available in the literature.
Numerical study on gas flow through a micro–nano porous medium based on finite difference schemes on quasi-uniform grids
Fazio, Riccardo;Jannelli, Alessandra
;ROTONDO, TIZIANA
2018-01-01
Abstract
In this paper, the unsteady isothermal flow of a gas through a semi-infinite micro nano porous medium described by a non-linear two-point boundary value problem on a semi-infinite interval has been considered. We solve this problem by a nonstandard finite difference method defined on quasi-uniform grids in order to derive a new numerical approximation. By introducing a stencil that is constructed in such a way that the boundary conditions at infinity are exactly assigned, the proposed method is effectively used to determine the numerical solution. In addition, a mesh refinement and the Richardson's extrapolation allow to improve the accuracy of the numerical solution and to define a posteriori estimator for the global error of the proposed numerical scheme. We determine the accurate initial slope du/dx(0) = 1.1917906497194208 calculated for a = 0.5 with good capturing the essential behavior of u(x). This clearly demonstrates that the numerical solutions presented in this paper result highly accurate and in excellent agreement with the existing solutions available in the literature.File | Dimensione | Formato | |
---|---|---|---|
3129195.pdf
solo gestori archivio
Descrizione: Articolo principale
Tipologia:
Versione Editoriale (PDF)
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
656.9 kB
Formato
Adobe PDF
|
656.9 kB | Adobe PDF | Visualizza/Apri Richiedi una copia |
Pubblicazioni consigliate
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.