Diffusion Physics syntax

The [Diffusion] block is a sub-block of the [Physics] block. DiffusionPhysicsBase-derived classes, such as DiffusionCG can be created, nested under the relevant discretization sub-block.

For example, a DiffusionCG can be created inside the Physics/Diffusion/ContinuousGalerkin block.

[Physics<<<{"href": "../index.html"}>>>]
  [Diffusion<<<{"href": "index.html"}>>>]
    [ContinuousGalerkin<<<{"href": "ContinuousGalerkin/index.html"}>>>]
      [diff]
        source_functor<<<{"description": "Source term in the diffusion problem. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number."}>>> = 2
        diffusivity_matprop<<<{"description": "Material property defining the diffusion coefficient"}>>> = '1'

        # Test all the ways of setting the boundary conditions
        neumann_boundaries<<<{"description": "Boundaries on which to apply a diffusive flux"}>>> = 'left_to_0 right_to_0 top_to_0 bottom_to_0'
        boundary_fluxes<<<{"description": "Functors to compute the diffusive flux on each Neumann boundary'"}>>> = '1 flux_pp flux_function flux_variable'
        dirichlet_boundaries<<<{"description": "Boundaries on which to apply a fixed value"}>>> = 'left_to_1 right_to_1 top_to_1 bottom_to_1'
        boundary_values<<<{"description": "Functors to compute the diffusive flux on each Dirichlet boundary'"}>>> = '2 value_pp value_function value_variable'
      []
    []
  []
[]

(moose/test/tests/physics/diffusion_cg.i)

[Mesh]
  [cmg]
    type = CartesianMeshGenerator
    dim = 2
    dx = '1 2'
    dy = '2 1'
    ix = '2 3'
    iy = '3 2'
    subdomain_id = '0 1
                    1 0'
  []
  [split_boundaries]
    type = BreakBoundaryOnSubdomainGenerator
    input = cmg
  []
  allow_renumbering = false
[]

[Physics]
  [Diffusion]
    [ContinuousGalerkin]
      [diff]
        source_functor = 2
        diffusivity_matprop = '1'

        # Test all the ways of setting the boundary conditions
        neumann_boundaries = 'left_to_0 right_to_0 top_to_0 bottom_to_0'
        boundary_fluxes = '1 flux_pp flux_function flux_variable'
        dirichlet_boundaries = 'left_to_1 right_to_1 top_to_1 bottom_to_1'
        boundary_values = '2 value_pp value_function value_variable'
      []
    []
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  num_steps = 10
  # Output the setup
  verbose = true
[]

# To test setting up a boundary condition with a postprocessor
[Postprocessors]
  [flux_pp]
    type = Receiver
    default = 1
    outputs = 'none'
  []
  [value_pp]
    type = Receiver
    default = 2
    outputs = 'none'
  []
[]

# To test setting up a boundary condition with a function
[Functions]
  [flux_function]
    type = ConstantFunction
    value = 1
  []
  [value_function]
    type = ConstantFunction
    value = 2
  []
[]

# To test setting up a boundary condition with a variable
[AuxVariables]
  [flux_variable]
    initial_condition = 1
  []
  [value_variable]
    initial_condition = 2
  []
[]

# Form output for testing
[VectorPostprocessors]
  [sample]
    type = NodalValueSampler
    variable = 'u'
    sort_by = 'id'
  []
[]

[Outputs]
  csv = true
[]