diff --git a/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.cpp b/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.cpp
new file mode 100644
index 0000000000000000000000000000000000000000..bdd0fa4ec06a2d8c92bd76c005d22679db6cc4d5
--- /dev/null
+++ b/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.cpp
@@ -0,0 +1,28 @@
+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ * Created on April 20, 2020, 10:47 AM
+ */
+
+#include "CheckVanGenuchtenExponentRange.h"
+
+#include "BaseLib/Error.h"
+
+namespace MaterialPropertyLib
+{
+void checkVanGenuchtenExponentRange(const double m)
+{
+    if (m <= 0 || m >= 1)
+    {
+        OGS_FATAL(
+            "The exponent value m = {:e} of van Genuchten saturation "
+            "model, is out of its range of(0, 1) ",
+            m);
+    }
+}
+}  // namespace MaterialPropertyLib
diff --git a/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.h b/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.h
new file mode 100644
index 0000000000000000000000000000000000000000..e93f7691b43a62d985c0c4002888f2439436caf4
--- /dev/null
+++ b/MaterialLib/MPL/Utils/CheckVanGenuchtenExponentRange.h
@@ -0,0 +1,17 @@
+/**
+ * \file
+ * \copyright
+ * Copyright (c) 2012-2020, OpenGeoSys Community (http://www.opengeosys.org)
+ *            Distributed under a Modified BSD License.
+ *              See accompanying file LICENSE.txt or
+ *              http://www.opengeosys.org/project/license
+ *
+ * Created on April 20, 2020, 10:47 AM
+ */
+
+#pragma once
+
+namespace MaterialPropertyLib
+{
+void checkVanGenuchtenExponentRange(const double m);
+}  // namespace MaterialPropertyLib