<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p><br>
</p>
<p>Dear Steffen,</p>
<p><br>
</p>
<p>I believe this issue may have been resolved in an update to the git repo</p>
<p>last week.</p>
<p><br>
</p>
<p>Best,</p>
<p><br>
</p>
<p>Paul</p>
<p><br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Thursday, September 28, 2017 7:40:27 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] proj_ortho error in pipe flow</font>
<div> </div>
</div>
<div>
<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>
</div>
</body>
</html>