SpecifiedVaporPressureBC

Prescribed vapor pressure boundary condition for moisture transport in concrete.

Description

This boundary condition provides the ability to prescribe the vapor pressure at a location to influence the solution for the relative humidity in a concrete moisture transport simulation. It computes a saturated vapor pressure, $var element = document.getElementById("moose-equation-52faac87-927e-4ad5-9d4d-cca6bf634a46");katex.render("P_{vs}", element, {displayMode:false,throwOnError:false});$ , defined by Bary et al. (2012) as:

(1)var element = document.getElementById("moose-equation-76d2e523-6f17-4a57-ad82-75576b06da44");katex.render("P_{vs} = P_{atm}\\exp\\left(4871.3\\frac{T-100}{373.15T}\\right)", element, {displayMode:true,throwOnError:false});

where

$var element = document.getElementById("moose-equation-250967c9-871c-466d-894b-ac3e07c38d6a");katex.render("T", element, {displayMode:false,throwOnError:false});$ = temperature in K
$var element = document.getElementById("moose-equation-ad1e2842-a2fd-4026-b3ee-9ee4574d7511");katex.render("P_{atm}", element, {displayMode:false,throwOnError:false});$ = standard atmospheric pressure (101325 Pa)

Note that although $var element = document.getElementById("moose-equation-6e202a5b-0837-4fda-a6c6-2c519c334857");katex.render("T", element, {displayMode:false,throwOnError:false});$ in the above formula is in K, the model is designed to accept T in units of $var element = document.getElementById("moose-equation-c6753431-591f-4b8c-a445-652f58008539");katex.render("\\degree", element, {displayMode:false,throwOnError:false});$ C and convert to K internally. This saturated vapor pressure is then used to compute a prescribed relative humidity, $var element = document.getElementById("moose-equation-c2904106-2fcc-4f44-b49b-ad9a1cfd2382");katex.render("P_{v,pres}", element, {displayMode:false,throwOnError:false});$ , which is prescribed using a Dirichlet-style boundary condition:

(2)var element = document.getElementById("moose-equation-dc67299d-f23c-4009-bfa3-b299bbd09102");katex.render("rh_{pres}=P_{v,pres}/P_{vs}", element, {displayMode:true,throwOnError:false});

To avoid sudden changes, this is ramped from initial prescribed conditions over a user-specified duration.

Input Parameters

boundaryThe list of boundary IDs from the mesh where this boundary condition applies
C++ Type:std::vector<BoundaryName>
Options:
Description:The list of boundary IDs from the mesh where this boundary condition applies
vapor_pressurein Pa
C++ Type:double
Options:
Description:in Pa
variableThe name of the variable that this residual object operates on
C++ Type:NonlinearVariableName
Options:
Description:The name of the variable that this residual object operates on

Required Parameters

T_ref20Initial temperature in C
Default:20
C++ Type:double
Options:
Description:Initial temperature in C
diag_save_inThe name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Options:
Description:The name of auxiliary variables to save this BC's diagonal jacobian contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
displacementsThe displacements
C++ Type:std::vector<VariableName>
Options:
Description:The displacements
duration0
Default:0
C++ Type:double
Options:
rh_ref0.96initial relative humidity
Default:0.96
C++ Type:double
Options:
Description:initial relative humidity
save_inThe name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
C++ Type:std::vector<AuxVariableName>
Options:
Description:The name of auxiliary variables to save this BC's residual contributions to. Everything about that variable must match everything about this variable (the type, what blocks it's on, etc.)
temperaturenonlinear variable name holding temperature field
C++ Type:std::vector<VariableName>
Options:
Description:nonlinear variable name holding temperature field

Optional Parameters

control_tagsAdds user-defined labels for accessing object parameters via control logic.
C++ Type:std::vector<std::string>
Options:
Description:Adds user-defined labels for accessing object parameters via control logic.
enableTrueSet the enabled status of the MooseObject.
Default:True
C++ Type:bool
Options:
Description:Set the enabled status of the MooseObject.
implicitTrueDetermines whether this object is calculated using an implicit or explicit form
Default:True
C++ Type:bool
Options:
Description:Determines whether this object is calculated using an implicit or explicit form
seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Options:
Description:The seed for the master random number generator
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
Options:
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

extra_matrix_tagsThe extra tags for the matrices this Kernel should fill
C++ Type:std::vector<TagName>
Options:
Description:The extra tags for the matrices this Kernel should fill
extra_vector_tagsThe extra tags for the vectors this Kernel should fill
C++ Type:std::vector<TagName>
Options:
Description:The extra tags for the vectors this Kernel should fill
matrix_tagssystem timeThe tag for the matrices this Kernel should fill
Default:system time
C++ Type:MultiMooseEnum
Options:nontime, system, time
Description:The tag for the matrices this Kernel should fill
vector_tagsresidualThe tag for the vectors this Kernel should fill
Default:residual
C++ Type:MultiMooseEnum
Options:nontime, time, residual
Description:The tag for the vectors this Kernel should fill

Tagging Parameters

Input Files

(test/tests/concrete_moisture_heat_transfer/maqbeth_1d.i)

References

Beno\^ıt Bary, Marcus V.G. de Morais, Stéphane Poyet, and Sabine Durand. Simulations of the thermo-hydro-mechanical behaviour of an annular reinforced concrete structure heated up to 200°c. Engineering Structures, 36:302–315, March 2012.

@article{Bary2012,
    author = "Bary, Beno{\^{\i}}t and de Morais, Marcus V.G. and Poyet, St{\'{e}}phane and Durand, Sabine",
    journal = "Engineering Structures",
    month = "March",
    pages = "302--315",
    title = "Simulations of the thermo-hydro-mechanical behaviour of an annular reinforced concrete structure heated up to 200{\textdegree}C",
    volume = "36",
    year = "2012"
}

(test/tests/concrete_moisture_heat_transfer/maqbeth_1d.i)

[Problem]
  type = ReferenceResidualProblem
  reference_vector = 'ref'
  extra_tag_vectors = 'ref'
  coord_type = RZ
[]

[Mesh]
  file = maqbeth_1d.e
  construct_side_list_from_node_list = true
[]

[Functions]
  [temp_hist]
    type = PiecewiseLinear
    x = '0.0  36000.0 114120.0 138600.0  213120'
    y = '20.0  20.0      150.0     80.0   200.0'
  []
[]

[Variables]
  [T]
    order = FIRST
    family = LAGRANGE
    initial_condition = 20.0
  []
  [rh]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.96
  []
[]

[AuxVariables]
  [thermal_conductivity]
    order = CONSTANT
    family = Monomial
  []

  [thermal_capacity]
    order = CONSTANT
    family = Monomial
  []

  [moisture_capacity]
    order = CONSTANT
    family = Monomial
  []

  [humidity_diffusivity]
    order = CONSTANT
    family = Monomial
  []

  [water_content]
    order = CONSTANT
    family = Monomial
  []
  [water_hydrated]
    order = CONSTANT
    family = Monomial
  []
[]

[Kernels]
  [T_td]
    type     = ConcreteThermalTimeIntegration
    variable = T
    extra_vector_tags = 'ref'
  []
  [T_diff]
    type     = ConcreteThermalConduction
    variable = T
    extra_vector_tags = 'ref'
  []

  [T_conv]
    type     = ConcreteThermalConvection
    variable = T
    relative_humidity = rh
    extra_vector_tags = 'ref'
  []

  [T_adsorption]
    type     = ConcreteLatentHeat
    variable = T
    H = rh
    extra_vector_tags = 'ref'
  []

  [rh_td]
    type     = ConcreteMoistureTimeIntegration
    variable = rh
    extra_vector_tags = 'ref'
  []

  [h_diff]
    type     = ConcreteMoistureDiffusion
    variable = rh
    temperature = T
    extra_vector_tags = 'ref'
  []
[]

[AuxKernels]
  [k]
    type = MaterialRealAux
    variable = thermal_conductivity
    property = thermal_conductivity
    execute_on = 'timestep_end'
  [../]
  [./capacity]
    type = MaterialRealAux
    variable = thermal_capacity
    property = thermal_capacity
    execute_on = 'timestep_end'
  [../]

  [./rh_capacity]
    type = MaterialRealAux
    variable = moisture_capacity
    property = moisture_capacity
    execute_on = 'timestep_end'
  [../]
  [./rh_diff]
    type = MaterialRealAux
    variable = humidity_diffusivity
    property = humidity_diffusivity
    execute_on = 'timestep_end'
  [../]
  [./wc_diff]
    type = MaterialRealAux
    variable = water_content
    property = moisture_content
    execute_on = 'timestep_end'
  [../]
  [./hydrw_duff]
    type = MaterialRealAux
    variable = water_hydrated
    property = hydrated_water
    execute_on = 'timestep_end'
  [../]
[]

[Materials]
  [./concrete]
    type = ConcreteThermalMoisture
    block = 1
    # options available: CONSTANT ASCE-1992 KODUR-2004 EUROCODE-2004 KIM-2003
    thermal_conductivity_model =  KODUR-2004
    thermal_capacity_model     =  KODUR-2004
    aggregate_type = Siliceous               #options: Siliceous Carbonate

    ref_density_of_concrete = 2231.0         # in kg/m^3
    ref_specific_heat_of_concrete = 1100.0   # in J/(Kg.0C)
    ref_thermal_conductivity_of_concrete = 3 # in W/(m.0C)


    # setup moisture capacity and humidity diffusivity models
    aggregate_pore_type = dense              #options: dense porous
    aggregate_mass = 1877.0                  #mass of aggregate (kg) per m^3 of concrete
    cement_type = 1                          #options: 1 2 3 4
    cement_mass = 354.0                      #mass of cement (kg) per m^3 of concrete
    water_to_cement_ratio       = 0.43
    concrete_cure_time          = 23.0       #curing time in (days)

    # options available for humidity diffusivity:
    moisture_diffusivity_model =  Bazant      #options: Bazant Xi Mensi
    D1 = 3.0e-10

    coupled_moisture_diffusivity_factor = 1.0e-3  # factor for mositure diffusivity due to heat

    # coupled nonlinear variables
    relative_humidity = rh
    temperature = T
  [../]
[]

[BCs]
  [./T_left]
    type = FunctionDirichletBC
    variable = T
    boundary = '1'
    function = temp_hist
  [../]

  [./T_right]
    type = ConvectiveFluxBC
    variable = T
    boundary = '2'
    final = 20.0
    rate = 10.0
  [../]

  [./rh_left]
    type = SpecifiedVaporPressureBC
    variable = rh
    boundary = '1'
    duration = 3600
    vapor_pressure = 2500.0
    temperature = T
  [../]
  [./rh_right]
    type = SpecifiedVaporPressureBC
    variable = rh
    boundary = '2'
    duration = 3600
    vapor_pressure = 2500.0
    temperature = T
  [../]
[]

[VectorPostprocessors]
  [profiles]
    type = LineValueSampler
    start_point = '0.5 0 0'
    end_point = '1.1 0 0'
    num_points = 100
    variable = 'T rh'
    sort_by = id
    outputs = csv
  []
[]

[Preconditioning]
  [smp]
    type = SMP
    full = true
  []
[]

[Executioner]
  type       = Transient
  solve_type = 'PJFNK'

  petsc_options_iname = '-pc_type -pc_factor_mat_solver_package'
  petsc_options_value = 'lu superlu_dist'

  automatic_scaling = true

  dt = 10000

  end_time = 900000.0

  l_max_its  = 50
  l_tol      = 1e-6
  nl_max_its = 10
  nl_rel_tol = 1e-6
  nl_abs_tol = 1e-10

  [Predictor]
    type = SimplePredictor
    scale = 1.0
    skip_times_old = '0.0  36000.0 114120.0 138600.0  213120'
  []
[]

[Outputs]
  exodus     = true
  [csv]
    type = CSV
    file_base = 'csv/out'
    time_data = true
  []
[]

Description
Input Parameters
Input Files
References