added binary diffusion calc - uses molecular weight of binary pair

2025-11-28 03:28:01 +00:00 · 2009-09-01 19:35:59 +01:00
parent ab9317a589
commit 49faec23f9
62 changed files with 287 additions and 0 deletions
--- a/src/thermophysicalModels/liquids/Ar/Ar.H
+++ b/src/thermophysicalModels/liquids/Ar/Ar.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/Ar/ArI.H
+++ b/src/thermophysicalModels/liquids/Ar/ArI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::Ar::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::Ar::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C10H22/C10H22.H
+++ b/src/thermophysicalModels/liquids/C10H22/C10H22.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C10H22/C10H22I.H
+++ b/src/thermophysicalModels/liquids/C10H22/C10H22I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C10H22::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C10H22::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C12H26/C12H26.H
+++ b/src/thermophysicalModels/liquids/C12H26/C12H26.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C12H26/C12H26I.H
+++ b/src/thermophysicalModels/liquids/C12H26/C12H26I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C12H26::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C12H26::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C13H28/C13H28.H
+++ b/src/thermophysicalModels/liquids/C13H28/C13H28.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C13H28/C13H28I.H
+++ b/src/thermophysicalModels/liquids/C13H28/C13H28I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C13H28::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C13H28::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C14H30/C14H30.H
+++ b/src/thermophysicalModels/liquids/C14H30/C14H30.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C14H30/C14H30I.H
+++ b/src/thermophysicalModels/liquids/C14H30/C14H30I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C14H30::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C14H30::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C16H34/C16H34.H
+++ b/src/thermophysicalModels/liquids/C16H34/C16H34.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C16H34/C16H34I.H
+++ b/src/thermophysicalModels/liquids/C16H34/C16H34I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C16H34::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C16H34::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C2H5OH/C2H5OH.H
+++ b/src/thermophysicalModels/liquids/C2H5OH/C2H5OH.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C2H5OH/C2H5OHI.H
+++ b/src/thermophysicalModels/liquids/C2H5OH/C2H5OHI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C2H5OH::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C2H5OH::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C2H6/C2H6.H
+++ b/src/thermophysicalModels/liquids/C2H6/C2H6.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C2H6/C2H6I.H
+++ b/src/thermophysicalModels/liquids/C2H6/C2H6I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C2H6::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C2H6::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C2H6O/C2H6O.H
+++ b/src/thermophysicalModels/liquids/C2H6O/C2H6O.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C2H6O/C2H6OI.H
+++ b/src/thermophysicalModels/liquids/C2H6O/C2H6OI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C2H6O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C2H6O::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C3H6O/C3H6O.H
+++ b/src/thermophysicalModels/liquids/C3H6O/C3H6O.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C3H6O/C3H6OI.H
+++ b/src/thermophysicalModels/liquids/C3H6O/C3H6OI.H
@ -101,4 +101,10 @@ inline Foam::scalar Foam::C3H6O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C3H6O::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C3H8/C3H8.H
+++ b/src/thermophysicalModels/liquids/C3H8/C3H8.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C3H8/C3H8I.H
+++ b/src/thermophysicalModels/liquids/C3H8/C3H8I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C3H8::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C3H8::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C4H10O/C4H10O.H
+++ b/src/thermophysicalModels/liquids/C4H10O/C4H10O.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C4H10O/C4H10OI.H
+++ b/src/thermophysicalModels/liquids/C4H10O/C4H10OI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C4H10O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C4H10O::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C6H14/C6H14.H
+++ b/src/thermophysicalModels/liquids/C6H14/C6H14.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C6H14/C6H14I.H
+++ b/src/thermophysicalModels/liquids/C6H14/C6H14I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C6H14::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C6H14::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C6H6/C6H6.H
+++ b/src/thermophysicalModels/liquids/C6H6/C6H6.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C6H6/C6H6I.H
+++ b/src/thermophysicalModels/liquids/C6H6/C6H6I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C6H6::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C6H6::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C7H16/C7H16.H
+++ b/src/thermophysicalModels/liquids/C7H16/C7H16.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C7H16/C7H16I.H
+++ b/src/thermophysicalModels/liquids/C7H16/C7H16I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C7H16::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C7H16::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C7H8/C7H8.H
+++ b/src/thermophysicalModels/liquids/C7H8/C7H8.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C7H8/C7H8I.H
+++ b/src/thermophysicalModels/liquids/C7H8/C7H8I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C7H8::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C7H8::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C8H10/C8H10.H
+++ b/src/thermophysicalModels/liquids/C8H10/C8H10.H
@ -152,6 +152,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C8H10/C8H10I.H
+++ b/src/thermophysicalModels/liquids/C8H10/C8H10I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C8H10::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C8H10::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C8H18/C8H18.H
+++ b/src/thermophysicalModels/liquids/C8H18/C8H18.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C8H18/C8H18I.H
+++ b/src/thermophysicalModels/liquids/C8H18/C8H18I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C8H18::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C8H18::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/C9H20/C9H20.H
+++ b/src/thermophysicalModels/liquids/C9H20/C9H20.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/C9H20/C9H20I.H
+++ b/src/thermophysicalModels/liquids/C9H20/C9H20I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::C9H20::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::C9H20::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/CH3OH/CH3OH.H
+++ b/src/thermophysicalModels/liquids/CH3OH/CH3OH.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/CH3OH/CH3OHI.H
+++ b/src/thermophysicalModels/liquids/CH3OH/CH3OHI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::CH3OH::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::CH3OH::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/CH4N2O/CH4N2O.H
+++ b/src/thermophysicalModels/liquids/CH4N2O/CH4N2O.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/CH4N2O/CH4N2OI.H
+++ b/src/thermophysicalModels/liquids/CH4N2O/CH4N2OI.H
@ -102,5 +102,11 @@ inline Foam::scalar Foam::CH4N2O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::CH4N2O::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //

--- a/src/thermophysicalModels/liquids/H2O/H2O.H
+++ b/src/thermophysicalModels/liquids/H2O/H2O.H
@ -152,6 +152,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/H2O/H2OI.H
+++ b/src/thermophysicalModels/liquids/H2O/H2OI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::H2O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::H2O::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/IC8H18/IC8H18.H
+++ b/src/thermophysicalModels/liquids/IC8H18/IC8H18.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/IC8H18/IC8H18I.H
+++ b/src/thermophysicalModels/liquids/IC8H18/IC8H18I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::IC8H18::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::IC8H18::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/IDEA/IDEA.H
+++ b/src/thermophysicalModels/liquids/IDEA/IDEA.H
@ -175,6 +175,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/IDEA/IDEAI.H
+++ b/src/thermophysicalModels/liquids/IDEA/IDEAI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::IDEA::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::IDEA::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/MB/MB.H
+++ b/src/thermophysicalModels/liquids/MB/MB.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/MB/MBI.H
+++ b/src/thermophysicalModels/liquids/MB/MBI.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::MB::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::MB::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/N2/N2.H
+++ b/src/thermophysicalModels/liquids/N2/N2.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/N2/N2I.H
+++ b/src/thermophysicalModels/liquids/N2/N2I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::N2::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::N2::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/aC10H7CH3/aC10H7CH3.H
+++ b/src/thermophysicalModels/liquids/aC10H7CH3/aC10H7CH3.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/aC10H7CH3/aC10H7CH3I.H
+++ b/src/thermophysicalModels/liquids/aC10H7CH3/aC10H7CH3I.H
@ -102,4 +102,12 @@ inline Foam::scalar Foam::aC10H7CH3::D(scalar p, scalar T) const
 }


+
+
+inline Foam::scalar Foam::aC10H7CH3::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/bC10H7CH3/bC10H7CH3.H
+++ b/src/thermophysicalModels/liquids/bC10H7CH3/bC10H7CH3.H
@ -153,6 +153,9 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/bC10H7CH3/bC10H7CH3I.H
+++ b/src/thermophysicalModels/liquids/bC10H7CH3/bC10H7CH3I.H
@ -102,4 +102,10 @@ inline Foam::scalar Foam::bC10H7CH3::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::bC10H7CH3::D(scalar p, scalar T, scalar Wb) const
+{
+    return D_.f(p, T, Wb);
+}
+
+
 // ************************************************************************* //
--- a/src/thermophysicalModels/liquids/iC3H8O/iC3H8O.H
+++ b/src/thermophysicalModels/liquids/iC3H8O/iC3H8O.H
@ -151,6 +151,10 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        // Note: behaves the same as D(p, T)
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/iC3H8O/iC3H8OI.H
+++ b/src/thermophysicalModels/liquids/iC3H8O/iC3H8OI.H
@ -102,5 +102,12 @@ inline Foam::scalar Foam::iC3H8O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::iC3H8O::D(scalar p, scalar T, scalar) const
+{
+    // Use 'standard' NSRDS D(p, T) function
+    return D(p, T);
+}
+
+
 // ************************************************************************* //

--- a/src/thermophysicalModels/liquids/liquid/liquid.H
+++ b/src/thermophysicalModels/liquids/liquid/liquid.H
@ -243,6 +243,9 @@ public:
            //- Vapour diffussivity [m2/s]
            virtual scalar D(scalar p, scalar T) const = 0;

+            //- Vapour diffussivity [m2/s] with specified binary pair
+            virtual scalar D(scalar p, scalar T, scalar Wb) const = 0;
+

    // I-O

--- a/src/thermophysicalModels/liquids/nC3H8O/nC3H8O.H
+++ b/src/thermophysicalModels/liquids/nC3H8O/nC3H8O.H
@ -151,6 +151,10 @@ public:
        //- Vapour diffussivity [m2/s]
        inline scalar D(scalar p, scalar T) const;

+        //- Vapour diffussivity [m2/s] with specified binary pair
+        // Note: behaves the same as D(p, T)
+        inline scalar D(scalar p, scalar T, scalar Wb) const;
+

    // I-O

--- a/src/thermophysicalModels/liquids/nC3H8O/nC3H8OI.H
+++ b/src/thermophysicalModels/liquids/nC3H8O/nC3H8OI.H
@ -102,5 +102,12 @@ inline Foam::scalar Foam::nC3H8O::D(scalar p, scalar T) const
 }


+inline Foam::scalar Foam::nC3H8O::D(scalar p, scalar T, scalar) const
+{
+    // Use 'standard' NSRDS D(p, T) function
+    return D(p, T);
+}
+
+
 // ************************************************************************* //

--- a/src/thermophysicalModels/thermophysicalFunctions/APIfunctions/APIdiffCoefFunc/APIdiffCoefFunc.H
+++ b/src/thermophysicalModels/thermophysicalFunctions/APIfunctions/APIdiffCoefFunc/APIdiffCoefFunc.H
@ -103,6 +103,14 @@ public:
            return 3.6059e-3*(pow(1.8*T, 1.75))*alpha_/(p*beta_);
        }

+        //- API vapour mass diffusivity function using properties from
+        //  construction - with specified binary pair
+        scalar f(scalar p, scalar T, scalar  Wa) const
+        {
+            const scalar alphaBinary = sqrt(1/wf_ + 1/Wa);
+            return 3.6059e-3*(pow(1.8*T, 1.75))*alphaBinary/(p*beta_);
+        }
+
        //- Write the function coefficients
        void writeData(Ostream& os) const
        {