ShearWaveIndicator

Computes the minimum element size based on the shear wave speed.

Description

An indicator that computes the minimum element size based on the shear-wave speed.

The minimum element length ( $var element = document.getElementById("moose-equation-f1586f3b-ef24-4172-b154-2c9c94fcf060");katex.render("\\lambda_{min}", element, {displayMode:false,throwOnError:false});$ ), which this indicator computes, is computed based on the cutoff frequency ( $var element = document.getElementById("moose-equation-29de804d-0d89-4822-9803-81c355dda72d");katex.render("f", element, {displayMode:false,throwOnError:false});$ ) and shear-wave speed ( $var element = document.getElementById("moose-equation-09e439cb-632b-4842-a2ec-bc7a2c939378");katex.render("v_s", element, {displayMode:false,throwOnError:false});$ ) defined below as:

(1)var element = document.getElementById("moose-equation-78f1d827-3834-4937-ad59-2eec8d1e2ce8");katex.render("\\lambda_{min} = \\frac{v_s}{f}", element, {displayMode:true,throwOnError:false});

The shear wave speed is supplied as a material property and the cutoff frequency is defined as a constant value via the input parameters.

Input Parameters

cutoff_frequencyThe cutoff frequency in Hertz.
C++ Type:double
Unit:(no unit assumed)
Controllable:No
Description:The cutoff frequency in Hertz.

Required 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
shear_wave_speedshear_wave_speedThe name of the material properties (type Real) that computes the shear wave velocity.
Default:shear_wave_speed
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:The name of the material properties (type Real) that computes the shear wave velocity.

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
use_displaced_meshFalseWhether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.
Default:False
C++ Type:bool
Controllable:No
Description:Whether or not this object should use the displaced mesh for computation. Note that in the case this is true but no displacements are provided in the Mesh block the undisplaced mesh will still be used.

Advanced Parameters

prop_getter_suffixAn optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
C++ Type:MaterialPropertyName
Unit:(no unit assumed)
Controllable:No
Description:An optional suffix parameter that can be appended to any attempt to retrieve/get material properties. The suffix will be prepended with a '_' character.
use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.

Material Property Retrieval Parameters

Input Files

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

(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
[]

(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/indicators/shear_wave_indicator/shear_wave_indicator.i)

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 30
  ny = 30
  xmax = 300
  ymax = 300
[]

[Adaptivity]
  [./Indicators]
    [./minimum_element_size]
      type = ShearWaveIndicator
      cutoff_frequency = 25
    [../]
  [../]

[]

[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
  [../]
[]

[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
[]

[Outputs]
  execute_on = 'TIMESTEP_END'
  hide = 'u'
  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
[]

(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
[]

(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
[]

Description
Input Parameters
Input Files

Usage

Development

Demonstration