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); }