void
AEFVUpwindInternalSideFlux::calcFlux(unsigned int /*iside*/,
dof_id_type /*ielem*/,
dof_id_type /*ineig*/,
const std::vector<Real> & uvec1,
const std::vector<Real> & uvec2,
const RealVectorValue & dwave,
std::vector<Real> & flux) const
{
mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
// assign the size of flux vector, e.g. = 1 for the advection equation
flux.resize(1);
// assume a constant velocity on the left
RealVectorValue uadv1(1.0, 1.0, 1.0);
// assume a constant velocity on the right
RealVectorValue uadv2(1.0, 1.0, 1.0);
// normal velocity on the left and right
Real vdon1 = uadv1 * dwave;
Real vdon2 = uadv2 * dwave;
// calculate the so-called a^plus and a^minus
Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
Real amins = 0.5 * (vdon2 - std::abs(vdon2));
// finally calculate the flux
flux[0] = aplus * uvec1[0] + amins * uvec2[0];
}
void
AEFVUpwindInternalSideFlux::calcJacobian(unsigned int /*iside*/,
dof_id_type /*ielem*/,
dof_id_type /*ineig*/,
const std::vector<Real> & libmesh_dbg_var(uvec1),
const std::vector<Real> & libmesh_dbg_var(uvec2),
const RealVectorValue & dwave,
DenseMatrix<Real> & jac1,
DenseMatrix<Real> & jac2) const
{
mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
mooseAssert(uvec2.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
// assign the size of Jacobian matrix, e.g. = (1, 1) for the advection equation
jac1.resize(1, 1);
jac2.resize(1, 1);
// assume a constant velocity on the left
RealVectorValue uadv1(1.0, 1.0, 1.0);
// assume a constant velocity on the right
RealVectorValue uadv2(1.0, 1.0, 1.0);
// normal velocity on the left and right
Real vdon1 = uadv1 * dwave;
Real vdon2 = uadv2 * dwave;
// calculate the so-called a^plus and a^minus
Real aplus = 0.5 * (vdon1 + std::abs(vdon1));
Real amins = 0.5 * (vdon2 - std::abs(vdon2));
// finally calculate the Jacobian matrix
jac1(0, 0) = aplus;
jac2(0, 0) = amins;
}
void
AEFVFreeOutflowBoundaryFlux::calcFlux(unsigned int /*iside*/,
dof_id_type /*ielem*/,
const std::vector<Real> & uvec1,
const RealVectorValue & dwave,
std::vector<Real> & flux) const
{
mooseAssert(uvec1.size() == 1, "Invalid size for uvec1. Must be single variable coupling.");
// assume the velocity vector is constant, e.g. = (1., 1., 1.)
RealVectorValue uadv1(1.0, 1.0, 1.0);
// assign the size of flux vector, e.g. = 1 for the advection equation
flux.resize(1);
// finally calculate the flux
flux[0] = (uadv1 * dwave) * uvec1[0];
}
