[Nek5000-users] Instantaneous velocity field

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Feb 21 14:43:02 CST 2018


this is not really instantaneous data, but rather statistics. Then you 
would need to run some kind of averaging (in time, space etc.) and 
proper postprocessing afterwards, which strongly depends on the actual 
flow case. I am sure that there are people in this group that have done 
this type of analysis for a variety of flow cases/geometries.

On 2018-02-21 19:34, nek5000-users at lists.mcs.anl.gov wrote:
> Hello,
> 
> Thank you for the reply. Yes, outpost() in fact does yield the 
> information I am looking for.  But I hope I could generate the field 
> information in .dat file similar to following pattern:
> 
>           open(unit=56,file='vel_fluc_prof.dat')
> 
>           write(56,'(A,1pe14.7)') '#time = ', time
> 
>           write(56,'(A)')
> 
>       &   '#  y     y+     uu     vv     ww     uv'
> 
> Could you please teach me how to do this?
> 
> Regards,
> 
> Tim
> 
> ------------------------------------------------------------------------
> 
> I forgot to add 'call' before outpost. Please read the block as:
> 
>    if (istep .eq. 10) then
> 
>         call outpost(vx,vy,vz,pr,t,'vel')
> 
>     endif
> 
> On 02/21/2018 04:33 PM, Navrose Navrose wrote:
> 
>>//
> 
>>/Doesn't outpost() (userchk subroutine) inside a conditional statement /
> 
>>/address your issue?/
> 
>>/Something like (if you have to flow information at istep=10):/
> 
>>//
> 
>>/  if (istep .eq. 10) then/
> 
>>/       outpost(vx,vy,vz,pr,t,'vel')/
> 
>>/   endif/
> 
>>//
> 
>>//
> 
>>/On 02/21/2018 04:22 PM, nek5000-users at lists.mcs.anl.gov 
> <https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users> wrote:/
> 
>>>//
> 
>>>/Hello,/
> 
>>>//
> 
>>>/Thank you for the reply. I meant the instantaneous velocity field for /
> 
>>>/the domain (u(x,y,z),v(x,y,z),w(x,y,z)) at some certain time step /
> 
>>>/(istep)./
> 
>>>//
> 
>>>/Regards,/
> 
>>>//
> 
>>>/Tim/
> 
> 
> 
> _______________________________________________
> 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