[Nek5000-users] Time derivative of Scalar field

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

Thanks Ketan.

Saikat,

You may also be able to just use the values in vxlag() and thus avoid the need
for extra storage.

A simple test would verify..   Of course, it's not that expensive to just store your
own copies for purposes of your own analysis.

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: Tuesday, April 25, 2017 1:06 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Time derivative of Scalar field

Hi Saikat,

I have attached a usr file that shows how to calculate finite difference based first derivative of the first passive scalar t(1,1,1,1,2) using the current time step and two previous time steps. The coefficients for the derivative are obtained using a call to fd_weights_full, and I am using a common block for storing the values of the scalar at two previous time-steps.

Ketan








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