MAST
heat_conduction_elem_base.h
Go to the documentation of this file.
1 /*
2  * MAST: Multidisciplinary-design Adaptation and Sensitivity Toolkit
3  * Copyright (C) 2013-2019 Manav Bhatia
4  *
5  * This library is free software; you can redistribute it and/or
6  * modify it under the terms of the GNU Lesser General Public
7  * License as published by the Free Software Foundation; either
8  * version 2.1 of the License, or (at your option) any later version.
9  *
10  * This library is distributed in the hope that it will be useful,
11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13  * Lesser General Public License for more details.
14  *
15  * You should have received a copy of the GNU Lesser General Public
16  * License along with this library; if not, write to the Free Software
17  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
18  */
19 
20 
21 #ifndef __mast__heat_conduction_elem_base__
22 #define __mast__heat_conduction_elem_base__
23 
24 // MAST includes
25 #include "base/elem_base.h"
26 
27 namespace MAST {
28 
29  // Forward declerations
30  class ElementPropertyCardBase;
31  class BoundaryConditionBase;
32  class FEMOperatorMatrix;
33  template <typename ValType> class FieldFunction;
34 
35 
53  class HeatConductionElementBase:
54  public MAST::ElementBase
55  {
56  public:
60  HeatConductionElementBase(MAST::SystemInitialization& sys,
61  MAST::AssemblyBase& assembly,
62  const MAST::GeomElem& elem,
63  const MAST::ElementPropertyCardBase& p);
64 
65 
66  virtual ~HeatConductionElementBase();
67 
68 
72  const MAST::ElementPropertyCardBase&
73  elem_property() {
74  return _property;
75  }
76 
77 
81  virtual void
82  internal_residual (bool request_jacobian,
83  RealVectorX& f,
84  RealMatrixX& jac);
85 
86 
90  virtual void
91  velocity_residual (bool request_jacobian,
92  RealVectorX& f,
93  RealMatrixX& jac_xdot,
94  RealMatrixX& jac);
95 
99  void
100  side_external_residual (bool request_jacobian,
101  RealVectorX& f,
102  RealMatrixX& jac,
103  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
104 
108  void
109  volume_external_residual (bool request_jacobian,
110  RealVectorX& f,
111  RealMatrixX& jac,
112  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
113 
117  virtual void
118  internal_residual_sensitivity (const MAST::FunctionBase& p,
119  RealVectorX& f);
120 
124  virtual void
125  internal_residual_boundary_velocity (const MAST::FunctionBase& p,
126  RealVectorX& f,
127  const unsigned int s,
128  const MAST::FieldFunction<RealVectorX>& vel_f);
129 
133  virtual void
134  velocity_residual_sensitivity (const MAST::FunctionBase& p,
135  RealVectorX& f);
136 
140  virtual void
141  velocity_residual_boundary_velocity (const MAST::FunctionBase& p,
142  RealVectorX& f,
143  const unsigned int s,
144  const MAST::FieldFunction<RealVectorX>& vel_f);
145 
149  void
150  side_external_residual_sensitivity (const MAST::FunctionBase& p,
151  RealVectorX& f,
152  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
153 
157  void
158  volume_external_residual_sensitivity (const MAST::FunctionBase& p,
159  RealVectorX& f,
160  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
161 
165  void
166  volume_external_residual_boundary_velocity(const MAST::FunctionBase& p,
167  RealVectorX& f,
168  const unsigned int s,
169  const MAST::FieldFunction<RealVectorX>& vel_f,
170  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
171 
172  protected:
173 
174 
175 
180  virtual void surface_flux_residual(bool request_jacobian,
181  RealVectorX& f,
182  RealMatrixX& jac,
183  const unsigned int s,
184  MAST::BoundaryConditionBase& bc);
185 
191  virtual void surface_flux_residual(bool request_jacobian,
192  RealVectorX& f,
193  RealMatrixX& jac,
194  MAST::BoundaryConditionBase& bc);
195 
200  virtual void surface_flux_residual_sensitivity(const MAST::FunctionBase& p,
201  RealVectorX& f,
202  const unsigned int s,
203  MAST::BoundaryConditionBase& bc);
204 
210  virtual void surface_flux_residual_sensitivity(const MAST::FunctionBase& p,
211  RealVectorX& f,
212  MAST::BoundaryConditionBase& bc);
213 
219  virtual void surface_flux_boundary_velocity(const MAST::FunctionBase& p,
220  RealVectorX& f,
221  const unsigned int s,
222  const MAST::FieldFunction<RealVectorX>& vel_f,
223  MAST::BoundaryConditionBase& bc);
224 
225 
230  virtual void surface_convection_residual(bool request_jacobian,
231  RealVectorX& f,
232  RealMatrixX& jac,
233  const unsigned int s,
234  MAST::BoundaryConditionBase& bc);
235 
236 
242  virtual void surface_convection_residual(bool request_jacobian,
243  RealVectorX& f,
244  RealMatrixX& jac,
245  MAST::BoundaryConditionBase& bc);
246 
247 
252  virtual void
253  surface_convection_residual_sensitivity(const MAST::FunctionBase& p,
254  RealVectorX& f,
255  const unsigned int s,
256  MAST::BoundaryConditionBase& bc);
257 
263  virtual void
264  surface_convection_residual_sensitivity(const MAST::FunctionBase& p,
265  RealVectorX& f,
266  MAST::BoundaryConditionBase& bc);
267 
273  virtual void
274  surface_convection_boundary_velocity(const MAST::FunctionBase& p,
275  RealVectorX& f,
276  const unsigned int s,
277  const MAST::FieldFunction<RealVectorX>& vel_f,
278  MAST::BoundaryConditionBase& bc);
279 
284  virtual void surface_radiation_residual(bool request_jacobian,
285  RealVectorX& f,
286  RealMatrixX& jac,
287  const unsigned int s,
288  MAST::BoundaryConditionBase& bc);
289 
290 
295  virtual void surface_radiation_residual(bool request_jacobian,
296  RealVectorX& f,
297  RealMatrixX& jac,
298  MAST::BoundaryConditionBase& bc);
299 
300 
305  virtual void
306  surface_radiation_residual_sensitivity(const MAST::FunctionBase& p,
307  RealVectorX& f,
308  const unsigned int s,
309  MAST::BoundaryConditionBase& bc);
310 
311 
316  virtual void
317  surface_radiation_residual_sensitivity(const MAST::FunctionBase& p,
318  RealVectorX& f,
319  MAST::BoundaryConditionBase& bc);
320 
326  virtual void
327  surface_radiation_boundary_velocity(const MAST::FunctionBase& p,
328  RealVectorX& f,
329  const unsigned int s,
330  const MAST::FieldFunction<RealVectorX>& vel_f,
331  MAST::BoundaryConditionBase& bc);
332 
337  virtual void volume_heat_source_residual(bool request_jacobian,
338  RealVectorX& f,
339  RealMatrixX& jac,
340  MAST::BoundaryConditionBase& bc);
341 
342 
347  virtual void
348  volume_heat_source_residual_sensitivity(const MAST::FunctionBase& p,
349  RealVectorX& f,
350  MAST::BoundaryConditionBase& bc);
351 
356  virtual void
357  volume_heat_source_boundary_velocity(const MAST::FunctionBase& p,
358  RealVectorX& f,
359  const unsigned int s,
360  const MAST::FieldFunction<RealVectorX>& vel_f,
361  MAST::BoundaryConditionBase& bc);
362 
368  void _initialize_mass_fem_operator(const unsigned int qp,
369  const MAST::FEBase& fe,
370  MAST::FEMOperatorMatrix& Bmat);
371 
372 
380  void _initialize_fem_gradient_operator(const unsigned int qp,
381  const unsigned int dim,
382  const MAST::FEBase& fe,
383  std::vector<MAST::FEMOperatorMatrix>& dBmat);
384 
388  const MAST::ElementPropertyCardBase& _property;
389  };
390 
391 }
392 
393 
394 #endif // __mast__heat_conduction_elem_base__
MAST::HeatConductionElementBase::surface_convection_residual_sensitivity
virtual void surface_convection_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector sensitivity and Jacobian due to surface convection.
Definition: heat_conduction_elem_base.cpp:1144
MAST::HeatConductionElementBase::_initialize_mass_fem_operator
void _initialize_mass_fem_operator(const unsigned int qp, const MAST::FEBase &fe, MAST::FEMOperatorMatrix &Bmat)
When mass = false, initializes the FEM operator matrix to the shape functions as .
Definition: heat_conduction_elem_base.cpp:1759
MAST::HeatConductionElementBase::side_external_residual
void side_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase * > &bc)
side external force contribution to system residual
Definition: heat_conduction_elem_base.cpp:230
MAST::HeatConductionElementBase::HeatConductionElementBase
HeatConductionElementBase(MAST::SystemInitialization &sys, MAST::AssemblyBase &assembly, const MAST::GeomElem &elem, const MAST::ElementPropertyCardBase &p)
Constructor.
Definition: heat_conduction_elem_base.cpp:36
MAST::HeatConductionElementBase::volume_heat_source_residual_sensitivity
virtual void volume_heat_source_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to volume heat source.
Definition: heat_conduction_elem_base.cpp:1661
MAST::HeatConductionElementBase::internal_residual
virtual void internal_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac)
internal force contribution to system residual
Definition: heat_conduction_elem_base.cpp:55
MAST::HeatConductionElementBase::surface_flux_residual_sensitivity
virtual void surface_flux_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector sensitivity and Jacobian due to surface flux on element side s...
Definition: heat_conduction_elem_base.cpp:892
MAST::HeatConductionElementBase::surface_flux_residual
virtual void surface_flux_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface flux on element side s.
Definition: heat_conduction_elem_base.cpp:812
MAST::HeatConductionElementBase::side_external_residual_sensitivity
void side_external_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, std::multimap< libMesh::boundary_id_type, MAST::BoundaryConditionBase * > &bc)
sensitivity of the side external force contribution to system residual
Definition: heat_conduction_elem_base.cpp:349
MAST::ElementBase::elem
const MAST::GeomElem & elem() const
Definition: elem_base.h:117
MAST::HeatConductionElementBase::volume_heat_source_boundary_velocity
virtual void volume_heat_source_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to volume heat source on element volumetric domain...
Definition: heat_conduction_elem_base.cpp:1699
MAST::HeatConductionElementBase::volume_external_residual_boundary_velocity
void volume_external_residual_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase * > &bc)
boundary velocity contribution of volume external force.
Definition: heat_conduction_elem_base.cpp:466
elem_base.h
MAST::ElementBase::assembly
MAST::AssemblyBase & assembly()
Definition: elem_base.h:103
MAST::HeatConductionElementBase::internal_residual_sensitivity
virtual void internal_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f)
sensitivity of the internal force contribution to system residual
Definition: heat_conduction_elem_base.cpp:534
MAST::HeatConductionElementBase::_initialize_fem_gradient_operator
void _initialize_fem_gradient_operator(const unsigned int qp, const unsigned int dim, const MAST::FEBase &fe, std::vector< MAST::FEMOperatorMatrix > &dBmat)
For mass = true, the FEM operator matrix is initilized to the weak form of the Laplacian ...
Definition: heat_conduction_elem_base.cpp:1782
MAST::HeatConductionElementBase::surface_radiation_boundary_velocity
virtual void surface_radiation_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface flux on element volumetric domain...
Definition: heat_conduction_elem_base.cpp:1547
RealMatrixX
Matrix< Real, Dynamic, Dynamic > RealMatrixX
Definition: mast_data_types.h:44
MAST::HeatConductionElementBase::surface_convection_residual
virtual void surface_convection_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface convection.
Definition: heat_conduction_elem_base.cpp:1028
RealVectorX
Matrix< Real, Dynamic, 1 > RealVectorX
Definition: mast_data_types.h:35
MAST::HeatConductionElementBase::surface_convection_boundary_velocity
virtual void surface_convection_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface flux on element volumetric domain...
Definition: heat_conduction_elem_base.cpp:1252
MAST::HeatConductionElementBase::volume_heat_source_residual
virtual void volume_heat_source_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to volume heat source.
Definition: heat_conduction_elem_base.cpp:1622
MAST::GeomElem
This class acts as a wrapper around libMesh::Elem for the purpose of providing a uniform interface fo...
Definition: geom_elem.h:59
MAST::HeatConductionElementBase::elem_property
const MAST::ElementPropertyCardBase & elem_property()
returns a constant reference to the finite element object
Definition: heat_conduction_elem_base.h:73
MAST::HeatConductionElementBase::volume_external_residual
void volume_external_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase * > &bc)
volume external force contribution to system residual
Definition: heat_conduction_elem_base.cpp:292
MAST::HeatConductionElementBase::surface_flux_boundary_velocity
virtual void surface_flux_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface flux on element volumetric domain...
Definition: heat_conduction_elem_base.cpp:971
MAST::HeatConductionElementBase
This element implements the Galerkin discretization of the heat conduction problem with the flux pro...
Definition: heat_conduction_elem_base.h:53
MAST::HeatConductionElementBase::velocity_residual
virtual void velocity_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac_xdot, RealMatrixX &jac)
inertial force contribution to system residual
Definition: heat_conduction_elem_base.cpp:157
MAST::HeatConductionElementBase::surface_radiation_residual_sensitivity
virtual void surface_radiation_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector sensitivity and Jacobian due to surface radiation flux on element side...
Definition: heat_conduction_elem_base.cpp:1443
MAST::HeatConductionElementBase::velocity_residual_sensitivity
virtual void velocity_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f)
sensitivity of the damping force contribution to system residual
Definition: heat_conduction_elem_base.cpp:689
MAST::HeatConductionElementBase::internal_residual_boundary_velocity
virtual void internal_residual_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f)
sensitivity of the internal force contribution to system residual
Definition: heat_conduction_elem_base.cpp:602
MAST::HeatConductionElementBase::surface_radiation_residual
virtual void surface_radiation_residual(bool request_jacobian, RealVectorX &f, RealMatrixX &jac, const unsigned int s, MAST::BoundaryConditionBase &bc)
Calculates the residual vector and Jacobian due to surface radiation flux on side s...
Definition: heat_conduction_elem_base.cpp:1323
MAST::HeatConductionElementBase::_property
const MAST::ElementPropertyCardBase & _property
element property
Definition: heat_conduction_elem_base.h:388
MAST::HeatConductionElementBase::volume_external_residual_sensitivity
void volume_external_residual_sensitivity(const MAST::FunctionBase &p, RealVectorX &f, std::multimap< libMesh::subdomain_id_type, MAST::BoundaryConditionBase * > &bc)
sensitivity of the volume external force contribution to system residual
Definition: heat_conduction_elem_base.cpp:409
MAST::HeatConductionElementBase::velocity_residual_boundary_velocity
virtual void velocity_residual_boundary_velocity(const MAST::FunctionBase &p, RealVectorX &f, const unsigned int s, const MAST::FieldFunction< RealVectorX > &vel_f)
sensitivity of the internal force contribution to system residual
Definition: heat_conduction_elem_base.cpp:741
MAST::ElementBase
This is the base class for elements that implement calculation of finite element quantities over the ...
Definition: elem_base.h:72