|
| 1 | +//****** diffusion coeff *********** |
| 2 | +D_H2 1.2097e-8; //Ngu et al., Chem. Engg. Sci., 251, 2022, 117478 |
| 3 | +D_CO2 5.3765e-9; //Ngu et al., Chem. Engg. Sci., 251, 2022, 117478 |
| 4 | +D_CO 3.9518e-9; //Ngu et al., Chem. Engg. Sci., 251, 2022, 117478 |
| 5 | +//****** Henry coeff *************** |
| 6 | +He_H2 0.02; |
| 7 | +He_CO 0.018; |
| 8 | +He_CO2 0.64; |
| 9 | +//****** molecular weight *************** |
| 10 | +Mw_CO2 0.044;//kg/mol |
| 11 | +Mw_CO 0.028; |
| 12 | +Mw_H2 0.002; |
| 13 | +//****Liquid properties************** |
| 14 | +CpMixLiq 4181; |
| 15 | +muMixLiq 0.001; // 1 cP |
| 16 | +kThermLiq 0.62; // W/m-K |
| 17 | +rho0MixLiq 1000; // kg/m^3 |
| 18 | +sigmaLiq 0.07; //surface tension N/m |
| 19 | +//*******inlet gas frac************* |
| 20 | +f_H2 0.07; |
| 21 | +f_CO2 0.70; |
| 22 | +f_CO 0.23; |
| 23 | +//********************************* |
| 24 | +LeLiqH2 #calc "$kThermLiq / $rho0MixLiq / $D_H2 / $CpMixLiq"; |
| 25 | +LeLiqCO #calc "$kThermLiq / $rho0MixLiq / $D_CO / $CpMixLiq"; |
| 26 | +LeLiqCO2 #calc "$kThermLiq / $rho0MixLiq / $D_CO2 / $CpMixLiq"; |
| 27 | +LeLiqMix #calc "$f_H2*$LeLiqH2+$f_CO2*$LeLiqCO2+$f_CO*$LeLiqCO"; |
| 28 | +//********************************* |
| 29 | +kH2 #calc "$D_H2*$rho0MixLiq*$CpMixLiq*$LeLiqMix"; |
| 30 | +PrH2 #calc "$muMixLiq*$CpMixLiq / $kH2"; |
| 31 | + |
| 32 | +kCO #calc "$D_CO*$rho0MixLiq*$CpMixLiq*$LeLiqMix"; |
| 33 | +PrCO #calc "$muMixLiq*$CpMixLiq / $kCO"; |
| 34 | + |
| 35 | +kCO2 #calc "$D_CO2*$rho0MixLiq*$CpMixLiq*$LeLiqMix"; |
| 36 | +PrCO2 #calc "$muMixLiq*$CpMixLiq / $kCO2"; |
| 37 | +//********************************* |
| 38 | +HtBcol 10.0; |
| 39 | +DiaBcol 0.72; |
| 40 | +LiqHt 7.0; |
| 41 | +NPS 8; |
| 42 | +NPD 8; |
| 43 | +NPY 200; |
| 44 | +//********************************* |
| 45 | +T0 303.15; |
| 46 | +//P0 101325.0; |
| 47 | +P0 1.0E5; |
| 48 | +Pbot #calc "$P0+1000.0*$LiqHt*9.8"; |
| 49 | +Pmid #calc "$P0+1000.0*0.5*$LiqHt*9.8"; |
| 50 | + |
| 51 | +//bubble column cross-section divided by sparger area |
| 52 | +ArbyAs 2.0; |
| 53 | +uSupVel 0.05; |
| 54 | +uGas #calc "$Pmid / $Pbot * $ArbyAs * $uSupVel"; |
| 55 | +A_Bcol #calc "3.14159 * $DiaBcol * $DiaBcol / 4.0"; |
| 56 | +rho_gas #calc "$Pmid / 287.0 / $T0"; |
| 57 | +mflowRate #calc "$uSupVel * $A_Bcol * $rho_gas"; |
| 58 | +//mflowRate 4.5e-4; |
| 59 | +//********************************* |
0 commit comments