40 #include "libmesh/quadrature.h"
41 #include "libmesh/fem_system.h"
81 this->_dim = system->get_mesh().mesh_dimension();
83 _temp_vars.init(system);
92 system->time_evolving(_temp_vars.T());
103 context.get_element_fe(_temp_vars.T())->get_JxW();
104 context.get_element_fe(_temp_vars.T())->get_phi();
105 context.get_element_fe(_temp_vars.T())->get_dphi();
106 context.get_element_fe(_temp_vars.T())->get_xyz();
108 context.get_side_fe(_temp_vars.T())->get_JxW();
109 context.get_side_fe(_temp_vars.T())->get_phi();
110 context.get_side_fe(_temp_vars.T())->get_dphi();
111 context.get_side_fe(_temp_vars.T())->get_xyz();
122 const unsigned int n_T_dofs = context.get_dof_indices(_temp_vars.T()).size();
128 const std::vector<libMesh::Real> &JxW =
129 context.get_element_fe(_temp_vars.T())->get_JxW();
133 const std::vector<std::vector<libMesh::RealGradient> >& T_gradphi =
134 context.get_element_fe(_temp_vars.T())->get_dphi();
142 libMesh::DenseSubMatrix<libMesh::Number> &KTT = context.get_elem_jacobian(_temp_vars.T(), _temp_vars.T());
144 libMesh::DenseSubVector<libMesh::Number> &FT = context.get_elem_residual(_temp_vars.T());
152 unsigned int n_qpoints = context.get_element_qrule().n_points();
154 for (
unsigned int qp=0; qp != n_qpoints; qp++)
157 libMesh::Gradient grad_T;
158 grad_T = context.interior_gradient(_temp_vars.T(), qp);
161 libMesh::Real _k_qp = this->_k(context, qp);
164 for (
unsigned int i=0; i != n_T_dofs; i++)
166 FT(i) += JxW[qp]*(-_k_qp*(T_gradphi[i][qp]*grad_T));
168 if (compute_jacobian)
170 for (
unsigned int j=0; j != n_T_dofs; j++)
175 KTT(i,j) += JxW[qp] * context.get_elem_solution_derivative() *
176 ( -_k_qp*(T_gradphi[i][qp]*T_gradphi[j][qp]) );
196 const std::vector<libMesh::Real> &JxW =
197 context.get_element_fe(_temp_vars.T())->get_JxW();
200 const std::vector<std::vector<libMesh::Real> >& phi =
201 context.get_element_fe(_temp_vars.T())->get_phi();
204 const unsigned int n_T_dofs = context.get_dof_indices(_temp_vars.T()).size();
207 libMesh::DenseSubVector<libMesh::Real> &F =
208 context.get_elem_residual(_temp_vars.T());
210 libMesh::DenseSubMatrix<libMesh::Real> &M =
211 context.get_elem_jacobian(_temp_vars.T(), _temp_vars.T());
213 unsigned int n_qpoints = context.get_element_qrule().n_points();
215 for (
unsigned int qp = 0; qp != n_qpoints; ++qp)
223 context.interior_rate(_temp_vars.T(), qp, T_dot);
225 for (
unsigned int i = 0; i != n_T_dofs; ++i)
227 F(i) -= JxW[qp]*(_rho*_Cp*T_dot*phi[i][qp] );
229 if( compute_jacobian )
231 for (
unsigned int j=0; j != n_T_dofs; j++)
235 context.get_elem_solution_rate_derivative()
236 * JxW[qp]*_rho*_Cp*phi[j][qp]*phi[i][qp] ;
249 (
const std::string & param_name,
254 _k.register_parameter(param_name, param_pointer);
virtual void mass_residual(bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
Mass matrix part(s) for element interiors. All boundary terms lie within the time_derivative part...
GRINS::ICHandlingBase * _ic_handler
void register_variables()
Physics abstract base class. Defines API for physics to be added to MultiphysicsSystem.
static std::string temperature_section()
INSTANTIATE_HEAT_TRANSFER_SUBCLASS(HeatConduction)
virtual void element_time_derivative(bool compute_jacobian, AssemblyContext &context, CachedValues &cache)
Time dependent part(s) of physics for element interiors.
static void check_and_register_variable(const std::string &var_name, const FEVariablesBase &variable)
First check if var_name is registered and then register.
virtual void register_parameter(const std::string ¶m_name, libMesh::ParameterMultiAccessor< libMesh::Number > ¶m_pointer) const
Each subclass will register its copy of an independent.
Base class for reading and handling initial conditions for physics classes.
static void read_density(const std::string &core_physics_name, const GetPot &input, ParameterUser ¶ms, libMesh::Real &rho)
Helper function to reading density from input.
static void read_specific_heat(const std::string &core_physics_name, const GetPot &input, ParameterUser ¶ms, libMesh::Real &cp)
Helper function to reading scalar specific heat from input.
Helper functions for parsing material properties.
virtual void init_context(AssemblyContext &context)
Initialize context for added physics variables.
virtual void set_time_evolving_vars(libMesh::FEMSystem *system)
Set which variables are time evolving.
static PhysicsName heat_conduction()
virtual void register_parameter(const std::string ¶m_name, libMesh::ParameterMultiAccessor< libMesh::Number > ¶m_pointer) const
Each subclass will register its copy of an independent.
virtual void init_variables(libMesh::FEMSystem *system)
Initialize variables for this physics.