DashpotBC::DashpotBC(const std::string & name, InputParameters parameters) :
IntegratedBC(name, parameters),
_component(getParam<unsigned int>("component")),
_coefficient(getParam<Real>("coefficient")),
_disp_x_var(coupled("disp_x")),
_disp_y_var(isCoupled("disp_y") ? coupled("disp_y") : 0),
_disp_z_var(isCoupled("disp_z") ? coupled("disp_z") : 0),
_disp_x_dot(coupledDot("disp_x")),
_disp_y_dot(isCoupled("disp_y") ? coupledDot("disp_y") : _zero),
_disp_z_dot(isCoupled("disp_z") ? coupledDot("disp_z") : _zero)
{}
PyrolysisGasSource::PyrolysisGasSource(const InputParameters & parameters) :
Kernel(parameters),
_rho_dt(coupledDot("rho")),
_rhoDot_dRho(coupledDotDu("rho")),
_rho_num(coupled("rho"))
{
}
ElementIntegralVariablePostprocessor::ElementIntegralVariablePostprocessor(const InputParameters & parameters) :
ElementIntegralPostprocessor(parameters),
MooseVariableInterface(parameters, false),
_u(coupledValue("variable")),
_grad_u(coupledGradient("variable")),
_u_dot(coupledDot("variable"))
{
addMooseVariableDependency(mooseVariable());
}
ElementVariablePostprocessor::ElementVariablePostprocessor(const std::string & name, InputParameters parameters) :
ElementPostprocessor(name, parameters),
MooseVariableInterface(parameters, false),
_u(coupledValue("variable")),
_grad_u(coupledGradient("variable")),
_u_dot(coupledDot("variable")),
_qp(0)
{
addMooseVariableDependency(mooseVariable());
}
ElementVariablePostprocessor::ElementVariablePostprocessor(const InputParameters & parameters)
: ElementPostprocessor(parameters),
MooseVariableInterface<Real>(this,
false,
"variable",
Moose::VarKindType::VAR_ANY,
Moose::VarFieldType::VAR_FIELD_STANDARD),
_u(coupledValue("variable")),
_grad_u(coupledGradient("variable")),
_u_dot(_is_transient ? coupledDot("variable") : _zero),
_qp(0)
{
addMooseVariableDependency(mooseVariable());
}
OptionallyCoupledForce::OptionallyCoupledForce(const InputParameters & parameters)
: Kernel(parameters),
_v_var(coupled("v")),
_v(coupledValue("v")),
_grad_v(coupledGradient("v")),
_second_v(coupledSecond("v")),
_v_dot(_is_transient ? coupledDot("v") : _zero),
_v_dot_du(_is_transient ? coupledDotDu("v") : _zero),
_v_coupled(isCoupled("v"))
{
if (!_v_coupled && _v_var < 64)
mooseError("Something is wrong with the coupling system. It should be producing really huge "
"numbers for coupled('v') But instead it generated: ",
_v_var);
}
DotCouplingKernel::DotCouplingKernel(const InputParameters & parameters)
: Kernel(parameters), _v_dot(coupledDot("v")), _dv_dot_dv(coupledDotDu("v"))
{
}
stzplasticrateAux_chidot::stzplasticrateAux_chidot(const InputParameters & parameters)
: AuxKernel(parameters),
_coupled_chidot(coupledDot("couple_chidot"))
{
}
DotCouplingKernel::DotCouplingKernel(const std::string & name, InputParameters parameters) :
Kernel(name, parameters),
_v_dot(coupledDot("v")),
_dv_dot_dv(coupledDotDu("v"))
{
}
// DEPRECATED CONSTRUCTOR
NavierStokesMaterial::NavierStokesMaterial(const std::string & deprecated_name, InputParameters parameters)
:
Material(deprecated_name, parameters),
_mesh_dimension(_mesh.dimension()),
_grad_u(coupledGradient("u")),
_grad_v(coupledGradient("v")),
_grad_w(_mesh_dimension == 3 ? coupledGradient("w") : _grad_zero),
_viscous_stress_tensor(declareProperty<RealTensorValue>("viscous_stress_tensor")),
_thermal_conductivity(declareProperty<Real>("thermal_conductivity")),
// Declared here but _not_ calculated here
// (See e.g. derived class, bighorn/include/materials/FluidTC1.h)
_dynamic_viscosity(declareProperty<Real>("dynamic_viscosity")),
// The momentum components of the inviscid flux Jacobians.
_calA(declareProperty<std::vector<RealTensorValue> >("calA")),
// "Velocity column" matrices
_calC(declareProperty<std::vector<RealTensorValue> >("calC")),
// Energy equation inviscid flux matrices, "cal E_{kl}" in the notes.
_calE(declareProperty<std::vector<std::vector<RealTensorValue> > >("calE")),
// Parameter values read in from input file
_R(getParam<Real>("R")),
_gamma(getParam<Real>("gamma")),
_Pr(getParam<Real>("Pr")),
// Coupled solution values needed for computing SUPG stabilization terms
_u_vel(coupledValue("u")),
_v_vel(coupledValue("v")),
_w_vel(_mesh.dimension() == 3 ? coupledValue("w") : _zero),
_temperature(coupledValue("temperature")),
_enthalpy(coupledValue("enthalpy")),
// Coupled solution values
_rho(coupledValue("rho")),
_rho_u(coupledValue("rhou")),
_rho_v(coupledValue("rhov")),
_rho_w(_mesh.dimension() == 3 ? coupledValue("rhow") : _zero),
_rho_e(coupledValue("rhoe")),
// Time derivative values
_drho_dt(coupledDot("rho")),
_drhou_dt(coupledDot("rhou")),
_drhov_dt(coupledDot("rhov")),
_drhow_dt( _mesh.dimension() == 3 ? coupledDot("rhow") : _zero),
_drhoe_dt(coupledDot("rhoe")),
// Gradients
_grad_rho(coupledGradient("rho")),
_grad_rho_u(coupledGradient("rhou")),
_grad_rho_v(coupledGradient("rhov")),
_grad_rho_w( _mesh.dimension() == 3 ? coupledGradient("rhow") : _grad_zero),
_grad_rho_e(coupledGradient("rhoe")),
// Material properties for stabilization
_hsupg(declareProperty<Real>("hsupg")),
_tauc(declareProperty<Real>("tauc")),
_taum(declareProperty<Real>("taum")),
_taue(declareProperty<Real>("taue")),
_strong_residuals(declareProperty<std::vector<Real> >("strong_residuals"))
{
// Load the velocity gradients up into a single vector for convenience
_vel_grads.resize(3);
_vel_grads[0] = &_grad_u;
_vel_grads[1] = &_grad_v;
_vel_grads[2] = &_grad_w;
}
NavierStokesMaterial::NavierStokesMaterial(const InputParameters & parameters)
: Material(parameters),
_mesh_dimension(_mesh.dimension()),
_grad_u(coupledGradient(NS::velocity_x)),
_grad_v(_mesh_dimension >= 2 ? coupledGradient(NS::velocity_y) : _grad_zero),
_grad_w(_mesh_dimension == 3 ? coupledGradient(NS::velocity_z) : _grad_zero),
_viscous_stress_tensor(declareProperty<RealTensorValue>("viscous_stress_tensor")),
_thermal_conductivity(declareProperty<Real>("thermal_conductivity")),
// Declared here but _not_ calculated here
// (See e.g. derived class, bighorn/include/materials/FluidTC1.h)
_dynamic_viscosity(declareProperty<Real>("dynamic_viscosity")),
// The momentum components of the inviscid flux Jacobians.
_calA(declareProperty<std::vector<RealTensorValue>>("calA")),
// "Velocity column" matrices
_calC(declareProperty<std::vector<RealTensorValue>>("calC")),
// Energy equation inviscid flux matrices, "cal E_{kl}" in the notes.
_calE(declareProperty<std::vector<std::vector<RealTensorValue>>>("calE")),
_vel_grads({&_grad_u, &_grad_v, &_grad_w}),
// Coupled solution values needed for computing SUPG stabilization terms
_u_vel(coupledValue(NS::velocity_x)),
_v_vel(_mesh.dimension() >= 2 ? coupledValue(NS::velocity_y) : _zero),
_w_vel(_mesh.dimension() == 3 ? coupledValue(NS::velocity_z) : _zero),
_temperature(coupledValue(NS::temperature)),
_enthalpy(coupledValue(NS::enthalpy)),
// Coupled solution values
_rho(coupledValue(NS::density)),
_rho_u(coupledValue(NS::momentum_x)),
_rho_v(_mesh.dimension() >= 2 ? coupledValue(NS::momentum_y) : _zero),
_rho_w(_mesh.dimension() == 3 ? coupledValue(NS::momentum_z) : _zero),
_rho_E(coupledValue(NS::total_energy)),
// Time derivative values
_drho_dt(coupledDot(NS::density)),
_drhou_dt(coupledDot(NS::momentum_x)),
_drhov_dt(_mesh.dimension() >= 2 ? coupledDot(NS::momentum_y) : _zero),
_drhow_dt(_mesh.dimension() == 3 ? coupledDot(NS::momentum_z) : _zero),
_drhoE_dt(coupledDot(NS::total_energy)),
// Gradients
_grad_rho(coupledGradient(NS::density)),
_grad_rho_u(coupledGradient(NS::momentum_x)),
_grad_rho_v(_mesh.dimension() >= 2 ? coupledGradient(NS::momentum_y) : _grad_zero),
_grad_rho_w(_mesh.dimension() == 3 ? coupledGradient(NS::momentum_z) : _grad_zero),
_grad_rho_E(coupledGradient(NS::total_energy)),
// Material properties for stabilization
_hsupg(declareProperty<Real>("hsupg")),
_tauc(declareProperty<Real>("tauc")),
_taum(declareProperty<Real>("taum")),
_taue(declareProperty<Real>("taue")),
_strong_residuals(declareProperty<std::vector<Real>>("strong_residuals")),
_fp(getUserObject<IdealGasFluidProperties>("fluid_properties"))
{
}
DotCouplingAux::DotCouplingAux(const InputParameters & parameters)
: AuxKernel(parameters), _v_dot(coupledDot("v"))
{
}
TestNewmarkTI::TestNewmarkTI(const InputParameters & parameters)
: AuxKernel(parameters),
_first(getParam<bool>("first")),
_value(_first ? coupledDot("displacement") : coupledDotDot("displacement"))
{
}
CoupledImplicitEuler::CoupledImplicitEuler(const std::string & name, InputParameters parameters)
:Kernel(name, parameters),
_v_dot(coupledDot("v")),
_dv_dot(coupledDotDu("v")),
_v_var(coupled("v"))
{}
combust::combust(const InputParameters & parameters) :
Kernel(parameters),
_coupled_val(coupledValue("coupled")),
_dcoupled_val_dt(coupledDot("coupled")),
_func(getFunction("function"))
{}
