FileMeshPhysicsComponent

This component first loads a mesh from an ExodusII file. It is equivalent to the FileMeshGenerator but can be used within a simulation with the geometry described with Components instead of a Mesh block.

This component then adds its block to the domain of definition of Physics actions. The Physics must have implemented the ::addBlocks routine.

Loading the mesh file

See the FileMeshComponent for explanations on how to load the mesh.

Defining Physics

The Physics active on the mesh loaded by this component are specified with the "physics" parameter.

Input Parameters

fileThe ExodusII mesh file name
C++ Type:FileName
Controllable:No
Description:The ExodusII mesh file name
positionTranslation vector for the file mesh [m]
C++ Type:libMesh::Point
Controllable:No
Description:Translation vector for the file mesh [m]

Required Parameters

physicsPhysics object(s) active on the Component
C++ Type:std::vector<PhysicsName>
Controllable:No
Description:Physics object(s) active on the Component

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.

Advanced Parameters

Input Files

(modules/thermal_hydraulics/test/tests/components/flow_component_ns/flow_component_ns_pinc.i)
(modules/thermal_hydraulics/test/tests/components/physics_component/single_physics.i)

(modules/thermal_hydraulics/test/tests/components/flow_component_ns/flow_component_ns_pinc.i)

# Operating conditions
u_inlet = 1
p_outlet = 0

[AuxVariables]
  [porosity]
    type = MooseVariableFVReal
    initial_condition = 0.5
  []
  [velocity_norm]
    type = MooseVariableFVReal
  []
[]

[Physics]
  [NavierStokes]
    [Flow]
      [flow]
        compressibility = 'incompressible'
        porous_medium_treatment = true

        density = 'rho'
        dynamic_viscosity = 'mu'

        initial_velocity = '${u_inlet} 0 0'
        initial_pressure = '${p_outlet}'

        mass_advection_interpolation = 'upwind'
        momentum_advection_interpolation = 'upwind'

        inlet_boundaries = 'comp1:left comp2:left'
        momentum_inlet_types = 'fixed-velocity fixed-velocity'
        momentum_inlet_functors = '${u_inlet} 0; ${u_inlet} 0'

        wall_boundaries = 'comp1:top comp1:bottom comp2:top comp2:bottom'
        momentum_wall_types = 'noslip symmetry noslip symmetry'

        outlet_boundaries = 'comp1:right comp2:right'
        momentum_outlet_types = 'fixed-pressure fixed-pressure'
        pressure_functors = '${p_outlet} ${p_outlet}'
      []
    []
  []
[]

[Components]
  [comp1]
    type = FileMeshPhysicsComponent
    file = rectangle.e
    position = '0 0 0'
    physics = 'flow'
  []

  [comp2]
    type = FileMeshPhysicsComponent
    file = rectangle.e
    position = '0 2 0'
    physics = 'flow'
  []
[]

[FunctorMaterials]
  [const_functor]
    type = ADGenericFunctorMaterial
    prop_names = 'rho mu'
    prop_values = '1  1'
  []
[]

[AuxKernels]
  [speed]
    type = ParsedAux
    variable = 'velocity_norm'
    coupled_variables = 'superficial_vel_x superficial_vel_y porosity'
    expression = 'sqrt(superficial_vel_x*superficial_vel_x + superficial_vel_y*superficial_vel_y) / porosity'
  []
[]

[Executioner]
  type = Steady
  solve_type = 'NEWTON'
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu       NONZERO               mumps'
  line_search = 'none'
  nl_rel_tol = 1e-12
  automatic_scaling = true
  off_diagonals_in_auto_scaling = true
  verbose = true
  scaling_group_variables = 'superficial_vel_x superficial_vel_y'
[]

[Debug]
  show_var_residual_norms = true
[]

# Some basic Postprocessors to examine the solution
[Postprocessors]
  [inlet-p]
    type = SideAverageValue
    variable = pressure
    boundary = 'comp1:left'
  []
  [outlet-u]
    type = SideAverageValue
    variable = superficial_vel_x
    boundary = 'comp1:right'
  []
[]

[Outputs]
  exodus = true
[]

(modules/thermal_hydraulics/test/tests/components/physics_component/single_physics.i)

# Operating conditions
u_inlet = 1
p_outlet = 0

[AuxVariables]
  [porosity]
    type = MooseVariableFVReal
    initial_condition = 0.5
  []
  [velocity_norm]
    type = MooseVariableFVReal
  []
[]

[Physics]
  [NavierStokes]
    [Flow]
      [flow]
        compressibility = 'incompressible'
        porous_medium_treatment = true

        density = 'rho'
        dynamic_viscosity = 'mu'

        initial_velocity = '${u_inlet} 0 0'
        initial_pressure = '${p_outlet}'

        mass_advection_interpolation = 'upwind'
        momentum_advection_interpolation = 'upwind'

        inlet_boundaries = 'comp1:left'
        momentum_inlet_types = 'fixed-velocity'
        momentum_inlet_functors = 'f1 f1'

        wall_boundaries = 'comp1:top comp1:bottom'
        momentum_wall_types = 'noslip symmetry'

        outlet_boundaries = 'comp1:right'
        momentum_outlet_types = 'fixed-pressure'
        pressure_functors = 'f1'
      []
    []
  []
[]

[Functions]
  [f1]
    type = ConstantFunction
    value = 1
  []
[]

[Components]
  [comp1]
    type = FileMeshPhysicsComponent
    file = rectangle.e
    position = '0 0 0'
    physics = 'flow'
  []
[]

[FunctorMaterials]
  [const_functor]
    type = ADGenericFunctorMaterial
    prop_names = 'rho mu'
    prop_values = '1  1'
  []
[]

[AuxKernels]
  [speed]
    type = ParsedAux
    variable = 'velocity_norm'
    coupled_variables = 'superficial_vel_x superficial_vel_y porosity'
    expression = 'sqrt(superficial_vel_x*superficial_vel_x + superficial_vel_y*superficial_vel_y) / porosity'
  []
[]

[Executioner]
  type = Steady
  solve_type = 'NEWTON'
  petsc_options_iname = '-pc_type -pc_factor_shift_type -pc_factor_mat_solver_type'
  petsc_options_value = 'lu       NONZERO               mumps'
  line_search = 'none'
  nl_rel_tol = 1e-12
  automatic_scaling = true
  off_diagonals_in_auto_scaling = true
  verbose = true
  scaling_group_variables = 'superficial_vel_x superficial_vel_y'
[]

[Debug]
  show_var_residual_norms = true
[]

# Some basic Postprocessors to examine the solution
[Postprocessors]
  [inlet-p]
    type = SideAverageValue
    variable = pressure
    boundary = 'comp1:left'
  []
  [outlet-u]
    type = SideAverageValue
    variable = superficial_vel_x
    boundary = 'comp1:right'
  []
[]

[Outputs]
  exodus = true
[]

Loading the mesh file
Defining Physics
Input Parameters
Input Files