<div dir="ltr">Thank You, Paul.  <br><br>Why does the conductivity/diffusivity (i.e. param(08)) have to be small and non-zero positive?  Can nek5k solve the linear advection equation without the discrete diffusion term?    <br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Feb 28, 2015 at 10:31 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>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>Set ifflow to false and ifheat to T and ifadvc = F T F F F F...</div>
<div><br>
</div>
<div>in the .rea file.</div>
<div><br>
</div>
<div>Then, supply your velocity field in useric  (ux, uy, uz), or, if you want time-varying velocity,</div>
<div>do so in usrchk by over-writing the velocity field arrays, vx(), vy(), and vz() at each step.</div>
<div><br>
</div>
<div><br>
</div>
<div>Paul<br>
<div style="font-family:Times New Roman;color:#000000;font-size:16px">
<hr>
<div style="direction:ltr"><font color="#000000" face="Tahoma"><b>From:</b> <a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users-bounces@lists.mcs.anl.gov" target="_blank">nek5000-users-bounces@lists.mcs.anl.gov</a>] on behalf of <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> [<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>]<br>
<b>Sent:</b> Saturday, February 28, 2015 9:24 AM<br>
<b>To:</b> <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a><br>
<b>Subject:</b> [Nek5000-users] Linear Advection Equation<br>
</font><br>
</div><div><div class="h5">
<div></div>
<div>
<div dir="ltr">
<div>
<div>Dear Users,<br>
<br>
</div>
Is there a way to solve a linear advection equation in nek5k without computing the Navier-Stokes equations?  I am simply interested in advecting a profile according to the following:<br>
<br>
</div>
D(phi)/Dt = d (phi) /dt +  v<b>.</b>grad(phi)= 0<br>
<div><br>
where v<b>.</b>grad(phi) = u.d(phi)/dx + v.d(phi)/dy + w.d(phi)/dz<br>
<br>
</div>
<div>and v is a prescribed velocity (for example v = (1,1) in 2D) <br>
</div>
<div><br>
</div>
<div>Best,<br>
<br>
</div>
<div>Saumil <br clear="all">
</div>
<div>
<div>
<div><br>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
</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><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr">Saumil S. Patel<br>Ph.D. Candidate<br>Department of Mechanical Engineerin<font>g<br><a href="https://sites.google.com/site/leeccny/" dir="ltr" target="_blank">Computational Multiphase Fluid Dynamics Group</a></font><br>The City College of New York<br>Telephone: <a value="+12126507134">1 (212) 650-7134</a></div></div>
</div>