MAST
Multidisciplinary-design Adaptation and Sensitivity Toolkit (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-2020 Manav Bhatia and MAST authors
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  const MAST::GeomElem& elem,
62  const MAST::ElementPropertyCardBase& p);
63 
64 
65  virtual ~HeatConductionElementBase();
66 
67 
71  const MAST::ElementPropertyCardBase&
72  elem_property() {
73  return _property;
74  }
75 
76 
80  virtual void
81  internal_residual (bool request_jacobian,
82  RealVectorX& f,
83  RealMatrixX& jac);
84 
85 
89  virtual void
90  velocity_residual (bool request_jacobian,
91  RealVectorX& f,
92  RealMatrixX& jac_xdot,
93  RealMatrixX& jac);
94 
98  void
99  side_external_residual (bool request_jacobian,
100  RealVectorX& f,
101  RealMatrixX& jac,
102  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
103 
107  void
108  volume_external_residual (bool request_jacobian,
109  RealVectorX& f,
110  RealMatrixX& jac,
111  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
112 
116  virtual void
117  internal_residual_sensitivity (const MAST::FunctionBase& p,
118  RealVectorX& f);
119 
123  virtual void
124  internal_residual_boundary_velocity (const MAST::FunctionBase& p,
125  RealVectorX& f,
126  const unsigned int s,
127  const MAST::FieldFunction<RealVectorX>& vel_f);
128 
132  virtual void
133  velocity_residual_sensitivity (const MAST::FunctionBase& p,
134  RealVectorX& f);
135 
139  virtual void
140  velocity_residual_boundary_velocity (const MAST::FunctionBase& p,
141  RealVectorX& f,
142  const unsigned int s,
143  const MAST::FieldFunction<RealVectorX>& vel_f);
144 
148  void
149  side_external_residual_sensitivity (const MAST::FunctionBase& p,
150  RealVectorX& f,
151  std::multimap<libMesh::boundary_id_type, MAST::BoundaryConditionBase*>& bc);
152 
156  void
157  volume_external_residual_sensitivity (const MAST::FunctionBase& p,
158  RealVectorX& f,
159  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
160 
164  void
165  volume_external_residual_boundary_velocity(const MAST::FunctionBase& p,
166  RealVectorX& f,
167  const unsigned int s,
168  const MAST::FieldFunction<RealVectorX>& vel_f,
169  std::multimap<libMesh::subdomain_id_type, MAST::BoundaryConditionBase*>& bc);
170 
171  protected:
172 
173 
174 
179  virtual void surface_flux_residual(bool request_jacobian,
180  RealVectorX& f,
181  RealMatrixX& jac,
182  const unsigned int s,
183  MAST::BoundaryConditionBase& bc);
184 
190  virtual void surface_flux_residual(bool request_jacobian,
191  RealVectorX& f,
192  RealMatrixX& jac,
193  MAST::BoundaryConditionBase& bc);
194 
199  virtual void surface_flux_residual_sensitivity(const MAST::FunctionBase& p,
200  RealVectorX& f,
201  const unsigned int s,
202  MAST::BoundaryConditionBase& bc);
203 
209  virtual void surface_flux_residual_sensitivity(const MAST::FunctionBase& p,
210  RealVectorX& f,
211  MAST::BoundaryConditionBase& bc);
212 
218  virtual void surface_flux_boundary_velocity(const MAST::FunctionBase& p,
219  RealVectorX& f,
220  const unsigned int s,
221  const MAST::FieldFunction<RealVectorX>& vel_f,
222  MAST::BoundaryConditionBase& bc);
223 
224 
229  virtual void surface_convection_residual(bool request_jacobian,
230  RealVectorX& f,
231  RealMatrixX& jac,
232  const unsigned int s,
233  MAST::BoundaryConditionBase& bc);
234 
235 
241  virtual void surface_convection_residual(bool request_jacobian,
242  RealVectorX& f,
243  RealMatrixX& jac,
244  MAST::BoundaryConditionBase& bc);
245 
246 
251  virtual void
252  surface_convection_residual_sensitivity(const MAST::FunctionBase& p,
253  RealVectorX& f,
254  const unsigned int s,
255  MAST::BoundaryConditionBase& bc);
256 
262  virtual void
263  surface_convection_residual_sensitivity(const MAST::FunctionBase& p,
264  RealVectorX& f,
265  MAST::BoundaryConditionBase& bc);
266 
272  virtual void
273  surface_convection_boundary_velocity(const MAST::FunctionBase& p,
274  RealVectorX& f,
275  const unsigned int s,
276  const MAST::FieldFunction<RealVectorX>& vel_f,
277  MAST::BoundaryConditionBase& bc);
278 
283  virtual void surface_radiation_residual(bool request_jacobian,
284  RealVectorX& f,
285  RealMatrixX& jac,
286  const unsigned int s,
287  MAST::BoundaryConditionBase& bc);
288 
289 
294  virtual void surface_radiation_residual(bool request_jacobian,
295  RealVectorX& f,
296  RealMatrixX& jac,
297  MAST::BoundaryConditionBase& bc);
298 
299 
304  virtual void
305  surface_radiation_residual_sensitivity(const MAST::FunctionBase& p,
306  RealVectorX& f,
307  const unsigned int s,
308  MAST::BoundaryConditionBase& bc);
309 
310 
315  virtual void
316  surface_radiation_residual_sensitivity(const MAST::FunctionBase& p,
317  RealVectorX& f,
318  MAST::BoundaryConditionBase& bc);
319 
325  virtual void
326  surface_radiation_boundary_velocity(const MAST::FunctionBase& p,
327  RealVectorX& f,
328  const unsigned int s,
329  const MAST::FieldFunction<RealVectorX>& vel_f,
330  MAST::BoundaryConditionBase& bc);
331 
336  virtual void volume_heat_source_residual(bool request_jacobian,
337  RealVectorX& f,
338  RealMatrixX& jac,
339  MAST::BoundaryConditionBase& bc);
340 
341 
346  virtual void
347  volume_heat_source_residual_sensitivity(const MAST::FunctionBase& p,
348  RealVectorX& f,
349  MAST::BoundaryConditionBase& bc);
350 
355  virtual void
356  volume_heat_source_boundary_velocity(const MAST::FunctionBase& p,
357  RealVectorX& f,
358  const unsigned int s,
359  const MAST::FieldFunction<RealVectorX>& vel_f,
360  MAST::BoundaryConditionBase& bc);
361 
367  void _initialize_mass_fem_operator(const unsigned int qp,
368  const MAST::FEBase& fe,
369  MAST::FEMOperatorMatrix& Bmat);
370 
371 
379  void _initialize_fem_gradient_operator(const unsigned int qp,
380  const unsigned int dim,
381  const MAST::FEBase& fe,
382  std::vector<MAST::FEMOperatorMatrix>& dBmat);
383 
387  const MAST::ElementPropertyCardBase& _property;
388  };
389 
390 }
391 
392 
393 #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:1218
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:1857
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:1751
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:54
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:958
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:877
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:1796
elem_base.h
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:554
MAST::HeatConductionElementBase::HeatConductionElementBase
HeatConductionElementBase(MAST::SystemInitialization &sys, const MAST::GeomElem &elem, const MAST::ElementPropertyCardBase &p)
Constructor.
Definition: heat_conduction_elem_base.cpp:36
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:1880
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:250
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:1636
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:1101
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:1333
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:1711
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:369
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:72
MAST::ElementBase::elem
const MAST::GeomElem & elem() const
Definition: elem_base.h:109
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:312
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:429
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:486
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:1044
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:156
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:1525
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:709
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:622
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:1404
MAST::HeatConductionElementBase::_property
const MAST::ElementPropertyCardBase & _property
element property
Definition: heat_conduction_elem_base.h:387
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:783
MAST::ElementBase
This is the base class for elements that implement calculation of finite element quantities over the ...
Definition: elem_base.h:72