v0.14.0
Classes | Public Types | Static Public Member Functions | List of all members
NavierStokesElement Struct Reference

Element for simulating viscous fluid flow. More...

#include <users_modules/basic_finite_elements/src/NavierStokesElement.hpp>

Classes

struct  BlockData
 
struct  CommonData
 
struct  FaceRule
 
struct  LoadScale
 
struct  OpAssembleLhs
 Base class for operators assembling LHS. More...
 
struct  OpAssembleLhsDiagLin
 Assemble linear (symmetric) part of the diagonal block of the LHS Operator for assembling linear (symmetric) part of the diagonal block of the LHS. More...
 
struct  OpAssembleLhsDiagNonLin
 Assemble non-linear (non-symmetric) part of the diagonal block of the LHS Operator for assembling non-linear (non-symmetric) part of the diagonal block of the LHS. More...
 
struct  OpAssembleLhsOffDiag
 Assemble off-diagonal block of the LHS Operator for assembling off-diagonal block of the LHS. More...
 
struct  OpAssembleRhs
 Base class for operators assembling RHS. More...
 
struct  OpAssembleRhsPressure
 Assemble the pressure component of the RHS vector. More...
 
struct  OpAssembleRhsVelocityLin
 Assemble linear part of the velocity component of the RHS vector. More...
 
struct  OpAssembleRhsVelocityNonLin
 Assemble non-linear part of the velocity component of the RHS vector. More...
 
struct  OpCalcDragForce
 Calculate drag force on the fluid-solid interface. More...
 
struct  OpCalcDragTraction
 Calculate drag traction on the fluid-solid interface. More...
 
struct  OpCalcVolumeFlux
 calculating volumetric flux More...
 
struct  OpPostProcDrag
 Post processing output of drag traction on the fluid-solid interface. More...
 
struct  OpPostProcVorticity
 Post processing output of the vorticity criterion levels. More...
 
struct  VolRule
 Set integration rule to volume elements. More...
 

Public Types

using PostProcVol = PostProcBrokenMeshInMoab< VolumeElementForcesAndSourcesCore >
 
using PostProcFace = PostProcBrokenMeshInMoab< FaceElementForcesAndSourcesCore >
 
using UserDataOperator = MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator
 
using FaceUserDataOperator = MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator
 
using EntData = EntitiesFieldData::EntData
 

Static Public Member Functions

static MoFEMErrorCode addElement (MoFEM::Interface &m_field, const string element_name, const string velocity_field_name, const string pressure_field_name, const string mesh_field_name, const int dim=3, Range *ents=nullptr)
 Setting up elements. More...
 
static MoFEMErrorCode setNavierStokesOperators (boost::shared_ptr< VolumeElementForcesAndSourcesCore > feRhs, boost::shared_ptr< VolumeElementForcesAndSourcesCore > feLhs, const std::string velocity_field, const std::string pressure_field, boost::shared_ptr< CommonData > common_data, const EntityType type=MBTET)
 Setting up operators for solving Navier-Stokes equations. More...
 
static MoFEMErrorCode setStokesOperators (boost::shared_ptr< VolumeElementForcesAndSourcesCore > feRhs, boost::shared_ptr< VolumeElementForcesAndSourcesCore > feLhs, const std::string velocity_field, const std::string pressure_field, boost::shared_ptr< CommonData > common_data, const EntityType type=MBTET)
 Setting up operators for solving Stokes equations. More...
 
static MoFEMErrorCode setCalcDragOperators (boost::shared_ptr< FaceElementForcesAndSourcesCore > dragFe, boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide > sideDragFe, std::string side_fe_name, const std::string velocity_field, const std::string pressure_field, boost::shared_ptr< CommonData > common_data)
 Setting up operators for calculating drag force on the solid surface. More...
 
static MoFEMErrorCode setPostProcDragOperators (boost::shared_ptr< PostProcFace > postProcDragPtr, boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide > sideDragFe, std::string side_fe_name, const std::string velocity_field, const std::string pressure_field, boost::shared_ptr< CommonData > common_data)
 Setting up operators for post processing output of drag traction. More...
 
static MoFEMErrorCode setCalcVolumeFluxOperators (boost::shared_ptr< VolumeElementForcesAndSourcesCore > fe_flux_ptr, const std::string velocity_field, boost::shared_ptr< CommonData > common_data, const EntityType type=MBTET)
 Setting up operators for calculation of volume flux. More...
 

Detailed Description

Element for simulating viscous fluid flow.

Examples
NavierStokesElement.hpp.

Definition at line 25 of file NavierStokesElement.hpp.

Member Typedef Documentation

◆ EntData

using NavierStokesElement::EntData = EntitiesFieldData::EntData

Definition at line 36 of file NavierStokesElement.hpp.

◆ FaceUserDataOperator

using NavierStokesElement::FaceUserDataOperator = MoFEM::FaceElementForcesAndSourcesCore::UserDataOperator

Definition at line 35 of file NavierStokesElement.hpp.

◆ PostProcFace

using NavierStokesElement::PostProcFace = PostProcBrokenMeshInMoab<FaceElementForcesAndSourcesCore>

Definition at line 30 of file NavierStokesElement.hpp.

◆ PostProcVol

using NavierStokesElement::PostProcVol = PostProcBrokenMeshInMoab<VolumeElementForcesAndSourcesCore>

Definition at line 28 of file NavierStokesElement.hpp.

◆ UserDataOperator

using NavierStokesElement::UserDataOperator = MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator

Definition at line 33 of file NavierStokesElement.hpp.

Member Function Documentation

◆ addElement()

static MoFEMErrorCode NavierStokesElement::addElement ( MoFEM::Interface &  m_field,
const string  element_name,
const string  velocity_field_name,
const string  pressure_field_name,
const string  mesh_field_name,
const int  dim = 3,
Range *  ents = nullptr 
)
inlinestatic

Setting up elements.

This functions adds element to the database, adds provided fields to rows and columns of the element, provides access of the element to the fields data and adds entities of particular dimension (or a given range of entities to the element)

Parameters
m_fieldMoFEM interface
element_nameName of the element
velocity_field_nameName of the velocity field
pressure_field_nameName of the pressure field
mesh_field_nameName for mesh node positions field
entsRange of entities to be added to element
Returns
Error code
Examples
navier_stokes.cpp.

Definition at line 129 of file NavierStokesElement.hpp.

135  {
136  MoFEMFunctionBegin;
137 
138  CHKERR m_field.add_finite_element(element_name);
139 
140  CHKERR m_field.modify_finite_element_add_field_row(element_name,
141  velocity_field_name);
142  CHKERR m_field.modify_finite_element_add_field_col(element_name,
143  velocity_field_name);
144  CHKERR m_field.modify_finite_element_add_field_data(element_name,
145  velocity_field_name);
146 
147  CHKERR m_field.modify_finite_element_add_field_row(element_name,
148  pressure_field_name);
149  CHKERR m_field.modify_finite_element_add_field_col(element_name,
150  pressure_field_name);
151  CHKERR m_field.modify_finite_element_add_field_data(element_name,
152  pressure_field_name);
153 
154  CHKERR m_field.modify_finite_element_add_field_data(element_name,
155  mesh_field_name);
156 
157  if (ents != nullptr) {
158  CHKERR m_field.add_ents_to_finite_element_by_dim(*ents, dim,
159  element_name);
160  } else {
161  CHKERR m_field.add_ents_to_finite_element_by_dim(0, dim, element_name);
162  }
163 
164  MoFEMFunctionReturn(0);
165  }

◆ setCalcDragOperators()

MoFEMErrorCode NavierStokesElement::setCalcDragOperators ( boost::shared_ptr< FaceElementForcesAndSourcesCore >  dragFe,
boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide >  sideDragFe,
std::string  side_fe_name,
const std::string  velocity_field,
const std::string  pressure_field,
boost::shared_ptr< CommonData >  common_data 
)
static

Setting up operators for calculating drag force on the solid surface.

Pushes operators for caluclating drag force components on the fluid-solid interface

Parameters
dragFepointer to face element instance
sideDragFepointer to volume on side element instance
velocity_fieldname of the velocity field
pressure_fieldname of the pressure field
common_datapointer to common data object
Returns
Error code
Examples
navier_stokes.cpp, and NavierStokesElement.cpp.

Definition at line 150 of file NavierStokesElement.cpp.

155  {
156  MoFEMFunctionBegin;
157 
158  auto det_ptr = boost::make_shared<VectorDouble>();
159  auto jac_ptr = boost::make_shared<MatrixDouble>();
160  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
161 
162  for (auto &sit : common_data->setOfFacesData) {
163  sideDragFe->getOpPtrVector().push_back(
164  new OpCalculateVectorFieldGradient<3, 3>(velocity_field,
165  common_data->gradVelPtr));
166  dragFe->getOpPtrVector().push_back(new OpCalculateHOJac<2>(jac_ptr));
167  dragFe->getOpPtrVector().push_back(
168  new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
169  dragFe->getOpPtrVector().push_back(
170  new OpSetHOInvJacToScalarBases<2>(H1, inv_jac_ptr));
171 
172  dragFe->getOpPtrVector().push_back(new OpCalculateScalarFieldValues(
173  pressure_field, common_data->pressPtr));
174 
175  dragFe->getOpPtrVector().push_back(
176  new NavierStokesElement::OpCalcDragTraction(
177  velocity_field, sideDragFe, side_fe_name, common_data, sit.second));
178  dragFe->getOpPtrVector().push_back(new NavierStokesElement::OpCalcDragForce(
179  velocity_field, common_data, sit.second));
180  }
181  MoFEMFunctionReturn(0);
182 };

◆ setCalcVolumeFluxOperators()

MoFEMErrorCode NavierStokesElement::setCalcVolumeFluxOperators ( boost::shared_ptr< VolumeElementForcesAndSourcesCore >  fe_flux_ptr,
const std::string  velocity_field,
boost::shared_ptr< CommonData >  common_data,
const EntityType  type = MBTET 
)
static

Setting up operators for calculation of volume flux.

Examples
NavierStokesElement.cpp.

Definition at line 124 of file NavierStokesElement.cpp.

127  {
128  MoFEMFunctionBegin;
129 
130  for (auto &sit : common_data->setOfBlocksData) {
131 
132  if (type == MBPRISM) {
133  boost::shared_ptr<MatrixDouble> inv_jac_ptr(new MatrixDouble);
134  fe_flux_ptr->getOpPtrVector().push_back(
135  new OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms());
136  fe_flux_ptr->getOpPtrVector().push_back(
137  new OpCalculateInvJacForFatPrism(inv_jac_ptr));
138  fe_flux_ptr->getOpPtrVector().push_back(
139  new OpSetInvJacH1ForFatPrism(inv_jac_ptr));
140  }
141  fe_flux_ptr->getOpPtrVector().push_back(new OpCalculateVectorFieldValues<3>(
142  velocity_field, common_data->velPtr));
143  fe_flux_ptr->getOpPtrVector().push_back(
144  new OpCalcVolumeFlux(velocity_field, common_data, sit.second));
145  }
146 
147  MoFEMFunctionReturn(0);
148 };

◆ setNavierStokesOperators()

MoFEMErrorCode NavierStokesElement::setNavierStokesOperators ( boost::shared_ptr< VolumeElementForcesAndSourcesCore >  feRhs,
boost::shared_ptr< VolumeElementForcesAndSourcesCore >  feLhs,
const std::string  velocity_field,
const std::string  pressure_field,
boost::shared_ptr< CommonData >  common_data,
const EntityType  type = MBTET 
)
static

Setting up operators for solving Navier-Stokes equations.

Pushes operators for solving Navier-Stokes equations to pipelines of RHS and LHS element instances

Parameters
feRhspointer to RHS element instance
feLhspointer to LHS element instance
velocity_fieldname of the velocity field
pressure_fieldname of the pressure field
common_datapointer to common data object
typetype of entities in the domain
Returns
Error code
Examples
navier_stokes.cpp, and NavierStokesElement.cpp.

Definition at line 22 of file NavierStokesElement.cpp.

26  {
27  MoFEMFunctionBegin;
28 
29  for (auto &sit : common_data->setOfBlocksData) {
30 
31  if (type == MBPRISM) {
32  boost::shared_ptr<MatrixDouble> inv_jac_ptr(new MatrixDouble);
33  fe_lhs_ptr->getOpPtrVector().push_back(
34  new OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms());
35  fe_lhs_ptr->getOpPtrVector().push_back(
36  new OpCalculateInvJacForFatPrism(inv_jac_ptr));
37  fe_lhs_ptr->getOpPtrVector().push_back(
38  new OpSetInvJacH1ForFatPrism(inv_jac_ptr));
39  fe_rhs_ptr->getOpPtrVector().push_back(
40  new OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms());
41  fe_rhs_ptr->getOpPtrVector().push_back(
42  new OpCalculateInvJacForFatPrism(inv_jac_ptr));
43  fe_rhs_ptr->getOpPtrVector().push_back(
44  new OpSetInvJacH1ForFatPrism(inv_jac_ptr));
45  }
46 
47  fe_lhs_ptr->getOpPtrVector().push_back(new OpCalculateVectorFieldValues<3>(
48  velocity_field, common_data->velPtr));
49  fe_lhs_ptr->getOpPtrVector().push_back(
50  new OpCalculateVectorFieldGradient<3, 3>(velocity_field,
51  common_data->gradVelPtr));
52 
53  fe_lhs_ptr->getOpPtrVector().push_back(new OpAssembleLhsDiagLin(
54  velocity_field, velocity_field, common_data, sit.second));
55  fe_lhs_ptr->getOpPtrVector().push_back(new OpAssembleLhsDiagNonLin(
56  velocity_field, velocity_field, common_data, sit.second));
57  fe_lhs_ptr->getOpPtrVector().push_back(new OpAssembleLhsOffDiag(
58  velocity_field, pressure_field, common_data, sit.second));
59 
60  fe_rhs_ptr->getOpPtrVector().push_back(new OpCalculateVectorFieldValues<3>(
61  velocity_field, common_data->velPtr));
62  fe_rhs_ptr->getOpPtrVector().push_back(
63  new OpCalculateVectorFieldGradient<3, 3>(velocity_field,
64  common_data->gradVelPtr));
65  fe_rhs_ptr->getOpPtrVector().push_back(new OpCalculateScalarFieldValues(
66  pressure_field, common_data->pressPtr));
67 
68  fe_rhs_ptr->getOpPtrVector().push_back(
69  new OpAssembleRhsVelocityLin(velocity_field, common_data, sit.second));
70  fe_rhs_ptr->getOpPtrVector().push_back(new OpAssembleRhsVelocityNonLin(
71  velocity_field, common_data, sit.second));
72  fe_rhs_ptr->getOpPtrVector().push_back(
73  new OpAssembleRhsPressure(pressure_field, common_data, sit.second));
74  }
75 
76  MoFEMFunctionReturn(0);
77 };

◆ setPostProcDragOperators()

MoFEMErrorCode NavierStokesElement::setPostProcDragOperators ( boost::shared_ptr< PostProcFace >  postProcDragPtr,
boost::shared_ptr< VolumeElementForcesAndSourcesCoreOnSide >  sideDragFe,
std::string  side_fe_name,
const std::string  velocity_field,
const std::string  pressure_field,
boost::shared_ptr< CommonData >  common_data 
)
static

Setting up operators for post processing output of drag traction.

Pushes operators for post processing ouput of drag traction components on the fluid-solid interface

Parameters
dragFepointer to face element instance
sideDragFepointer to volume on side element instance
velocity_fieldname of the velocity field
pressure_fieldname of the pressure field
common_datapointer to common data object
Returns
Error code
Examples
navier_stokes.cpp, and NavierStokesElement.cpp.

Definition at line 184 of file NavierStokesElement.cpp.

189  {
190  MoFEMFunctionBegin;
191 
192  auto det_ptr = boost::make_shared<VectorDouble>();
193  auto jac_ptr = boost::make_shared<MatrixDouble>();
194  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
195 
196  for (auto &sit : common_data->setOfFacesData) {
197  sideDragFe->getOpPtrVector().push_back(
198  new OpCalculateVectorFieldGradient<3, 3>(velocity_field,
199  common_data->gradVelPtr));
200  postProcDragPtr->getOpPtrVector().push_back(
201  new OpCalculateHOJac<2>(jac_ptr));
202  postProcDragPtr->getOpPtrVector().push_back(
203  new OpInvertMatrix<2>(jac_ptr, det_ptr, inv_jac_ptr));
204  postProcDragPtr->getOpPtrVector().push_back(
205  new OpSetHOInvJacToScalarBases<2>(H1, inv_jac_ptr));
206 
207  postProcDragPtr->getOpPtrVector().push_back(
208  new OpCalculateVectorFieldValues<3>(velocity_field,
209  common_data->velPtr));
210  postProcDragPtr->getOpPtrVector().push_back(
211  new OpCalculateScalarFieldValues(pressure_field,
212  common_data->pressPtr));
213 
214  postProcDragPtr->getOpPtrVector().push_back(
215  new NavierStokesElement::OpCalcDragTraction(
216  velocity_field, sideDragFe, side_fe_name, common_data, sit.second));
217  postProcDragPtr->getOpPtrVector().push_back(
218  new NavierStokesElement::OpPostProcDrag(
219  velocity_field, postProcDragPtr->getPostProcMesh(),
220  postProcDragPtr->getMapGaussPts(), common_data, sit.second));
221  }
222  MoFEMFunctionReturn(0);
223 };

◆ setStokesOperators()

MoFEMErrorCode NavierStokesElement::setStokesOperators ( boost::shared_ptr< VolumeElementForcesAndSourcesCore >  feRhs,
boost::shared_ptr< VolumeElementForcesAndSourcesCore >  feLhs,
const std::string  velocity_field,
const std::string  pressure_field,
boost::shared_ptr< CommonData >  common_data,
const EntityType  type = MBTET 
)
static

Setting up operators for solving Stokes equations.

Pushes operators for solving Stokes equations to pipelines of RHS and LHS element instances

Parameters
feRhspointer to RHS element instance
feLhspointer to LHS element instance
velocity_fieldname of the velocity field
pressure_fieldname of the pressure field
common_datapointer to common data object
typetype of entities in the domain
Returns
Error code
Examples
navier_stokes.cpp, and NavierStokesElement.cpp.

Definition at line 79 of file NavierStokesElement.cpp.

83  {
84  MoFEMFunctionBegin;
85 
86  for (auto &sit : common_data->setOfBlocksData) {
87 
88  if (type == MBPRISM) {
89  boost::shared_ptr<MatrixDouble> inv_jac_ptr(new MatrixDouble);
90  fe_lhs_ptr->getOpPtrVector().push_back(
91  new OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms());
92  fe_lhs_ptr->getOpPtrVector().push_back(
93  new OpCalculateInvJacForFatPrism(inv_jac_ptr));
94  fe_lhs_ptr->getOpPtrVector().push_back(
95  new OpSetInvJacH1ForFatPrism(inv_jac_ptr));
96  fe_rhs_ptr->getOpPtrVector().push_back(
97  new OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms());
98  fe_rhs_ptr->getOpPtrVector().push_back(
99  new OpCalculateInvJacForFatPrism(inv_jac_ptr));
100  fe_rhs_ptr->getOpPtrVector().push_back(
101  new OpSetInvJacH1ForFatPrism(inv_jac_ptr));
102  }
103 
104  fe_lhs_ptr->getOpPtrVector().push_back(new OpAssembleLhsDiagLin(
105  velocity_field, velocity_field, common_data, sit.second));
106  fe_lhs_ptr->getOpPtrVector().push_back(new OpAssembleLhsOffDiag(
107  velocity_field, pressure_field, common_data, sit.second));
108 
109  fe_rhs_ptr->getOpPtrVector().push_back(
110  new OpCalculateVectorFieldGradient<3, 3>(velocity_field,
111  common_data->gradVelPtr));
112 
113  fe_rhs_ptr->getOpPtrVector().push_back(new OpCalculateScalarFieldValues(
114  pressure_field, common_data->pressPtr));
115  fe_rhs_ptr->getOpPtrVector().push_back(
116  new OpAssembleRhsVelocityLin(velocity_field, common_data, sit.second));
117  fe_rhs_ptr->getOpPtrVector().push_back(
118  new OpAssembleRhsPressure(pressure_field, common_data, sit.second));
119  }
120 
121  MoFEMFunctionReturn(0);
122 };
The documentation for this struct was generated from the following files:
NavierStokesElement::OpPostProcDrag
Post processing output of drag traction on the fluid-solid interface.
Definition: NavierStokesElement.hpp:552
NavierStokesElement::OpCalcDragTraction
Calculate drag traction on the fluid-solid interface.
Definition: NavierStokesElement.hpp:521
H1
@ H1
continuous field
Definition: definitions.h:85
MoFEM::CoreInterface::modify_finite_element_add_field_row
virtual MoFEMErrorCode modify_finite_element_add_field_row(const std::string &fe_name, const std::string name_row)=0
set field row which finite element use
MoFEM::OpCalculateVectorFieldValues
Get values at integration pts for tensor filed rank 1, i.e. vector field.
Definition: UserDataOperators.hpp:466
MoFEM::Types::MatrixDouble
UBlasMatrix< double > MatrixDouble
Definition: Types.hpp:77
MoFEM::OpSetHOInvJacToScalarBases< 2 >
Definition: HODataOperators.hpp:78
MoFEM::OpCalculateHOJac< 2 >
Definition: HODataOperators.hpp:273
MoFEM::OpCalculateInvJacForFatPrism
Calculate inverse of jacobian for face element.
Definition: UserDataOperators.hpp:3145
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:548
MoFEM::CoreInterface::add_finite_element
virtual MoFEMErrorCode add_finite_element(const std::string &fe_name, enum MoFEMTypes bh=MF_EXCL, int verb=DEFAULT_VERBOSITY)=0
add finite element
MoFEM::CoreInterface::modify_finite_element_add_field_col
virtual MoFEMErrorCode modify_finite_element_add_field_col(const std::string &fe_name, const std::string name_row)=0
set field col which finite element use
MoFEM::OpSetInvJacH1ForFatPrism
Transform local reference derivatives of shape functions to global derivatives.
Definition: UserDataOperators.hpp:3172
convert.type
type
Definition: convert.py:64
MoFEM::OpCalculateScalarFieldValues
Get value at integration points for scalar field.
Definition: UserDataOperators.hpp:82
MoFEM::CoreInterface::add_ents_to_finite_element_by_dim
virtual MoFEMErrorCode add_ents_to_finite_element_by_dim(const EntityHandle entities, const int dim, const std::string &name, const bool recursive=true)=0
add entities to finite element
MoFEM::CoreInterface::modify_finite_element_add_field_data
virtual MoFEMErrorCode modify_finite_element_add_field_data(const std::string &fe_name, const std::string name_field)=0
set finite element field data
MoFEM::OpMultiplyDeterminantOfJacobianAndWeightsForFatPrisms
Operator for fat prism element updating integration weights in the volume.
Definition: UserDataOperators.hpp:3124
MoFEM::OpCalculateVectorFieldGradient
Get field gradients at integration pts for scalar filed rank 0, i.e. vector field.
Definition: UserDataOperators.hpp:1535
MoFEM::OpInvertMatrix
Definition: UserDataOperators.hpp:3249
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
NavierStokesElement::OpCalcDragForce
Calculate drag force on the fluid-solid interface.
Definition: NavierStokesElement.hpp:490
Generated by Doxygen 1.8.17 and hosted at University of Glasgow