<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<p>Hi Neks,<br>
<br>
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. <br>
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).<br>
<br>
I tried to restrict the projection only for velocity and
temperature in hsolve subroutine by adding the following code
snippet:<br>
</p>
<p>...<br>
if (ifstdh) then<br>
call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)<br>
! My restriction for using projection only on ifield 1 and 2 but
not on 3 (first passive scalar)<br>
elseif (ifield.gt.2) then<br>
call hmholtz(name,u,r,h1,h2,vmk,vml,imsh,tol,maxit,isd)<br>
else<br>
...<br>
<br>
<br>
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. <br>
I have tested this in a simple test case which is uploaded here:<br>
</p>
<p><a class="moz-txt-link-freetext" href="https://bwsyncandshare.kit.edu/dl/fiVn2UB9QqWZN7awPsJHJSVF">https://bwsyncandshare.kit.edu/dl/fiVn2UB9QqWZN7awPsJHJSVF</a><br>
<br>
<br>
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.<br>
The adjustments from my previous message (to extend approxt and
napproxt) only postponed the projection errors to later timesteps.<br>
</p>
<p>Maybe you can point me in the right direction?<br>
<br>
Best Regards,<br>
Steffen Straub<br>
<br>
<br>
<br>
<br>
</p>
<br>
<div class="moz-cite-prefix">On 09/28/2017 07:00 PM,
<a class="moz-txt-link-abbreviated" href="mailto:nek5000-users-request@lists.mcs.anl.gov">nek5000-users-request@lists.mcs.anl.gov</a> wrote:<br>
</div>
<blockquote type="cite"
cite="mid:mailman.23.1506618014.75224.nek5000-users@lists.mcs.anl.gov">
<pre wrap="">Message: 1
Date: Thu, 28 Sep 2017 14:40:27 +0200
From: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">nek5000-users@lists.mcs.anl.gov</a>
To: <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users@lists.mcs.anl.gov></a>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
Message-ID:
<a class="moz-txt-link-rfc2396E" href="mailto:mailman.558.1506602432.6921.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.558.1506602432.6921.nek5000-users@lists.mcs.anl.gov></a>
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 <i class="moz-txt-slash"><span class="moz-txt-tag">/</span>trthoi<span class="moz-txt-tag">/</span></i> 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: <a class="moz-txt-link-abbreviated" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">nek5000-users@lists.mcs.anl.gov</a>
To: <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true">"nek5000-users@lists.mcs.anl.gov"</a>
<a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users@lists.mcs.anl.gov></a>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
Message-ID:
<a class="moz-txt-link-rfc2396E" href="mailto:mailman.565.1506607685.6921.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.565.1506607685.6921.nek5000-users@lists.mcs.anl.gov></a>
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</pre>
</blockquote>
<br>
</body>
</html>