FluidPressure Struct Reference

Fluid pressure forces. More...

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

Collaboration diagram for FluidPressure:

Classes
struct	FluidData
struct	MyTriangleFE
struct	OpCalculatePressure

Public Types
typedef int	MeshSetId

Public Member Functions
MyTriangleFE &	getLoopFe ()
	FluidPressure (MoFEM::Interface &m_field)
MoFEMErrorCode	addNeumannFluidPressureBCElements (const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")
MoFEMErrorCode	setNeumannFluidPressureFiniteElementOperators (string field_name, Vec F, bool allow_negative_pressure=true, bool ho_geometry=false)

Public Attributes
MoFEM::Interface &	mField
MyTriangleFE	fe
std::map< MeshSetId, FluidData >	setOfFluids
boost::ptr_vector< MethodForForceScaling >	methodsOp

Detailed Description

Fluid pressure forces.

Todo:

Implementation for large displacements

Examples

elasticity.cpp, and nonlinear_dynamics.cpp.

Definition at line 20 of file FluidPressure.hpp.

Member Typedef Documentation

MeshSetId

typedef int FluidPressure::MeshSetId

Definition at line 39 of file FluidPressure.hpp.

Constructor & Destructor Documentation

FluidPressure()

FluidPressure::FluidPressure ( MoFEM::Interface &  m_field )

inline

Definition at line 37 of file FluidPressure.hpp.

: mField(m_field), fe(mField) {}

Member Function Documentation

addNeumannFluidPressureBCElements()

MoFEMErrorCode FluidPressure::addNeumannFluidPressureBCElements	(	const std::string	field_name,
		const std::string	mesh_nodals_positions = "MESH_NODE_POSITIONS"
	)

Definition at line 67 of file FluidPressure.cpp.

                                                                       {
  MoFEMFunctionBegin;
 
  CHKERR mField.add_finite_element("FLUID_PRESSURE_FE", MF_ZERO);
  CHKERR mField.modify_finite_element_add_field_row("FLUID_PRESSURE_FE",
                                                    field_name);
  CHKERR mField.modify_finite_element_add_field_col("FLUID_PRESSURE_FE",
                                                    field_name);
  CHKERR mField.modify_finite_element_add_field_data("FLUID_PRESSURE_FE",
                                                     field_name);
  if (mField.check_field(mesh_nodals_positions)) {
    CHKERR mField.modify_finite_element_add_field_data("FLUID_PRESSURE_FE",
                                                       mesh_nodals_positions);
  }
 
  // takes skin of block of entities
  Skinner skin(&mField.get_moab());
  // loop over all blocksets and get data which name is FluidPressure
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(mField, BLOCKSET, bit)) {
 
    if (bit->getName().compare(0, 14, "FLUID_PRESSURE") == 0) {
 
      // get block attributes
      std::vector<double> attributes;
      CHKERR bit->getAttributes(attributes);
      if (attributes.size() < 7) {
        SETERRQ1(PETSC_COMM_SELF, 1,
                 "not enough block attributes to deffine fluid pressure "
                 "element, attributes.size() = %d ",
                 attributes.size());
      }
      setOfFluids[bit->getMeshsetId()].dEnsity = attributes[0];
      setOfFluids[bit->getMeshsetId()].aCCeleration.resize(3);
      setOfFluids[bit->getMeshsetId()].aCCeleration[0] = attributes[1];
      setOfFluids[bit->getMeshsetId()].aCCeleration[1] = attributes[2];
      setOfFluids[bit->getMeshsetId()].aCCeleration[2] = attributes[3];
      setOfFluids[bit->getMeshsetId()].zEroPressure.resize(3);
      setOfFluids[bit->getMeshsetId()].zEroPressure[0] = attributes[4];
      setOfFluids[bit->getMeshsetId()].zEroPressure[1] = attributes[5];
      setOfFluids[bit->getMeshsetId()].zEroPressure[2] = attributes[6];
      // get blok tetrahedron and triangles
      Range tets;
      CHKERR mField.get_moab().get_entities_by_type(bit->meshset, MBTET, tets,
                                                    true);
      Range tris;
      CHKERR mField.get_moab().get_entities_by_type(
          bit->meshset, MBTRI, setOfFluids[bit->getMeshsetId()].tRis, true);
      // this get triangles only on block surfaces
      Range tets_skin_tris;
      CHKERR skin.find_skin(0, tets, false, tets_skin_tris);
      setOfFluids[bit->getMeshsetId()].tRis.merge(tets_skin_tris);
      std::ostringstream ss;
      ss << setOfFluids[bit->getMeshsetId()] << std::endl;
      PetscPrintf(mField.get_comm(), ss.str().c_str());
 
      CHKERR mField.add_ents_to_finite_element_by_type(
          setOfFluids[bit->getMeshsetId()].tRis, MBTRI, "FLUID_PRESSURE_FE");
    }
  }
 
  MoFEMFunctionReturn(0);
}

getLoopFe()

MyTriangleFE& FluidPressure::getLoopFe ( )

inline

Definition at line 35 of file FluidPressure.hpp.

{ return fe; }

setNeumannFluidPressureFiniteElementOperators()

MoFEMErrorCode FluidPressure::setNeumannFluidPressureFiniteElementOperators	(	string	field_name,
		Vec	F,
		bool	allow_negative_pressure = true,
		bool	ho_geometry = false
	)

Definition at line 131 of file FluidPressure.cpp.

                                                                              {
  MoFEMFunctionBegin;
  std::map<MeshSetId, FluidData>::iterator sit = setOfFluids.begin();
  for (; sit != setOfFluids.end(); sit++) {
    // add finite element
    fe.getOpPtrVector().push_back(
        new OpCalculatePressure(field_name, F, sit->second, methodsOp,
                                allow_negative_pressure, ho_geometry));
  }
  MoFEMFunctionReturn(0);
}

Member Data Documentation

fe

MyTriangleFE FluidPressure::fe

Definition at line 34 of file FluidPressure.hpp.

methodsOp

boost::ptr_vector<MethodForForceScaling> FluidPressure::methodsOp

Definition at line 51 of file FluidPressure.hpp.

mField

MoFEM::Interface& FluidPressure::mField

Definition at line 22 of file FluidPressure.hpp.

setOfFluids

std::map<MeshSetId, FluidData> FluidPressure::setOfFluids

Definition at line 49 of file FluidPressure.hpp.

---

The documentation for this struct was generated from the following files:

• users_modules/basic_finite_elements/src/FluidPressure.hpp
• users_modules/basic_finite_elements/src/impl/FluidPressure.cpp