[Nek5000-users] Linearised Navier Stokes

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


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

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


More information about the Nek5000-users mailing list