KratosMultiphysics
KRATOS Multiphysics (Kratos) is a framework for building parallel, multi-disciplinary simulation software, aiming at modularity, extensibility, and high performance. Kratos is written in C++, and counts with an extensive Python interface.
conservative_element_fc.h
Go to the documentation of this file.
1 // | / |
2 // ' / __| _` | __| _ \ __|
3 // . \ | ( | | ( |\__ `
4 // _|\_\_| \__,_|\__|\___/ ____/
5 // Multi-Physics
6 //
7 // License: BSD License
8 // Kratos default license: kratos/license.txt
9 //
10 // Main authors: Miguel Maso Sotomayor
11 //
12 
13 #ifndef KRATOS_CONSERVATIVE_ELEMENT_FC_H_INCLUDED
14 #define KRATOS_CONSERVATIVE_ELEMENT_FC_H_INCLUDED
15 
16 // System includes
17 
18 
19 // External includes
20 
21 
22 // Project includes
23 #include "conservative_element.h"
24 
25 namespace Kratos
26 {
29 
32 
36 
40 
44 
48 
50 template<std::size_t TNumNodes>
51 class ConservativeElementFC : public ConservativeElement<TNumNodes>
52 {
53 public:
56 
57  typedef std::size_t IndexType;
58 
59  typedef Node NodeType;
60 
62 
64 
66 
68 
70 
72 
74 
78 
80 
84 
89 
93  ConservativeElementFC(IndexType NewId, const NodesArrayType& ThisNodes) : ConservativeElementType(NewId, ThisNodes){}
94 
98  ConservativeElementFC(IndexType NewId, GeometryType::Pointer pGeometry) : ConservativeElementType(NewId, pGeometry){}
99 
103  ConservativeElementFC(IndexType NewId, GeometryType::Pointer pGeometry, typename PropertiesType::Pointer pProperties) : ConservativeElementType(NewId, pGeometry, pProperties){}
104 
108  ~ ConservativeElementFC() override {};
109 
113 
121  Element::Pointer Create(IndexType NewId, NodesArrayType const& ThisNodes, typename PropertiesType::Pointer pProperties) const override
122  {
123  return Kratos::make_intrusive<ConservativeElementFC<TNumNodes>>(NewId, this->GetGeometry().Create(ThisNodes), pProperties);
124  }
125 
133  Element::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, typename PropertiesType::Pointer pProperties) const override
134  {
135  return Kratos::make_intrusive<ConservativeElementFC<TNumNodes>>(NewId, pGeom, pProperties);
136  }
137 
144  Element::Pointer Clone(IndexType NewId, NodesArrayType const& ThisNodes) const override
145  {
146  Element::Pointer p_new_elem = Create(NewId, this->GetGeometry().Create(ThisNodes), this->pGetProperties());
147  p_new_elem->SetData(this->GetData());
148  p_new_elem->Set(Flags(*this));
149  return p_new_elem;
150  }
151 
155 
165  MatrixType& rDampingMatrix,
166  const ProcessInfo& rCurrentProcessInfo) override;
167 
171 
175  std::string Info() const override
176  {
177  return "ConservativeElementFC";
178  }
179 
181 
182 protected:
185 
187 
191 
193  BoundedMatrix<double,3,3>& rViscosity,
194  BoundedMatrix<double,2,2>& rDiffusion,
195  const ElementData& rData,
196  const array_1d<double,TNumNodes>& rN,
197  const BoundedMatrix<double,TNumNodes,2>& rDN_DX) override;
198 
200 
201 private:
204 
205  friend class Serializer;
206 
207  void save(Serializer& rSerializer) const override
208  {
210  }
211 
212  void load(Serializer& rSerializer) override
213  {
215  }
216 
220 
221 
225 
226 
228 
229 }; // Class ConservativeElementFC
230 
234 
235 
237 
239 
240 } // namespace Kratos.
241 
242 #endif // KRATOS_CONSERVATIVE_ELEMENT_FC_H_INCLUDED defined
Implementation of a linear element for shallow water problems.
Definition: conservative_element_fc.h:52
static constexpr IndexType mLocalSize
Definition: conservative_element_fc.h:186
ConservativeElementType::NodesArrayType NodesArrayType
Definition: conservative_element_fc.h:67
~ ConservativeElementFC() override
Destructor.
Definition: conservative_element_fc.h:108
ConservativeElementFC()
Default constructor.
Definition: conservative_element_fc.h:88
void CalculateDampingMatrix(MatrixType &rDampingMatrix, const ProcessInfo &rCurrentProcessInfo) override
Definition: conservative_element_fc.cpp:31
ConservativeElementType::MatrixType MatrixType
Definition: conservative_element_fc.h:73
std::size_t IndexType
Definition: conservative_element_fc.h:57
void CalculateArtificialViscosity(BoundedMatrix< double, 3, 3 > &rViscosity, BoundedMatrix< double, 2, 2 > &rDiffusion, const ElementData &rData, const array_1d< double, TNumNodes > &rN, const BoundedMatrix< double, TNumNodes, 2 > &rDN_DX) override
Definition: conservative_element_fc.cpp:75
ConservativeElementFC(IndexType NewId, GeometryType::Pointer pGeometry)
Constructor using Geometry.
Definition: conservative_element_fc.h:98
ConservativeElementFC(IndexType NewId, GeometryType::Pointer pGeometry, typename PropertiesType::Pointer pProperties)
Constructor using Geometry and Properties.
Definition: conservative_element_fc.h:103
ConservativeElementType::WaveElementType WaveElementType
Definition: conservative_element_fc.h:65
KRATOS_CLASS_INTRUSIVE_POINTER_DEFINITION(ConservativeElementFC)
ConservativeElement< TNumNodes > ConservativeElementType
Definition: conservative_element_fc.h:63
Element::Pointer Create(IndexType NewId, GeometryType::Pointer pGeom, typename PropertiesType::Pointer pProperties) const override
Create a new element pointer.
Definition: conservative_element_fc.h:133
Node NodeType
Definition: conservative_element_fc.h:59
ConservativeElementType::ElementData ElementData
Definition: conservative_element_fc.h:71
Element::Pointer Create(IndexType NewId, NodesArrayType const &ThisNodes, typename PropertiesType::Pointer pProperties) const override
Create a new element pointer.
Definition: conservative_element_fc.h:121
Element::Pointer Clone(IndexType NewId, NodesArrayType const &ThisNodes) const override
Create a new element pointer and clone the previous element data.
Definition: conservative_element_fc.h:144
ConservativeElementFC(IndexType NewId, const NodesArrayType &ThisNodes)
Constructor using an array of nodes.
Definition: conservative_element_fc.h:93
ConservativeElementType::PropertiesType PropertiesType
Definition: conservative_element_fc.h:69
Geometry< NodeType > GeometryType
Definition: conservative_element_fc.h:61
std::string Info() const override
Turn back information as a string.
Definition: conservative_element_fc.h:175
Implementation of a linear element for shallow water problems.
Definition: conservative_element.h:52
static constexpr IndexType mLocalSize
Definition: conservative_element.h:182
Base class for all Elements.
Definition: element.h:60
PropertiesType::Pointer pGetProperties()
returns the pointer to the property of the element. Does not throw an error, to allow copying of elem...
Definition: element.h:1014
std::size_t IndexType
Definition: flags.h:74
Flags()
Default constructor.
Definition: flags.h:119
GeometryType & GetGeometry()
Returns the reference of the geometry.
Definition: geometrical_object.h:158
DataValueContainer & GetData()
Definition: geometrical_object.h:212
Geometry base class.
Definition: geometry.h:71
virtual Pointer Create(PointsArrayType const &rThisPoints) const
Creates a new geometry pointer.
Definition: geometry.h:813
This class defines the node.
Definition: node.h:65
PointerVector is a container like stl vector but using a vector to store pointers to its data.
Definition: pointer_vector.h:72
ProcessInfo holds the current value of different solution parameters.
Definition: process_info.h:59
Properties encapsulates data shared by different Elements or Conditions. It can store any type of dat...
Definition: properties.h:69
The serialization consists in storing the state of an object into a storage format like data file or ...
Definition: serializer.h:123
Implementation of a linear element for shallow water problems.
Definition: wave_element.h:55
#define KRATOS_SERIALIZE_SAVE_BASE_CLASS(Serializer, BaseType)
Definition: define.h:812
#define KRATOS_SERIALIZE_LOAD_BASE_CLASS(Serializer, BaseType)
Definition: define.h:815
REF: G. R. Cowper, GAUSSIAN QUADRATURE FORMULAS FOR TRIANGLES.
Definition: mesh_condition.cpp:21
def load(f)
Definition: ode_solve.py:307
Definition: wave_element.h:299