During pressure equilibration in conjunction with JWL++, the bisection algorithm generates high pressures. When trying to intersect these pressures with the tabular hydrostat, large values of strain are generated causing the value of J = detF to go to zero and the value of rho_cur to go to infinity.
During pressure equilibration in conjunction with JWL++, the bisection algorithm generates high pressures. When trying to intersect these pressures with the tabular hydrostat, large values of strain are generated causing the value of J = detF to go to zero and the value of rho_cur to go to infinity.