<div dir="ltr">Hello Ammar<div style>The initial c lies between 0 and 1 and I want it to stay in this range. The problem is that there are narrow interfaces in which c changes from 0 to 1 and the change will not be monotone causing c to takes values outside [0,1].</div>
<div style>Thanks</div><div style>praveen</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Oct 13, 2013 at 7:42 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"><div dir="auto"><div>I may have misunderstood your Q</div><div>Did you mean that c initially has a max = 1 and that mixing will only cause it to go down? </div>
<div>In this case check mass conservation</div><div>And make sure the initial density is accurately specified and that it gets updated correctly</div><div class="im"><div>Ammar</div><div><br>Sent from my iPhone</div></div>
<div><div class="h5"><div><br>On Oct 13, 2013, at 9:01 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br><br></div><blockquote type="cite"><div><div><div style="direction:ltr">
<span>I am not sure that what you are solving for c is a normalized scalar like a mass fraction Y for example </span></div><div style="direction:ltr"><span>If c is not normalized May be the formulation is equivalent to</span></div>
<div style="direction:ltr"><span>C= rho* Y where y is a normalized scalar and rho is the total density..</span></div><div style="direction:ltr"><span><br></span></div><div style="direction:ltr"><span>Finally check mass conservation</span></div>
<div style="direction:ltr"><span>By integrating over the volume of the channel, the density equation</span></div><div style="direction:ltr"><span>d/dt integral(rho) + integral (u. Grad rho) =error </span></div><div style="direction:ltr">
<span>And see how error changes over time</span></div><div style="direction:ltr"><span><br></span></div><div style="direction:ltr"><span>Ammar</span></div><br><span>Sent from my iPhone</span></div><div><br>On Oct 13, 2013, at 1:46 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
<br></div><blockquote type="cite"><div><div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Oct 13, 2013 at 10:35 AM, <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"><div>You can also try error function like profile with 0 value at the inlet and increases to a peak somewhere and stays constant up to the exit so gradient is0 at the exit<br>
Insulation normally refers to the diffusive flux (first form)</div><div><div>Ammar</div></div></blockquote></div><br>Thanks for the tips. I am using a smoothed initial condition for scalar now. However after some time the scalar does not stay within [0,1] which is important. My time step is already small with cfl ~= 0.1. What can I do to get better behaviour of the scalar ?</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">The diffusion coefficient for my scalar is quite small, D = 5*10^(-5).</div><div class="gmail_extra"><br></div><div class="gmail_extra">Regards</div>
<div class="gmail_extra">praveen</div></div>
</div></blockquote><blockquote type="cite"><div><span>_______________________________________________</span><br><span>Nek5000-users mailing list</span><br><span><a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a></span><br>
<span><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a></span><br></div></blockquote></div></blockquote><blockquote type="cite">
<div><span>_______________________________________________</span><br><span>Nek5000-users mailing list</span><br><span><a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a></span><br>
<span><a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a></span><br></div></blockquote></div></div></div><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>
<br></blockquote></div><br></div>