In this work we present finite element approximations of relaxed systems for nonlinear diffusion problems, which can also tackle the cases of degenerate and strongly degenerate diffusion equations. Relaxation schemes take advantage of the replacement of the original partial differential equation (PDE) with a semilinear hyperbolic system of equations, with a stiff source term, tuned by a relaxation parameter ε. When ε → 0+, the system relaxes onto the original PDE: in this way, a consistent discretization of the relaxation system for vanishing ε yields a consistent discretization of the original PDE. The numerical schemes obtained with this procedure do not require solving implicit nonlinear problems and possess the robustness of upwind discretizations. The proposed approximations are based on a discontinuous Galerkin method in space and on suitable implicitexplicit integration in time. Then, in principle, we can achieve any order of accuracy and obtain stable solutions, even when the diffusion equation becomes degenerate and solution singularities develop. Moreover, when needed, we can easily incorporate slope limiters within our schemes in order to handle spurious oscillatory phenomena. Some preliminary theoretical results are given, along with several numerical tests in one and two space dimensions, both for linear and nonlinear diffusion problems, including a degenerate diffusion equation, that provide numerical evidence of the properties of the presented approach.

Cavalli, F., Naldi, G., Perugia, I., Discontinuous Galerkin Approximation of Relaxation Models for Linear and Nonlinear Diffusion Equations, <<SIAM JOURNAL ON SCIENTIFIC COMPUTING>>, 2012; 34 (1): A105-A136. [doi:10.1137/110827752] [http://hdl.handle.net/10807/85480]

Discontinuous Galerkin Approximation of Relaxation Models for Linear and Nonlinear Diffusion Equations

Cavalli, Fausto
Primo
;
2012

Abstract

In this work we present finite element approximations of relaxed systems for nonlinear diffusion problems, which can also tackle the cases of degenerate and strongly degenerate diffusion equations. Relaxation schemes take advantage of the replacement of the original partial differential equation (PDE) with a semilinear hyperbolic system of equations, with a stiff source term, tuned by a relaxation parameter ε. When ε → 0+, the system relaxes onto the original PDE: in this way, a consistent discretization of the relaxation system for vanishing ε yields a consistent discretization of the original PDE. The numerical schemes obtained with this procedure do not require solving implicit nonlinear problems and possess the robustness of upwind discretizations. The proposed approximations are based on a discontinuous Galerkin method in space and on suitable implicitexplicit integration in time. Then, in principle, we can achieve any order of accuracy and obtain stable solutions, even when the diffusion equation becomes degenerate and solution singularities develop. Moreover, when needed, we can easily incorporate slope limiters within our schemes in order to handle spurious oscillatory phenomena. Some preliminary theoretical results are given, along with several numerical tests in one and two space dimensions, both for linear and nonlinear diffusion problems, including a degenerate diffusion equation, that provide numerical evidence of the properties of the presented approach.
Inglese
Cavalli, F., Naldi, G., Perugia, I., Discontinuous Galerkin Approximation of Relaxation Models for Linear and Nonlinear Diffusion Equations, <>, 2012; 34 (1): A105-A136. [doi:10.1137/110827752] [http://hdl.handle.net/10807/85480]
File in questo prodotto:
Non ci sono file associati a questo prodotto.

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: http://hdl.handle.net/10807/85480
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 3
  • ???jsp.display-item.citation.isi??? 5
social impact