PNSInitialCondition

PNSInitialCondition can be used to initialize variables by evaluating fluid properties for porous flow calculations. It can accept both interstitial and superficial velocities for initializing superficial or interstitial velocity variables and energy.

Example input file syntax

In this example, we initialize the first component of the superficial momentum with a specified superficial velocity. The superficial velocity is specified as a vector. We could have alternatively specified an interstitial velocity and the porosity.

[superficial_rhou]
  type = PNSInitialCondition
  fluid_properties = 'fp'
  initial_pressure = ${p_initial}
  initial_temperature = ${T}
  initial_superficial_velocity = '1 0.2 18'
  porosity = porosity
  variable = 'superficial_rhou'
[]

Input Parameters

fluid_propertiesThe name of the user object for fluid properties
C++ Type:UserObjectName
Controllable:No
Description:The name of the user object for fluid properties
initial_pressureThe initial pressure, assumed constant everywhere
C++ Type:double
Controllable:No
Description:The initial pressure, assumed constant everywhere
initial_temperatureThe initial temperature, assumed constant everywhere
C++ Type:double
Controllable:No
Description:The initial temperature, assumed constant everywhere
variableThe variable this initial condition is supposed to provide values for.
C++ Type:VariableName
Controllable:No
Description:The variable this initial condition is supposed to provide values for.

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
boundaryThe list of boundaries (ids or names) from the mesh where this object applies
C++ Type:std::vector<BoundaryName>
Controllable:No
Description:The list of boundaries (ids or names) from the mesh where this object applies
initial_interstitial_velocityThe initial interstitial velocity, assumed constant everywhere
C++ Type:libMesh::VectorValue<double>
Controllable:No
Description:The initial interstitial velocity, assumed constant everywhere
initial_superficial_velocityThe initial superficial velocity, assumed constant everywhere
C++ Type:libMesh::VectorValue<double>
Controllable:No
Description:The initial superficial velocity, assumed constant everywhere
porosityPorosity variable (defaults to porosity material property).
C++ Type:std::vector<VariableName>
Controllable:No
Description:Porosity variable (defaults to porosity material property).
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
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.
variable_typeSpecifies what this variable is in the Navier Stokes namespace of variables
C++ Type:MooseEnum
Options:ht, e, Mach, pressure, rho, rhou, rhov, rhow, rho_et, specific_volume, temperature, vel_x, vel_y, vel_z, superficial_vel_x, superficial_vel_y, superficial_vel_z, superficial_rho, superficial_rhou, superficial_rhov, superficial_rhow, superficial_rho_et, superficial_rho_ht
Controllable:No
Description:Specifies what this variable is in the Navier Stokes namespace of variables

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.
ignore_uo_dependencyFalseWhen set to true, a UserObject retrieved by this IC will not be executed before the this IC
Default:False
C++ Type:bool
Controllable:No
Description:When set to true, a UserObject retrieved by this IC will not be executed before the this IC

Advanced Parameters

(moose/modules/navier_stokes/test/tests/ics/pns_test.i)

p_initial=1.01e5
T=273.15

[Mesh]
  [cartesian]
    type = GeneratedMeshGenerator
    dim = 2
    xmin = 0
    xmax = 1
    ymin = 1
    ymax = 2
    nx = 4
    ny = 4
  []
[]

[FluidProperties]
  [fp]
    type = IdealGasFluidProperties
  []
[]

[Problem]
  kernel_coverage_check = false
  solve = false
  skip_nl_system_check = true
[]

[AuxVariables]
  [porosity]
    initial_condition = 0.2
  []
  [pressure]
    type = MooseVariableFVReal
  []
  [superficial_vel_x]
    type = MooseVariableFVReal
  []
  [superficial_vel_y]
    type = MooseVariableFVReal
  []
  [superficial_vel_z]
    type = MooseVariableFVReal
  []
  [temperature]
    type = MooseVariableFVReal
  []
  [vel_x]
    type = MooseVariableFVReal
  []
  [vel_y]
    type = MooseVariableFVReal
  []
  [vel_z]
    type = MooseVariableFVReal
  []
  [superficial_rho_ht]
    type = MooseVariableFVReal
  []
  [ht]
    type = MooseVariableFVReal
  []
  [e]
    type = MooseVariableFVReal
  []
  [Mach]
    type = MooseVariableFVReal
  []
  [superficial_rho]
    type = MooseVariableFVReal
  []
  [superficial_rhou]
    type = MooseVariableFVReal
  []
  [superficial_rhov]
    type = MooseVariableFVReal
  []
  [superficial_rhow]
    type = MooseVariableFVReal
  []
  [superficial_rho_et]
    type = MooseVariableFVReal
  []
  [rho]
    type = MooseVariableFVReal
  []
  [rhou]
    type = MooseVariableFVReal
  []
  [rhov]
    type = MooseVariableFVReal
  []
  [rhow]
    type = MooseVariableFVReal
  []
  [rho_et]
    type = MooseVariableFVReal
  []
  [specific_volume]
    type = MooseVariableFVReal
  []
  [pressure_2]
  []
  [superficial_vel_x_2]
  []
  [superficial_vel_y_2]
  []
  [superficial_vel_z_2]
  []
  [vel_x_2]
  []
  [vel_y_2]
  []
  [vel_z_2]
  []
  [temperature_2]
  []
  [ht_2]
  []
  [superficial_rho_ht_2]
  []
  [e_2]
  []
  [Mach_2]
  []
  [superficial_rho_2]
  []
  [superficial_rhou_2]
  []
  [superficial_rhov_2]
  []
  [superficial_rhow_2]
  []
  [superficial_rho_et_2]
  []
  [rho_2]
  []
  [rhou_2]
  []
  [rhov_2]
  []
  [rhow_2]
  []
  [rho_et_2]
  []
  [specific_volume_2]
  []
[]

[GlobalParams]
  fluid_properties = 'fp'
  initial_pressure = ${p_initial}
  initial_temperature = ${T}
  initial_superficial_velocity = '1 0.2 18'
  porosity = porosity
[]

[ICs]
  [p]
    type = PNSInitialCondition
    variable = 'pressure'
  []
  [vel_x]
    type = PNSInitialCondition
    variable = 'vel_x'
  []
  [vel_y]
    type = PNSInitialCondition
    variable = 'vel_y'
  []
  [vel_z]
    type = PNSInitialCondition
    variable = 'vel_z'
  []
  [superficial_vel_x]
    type = PNSInitialCondition
    variable = 'superficial_vel_x'
  []
  [superficial_vel_y]
    type = PNSInitialCondition
    variable = 'superficial_vel_y'
  []
  [superficial_vel_z]
    type = PNSInitialCondition
    variable = 'superficial_vel_z'
  []
  [temperature]
    type = PNSInitialCondition
    variable = 'temperature'
  []
  [ht]
    type = PNSInitialCondition
    variable = 'ht'
  []
  [superficial_rho_ht]
    type = PNSInitialCondition
    variable = 'superficial_rho_ht'
  []
  [e]
    type = PNSInitialCondition
    variable = 'e'
  []
  [Mach]
    type = PNSInitialCondition
    variable = 'Mach'
  []
  [superficial_rho]
    type = PNSInitialCondition
    variable = 'superficial_rho'
  []
  [superficial_rhou]
    type = PNSInitialCondition
    fluid_properties = 'fp'
    initial_pressure = ${p_initial}
    initial_temperature = ${T}
    initial_superficial_velocity = '1 0.2 18'
    porosity = porosity
    variable = 'superficial_rhou'
  []
  [superficial_rhov]
    type = PNSInitialCondition
    variable = 'superficial_rhov'
  []
  [superficial_rhow]
    type = PNSInitialCondition
    variable = 'superficial_rhow'
  []
  [rho]
    type = PNSInitialCondition
    variable = 'rho'
  []
  [rhou]
    type = PNSInitialCondition
    variable = 'rhou'
  []
  [rhov]
    type = PNSInitialCondition
    variable = 'rhov'
  []
  [rhow]
    type = PNSInitialCondition
    variable = 'rhow'
  []
  [rho_et]
    type = PNSInitialCondition
    variable = 'rho_et'
  []
  [superficial_rho_et]
    type = PNSInitialCondition
    variable = 'superficial_rho_et'
  []
  [specific_volume]
    type = PNSInitialCondition
    variable = 'specific_volume'
  []
  [p_2]
    type = PNSInitialCondition
    variable = 'pressure_2'
    variable_type = 'pressure'
  []
  [superficial_vel_x_2]
    type = PNSInitialCondition
    variable = 'superficial_vel_x_2'
    variable_type = 'superficial_vel_x'
  []
  [superficial_vel_y_2]
    type = PNSInitialCondition
    variable = 'superficial_vel_y_2'
    variable_type = 'superficial_vel_y'
  []
  [superficial_vel_z_2]
    type = PNSInitialCondition
    variable = 'superficial_vel_z_2'
    variable_type = 'superficial_vel_z'
  []
  [vel_x_2]
    type = PNSInitialCondition
    variable = 'vel_x_2'
    variable_type = 'vel_x'
  []
  [vel_y_2]
    type = PNSInitialCondition
    variable = 'vel_y_2'
    variable_type = 'vel_y'
  []
  [vel_z_2]
    type = PNSInitialCondition
    variable = 'vel_z_2'
    variable_type = 'vel_z'
  []
  [temperature_2]
    type = PNSInitialCondition
    variable = 'temperature_2'
    variable_type = 'temperature'
  []
  [superficial_ht_2]
    type = PNSInitialCondition
    variable = 'superficial_rho_ht_2'
    variable_type = 'superficial_rho_ht'
  []
  [ht_2]
    type = PNSInitialCondition
    variable = 'ht_2'
    variable_type = 'ht'
  []
  [e_2]
    type = PNSInitialCondition
    variable = 'e_2'
    variable_type = 'e'
  []
  [Mach_2]
    type = PNSInitialCondition
    variable = 'Mach_2'
    variable_type = 'Mach'
  []
  [superficial_rho_2]
    type = PNSInitialCondition
    variable = 'superficial_rho_2'
    variable_type = 'superficial_rho'
  []
  [superficial_rhou_2]
    type = PNSInitialCondition
    variable = 'superficial_rhou_2'
    variable_type = 'superficial_rhou'
  []
  [superficial_rhov_2]
    type = PNSInitialCondition
    variable = 'superficial_rhov_2'
    variable_type = 'superficial_rhov'
  []
  [superficial_rhow_2]
    type = PNSInitialCondition
    variable = 'superficial_rhow_2'
    variable_type = 'superficial_rhow'
  []
  [superficial_rho_et_2]
    type = PNSInitialCondition
    variable = 'superficial_rho_et_2'
    variable_type = 'superficial_rho_et'
  []
  [rho_2]
    type = PNSInitialCondition
    variable = 'rho_2'
    variable_type = 'rho'
  []
  [rhou_2]
    type = PNSInitialCondition
    variable = 'rhou_2'
    variable_type = 'rhou'
  []
  [rhov_2]
    type = PNSInitialCondition
    variable = 'rhov_2'
    variable_type = 'rhov'
  []
  [rhow_2]
    type = PNSInitialCondition
    variable = 'rhow_2'
    variable_type = 'rhow'
  []
  [rho_et_2]
    type = PNSInitialCondition
    variable = 'rho_et_2'
    variable_type = 'rho_et'
  []
  [specific_volume_2]
    type = PNSInitialCondition
    variable = 'specific_volume_2'
    variable_type = 'specific_volume'
  []
[]

[Executioner]
  type = Steady
[]

[Outputs]
  exodus = true
[]

Example input file syntax
Input Parameters