lorenzo alessio botti <[email protected]> writes:

> I almost never tune those parameters.
> Interior penalty can be problematic in case of very stretched elements but
> BR2 works great in 2D and 3D with stability parameter \eta= 3, 4 or number
> of element faces (as the theory suggests).

How do you define the stabilization when the viscosity/diffusivity has a
discontinuity across the interface?  Make it a healthy jump of ten
orders of magnitude.

> I think that the major difference is the convergence rates of iterative
> solvers, with cG you need far less iterations. That's the reason why I use
> cG for the pressure solvers in operator splitting algorithms for
> incompressible flows.

Okay, but then your discrete velocity field is only divergence-free up
to discretization error instead of up to iterative solver tolerance.

------------------------------------------------------------------------------
Precog is a next-generation analytics platform capable of advanced
analytics on semi-structured data. The platform includes APIs for building
apps and a phenomenal toolset for data science. Developers can use
our toolset for easy data analysis & visualization. Get a free account!
http://www2.precog.com/precogplatform/slashdotnewsletter
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to