SubdomainInterceptedGenerator

Classifies mesh elements into inside, outside, and optionally intercepted subdomains based on an analytical signed distance function expression.

Overview

The SubdomainInterceptedGenerator classifies and assigns subdomain IDs to elements in an input mesh based on an implicit interface defined by an analytical signed distance function expression $var element = document.getElementById("moose-equation-4f98c3d8-a189-4e58-bb2a-5efad9acbb81");katex.render("\\phi(x, y, z)", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ .

This generator evaluates the signed distance function across the mesh to assign elements to: - An inside** subdomain ("subdomain_id_inside"), - An outside** subdomain ("subdomain_id_outside"), and - Optionally, a dedicated intercepted** subdomain ("subdomain_id_intercepted") for elements cut by the interface.

Classification Method

The interface is defined by the level-set condition: $var element = document.getElementById("moose-equation-c2f1f660-135b-480e-910d-47c316611c0f");katex.render("\\phi(x, y, z) = \\text{threshold}", element, {displayMode:true,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ where $var element = document.getElementById("moose-equation-a95fc649-aef0-47af-8343-b3f1c16830b5");katex.render("\\text{threshold}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ is set by "threshold" (default is $var element = document.getElementById("moose-equation-a172fd17-8d91-4160-8e49-b42f02c5a57f");katex.render("0.0", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ ).

Domain Orientation

The orientation of the inside domain relative to the signed distance function is controlled by "is_domain_inside_surface": - is_domain_inside_surface = true: Follows the standard convention where negative signed distance values represent the interior of an enclosed region. Points satisfying $var element = document.getElementById("moose-equation-c9961707-bf81-498a-9d6e-b652516da423");katex.render("\\phi(x, y, z) < \\text{threshold}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ are categorized as inside the domain, and points satisfying $var element = document.getElementById("moose-equation-04bfc58b-2a83-46dc-84b5-0bfd1d3b598d");katex.render("\\phi(x, y, z) > \\text{threshold}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ are categorized as outside. - is_domain_inside_surface = false: Inverts the convention. Points satisfying $var element = document.getElementById("moose-equation-42966b68-9ba1-43b9-8452-f2012c21f849");katex.render("\\phi(x, y, z) > \\text{threshold}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ are categorized as inside the domain, and points satisfying $var element = document.getElementById("moose-equation-190e1140-de00-4e80-97f4-67adfd763150");katex.render("\\phi(x, y, z) < \\text{threshold}", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ are categorized as outside.

Element Occupancy and Intercepted Elements

Each element in the input mesh is evaluated at its vertices (nodes) and sampled across its volume using a Gauss quadrature rule of order "qrule_order" (default is 9). The generator measures: 1. Whether all element nodes lie inside the domain (all_nodes_in_domain). 2. Whether all element nodes lie outside the domain (all_nodes_outside_domain). 3. The fraction of the element volume occupied by the inside region (fraction_of_element_in_domain), evaluated via numerical quadrature.

Elements are then categorized as follows: - Fully inside**: When all nodes lie inside and fraction_of_element_in_domain == 1.0, the element is assigned "subdomain_id_inside". - Fully outside**: When all nodes lie outside and fraction_of_element_in_domain == 0.0, the element is assigned "subdomain_id_outside". - Intercepted (cut) elements**: Elements that contain both inside and outside points are cut by the implicit boundary: - If "mark_intercepted" is true, these elements are assigned a dedicated subdomain ID specified by "subdomain_id_intercepted". - If "mark_intercepted" is false (the default), the element is resolved to either inside or outside based on the volume fraction threshold "lambda" (default is $var element = document.getElementById("moose-equation-be83c90f-e592-4d5c-9cc9-8a87b028717f");katex.render("0.5", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ ). If the fraction of the element volume outside the domain exceeds $var element = document.getElementById("moose-equation-58bf89eb-e62f-48e3-91e8-26e371051e03");katex.render("\\lambda", element, {displayMode:false,throwOnError:false,macros:{"\\eqc":"\\,,","\\eqp":"\\,.","\\pd":"\\frac{\\partial #1}{\\partial #2}","\\pr":"\\left(#1\\right)","\\ddt":"\\frac{d #1}{d t}"}});$ , the element is assigned to outside; otherwise, it is assigned to inside.

Multi-Geometry Workflows

When constructing domains with multiple geometric features, inclusions, or composite shapes, multiple SubdomainInterceptedGenerator instances can be chained in sequence. To prevent subsequent generators from unintentionally overwriting previously assigned subdomains, several parameters are available:

"multi_geo": Enables multi-geometry retention checks.
"keep_inside_as_inside": When multi_geo = true, any element whose existing subdomain ID already matches "subdomain_id_inside" is skipped and preserved, preventing later generators from overwriting it.
"keep_outside_as_outside": When multi_geo = true, any element whose existing subdomain ID already matches "subdomain_id_outside" is skipped and preserved.
"modify_inside_only": When true, only elements categorized as inside receive "subdomain_id_inside". Elements categorized as outside retain their pre-existing subdomain IDs.
"modify_outside_only": When true, only elements categorized as outside receive "subdomain_id_outside". Elements categorized as inside retain their pre-existing subdomain IDs.

Example Input File Syntax

Inclined Interface Partitioning

The following example demonstrates using SubdomainInterceptedGenerator to partition a 2D mesh across an inclined plane into subdomains 1 and 2, with is_domain_inside_surface = true:

[Mesh<<<{"href": "../../syntax/Mesh/index.html"}>>>]
  [gen]
    type = CartesianMeshGenerator<<<{"description": "This CartesianMeshGenerator creates a non-uniform Cartesian mesh.", "href": "CartesianMeshGenerator.html"}>>>
    dim<<<{"description": "The dimension of the mesh to be generated"}>>> = 2
    dx<<<{"description": "Intervals in the X direction"}>>> = '${x0} ${x0}'
    dy<<<{"description": "Intervals in the Y direction (required when dim>1 otherwise ignored)"}>>> = 1
    ix<<<{"description": "Number of grids in all intervals in the X direction (default to all one)"}>>> = '${nx_half} ${nx_half}'
    iy<<<{"description": "Number of grids in all intervals in the Y direction (default to all one)"}>>> = '${nx}'
    subdomain_id<<<{"description": "Block IDs (default to all zero)"}>>> = '1 2'
  []
  [intercepted]
    type = SubdomainInterceptedGenerator<<<{"description": "Classifies mesh elements into inside, outside, and optionally intercepted subdomains based on an analytical signed distance function expression.", "href": "SubdomainInterceptedGenerator.html"}>>>
    input<<<{"description": "The mesh to modify"}>>> = gen
    subdomain_id_inside<<<{"description": "Subdomain ID to assign to elements categorized inside the domain."}>>> = 2
    subdomain_id_outside<<<{"description": "Subdomain ID to assign to elements categorized outside the domain."}>>> = 1
    lambda<<<{"description": "Threshold applied to the fraction outside the retained domain for a partial element. This parameter cannot be set when 'mark_intercepted' is true."}>>> = 0.5
    is_domain_inside_surface<<<{"description": "Whether the inside domain corresponds to the negative signed distance region (< threshold). If true, points with signed distance < threshold are categorized as inside and > threshold as outside. If false, points with signed distance > threshold are categorized as inside."}>>> = true
    signed_dist_function<<<{"description": "Parsed function expression (in terms of x, y, and z) defining the signed distance function to the interface."}>>> = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
  [break]
    type = BreakMeshByBlockGenerator<<<{"description": "Break the mesh at interfaces between blocks. New nodes will be generated so elements on each side of the break are no longer connected.", "href": "BreakMeshByBlockGenerator.html"}>>>
    input<<<{"description": "The mesh we want to modify"}>>> = intercepted
    split_interface<<<{"description": "Interfaces are always created on both sides and split by block pairs."}>>> = true
    add_interface_on_two_sides<<<{"description": "Interfaces are always created on both sides and split by block pairs."}>>> = true
  []
[]

Inverted Domain Orientation

When is_domain_inside_surface = false, points on the positive side of the signed distance function are categorized into the inside subdomain:

  [subdomain_intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 1
    subdomain_id_outside = 2
    lambda = 0.5
    is_domain_inside_surface = false
    signed_dist_function = '-x+${x0}'
  []

Input Parameters

inputThe mesh to modify
C++ Type:MeshGeneratorName
Controllable:No
Description:The mesh to modify
is_domain_inside_surfaceFalseWhether the inside domain corresponds to the negative signed distance region (< threshold). If true, points with signed distance < threshold are categorized as inside and > threshold as outside. If false, points with signed distance > threshold are categorized as inside.
Default:False
C++ Type:bool
Controllable:No
Description:Whether the inside domain corresponds to the negative signed distance region (< threshold). If true, points with signed distance < threshold are categorized as inside and > threshold as outside. If false, points with signed distance > threshold are categorized as inside.
signed_dist_functionParsed function expression (in terms of x, y, and z) defining the signed distance function to the interface.
C++ Type:std::string
Controllable:No
Description:Parsed function expression (in terms of x, y, and z) defining the signed distance function to the interface.
subdomain_id_insideSubdomain ID to assign to elements categorized inside the domain.
C++ Type:unsigned short
Controllable:No
Description:Subdomain ID to assign to elements categorized inside the domain.
subdomain_id_outsideSubdomain ID to assign to elements categorized outside the domain.
C++ Type:unsigned short
Controllable:No
Description:Subdomain ID to assign to elements categorized outside the domain.

Required Parameters

epsilon0Fuzzy comparison tolerance
Default:0
C++ Type:Real
Unit:(no unit assumed)
Controllable:No
Description:Fuzzy comparison tolerance
keep_inside_as_insideFalseWhen 'multi_geo' is true, elements whose existing subdomain ID matches 'subdomain_id_inside' are skipped, preserving their existing subdomain assignment.
Default:False
C++ Type:bool
Controllable:No
Description:When 'multi_geo' is true, elements whose existing subdomain ID matches 'subdomain_id_inside' are skipped, preserving their existing subdomain assignment.
keep_outside_as_outsideFalseWhen 'multi_geo' is true, elements whose existing subdomain ID matches 'subdomain_id_outside' are skipped, preserving their existing subdomain assignment.
Default:False
C++ Type:bool
Controllable:No
Description:When 'multi_geo' is true, elements whose existing subdomain ID matches 'subdomain_id_outside' are skipped, preserving their existing subdomain assignment.
lambda0.5Threshold applied to the fraction outside the retained domain for a partial element. This parameter cannot be set when 'mark_intercepted' is true.
Default:0.5
C++ Type:Real
Unit:(no unit assumed)
Range:lambda >= 0 & lambda <= 1
Controllable:No
Description:Threshold applied to the fraction outside the retained domain for a partial element. This parameter cannot be set when 'mark_intercepted' is true.
mark_interceptedFalseWhether to assign a dedicated subdomain ID when intercepted.
Default:False
C++ Type:bool
Controllable:No
Description:Whether to assign a dedicated subdomain ID when intercepted.
modify_inside_onlyFalseIf true, only elements categorized as inside are updated to 'subdomain_id_inside'; elements categorized as outside retain their existing subdomain IDs.
Default:False
C++ Type:bool
Controllable:No
Description:If true, only elements categorized as inside are updated to 'subdomain_id_inside'; elements categorized as outside retain their existing subdomain IDs.
modify_outside_onlyFalseIf true, only elements categorized as outside are updated to 'subdomain_id_outside'; elements categorized as inside retain their existing subdomain IDs.
Default:False
C++ Type:bool
Controllable:No
Description:If true, only elements categorized as outside are updated to 'subdomain_id_outside'; elements categorized as inside retain their existing subdomain IDs.
multi_geoFalseWhether to enable multi-geometry handling when chaining generators, allowing elements with specific existing subdomain IDs to be preserved without reclassification (see 'keep_inside_as_inside' and 'keep_outside_as_outside').
Default:False
C++ Type:bool
Controllable:No
Description:Whether to enable multi-geometry handling when chaining generators, allowing elements with specific existing subdomain IDs to be preserved without reclassification (see 'keep_inside_as_inside' and 'keep_outside_as_outside').
qrule_order9Quadrature order used to estimate the retained-domain fraction.
Default:9
C++ Type:int
Range:qrule_order >= 0 & qrule_order <= 10
Controllable:No
Description:Quadrature order used to estimate the retained-domain fraction.
subdomain_id_intercepted65535Subdomain ID assigned to intercepted elements.
Default:65535
C++ Type:unsigned short
Controllable:No
Description:Subdomain ID assigned to intercepted elements.
threshold0Level-set threshold value defining the interface boundary (default 0.0).
Default:0
C++ Type:Real
Unit:(no unit assumed)
Controllable:No
Description:Level-set threshold value defining the interface boundary (default 0.0).

Optional Parameters

disable_fpoptimizerFalseDisable the function parser algebraic optimizer
Default:False
C++ Type:bool
Controllable:No
Description:Disable the function parser algebraic optimizer
enable_ad_cacheTrueEnable caching of function derivatives for faster startup time
Default:True
C++ Type:bool
Controllable:No
Description:Enable caching of function derivatives for faster startup time
enable_auto_optimizeTrueEnable automatic immediate optimization of derivatives
Default:True
C++ Type:bool
Controllable:No
Description:Enable automatic immediate optimization of derivatives
enable_jitTrueEnable just-in-time compilation of function expressions for faster evaluation
Default:True
C++ Type:bool
Controllable:No
Description:Enable just-in-time compilation of function expressions for faster evaluation
evalerror_behaviornanWhat to do if evaluation error occurs. Options are to pass a nan, pass a nan with a warning, throw a error, or throw an exception
Default:nan
C++ Type:MooseEnum
Options:nan, nan_warning, error, exception
Controllable:No
Description:What to do if evaluation error occurs. Options are to pass a nan, pass a nan with a warning, throw a error, or throw an exception

Parsed Expression Advanced Parameters

enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.
save_with_nameKeep the mesh from this mesh generator in memory with the name specified
C++ Type:std::string
Controllable:No
Description:Keep the mesh from this mesh generator in memory with the name specified

Advanced Parameters

nemesisFalseWhether or not to output the mesh file in the nemesisformat (only if output = true)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file in the nemesisformat (only if output = true)
outputFalseWhether or not to output the mesh file after generating the mesh
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to output the mesh file after generating the mesh
show_infoFalseWhether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not to show mesh info after generating the mesh (bounding box, element types, sidesets, nodesets, subdomains, etc)

Debugging Parameters

Input Files

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_ad_elastic.i)
(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_common.i)
(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_sczm_tri_material.i)
(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_czm_only_quad_material_mimic.i)
(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_nonad_salehani.i)
(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/sczm_triangle_mms.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_common.i)

nx = 8
nx_half = '${fparse nx / 2}'
E = 1e3
poi = 0.3
theta = 60
theta_rad = '${fparse theta * pi / 180}'
x0 = 0.5

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [gen]
    type = CartesianMeshGenerator
    dim = 2
    dx = '${x0} ${x0}'
    dy = 1
    ix = '${nx_half} ${nx_half}'
    iy = '${nx}'
    subdomain_id = '1 2'
  []
  [intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 2
    subdomain_id_outside = 1
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
  [break]
    type = BreakMeshByBlockGenerator
    input = intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = signed_dist_func
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = INITIAL
  []
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
[]

[AuxVariables]
  [react_x]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = ref
    v = disp_x
    variable = react_x
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [pull]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [react_x]
    type = NodalSum
    variable = react_x
    boundary = right
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = none
  nl_max_its = 50
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-10
  dt = 1
  end_time = 5
[]

[Outputs]
  [csv]
    type = CSV
    precision = 15
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/sczm_triangle_mms.i)

E1 = 0.1
E2 = 1
poi = 0
x0 = 0.25
a = '${fparse (pi * cos(pi * x0) * (1 - E1 / E2)) / (2 * x0)}'
b = '${fparse -pi * cos(pi * x0) * (E1 + (x0 / 2) * (1 - E1 / E2))}'

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Mesh]
  [gen]
    type = FileMeshGenerator
    file = sczm_triangle_mesh.xda
  []

  [subdomain_intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 1
    subdomain_id_outside = 2
    lambda = 0.5
    is_domain_inside_surface = false
    signed_dist_function = '-x+${x0}'
  []

  [break]
    type = BreakMeshByBlockGenerator
    block_pairs = '1 2'
    input = subdomain_intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []

  [ss]
    type = SideSetsFromNormalsGenerator
    input = break
    normals = '-1 0 0
                1 0 0'
    fixed_normal = true
    new_boundary = 'left right'
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        use_automatic_differentiation = true
        generate_output = 'stress_xx stress_xy stress_yy stress_zz strain_xx strain_xy strain_yy strain_zz'
      []
    []
  []
[]

[Physics/SolidMechanics/ShiftedCohesiveZone]
  [czm_ik]
    use_automatic_differentiation = true
    boundary = 'Block1_Block2'
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = 'signed_dist_func'
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = 'INITIAL'
  []
[]

[Materials]
  [elastic_stress]
    type = ADComputeLinearElasticStress
  []
  [elasticity_tensor1]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E1}
    block = 1
  []
  [elasticity_tensor2]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E2}
    block = 2
  []
  [interface_traction]
    type = ADPureElasticTractionSeparation
    normal_stiffness = 1
    tangent_stiffness = 1
    boundary = 'Block1_Block2'
  []
[]

[AuxVariables]
  [a]
    family = MONOMIAL
  []
[]

[AuxKernels]
  [a]
    type = ConstantAux
    variable = a
    value = ${a}
  []
[]

[Kernels]
  [source_x_left]
    type = BodyForce
    variable = disp_x
    function = bx_func_left
    block = 1
  []
  [source_x_right]
    type = BodyForce
    variable = disp_x
    function = bx_func_right
    block = 2
  []
[]

[Functions]
  [signed_dist_func]
    type = ParsedFunction
    expression = '-x+${x0}'
  []
  [ux_func_left]
    type = ParsedFunction
    expression = '-sin(pi*x)'
  []

  [ux_func_right]
    type = ParsedFunction
    expression = '-sin(pi*x)+${a}*x^2+${b}'
  []

  [bx_func_left]
    type = ParsedFunction
    expression = '-pi^2*${E1}*sin(pi*x)'
  []

  [bx_func_right]
    type = ParsedFunction
    expression = '-${E2}*(2*${a} + pi^2*sin(pi*x))'
  []
[]

[BCs]
  [y_anchor]
    type = DirichletBC
    variable = disp_y
    boundary = 'left right'
    value = 0.0
  []

  [x_left]
    type = FunctionDirichletBC
    boundary = 'left'
    function = ux_func_left
    variable = disp_x
  []

  [x_right]
    type = FunctionDirichletBC
    boundary = 'right'
    function = ux_func_right
    variable = disp_x
  []
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = 'none'
[]

[Postprocessors]
  [l2_error_left_x]
    type = ElementL2Error
    variable = disp_x
    function = ux_func_left
    block = '1'
  []

  [l2_error_right_x]
    type = ElementL2Error
    variable = disp_x
    function = ux_func_right
    block = '2'
  []

  [l2_sum]
    type = SumPostprocessor
    values = 'l2_error_left_x l2_error_right_x'
  []

  [n_elements]
    type = NumElements
  []
[]

[Outputs]
  exodus = true
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_ad_elastic.i)

nx = 8
nx_half = '${fparse nx / 2}'
E = 1e3
poi = 0.3
theta = 60
theta_rad = '${fparse theta * pi / 180}'
x0 = 0.5

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [gen]
    type = CartesianMeshGenerator
    dim = 2
    dx = '${x0} ${x0}'
    dy = 1
    ix = '${nx_half} ${nx_half}'
    iy = '${nx}'
    subdomain_id = '1 2'
  []
  [intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 2
    subdomain_id_outside = 1
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
  [break]
    type = BreakMeshByBlockGenerator
    input = intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = signed_dist_func
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = INITIAL
  []
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
[]

[AuxVariables]
  [react_x]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = ref
    v = disp_x
    variable = react_x
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [pull]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [react_x]
    type = NodalSum
    variable = react_x
    boundary = right
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = none
  nl_max_its = 50
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-10
  dt = 1
  end_time = 5
[]

[Outputs]
  [csv]
    type = CSV
    precision = 15
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        extra_vector_tags = ref
        add_variables = true
        use_automatic_differentiation = true
      []
    []
  []
[]

[Physics/SolidMechanics/ShiftedCohesiveZone]
  [czm]
    use_automatic_differentiation = true
    boundary = Block1_Block2
  []
[]

[Materials]
  [elastic_stress]
    type = ADComputeLinearElasticStress
  []
  [elasticity_tensor_1]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E}
    block = 1
  []
  [elasticity_tensor_2]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E}
    block = 2
  []
  [czm]
    type = ADPureElasticTractionSeparation
    normal_stiffness = 5000
    tangent_stiffness = 5000
    boundary = Block1_Block2
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_common.i)

nx = 8
nx_half = '${fparse nx / 2}'
E = 1e3
poi = 0.3
theta = 60
theta_rad = '${fparse theta * pi / 180}'
x0 = 0.5

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [gen]
    type = CartesianMeshGenerator
    dim = 2
    dx = '${x0} ${x0}'
    dy = 1
    ix = '${nx_half} ${nx_half}'
    iy = '${nx}'
    subdomain_id = '1 2'
  []
  [intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 2
    subdomain_id_outside = 1
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
  [break]
    type = BreakMeshByBlockGenerator
    input = intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = signed_dist_func
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = INITIAL
  []
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
[]

[AuxVariables]
  [react_x]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = ref
    v = disp_x
    variable = react_x
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [pull]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [react_x]
    type = NodalSum
    variable = react_x
    boundary = right
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = none
  nl_max_its = 50
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-10
  dt = 1
  end_time = 5
[]

[Outputs]
  [csv]
    type = CSV
    precision = 15
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_sczm_tri_material.i)

E1 = 1e3
E2 = 1e3
poi = 0.3

desired_area = 0.1
number_of_point = '${fparse int(sqrt(1/(2*desired_area))) -1}'

x0 = 0.5

dt = 1

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [square_boundary]
    type = PolyLineMeshGenerator
    points = '0.0 0.0 0.0
            1.0 0.0 0.0
            1.0 1.0 0.0
            0.0 1.0 0.0'
    loop = true
  []

  [gen]
    type = XYDelaunayGenerator
    boundary = 'square_boundary'
    desired_area = ${desired_area}
    add_nodes_per_boundary_segment = ${number_of_point}
    refine_boundary = false
  []

  [subdomain_intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 1
    subdomain_id_outside = 2
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = 'x-${x0}'
  []

  [break]
    type = BreakMeshByBlockGenerator
    input = subdomain_intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []

  [ss]
    type = SideSetsFromNormalsGenerator
    input = break
    normals = '-1 0 0
                1 0 0
                0 -1 0'
    fixed_normal = true
    new_boundary = 'left right bottom'
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        extra_vector_tags = 'ref'
        add_variables = true
        use_automatic_differentiation = false
        generate_output = 'stress_xx stress_xy stress_yy stress_zz strain_xx strain_xy strain_yy strain_zz'
      []
    []
  []
[]
[Physics/SolidMechanics/ShiftedCohesiveZone]
  [czm_ik]
    use_automatic_differentiation = false
    boundary = 'Block1_Block2'
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = 'signed_dist_func'
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = 'INITIAL'
    # The coarse test mesh yields surrogate distances comparable to the element size.
    suppress_distance_warning = true
  []
[]

[Materials]
  [elastic_stress]
    type = ComputeLinearElasticStress
  []
  [elasticity_tensor_in]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E1}
    block = 1
  []
  [elasticity_tensor_out]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E2}
    block = 2
  []
  [czm]
    type = SalehaniIrani3DCTraction
    normal_gap_at_maximum_normal_traction = 1
    tangential_gap_at_maximum_shear_traction = 0.5
    maximum_normal_traction = 500
    maximum_shear_traction = 300
    boundary = 'Block1_Block2'
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = 'none'

  nl_max_its = 500
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-50
  dt = ${dt}
  end_time = 5
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = 'x-${x0}'
  []
  [const_1]
    type = ConstantFunction
    value = 1
  []
[]

[AuxVariables]
  [react_x]
  []
  [react_y]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = 'ref'
    v = 'disp_x'
    variable = 'react_x'
    remove_variable_scaling = true
  []
  [react_y]
    type = TagVectorAux
    vector_tag = 'ref'
    v = 'disp_y'
    variable = 'react_y'
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
  [displacement_x_right]
    #Anchors the left side against deformation in the x-direction
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 'right'
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [area]
    type = ElementIntegralFunctorPostprocessor
    functor = 1
    block = '1 2'
    execute_on = 'final'
  []

  [n_elements]
    type = NumElements
  []

  [react_x]
    type = NodalSum
    variable = 'react_x'
    boundary = 'right'
  []
  [react_y]
    type = NodalSum
    variable = 'react_y'
    boundary = 'right'
  []
  [length]
    type = AreaPostprocessor
    boundary = 'right'
  []
[]

[Outputs]
  exodus = true
  [csv]
    type = CSV
    precision = 15
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_czm_only_quad_material_mimic.i)

E1 = 1e3
E2 = 1e3
poi = 0.3

x0 = 0.5

dt = 1

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2
    nx = 6
    ny = 6
    elem_type = QUAD4
  []

  [subdomain_intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 1
    subdomain_id_outside = 2
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = 'x-${x0}'
  []

  [break]
    type = BreakMeshByBlockGenerator
    input = subdomain_intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []

  [ss]
    type = SideSetsFromNormalsGenerator
    input = break
    normals = '-1 0 0
                1 0 0
                0 -1 0'
    fixed_normal = true
    new_boundary = 'left right bottom'
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        extra_vector_tags = 'ref'
        add_variables = true
        use_automatic_differentiation = false
        generate_output = 'stress_xx stress_xy stress_yy stress_zz strain_xx strain_xy strain_yy strain_zz'
      []
    []
  []
[]
[Physics/SolidMechanics/CohesiveZone]
  [czm_ik]
    use_automatic_differentiation = false
    boundary = 'Block1_Block2'
  []
[]

[Materials]
  [elastic_stress]
    type = ComputeLinearElasticStress
  []
  [elasticity_tensor_in]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E1}
    block = 1
  []
  [elasticity_tensor_out]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E2}
    block = 2
  []
  [czm]
    type = SalehaniIrani3DCTraction
    normal_gap_at_maximum_normal_traction = 1
    tangential_gap_at_maximum_shear_traction = 0.5
    maximum_normal_traction = 500
    maximum_shear_traction = 300
    boundary = 'Block1_Block2'
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = 'none'

  nl_max_its = 500
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-50
  dt = ${dt}
  end_time = 5
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = 'x-${x0}'
  []
  [const_1]
    type = ConstantFunction
    value = 1
  []
[]

[AuxVariables]
  [react_x]
  []
  [react_y]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = 'ref'
    v = 'disp_x'
    variable = 'react_x'
    remove_variable_scaling = true
  []
  [react_y]
    type = TagVectorAux
    vector_tag = 'ref'
    v = 'disp_y'
    variable = 'react_y'
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = 'left'
    value = 0.0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = 'bottom'
    value = 0.0
  []
  [displacement_x_right]
    #Anchors the left side against deformation in the x-direction
    type = FunctionDirichletBC
    variable = disp_x
    boundary = 'right'
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [area]
    type = ElementIntegralFunctorPostprocessor
    functor = 1
    block = '1 2'
    execute_on = 'final'
  []

  [n_elements]
    type = NumElements
  []

  [react_x]
    type = NodalSum
    variable = 'react_x'
    boundary = 'right'
  []
  [react_y]
    type = NodalSum
    variable = 'react_y'
    boundary = 'right'
  []
  [length]
    type = AreaPostprocessor
    boundary = 'right'
  []
[]

[Outputs]
  exodus = true
  [csv]
    type = CSV
    precision = 15
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_nonad_salehani.i)

nx = 8
nx_half = '${fparse nx / 2}'
E = 1e3
poi = 0.3
theta = 60
theta_rad = '${fparse theta * pi / 180}'
x0 = 0.5

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Problem]
  extra_tag_vectors = 'ref'
[]

[Mesh]
  [gen]
    type = CartesianMeshGenerator
    dim = 2
    dx = '${x0} ${x0}'
    dy = 1
    ix = '${nx_half} ${nx_half}'
    iy = '${nx}'
    subdomain_id = '1 2'
  []
  [intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 2
    subdomain_id_outside = 1
    lambda = 0.5
    is_domain_inside_surface = true
    signed_dist_function = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
  [break]
    type = BreakMeshByBlockGenerator
    input = intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = signed_dist_func
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = INITIAL
  []
[]

[Functions]
  [displacement_with_time]
    type = ParsedFunction
    expression = '1e-2*t'
  []
  [signed_dist_func]
    type = ParsedFunction
    expression = '(x-${x0})*sin(${theta_rad}) - cos(${theta_rad})*(y-${x0})'
  []
[]

[AuxVariables]
  [react_x]
  []
[]

[AuxKernels]
  [react_x]
    type = TagVectorAux
    vector_tag = ref
    v = disp_x
    variable = react_x
    remove_variable_scaling = true
  []
[]

[BCs]
  [anchor_x]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0
  []
  [anchor_y]
    type = DirichletBC
    variable = disp_y
    boundary = bottom
    value = 0
  []
  [pull]
    type = FunctionDirichletBC
    variable = disp_x
    boundary = right
    function = displacement_with_time
    preset = false
  []
[]

[Postprocessors]
  [react_x]
    type = NodalSum
    variable = react_x
    boundary = right
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = none
  nl_max_its = 50
  nl_abs_tol = 1e-12
  nl_rel_tol = 1e-10
  dt = 1
  end_time = 5
[]

[Outputs]
  [csv]
    type = CSV
    precision = 15
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        extra_vector_tags = ref
        add_variables = true
        use_automatic_differentiation = false
      []
    []
  []
[]

[Physics/SolidMechanics/ShiftedCohesiveZone]
  [czm]
    use_automatic_differentiation = false
    boundary = Block1_Block2
  []
[]

[Materials]
  [elastic_stress]
    type = ComputeLinearElasticStress
  []
  [elasticity_tensor_1]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E}
    block = 1
  []
  [elasticity_tensor_2]
    type = ComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E}
    block = 2
  []
  [czm]
    type = SalehaniIrani3DCTraction
    normal_gap_at_maximum_normal_traction = 0.1
    tangential_gap_at_maximum_shear_traction = 0.1
    # Match the AD law's normal peak and initial tangential stiffness for comparison.
    maximum_normal_traction = 183.939720585721
    maximum_shear_traction = 303.265329856317
    boundary = Block1_Block2
  []
[]

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/sczm_triangle_mms.i)

E1 = 0.1
E2 = 1
poi = 0
x0 = 0.25
a = '${fparse (pi * cos(pi * x0) * (1 - E1 / E2)) / (2 * x0)}'
b = '${fparse -pi * cos(pi * x0) * (E1 + (x0 / 2) * (1 - E1 / E2))}'

[GlobalParams]
  displacements = 'disp_x disp_y'
  use_displaced_mesh = false
  sbm_distance_uo = sbm_distance_uo
[]

[Mesh]
  [gen]
    type = FileMeshGenerator
    file = sczm_triangle_mesh.xda
  []

  [subdomain_intercepted]
    type = SubdomainInterceptedGenerator
    input = gen
    subdomain_id_inside = 1
    subdomain_id_outside = 2
    lambda = 0.5
    is_domain_inside_surface = false
    signed_dist_function = '-x+${x0}'
  []

  [break]
    type = BreakMeshByBlockGenerator
    block_pairs = '1 2'
    input = subdomain_intercepted
    split_interface = true
    add_interface_on_two_sides = true
  []

  [ss]
    type = SideSetsFromNormalsGenerator
    input = break
    normals = '-1 0 0
                1 0 0'
    fixed_normal = true
    new_boundary = 'left right'
  []
[]

[Variables]
  [disp_x]
  []
  [disp_y]
  []
[]

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        use_automatic_differentiation = true
        generate_output = 'stress_xx stress_xy stress_yy stress_zz strain_xx strain_xy strain_yy strain_zz'
      []
    []
  []
[]

[Physics/SolidMechanics/ShiftedCohesiveZone]
  [czm_ik]
    use_automatic_differentiation = true
    boundary = 'Block1_Block2'
  []
[]

[UserObjects]
  [sbm_distance_uo]
    type = BoundaryShortestDistanceToSurface
    surfaces = 'signed_dist_func'
    boundary = 'Block1_Block2'
    execution_order_group = 0
    execute_on = 'INITIAL'
  []
[]

[Materials]
  [elastic_stress]
    type = ADComputeLinearElasticStress
  []
  [elasticity_tensor1]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E1}
    block = 1
  []
  [elasticity_tensor2]
    type = ADComputeIsotropicElasticityTensor
    poissons_ratio = ${poi}
    youngs_modulus = ${E2}
    block = 2
  []
  [interface_traction]
    type = ADPureElasticTractionSeparation
    normal_stiffness = 1
    tangent_stiffness = 1
    boundary = 'Block1_Block2'
  []
[]

[AuxVariables]
  [a]
    family = MONOMIAL
  []
[]

[AuxKernels]
  [a]
    type = ConstantAux
    variable = a
    value = ${a}
  []
[]

[Kernels]
  [source_x_left]
    type = BodyForce
    variable = disp_x
    function = bx_func_left
    block = 1
  []
  [source_x_right]
    type = BodyForce
    variable = disp_x
    function = bx_func_right
    block = 2
  []
[]

[Functions]
  [signed_dist_func]
    type = ParsedFunction
    expression = '-x+${x0}'
  []
  [ux_func_left]
    type = ParsedFunction
    expression = '-sin(pi*x)'
  []

  [ux_func_right]
    type = ParsedFunction
    expression = '-sin(pi*x)+${a}*x^2+${b}'
  []

  [bx_func_left]
    type = ParsedFunction
    expression = '-pi^2*${E1}*sin(pi*x)'
  []

  [bx_func_right]
    type = ParsedFunction
    expression = '-${E2}*(2*${a} + pi^2*sin(pi*x))'
  []
[]

[BCs]
  [y_anchor]
    type = DirichletBC
    variable = disp_y
    boundary = 'left right'
    value = 0.0
  []

  [x_left]
    type = FunctionDirichletBC
    boundary = 'left'
    function = ux_func_left
    variable = disp_x
  []

  [x_right]
    type = FunctionDirichletBC
    boundary = 'right'
    function = ux_func_right
    variable = disp_x
  []
[]

[Executioner]
  type = Steady
  solve_type = NEWTON
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu superlu_dist'
  automatic_scaling = true
  line_search = 'none'
[]

[Postprocessors]
  [l2_error_left_x]
    type = ElementL2Error
    variable = disp_x
    function = ux_func_left
    block = '1'
  []

  [l2_error_right_x]
    type = ElementL2Error
    variable = disp_x
    function = ux_func_right
    block = '2'
  []

  [l2_sum]
    type = SumPostprocessor
    values = 'l2_error_left_x l2_error_right_x'
  []

  [n_elements]
    type = NumElements
  []
[]

[Outputs]
  exodus = true
[]

Overview
Classification Method
Multi-Geometry Workflows
Example Input File Syntax
Input Parameters
Input Files

SubdomainInterceptedGenerator

Overview

Classification Method

Domain Orientation

Element Occupancy and Intercepted Elements

Multi-Geometry Workflows

Example Input File Syntax

Inclined Interface Partitioning

Inverted Domain Orientation

Input Parameters

Required Parameters

Optional Parameters

Parsed Expression Advanced Parameters

Advanced Parameters

Debugging Parameters

Input Files

subdomain_id_inside

subdomain_id_outside

subdomain_id_intercepted

threshold

is_domain_inside_surface

qrule_order

subdomain_id_inside

subdomain_id_outside

mark_intercepted

subdomain_id_intercepted

mark_intercepted

lambda

multi_geo

keep_inside_as_inside

subdomain_id_inside

keep_outside_as_outside

subdomain_id_outside

modify_inside_only

subdomain_id_inside

modify_outside_only

subdomain_id_outside

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_common.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/sczm_triangle_mms.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_ad_elastic.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_common.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_sczm_tri_material.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/rectangle_czm_only_quad_material_mimic.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/angled_sczm_nonad_salehani.i)

(modules/shifted_boundary_method/test/tests/shifted_cohesive_zone/sczm_triangle_mms.i)