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

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 More...

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 More...

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 More...

double	W = 0.25

double	h = 0.01 / 8

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. More...

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 More...

auto	get_dofs_ents_all
	get entities of dofs in the problem - used for debugging More...

static boost::weak_ptr< TSPrePostProc >	tsPrePostProc

Typedef Documentation

◆ AssemblyBoundaryEleOp

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

Definition at line 104 of file free_surface.cpp.

◆ AssemblyDomainEleOp

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

Definition at line 103 of file free_surface.cpp.

◆ BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 91 of file free_surface.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 93 of file free_surface.cpp.

◆ BoundaryNaturalBC

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

Definition at line 128 of file free_surface.cpp.

◆ BoundaryParentEle

using BoundaryParentEle = ElementsAndOps<SPACE_DIM>::BoundaryParentEle

Examples: free_surface.cpp.

Definition at line 92 of file free_surface.cpp.

◆ DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 88 of file free_surface.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 90 of file free_surface.cpp.

◆ DomianParentEle

using DomianParentEle = ElementsAndOps<SPACE_DIM>::DomianParentEle

Examples: free_surface.cpp.

Definition at line 89 of file free_surface.cpp.

◆ ElementsAndOps

template<int DIM>

using ElementsAndOps = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>

Definition at line 86 of file free_surface.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Definition at line 101 of file free_surface.cpp.

◆ OpBoundaryAssembleScalar

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

Examples: free_surface.cpp.

Definition at line 120 of file free_surface.cpp.

◆ OpBoundaryMassL

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

Examples: free_surface.cpp.

Definition at line 113 of file free_surface.cpp.

◆ OpDomainAssembleScalar

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

Examples: free_surface.cpp.

Definition at line 118 of file free_surface.cpp.

◆ OpDomainAssembleVector

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

Examples: free_surface.cpp.

Definition at line 116 of file free_surface.cpp.

◆ OpDomainMassG

using OpDomainMassG = OpDomainMassH

Examples: free_surface.cpp.

Definition at line 112 of file free_surface.cpp.

◆ OpDomainMassH

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

Examples: free_surface.cpp.

Definition at line 109 of file free_surface.cpp.

◆ OpDomainMassP

using OpDomainMassP = OpDomainMassH

Examples: free_surface.cpp.

Definition at line 111 of file free_surface.cpp.

◆ OpDomainMassU

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

Examples: free_surface.cpp.

Definition at line 107 of file free_surface.cpp.

◆ OpFluidFlux

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

Definition at line 129 of file free_surface.cpp.

◆ OpMixScalarTimesDiv

template<CoordinateTypes COORD_TYPE>

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

Definition at line 124 of file free_surface.cpp.

◆ PostProcEleBdyCont

using PostProcEleBdyCont = PostProcBrokenMeshInMoabBaseCont<BoundaryEle>

Definition at line 99 of file free_surface.cpp.

◆ PostProcEleDomain

using PostProcEleDomain = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 97 of file free_surface.cpp.

◆ PostProcEleDomainCont

using PostProcEleDomainCont = PostProcBrokenMeshInMoabBaseCont<DomainEle>

Definition at line 98 of file free_surface.cpp.

◆ SideEle

using SideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle

Definition at line 94 of file free_surface.cpp.

◆ SideOp

using SideOp = SideEle::UserDataOperator

Definition at line 95 of file free_surface.cpp.

Enumeration Type Documentation

◆ FR

enum FR

Enumerator
F
R

Definition at line 394 of file free_surface.cpp.

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

R

@ R

Definition: free_surface.cpp:394

F

@ F

Definition: free_surface.cpp:394

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]

Definition at line 2406 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 81 of file free_surface.cpp.

◆ a0

double a0 = 980

Examples: free_surface.cpp.

Definition at line 160 of file free_surface.cpp.

◆ BASE_DIM

constexpr int BASE_DIM = 1

constexpr

Definition at line 77 of file free_surface.cpp.

◆ bit

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

Examples: free_surface.cpp.

Definition at line 312 of file free_surface.cpp.

◆ 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 284 of file free_surface.cpp.

◆ 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 278 of file free_surface.cpp.

◆ coord_type

int coord_type = EXECUTABLE_COORD_TYPE

Examples: free_surface.cpp.

Definition at line 75 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 203 of file free_surface.cpp.

◆ 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 187 of file free_surface.cpp.

◆ 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 204 of file free_surface.cpp.

◆ 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 215 of file free_surface.cpp.

◆ debug

constexpr bool debug = true

constexpr

Examples: free_surface.cpp.

Definition at line 144 of file free_surface.cpp.

◆ eps

double eps = 1e-12

Examples: free_surface.cpp.

Definition at line 175 of file free_surface.cpp.

◆ eta

double eta = h

Examples: free_surface.cpp, and plot_base.cpp.

Definition at line 170 of file free_surface.cpp.

◆ eta2

double eta2 = eta * eta

Examples: free_surface.cpp.

Definition at line 171 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 149 of file free_surface.cpp.

◆ 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 252 of file free_surface.cpp.

◆ 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 368 of file free_surface.cpp.

◆ 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 340 of file free_surface.cpp.

◆ get_f

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

Examples: free_surface.cpp.

Definition at line 221 of file free_surface.cpp.

◆ get_f_dh

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

Examples: free_surface.cpp.

Definition at line 222 of file free_surface.cpp.

◆ get_fe_bit

auto get_fe_bit

Initial value:

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

Examples: free_surface.cpp.

Definition at line 314 of file free_surface.cpp.

◆ 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 318 of file free_surface.cpp.

◆ get_M

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

Examples: free_surface.cpp.

Definition at line 249 of file free_surface.cpp.

◆ get_M0

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

Examples: free_surface.cpp.

Definition at line 224 of file free_surface.cpp.

◆ get_M0_dh

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

Examples: free_surface.cpp.

Definition at line 225 of file free_surface.cpp.

◆ get_M2

auto get_M2

Initial value:

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

Examples: free_surface.cpp.

Definition at line 227 of file free_surface.cpp.

◆ get_M2_dh

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

Examples: free_surface.cpp.

Definition at line 231 of file free_surface.cpp.

◆ 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 233 of file free_surface.cpp.

◆ 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 242 of file free_surface.cpp.

◆ get_M_dh

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

Examples: free_surface.cpp.

Definition at line 250 of file free_surface.cpp.

◆ get_projection_bit

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

Examples: free_surface.cpp.

Definition at line 154 of file free_surface.cpp.

◆ get_skin_child_bit

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

Examples: free_surface.cpp.

Definition at line 153 of file free_surface.cpp.

◆ get_skin_parent_bit

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

Examples: free_surface.cpp.

Definition at line 152 of file free_surface.cpp.

◆ get_skin_projection_bit

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

Examples: free_surface.cpp.

Definition at line 155 of file free_surface.cpp.

◆ get_start_bit

auto get_start_bit

Initial value:

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

Examples: free_surface.cpp.

Definition at line 146 of file free_surface.cpp.

◆ h

double h = 0.01 / 8

Examples: free_surface.cpp.

Definition at line 169 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 82 of file free_surface.cpp.

◆ i

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

Definition at line 132 of file free_surface.cpp.

◆ init_h

auto init_h

Initial value:

= [](double r, double y, double theta) {
 
  return bubble_device(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 295 of file free_surface.cpp.

◆ integration_rule

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

Examples: 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_divergence_operator.cpp, heat_equation.cpp, helmholtz.cpp, nonlinear_elastic.cpp, photon_diffusion.cpp, scalar_check_approximation.cpp, shallow_wave.cpp, thermo_elastic.cpp, and wave_equation.cpp.

Definition at line 185 of file free_surface.cpp.

◆ j

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

Definition at line 133 of file free_surface.cpp.

◆ k

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

Definition at line 134 of file free_surface.cpp.

◆ kappa

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

Examples: free_surface.cpp.

Definition at line 183 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 270 of file free_surface.cpp.

◆ 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 260 of file free_surface.cpp.

◆ l

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

Examples: free_surface.cpp.

Definition at line 135 of file free_surface.cpp.

◆ lambda

double lambda = 73

surface tension

Examples: free_surface.cpp.

Definition at line 165 of file free_surface.cpp.

◆ marker

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

Definition at line 313 of file free_surface.cpp.

◆ md

double md = 1e-2

Examples: free_surface.cpp.

Definition at line 174 of file free_surface.cpp.

◆ mu_ave

double mu_ave = (mu_p + mu_m) / 2

Examples: free_surface.cpp.

Definition at line 180 of file free_surface.cpp.

◆ mu_diff

double mu_diff = (mu_p - mu_m) / 2

Examples: free_surface.cpp.

Definition at line 181 of file free_surface.cpp.

◆ mu_m

double mu_m = 0.010101 * 1e2

Examples: free_surface.cpp.

Definition at line 162 of file free_surface.cpp.

◆ mu_p

double mu_p = 0.000182 * 1e2

Examples: free_surface.cpp.

Definition at line 164 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 201 of file free_surface.cpp.

◆ my_min

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

Examples: free_surface.cpp.

Definition at line 202 of file free_surface.cpp.

◆ nb_levels

int nb_levels = 4

Examples: free_surface.cpp.

Definition at line 141 of file free_surface.cpp.

◆ order

int order = 3

approximation order

Examples: free_surface.cpp.

Definition at line 140 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 211 of file free_surface.cpp.

◆ refine_overlap

int refine_overlap = 4

Examples: free_surface.cpp.

Definition at line 142 of file free_surface.cpp.

◆ rho_ave

double rho_ave = (rho_p + rho_m) / 2

Examples: free_surface.cpp.

Definition at line 178 of file free_surface.cpp.

◆ rho_diff

double rho_diff = (rho_p - rho_m) / 2

Examples: free_surface.cpp.

Definition at line 179 of file free_surface.cpp.

◆ rho_m

double rho_m = 0.998

Examples: free_surface.cpp.

Definition at line 161 of file free_surface.cpp.

◆ rho_p

double rho_p = 0.0012

Examples: free_surface.cpp.

Definition at line 163 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 324 of file free_surface.cpp.

◆ SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

Examples: dynamic_first_order_con_law.cpp, and free_surface.cpp.

Definition at line 78 of file free_surface.cpp.

◆ t_kd

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

constexpr

◆ tol

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

Definition at line 176 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 467 of file free_surface.cpp.

◆ U_FIELD_DIM

constexpr int U_FIELD_DIM = SPACE_DIM

constexpr

Examples: free_surface.cpp.

Definition at line 79 of file free_surface.cpp.

◆ W

double W = 0.25

Examples: free_surface.cpp.

Definition at line 166 of file free_surface.cpp.

◆ wetting_angle

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

Examples: free_surface.cpp.

Definition at line 310 of file free_surface.cpp.

◆ 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 194 of file free_surface.cpp.

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ AssemblyBoundaryEleOp

◆ AssemblyDomainEleOp

◆ BoundaryEle

◆ BoundaryEleOp

◆ BoundaryNaturalBC

◆ BoundaryParentEle

◆ DomainEle

◆ DomainEleOp

◆ DomianParentEle

◆ ElementsAndOps

◆ EntData

◆ OpBoundaryAssembleScalar

◆ OpBoundaryMassL

◆ OpDomainAssembleScalar

◆ OpDomainAssembleVector

◆ OpDomainMassG

◆ OpDomainMassH

◆ OpDomainMassP

◆ OpDomainMassU

◆ OpFluidFlux

◆ OpMixScalarTimesDiv

◆ PostProcEleBdyCont

◆ PostProcEleDomain

◆ PostProcEleDomainCont

◆ SideEle

◆ SideOp

Enumeration Type Documentation

◆ FR

Function Documentation

◆ main()

Variable Documentation

◆ A

◆ a0

◆ BASE_DIM

◆ bit

◆ bubble_device

◆ capillary_tube

◆ coord_type

◆ cut_off

◆ cylindrical

◆ d_cut_off

◆ d_phase_function_h

◆ debug

◆ eps

◆ eta

◆ eta2

◆ get_current_bit

◆ get_D

◆ get_dofs_ents_all

◆ get_dofs_ents_by_field_name

◆ get_f

◆ get_f_dh

◆ get_fe_bit

◆ get_global_size

◆ get_M

◆ get_M0

◆ get_M0_dh

◆ get_M2

◆ get_M2_dh

◆ get_M3

◆ get_M3_dh

◆ get_M_dh

◆ get_projection_bit

◆ get_skin_child_bit

◆ get_skin_parent_bit

◆ get_skin_projection_bit

◆ get_start_bit

◆ h

◆ help

◆ I

◆ i

◆ init_h

◆ integration_rule

◆ j