Navier Stokes Solid Heat Transfer / PNSFVSolidHeatTransfer

Define the Navier Stokes porous media solid energy equation

Equation

This Physics object creates the kernels and boundary conditions to solve the heat transfer equation for the solid phase in a porous media flow problem.

(1ϵ)ρhst(κsTs)α(TfTs)(1ϵ)Q˙=0\dfrac{\partial (1-\epsilon) \rho h_s}{\partial t} - \nabla \cdot (\kappa_s \nabla T_s) - \alpha (T_f - T_s) - (1-\epsilon) \dot{Q} = 0

where:

  • hsh_s is the solid specific enthalpy, computed from the specific heat cpsc_{ps} and the solid temperature

  • ρ\rho is the solid density

  • ϵ\epsilon is the porosity

  • TsT_s is the solid temperature

  • TfT_f is the fluid temperature

  • ksk_s the solid thermal conductivity

  • (1ϵ)Q(1-\epsilon) Q is the source term, corresponding to energy deposited directly in the solid phase

  • α\alpha is the ambient convection volumetric heat transfer coefficient

The enthalpy is used in lieu of ρcpT\rho c_p T to be able to model solids with temperature dependent specific heat.

The kernels potentially created for this equation are:

commentnote

Additional details on porous media flow equations can be found on this page.

Automatically created variables

The PNSFVSolidHeatTransferPhysics by default will automatically create a nonlinear variable for the solid phase temperature. It should be named as follows:

static const std::string T_solidus = "T_solidus";
(contrib/moose/modules/navier_stokes/include/base/NS.h)

Coupling with other Physics

In the presence of fluid flow, a Navier Stokes Fluid Heat Transfer / WCNSFVFluidHeatTransferPhysics should be created using the [Physics/NavierStokes/FluidHeatTransfer/<name>] syntax. The following input performs the coupling between the fluid equations and the solid temperature equations. The coupling between the fluid and solid domain is performed through a volumetric ambient convection term.

[Physics]
  [NavierStokes]
    [Flow]
      [flow]
        compressibility = 'weakly-compressible'
        porous_medium_treatment = true
        define_variables = true

        pressure_variable = 'pressure'

        density = 'rho'
        dynamic_viscosity = 'mu'

        initial_velocity = '${u_inlet} 1e-6 0'
        initial_pressure = '${p_outlet}'

        inlet_boundaries = 'left'
        momentum_inlet_types = 'fixed-velocity'
        momentum_inlet_function = '${u_inlet} 0'

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

        outlet_boundaries = 'right'
        momentum_outlet_types = 'fixed-pressure'
        pressure_function = '${p_outlet}'

        mass_advection_interpolation = 'average'
        momentum_advection_interpolation = 'average'
      []
    []
    [FluidHeatTransfer]
      [fluid]
        thermal_conductivity = 'k'
        effective_conductivity = true
        specific_heat = 'cp'

        initial_temperature = '${T_inlet}'

        # See 'flow' for inlet boundaries
        energy_inlet_types = 'fixed-temperature'
        energy_inlet_function = '${T_inlet}'

        # See 'flow' for wall boundaries
        energy_wall_types = 'heatflux heatflux'
        energy_wall_function = '0 0'

        ambient_convection_alpha = 'h_cv'
        ambient_temperature = 'T_solid'

        energy_advection_interpolation = 'average'
      []
    []

    [SolidHeatTransfer]
      [solid]
        block = 0
        initial_temperature = 100
        transient = true
        # To match the previous test results
        solid_temperature_two_term_bc_expansion = true

        thermal_conductivity_solid = '${k_s}'
        cp_solid = ${cp_s}
        rho_solid = ${rho_s}

        fixed_temperature_boundaries = 'top'
        boundary_temperatures = '${top_side_temperature}'

        ambient_convection_alpha = 'h_cv'
        ambient_convection_temperature = 'T_fluid'
        verbose = true
      []
    []
  []
[]
(contrib/moose/modules/navier_stokes/test/tests/finite_volume/pwcns/channel-flow/2d-transient-physics.i)
warningwarning

Conjugate heat transfer on a surface on the boundary of the fluid domain is not currently implemented with the Physics syntax. Please use a FVConvectionCorrelationInterface for that purpose.

Input Parameters

  • blockBlocks (subdomains) that this Physics is active on.

    C++ Type:std::vector<SubdomainName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Blocks (subdomains) that this Physics is active on.

  • fluid_temperature_variableT_fluidName of the fluid temperature variable

    Default:T_fluid

    C++ Type:NonlinearVariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:Name of the fluid temperature variable

  • heat_source_functorFunctor providing the heat source. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Functor providing the heat source. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • initial_temperature300Initial value of the temperature variable

    Default:300

    C++ Type:FunctionName

    Unit:(no unit assumed)

    Controllable:No

    Description:Initial value of the temperature variable

  • porosityporosityThe name of the auxiliary variable for the porosity field. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:porosity

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The name of the auxiliary variable for the porosity field. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • solid_temperature_variableT_solidName of the solid phase temperature variable

    Default:T_solid

    C++ Type:VariableName

    Unit:(no unit assumed)

    Controllable:No

    Description:Name of the solid phase temperature variable

  • temperature_scaling1Scaling factor for the heat conduction equation

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Scaling factor for the heat conduction equation

  • transientsame_as_problemWhether the physics is to be solved as a transient

    Default:same_as_problem

    C++ Type:MooseEnum

    Unit:(no unit assumed)

    Options:true, false, same_as_problem

    Controllable:No

    Description:Whether the physics is to be solved as a transient

  • verboseFalseFlag to facilitate debugging a Physics

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Flag to facilitate debugging a Physics

Optional Parameters

  • active__all__ If specified only the blocks named will be visited and made active

    Default:__all__

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:If specified only the blocks named will be visited and made active

  • control_tagsAdds user-defined labels for accessing object parameters via control logic.

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:Adds user-defined labels for accessing object parameters via control logic.

  • ghost_layers2Number of layers of elements to ghost near process domain boundaries

    Default:2

    C++ Type:unsigned short

    Unit:(no unit assumed)

    Controllable:No

    Description:Number of layers of elements to ghost near process domain boundaries

  • inactiveIf specified blocks matching these identifiers will be skipped.

    C++ Type:std::vector<std::string>

    Unit:(no unit assumed)

    Controllable:No

    Description:If specified blocks matching these identifiers will be skipped.

Advanced Parameters

  • ambient_convection_alphaThe heat exchange coefficients for each block in 'ambient_convection_blocks'.

    C++ Type:std::vector<MooseFunctorName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The heat exchange coefficients for each block in 'ambient_convection_blocks'.

  • ambient_convection_blocksThe blocks where the ambient convection is present.

    C++ Type:std::vector<std::vector<SubdomainName>>

    Unit:(no unit assumed)

    Controllable:No

    Description:The blocks where the ambient convection is present.

  • ambient_convection_temperatureT_fluid The fluid temperature for each block in 'ambient_convection_blocks'.

    Default:T_fluid

    C++ Type:std::vector<MooseFunctorName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The fluid temperature for each block in 'ambient_convection_blocks'.

Ambient Convection Parameters

  • boundary_heat_fluxesFunctors to compute the heat flux on each boundary in 'heat_flux_boundaries'

    C++ Type:std::vector<MooseFunctorName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Functors to compute the heat flux on each boundary in 'heat_flux_boundaries'

  • boundary_temperaturesFunctors to compute the heat flux on each boundary in 'fixed_temperature_boundaries'

    C++ Type:std::vector<MooseFunctorName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Functors to compute the heat flux on each boundary in 'fixed_temperature_boundaries'

  • fixed_temperature_boundariesBoundaries on which to apply a fixed temperature

    C++ Type:std::vector<BoundaryName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Boundaries on which to apply a fixed temperature

  • heat_flux_boundariesBoundaries on which to apply a heat flux

    C++ Type:std::vector<BoundaryName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Boundaries on which to apply a heat flux

  • insulated_boundariesBoundaries on which to apply a zero heat flux

    C++ Type:std::vector<BoundaryName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Boundaries on which to apply a zero heat flux

Thermal Boundaries Parameters

  • cp_solidcp_solidSpecific heat functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:cp_solid

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Specific heat functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • rho_solidrho_solidDensity functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    Default:rho_solid

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:Density functor. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • thermal_conductivity_blocksBlocks which each thermal conductivity is defined

    C++ Type:std::vector<std::vector<SubdomainName>>

    Unit:(no unit assumed)

    Controllable:No

    Description:Blocks which each thermal conductivity is defined

  • thermal_conductivity_solidThermal conductivity, which may have different names depending on the subdomain

    C++ Type:std::vector<MooseFunctorName>

    Unit:(no unit assumed)

    Controllable:No

    Description:Thermal conductivity, which may have different names depending on the subdomain

  • use_external_enthalpy_materialFalseTo indicate if the enthalpy material is set up outside of the action.

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:To indicate if the enthalpy material is set up outside of the action.

Material Properties Parameters

  • external_heat_sourceThe name of a functor which contains the external heat source for the energy equation. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

    C++ Type:MooseFunctorName

    Unit:(no unit assumed)

    Controllable:No

    Description:The name of a functor which contains the external heat source for the energy equation. A functor is any of the following: a variable, a functor material property, a function, a post-processor, or a number.

  • external_heat_source_blocksThe blocks where the heat source is present.

    C++ Type:std::vector<SubdomainName>

    Unit:(no unit assumed)

    Controllable:No

    Description:The blocks where the heat source is present.

  • external_heat_source_coeff1Multiplier for the coupled heat source term.

    Default:1

    C++ Type:double

    Unit:(no unit assumed)

    Controllable:No

    Description:Multiplier for the coupled heat source term.

Solid Porous Medium Heat Source Parameters

  • initial_from_file_timestepLATESTGives the time step number (or "LATEST") for which to read the Exodus solution

    Default:LATEST

    C++ Type:std::string

    Unit:(no unit assumed)

    Controllable:No

    Description:Gives the time step number (or "LATEST") for which to read the Exodus solution

  • initialize_variables_from_mesh_fileFalseDetermines if the variables that are added by the action are initializedfrom the mesh file (only for Exodus format)

    Default:False

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:Determines if the variables that are added by the action are initializedfrom the mesh file (only for Exodus format)

Restart From Exodus Parameters

  • solid_temperature_face_interpolationaverageThe numerical scheme to interpolate the temperature/energy to the face for conduction (separate from the advected quantity interpolation).

    Default:average

    C++ Type:MooseEnum

    Unit:(no unit assumed)

    Options:average, skewness-corrected

    Controllable:No

    Description:The numerical scheme to interpolate the temperature/energy to the face for conduction (separate from the advected quantity interpolation).

  • solid_temperature_two_term_bc_expansionTrueIf a two-term Taylor expansion is needed for the determination of the boundary valuesof the temperature/energy.

    Default:True

    C++ Type:bool

    Unit:(no unit assumed)

    Controllable:No

    Description:If a two-term Taylor expansion is needed for the determination of the boundary valuesof the temperature/energy.

Numerical Scheme Parameters