MFiX中后处理输出雷诺数Re的方法

操作方法

废话少说，直接上操作方法
从\legacy_tests\dem-tests\evaporation\usr_rates_des.f中复制代码如下

!----------------------------------------------------------------------!
! Function: calc_NRe(M)                                                !
!                                                                      !
! Purpose: Calculate the Reynolds number.                              !
!                                                                      !
!''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''''!
      DOUBLE PRECISION FUNCTION cal_NRe(M)

      INTEGER M

! Various fluid cell indicies
      INTEGER I, IMJK, IJMK, IJKM
! Gas Velocity - cell centered
      DOUBLE PRECISION UGC, VGC, WGC
! Solids Velocity - cell centered
      DOUBLE PRECISION USCM, VSCM, WSCM
! Relative velocity.
      DOUBLE PRECISION VREL

! Initialize fluid cell variables
      I =  I_OF(IJK)
      IMJK  = IM_OF(IJK)
      IJMK  = JM_OF(IJK)
      IJKM  = KM_OF(IJK)

! Calculate velocity components at i, j, k
! Gas
      UGC = AVG_X_E(U_G(IMJK),U_G(IJK),I)
      VGC = AVG_Y_N(V_G(IJMK),V_G(IJK))
      WGC = AVG_Z_T(W_G(IJKM),W_G(IJK))
! Solids
      USCM = DES_VEL_NEW(NP,1)
      VSCM = DES_VEL_NEW(NP,2)
      WSCM = DES_VEL_NEW(NP,3)

! magnitude of gas-solids relative velocity
      VREL = SQRT((UGC - USCM)**2 + (VGC - VSCM)**2 + (WGC - WSCM)**2)

! Reynods Number
      IF(MU_g(IJK) > ZERO) THEN
         cal_NRe = 2.0d0 * DES_RADIUS(NP) * VREL * RO_g(IJK) / MU_g(IJK)
      ELSE
         cal_NRe = LARGE_NUMBER
      ENDIF

      RETURN
      END FUNCTION cal_NRe

粘贴到自己文件夹下面的des_usr_rates.f中
位置如图
并且在粘贴的代码之前加上如下代码（如上图所示：

! Reynolds Number
      N_Re = cal_NRe(1)
      des_usr_var(12,NP)=N_Re

       IF (PM==1) THEN
          write(*,*) 'N_Re=', N_Re
       ENDIF 
      
      RETURN
      
      contains

注意别忘了contains这个关键字
其中RETURN是des_usr_rates本来就有的

这样在des_usr_var(12,:)这个变量里就能输出Re数了。（des_usr_var用法见前面的博客）
并且屏幕上也会打印出雷诺数
流化床里雷诺数的量级在几十上下

之前也写过一个博客，删了。那个没什么用，看这个。

完毕

---