TPZMHMixedHybridMeshControl.h

Go to the documentation of this file.

//
//  TPZMHMixedHybridMeshControl.hpp
//  PZ
//
//  Created by Philippe Devloo on 09/10/16.
//
//

#ifndef TPZMHMixedHybridMeshControl_hpp
#define TPZMHMixedHybridMeshControl_hpp

#include <stdio.h>

#include "TPZMHMixedMeshControl.h"
class TPZMultiphysicsElement;

class TPZMHMixedHybridMeshControl : public TPZMHMixedMeshControl
{
    
    

public:
    
    TPZMHMixedHybridMeshControl() : TPZMHMixedMeshControl()
    {
        
    }
    
    TPZMHMixedHybridMeshControl(int dimension) : TPZMHMixedMeshControl(dimension)
    {
        
    }

    //    TPZMHMixedHybridMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh, std::set<int64_t> &coarseindices);
    
    TPZMHMixedHybridMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh, TPZVec<int64_t> &coarseindices);
    
    TPZMHMixedHybridMeshControl(TPZAutoPointer<TPZGeoMesh> gmesh);
    
    
    TPZMHMixedHybridMeshControl(const TPZMHMixedHybridMeshControl &copy) : TPZMHMixedMeshControl(copy), fHDivWrapperMatId(copy.fHDivWrapperMatId)
    {
    }
    
    virtual ~TPZMHMixedHybridMeshControl();
    
    TPZMHMixedHybridMeshControl &operator=(const TPZMHMixedHybridMeshControl &cp)
    {
        TPZMHMixedMeshControl::operator=(cp);
        fHDivWrapperMatId = cp.fHDivWrapperMatId;
        return *this;
    }
    
    virtual void InsertPeriferalMaterialObjects();
    
    virtual void InsertPeriferalHdivMaterialObjects();
    
    virtual void InsertPeriferalPressureMaterialObjects();
    
    virtual void BuildComputationalMesh(bool usersubstructure);
    
public:
    
    int fHDivWrapperMatId = 505;
    
    int fPressureDim1MatId = 506;
    
    int fPressureDim2MatId = 507;
    
    std::set<int> fFractureFlowDim1MatId;
    
    virtual bool IsSkeletonMatid(int matid)
    {
        return matid == fSkeletonMatId || fSkeletonWithFlowMatId.find(matid) != fSkeletonWithFlowMatId.end();
    }

    int fHomogeneousNeumannBcMatId = 509;
    
    std::set<int> fSkeletonWithFlowMatId;
    
    int fSkeletonWithFlowPressureMatId = 12;
    
protected:
    // create the elements domain per domain with approximation spaces disconnected from each other
    virtual void CreateInternalFluxElements();

    // create the approximation space associated with the skeleton and restrain the connects
    virtual void CreateSkeleton();
        
    virtual void CreateMultiPhysicsInterfaceElements(int dim);


public:
    
    void InsertFractureFlowMaterial(int matid)
    {
        if (fFractureFlowDim1MatId.find(matid) != fFractureFlowDim1MatId.end()) {
            DebugStop();
        }
        fFractureFlowDim1MatId.insert(matid);
    }
    
    // create HDiv approximation spaces associated with the pressure skeleton elements
    // initially we will create only HDiv elements associated with the coarse skeleton elements
    void CreateSkeletonAxialFluxes();
    
    // create HDiv approximation spaces associated with the pressure elements on internal side
    // initially we will create only HDiv elements associated with the coarse skeleton elements
    void CreateInternalAxialFluxes();
    
    // the material id of the pressure element can be either fPressureSkeletonMatId or fPressureDim1MatId
    // An H(div) and a pressure element will be created with material Id of the geometric element
    // HDivWrapper boundary elements will also be created
    void CreateAxialFluxElement(TPZInterpolatedElement *PressureElement, int gelfluxmatid);
    
    void SplitFluxElementsAroundFractures();
    

    // if the element is neighbour of a pressure boundary condition, maintain the pressure bc
    // if the element is neighbour of a flux boundary condition, apply a homogeneous flux condition
    void AdjustBoundaryConditionsOfFractures();
    
    //  only implemented in the hybrid Hdiv version
    virtual void CheckMeshConsistency();
    
    virtual void PrintFriendly(std::ostream &out);

    virtual void SetHybridize(bool flag)
    {
        fHybridize = flag;
    }

protected:

    virtual void CreatePressureMHMMesh();
    

    // create the dim-1 pressure elements between the hdiv elements
    // this method now relies on the existence of the geometric elements
    void CreatePressureInterfaces();
    
    virtual void HybridizeSkeleton(int skeletonmatid, int pressurematid);
    
    // when dim==2 create a point pressure element for each skeleton element
    // only create internal elements - we do not hybridize the contours
    // associate an L2 do nothing material
    void CreateLowerDimensionPressureElements();
    
    // the method will look if there is a neighbouring geometric element with material id fPressureDim1MatId
    bool NeedsHDivWrapper(TPZGeoElSide gelside);
    
    // create a boundary element around each flux element to optimize calculations
    void CreateHDivWrappers();
    
    // Find the two connected flux elements to a pressure element
    void FindConnectedElements(TPZGeoElSide &pressureindex, int domain, TPZVec<TPZCompElSide> &fluxconnected);
    
    virtual void GroupandCondenseElements();

    void GroupElements(TPZCompMesh *cmesh);
    
    void ApplyNeighbourBoundaryCondition(TPZGeoEl *gel);
    
};

#endif /* TPZMHMixedHybridMeshControl_hpp */