https://mooseframework.inl.gov
Loading...
Searching...
No Matches
PorousFlowPreDis.h
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
10#pragma once
11
13#include "PorousFlowDictator.h"
14
23template <bool is_ad>
25{
26public:
28
30
31protected:
32 virtual GenericReal<is_ad> computeQpResidual() override;
33 virtual Real computeQpJacobian() override;
34 virtual Real computeQpOffDiagJacobian(unsigned int jvar) override;
35
37 Real computeQpJac(unsigned int pvar);
38
40 const std::vector<Real> _mineral_density;
41
44
46 const unsigned int _aq_ph;
47
50
53
56
59
62
64 const std::vector<Real> _stoichiometry;
65
67};
68
typename GenericMaterialPropertyStruct< T, is_ad >::type GenericMaterialProperty
Moose::GenericType< Real, is_ad > GenericReal
PorousFlowPreDisTempl< false > PorousFlowPreDis
PorousFlowPreDisTempl< true > ADPorousFlowPreDis
const InputParameters & parameters() const
This holds maps between the nonlinear variables used in a PorousFlow simulation and the variable numb...
Base class for PorousFlow kernels that use mass-lumped (nodal) material properties.
Kernel = sum (stoichiometry * density * porosity_old * saturation * reaction_rate) where the sum is o...
static InputParameters validParams()
virtual Real computeQpJacobian() override
const MaterialProperty< std::vector< std::vector< Real > > > *const _dreaction_rate_dvar
d(reaction rate)/d(PorousFlow variable) – null for AD path
const PorousFlowDictator & _dictator
PorousFlowDictator UserObject.
Real computeQpJac(unsigned int pvar)
Derivative of residual wrt PorousFlow variable pvar (non-AD path only)
const std::vector< Real > _stoichiometry
Stoichiometric coefficients.
const std::vector< Real > _mineral_density
Density of the mineral species.
const unsigned int _aq_ph
Aqueous phase number.
const MaterialProperty< std::vector< std::vector< Real > > > *const _dsaturation_dvar
d(saturation)/d(PorousFlow var) – null for AD path
virtual Real computeQpOffDiagJacobian(unsigned int jvar) override
const MaterialProperty< Real > & _porosity_old
Old value of porosity (always non-AD: old values carry no current-timestep derivatives)
virtual GenericReal< is_ad > computeQpResidual() override
const GenericMaterialProperty< std::vector< Real >, is_ad > & _reaction_rate
Reaction rate of the yielding the secondary species (AD or non-AD)
const GenericMaterialProperty< std::vector< Real >, is_ad > & _saturation
Saturation (AD or non-AD)