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

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

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.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20120111/7b6f305a/attachment.html>