Beispiel #1
0
Real
LevelSetAdvectionSUPG::computeQpJacobian()
{
  computeQpVelocity();
  Real tau = _current_elem->hmin() / (2 * _velocity.norm());
  return (tau * _velocity * _grad_test[_i][_qp]) * (_velocity * _grad_phi[_j][_qp]);
}
Real
LevelSetAdvection::computeQpJacobian()
{
  computeQpVelocity();
  return _test[_i][_qp] * (_velocity * _grad_phi[_j][_qp]);
}
Real
LevelSetAdvection::computeQpResidual()
{
  computeQpVelocity();
  return _test[_i][_qp] * (_velocity * _grad_u[_qp]);
}