> There is another strange problem, I change my RKDG to MOMOMIAL but I can
> only
>
> ran succesful when I use CONSTANT, SECOND and FOURTH basis, FIRST and THIRD
>
> is unsteady, can your epplitic code works well with FIRST or THIRD?
Yes, it works. I have tested it till P8 and as you can see libMesh
does both h, p, and hp refinement.
> My code solve euler equation and I am trying to solve NS equaiton later.
I'm interested in solving NS too.
Actually I have also some code for the Stokes equation but I miss the
convective flux treatment
and so I would like to implement an Euler solver in the near future.
I wish I had time to take a look at the matrix free approach.
My Stokes code uses an artificial compressibility flux approach but
I'm not happy of the convergence I
can achive with GMRES in large 3D problems. I add a couple of
screenshot of the stokes solution
of the lid-driven cavity with P4 and h_ref, Re is 400.
I'm very interested to know what kind of scheme you are using and how
it performs with libMesh.
Thank you very much for the reply.
Lorenzo
-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
