Dear friends,

I work out the things further. So, finally, I realized that the code is not taking the base flow component of the scalar field. 

I mean to say that for a scalar field T, 
T = T_base + T_ fluctuation, 

but somehow in the code it is

T = T_ fluctuation.

I would be grateful if David or any other developer could check this.

with regards
Abhishek 

---------------------------------------------------------------------------------------
Abhishek Kumar
Post Doctoral Researcher
Applied Mathematics Research Centre
Coventry University, Coventry CV15FB
The United Kingdom

Email # ac7600@coventry.ac.uk, abhishek.kir@gmail.com
---------------------------------------------------------------------------------------


 


On Tue, Jan 30, 2018 at 5:12 PM, Abhishek Kumar <abhishek.kir@gmail.com> wrote:
Dear friends,

I am trying to do direct stability analysis of Rayleigh Benard Convection (RBC) in Nektar++. 

I observe positive eigenvalues even for Ra=1000, which is below critical. I have attached my XML file. It would be great if someone would give comment on this issue.


With regards
Abhishek Kumar

---------------------------------------------------------------------------------------
Abhishek Kumar
Post Doctoral Researcher
Applied Mathematics Research Centre
Coventry University, Coventry CV15FB
The United Kingdom

Email # ac7600@coventry.ac.uk, abhishek.kir@gmail.com
---------------------------------------------------------------------------------------



--
With regards
Abhishek Kumar