Set of high-level function declaring elements and setting operators to apply forces/fluxes. More...

#include "users_modules/basic_finite_elements/src/SurfacePressure.hpp"

Static Public Member Functions
static MoFEMErrorCode	addNeumannBCElements (MoFEM::Interface &m_field, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS", Range *intersect_ptr=NULL)
	Declare finite element.

static MoFEMErrorCode	setMomentumFluxOperators (MoFEM::Interface &m_field, boost::ptr_map< std::string, NeumannForcesSurface > &neumann_forces, Vec F, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")
	Set operators to finite elements calculating right hand side vector.

static MoFEMErrorCode	addNeumannFluxBCElements (MoFEM::Interface &m_field, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")

static MoFEMErrorCode	setMassFluxOperators (MoFEM::Interface &m_field, boost::ptr_map< std::string, NeumannForcesSurface > &neumann_forces, Vec F, const std::string field_name, const std::string mesh_nodals_positions="MESH_NODE_POSITIONS")

Detailed Description

Set of high-level function declaring elements and setting operators to apply forces/fluxes.

Definition at line 1261 of file SurfacePressure.hpp.

Member Function Documentation

◆ addNeumannBCElements()

MoFEMErrorCode MetaNeumannForces::addNeumannBCElements	(	MoFEM::Interface &	m_field,
		const std::string	field_name,
		const std::string	mesh_nodals_positions = `"MESH_NODE_POSITIONS"`,
		Range *	intersect_ptr = `NULL`
	)

static

Declare finite element.

Search cubit sidesets and blocksets with pressure bc and declare surface elemen

Block set has to have name “PRESSURE”. Can have name “PRESSURE_01” or any other name with prefix. The first attribute of block set is pressure value.

Parameters

m_field	Interface insurance
field_name	Field name (f.e. DISPLACEMENT)
mesh_nodals_positions	Name of field on which ho-geometry is defined
intersect_ptr	Pointer to range to interect meshset entities

Returns: Error code

Examples: mofem/tutorials/cor-10_navier_stokes/navier_stokes.cpp, mofem/tutorials/cor-7_elasticity_mixed_formulation/elasticity_mixed_formulation.cpp, mofem/users_modules/basic_finite_elements/elasticity/elasticity.cpp, and mofem/users_modules/basic_finite_elements/nonlinear_elasticity/nonlinear_dynamics.cpp.

Definition at line 1974 of file SurfacePressure.cpp.

                                                                 {
  MoFEMFunctionBegin;
 
  // Define boundary element that operates on rows, columns and data of a
  // given field
  CHKERR m_field.add_finite_element("FORCE_FE", MF_ZERO);
  CHKERR m_field.modify_finite_element_add_field_row("FORCE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_col("FORCE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_data("FORCE_FE", field_name);
  if (m_field.check_field(mesh_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("FORCE_FE",
                                                        mesh_nodals_positions);
  }
  // Add entities to that element, here we add all triangles with FORCESET
  // from cubit
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                  it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    if (intersect_ptr)
      tris = intersect(tris, *intersect_ptr);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FORCE_FE");
  }
 
  CHKERR m_field.add_finite_element("PRESSURE_FE", MF_ZERO);
  CHKERR m_field.modify_finite_element_add_field_row("PRESSURE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_col("PRESSURE_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_data("PRESSURE_FE",
                                                      field_name);
  if (m_field.check_field(mesh_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("PRESSURE_FE",
                                                        mesh_nodals_positions);
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    if (intersect_ptr)
      tris = intersect(tris, *intersect_ptr);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
  }
 
  // Reading forces from BLOCKSET
 
  const string block_set_force_name("FORCE");
  // search for block named FORCE and add its attributes to FORCE_FE element
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_force_name.length(),
                              block_set_force_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FORCE_FE");
    }
  }
  // search for block named PRESSURE and add its attributes to PRESSURE_FE
  // element
  const string block_set_pressure_name("PRESSURE");
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_pressure_name.length(),
                              block_set_pressure_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
    }
  }
 
  // search for block named LINEAR_PRESSURE and add its attributes to
  // PRESSURE_FE element
  const string block_set_linear_pressure_name("LINEAR_PRESSURE");
  for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
    if (it->getName().compare(0, block_set_linear_pressure_name.length(),
                              block_set_linear_pressure_name) == 0) {
      Range tris;
      CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                          true);
      if (intersect_ptr)
        tris = intersect(tris, *intersect_ptr);
      CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "PRESSURE_FE");
    }
  }
 
  MoFEMFunctionReturn(0);
}

◆ addNeumannFluxBCElements()

MoFEMErrorCode MetaNeumannForces::addNeumannFluxBCElements	(	MoFEM::Interface &	m_field,
		const std::string	field_name,
		const std::string	mesh_nodals_positions = `"MESH_NODE_POSITIONS"`
	)

static

Definition at line 2143 of file SurfacePressure.cpp.

                                           {
  MoFEMFunctionBegin;
 
  CHKERR m_field.add_finite_element("FLUX_FE", MF_ZERO);
  CHKERR m_field.modify_finite_element_add_field_row("FLUX_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_col("FLUX_FE", field_name);
  CHKERR m_field.modify_finite_element_add_field_data("FLUX_FE", field_name);
  if (m_field.check_field(mesh_nodals_positions)) {
    CHKERR m_field.modify_finite_element_add_field_data("FLUX_FE",
                                                        mesh_nodals_positions);
  }
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    Range tris;
    CHKERR m_field.get_moab().get_entities_by_dimension(it->meshset, 2, tris,
                                                        true);
    CHKERR m_field.add_ents_to_finite_element_by_dim(tris, 2, "FLUX_FE");
  }
 
  MoFEMFunctionReturn(0);
}

◆ setMassFluxOperators()

MoFEMErrorCode MetaNeumannForces::setMassFluxOperators	(	MoFEM::Interface &	m_field,
		boost::ptr_map< std::string, NeumannForcesSurface > &	neumann_forces,
		Vec	F,
		const std::string	field_name,
		const std::string	mesh_nodals_positions = `"MESH_NODE_POSITIONS"`
	)

static

Definition at line 2168 of file SurfacePressure.cpp.

                                                                       {
  MoFEMFunctionBegin;
  bool ho_geometry = m_field.check_field(mesh_nodals_positions);
 
  string fe_name;
  fe_name = "FLUX_FE";
  neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
 
  auto &nf = neumann_forces.at(fe_name);
  auto &fe = nf.getLoopFe();
 
  if (ho_geometry)
    CHKERR AddHOOps<2, 3, 3>::add(fe.getOpPtrVector(), {},
                                  mesh_nodals_positions);
 
  for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field,
                                                  SIDESET | PRESSURESET, it)) {
    CHKERR nf.addFlux(field_name, F, it->getMeshsetId(), ho_geometry);
  }
  MoFEMFunctionReturn(0);
}

◆ setMomentumFluxOperators()

MoFEMErrorCode MetaNeumannForces::setMomentumFluxOperators	(	MoFEM::Interface &	m_field,
		boost::ptr_map< std::string, NeumannForcesSurface > &	neumann_forces,
		Vec	F,
		const std::string	field_name,
		const std::string	mesh_nodals_positions = `"MESH_NODE_POSITIONS"`
	)

static

Set operators to finite elements calculating right hand side vector.

Parameters

m_field	Interface
neumann_forces	Map of pointers to force/pressure elements
F	Right hand side vector
field_name	Field name (f.e. DISPLACEMENT)
mesh_nodals_positions	Name of field on which ho-geometry is defined

Returns: Error code

Examples: mofem/tutorials/cor-10_navier_stokes/navier_stokes.cpp, mofem/tutorials/cor-7_elasticity_mixed_formulation/elasticity_mixed_formulation.cpp, and mofem/users_modules/basic_finite_elements/elasticity/elasticity.cpp.

Definition at line 2069 of file SurfacePressure.cpp.

                                                                       {
  MoFEMFunctionBegin;
  bool ho_geometry = m_field.check_field(mesh_nodals_positions);
  // Add forces
  {
    std::string fe_name = "FORCE_FE";
    neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
    auto &nf = neumann_forces.at(fe_name);
    auto &fe = nf.getLoopFe();
 
    if (ho_geometry)
      CHKERR AddHOOps<2, 3, 3>::add(fe.getOpPtrVector(), {},
                                    mesh_nodals_positions);
 
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(m_field, NODESET | FORCESET,
                                                    it)) {
      CHKERR nf.addForce(field_name, F, it->getMeshsetId(), ho_geometry, false);
    }
    // Reading forces from BLOCKSET
    const string block_set_force_name("FORCE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_force_name.length(),
                                block_set_force_name) == 0) {
        CHKERR nf.addForce(field_name, F, it->getMeshsetId(), ho_geometry,
                           true);
      }
    }
  }
 
  // Add pressures
  {
    std::string fe_name = "PRESSURE_FE";
    neumann_forces.insert(fe_name, new NeumannForcesSurface(m_field));
 
    auto &nf = neumann_forces.at(fe_name);
    auto &fe = nf.getLoopFe();
 
    if (ho_geometry)
      CHKERR AddHOOps<2, 3, 3>::add(fe.getOpPtrVector(), {},
                                    mesh_nodals_positions);
 
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             m_field, SIDESET | PRESSURESET, it)) {
      CHKERR nf.addPressure(field_name, F, it->getMeshsetId(), ho_geometry,
                            false);
    }
 
    // Reading pressures from BLOCKSET
    const string block_set_pressure_name("PRESSURE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_pressure_name.length(),
                                block_set_pressure_name) == 0) {
        CHKERR nf.addPressure(field_name, F, it->getMeshsetId(), ho_geometry,
                              true);
      }
    }
 
    // Reading pressures from BLOCKSET
    const string block_set_linear_pressure_name("LINEAR_PRESSURE");
    for (_IT_CUBITMESHSETS_BY_SET_TYPE_FOR_LOOP_(m_field, BLOCKSET, it)) {
      if (it->getName().compare(0, block_set_linear_pressure_name.length(),
                                block_set_linear_pressure_name) == 0) {
        CHKERR nf.addLinearPressure(field_name, F, it->getMeshsetId(),
                                    ho_geometry);
      }
    }
  }
 
  MoFEMFunctionReturn(0);
}

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

users_modules/basic_finite_elements/src/SurfacePressure.hpp
users_modules/basic_finite_elements/src/impl/SurfacePressure.cpp

Static Public Member Functions

Detailed Description

Member Function Documentation

◆ addNeumannBCElements()

◆ addNeumannFluxBCElements()

◆ setMassFluxOperators()

◆ setMomentumFluxOperators()