[Nek5000-users] Linear Advection Equation

[nek5000-users at lists.mcs.anl.gov]

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20150228/342ef106/attachment.html>