化学反应新鲜出炉

rainman08

第十三例，酸碱反应

[ 本帖最后由 rainman08 于 2008-11-11 09:45 编辑 ]

rainman08

LIBRARY:
  CEL:
    EXPRESSIONS:
      Rate=4.*density*ed/ke
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      HeatReaction=460. [kJ kg^-1]
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      i=2.7905 [kg kg^-1]
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      AcidSource=-Rate*min(acid.mf,alkali.mf/i)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      AcidSourceCoeff=AcidSource/max(acid.mf,1.0e-16)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      AlkaliSource=-Rate*min(alkali.mf,acid.mf*i)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      AlkaliSourceCoeff=AlkaliSource/max(alkali.mf,1.0e-16)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      ProductSource=(1.0+i)*Rate*min(acid.mf,alkali.mf/i)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      HeatSource=HeatReaction*Rate*min(acid.mf,alkali.mf/i)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      alpha=1.0E-05 [kmole kg^-1]
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      Xions=alpha*density*(acid.mf+product.mf/(1+i))
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      Yions=alpha*density*(alkali.mf/i + product.mf/(1+i))
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      a=1. [kg kg^-1]
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      b=Yions - Xions
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      c=-1.0E-14 [kmole^2 m^-6]
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      Hions=(-b+sqrt(max(b*b-4*a*c,0.0 [ kmole^2 m^-6 ])))/(2*a)
    END
  END
END
LIBRARY:
  CEL:
    EXPRESSIONS:
      pH=-log10(max( (Hions/ 1. [kmol m^-3]), 1.0E-15))
    END
  END
END

ustbdynamic

呵呵  不错啊