[Nek5000-users] Linear Advection Equation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sun Mar 8 09:17:40 CDT 2015
Thank You, Paul.
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?
On Sat, Feb 28, 2015 at 10:31 AM, <nek5000-users at lists.mcs.anl.gov> wrote:
>
>
>
> Set ifflow to false and ifheat to T and ifadvc = F T F F F F...
>
> in the .rea file.
>
> Then, supply your velocity field in useric (ux, uy, uz), or, if you
> want time-varying velocity,
> do so in usrchk by over-writing the velocity field arrays, vx(), vy(), and
> vz() at each step.
>
>
> Paul
> ------------------------------
> *From:* nek5000-users-bounces at lists.mcs.anl.gov [
> nek5000-users-bounces at lists.mcs.anl.gov] on behalf of
> nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
> *Sent:* Saturday, February 28, 2015 9:24 AM
> *To:* nek5000-users at lists.mcs.anl.gov
> *Subject:* [Nek5000-users] Linear Advection Equation
>
> Dear Users,
>
> 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:
>
> D(phi)/Dt = d (phi) /dt + v*.*grad(phi)= 0
>
> where v*.*grad(phi) = u.d(phi)/dx + v.d(phi)/dy + w.d(phi)/dz
>
> and v is a prescribed velocity (for example v = (1,1) in 2D)
>
> Best,
>
> Saumil
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>
--
Saumil S. Patel
Ph.D. Candidate
Department of Mechanical Engineering
Computational Multiphase Fluid Dynamics Group
<https://sites.google.com/site/leeccny/>
The City College of New York
Telephone: 1 (212) 650-7134
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20150308/caa65f30/attachment.html>
More information about the Nek5000-users
mailing list