seepage.cpp File Reference

#include <MoFEM.hpp>
#include <MatrixFunction.hpp>
#include <IntegrationRules.hpp>
#include <SeepageOps.hpp>

Go to the source code of this file.

Classes
struct	Seepage
struct	Seepage::BlockedParameters

Macros
#define	EXECUTABLE_DIMENSION   2

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	BoundaryEle
using	PostProcEle = PostProcBrokenMeshInMoab<DomainEle>
using	AssemblyDomainEleOp
using	OpKCauchy
	[Only used with Hooke equation (linear material model)]
using	OpInternalForceCauchy
using	OpMass
	[Only used with Hooke equation (linear material model)]
using	OpInertiaForce
using	OpKPiola
	[Only used for dynamics]
using	OpInternalForcePiola
using	OpBoundaryMass
	[Only used with Hencky/nonlinear material]
using	OpBoundaryVec
using	OpBoundaryInternal
using	OpBaseDivU
	[Essential boundary conditions]
using	OpHdivHdiv
	Integrate Lhs base of flux (1/k) base of flux (FLUX x FLUX)
using	OpHdivQ
	Integrate Lhs div of base of flux time base of temperature (FLUX x T) and transpose of it, i.e. (T x FLAX)
using	OpCapacity
	Integrate Lhs base of temperature times (heat capacity) times base of temperature (T x T)
using	OpHdivFlux
	Integrating Rhs flux base (1/k) flux (FLUX)
using	OpHDivH
	Integrate Rhs div flux base times temperature (T)
using	OpBaseDotH
	Integrate Rhs base of temperature time heat capacity times heat rate (T)
using	OpBaseDivFlux = OpBaseDotH
	Integrate Rhs base of temperature times divergent of flux (T)
using	DomainNaturalBCRhs
	[Body and heat source]
using	OpBodyForce
using	OpHeatSource
using	DomainNaturalBCLhs
using	BoundaryNaturalBC
	[Body and heat source]
using	OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>
using	OpTemperatureBC
using	OpEssentialFluxRhs
	[Natural boundary conditions]
using	OpEssentialFluxLhs

Functions
int	main (int argc, char *argv[])

Variables
constexpr int	SPACE_DIM
double	scale = 1.
	[Essential boundary conditions (Least square approach)]
double	default_young_modulus = 1
double	default_poisson_ratio = 0.25
double	default_conductivity = 1
double	fluid_density = 10
auto	save_range
static char	help [] = "...\n\n"
	[Solve]

Macro Definition Documentation

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION   2

Definition at line 24 of file seepage.cpp.

Typedef Documentation

AssemblyDomainEleOp

using AssemblyDomainEleOp

Initial value:

 
    FormsIntegrators<DomainEleOp>::Assembly<PETSC>::OpBase

Definition at line 45 of file seepage.cpp.

BoundaryEle

using BoundaryEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 40 of file seepage.cpp.

BoundaryNaturalBC

using BoundaryNaturalBC

Initial value:

 
    NaturalBC<BoundaryEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

[Body and heat source]

[Natural boundary conditions]

Definition at line 144 of file seepage.cpp.

DomainEle

using DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle

Definition at line 38 of file seepage.cpp.

DomainNaturalBCLhs

using DomainNaturalBCLhs

Initial value:

 
    NaturalBC<DomainEleOp>::Assembly<PETSC>::BiLinearForm<GAUSS>

Definition at line 139 of file seepage.cpp.

DomainNaturalBCRhs

using DomainNaturalBCRhs

Initial value:

 
    NaturalBC<DomainEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

[Body and heat source]

Definition at line 133 of file seepage.cpp.

OpBaseDivFlux

using OpBaseDivFlux = OpBaseDotH

Integrate Rhs base of temperature times divergent of flux (T)

Examples

seepage.cpp.

Definition at line 130 of file seepage.cpp.

OpBaseDivU

using OpBaseDivU

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpMixScalarTimesDiv<SPACE_DIM>

[Essential boundary conditions]

Examples

seepage.cpp.

Definition at line 78 of file seepage.cpp.

OpBaseDotH

using OpBaseDotH

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm<
    GAUSS>::OpBaseTimesScalarField<1>

Integrate Rhs base of temperature time heat capacity times heat rate (T)

Examples

seepage.cpp.

Definition at line 123 of file seepage.cpp.

OpBodyForce

using OpBodyForce

Initial value:

 
    DomainNaturalBCRhs::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, SPACE_DIM>

Definition at line 135 of file seepage.cpp.

OpBoundaryInternal

using OpBoundaryInternal

Initial value:

 FormsIntegrators<BoundaryEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 1>

Definition at line 74 of file seepage.cpp.

OpBoundaryMass

using OpBoundaryMass

Initial value:

 FormsIntegrators<BoundaryEleOp>::Assembly<
    PETSC>::BiLinearForm<GAUSS>::OpMass<1, SPACE_DIM>

[Only used with Hencky/nonlinear material]

[Essential boundary conditions]

Examples

helmholtz.cpp.

Definition at line 70 of file seepage.cpp.

OpBoundaryVec

using OpBoundaryVec

Initial value:

 FormsIntegrators<BoundaryEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 0>

Definition at line 72 of file seepage.cpp.

OpCapacity

using OpCapacity

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpMass<1, 1>

Integrate Lhs base of temperature times (heat capacity) times base of temperature (T x T)

Definition at line 102 of file seepage.cpp.

OpEssentialFluxLhs

using OpEssentialFluxLhs

Initial value:

 
    EssentialBC<BoundaryEleOp>::Assembly<PETSC>::BiLinearForm<
        GAUSS>::OpEssentialLhs<HeatFluxCubitBcData, 3, SPACE_DIM>

Definition at line 155 of file seepage.cpp.

OpEssentialFluxRhs

using OpEssentialFluxRhs

Initial value:

 
    EssentialBC<BoundaryEleOp>::Assembly<PETSC>::LinearForm<
        GAUSS>::OpEssentialRhs<HeatFluxCubitBcData, 3, SPACE_DIM>

[Natural boundary conditions]

[Essential boundary conditions (Least square approach)]

Definition at line 152 of file seepage.cpp.

OpForce

using OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>

Definition at line 146 of file seepage.cpp.

OpHdivFlux

using OpHdivFlux

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm<
    GAUSS>::OpBaseTimesVector<3, 3, 1>

Integrating Rhs flux base (1/k) flux (FLUX)

Examples

seepage.cpp.

Definition at line 108 of file seepage.cpp.

OpHDivH

using OpHDivH

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm<
    GAUSS>::OpMixDivTimesU<3, 1, 2>

Integrate Rhs div flux base times temperature (T)

Examples

seepage.cpp, and test_broken_space.cpp.

Definition at line 115 of file seepage.cpp.

OpHdivHdiv

using OpHdivHdiv

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpMass<3, 3>

Integrate Lhs base of flux (1/k) base of flux (FLUX x FLUX)

Examples

seepage.cpp.

Definition at line 86 of file seepage.cpp.

OpHdivQ

using OpHdivQ

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpMixDivTimesScalar<SPACE_DIM>

Integrate Lhs div of base of flux time base of temperature (FLUX x T) and transpose of it, i.e. (T x FLAX)

Examples

seepage.cpp.

Definition at line 94 of file seepage.cpp.

OpHeatSource

using OpHeatSource

Initial value:

 
    DomainNaturalBCRhs::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, 1>

Definition at line 137 of file seepage.cpp.

OpInertiaForce

using OpInertiaForce

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 1>

Definition at line 58 of file seepage.cpp.

OpInternalForceCauchy

using OpInternalForceCauchy

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpGradTimesSymTensor<1, SPACE_DIM, SPACE_DIM>

Examples

seepage.cpp.

Definition at line 51 of file seepage.cpp.

OpInternalForcePiola

using OpInternalForcePiola

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpGradTimesTensor<1, SPACE_DIM, SPACE_DIM>

Examples

HenckyOps.hpp, PlasticOps.hpp, and adolc_plasticity.cpp.

Definition at line 65 of file seepage.cpp.

OpKCauchy

using OpKCauchy

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpGradSymTensorGrad<1, SPACE_DIM, SPACE_DIM, 0>

[Only used with Hooke equation (linear material model)]

Examples

seepage.cpp.

Definition at line 49 of file seepage.cpp.

OpKPiola

using OpKPiola

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpGradTensorGrad<1, SPACE_DIM, SPACE_DIM, 1>

[Only used for dynamics]

[Only used with Hencky/nonlinear material]

Examples

HenckyOps.hpp, PlasticOps.hpp, and adolc_plasticity.cpp.

Definition at line 63 of file seepage.cpp.

OpMass

using OpMass

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm<
    GAUSS>::OpMass<1, SPACE_DIM>

[Only used with Hooke equation (linear material model)]

[Only used for dynamics]

Examples

eigen_elastic.cpp, heat_method.cpp, plastic.cpp, and scalar_check_approximation.cpp.

Definition at line 56 of file seepage.cpp.

OpTemperatureBC

using OpTemperatureBC

Initial value:

 
    BoundaryNaturalBC::OpFlux<NaturalTemperatureMeshsets, 3, SPACE_DIM>

Definition at line 147 of file seepage.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 43 of file seepage.cpp.

Function Documentation

main()

int main	(	int	argc,
		char *	argv[] )

[Register MoFEM discrete manager in PETSc]

[Register MoFEM discrete manager in PETSc

[Create MoAB]

< mesh database

< mesh database interface

[Create MoAB]

[Create MoFEM]

< finite element database

< finite element database insterface

[Create MoFEM]

[Load mesh]

[Load mesh]

[Seepage]

[Seepage]

Definition at line 1110 of file seepage.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for example
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "Seepage"));
  LogManager::setLog("Seepage");
  MOFEM_LOG_TAG("Seepage", "seepage");
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmSync(), "SeepageSync"));
  LogManager::setLog("SeepageSync");
  MOFEM_LOG_TAG("SeepageSync", "seepage");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    //! [Register MoFEM discrete manager in PETSc
 
    //! [Create MoAB]
    moab::Core mb_instance;              ///< mesh database
    moab::Interface &moab = mb_instance; ///< mesh database interface
    //! [Create MoAB]
 
    //! [Create MoFEM]
    MoFEM::Core core(moab);           ///< finite element database
    MoFEM::Interface &m_field = core; ///< finite element database insterface
    //! [Create MoFEM]
 
    //! [Load mesh]
    Simple *simple = m_field.getInterface<Simple>();
    CHKERR simple->getOptions();
    CHKERR simple->loadFile();
    //! [Load mesh]
 
    //! [Seepage]
    Seepage ex(m_field);
    CHKERR ex.runProblem();
    //! [Seepage]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

default_conductivity

double default_conductivity = 1

Examples

seepage.cpp.

Definition at line 164 of file seepage.cpp.

default_poisson_ratio

double default_poisson_ratio = 0.25

Examples

seepage.cpp.

Definition at line 163 of file seepage.cpp.

default_young_modulus

double default_young_modulus = 1

Examples

seepage.cpp.

Definition at line 162 of file seepage.cpp.

fluid_density

double fluid_density = 10

Examples

seepage.cpp.

Definition at line 165 of file seepage.cpp.

help

char help[] = "...\n\n"

static

[Solve]

Definition at line 1108 of file seepage.cpp.

save_range

auto save_range

Initial value:

= [](moab::Interface &moab, const std::string name,
                     const Range r) {
  MoFEMFunctionBegin;
  auto out_meshset = get_temp_meshset_ptr(moab);
  CHKERR moab.add_entities(*out_meshset, r);
  CHKERR moab.write_file(name.c_str(), "VTK", "", out_meshset->get_ptr(), 1);
  MoFEMFunctionReturn(0);
}

Examples

seepage.cpp.

Definition at line 170 of file seepage.cpp.

                                    {
  MoFEMFunctionBegin;
  auto out_meshset = get_temp_meshset_ptr(moab);
  CHKERR moab.add_entities(*out_meshset, r);
  CHKERR moab.write_file(name.c_str(), "VTK", "", out_meshset->get_ptr(), 1);
  MoFEMFunctionReturn(0);
};

scale

double scale = 1.

[Essential boundary conditions (Least square approach)]

Definition at line 160 of file seepage.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Examples

seepage.cpp.

Definition at line 34 of file seepage.cpp.