diff --git a/Source/Constraint/PbdConstraints/imstkPbdConstraint.h b/Source/Constraint/PbdConstraints/imstkPbdConstraint.h
index 4fe9b6a5a7d5dbdf3a5096dd199b2191c24eebe3..ccb0dfc39cafcaa7ab6817a43445e4b191550ed6 100644
--- a/Source/Constraint/PbdConstraints/imstkPbdConstraint.h
+++ b/Source/Constraint/PbdConstraints/imstkPbdConstraint.h
@@ -83,6 +83,6 @@ public:
 
 protected:
     std::vector<size_t> m_vertexIds;   ///> index of points for the constraint
-    double m_epsilon = 1.0e-6;         ///> Tolerance used for the costraints
+    double m_epsilon = 1.0e-16;        ///> Tolerance used for the costraints
 };
 }
diff --git a/Source/DynamicalModels/ObjectModels/imstkPbdModel.h b/Source/DynamicalModels/ObjectModels/imstkPbdModel.h
index 77aa54a2e37196d55b7a2a854b86f33564b24616..7c491b930cd80a69041dc2a76f0d05bcaef85b59 100644
--- a/Source/DynamicalModels/ObjectModels/imstkPbdModel.h
+++ b/Source/DynamicalModels/ObjectModels/imstkPbdModel.h
@@ -42,20 +42,20 @@ struct PBDModelConfig
     double m_viscousDampingCoeff = 0.01;     ///> Viscous damping coefficient [0, 1]
 
     double m_contactStiffness = 1.0;         ///> Contact stiffness for collisions
-    double m_proximity;                      ///> Proximity for collisions
+    double m_proximity = 0.1;                ///> Proximity for collisions
 
-    unsigned int m_maxIter;                 ///> Max. pbd iterations
-    double m_dt;                            ///> Time step size
-    double m_DefaultDt;                     ///> Default Time step size
+    unsigned int m_maxIter = 10;                 ///> Max. pbd iterations
+    double m_dt = 0.01;                          ///> Time step size
+    double m_DefaultDt = 0.01;                   ///> Default Time step size
 
     std::vector<std::size_t> m_fixedNodeIds; ///> Nodal IDs of the nodes that are fixed
-    Vec3r m_gravity;                         ///> Gravity
+    Vec3r m_gravity = Vec3r(0, -9.81, 0);                         ///> Gravity
 
-    double m_mu;           ///> Lame constant, if constraint type is FEM
-    double m_lambda;       ///> Lame constant, if constraint type is FEM
+    double m_mu = 0;           ///> Lame constant, if constraint type is FEM
+    double m_lambda = 0;       ///> Lame constant, if constraint type is FEM
 
-    double m_YoungModulus; ///> FEM parameter, if constraint type is FEM
-    double m_PoissonRatio; ///> FEM parameter, if constraint type is FEM
+    double m_YoungModulus = 1000; ///> FEM parameter, if constraint type is FEM
+    double m_PoissonRatio = 0.2;  ///> FEM parameter, if constraint type is FEM
 
     std::vector<std::pair<PbdConstraint::Type, double>>    m_RegularConstraints; ///> Constraints except FEM
     std::vector<std::pair<PbdConstraint::Type,