<div dir="ltr"><div class="gmail_extra">These tips are most useful. I am able to run some high Re and high Pe cases but I still need really large meshes to ensure that the scalar stays bounded between [0,1].</div><div class="gmail_extra">
<br></div><div class="gmail_extra">Is the characteristic time stepping going to be useful for these cases ? I have read that it allows larger time steps. Since "characteristic" means some upwind effect, I was wondering does it help to maintain better monotonicity of advected scalars ? </div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Thanks</div><div class="gmail_extra">praveen</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 14, 2013 at 12:37 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 id=":1gg" style="overflow:hidden">Hi Praveen,<br>
<br>
Your clipping implementation will as you have it -- I normally would<br>
do it in a single loop:<br>
<br>
      n=nx1*ny1*nz1*nelt<br>
      do i=1,n<br>
         t(i,1,1,1,2)=min(1.,t(i,1,1,1,2))<br>
         t(i,1,1,1,2)=max(0.,t(i,1,1,1,2))<br>
      enddo<br>
<br>
You can turn on filtering by setting param(103) in your .rea file.<br>
Setting p103=.05, say, would reduce the energy in the Nth mode by<br>
5% each timestep --- this is usually a good choice.  Sometimes it<br>
pays to bump this up to, say, .1, but usually a better strategy is<br>
to add another mode to the filter, by setting p101 to 1, or 2.<br>
In this case, the filter amplitude for the N-1 mode is alpha/4<br>
when p101=1 and when p101=2, we have amplitudes as:<br>
<br>
   a(N) = alpha := p103<br>
   a(N-1) = alpha*(2/3)**2 = 4*alpha/9<br>
   a(N-2) = alpha*(1/3)**2 = alpha/9<br>
<br>
which amounts to fitting a parabola through the filter amplitudes so<br>
that the filter strength decays smoothly to zero.<br>
<br>
p103 = .05<br>
p101 = 1<br>
<br>
are my common choices when lx1=8.  If lx1 < 8, I typically set p101=0<br>
as you're otherwise cutting too deeply into low end of the spectrum.<br>
<br>
All this being said, you'll no doubt still need clipping if Sc=100,000.<br>
<div class="im"><br>
Paul</div></div></blockquote></div><br><br></div></div>