On Tue, Nov 21, 2017 at 7:10 AM, Tobias Moehle <tobias.moe...@uni-rostock.de
> wrote:

> Dear all,
>
> I have some problem with Lagrange multipliers: Maybe I misinterpret my
> results, but somehow they don't have the effect I expected.
>
> Looking at the solution of systems_of_equations_ex3, (using paraview) I
> see the (real part of) pressure r_p to be 0 at time step 0 and than it
> jumps to  ~10, where it stays for the rest of the simulation. But isn't it
> actually enforced to stay at 0 "forever"?
>


This is beacuse the following lines:

  // We can set the mean of the pressure by setting Falpha.  Typically
  // a value of zero is chosen, but the value should be arbitrary.
  navier_stokes_system.rhs->add(navier_stokes_system.rhs->size()-1, 10.);

So in this case we set the mean pressure to be 10. We do this just for
illustration purposes, to show that the Lagrange multiplier can be used to
set the mean pressure to any value you like (but as mentioned in the
comment normally one would just set it to zero). If you want to change the
value, just change the "10" in the line above.

David
------------------------------------------------------------------------------
Check out the vibrant tech community on one of the world's most
engaging tech sites, Slashdot.org! http://sdm.link/slashdot
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to