[Nek5000-users] "usrdiv" not working with Pn/Pn-2 formulation

[nek5000-users at lists.mcs.anl.gov]

Thanks Harish,

We have changed induct.f to include the usrdiv call.  The current revision
will now reflect this change.

-Katie

On Wed, Jan 11, 2012 at 10:57 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hello,
>
> I was trying to use the "subroutine fill_div(usrdiv)" with the Pn/Pn-2
> formulation for a backward facing step simulation. This used to work
> before, now it is not working in my simulation.
>
> I did a little bit of digging by reverting to old Nek revisions and found
> the following:
>
> 1. Until revision 739, "usrdiv" was called in the main code with the
> following lines for Pn/Pn-2 formulation in file navier4.f.
>
> subroutine incompr()
> .......
> .......
> .......
> call add2col2(respr,bm2,usrdiv,ntot2) ! User-defined divergence
> .......
> .......
>
> 2. The "subroutine incompr()" is called from "subroutine plan3()", which
> is in planx.f
>
> 3. From revision 740 "subroutine plan3()" calls "subroutine incomprn()"
> and not "subroutine incompr()". And incompr() is deleted from the
> navier4.f.
>
> 4. "usrdiv" is missing in the new incomprn() subroutine. Hence, the value
> of "usrdiv" computed in *.usr file is not added to the main solution.
>
> Just wanted to point this problem out. This may be useful for others in
> turbulent simulations.
>
> I got around this problem by interpolating an old solution onto the new
> grid.
>
> Regards,
>
> Harish.
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20120112/bc3d44e1/attachment.html>