1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Subroutine: CALC_K_g C 4 ! Purpose: Calculate the effective conductivity of fluid phase C 5 ! C 6 ! C 7 ! Comments: C 8 ! This routine will not be called if k_g0 is defined C 9 ! C 10 ! C 11 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 12 SUBROUTINE CALC_K_G() 13 14 ! Modules 15 !---------------------------------------------------------------------// 16 USE param1, only: undefined 17 USE physprop, only: k_g0, k_g 18 USE sendrecv, only: send_recv 19 ! invoke user defined quantity 20 USE usr_prop, only: usr_kg, calc_usr_prop 21 USE usr_prop, only: gas_conductivity 22 IMPLICIT NONE 23 !---------------------------------------------------------------------// 24 25 IF (USR_Kg) THEN 26 CALL CALC_USR_PROP(Gas_Conductivity,lm=0) 27 ELSEIF (K_g0 == UNDEFINED) THEN 28 ! unncessary check but included for clarity 29 CALL CALC_DEFAULT_Kg 30 ENDIF 31 32 CALL send_recv(K_G, 2) 33 34 RETURN 35 END SUBROUTINE CALC_K_G 36 37 38 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 39 ! C 40 ! Purpose: Compute the default value for gas conductivity where the C 41 ! gas phase is assumed to be air C 42 ! Author:M. Syamlal Date: 24-APR-96 C 43 ! C 44 ! Literature/Document References: C 45 ! Bird, Stewart, and Lightfoot (1960) -- C 46 ! Temperature dependence from formula 8.3-12 on p. 255 and C 47 ! conductivity value at 300 K from p. 263 C 48 ! C 49 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 50 SUBROUTINE CALC_DEFAULT_Kg 51 52 ! Modules 53 !---------------------------------------------------------------------// 54 use compar, only: ijkstart3, ijkend3 55 USE fldvar, only: T_g 56 USE functions, only: fluid_at 57 USE param1, only: zero 58 USE physprop, only: K_g 59 USE run, only: units 60 IMPLICIT NONE 61 62 ! Local variables 63 !---------------------------------------------------------------------// 64 ! Indices 65 INTEGER :: IJK 66 !---------------------------------------------------------------------// 67 68 !!$omp parallel do private(ijk) & 69 !!$omp& schedule(dynamic,chunk_size) 70 DO IJK = ijkstart3, ijkend3 71 IF (FLUID_AT(IJK)) THEN 72 ! Gas conductivity (air) in cal/(s.cm.K) 73 K_G(IJK) = 6.02D-5*SQRT(T_G(IJK)/300.D0) 74 ELSE 75 K_G(IJK) = ZERO 76 ENDIF 77 78 ! 1 cal = 4.183925D0 J 79 IF (UNITS == 'SI') K_G(IJK) = 418.3925D0*K_G(IJK) !J/s.m.K 80 81 ENDDO 82 83 RETURN 84 END SUBROUTINE CALC_DEFAULT_Kg 85