<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1">
<META NAME="Generator" CONTENT="MS Exchange Server version 6.5.7638.1">
<TITLE>RE : [dune-pdelab] Instationary Navier-Stokes equations with Finite Volume method</TITLE>
</HEAD>
<BODY>
<!-- Converted from text/plain format -->
<P><FONT SIZE=2>Hello,<BR>
<BR>
First of all, thanks for your prompt answer !<BR>
<BR>
Once I have written my own operator and run my simulation, do you have any tools (functions ?) for debugging ? By debugging, I mean dumping the values of a matrix in case the maths are wrong<BR>
<BR>
Thanks again<BR>
Samuel<BR>
<BR>
<BR>
-------- Message d'origine--------<BR>
De: Christian Engwer [<A HREF="mailto:christian.engwer@uni-muenster.de">mailto:christian.engwer@uni-muenster.de</A>]<BR>
Date: jeu. 11/09/2014 14:36<BR>
À: CERRITO Samuel<BR>
Cc: dune-pdelab@dune-project.org; Xavier.KUHN@cea.fr<BR>
Objet : Re: [dune-pdelab] Instationary Navier-Stokes equations with Finite Volume method<BR>
<BR>
Hi,<BR>
<BR>
> I use classes from stokesdg.hh as local operators and P0LocalFiniteElementMap as finite element map for pressure and velocity components. My idea was to turn discontinuous Galerkin case into finite volume method by using a zero order polynomial (P0LocalFiniteElementMap).<BR>
><BR>
> This method doesn't seem to work, do you think the classes I use are appropriate to solve these equations ? I am aware of the existence of StokesDGVectorFEM but I'm not quite sure about how it works.<BR>
<BR>
the problem is not DUNE related, but due to the maths. The p0/p0 pair<BR>
does not fulfill the inf-sup condition and thus you get wrong<BR>
results. Furthermore the DG method is not a FV method, although they<BR>
have a lot of similarities. In particular, you have to be careful<BR>
regarding the definition of your fluxes.<BR>
<BR>
> Are there any other local operators I can try ?<BR>
<BR>
you should write your own operator. And! more important first think<BR>
about which discretization you want to use.<BR>
<BR>
You'll have to either use a pair like P1/P0, or use a stabilization,<BR>
like LPS, or other techniques.<BR>
<BR>
Ciao<BR>
Christian<BR>
<BR>
<BR>
</FONT>
</P>
</BODY>
</HTML>