disjoint_neighbor_test.C File Reference

Go to the source code of this file.

Classes
class	DisjointNeighborTest

Functions
Number	heat_exact (const Point &p, const Parameters &, const std::string &, const std::string &)
Number	left_solution_fn (const Point &, const Parameters &, const std::string &, const std::string &)
Number	right_solution_fn (const Point &, const Parameters &params, const std::string &, const std::string &)
void	assemble_temperature_jump (EquationSystems &es, const std::string &)
	CPPUNIT_TEST_SUITE_REGISTRATION (DisjointNeighborTest)

Variables
static const Real	b = 1.0
static const Real	Cgap = 1.0
static const boundary_id_type	left_id = 3
static const boundary_id_type	right_id = 1
static const boundary_id_type	interface_left_id = 5
static const boundary_id_type	interface_right_id = 6
static const boundary_id_type	interface1_left_id = 5
static const boundary_id_type	interface1_right_id = 6
static const boundary_id_type	interface2_left_id = 7
static const boundary_id_type	interface2_right_id = 8
static const boundary_id_type	interface3_left_id = 9
static const boundary_id_type	interface3_right_id = 10
static const int	Nx = 2
static const int	Ny = 2

Function Documentation

assemble_temperature_jump()

void assemble_temperature_jump	(	EquationSystems &	es,
		const std::string &
	)

Definition at line 69 of file disjoint_neighbor_test.C.

References libMesh::SparseMatrix< T >::add_matrix(), libMesh::NumericVector< T >::add_vector(), b, libMesh::FEGenericBase< OutputType >::build(), Cgap, libMesh::NumericVector< T >::close(), libMesh::SparseMatrix< T >::close(), libMesh::FEType::default_quadrature_order(), libMesh::MeshBase::get_boundary_info(), libMesh::System::get_dof_map(), libMesh::EquationSystems::get_mesh(), libMesh::EquationSystems::get_system(), interface_left_id, interface_right_id, libMesh::invalid_uint, libMesh::libmesh_assert(), libMesh::ImplicitSystem::matrix, mesh, libMesh::MeshBase::mesh_dimension(), libMesh::Elem::neighbor_ptr(), libMesh::DenseVector< T >::resize(), libMesh::DenseMatrix< T >::resize(), libMesh::ExplicitSystem::rhs, libMesh::BoundaryInfo::side_with_boundary_id(), libMesh::Elem::which_neighbor_am_i(), libMesh::DenseMatrix< T >::zero(), and libMesh::DenseVector< T >::zero().

Referenced by DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().

  {
    const MeshBase &mesh = es.get_mesh();
    LinearImplicitSystem &system = es.get_system<LinearImplicitSystem>("TempJump");
    const DofMap &dof_map = system.get_dof_map();
    FEType fe_type = dof_map.variable_type(0);

    // FE objects for volume and face integration
    auto fe = FEBase::build(mesh.mesh_dimension(), fe_type);
    auto fe_face_L = FEBase::build(mesh.mesh_dimension(), fe_type);
    auto fe_face_R = FEBase::build(mesh.mesh_dimension(), fe_type);

    // Quadrature rules
    QGauss qrule(mesh.mesh_dimension(), fe_type.default_quadrature_order());
    QGauss qface(mesh.mesh_dimension() - 1, fe_type.default_quadrature_order());
    fe->attach_quadrature_rule(&qrule);
    fe_face_L->attach_quadrature_rule(&qface);
    fe_face_R->attach_quadrature_rule(&qface);

    const auto &JxW_vol = fe->get_JxW();
    const auto &dphi_vol = fe->get_dphi();
    const auto &phi_L = fe_face_L->get_phi();
    const auto &phi_R = fe_face_R->get_phi();
    const auto &JxW_face = fe_face_L->get_JxW();

    DenseMatrix<Number> Ke;
    DenseVector<Number> Fe;
    std::vector<dof_id_type> dof_indices;

    const BoundaryInfo &boundary = mesh.get_boundary_info();
    const double conductance = Cgap * b; // interface conductance

    for (const Elem *elem : mesh.active_local_element_ptr_range())
    {
      dof_map.dof_indices(elem, dof_indices);
      const unsigned int n_dofs = dof_indices.size();

      Ke.resize(n_dofs, n_dofs);
      Fe.resize(n_dofs);
      Ke.zero();
      Fe.zero();

      // --- Volume contribution ---
      fe->reinit(elem);
      for (unsigned int qp = 0; qp < qrule.n_points(); qp++)
        for (unsigned int i = 0; i < n_dofs; i++)
          for (unsigned int j = 0; j < n_dofs; j++)
            Ke(i, j) += conductance * JxW_vol[qp] * (dphi_vol[i][qp] * dphi_vol[j][qp]);

      // --- Left-side interface ---
      unsigned int side = boundary.side_with_boundary_id(elem, interface_left_id);
      if (side != libMesh::invalid_uint)
        {
          fe_face_L->reinit(elem, side);
          const Elem *neighbor = elem->neighbor_ptr(side);
          libmesh_assert(neighbor);
          unsigned int n_side = neighbor->which_neighbor_am_i(elem);
          fe_face_R->reinit(neighbor, n_side);

          std::vector<dof_id_type> dofs_L, dofs_R;
          dof_map.dof_indices(elem, dofs_L);
          dof_map.dof_indices(neighbor, dofs_R);

          DenseMatrix<Number> Ke_ll(n_dofs, n_dofs);
          DenseMatrix<Number> Ke_lr(n_dofs, dofs_R.size());

          for (unsigned int qp = 0; qp < qface.n_points(); qp++)
            {
              for (unsigned int i = 0; i < n_dofs; i++)
                for (unsigned int j = 0; j < n_dofs; j++)
                  Ke_ll(i, j) += Cgap * phi_L[i][qp] * phi_L[j][qp] * JxW_face[qp];

              for (unsigned int i = 0; i < n_dofs; i++)
                for (unsigned int j = 0; j < dofs_R.size(); j++)
                  Ke_lr(i, j) += -Cgap * phi_L[i][qp] * phi_R[j][qp] * JxW_face[qp];
            }

          Ke += Ke_ll;
          system.matrix->add_matrix(Ke_lr, dofs_L, dofs_R);
        }

      // --- Right-side interface ---
      side = boundary.side_with_boundary_id(elem, interface_right_id);
      if (side != libMesh::invalid_uint)
        {
          fe_face_R->reinit(elem, side);
          const Elem *neighbor = elem->neighbor_ptr(side);
          libmesh_assert(neighbor);
          unsigned int n_side = neighbor->which_neighbor_am_i(elem);
          fe_face_L->reinit(neighbor, n_side);

          std::vector<dof_id_type> dofs_R, dofs_L;
          dof_map.dof_indices(elem, dofs_R);
          dof_map.dof_indices(neighbor, dofs_L);

          DenseMatrix<Number> Ke_rr(n_dofs, n_dofs);
          DenseMatrix<Number> Ke_rl(n_dofs, dofs_L.size());

          for (unsigned int qp = 0; qp < qface.n_points(); qp++)
            {
              for (unsigned int i = 0; i < n_dofs; i++)
                for (unsigned int j = 0; j < n_dofs; j++)
                  Ke_rr(i, j) += Cgap * phi_R[i][qp] * phi_R[j][qp] * JxW_face[qp];

              for (unsigned int i = 0; i < n_dofs; i++)
                for (unsigned int j = 0; j < dofs_L.size(); j++)
                  Ke_rl(i, j) += -Cgap * phi_R[i][qp] * phi_L[j][qp] * JxW_face[qp];
            }

          Ke += Ke_rr;
          system.matrix->add_matrix(Ke_rl, dofs_R, dofs_L);
        }

      // Apply constraints and add local contributions
      dof_map.heterogenously_constrain_element_matrix_and_vector(Ke, Fe, dof_indices);
      system.matrix->add_matrix(Ke, dof_indices);
      system.rhs->add_vector(Fe, dof_indices);
    }

    system.matrix->close();
    system.rhs->close();
  }

CPPUNIT_TEST_SUITE_REGISTRATION()

CPPUNIT_TEST_SUITE_REGISTRATION

(

DisjointNeighborTest

)

heat_exact()

Number heat_exact	(	const Point &	p,
		const Parameters &	,
		const std::string &	,
		const std::string &
	)

Definition at line 43 of file disjoint_neighbor_test.C.

References b.

Referenced by DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().

{
  return (p(0) < 0.5) ? p(0) : p(0) + b;
}

left_solution_fn()

Number left_solution_fn	(	const Point &	,
		const Parameters &	,
		const std::string &	,
		const std::string &
	)

Definition at line 51 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().

{
  return 0.0;
}

right_solution_fn()

Number right_solution_fn	(	const Point &	,
		const Parameters &	params,
		const std::string &	,
		const std::string &
	)

Definition at line 59 of file disjoint_neighbor_test.C.

References b, libMesh::Parameters::get(), and libMesh::Real.

Referenced by DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().

{
  const Real b = params.get<Real>("b");
  return 1.0 + b;
}

Variable Documentation

b

const Real b = 1.0

static

Definition at line 24 of file disjoint_neighbor_test.C.

Referenced by libMesh::DenseMatrix< Real >::_cholesky_back_substitute(), libMesh::DenseMatrix< Real >::_lu_back_substitute(), libMesh::DenseMatrix< Real >::_lu_back_substitute_lapack(), assemble_elasticity(), assemble_shell(), assemble_temperature_jump(), libMesh::Parallel::BinSorter< KeyType, IdxType >::binsort(), libMesh::System::boundary_project_solution(), libMesh::System::boundary_project_vector(), libMesh::Parallel::Histogram< KeyType, IdxType >::build_histogram(), NonlinearNeoHookeCurrentConfig::calculate_stress(), libMesh::ExodusII_IO_Helper::check_existing_vars(), libMesh::DenseMatrix< Real >::cholesky_solve(), libMesh::DistributedVector< T >::close(), libMesh::VariationalSmootherConstraint::constrain_node_to_line(), ContainsPointTest::containsPointTri3Helper(), libMesh::cross_norm(), libMesh::cross_norm_sq(), DMCreateDomainDecomposition_libMesh(), DMlibMeshCreateDomainDecompositionDM(), libMesh::QGauss::dunavant_rule2(), libMesh::VariationalSmootherSystem::element_time_derivative(), libMesh::ThreadBufferedSyncbuf::emit_from_thread_local_buffer(), libMesh::DofMap::get_variable_array(), libMesh::Edge3::has_invertible_map(), libMesh::Edge4::has_invertible_map(), libMesh::Utility::hashword(), libMesh::Utility::hashword2(), heat_exact(), libMesh::QGauss::init_2D(), libMesh::QGauss::init_3D(), libMesh::TensorTools::inner_product(), libMesh::Poly2TriTriangulator::insert_refinement_points(), libMesh::TypeNTensor< N, T >::int_pow(), libMesh::intersect_constraints(), libMesh::TypeTensor< T >::inverse(), libMesh::QGauss::keast_rule(), libMesh::DenseMatrix< Real >::lu_solve(), libMesh::Parallel::Histogram< KeyType, IdxType >::make_histogram(), libMesh::PetscVector< libMesh::Number >::map_global_to_local_index(), libMesh::TriangulatorInterface::MeshedHole::MeshedHole(), libMesh::PeriodicBoundaries::neighbor(), libMesh::CompareElemIdsByLevel::operator()(), std::plus< boost::multiprecision::float128 >::operator()(), std::multiplies< boost::multiprecision::float128 >::operator()(), libMesh::BoundaryProjectSolution::operator()(), libMesh::Utility::ReverseBytes::operator()(), libMesh::Utility::CompareUnderlying::operator()(), libMesh::CompareDofObjectsByID::operator()(), libMesh::CompareDofObjectsByPIDAndThenID::operator()(), libMesh::operator*(), libMesh::DenseMatrix< Real >::outer_product(), libMesh::outer_product(), libMesh::RadialBasisInterpolation< KDDim, RBF >::prepare_for_use(), libMesh::Face::quasicircumcenter(), right_solution_fn(), libMesh::VTKIO::set_compression(), set_system_parameters(), libMesh::SIGN(), libMesh::TypeTensor< T >::solve(), libMesh::UnstructuredMesh::stitching_helper(), AllTriTest::test_helper_c0polyhedron(), ContainsPointTest::testContainsPointTri3(), InfFERadialTest::testInfQuants(), DenseMatrixTest::testOuterProduct(), TypeTensorTest::testOuterProduct(), DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch(), DisjointNeighborTest::testStitchCrossMesh(), DisjointNeighborTest::testTempJump(), DisjointNeighborTest::testTempJumpRefine(), MeshSmootherTest::testVariationalSmoother(), triangulate_domain(), libMesh::triple_product(), libMesh::Edge3::volume(), and libMesh::Tet4::volume().

Cgap

const Real Cgap = 1.0

static

Definition at line 25 of file disjoint_neighbor_test.C.

Referenced by assemble_temperature_jump().

interface1_left_id

const boundary_id_type interface1_left_id = 5

static

Definition at line 33 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface1_right_id

const boundary_id_type interface1_right_id = 6

static

Definition at line 34 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface2_left_id

const boundary_id_type interface2_left_id = 7

static

Definition at line 35 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface2_right_id

const boundary_id_type interface2_right_id = 8

static

Definition at line 36 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface3_left_id

const boundary_id_type interface3_left_id = 9

static

Definition at line 37 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface3_right_id

const boundary_id_type interface3_right_id = 10

static

Definition at line 38 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), and DisjointNeighborTest::testPreserveDisjointNeighborPairsAfterStitch().

interface_left_id

const boundary_id_type interface_left_id = 5

static

Definition at line 30 of file disjoint_neighbor_test.C.

Referenced by assemble_temperature_jump(), DisjointNeighborTest::build_split_mesh_with_interface(), DisjointNeighborTest::build_two_disjoint_elems(), DisjointNeighborTest::testDisjointNeighborConflictError(), DisjointNeighborTest::testStitchCrossMesh(), and DisjointNeighborTest::testStitchingDiscontinuousBoundaries().

interface_right_id

const boundary_id_type interface_right_id = 6

static

Definition at line 31 of file disjoint_neighbor_test.C.

Referenced by assemble_temperature_jump(), DisjointNeighborTest::build_split_mesh_with_interface(), DisjointNeighborTest::build_two_disjoint_elems(), DisjointNeighborTest::testDisjointNeighborConflictError(), DisjointNeighborTest::testStitchCrossMesh(), and DisjointNeighborTest::testStitchingDiscontinuousBoundaries().

left_id

const boundary_id_type left_id = 3

static

Definition at line 28 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), DisjointNeighborTest::build_split_mesh_with_interface(), DisjointNeighborTest::build_two_disjoint_elems(), BoundaryMeshTest::sync_and_test_meshes(), DisjointNeighborTest::testDisjointNeighborConflictError(), DisjointNeighborTest::testStitchCrossMesh(), DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().

Nx

const int Nx = 2

static

Definition at line 40 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::testTempJumpLocalRefineFail(), and DisjointNeighborTest::testTempJumpRefine().

Ny

const int Ny = 2

static

Definition at line 40 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::testTempJumpLocalRefineFail(), and DisjointNeighborTest::testTempJumpRefine().

right_id

const boundary_id_type right_id = 1

static

Definition at line 29 of file disjoint_neighbor_test.C.

Referenced by DisjointNeighborTest::build_four_disjoint_elems(), DisjointNeighborTest::build_split_mesh_with_interface(), DisjointNeighborTest::build_two_disjoint_elems(), BoundaryMesh0DTest::sync_and_test_meshes(), BoundaryMeshTest::sync_and_test_meshes(), DisjointNeighborTest::testDisjointNeighborConflictError(), DisjointNeighborTest::testStitchCrossMesh(), DisjointNeighborTest::testTempJump(), and DisjointNeighborTest::testTempJumpRefine().