diff --git a/src/thermophysicalModels/specie/thermo/thermo/thermo.H b/src/thermophysicalModels/specie/thermo/thermo/thermo.H
index f545b798..2ea2127d 100644
--- a/src/thermophysicalModels/specie/thermo/thermo/thermo.H
+++ b/src/thermophysicalModels/specie/thermo/thermo/thermo.H
@@ -333,6 +333,22 @@ public:
             ) const;
 
 
+        // Electric charge function
+
+
+            // Mole specific properties
+
+                //- Electric charge [C/kmol]
+                inline scalar qc() const;
+
+
+            // Mass specific properties
+
+                //- Electric charge [C/kg]
+                inline scalar QC() const;
+
+
+
         // I-O
 
             //- Write to Ostream
diff --git a/src/thermophysicalModels/specie/thermo/thermo/thermoI.H b/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
index 1de31499..58245095 100644
--- a/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
+++ b/src/thermophysicalModels/specie/thermo/thermo/thermoI.H
@@ -539,4 +539,20 @@ inline Foam::species::thermo<Thermo, Type> Foam::species::operator==
 }
 
 
+template<class Thermo, template<class> class Type>
+inline Foam::scalar
+Foam::species::thermo<Thermo, Type>::qc() const
+{
+    return this->z() * this->NA * this->e;
+}
+
+
+template<class Thermo, template<class> class Type>
+inline Foam::scalar
+Foam::species::thermo<Thermo, Type>::QC() const
+{
+    return this->z() * this->NA * this->e / this->W();
+}
+
+
 // ************************************************************************* //