예제 #1
0
파일: DashpotBC.C 프로젝트: DarinReid/moose
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)
{}
예제 #2
0
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());
}
예제 #4
0
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());
}
예제 #5
0
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());
}
예제 #6
0
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);
}
예제 #7
0
DotCouplingKernel::DotCouplingKernel(const InputParameters & parameters)
  : Kernel(parameters), _v_dot(coupledDot("v")), _dv_dot_dv(coupledDotDu("v"))
{
}
예제 #8
0
stzplasticrateAux_chidot::stzplasticrateAux_chidot(const InputParameters & parameters)
  : AuxKernel(parameters),
    _coupled_chidot(coupledDot("couple_chidot"))

{
}
예제 #9
0
DotCouplingKernel::DotCouplingKernel(const std::string & name, InputParameters parameters) :
    Kernel(name, parameters),
    _v_dot(coupledDot("v")),
    _dv_dot_dv(coupledDotDu("v"))
{
}
예제 #10
0
// 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;
  }
예제 #11
0
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"))
{
}
예제 #12
0
DotCouplingAux::DotCouplingAux(const InputParameters & parameters)
  : AuxKernel(parameters), _v_dot(coupledDot("v"))
{
}
예제 #13
0
TestNewmarkTI::TestNewmarkTI(const InputParameters & parameters)
  : AuxKernel(parameters),
    _first(getParam<bool>("first")),
    _value(_first ? coupledDot("displacement") : coupledDotDot("displacement"))
{
}
예제 #14
0
CoupledImplicitEuler::CoupledImplicitEuler(const std::string & name, InputParameters parameters)
  :Kernel(name, parameters),
   _v_dot(coupledDot("v")),
   _dv_dot(coupledDotDu("v")),
   _v_var(coupled("v"))
{}
예제 #15
0
파일: combust.C 프로젝트: jgreen9797/raptor
combust::combust(const InputParameters & parameters) :
    Kernel(parameters),
    _coupled_val(coupledValue("coupled")),
    _dcoupled_val_dt(coupledDot("coupled")),
    _func(getFunction("function"))
{}