<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 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.<br>
Can you please confirm that my adjustments are reasonable?<br>
<br>
In conduct.f when hsolve is called, I have adjusted the argument
list like this:</p>
<p>- napproxt(1) = laxtt ! Fix this... pff 10/10/15<br>
+ napproxt(1,ifield-1) = laxtt ! Fix this... pff 10/10/15<br>
<br>
call hsolve (name4t,TA,TB,H1,H2 <br>
$ ,tmask(1,1,1,1,ifield-1)<br>
$ ,tmult(1,1,1,1,ifield-1)<br>
$ ,imesh,tolht(ifield),nmxh,1<br>
- $ ,approxt,napproxt,bintm1)<br>
+ $
,approxt(1,1,ifield-1),napproxt(1,ifield-1)<br>
+ $ ,bintm1)<br>
<br>
such that approxt and napproxt refers to each ifield, ie.
temperature, passive scalar 1, passive scalar 2, ... similar to
velx, vely, velz.</p>
<p><br>
</p>
<p>Accordingly, I have adjusted ORTHOT:</p>
<p>- integer napproxt(2)<br>
+<br>
+C Add another dimension to (n)approxt for each passive scalar and
temp.<br>
+ integer napproxt(2,ldimt)<br>
common /trthoi/ napproxt<br>
<br>
- real approxt(ktott,0:laxtt)<br>
+ real approxt(ktott,0:laxtt,ldimt)<br>
<br>
</p>
<p>and induct.f</p>
<p>- napproxt(1) = laxtt<br>
+ napproxt(1,ifield-1) = laxtt<br>
<br>
</p>
<p><br>
</p>
<pre wrap="">
</pre>
<br>
<div class="moz-cite-prefix">On 09/14/2017 06:29 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.780.1505406553.10190.nek5000-users@lists.mcs.anl.gov">
<pre wrap="">Date: Thu, 14 Sep 2017 15:34:39 +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.770.1505403282.10190.nek5000-users@lists.mcs.anl.gov" moz-do-not-send="true"><mailman.770.1505403282.10190.nek5000-users@lists.mcs.anl.gov></a>
Content-Type: text/plain; charset="windows-1252"
Hi Steffen,
Thanks for bringing this to our attention.
Adam and I have been looking into this... (no doubt, it's something
stupid done by me awhile back).
Hope to have it resolved by next week.
Best,
Paul
________________________________
From: Nek5000-users <a class="moz-txt-link-rfc2396E" href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" moz-do-not-send="true"><nek5000-users-bounces@lists.mcs.anl.gov></a> on behalf of <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> <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>
Sent: Thursday, September 14, 2017 8:18:27 AM
To: <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>
Subject: Re: [Nek5000-users] proj_ortho error in pipe flow
Hi,
I have made some more tests to find the cause of this issue.
Attached is a laminar pipe flow simulation with temperature and one
passive scalar. The original setup is from Paul's "nusselt.tgz", where
he answered another question of mine.
I found that I can reproduce the same error I see in my turbulent
simulations when there is an additional passive scalar with a different
type of boundary condition, i.e. Neumann boundary condition at the pipe
wall for temperature and Dirichlet boundary condition for the passive
scalar.
Can someone please check my setup and tell me why this causes the
projection to fail?
I hope to save some computational resources by using the projection.
However, if it is not possible, I will keep param(94) = 0.
Kind Regards,
Steffen Straub</pre>
</blockquote>
<br>
</body>
</html>