free_surface.cpp File Reference

#include <MoFEM.hpp>
#include <petsc/private/petscimpl.h>
#include <FreeSurfaceOps.hpp>

Go to the source code of this file.

Classes
struct	TSPrePostProc
	Set of functions called by PETSc solver used to refine and update mesh. More...
struct	FreeSurface
struct	Monitor
	[Push operators to pipeline] More...

Typedefs
template<int DIM>
using	ElementsAndOps = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >
using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle
using	DomianParentEle = ElementsAndOps< SPACE_DIM >::DomianParentEle
using	DomainEleOp = DomainEle::UserDataOperator
using	BoundaryEle = ElementsAndOps< SPACE_DIM >::BoundaryEle
using	BoundaryParentEle = ElementsAndOps< SPACE_DIM >::BoundaryParentEle
using	BoundaryEleOp = BoundaryEle::UserDataOperator
using	SideEle = ElementsAndOps< SPACE_DIM >::FaceSideEle
using	SideOp = SideEle::UserDataOperator
using	PostProcEleDomain = PostProcBrokenMeshInMoab< DomainEle >
using	PostProcEleDomainCont = PostProcBrokenMeshInMoabBaseCont< DomainEle >
using	PostProcEleBdyCont = PostProcBrokenMeshInMoabBaseCont< BoundaryEle >
using	EntData = EntitiesFieldData::EntData
using	AssemblyDomainEleOp = FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase
using	AssemblyBoundaryEleOp = FormsIntegrators< BoundaryEleOp >::Assembly< A >::OpBase
using	OpDomainMassU = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< I >::OpMass< BASE_DIM, U_FIELD_DIM >
using	OpDomainMassH = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< I >::OpMass< BASE_DIM, 1 >
using	OpDomainMassP = OpDomainMassH
using	OpDomainMassG = OpDomainMassH
using	OpBoundaryMassL = FormsIntegrators< BoundaryEleOp >::Assembly< A >::BiLinearForm< I >::OpMass< BASE_DIM, 1 >
using	OpDomainAssembleVector = FormsIntegrators< DomainEleOp >::Assembly< A >::LinearForm< I >::OpBaseTimesVector< BASE_DIM, SPACE_DIM, 1 >
using	OpDomainAssembleScalar = FormsIntegrators< DomainEleOp >::Assembly< A >::LinearForm< I >::OpBaseTimesScalar< BASE_DIM >
using	OpBoundaryAssembleScalar = FormsIntegrators< BoundaryEleOp >::Assembly< A >::LinearForm< I >::OpBaseTimesScalar< BASE_DIM >
template<CoordinateTypes COORD_TYPE>
using	OpMixScalarTimesDiv = FormsIntegrators< DomainEleOp >::Assembly< A >::BiLinearForm< I >::OpMixScalarTimesDiv< SPACE_DIM, COORD_TYPE >
using	BoundaryNaturalBC = NaturalBC< BoundaryEleOp >::Assembly< A >::LinearForm< I >
using	OpFluidFlux = BoundaryNaturalBC::OpFlux< NaturalMeshsetType< BLOCKSET >, 1, 1 >

Enumerations
enum	FR { F , R }

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

Variables
static char	help [] = "...\n\n"
int	coord_type = EXECUTABLE_COORD_TYPE
constexpr int	BASE_DIM = 1
constexpr int	SPACE_DIM = 2
constexpr int	U_FIELD_DIM = SPACE_DIM
constexpr AssemblyType	A = AssemblyType::PETSC
constexpr IntegrationType	I
FTensor::Index< 'i', SPACE_DIM >	i
FTensor::Index< 'j', SPACE_DIM >	j
FTensor::Index< 'k', SPACE_DIM >	k
FTensor::Index< 'l', SPACE_DIM >	l
constexpr auto	t_kd = FTensor::Kronecker_Delta_symmetric<int>()
int	order = 3
	approximation order
int	nb_levels = 4
int	refine_overlap = 4
constexpr bool	debug = true
auto	get_start_bit
auto	get_current_bit
	dofs bit used to do calculations
auto	get_skin_parent_bit = []() { return 2 * get_start_bit() + 2; }
auto	get_skin_child_bit = []() { return 2 * get_start_bit() + 3; }
auto	get_projection_bit = []() { return 2 * get_start_bit() + 4; }
auto	get_skin_projection_bit = []() { return 2 * get_start_bit() + 5; }
double	a0 = 980
double	rho_m = 0.998
double	mu_m = 0.010101 * 1e2
double	rho_p = 0.0012
double	mu_p = 0.000182 * 1e2
double	lambda = 73
	surface tension
double	W = 0.25
double	h = 0.03
double	eta = h
double	eta2 = eta * eta
double	md = 1e-2
double	eps = 1e-12
double	tol = std::numeric_limits<float>::epsilon()
double	rho_ave = (rho_p + rho_m) / 2
double	rho_diff = (rho_p - rho_m) / 2
double	mu_ave = (mu_p + mu_m) / 2
double	mu_diff = (mu_p - mu_m) / 2
double	kappa = (3. / (4. * std::sqrt(2. * W))) * (lambda / eta)
auto	integration_rule = [](int, int, int) { return 2 * order + 1; }
auto	cylindrical
auto	wetting_angle_sub_stepping
auto	my_max = [](const double x) { return (x - 1 + std::abs(x + 1)) / 2; }
auto	my_min = [](const double x) { return (x + 1 - std::abs(x - 1)) / 2; }
auto	cut_off = [](const double h) { return my_max(my_min(h)); }
auto	d_cut_off
auto	phase_function
auto	d_phase_function_h
auto	get_f = [](const double h) { return 4 * W * h * (h * h - 1); }
auto	get_f_dh = [](const double h) { return 4 * W * (3 * h * h - 1); }
auto	get_M0 = [](auto h) { return md; }
auto	get_M0_dh = [](auto h) { return 0; }
auto	get_M2
auto	get_M2_dh = [](auto h) { return -md * 2 * h; }
auto	get_M3
auto	get_M3_dh
auto	get_M = [](auto h) { return get_M0(h); }
auto	get_M_dh = [](auto h) { return get_M0_dh(h); }
auto	get_D
auto	kernel_oscillation
auto	kernel_eye
auto	capillary_tube
auto	bubble_device
auto	init_h
	Initialisation function.
auto	wetting_angle = [](double water_level) { return water_level; }
auto	bit = [](auto b) { return BitRefLevel().set(b); }
auto	marker = [](auto b) { return BitRefLevel().set(BITREFLEVEL_SIZE - b); }
auto	get_fe_bit
auto	get_global_size
auto	save_range
auto	get_dofs_ents_by_field_name
	get entities of dofs in the problem - used for debugging
auto	get_dofs_ents_all
	get entities of dofs in the problem - used for debugging
static boost::weak_ptr< TSPrePostProc >	tsPrePostProc

Typedef Documentation

AssemblyBoundaryEleOp

using AssemblyBoundaryEleOp = FormsIntegrators<BoundaryEleOp>::Assembly<A>::OpBase

Definition at line 106 of file free_surface.cpp.

AssemblyDomainEleOp

using AssemblyDomainEleOp = FormsIntegrators<DomainEleOp>::Assembly<A>::OpBase

Definition at line 105 of file free_surface.cpp.

BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 93 of file free_surface.cpp.

BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 95 of file free_surface.cpp.

BoundaryNaturalBC

using BoundaryNaturalBC = NaturalBC<BoundaryEleOp>::Assembly<A>::LinearForm<I>

Definition at line 130 of file free_surface.cpp.

BoundaryParentEle

using BoundaryParentEle = ElementsAndOps<SPACE_DIM>::BoundaryParentEle

Examples

free_surface.cpp.

Definition at line 94 of file free_surface.cpp.

DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 90 of file free_surface.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 92 of file free_surface.cpp.

DomianParentEle

using DomianParentEle = ElementsAndOps<SPACE_DIM>::DomianParentEle

Examples

free_surface.cpp.

Definition at line 91 of file free_surface.cpp.

ElementsAndOps

template<int DIM>

using ElementsAndOps = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>

Definition at line 88 of file free_surface.cpp.

EntData

using EntData = EntitiesFieldData::EntData

Definition at line 103 of file free_surface.cpp.

OpBoundaryAssembleScalar

using OpBoundaryAssembleScalar = FormsIntegrators<BoundaryEleOp>::Assembly< A>::LinearForm<I>::OpBaseTimesScalar<BASE_DIM>

Examples

free_surface.cpp.

Definition at line 122 of file free_surface.cpp.

OpBoundaryMassL

using OpBoundaryMassL = FormsIntegrators<BoundaryEleOp>::Assembly< A>::BiLinearForm<I>::OpMass<BASE_DIM, 1>

Examples

free_surface.cpp.

Definition at line 115 of file free_surface.cpp.

OpDomainAssembleScalar

using OpDomainAssembleScalar = FormsIntegrators<DomainEleOp>::Assembly< A>::LinearForm<I>::OpBaseTimesScalar<BASE_DIM>

Examples

free_surface.cpp.

Definition at line 120 of file free_surface.cpp.

OpDomainAssembleVector

using OpDomainAssembleVector = FormsIntegrators<DomainEleOp>::Assembly< A>::LinearForm<I>::OpBaseTimesVector<BASE_DIM, SPACE_DIM, 1>

Examples

free_surface.cpp.

Definition at line 118 of file free_surface.cpp.

OpDomainMassG

using OpDomainMassG = OpDomainMassH

Examples

free_surface.cpp.

Definition at line 114 of file free_surface.cpp.

OpDomainMassH

using OpDomainMassH = FormsIntegrators<DomainEleOp>::Assembly<A>::BiLinearForm< I>::OpMass<BASE_DIM, 1>

Examples

free_surface.cpp.

Definition at line 111 of file free_surface.cpp.

OpDomainMassP

using OpDomainMassP = OpDomainMassH

Examples

free_surface.cpp.

Definition at line 113 of file free_surface.cpp.

OpDomainMassU

using OpDomainMassU = FormsIntegrators<DomainEleOp>::Assembly<A>::BiLinearForm< I>::OpMass<BASE_DIM, U_FIELD_DIM>

Examples

free_surface.cpp.

Definition at line 109 of file free_surface.cpp.

OpFluidFlux

using OpFluidFlux = BoundaryNaturalBC::OpFlux<NaturalMeshsetType<BLOCKSET>, 1, 1>

Definition at line 131 of file free_surface.cpp.

OpMixScalarTimesDiv

template<CoordinateTypes COORD_TYPE>

using OpMixScalarTimesDiv = FormsIntegrators<DomainEleOp>::Assembly< A>::BiLinearForm<I>::OpMixScalarTimesDiv<SPACE_DIM, COORD_TYPE>

Definition at line 126 of file free_surface.cpp.

PostProcEleBdyCont

using PostProcEleBdyCont = PostProcBrokenMeshInMoabBaseCont<BoundaryEle>

Definition at line 101 of file free_surface.cpp.

PostProcEleDomain

using PostProcEleDomain = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 99 of file free_surface.cpp.

PostProcEleDomainCont

using PostProcEleDomainCont = PostProcBrokenMeshInMoabBaseCont<DomainEle>

Definition at line 100 of file free_surface.cpp.

SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle

Definition at line 96 of file free_surface.cpp.

SideOp

using SideOp = SideEle::UserDataOperator

Definition at line 97 of file free_surface.cpp.

Enumeration Type Documentation

FR

enum FR

Enumerator
F
R

Definition at line 396 of file free_surface.cpp.

{ F, R }; // F - forward, and reverse

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]

[FreeSurface]

[FreeSurface]

Definition at line 2403 of file free_surface.cpp.

                                 {
 
#ifdef PYTHON_INIT_SURFACE
  Py_Initialize();
#endif
 
  // 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(), "FS"));
  LogManager::setLog("FS");
  MOFEM_LOG_TAG("FS", "free surface");
 
  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]
 
    //! [FreeSurface]
    FreeSurface ex(m_field);
    CHKERR ex.runProblem();
    //! [FreeSurface]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
 
#ifdef PYTHON_INIT_SURFACE
  if (Py_FinalizeEx() < 0) {
    exit(120);
  }
#endif
}

Variable Documentation

A

constexpr AssemblyType A = AssemblyType::PETSC

constexpr

Examples

free_surface.cpp.

Definition at line 83 of file free_surface.cpp.

a0

double a0 = 980

Examples

free_surface.cpp.

Definition at line 162 of file free_surface.cpp.

BASE_DIM

constexpr int BASE_DIM = 1

constexpr

Definition at line 79 of file free_surface.cpp.

bit

auto bit = [](auto b) { return BitRefLevel().set(b); }

Examples

free_surface.cpp.

Definition at line 314 of file free_surface.cpp.

{ return BitRefLevel().set(b); };

bubble_device

auto bubble_device

Initial value:

= [](double x, double y, double z) {
  return -tanh((-0.039 - x) / (eta * std::sqrt(2)));
}

Examples

free_surface.cpp.

Definition at line 286 of file free_surface.cpp.

                                                      {
  return -tanh((-0.039 - x) / (eta * std::sqrt(2)));
};

capillary_tube

auto capillary_tube

Initial value:

= [](double x, double y, double z) {
  constexpr double water_height = 0.;
  return tanh((water_height - y) / (eta * std::sqrt(2)));
  ;
}

Examples

free_surface.cpp.

Definition at line 280 of file free_surface.cpp.

                                                       {
  constexpr double water_height = 0.;
  return tanh((water_height - y) / (eta * std::sqrt(2)));
  ;
};

coord_type

int coord_type = EXECUTABLE_COORD_TYPE

Examples

free_surface.cpp.

Definition at line 77 of file free_surface.cpp.

cut_off

auto cut_off = [](const double h) { return my_max(my_min(h)); }

Examples

free_surface.cpp.

Definition at line 205 of file free_surface.cpp.

{ return my_max(my_min(h)); };

cylindrical

auto cylindrical

Initial value:

= [](const double r) {
  if (coord_type == CYLINDRICAL)
    return 2 * M_PI * r;
  else
    return 1.;
}

Examples

free_surface.cpp.

Definition at line 189 of file free_surface.cpp.

                                      {
  if (coord_type == CYLINDRICAL)
    return 2 * M_PI * r;
  else
    return 1.;
};

d_cut_off

auto d_cut_off

Initial value:

= [](const double h) {
  if (h >= -1 && h < 1)
    return 1.;
  else
    return 0.;
}

Examples

free_surface.cpp.

Definition at line 206 of file free_surface.cpp.

                                    {
  if (h >= -1 && h < 1)
    return 1.;
  else
    return 0.;
};

d_phase_function_h

auto d_phase_function_h

Initial value:

= [](const double h, const double diff) {
  return diff;
}

Examples

free_surface.cpp.

Definition at line 217 of file free_surface.cpp.

                                                                {
  return diff;
};

debug

constexpr bool debug = true

constexpr

Definition at line 146 of file free_surface.cpp.

eps

double eps = 1e-12

Examples

free_surface.cpp.

Definition at line 177 of file free_surface.cpp.

eta

double eta = h

Examples

free_surface.cpp, and plot_base.cpp.

Definition at line 172 of file free_surface.cpp.

eta2

double eta2 = eta * eta

Examples

free_surface.cpp.

Definition at line 173 of file free_surface.cpp.

get_current_bit

auto get_current_bit

Initial value:

= []() {
  return 2 * get_start_bit() + 1;
}

dofs bit used to do calculations

Examples

free_surface.cpp.

Definition at line 151 of file free_surface.cpp.

                            {
  return 2 * get_start_bit() + 1;
}; ///< dofs bit used to do calculations

get_D

auto get_D

Initial value:

= [](const double A) {
  FTensor::Ddg<double, SPACE_DIM, SPACE_DIM> t_D;
  t_D(i, j, k, l) = A * ((t_kd(i, k) ^ t_kd(j, l)) / 4.);
  return t_D;
}

Examples

free_surface.cpp.

Definition at line 254 of file free_surface.cpp.

                                {
  FTensor::Ddg<double, SPACE_DIM, SPACE_DIM> t_D;
  t_D(i, j, k, l) = A * ((t_kd(i, k) ^ t_kd(j, l)) / 4.);
  return t_D;
};

get_dofs_ents_all

auto get_dofs_ents_all

Initial value:

= [](auto dm) {
  auto prb_ptr = getProblemPtr(dm);
  std::vector<EntityHandle> ents_vec;
 
  MoFEM::Interface *m_field_ptr;
  CHKERR DMoFEMGetInterfacePtr(dm, &m_field_ptr);
 
  auto dofs = prb_ptr->numeredRowDofsPtr;
  ents_vec.reserve(dofs->size());
 
  for (auto d : *dofs) {
    ents_vec.push_back(d->getEnt());
  }
 
  std::sort(ents_vec.begin(), ents_vec.end());
  auto it = std::unique(ents_vec.begin(), ents_vec.end());
  Range r;
  r.insert_list(ents_vec.begin(), it);
  return r;
}

get entities of dofs in the problem - used for debugging

Examples

free_surface.cpp.

Definition at line 370 of file free_surface.cpp.

                                     {
  auto prb_ptr = getProblemPtr(dm);
  std::vector<EntityHandle> ents_vec;
 
  MoFEM::Interface *m_field_ptr;
  CHKERR DMoFEMGetInterfacePtr(dm, &m_field_ptr);
 
  auto dofs = prb_ptr->numeredRowDofsPtr;
  ents_vec.reserve(dofs->size());
 
  for (auto d : *dofs) {
    ents_vec.push_back(d->getEnt());
  }
 
  std::sort(ents_vec.begin(), ents_vec.end());
  auto it = std::unique(ents_vec.begin(), ents_vec.end());
  Range r;
  r.insert_list(ents_vec.begin(), it);
  return r;
};

get_dofs_ents_by_field_name

auto get_dofs_ents_by_field_name

Initial value:

= [](auto dm, auto field_name) {
  auto prb_ptr = getProblemPtr(dm);
  std::vector<EntityHandle> ents_vec;
 
  MoFEM::Interface *m_field_ptr;
  CHKERR DMoFEMGetInterfacePtr(dm, &m_field_ptr);
 
  auto bit_number = m_field_ptr->get_field_bit_number(field_name);
  auto dofs = prb_ptr->numeredRowDofsPtr;
  auto lo_it = dofs->lower_bound(FieldEntity::getLoBitNumberUId(bit_number));
  auto hi_it = dofs->upper_bound(FieldEntity::getHiBitNumberUId(bit_number));
  ents_vec.reserve(std::distance(lo_it, hi_it));
 
  for (; lo_it != hi_it; ++lo_it) {
    ents_vec.push_back((*lo_it)->getEnt());
  }
 
  std::sort(ents_vec.begin(), ents_vec.end());
  auto it = std::unique(ents_vec.begin(), ents_vec.end());
  Range r;
  r.insert_list(ents_vec.begin(), it);
  return r;
}

get entities of dofs in the problem - used for debugging

Examples

free_surface.cpp.

Definition at line 342 of file free_surface.cpp.

                                                                {
  auto prb_ptr = getProblemPtr(dm);
  std::vector<EntityHandle> ents_vec;
 
  MoFEM::Interface *m_field_ptr;
  CHKERR DMoFEMGetInterfacePtr(dm, &m_field_ptr);
 
  auto bit_number = m_field_ptr->get_field_bit_number(field_name);
  auto dofs = prb_ptr->numeredRowDofsPtr;
  auto lo_it = dofs->lower_bound(FieldEntity::getLoBitNumberUId(bit_number));
  auto hi_it = dofs->upper_bound(FieldEntity::getHiBitNumberUId(bit_number));
  ents_vec.reserve(std::distance(lo_it, hi_it));
 
  for (; lo_it != hi_it; ++lo_it) {
    ents_vec.push_back((*lo_it)->getEnt());
  }
 
  std::sort(ents_vec.begin(), ents_vec.end());
  auto it = std::unique(ents_vec.begin(), ents_vec.end());
  Range r;
  r.insert_list(ents_vec.begin(), it);
  return r;
};

get_f

auto get_f = [](const double h) { return 4 * W * h * (h * h - 1); }

Examples

free_surface.cpp.

Definition at line 223 of file free_surface.cpp.

{ return 4 * W * h * (h * h - 1); };

get_f_dh

auto get_f_dh = [](const double h) { return 4 * W * (3 * h * h - 1); }

Examples

free_surface.cpp.

Definition at line 224 of file free_surface.cpp.

{ return 4 * W * (3 * h * h - 1); };

get_fe_bit

auto get_fe_bit

Initial value:

= [](FEMethod *fe_ptr) {
  return fe_ptr->numeredEntFiniteElementPtr->getBitRefLevel();
}

Examples

free_surface.cpp.

Definition at line 316 of file free_surface.cpp.

                                       {
  return fe_ptr->numeredEntFiniteElementPtr->getBitRefLevel();
};

get_global_size

auto get_global_size

Initial value:

= [](int l_size) {
  int g_size;
  MPI_Allreduce(&l_size, &g_size, 1, MPI_INT, MPI_SUM, MPI_COMM_WORLD);
  return g_size;
}

Examples

free_surface.cpp.

Definition at line 320 of file free_surface.cpp.

                                      {
  int g_size;
  MPI_Allreduce(&l_size, &g_size, 1, MPI_INT, MPI_SUM, MPI_COMM_WORLD);
  return g_size;
};

get_M

auto get_M = [](auto h) { return get_M0(h); }

Examples

free_surface.cpp.

Definition at line 251 of file free_surface.cpp.

{ return get_M0(h); };

get_M0

auto get_M0 = [](auto h) { return md; }

Examples

free_surface.cpp.

Definition at line 226 of file free_surface.cpp.

{ return md; };

get_M0_dh

auto get_M0_dh = [](auto h) { return 0; }

Examples

free_surface.cpp.

Definition at line 227 of file free_surface.cpp.

{ return 0; };

get_M2

auto get_M2

Initial value:

= [](auto h) {
  return md * (1 - h * h);
}

Examples

free_surface.cpp.

Definition at line 229 of file free_surface.cpp.

                         {
  return md * (1 - h * h);
};

get_M2_dh

auto get_M2_dh = [](auto h) { return -md * 2 * h; }

Examples

free_surface.cpp.

Definition at line 233 of file free_surface.cpp.

{ return -md * 2 * h; };

get_M3

auto get_M3

Initial value:

= [](auto h) {
  const double h2 = h * h;
  const double h3 = h2 * h;
  if (h >= 0)
    return md * (2 * h3 - 3 * h2 + 1);
  else
    return md * (-2 * h3 - 3 * h2 + 1);
}

Examples

free_surface.cpp.

Definition at line 235 of file free_surface.cpp.

                         {
  const double h2 = h * h;
  const double h3 = h2 * h;
  if (h >= 0)
    return md * (2 * h3 - 3 * h2 + 1);
  else
    return md * (-2 * h3 - 3 * h2 + 1);
};

get_M3_dh

auto get_M3_dh

Initial value:

= [](auto h) {
  if (h >= 0)
    return md * (6 * h * (h - 1));
  else
    return md * (-6 * h * (h + 1));
}

Examples

free_surface.cpp.

Definition at line 244 of file free_surface.cpp.

                            {
  if (h >= 0)
    return md * (6 * h * (h - 1));
  else
    return md * (-6 * h * (h + 1));
};

get_M_dh

auto get_M_dh = [](auto h) { return get_M0_dh(h); }

Examples

free_surface.cpp.

Definition at line 252 of file free_surface.cpp.

{ return get_M0_dh(h); };

get_projection_bit

auto get_projection_bit = []() { return 2 * get_start_bit() + 4; }

Examples

free_surface.cpp.

Definition at line 156 of file free_surface.cpp.

{ return 2 * get_start_bit() + 4; };

get_skin_child_bit

auto get_skin_child_bit = []() { return 2 * get_start_bit() + 3; }

Examples

free_surface.cpp.

Definition at line 155 of file free_surface.cpp.

{ return 2 * get_start_bit() + 3; };

get_skin_parent_bit

auto get_skin_parent_bit = []() { return 2 * get_start_bit() + 2; }

Examples

free_surface.cpp.

Definition at line 154 of file free_surface.cpp.

{ return 2 * get_start_bit() + 2; };

get_skin_projection_bit

auto get_skin_projection_bit = []() { return 2 * get_start_bit() + 5; }

Examples

free_surface.cpp.

Definition at line 157 of file free_surface.cpp.

{ return 2 * get_start_bit() + 5; };

get_start_bit

auto get_start_bit

Initial value:

= []() {
  return nb_levels + 1;
}

Examples

free_surface.cpp.

Definition at line 148 of file free_surface.cpp.

                          {
  return nb_levels + 1;
}; //< first refinement level for computational mesh

h

double h = 0.03

Examples

free_surface.cpp.

Definition at line 171 of file free_surface.cpp.

help

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

static

Definition at line 17 of file free_surface.cpp.

I

constexpr IntegrationType I

constexpr

Initial value:

=
    IntegrationType::GAUSS

Definition at line 84 of file free_surface.cpp.

i

FTensor::Index<'i', SPACE_DIM> i

Definition at line 134 of file free_surface.cpp.

init_h

auto init_h

Initial value:

= [](double r, double y, double theta) {
 
 
 
 
 
 
 
 
  
  return capillary_tube(r, y, theta);
  
 
}

Initialisation function.

Note

If UMs are compiled with Python to initialise phase field "H" surface.py function is used, which has to be present in execution folder.

Examples

free_surface.cpp, and shallow_wave.cpp.

Definition at line 297 of file free_surface.cpp.

                                                   {
#ifdef PYTHON_INIT_SURFACE
  double s = 1;
  if (auto ptr = surfacePythonWeakPtr.lock()) {
    CHK_THROW_MESSAGE(ptr->evalSurface(r, y, theta, eta, s),
                      "error eval python");
  }
  return s;
#else
  // return bubble_device(r, y, theta);
  return capillary_tube(r, y, theta);
  // return kernel_eye(r, y, theta);
#endif
};

integration_rule

auto integration_rule = [](int, int, int) { return 2 * order + 1; }

Examples

ArcLengthExample.hpp, NonlinearElasticExample.hpp, approx_sphere.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, free_surface.cpp, hcurl_check_approx_in_2d.cpp, hcurl_curl_operator.cpp, hdiv_check_approx_in_3d.cpp, hdiv_divergence_operator.cpp, heat_equation.cpp, helmholtz.cpp, mofem/tutorials/vec-7/adjoint.cpp, photon_diffusion.cpp, scalar_check_approximation.cpp, seepage.cpp, shallow_wave.cpp, test_broken_space.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 187 of file free_surface.cpp.

{ return 2 * order + 1; };

j

FTensor::Index<'j', SPACE_DIM> j

Definition at line 135 of file free_surface.cpp.

k

FTensor::Index<'k', SPACE_DIM> k

Definition at line 136 of file free_surface.cpp.

kappa

double kappa = (3. / (4. * std::sqrt(2. * W))) * (lambda / eta)

Examples

free_surface.cpp.

Definition at line 185 of file free_surface.cpp.

kernel_eye

auto kernel_eye

Initial value:

= [](double r, double y, double) {
  constexpr double y0 = 0.5;
  constexpr double R = 0.4;
  y -= y0;
  const double d = std::sqrt(r * r + y * y);
  return tanh((R - d) / (eta * std::sqrt(2)));
}

Examples

free_surface.cpp.

Definition at line 272 of file free_surface.cpp.

                                                 {
  constexpr double y0 = 0.5;
  constexpr double R = 0.4;
  y -= y0;
  const double d = std::sqrt(r * r + y * y);
  return tanh((R - d) / (eta * std::sqrt(2)));
};

kernel_oscillation

auto kernel_oscillation

Initial value:

= [](double r, double y, double) {
  constexpr int n = 3;
  constexpr double R = 0.0125;
  constexpr double A = R * 0.2;
  const double theta = atan2(r, y);
  const double w = R + A * cos(n * theta);
  const double d = std::sqrt(r * r + y * y);
  return tanh((w - d) / (eta * std::sqrt(2)));
}

Examples

free_surface.cpp.

Definition at line 262 of file free_surface.cpp.

                                                         {
  constexpr int n = 3;
  constexpr double R = 0.0125;
  constexpr double A = R * 0.2;
  const double theta = atan2(r, y);
  const double w = R + A * cos(n * theta);
  const double d = std::sqrt(r * r + y * y);
  return tanh((w - d) / (eta * std::sqrt(2)));
};

l

FTensor::Index<'l', SPACE_DIM> l

Examples

free_surface.cpp.

Definition at line 137 of file free_surface.cpp.

lambda

double lambda = 73

surface tension

Examples

free_surface.cpp.

Definition at line 167 of file free_surface.cpp.

marker

auto marker = [](auto b) { return BitRefLevel().set(BITREFLEVEL_SIZE - b); }

Definition at line 315 of file free_surface.cpp.

{ return BitRefLevel().set(BITREFLEVEL_SIZE - b); };

md

double md = 1e-2

Examples

free_surface.cpp.

Definition at line 176 of file free_surface.cpp.

mu_ave

double mu_ave = (mu_p + mu_m) / 2

Examples

free_surface.cpp.

Definition at line 182 of file free_surface.cpp.

mu_diff

double mu_diff = (mu_p - mu_m) / 2

Examples

free_surface.cpp.

Definition at line 183 of file free_surface.cpp.

mu_m

double mu_m = 0.010101 * 1e2

Examples

free_surface.cpp.

Definition at line 164 of file free_surface.cpp.

mu_p

double mu_p = 0.000182 * 1e2

Examples

free_surface.cpp.

Definition at line 166 of file free_surface.cpp.

my_max

auto my_max = [](const double x) { return (x - 1 + std::abs(x + 1)) / 2; }

Examples

free_surface.cpp.

Definition at line 203 of file free_surface.cpp.

{ return (x - 1 + std::abs(x + 1)) / 2; };

my_min

auto my_min = [](const double x) { return (x + 1 - std::abs(x - 1)) / 2; }

Examples

free_surface.cpp.

Definition at line 204 of file free_surface.cpp.

{ return (x + 1 - std::abs(x - 1)) / 2; };

nb_levels

int nb_levels = 4

Examples

free_surface.cpp.

Definition at line 143 of file free_surface.cpp.

order

int order = 3

approximation order

Examples

free_surface.cpp.

Definition at line 142 of file free_surface.cpp.

phase_function

auto phase_function

Initial value:

= [](const double h, const double diff, const double ave) {
  return diff * h + ave;
}

Examples

free_surface.cpp.

Definition at line 213 of file free_surface.cpp.

                                                                              {
  return diff * h + ave;
};

refine_overlap

int refine_overlap = 4

Examples

free_surface.cpp.

Definition at line 144 of file free_surface.cpp.

rho_ave

double rho_ave = (rho_p + rho_m) / 2

Examples

free_surface.cpp.

Definition at line 180 of file free_surface.cpp.

rho_diff

double rho_diff = (rho_p - rho_m) / 2

Examples

free_surface.cpp.

Definition at line 181 of file free_surface.cpp.

rho_m

double rho_m = 0.998

Examples

free_surface.cpp.

Definition at line 163 of file free_surface.cpp.

rho_p

double rho_p = 0.0012

Examples

free_surface.cpp.

Definition at line 165 of file free_surface.cpp.

save_range

auto save_range

Initial value:

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

Examples

free_surface.cpp.

Definition at line 326 of file free_surface.cpp.

                                    {
  MoFEMFunctionBegin;
  if (get_global_size(r.size())) {
    auto out_meshset = get_temp_meshset_ptr(moab);
    CHKERR moab.add_entities(*out_meshset, r);
    CHKERR moab.write_file(name.c_str(), "MOAB", "PARALLEL=WRITE_PART",
                           out_meshset->get_ptr(), 1);
  }
  MoFEMFunctionReturn(0);
};

SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

Examples

free_surface.cpp.

Definition at line 80 of file free_surface.cpp.

t_kd

constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>()

constexpr

Examples

ADOLCPlasticityMaterialModels.hpp, EshelbianOperators.cpp, EshelbianPlasticity.cpp, HenckyOps.hpp, NeoHookean.hpp, PlasticOps.hpp, PlasticOpsGeneric.hpp, PlasticOpsLargeStrains.hpp, PlasticOpsSmallStrains.hpp, approx_sphere.cpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, elasticity.cpp, free_surface.cpp, matrix_function.cpp, plastic.cpp, plate.cpp, seepage.cpp, and shallow_wave.cpp.

Definition at line 139 of file free_surface.cpp.

tol

double tol = std::numeric_limits<float>::epsilon()

Definition at line 178 of file free_surface.cpp.

tsPrePostProc

boost::weak_ptr<TSPrePostProc> tsPrePostProc

static

Examples

dynamic_first_order_con_law.cpp, and free_surface.cpp.

Definition at line 469 of file free_surface.cpp.

U_FIELD_DIM

constexpr int U_FIELD_DIM = SPACE_DIM

constexpr

Examples

free_surface.cpp.

Definition at line 81 of file free_surface.cpp.

W

double W = 0.25

Examples

free_surface.cpp.

Definition at line 168 of file free_surface.cpp.

wetting_angle

auto wetting_angle = [](double water_level) { return water_level; }

Examples

free_surface.cpp.

Definition at line 312 of file free_surface.cpp.

{ return water_level; };

wetting_angle_sub_stepping

auto wetting_angle_sub_stepping

Initial value:

= [](auto ts_step) {
  constexpr int sub_stepping = 16;
  return std::min(1., static_cast<double>(ts_step) / sub_stepping);
}

Examples

free_surface.cpp.

Definition at line 196 of file free_surface.cpp.

                                                   {
  constexpr int sub_stepping = 16;
  return std::min(1., static_cast<double>(ts_step) / sub_stepping);
};