<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
Dear nek5000 team,
<div><br>
</div>
<div>Could you help me with the 2 questions below? Thank you very much !</div>
<div><br>
</div>
<div>1)</div>
<div>Concerning the routine “axhelm” in hmholtz.f. If I am right, one possibility of this routine can be the discretisation of the diffusion operator. However it is not clear to me if axhelm performs:</div>
<div><br>
</div>
<div>a) vdiff  x   Laplacian(phi)</div>
<div><br>
</div>
<div>or  </div>
<div><br>
</div>
<div>b) grad ( vdiff x grad (phi) )</div>
<div><br>
</div>
<div>phi being the variable solved for (T or PS for example).</div>
<div><br>
</div>
<div>So, is it actually a) or b) ?</div>
<div><br>
</div>
<div><br>
</div>
<div>2) </div>
<div>In the routine wlaplacian in navier1.f, I am not sure to understand the consequence of the line:</div>
<div><br>
</div>
<div>call sub2 (out,wrk,ntot)</div>
<div><br>
</div>
<div>after the lines:</div>
<div><br>
</div>
<div>
<div>call bcneusc(out,1)</div>
<div>call axhelm(wrk,a,diff,h2,imesh,1)</div>
</div>
<div><br>
</div>
<div>If I am right, this will change the sign in front of the diffusion operator, right? </div>
<div>Because in makeq.f, I am not sure that the line</div>
<div>call add2(bq(1,1,1,1,ifield-1),w1,ntot)</div>
<div><br>
</div>
<div>will actually assemble the RHS of the transport equation :  - (vtrans) x  U grad (phi)  +  grad ( vdiff x grad (phi) )</div>
<div><br>
</div>
<div><br>
</div>
<div>Thank you for your help !</div>
<div><br>
<div>
<div style="color: rgb(0, 0, 0); letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div>Best regards,</div>
<div>Emmanuel</div>
</div>
<br class="Apple-interchange-newline">
<br class="Apple-interchange-newline">
</div>
<br>
</div>
</body>
</html>