[Nek5000-users] Linearised Navier Stokes

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Aug 10 07:43:54 CDT 2017


that looks quite suspicious, indeed. I would "only" be a problem for 3D 
linearised temperature. Let me discuss with Adam next week.

Philipp

On 2017-08-10 14:33, nek5000-users at lists.mcs.anl.gov wrote:
> Hi all,
> 
> Speaking of which, I think there is a small mistake in subroutine 
> convabp in the perturb.f file.
> 
> 
> if (.not.ifadj) then
> 
> call opcopy(ua,ub,uc,vx,vy,vz)
> 
> call opcopy(vx,vy,vz,vxp(1,jp),vyp(1,jp),v*y*p(1,jp))
> 
> call convop(ta,t(1,1,1,1,ifield-1)) ! dU.grad T
> 
> call opcopy(vx,vy,vz,ua,ub,uc)
> 
> endif
> 
> 
> where  v*y*p(1,jp) should be v*z*p(1,jp)
> 
> best,
> 
> Pierre-Emmanuel
> 
> 
> 
> On 2017-08-10 14:21, nek5000-users at lists.mcs.anl.gov wrote:
>> Hi,
>> the linearised version is included in the main version of the code. 
>> You activate perturbations by setting NPERT different from zero. The 
>> relevant routines in the code are collected in the file perturb.f.
>>
>> Hope this helps,
>> philipp
>>
>> On 2017-08-10 12:56, nek5000-users at lists.mcs.anl.gov wrote:
>>> Hi Nek Forum,
>>>
>>>      How to obtain the linearised version of the Nek5000 code. Please 
>>> advice.
>>>
>>> Thanks
>>>   Swarandeep
>>> _______________________________________________
>>> Nek5000-users mailing list
>>> Nek5000-users at lists.mcs.anl.gov
>>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>>>
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> 
> 
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> 


More information about the Nek5000-users mailing list