diff --git a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
index 8e09f59a2460bc396ae0ec7e34cce233f84eac64..63a60f7737186c674ca6bd0c29ff4be954365929 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-impl.h
@@ -87,7 +87,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::getFlux(
     const double mu =
         _material_properties.getViscosity(pressure, _reference_temperature);
 
-    Eigen::Vector3d flux;
+    Eigen::Vector3d flux(0.0, 0.0, 0.0);
     flux.head<GlobalDim>() =
         -permeability / mu * shape_matrices.dNdx *
         Eigen::Map<const NodalVectorType>(local_x.data(), local_x.size());