<div dir="ltr">Hi<div style>My problem is flow in a 2d channel. The initial condition is zero velocity in the channel and there is an inlet flow. I have a passive scalar which is flowing in through inlet. Initial temperature is zero everywhere. I am having trouble getting initial iterations to converge. I am trying stress formulation and passive scalars for first time and I do not see what is going wrong.</div>
<div style><br></div><div style><br></div><div style><div> Solving for heat</div><div> Solving for fluid F F T</div><div> Temperature/Passive scalar solution</div><div> 1.0000000000000001E-009 p22 1 2</div>
<div> New CG1-tolerance (RINIT*epsm) = 8.4585499354425527E-033 2. 5514287912844953E-033</div><div> 1 1 Helmholtz TEMP F: 8.4585E-20 1.0000E-09 2.8571E-04 0.0000E+00</div><div> 1 Hmholtz TEMP: 0 8.4585E-20 8.4585E-20 1.0000E-09</div>
<div> 1.0000000000000001E-009 p22 1 2</div><div> New CG1-tolerance (RINIT*epsm) = 5.8174583141330990E-032 1. 7860001538991468E-032</div><div> 1 1 Helmholtz PS 1 F: 5.8175E-19 1.0000E-09 2.0000E-03 0.0000E+00</div>
<div> 1 Hmholtz PS 1: 0 5.8175E-19 5.8175E-19 1.0000E-09</div><div> 1 0.0000E+00 6.0799E-03 Heat done</div><div> New CG1(stress)-tolerance (RINIT*epsm) = 6.9206392325228944E-014 4. 2096427334076814E-033</div>
<div> 0 1.0000E-09 6.9206E-14 1.0000E-09 tol,matmod</div><div> 1 tx3: 1 0.3584E+00 0.1000E-08 0.3584E+00</div><div> 1 tx3: 2 0.2546E+00 0.1000E-08 0.3584E+00</div><div> 1 tx3: 3 0.1905E+00 0.1000E-08 0.3584E+00</div>
<div>.</div><div>.</div><div>.</div><div><div> 1 tx3: 198 0.1080E-02 0.1000E-08 0.3584E+00</div><div> 1 tx3: 199 0.1067E-02 0.1000E-08 0.3584E+00</div><div> 1 tx3: 200 0.1033E-02 0.1000E-08 0.3584E+00</div>
<div> 1 201 Unconverged Helmholtz3/Fluid: rbnorm = 0.103306E-02 0.100000E-08</div><div> 1 1.00000E-07 2.34297E-01 2.34297E-01 1.00000E+00 1 Divergence</div><div> New CG1(stress)-tolerance (RINIT*epsm) = 1.5438583947746548E-016 4. 2096427334076814E-033</div>
<div> 0 1.0000E-09 1.5439E-16 1.0000E-09 tol,matmod</div><div> 1 tx3: 1 0.8855E-02 0.1000E-08 0.8855E-02</div><div> 1 tx3: 2 0.5463E-02 0.1000E-08 0.8855E-02</div><div> 1 tx3: 3 0.5217E-02 0.1000E-08 0.8855E-02</div>
</div><div><br></div><div style>Thanks a lot</div><div style>praveen</div></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sat, Oct 12, 2013 at 10:04 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Hi Praveen,<br>
<br>
That looks right... As always, it's best to run a test<br>
to verify first.<br>
<div class="im HOEnZb"><br>
Paul<br>
<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
</div><div class="HOEnZb"><div class="h5">Sent: Saturday, October 12, 2013 11:17:06 AM<br>
Subject: Re: [Nek5000-users] Variable density flows<br>
<br>
<br>
<br>
Thanks to Paul and Ammar. I will try the variable density approach.<br>
<br>
<br>
I am first trying a case in which viscosity depends on temperature and a passive scalar. My uservp looks like this. Can I use "temp" and "PS" to get temperature and passive scalar as below ?<br>
<br>
<br>
<br>
subroutine uservp (ix,iy,iz,ieg)<br>
include 'SIZE'<br>
include 'TOTAL'<br>
include 'NEKUSE'<br>
common /mygeom/ xmin,xmax<br>
common /myparam/ r_k,Rey_no,Pr_no,Sc_no,visc_rat<br>
<br>
<br>
if(ifield.eq.1)then ! momentum equation<br>
utrans= 1.0<br>
udiff = exp(PS(1)*log(visc_rat) - temp)/Rey_no<br>
else if(ifield.eq.2)then ! temperature eqn<br>
utrans= 1.0<br>
udiff = PS(1)*(r_k - 1.0) + 1.0<br>
udiff = udiff/(Rey_no*Pr_no)<br>
else if(ifield.eq.3)then ! passive scalar<br>
utrans= 1.0<br>
udiff = 1.0/(Rey_no*Sc_no)<br>
else<br>
print*,'Only 3 fields'<br>
call exitt()<br>
endif<br>
<br>
<br>
return<br>
end<br>
<br>
<br>
Thanks<br>
praveen<br>
</div></div><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</div></div></blockquote></div><br></div>