[Nek5000-users] About the iobj in the drag calculation

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

Hi,
    I'm sorry to trouble. I find these lines in set_obj subroutine
            if (f.eq.1) iobj=1  ! lower wall
            if (f.eq.3) iobj=2  ! upper wall
    so this defines which wall to be represented. Another question comes that why the urms is calculate like
    do i=1,m
                write(56,3) yy(i)+1
     &                     ,(yy(i)+1)*Re_tau
     &                     , (urms_pl(i)-(uavg_pl(i))**2)/u_tau**2
     &                     , vrms_pl(i)/u_tau**2
     &                     , wrms_pl(i)/u_tau**2
     &                     , uvms_pl(i)/u_tau**2
                write(57,3)  yy(i) + 1.
     &                     , (yy(i)+1.)*Re_tau
     &                     , uavg_pl(i)/u_tau
    3           format(1p15e17.9)
            enddo

As can be seen,R_uu=(urms_pl(i)-(uavg_pl(i))**2)/u_tau**2, while R_vv=vrms_pl(i)/u_tau**2,why?






在 2014-05-29 02:24:34，nek5000-users at lists.mcs.anl.gov 写道：

Hi,all:
    I read Paul's reply in this mail https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2014-May/002756.html
    And got to know something about objects. As you say, the object is a collection of faces, here , I have a question about the iboj, how could I know which face it represent with iboj? As in turbChannel.usr,
     dragx_avg = alpha*dragx_avg + beta*0.5*(dragx(1)+dragx(2))
     I think it means to do the time-average of the average dragx on planey1 and planey2, as
     0.5*(dragx(1)+dragx(2))
 indicates, while how to calculate the dragx on the planex1, planex2, planez1 and planez2? for the maxobj in this example is only 4 and dragx is defined as
     dragx(0:maxobj)
   
Thank you
Xianbei




来自网易手机号码邮箱了解更多
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140529/7f3aaf87/attachment.html>