[Nek5000-users] Lyapunov Exponents

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

Hi Alireza,

this routine was donated by some of our users and is definitely
on the experimental side of things... My suggestion would be to
look carefully at it and see if you can understand it.   As I recall,
a key point was being able to come up with a meaningful norm when
both velocity and temperature are involved in the eigenvector.

Paul


On Mon, 23 Jan 2012, nek5000-users at lists.mcs.anl.gov wrote:

> Hi,
>
> I want to compute the perturbation vectors and Lyapunov exponents in a 
> simulation of Rayleigh-Benard convection. But I realized that the 
> "computelyap" subroutine in perturb.f file is never called anywhere in the 
> code. I wonder if you could help me how to handle the calculation of the 
> Lyapunov exponents and store the corresponding perturbation fld files and lyp 
> file.
>
> Regards,
> Alireza
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users