thermo_elastic.cpp File Reference

#include <MoFEM.hpp>
#include <ThermalConvection.hpp>
#include <ThermalRadiation.hpp>
#include <ThermoElasticOps.hpp>

Go to the source code of this file.

Classes
struct	ThermoElasticProblem
struct	ThermoElasticProblem::BlockedParameters

Macros
#define	EXECUTABLE_DIMENSION   3

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	BoundaryEle
using	PostProcEle = PostProcBrokenMeshInMoab<DomainEle>
using	SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>
using	SideEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::FaceSideEle
using	AssemblyDomainEleOp
using	OpKCauchy
	[Linear elastic problem]
using	OpInternalForceCauchy
using	OpHdivHdiv
	[Linear elastic problem]
using	OpHdivT
	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	OpHDivTemp
	Integrate Rhs div flux base times temperature (T)
using	OpBaseDotT
	Integrate Rhs base of temperature time heat capacity times heat rate (T)
using	OpBaseDivFlux = OpBaseDotT
	Integrate Rhs base of temperature times divergence of flux (T)
using	DomainNaturalBCRhs
	[Thermal problem]
using	OpBodyForce
using	OpHeatSource
using	DomainNaturalBCLhs
using	BoundaryNaturalBC
	[Body and heat source]
using	OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>
using	OpEssentialFluxRhs
	[Natural boundary conditions]
using	OpEssentialFluxLhs
using	OpSetTemperatureRhs
using	OpSetTemperatureLhs

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

Variables
constexpr int	SPACE_DIM
double	default_young_modulus = 1
	[Essential boundary conditions (Least square approach)]
double	default_poisson_ratio = 0.25
double	ref_temp = 0.0
double	init_temp = 0.0
PetscBool	is_plane_strain = PETSC_FALSE
double	default_coeff_expansion = 1
double	default_heat_conductivity
double	default_heat_capacity = 1
int	order_temp = 2
int	order_flux = 3
int	order_disp = 3
int	atom_test = 0
int	save_every = 1
PetscBool	do_output_domain
PetscBool	do_output_skin
auto	save_range
static char	help [] = "...\n\n"
	[Solve]

Macro Definition Documentation

EXECUTABLE_DIMENSION

#define EXECUTABLE_DIMENSION   3

Definition at line 10 of file thermo_elastic.cpp.

Typedef Documentation

AssemblyDomainEleOp

using AssemblyDomainEleOp

Initial value:

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

Definition at line 32 of file thermo_elastic.cpp.

BoundaryEle

using BoundaryEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 25 of file thermo_elastic.cpp.

BoundaryNaturalBC

using BoundaryNaturalBC

Initial value:

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

[Body and heat source]

[Natural boundary conditions]

Definition at line 110 of file thermo_elastic.cpp.

DomainEle

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

Definition at line 23 of file thermo_elastic.cpp.

DomainNaturalBCLhs

using DomainNaturalBCLhs

Initial value:

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

Definition at line 105 of file thermo_elastic.cpp.

DomainNaturalBCRhs

using DomainNaturalBCRhs

Initial value:

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

[Thermal problem]

[Body and heat source]

Definition at line 99 of file thermo_elastic.cpp.

OpBaseDivFlux

using OpBaseDivFlux = OpBaseDotT

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

Examples

thermo_elastic.cpp.

Definition at line 94 of file thermo_elastic.cpp.

OpBaseDotT

using OpBaseDotT

Initial value:

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

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

Examples

thermo_elastic.cpp.

Definition at line 87 of file thermo_elastic.cpp.

OpBodyForce

using OpBodyForce

Initial value:

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

Examples

dynamic_first_order_con_law.cpp.

Definition at line 101 of file thermo_elastic.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)

Examples

thermo_elastic.cpp.

Definition at line 66 of file thermo_elastic.cpp.

OpEssentialFluxLhs

using OpEssentialFluxLhs

Initial value:

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

Definition at line 120 of file thermo_elastic.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 117 of file thermo_elastic.cpp.

OpForce

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

Definition at line 112 of file thermo_elastic.cpp.

OpHdivFlux

using OpHdivFlux

Initial value:

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

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

Examples

test_broken_space.cpp, and thermo_elastic.cpp.

Definition at line 72 of file thermo_elastic.cpp.

OpHdivHdiv

using OpHdivHdiv

Initial value:

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

[Linear elastic problem]

[Thermal problem]

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

Examples

phase.cpp, test_broken_space.cpp, and thermo_elastic.cpp.

Definition at line 50 of file thermo_elastic.cpp.

OpHdivT

using OpHdivT

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

thermo_elastic.cpp.

Definition at line 58 of file thermo_elastic.cpp.

OpHDivTemp

using OpHDivTemp

Initial value:

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

Integrate Rhs div flux base times temperature (T)

Examples

thermo_elastic.cpp.

Definition at line 79 of file thermo_elastic.cpp.

OpHeatSource

using OpHeatSource

Initial value:

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

Definition at line 103 of file thermo_elastic.cpp.

OpInternalForceCauchy

using OpInternalForceCauchy

Initial value:

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

Examples

PlasticOps.hpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 39 of file thermo_elastic.cpp.

OpKCauchy

using OpKCauchy

Initial value:

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

[Linear elastic problem]

Examples

PlasticOps.hpp, and thermo_elastic.cpp.

Definition at line 36 of file thermo_elastic.cpp.

OpSetTemperatureLhs

using OpSetTemperatureLhs

Initial value:

 
    DomainNaturalBCLhs::OpFlux<SetTargetTemperature, 1, 1>

Definition at line 154 of file thermo_elastic.cpp.

OpSetTemperatureRhs

using OpSetTemperatureRhs

Initial value:

 
    DomainNaturalBCRhs::OpFlux<SetTargetTemperature, 1, 1>

Definition at line 152 of file thermo_elastic.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 28 of file thermo_elastic.cpp.

SideEle

using SideEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::FaceSideEle

Definition at line 30 of file thermo_elastic.cpp.

SkinPostProcEle

using SkinPostProcEle = PostProcBrokenMeshInMoab<BoundaryEle>

Definition at line 29 of file thermo_elastic.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 interface

[Create MoFEM]

[Load mesh]

[Load mesh]

[ThermoElasticProblem]

[ThermoElasticProblem]

Definition at line 1543 of file thermo_elastic.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(), "ThermoElastic"));
  LogManager::setLog("ThermoElastic");
  MOFEM_LOG_TAG("ThermoElastic", "ThermoElastic");
 
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmSync(), "ThermoElasticSync"));
  LogManager::setLog("ThermoElasticSync");
  MOFEM_LOG_TAG("ThermoElasticSync", "ThermoElasticSync");
 
  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 interface
    //! [Create MoFEM]
 
    //! [Load mesh]
    Simple *simple = m_field.getInterface<Simple>();
    CHKERR simple->getOptions();
    CHKERR simple->loadFile();
    //! [Load mesh]
 
    //! [ThermoElasticProblem]
    ThermoElasticProblem ex(m_field);
    CHKERR ex.runProblem();
    //! [ThermoElasticProblem]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

atom_test

int atom_test = 0

Examples

thermo_elastic.cpp.

Definition at line 143 of file thermo_elastic.cpp.

default_coeff_expansion

double default_coeff_expansion = 1

Examples

thermo_elastic.cpp.

Definition at line 132 of file thermo_elastic.cpp.

default_heat_capacity

double default_heat_capacity = 1

Examples

thermo_elastic.cpp.

Definition at line 136 of file thermo_elastic.cpp.

default_heat_conductivity

double default_heat_conductivity

Initial value:

=
    1

Examples

thermo_elastic.cpp.

Definition at line 133 of file thermo_elastic.cpp.

default_poisson_ratio

double default_poisson_ratio = 0.25

Examples

seepage.cpp, and thermo_elastic.cpp.

Definition at line 126 of file thermo_elastic.cpp.

default_young_modulus

double default_young_modulus = 1

[Essential boundary conditions (Least square approach)]

Examples

seepage.cpp, and thermo_elastic.cpp.

Definition at line 125 of file thermo_elastic.cpp.

do_output_domain

PetscBool do_output_domain

Examples

thermo_elastic.cpp.

Definition at line 146 of file thermo_elastic.cpp.

do_output_skin

PetscBool do_output_skin

Examples

thermo_elastic.cpp.

Definition at line 147 of file thermo_elastic.cpp.

help

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

static

[Solve]

Definition at line 1541 of file thermo_elastic.cpp.

init_temp

double init_temp = 0.0

Examples

mortar_contact_thermal.cpp, simple_contact_thermal.cpp, and thermo_elastic.cpp.

Definition at line 128 of file thermo_elastic.cpp.

is_plane_strain

PetscBool is_plane_strain = PETSC_FALSE

Examples

thermo_elastic.cpp.

Definition at line 130 of file thermo_elastic.cpp.

order_disp

int order_disp = 3

Examples

thermo_elastic.cpp.

Definition at line 141 of file thermo_elastic.cpp.

order_flux

int order_flux = 3

Examples

thermo_elastic.cpp.

Definition at line 140 of file thermo_elastic.cpp.

order_temp

int order_temp = 2

Examples

thermo_elastic.cpp.

Definition at line 139 of file thermo_elastic.cpp.

ref_temp

double ref_temp = 0.0

Examples

ThermoElasticOps.hpp, and thermo_elastic.cpp.

Definition at line 127 of file thermo_elastic.cpp.

save_every

int save_every = 1

Examples

thermo_elastic.cpp.

Definition at line 145 of file thermo_elastic.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

EshelbianPlasticity.cpp, free_surface.cpp, level_set.cpp, seepage.cpp, and thermo_elastic.cpp.

Definition at line 157 of file thermo_elastic.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);
};

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Examples

thermo_elastic.cpp.

Definition at line 19 of file thermo_elastic.cpp.