diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermBrooksCorey.cpp b/MaterialLib/MPL/Properties/RelativePermeability/RelPermBrooksCorey.cpp
index 4865dd8fdbc0144c43bf5c421f60775d7515256c..49b3a9f4824e439ed4daba9d7ee74484f3912cf0 100644
--- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermBrooksCorey.cpp
+++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermBrooksCorey.cpp
@@ -41,7 +41,7 @@ PropertyDataType RelPermBrooksCorey::value(
 {
     /// here, an extra computation of saturation is forced, guaranteeing a
     /// correct value. In order to speed up the computing time, saturation could
-    /// be insertred into the primary variable array after it is computed in the
+    /// be inserted into the primary variable array after it is computed in the
     /// FEM assembly.
     auto const s_L = std::visit(
         [&variable_array, &pos, t, dt](auto&& scale) -> double {
@@ -83,7 +83,7 @@ PropertyDataType RelPermBrooksCorey::dValue(
            " derivatives with respect to liquid saturation only.");
     /// here, an extra computation of saturation is forced, guaranteeing a
     /// correct value. In order to speed up the computing time, saturation could
-    /// be insertred into the primary variable array after it is computed in the
+    /// be inserted into the primary variable array after it is computed in the
     /// FEM assembly.
     auto const s_L = std::visit(
         [&variable_array, &pos, t, dt](auto&& scale) -> double {
diff --git a/MaterialLib/MPL/Properties/RelativePermeability/RelPermLiakopoulos.cpp b/MaterialLib/MPL/Properties/RelativePermeability/RelPermLiakopoulos.cpp
index 7fab447cf45e41d1925a9929ae48a727768feadc..efb0545f7d1d04d74372a184fa2e97093a02b419 100644
--- a/MaterialLib/MPL/Properties/RelativePermeability/RelPermLiakopoulos.cpp
+++ b/MaterialLib/MPL/Properties/RelativePermeability/RelPermLiakopoulos.cpp
@@ -32,7 +32,7 @@ PropertyDataType RelPermLiakopoulos::value(
 {
     /// here, an extra computation of saturation is forced, guaranteeing a
     /// correct value. In order to speed up the computing time, saturation could
-    /// be insertred into the primary variable array after it is computed in the
+    /// be inserted into the primary variable array after it is computed in the
     /// FEM assembly.
     auto const s_L = std::visit(
         [&variable_array, &pos, t, dt](auto&& scale) -> double {
@@ -71,7 +71,7 @@ PropertyDataType RelPermLiakopoulos::dValue(
            " derivatives with respect to liquid saturation only.");
     /// here, an extra computation of saturation is forced, guaranteeing a
     /// correct value. In order to speed up the computing time, saturation could
-    /// be insertred into the primary variable array after it is computed in the
+    /// be inserted into the primary variable array after it is computed in the
     /// FEM assembly.
     auto const s_L = std::visit(
         [&variable_array, &pos, t, dt](auto&& scale) -> double {