Should d2udt2 be set to 0 when Dirichlet boundary condition is applied?

[Wayne901]

Should d2udt2 be set to 0 when Dirichlet boundary condition is applied?

[Wayne901]

May someone give some hint on this problem? I'm still confused about it. The different result is not observed when t_final <= 4. I set d2udt2 to 0 after T_solver with

T_solver.Mult(z, d2udt2);
for (int i = 0; i < ess_tdof_list.Size(); i++)
  d2udt2[ess_tdof_list[i]] = 0.0;

[Wayne901]

I wrote another wave equation code based on the coupled first-order equations, like ex10.
The result at t=5 is same as the ex23 with d2udt2 set to 0.

[IdoAkkerman]

I have tried to address the issue in a small PR #4233. I added as step after solving for the acceleration where the acceleration is hardwired set to zero.

I think, this should technically not have been necessary, as the solve should have resulted in zero accelerations on those places. But apparently round-off and/or the approximate nature of the linear solves create for excessive slip of the solution.

[IdoAkkerman]

I have checked the BC values and after zeroing the accelerations the solution values at the BC were of the order e-105 which is effectively zero.