Open AccessElementos finitos em formulação mista de mínimos quadrados para a simulação da convecção-difusão em regime transiente
Author(s) -
Vitoriano Ruas,
Regina Célia P. Leal-Toledo,
Maurício Kischinhevsky
Publication year - 2013
Publication title -
revista internacional de métodos numéricos para cálculo y diseño en ingeniería
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.213
H-Index - 9
eISSN - 1886-158X
pISSN - 0213-1315
DOI - 10.1016/j.rimni.2012.12.002
Subject(s) - mathematics , piecewise , polynomial , boundary (topology) , divergence (linguistics) , discretization , mathematical analysis , finite element method , similarity solution , dimension (graph theory) , boundary layer , pure mathematics , physics , philosophy , linguistics , thermodynamics
AbstractA mixed least-squares finite element scheme designed for solving the transient convection-diffusion equations expressed in terms of both the primal unknown and its flux, incorporating or not a reaction term, is studied. Once a time discretization of the Crank-Nicholson type is performed, the resulting system of equations allows for a stable approximation of both fields, by means of classical Lagrange continuous piecewise polynomial functions of arbitrary degree, in both simplicial and non-simplicial geometry, in any space dimension. The scheme is also convergent in space in the mean-square sense as far as the primal unknown field, its gradient, the flux variable and its divergence are concerned, and in time in an appropriate sense for each one of these four fields. Numerical results certify that the scheme performs well for any Péclet number, thereby allowing to confirm theoretical predictions, at least in the case where there is no narrow boundary layer. In the latter case however the method fails to produce reliable results. The technique is also compared with three existing methods to solve the convection-diffusion equations in the transient case. These include two recent ones proposed by the first author and collaborators