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

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

You are welcome!

Harish.

On Thu, Jan 12, 2012 at 12:08 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

>
> 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
>>
>>
>
> _______________________________________________
> 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/fea70394/attachment.html>