InterfaceKernels System

Interface kernels are meant to assist in coupling different physics across sub-domains. The most straightforward example is the case in which one wants to set the flux of a specie A in subdomain 0 equal to the flux of a specie B in subdomain 1 at the boundary between subdomains 0 and 1. In mathematical terms, we might be interested in establishing the condition:

$var element = document.getElementById("moose-equation-dbcc0a24-b2da-4226-bf97-2f91810ccee4");katex.render("-D_0 \\frac{\\partial c_0}{\\partial x} = -D_1 \\frac{\\partial c_1}{\\partial x}", element, {displayMode:true,throwOnError:false});$

where $var element = document.getElementById("moose-equation-70005f8a-4a71-4d6a-9377-bb795e2dc15f");katex.render("D_i", element, {displayMode:false,throwOnError:false});$ is the diffusion coefficient of specie $var element = document.getElementById("moose-equation-886d2eb5-40c8-4e99-abb6-3ebb5c2cc972");katex.render("i", element, {displayMode:false,throwOnError:false});$ in subdomain $var element = document.getElementById("moose-equation-e535569b-d60b-4221-b764-4e7c58f5b70d");katex.render("i", element, {displayMode:false,throwOnError:false});$ , and $var element = document.getElementById("moose-equation-212853aa-9bae-485f-a023-66c0f99cee43");katex.render("c_i", element, {displayMode:false,throwOnError:false});$ is the concentration of specie $var element = document.getElementById("moose-equation-f4700fd6-6aaf-4e88-b4e5-c0485121836f");katex.render("i", element, {displayMode:false,throwOnError:false});$ in subdomain $var element = document.getElementById("moose-equation-9925a41b-894b-4e87-a648-1b9cc837064f");katex.render("i", element, {displayMode:false,throwOnError:false});$ . An example of this condition is shown in the MOOSE test directory; see files below:

(moose/test/tests/interfacekernels/2d_interface/coupled_value_coupled_flux.i)

(moose/framework/src/interfacekernels/InterfaceDiffusion.C)

(moose/framework/include/interfacekernels/InterfaceDiffusion.h)

Interface kernels can be used to provide any general flux condition at an interface, and even more generally can be used to impose any interfacial condition that requires access to values of different variables and gradients of different variables on either side of an interface. In an input file, the user will specify at a minimum the following parameters:

type: The type of interface kernel to be used
variable: This is the "primary" variable. Note that the primary variable must exist on the same subdomain as the sideset specified in the boundary parameter. The existence of a "primary" and "secondary" or "neighbor" variable ensures that the interface kernel residual and jacobian functions get called the correct number of times. variable could be $var element = document.getElementById("moose-equation-82224f02-eb8b-4f65-94f9-5fb482caa32b");katex.render("c_0", element, {displayMode:false,throwOnError:false});$ from our example above.
neighbor_var: The "secondary" variable. This could be $var element = document.getElementById("moose-equation-02af9932-48ba-43af-8508-d1687972fc81");katex.render("c_1", element, {displayMode:false,throwOnError:false});$ from our example above.
boundary: The interfacial boundary between the subdomains. Note that this must be a sideset and again must exist on the same subdomain as the primary variable. The fact that this boundary is a sideset allows access to variable gradients.

For additional information about the interface kernel system, don't hesitate to contact the MOOSE Discussion forum.

Available Objects

Moose App
ADPenaltyInterfaceDiffusionA penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.
ADVectorPenaltyInterfaceDiffusionA penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.
InterfaceDiffusionThe kernel is utilized to establish flux equivalence on an interface for variables.
InterfaceReactionImplements a reaction to establish ReactionRate=k_f*u-k_b*v at interface.
PenaltyInterfaceDiffusionA penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.
VectorPenaltyInterfaceDiffusionA penalty-based interface condition that forcesthe continuity of variables and the flux equivalence across an interface.
Tensor Mechanics App
ADCZMInterfaceKernelSmallStrainCZM Interface kernel to use when using the small strain kinematic formulation.
ADCZMInterfaceKernelTotalLagrangianCZM Interface kernel to use when using the total Lagrangian formulation.
CZMInterfaceKernelSmallStrainCZM Interface kernel to use when using the Small Strain kinematic formulation.
CZMInterfaceKernelTotalLagrangian

Available Actions

Moose App
AddInterfaceKernelActionAdd an InterfaceKernel object to the simulation.

(moose/test/tests/interfacekernels/2d_interface/coupled_value_coupled_flux.i)

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 2
    xmax = 2
    ny = 2
    ymax = 2
  []
  [./subdomain1]
    input = gen
    type = SubdomainBoundingBoxGenerator
    bottom_left = '0 0 0'
    top_right = '1 1 0'
    block_id = 1
  [../]
  [./interface]
    type = SideSetsBetweenSubdomainsGenerator
    input = subdomain1
    primary_block = '0'
    paired_block = '1'
    new_boundary = 'primary0_interface'
  [../]
  [./break_boundary]
    input = interface
    type = BreakBoundaryOnSubdomainGenerator
  [../]
[]

[Variables]
  [./u]
    order = FIRST
    family = LAGRANGE
    block = 0
  [../]

  [./v]
    order = FIRST
    family = LAGRANGE
    block = 1
  [../]
[]

[Kernels]
  [./diff_u]
    type = CoeffParamDiffusion
    variable = u
    D = 4
    block = 0
  [../]
  [./diff_v]
    type = CoeffParamDiffusion
    variable = v
    D = 2
    block = 1
  [../]
  [./source_u]
    type = BodyForce
    variable = u
    value = 1
  [../]
[]

[InterfaceKernels]
  [./interface]
    type = PenaltyInterfaceDiffusion
    variable = u
    neighbor_var = v
    boundary = primary0_interface
    penalty = 1e6
  [../]
[]

[BCs]
  [./u]
    type = VacuumBC
    variable = u
    boundary = 'left_to_0 bottom_to_0 right top'
  [../]
  [./v]
    type = VacuumBC
    variable = v
    boundary = 'left_to_1 bottom_to_1'
  [../]
[]

[Postprocessors]
  [./u_int]
    type = ElementIntegralVariablePostprocessor
    variable = u
    block = 0
  [../]
  [./v_int]
    type = ElementIntegralVariablePostprocessor
    variable = v
    block = 1
  [../]
[]

[Preconditioning]
  [./smp]
    type = SMP
    full = true
  [../]
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
[]

[Outputs]
  exodus = true
  print_linear_residuals = true
[]

(moose/framework/src/interfacekernels/InterfaceDiffusion.C)

// This file is part of the MOOSE framework
// https://www.mooseframework.org
//
// All rights reserved, see COPYRIGHT for full restrictions
// https://github.com/idaholab/moose/blob/master/COPYRIGHT
//
// Licensed under LGPL 2.1, please see LICENSE for details
// https://www.gnu.org/licenses/lgpl-2.1.html

#include "InterfaceDiffusion.h"

registerMooseObject("MooseApp", InterfaceDiffusion);

InputParameters
InterfaceDiffusion::validParams()
{
  InputParameters params = InterfaceKernel::validParams();
  params.addParam<MaterialPropertyName>("D", "D", "The diffusion coefficient.");
  params.addParam<MaterialPropertyName>(
      "D_neighbor", "D_neighbor", "The neighboring diffusion coefficient.");
  params.addClassDescription(
      "The kernel is utilized to establish flux equivalence on an interface for variables.");
  return params;
}

InterfaceDiffusion::InterfaceDiffusion(const InputParameters & parameters)
  : InterfaceKernel(parameters),
    _D(getMaterialProperty<Real>("D")),
    _D_neighbor(getNeighborMaterialProperty<Real>("D_neighbor"))
{
}

Real
InterfaceDiffusion::computeQpResidual(Moose::DGResidualType type)
{
  Real r = 0;

  switch (type)
  {
    case Moose::Element:
      r = _test[_i][_qp] * -_D_neighbor[_qp] * _grad_neighbor_value[_qp] * _normals[_qp];
      break;

    case Moose::Neighbor:
      r = _test_neighbor[_i][_qp] * _D[_qp] * _grad_u[_qp] * _normals[_qp];
      break;
  }

  return r;
}

Real
InterfaceDiffusion::computeQpJacobian(Moose::DGJacobianType type)
{
  Real jac = 0;

  switch (type)
  {
    case Moose::ElementElement:
    case Moose::NeighborNeighbor:
      break;

    case Moose::NeighborElement:
      jac = _test_neighbor[_i][_qp] * _D[_qp] * _grad_phi[_j][_qp] * _normals[_qp];
      break;

    case Moose::ElementNeighbor:
      jac = _test[_i][_qp] * -_D_neighbor[_qp] * _grad_phi_neighbor[_j][_qp] * _normals[_qp];
      break;
  }

  return jac;
}

(moose/framework/include/interfacekernels/InterfaceDiffusion.h)

// This file is part of the MOOSE framework
// https://www.mooseframework.org
//
// All rights reserved, see COPYRIGHT for full restrictions
// https://github.com/idaholab/moose/blob/master/COPYRIGHT
//
// Licensed under LGPL 2.1, please see LICENSE for details
// https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "InterfaceKernel.h"

/**
 * DG kernel for interfacing diffusion between two variables on adjacent blocks
 */
class InterfaceDiffusion : public InterfaceKernel
{
public:
  static InputParameters validParams();

  InterfaceDiffusion(const InputParameters & parameters);

protected:
  virtual Real computeQpResidual(Moose::DGResidualType type) override;
  virtual Real computeQpJacobian(Moose::DGJacobianType type) override;

  const MaterialProperty<Real> & _D;
  const MaterialProperty<Real> & _D_neighbor;
};