Ammonia viscosity and tests

Signed-off-by: Ian Bell <ian.h.bell@gmail.com>

dev/fluids/Ammonia.json

@@ -290,5 +290,141 @@
       "rhoVtriple_units": "mol/m^3"
     }
   ], 
-  "NAME": "Ammonia"
+  "NAME": "Ammonia", 
+  "TRANSPORT": {
+    "viscosity": {
+      "BibTeX": "Fenghour-JPCRD-1995", 
+      "_note": "There is a mysterious factor of 100 missing from the Fenghour equation for the dilute gas.  From REFPROP fluid file: !=0.021357*SQRT(MW)*(unknown factor of 100)  [Chapman-Enskog term]", 
+      "dilute": {
+        "C": 2.1357e-06, 
+        "a": [
+          4.9931822, 
+          -0.61122364, 
+          0.0, 
+          0.18535124, 
+          -0.11160946
+        ], 
+        "molar_mass": 0.01703026, 
+        "molar_mass_units": "kg/mol", 
+        "t": [
+          0, 
+          1, 
+          2, 
+          3, 
+          4
+        ], 
+        "type": "collision_integral"
+      }, 
+      "epsilon_over_k": 386, 
+      "epsilon_over_k_units": "K", 
+      "higher_order": {
+        "T_reduce": 405.4, 
+        "T_reduce_units": "K", 
+        "a": [
+          4.005040600989671e-06, 
+          -1.4107915123955129e-05, 
+          3.4760743039321816e-05, 
+          4.631310990138071e-06, 
+          -3.937374461785061e-06, 
+          -1.200075068367531e-05, 
+          1.9284977991745303e-06
+        ], 
+        "d1": [
+          3, 
+          3, 
+          2, 
+          4, 
+          4, 
+          2, 
+          4
+        ], 
+        "d2": [
+          1
+        ], 
+        "f": [
+          0.0
+        ], 
+        "g": [
+          1
+        ], 
+        "gamma": [
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0, 
+          0
+        ], 
+        "h": [
+          0
+        ], 
+        "l": [
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1, 
+          1
+        ], 
+        "p": [
+          1
+        ], 
+        "q": [
+          0
+        ], 
+        "rhomolar_reduce": 13211.8, 
+        "rhomolar_reduce_units": "mol/m^3", 
+        "t1": [
+          0, 
+          1, 
+          2, 
+          2, 
+          3, 
+          4, 
+          4
+        ], 
+        "t2": [
+          0
+        ], 
+        "type": "modified_Batschinski_Hildebrand"
+      }, 
+      "initial_density": {
+        "b": [
+          -1.7999496, 
+          46.692621, 
+          -534.60794, 
+          3360.4074, 
+          -13019.164, 
+          33414.23, 
+          -58711.743, 
+          71426.686, 
+          -59834.012, 
+          33652.741, 
+          -12027.35, 
+          2434.8205, 
+          -208.07957
+        ], 
+        "t": [
+          -0.0, 
+          -0.5, 
+          -1.0, 
+          -1.5, 
+          -2.0, 
+          -2.5, 
+          -3.0, 
+          -3.5, 
+          -4.0, 
+          -4.5, 
+          -5.0, 
+          -5.5, 
+          -6.0
+        ], 
+        "type": "Rainwater-Friend"
+      }, 
+      "sigma_eta": 2.957e-10, 
+      "sigma_eta_units": "m"
+    }
+  }
 }

doc/notebooks/transport_correlation_conversion.ipynb

@@ -190,6 +190,8 @@
      "cell_type": "code",
      "collapsed": false,
      "input": [
+      "# rho is units of mol/L, so convert the density to mol/L (poorly documented in paper)\n",
+      "\n",
       "D = [(3,0,0.17366936e-8),\n",
       " (3,1,-0.64250359e-8),\n",
       " (2,2,2.19664285e-7),\n",
@@ -198,8 +200,8 @@
       " (2,4,-0.83651107e-7),\n",
       " (4,4,0.77012274e-10)]\n",
       " \n",
-      "e_k = 386\n",
-      "rhoc = 225\n",
+      "e_k = 386.0\n",
+      "rhoc = 225.0/17.03026\n",
       "Tc = 405.4\n",
       "v = []\n",
       "for i,j,dij in D:\n",
@@ -213,11 +215,11 @@
        "output_type": "stream",
        "stream": "stdout",
        "text": [
-        "[0.0197820255375, -0.0696829732996963, 0.010081649810518732, 0.38957318686639053, -0.33120114806158074, -0.003480574788614132, 0.16221994918696753]\n"
+        "[4.005040600989671e-06, -1.4107915123955129e-05, 3.4760743039321816e-05, 4.631310990138071e-06, -3.937374461785061e-06, -1.200075068367531e-05, 1.9284977991745303e-06]\n"
        ]
       }
      ],
-     "prompt_number": 2
+     "prompt_number": 6
     },
     {
      "cell_type": "heading",
@@ -519,14 +521,6 @@
       }
      ],
      "prompt_number": 8
-    },
-    {
-     "cell_type": "code",
-     "collapsed": false,
-     "input": [],
-     "language": "python",
-     "metadata": {},
-     "outputs": []
     }
    ],
    "metadata": {}

src/Tests/CoolProp-Tests.cpp

@@ -66,13 +66,13 @@ vel("Hydrogen", "T", 35, "Dmass", 100, "V", 5.6658815444656859e-005, 1e-3),
 vel("DimethylEther", "T", 253.146, "Dmass", 734.28, "V", 0.20444e-3, 3e-3),
 vel("DimethylEther", "T", 373.132, "Dmass", 613.78, "V", 0.09991e-3, 3e-3),
 
-//// From Fenghour, JPCRD, 1995
-//vel("Ammonia", "T", 200, "Dmolar", 3.9, "V", 6.95e-6, 1e-3),
-//vel("Ammonia", "T", 200, "Dmolar", 42754.4, "V", 507.28e-6, 1e-3),
-//vel("Ammonia", "T", 398, "Dmolar", 7044.7, "V", 17.67e-6, 1e-3),
-//vel("Ammonia", "T", 398, "Dmolar", 21066.7, "V", 43.95e-6, 1e-3),
+// From Fenghour, JPCRD, 1995
+vel("Ammonia", "T", 200, "Dmolar", 3.9, "V", 6.95e-6, 1e-3),
+vel("Ammonia", "T", 200, "Dmolar", 42754.4, "V", 507.28e-6, 1e-3),
+vel("Ammonia", "T", 398, "Dmolar", 7044.7, "V", 17.67e-6, 1e-3),
+vel("Ammonia", "T", 398, "Dmolar", 21066.7, "V", 43.95e-6, 1e-3),
 
-//// From Lemmon and Jacobsen, JPCRD, 2004
+// From Lemmon and Jacobsen, JPCRD, 2004
 vel("Nitrogen", "T", 100, "Dmolar", 1e-14, "V", 6.90349e-6, 1e-3),
 vel("Nitrogen", "T", 300, "Dmolar", 1e-14, "V", 17.8771e-6, 1e-3),
 vel("Nitrogen", "T", 100, "Dmolar", 25000, "V", 79.7418e-6, 1e-3),