Real
NSEnergyInviscidSpecifiedPressureBC::computeJacobianHelper(unsigned var_number)
{
// Velocity vector object
RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
// Normal component of velocity
Real un = vel * _normals[_qp];
// For specified pressure, term "C" is zero, see base class for details.
return qpJacobianTermA(var_number, _specified_pressure) + qpJacobianTermB(var_number, un);
}
Real
NSEnergyInviscidUnspecifiedBC::computeJacobianHelper(unsigned var_number)
{
// Velocity vector object
RealVectorValue vel(_u_vel[_qp], _v_vel[_qp], _w_vel[_qp]);
// Normal component of velocity
Real un = vel * _normals[_qp];
// When both u.n and pressure are unspecified, all 3 Jacobian terms apply.
// See base class for details.
return qpJacobianTermA(var_number, _pressure[_qp]) + qpJacobianTermB(var_number, un) +
qpJacobianTermC(var_number, un);
}
