20 params.
addClassDescription(
"Conservative form of $\\nabla \\cdot \\vec{v} u$ (advection), using " 21 "the Flux Limited TVD scheme invented by Kuzmin and Turek");
23 "AdvectiveFluxCalculator UserObject");
36 mooseError(
"FluxLimitedTVDAdvection::computeQpResidual() called\n");
57 Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
58 for (
unsigned int i = 0; i <
_save_in.size(); i++)
79 std::vector<dof_id_type> idof_indices(
89 std::vector<dof_id_type> jdof_indices(derivs.size());
92 for (
const auto & node_j_deriv : derivs)
100 deriv_matrix(0,
j) = node_j_deriv.second / valence;
unsigned getValence(dof_id_type node_i) const
Returns the valence of the global node i Valence is the number of times the node is encountered in a ...
static InputParameters validParams()
void accumulateTaggedLocalResidual()
Real getFluxOut(dof_id_type node_i) const
Returns the flux out of lobal node id.
std::vector< MooseVariableFEBase *> _save_in
FluxLimitedTVDAdvection(const InputParameters ¶meters)
unsigned int number() const
virtual Real computeQpResidual() override
virtual void precalculateResidual()
void addJacobian(Assembly &assembly, const Residuals &residuals, const Indices &dof_indices, Real scaling_factor)
const AdvectiveFluxCalculatorBase & _fluo
The user object that computes Kuzmin and Turek's K_ij, R+ and R-, etc quantities. ...
virtual void computeResidual() override
virtual void precalculateJacobian()
Advection of the variable with velocity set in the AdvectiveFluxCalculator.
Base class to compute Advective fluxes.
const std::map< dof_id_type, Real > & getdFluxOutdu(dof_id_type node_i) const
Returns r where r[j] = d(flux out of global node i)/du(global node j) used in Jacobian computations...
unsigned int number() const
static InputParameters validParams()
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
DenseVector< Number > _local_re
void mooseError(Args &&... args) const
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
const Elem *const & _current_elem
registerMooseObject("PorousFlowApp", FluxLimitedTVDAdvection)
void prepareVectorTag(Assembly &assembly, unsigned int ivar)
void prepareMatrixTag(Assembly &assembly, unsigned int ivar, unsigned int jvar)
virtual void computeJacobian() override
void scalingFactor(const std::vector< Real > &factor)