محاسبات فلش با استفاده از معادله حالت Soave-Redlich-Kwong در Mathematica
کد:
(*CacheID: 232*)
(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[ 10777, 321]*)
(*NotebookOutlinePosition[ 11421, 343]*)
(* CellTagsIndexPosition[ 11377, 339]*)
(*WindowFrame->Normal*)
Notebook[{
Cell[CellGroupData[{
Cell[TextData[{
StyleBox["Flash Calculations using the Soave-Redlich-Kwong Equation of \
State", "Subsubtitle",
FontSize->18,
FontVariations->{"CompatibilityType"->0}],
StyleBox["\n",
FontWeight->"Plain",
FontVariations->{"CompatibilityType"->0}],
StyleBox["Authors' Data",
FontSize->14,
FontWeight->"Bold"],
StyleBox[": Zakia Nasri",
FontSize->14],
StyleBox[" and",
TextAlignment->Center,
FontSize->13,
Background->RGBColor[0.605478, 0.996109, 0.605478]],
StyleBox[" Housam Binous",
FontSize->14],
StyleBox["\n",
TextAlignment->Center,
FontFamily->"MS Shell Dlg",
FontSize->8.5,
Background->RGBColor[0.605478, 0.996109, 0.605478]],
StyleBox["Department of Chemical Engineering\nNational Institute of Applied \
Sciences and Technology\nTunis, TUNISIA\nEmail: binoushousam@yahoo.com ",
FontSize->14,
FontWeight->"Plain"]
}], "Title",
TextAlignment->Center,
Background->RGBColor[0.605478, 0.996109, 0.605478]],
Cell[TextData[{
StyleBox["Acknowledgement :",
FontWeight->"Bold"],
"\nThis problem was presented in the Master's report of R. Privat under the \
direction of J-P Corriou."
}], "Subsubtitle",
Background->RGBColor[0.773449, 0.996109, 0.996109]],
Cell[BoxData[
\(Off[General::"\<spell1\>"]\)], "Input"],
Cell[CellGroupData[{
Cell["Reduced Pressure and Temperature", "Subsubtitle",
Background->RGBColor[1, 1, 0.658824]],
Cell[BoxData[
\(Tre[i_] := T/Tc[i]\)], "Input"],
Cell[BoxData[
\(Pre[i_] := P/Pc[i]\)], "Input"]
}, Closed]],
Cell[CellGroupData[{
Cell[TextData[StyleBox["Equilibrium constants using the \[Phi]-\[Phi] methode \
and the Soave-Redlich-Kwong Equation of State", "Subsubtitle",
FontSize->18,
FontVariations->{"CompatibilityType"->0}]], "Subsubtitle",
Background->RGBColor[1, 1, 0.658824]],
Cell[BoxData[
\(Ap[i_] := 0.42747\ a[i]\ Pre[i]/Tre[i]^2\)], "Input"],
Cell[BoxData[
\(Bp[i_] := 0.08664\ Pre[i]/Tre[i]\)], "Input"],
Cell[BoxData[
\(\[Phi]\_v[i_] :=
Exp[\((Z\_v - 1)\)\ Bp[i]/B\_v - Log[Z\_v - B\_v] -
A\_v/B\_v\ \((2\ Ap[i]^0.5/A\_v^0.5 -
Bp[i]/B\_v)\)\ Log[\((Z\_v + B\_v)\)/Z\_v]]\)], "Input"],
Cell[BoxData[
\(\[Phi]\_l[i_] :=
Exp[\((Z\_l - 1)\)\ Bp[i]/B\_l - Log[Z\_l - B\_l] -
A\_l/B\_l\ \((2\ Ap[i]^0.5/A\_l^0.5 -
Bp[i]/B\_l)\)\ Log[\((Z\_l + B\_l)\)/Z\_l]]\)], "Input"],
Cell[BoxData[
\(a[i_] := \((1 + m[i]\ \((1 - Tre[i]^0.5)\))\)^2\)], "Input"],
Cell[BoxData[
\(m[i_] := 0.480 + 1.574\ \[Omega][i] - 0.176\ \[Omega][i]^2\)], "Input"],
*)
