[Nek5000-users] proj_ortho error in pipe flow

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

Hi Neks,

I saw that there was a new commit on Github "import Paul's projection 
fix"  on October, 4th so I tried the newest Nek5000 version with my setup.
Projection seems to be running fine if only temperature is solved 
additionally to pressure and velocity. However, when I want to solve for 
an additional passive scalar the projection fails (for both unperturbed 
initial conditions as well as chaotic fields).

I tried to restrict the projection only for velocity and temperature in 
hsolve subroutine by adding the following code snippet:

...
if (ifstdh) then
     call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)
! My restriction for using projection only on ifield 1 and 2 but not on 
3 (first passive scalar)
elseif (ifield.gt.2) then
     call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)
else
...


If I restrict the projection in such a way, I do not get errors for the 
projection of the temperature field but on the other hand no projection 
for passive scalar(s) is done either.
I have tested this in a simple test case which is uploaded here:

https://bwsyncandshare.kit.edu/dl/fiVn2UB9QqWZN7awPsJHJSVF


I believe there is a problem when more than one temperature field is 
present and one activates projection (p94 != 0). Unfortunately, my 
knowledge of the code is not sufficient enough to really dig deeper and 
solve this problem.
The adjustments from my previous message (to extend approxt and 
napproxt) only postponed the projection errors to later timesteps.

Maybe you can point me in the right direction?

Best Regards,
Steffen Straub





On 09/28/2017 07:00 PM, nek5000-users-request at lists.mcs.anl.gov wrote:
> Message: 1
> Date: Thu, 28 Sep 2017 14:40:27 +0200
> From:nek5000-users at lists.mcs.anl.gov
> To:<nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
> Message-ID:
> 	<mailman.558.1506602432.6921.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="utf-8"; Format="flowed"
>
> Hi Neks,
>
> I might have fixed my proj_ortho error issue. At least it is running now
> without errors on a small test case on my local machine. However, I am
> very much unsure if what I was changing is correct.
> Can you please confirm that my adjustments are reasonable?
>
> In conduct.f when hsolve is called, I have adjusted the argument list
> like this:
>
> -????? napproxt(1) = laxtt? ! Fix this... pff 10/10/15
> +????? napproxt(1,ifield-1) = laxtt? ! Fix this... pff 10/10/15
>
>   ??????????? call hsolve? (name4t,TA,TB,H1,H2
>   ????? $?????????????????? ,tmask(1,1,1,1,ifield-1)
>   ????? $?????????????????? ,tmult(1,1,1,1,ifield-1)
>   ????? $?????????????????? ,imesh,tolht(ifield),nmxh,1
> -???? $?????????????????? ,approxt,napproxt,bintm1)
> +???? $ ,approxt(1,1,ifield-1),napproxt(1,ifield-1)
> +???? $?????????????????? ,bintm1)
>
> such that approxt and napproxt refers to each ifield, ie. temperature,
> passive scalar 1, passive scalar 2, ... similar to velx, vely, velz.
>
>
> Accordingly, I have adjusted ORTHOT:
>
> -????? integer???????? napproxt(2)
> +
> +C Add another dimension to (n)approxt for each passive scalar and temp.
> +????? integer???????? napproxt(2,ldimt)
>   ?????? common/trthoi/  napproxt
>
> -????? real??????????? approxt(ktott,0:laxtt)
> +????? real??????????? approxt(ktott,0:laxtt,ldimt)
>
> and induct.f
>
> -????? napproxt(1) = laxtt
> +????? napproxt(1,ifield-1) = laxtt
>
>
>
>
>
> Message: 2
> Date: Thu, 28 Sep 2017 14:04:58 +0000
> From:nek5000-users at lists.mcs.anl.gov
> To:"nek5000-users at lists.mcs.anl.gov"
> 	<nek5000-users at lists.mcs.anl.gov>
> Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
> Message-ID:
> 	<mailman.565.1506607685.6921.nek5000-users at lists.mcs.anl.gov>
> Content-Type: text/plain; charset="us-ascii"
>
>
> Dear Steffen,
>
>
> I believe this issue may have been resolved in an update to the git repo
>
> last week.
>
>
> Best,
>
>
> Paul

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171010/773f2f28/attachment.html>