https://mooseframework.inl.gov
PorousFlowAdvectiveFluxCalculatorUnsaturatedMultiComponent.C
Go to the documentation of this file.
1 //* This file is part of the MOOSE framework
2 //* https://mooseframework.inl.gov
3 //*
4 //* All rights reserved, see COPYRIGHT for full restrictions
5 //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
6 //*
7 //* Licensed under LGPL 2.1, please see LICENSE for details
8 //* https://www.gnu.org/licenses/lgpl-2.1.html
9 
11 
13 
16 {
18  params.addClassDescription(
19  "Computes the advective flux of fluid of given phase and component. Hence this UserObject "
20  "is relevant to multi-phase, multi-component situations. Explicitly, the UserObject "
21  "computes (mass_fraction * density * relative_permeability / viscosity) * (- permeability * "
22  "(grad(P) - density * gravity)), using the Kuzmin-Turek FEM-TVD multidimensional "
23  "stabilization scheme");
24  return params;
25 }
26 
30  _relative_permeability(
31  getMaterialProperty<std::vector<Real>>("PorousFlow_relative_permeability_nodal")),
32  _drelative_permeability_dvar(getMaterialProperty<std::vector<std::vector<Real>>>(
33  "dPorousFlow_relative_permeability_nodal_dvar"))
34 {
35 }
36 
37 Real
39 {
40  return _relative_permeability[i][_phase] *
42 }
43 
44 Real
46  unsigned pvar) const
47 {
50  du += _relative_permeability[i][_phase] *
52  return du;
53 }
virtual Real computeU(unsigned i) const override
Computes the value of u at the local node id of the current element (_current_elem) ...
const MaterialProperty< std::vector< Real > > & _relative_permeability
Relative permeability of each phase.
const MaterialProperty< std::vector< std::vector< Real > > > & _drelative_permeability_dvar
Derivative of relative permeability of each phase wrt PorousFlow variables.
registerMooseObject("PorousFlowApp", PorousFlowAdvectiveFluxCalculatorUnsaturatedMultiComponent)
Computes the advective flux of fluid of given phase and fluid component.
virtual Real computeU(unsigned i) const override
Computes the value of u at the local node id of the current element (_current_elem) ...
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
virtual Real computedU_dvar(unsigned i, unsigned pvar) const override
Compute d(u)/d(porous_flow_variable)
void addClassDescription(const std::string &doc_string)
virtual Real computedU_dvar(unsigned i, unsigned pvar) const override
Compute d(u)/d(porous_flow_variable)