MinimumElementSizeMarker

Marks element for refinement based on the element size.

Description

A marker that flags elements for refinement based on element length. The element size is provided to this object either via the "indicator" parameter or an "element_size" value. If an element has a minimum dimension greater than the size defined by either the Indicator or the element_size input parameter, then the element is marked for refinement, otherwise it is not marked.

Input Parameters

blockThe list of blocks (ids or names) that this object will be applied
C++ Type:std::vector<SubdomainName>
Controllable:No
Description:The list of blocks (ids or names) that this object will be applied
element_sizeMinimum element size to be enforced in the mesh.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Minimum element size to be enforced in the mesh.
factor0Additional factor to add to the computed minimum element size.
Default:0
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Additional factor to add to the computed minimum element size.
indicatorThe name of the Indicator that this Marker uses.
C++ Type:IndicatorName
Controllable:No
Description:The name of the Indicator that this Marker uses.
scale1Additional scale term to multiply with the computed minimum element size.
Default:1
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:Additional scale term to multiply with the computed minimum element size.

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Controllable:No
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Controllable:No
Description:Set the enabled status of the MooseObject.
outputsVector of output names where you would like to restrict the output of variables(s) associated with this object
C++ Type:std::vector<OutputName>
Controllable:No
Description:Vector of output names where you would like to restrict the output of variables(s) associated with this object

Advanced Parameters

Input Files

(test/tests/problem/interface_adapt/interface.i)
(test/tests/markers/minimum_element_size/minimum_element_size.i)
(test/tests/markers/minimum_element_size/minimum_element_size_value.i)
(test/tests/indicators/interface_and_shear_wave_indicators/interface_and_shear_wave_indicators.i)
(examples/ex11b/3d_beam_static_refine_PP.i)
(test/tests/problem/auto_adapt/auto_adapt.i)

(test/tests/problem/interface_adapt/interface.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
  xmax = 20
  ymax = 20
[]

[Adaptivity]
  initial_marker = marker
  initial_steps = 5
  max_h_level = 3
  [./Indicators]
    [./minimum_element_size]
      type = ShearWaveIndicator
      cutoff_frequency = 25
    [../]
    [./layer_interface]
      type = LayerInterfaceIndicator
      variable = layer
    [../]
  [../]
  [./Markers]
    [./minimum_element_marker]
      type = MinimumElementSizeMarker
      indicator = minimum_element_size
    [../]
    [./interface_marker]
      type = ErrorToleranceMarker
      indicator = layer_interface
      refine = 0.5
    [../]
    [./marker]
      type = ComboMarker
      markers = 'minimum_element_marker interface_marker'
    [../]
  []
[]

[Variables]
  [./u]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    boundary = left
    variable = u
    value = 1
  [../]
  [./right]
    type = DirichletBC
    boundary = right
    variable = u
    value = 0
  [../]
[]

[AuxVariables]
  [./layer]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./hmin]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./hmin]
    type = ElementLengthAux
    variable = hmin
    method = min
    execute_on = 'initial timestep_end'
  [../]
[]

[Functions]
  [./image]
    type = ImageFunction
    file = ../../images/three_layer.png
    component = 0
  [../]
[]

[ICs]
  [./layer_ic]
    type = FunctionIC
    function = image
    variable = layer
  [../]
[]

[Materials]
  [./linear]
    type = LinearSoilMaterial
    layer_variable = layer
    layer_ids     = '32     64     128'
    shear_modulus = '1.35e6 4.25e8 5.623e9'
    density       = '1500   1700   2500'
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  hide = 'u'
  exodus = true
[]

(test/tests/markers/minimum_element_size/minimum_element_size.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
  xmax = 20
  ymax = 20
[]

[Adaptivity]
  initial_marker = marker
  initial_steps = 5
  [Indicators]
    [minimum_element_size]
      type = ShearWaveIndicator
      cutoff_frequency = 25
    []
  []
  [Markers]
    [marker]
      type = MinimumElementSizeMarker
      indicator = minimum_element_size
    []
  []
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    boundary = left
    variable = u
    value = 1
  []
  [right]
    type = DirichletBC
    boundary = right
    variable = u
    value = 0
  []
[]

[AuxVariables]
  [layer]
    family = MONOMIAL
    order = CONSTANT
  []
  [hmin]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [hmin]
    type = ElementLengthAux
    variable = hmin
    method = min
    execute_on = 'initial timestep_end'
  []
[]

[Functions]
  [image]
    type = ImageFunction
    file = ../../images/three_layer.png
    component = 0
  []
[]

[ICs]
  [layer_ic]
    type = FunctionIC
    function = image
    variable = layer
  []
[]

[Materials]
  [linear]
    type = LinearSoilMaterial
    layer_variable = layer
    layer_ids     = '32     64     128'
    shear_modulus = '1.35e6 4.25e8 5.623e9'
    density       = '1500   1700   2500'
    outputs = exodus
  []
[]

[Executioner]
  type = Transient
  num_steps = 1
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  hide = 'u'
  exodus = true
[]

(test/tests/markers/minimum_element_size/minimum_element_size_value.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 5
  ny = 5
[]

[Adaptivity]
  initial_marker = marker
  initial_steps = 5
  [Markers]
    [marker]
      type = MinimumElementSizeMarker
      element_size = 0.05
    []
  []
[]

[Variables]
  [u]
  []
[]

[Kernels]
  [diff]
    type = Diffusion
    variable = u
  []
[]

[BCs]
  [left]
    type = DirichletBC
    variable = u
    boundary = left
    value = 0
  []
  [right]
    type = DirichletBC
    variable = u
    boundary = right
    value = 1
  []
[]

[Executioner]
  type = Steady
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  exodus = true
[]

(test/tests/indicators/interface_and_shear_wave_indicators/interface_and_shear_wave_indicators.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  xmax = 12000
  ymax = 200
  nx = 120 # don't make this smaller, this was the smallest I could used to reproduce the MPI error
  ny = 20  # ditto
[]

[Adaptivity]
  initial_marker = marker
  initial_steps = 2
  max_h_level = 2
  [./Indicators]
    [./minimum_element_size]
      type = ShearWaveIndicator
      cutoff_frequency = 2
    [../]
    [./layer_interface]
       type = LayerInterfaceIndicator
       variable = layer
    [../]
  [../]
  [./Markers]
    [./minimum_element_marker]
      type = MinimumElementSizeMarker
      indicator = minimum_element_size
      factor = 0.1
    [../]
    [./interface_marker]
      type = ErrorToleranceMarker
      indicator = layer_interface
      refine = 0.5
    [../]
    [./marker]
      type = ComboMarker
      markers = 'minimum_element_marker interface_marker'
    [../]
  []
[]

[Variables]
  [./u]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    boundary = 1
    variable = u
    value = 1
  [../]
  [./right]
    type = DirichletBC
    boundary = 2
    variable = u
    value = 0
  [../]
[]

[AuxVariables]
  [./layer]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./hmin]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./hmin]
    type = ElementLengthAux
    variable = hmin
    method = min
    execute_on = 'initial timestep_end'
  [../]
[]

[Functions]
  [./image]
    type = ImageFunction
    file = input.png
    component = 0
  [../]
[]

[ICs]
  [./id_ic]
    type = FunctionIC
    function = image
    variable = layer
  [../]
[]

[Materials]
  [./linear]
    type = LinearSoilMaterial
    layer_variable = layer
    layer_ids = '204 179 153 128 102 77 51 26 0'
    shear_modulus = '8.585e7 4.924e8 2.689e8 1.471e8 7.266e8 4.245e8 8.3911e8 1.26e9 1.668e9'
    density = '1600.0 1600.0 1600.0 1600.0 1700.0 1600.0 1700.0 1700.0 1700.0'
    outputs = exodus
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre    boomeramg'
[]

[Outputs]
  hide = 'u'
  exodus = true
[]

(examples/ex11b/3d_beam_static_refine_PP.i)

[Mesh]
  [mesh_gen]
    type = GeneratedMeshGenerator
    dim = 3
    nx = 10
    ny = 2
    nz = 2
    xmin = 0.0
    xmax = 5.0
    ymin = 0.0
    ymax = 1.0
    zmin = 0.0
    zmax = 1.0
    elem_type = 'HEX20'
  []
  [nodeset_gen]
    type = ExtraNodesetGenerator
    new_boundary = 'center_node'
    coord = '5.0 0.0 0.5'
    input = 'mesh_gen'
  []
  [subdomain_near]
    type = SubdomainBoundingBoxGenerator
    bottom_left = '0.0 0.0 0.0'
    top_right = '1.0 1.0 1.0'
    location = INSIDE
    block_id = 1
    input = nodeset_gen
  []
  [subdomain_middle]
    type = SubdomainBoundingBoxGenerator
    bottom_left = '1.0 0.0 0.0'
    top_right = '3.0 1.0 1.0'
    location = INSIDE
    block_id = 2
    input = subdomain_near
  []
[]

[Adaptivity]
  marker = combo
  initial_steps = 5
  [Indicators]
    [shear_wave]
      type = ShearWaveIndicator
      cutoff_frequency = 1000
    []
  []
  [Markers]
    [marker1]
      type = MinimumElementSizeMarker
      element_size = 0.15
      block = 1
    []
    [marker2]
      type = MinimumElementSizeMarker
      indicator = shear_wave
      block = 2
    []
    [combo]
      type = ComboMarker
      markers = 'marker1 marker2'
      block = '1 2'
    []
  []
[]

[GlobalParams]
  displacements = 'disp_x disp_y disp_z'
[]

[Physics/SolidMechanics/QuasiStatic]
  [strain_calculator]
    strain = FINITE
    generate_output = 'stress_xx'
  []
[]

[Variables]
  [disp_x]
    order = SECOND
    family = LAGRANGE
  []
  [disp_y]
    order = SECOND
    family = LAGRANGE
  []
  [disp_z]
    order = SECOND
    family = LAGRANGE
  []
[]

[BCs]
  [fixx0]
    type = DirichletBC
    variable = disp_x
    boundary = left
    value = 0.0
  []
  [fixy1]
    type = DirichletBC
    variable = disp_y
    boundary = left
    value = 0.0
  []
  [fixy2]
    type = DirichletBC
    variable = disp_z
    boundary = left
    value = 0.0
  []
[]

[NodalKernels]
  [force_y2]
    type = ConstantRate
    variable = disp_y
    boundary = 'center_node'
    rate = -1000
  []
[]

[Materials]
  [elasticity]
    type = ComputeIsotropicElasticityTensor
    youngs_modulus = 10e9
    poissons_ratio = 0.3
  []
  [stress]
    type = ComputeFiniteStrainElasticStress
  []
  [shear_wave_speed]
    type = GenericConstantMaterial
    prop_names = shear_wave_speed
    prop_values =  250.0
  []
[]

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

[Executioner]
  type = Steady
  solve_type = NEWTON
  petsc_options = '-ksp_snes_ew'
  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu       superlu_dist'
  line_search = 'none'
  nl_rel_tol = 1e-6
  nl_abs_tol = 1e-10
[]

[Postprocessors]
  [moment_z]
    type = SidesetMoment
    stress_direction = '1 0 0'
    stress_tensor = stress
    boundary = 'left'
    reference_point = '0.0 0.0 0.5'
    moment_direction = '0 0 1'
  []
  [avg_stress_xx_side]
    type = SideAverageValue
    variable = stress_xx
    boundary = left
  []
  [avg_stress_xx_block]
    type = ElementAverageValue
    variable = stress_xx
    block = 1
  []
[]

[Outputs]
  file_base = 11_b_out
  exodus = true
  csv = true
  perf_graph = true
[]

(test/tests/problem/auto_adapt/auto_adapt.i)

# @requirement F3.2
[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 10
  ny = 10
  xmax = 20
  ymax = 20
[]

[Adaptivity]
  initial_marker = marker
  initial_steps = 5
  [./Indicators]
    [./minimum_element_size]
      type = ShearWaveIndicator
      cutoff_frequency = 25
    [../]
  [../]
  [./Markers]
    [./marker]
      type = MinimumElementSizeMarker
      indicator = minimum_element_size
    [../]
  []
[]

[Variables]
  [./u]
  [../]
[]

[Kernels]
  [./diff]
    type = Diffusion
    variable = u
  [../]
[]

[BCs]
  [./left]
    type = DirichletBC
    boundary = left
    variable = u
    value = 1
  [../]
  [./right]
    type = DirichletBC
    boundary = right
    variable = u
    value = 0
  [../]
[]

[AuxVariables]
  [./layer]
    family = MONOMIAL
    order = CONSTANT
  [../]
  [./hmin]
    family = MONOMIAL
    order = CONSTANT
  [../]
[]

[AuxKernels]
  [./hmin]
    type = ElementLengthAux
    variable = hmin
    method = min
    execute_on = 'initial timestep_end'
  [../]
[]

[Functions]
  [./image]
    type = ImageFunction
    file = ../../images/three_layer.png
    component = 0
  [../]
[]

[ICs]
  [./layer_ic]
    type = FunctionIC
    function = image
    variable = layer
  [../]
[]

[Materials]
  [./linear]
    type = LinearSoilMaterial
    layer_variable = layer
    layer_ids     = '32     64     128'
    shear_modulus = '1.35e6 4.25e8 5.623e5'
    density       = '1500   1700   2500'
  [../]
[]

[Executioner]
  type = Transient
  num_steps = 1
  solve_type = 'PJFNK'
  petsc_options_iname = '-pc_type -pc_hypre_type'
  petsc_options_value = 'hypre boomeramg'
[]

[Outputs]
  hide = 'u'
  exodus = true
[]

Description
Input Parameters
Input Files

Usage

Development

Demonstration