- variableThe name of the variable that this object applies to
C++ Type:AuxVariableName
Unit:(no unit assumed)
Controllable:No
Description:The name of the variable that this object applies to
PTFluidVelocityAux
buildconstruction:Undocumented Class
The PTFluidVelocityAux has not been documented. The content listed below should be used as a starting point for documenting the class, which includes the typical automatic documentation associated with a MooseObject; however, what is contained is ultimately determined by what is necessary to make the documentation clear for users.
Overview
Example Input File Syntax
Input 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
- check_boundary_restrictedTrueWhether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
Default:True
C++ Type:bool
Controllable:No
Description:Whether to check for multiple element sides on the boundary in the case of a boundary restricted, element aux variable. Setting this to false will allow contribution to a single element's elemental value(s) from multiple boundary sides on the same element (example: when the restricted boundary exists on two or more sides of an element, such as at a corner of a mesh
- component0component of the Darcy flux vector
Default:0
C++ Type:int
Controllable:No
Description:component of the Darcy flux vector
- execute_onLINEAR TIMESTEP_ENDThe list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
Default:LINEAR TIMESTEP_END
C++ Type:ExecFlagEnum
Options:NONE, INITIAL, LINEAR, LINEAR_CONVERGENCE, NONLINEAR, NONLINEAR_CONVERGENCE, POSTCHECK, TIMESTEP_END, TIMESTEP_BEGIN, MULTIAPP_FIXED_POINT_END, MULTIAPP_FIXED_POINT_BEGIN, MULTIAPP_FIXED_POINT_CONVERGENCE, FINAL, CUSTOM, PRE_DISPLACE
Controllable:No
Description:The list of flag(s) indicating when this object should be executed. For a description of each flag, see https://mooseframework.inl.gov/source/interfaces/SetupInterface.html.
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:Yes
Description:Set the enabled status of the MooseObject.
- seed0The seed for the master random number generator
Default:0
C++ Type:unsigned int
Controllable:No
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
Controllable:No
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
- 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
Unit:(no unit assumed)
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.
- use_interpolated_stateFalseFor the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Default:False
C++ Type:bool
Controllable:No
Description:For the old and older state use projected material properties interpolated at the quadrature points. To set up projection use the ProjectedStatefulMaterialStorageAction.
Material Property Retrieval Parameters
Input Files
- (test/tests/PT_H_stochastic_permeability/inp.i)
- (test/tests/PT_TH_elder/step2_wseos.i)
- (test/tests/PT_H_SideFluxAverage/inp.i)
- (test/tests/PT_H_permeability_function/test2.i)
- (test/tests/PT_H_gravity/inp.i)
- (test/tests/PT_TH_SGW2016_Reference_Case/inp3d_run.i)
- (test/tests/PT_TH_elder/step2_comp.i)
- (test/tests/PT_TH_SGW2016_Reference_Case/inp3d_ini.i)
- (test/tests/PT_H_SideFluxIntegral/inp.i)
- (test/tests/PT_H_permeability_function/test1.i)
(test/tests/PT_H_stochastic_permeability/inp.i)
### This is a primitive implementation of
### stochastic field for permeability.
### A structured hex mesh must be used to be compatible with
### the (x,y,z) location of cell-centered permeability.
############################################################
[Mesh]
file = mesh.e
[]
############################################################
[Variables]
[./P]
initial_condition = 1e6
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./k]
order = CONSTANT
family = MONOMIAL
[../]
[./stok]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./k]
type = PTPermeabilityAux
variable = k
[../]
[./stok]
type = StochasticFieldAux
variable = stok
file_name = IESE_Permeability.txt
[../]
[]
############################################################
[BCs]
[./P_right]
type = DirichletBC
variable = P
boundary = '1'
value = 10e6
[../]
[./P_left]
type = DirichletBC
variable = P
boundary = '2'
value = 1e6
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '1'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
stochastic_permeability = stok
permeability = 1e-12
porosity = 1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 0.0
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-3
nl_max_its = 500
nl_abs_tol = 1e-9
[]
############################################################
[Outputs]
[./Exodus]
type = Exodus
file_base = out
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
(test/tests/PT_TH_elder/step2_wseos.i)
#### FALCON test case: PT_TH_elder
#### Author: Yidong Xia (INL)
#### Description:
#### P-T based coupled mass & energy balance –--
#### thermally induced buoyant convection and instability
#### (using water-steam EOS)
############################################################
[Mesh]
file = out1.e
[]
############################################################
[Variables]
[./P]
initial_from_file_var = P
initial_from_file_timestep = 2
[../]
[./T]
initial_condition = 285.15
[../]
[]
############################################################
[Preconditioning]
[./SMP]
#type = FDP
type = SMP
off_diag_row = 'P'
off_diag_column = 'T'
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = FIRST
family = MONOMIAL
[../]
[./v_y]
order = FIRST
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_time_derivative]
type = PTMassTimeDerivative
variable = P
[../]
[./P_mass_residual]
type = PTMassResidual
variable = P
coupled_temperature = T
[../]
[./T_time_derivative]
type = PTEnergyTimeDerivative
variable = T
[../]
[./T_energy_residual]
type = PTEnergyResidual
variable = T
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[]
############################################################
[BCs]
#### 1 = bottom left half (20 Celsius)
#### 2 = bottom right half (12 Celsius)
#### 3 = top boundary (12 Celsius)
[./P_top]
type = DirichletBC
variable = P
boundary = '3'
value = 101325
[../]
[./T_top]
type = DirichletBC
variable = T
boundary = '3'
value = 285.15
[../]
[./T_bottom1]
type = DirichletBC
variable = T
boundary = '1'
value = 293.15
[../]
[./T_bottom2]
type = DirichletBC
variable = T
boundary = '2'
value = 285.15
[../]
[]
############################################################
[Materials]
[./porous]
type = PTGeothermal
block = 1
pressure = P
temperature = T
fluid_property_formulation = 'wseos'
stabilizer = 'none'
permeability = 1.0e-10 # default = 1.000e-12
porosity = 0.4 # default = 0.300e+00
compressibility = 1.0e-10 # default = 1.000e-05
density_rock = 2480 # default = 2.500e+03
density_water = 999.526088 # default = 1.000e+03
viscosity_water = 0.0012389 # default = 1.200e-04
specific_heat_rock = 920 # default = 0.920e+03
specific_heat_water = 4180 # default = 4.186e+03
thermal_conductivity = 1.0 # default = 2.500e+00
gravity = 9.80665 # default = 9.80665
gravity_direction = '0 -1 0' # default = '0 0 -1'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
num_steps = 1 # use 50000 for complete run
dt = 10000
end_time = 630720000
l_max_its = 30
l_tol = 1e-4 # could use 1e-3 for faster stepping
nl_max_its = 15
nl_rel_tol = 1e-4 # could use 1e-3 for faster stepping
nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
file_base = out2_wseos
# use 100 for complete run
exodus = true
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
(test/tests/PT_H_SideFluxAverage/inp.i)
############################################################
[Mesh]
type = GeneratedMesh
dim = 3
nx = 20
ny = 2
nz = 2
xmin = 0
xmax = 20
ymin = 0
ymax = 2
zmin = 0
zmax = 2
[]
############################################################
[Variables]
[./P]
initial_condition = 30e6
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./v_z]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./vz]
type = PTFluidVelocityAux
variable = v_z
component = 2
[../]
[]
############################################################
[BCs]
[./P_in]
type = PTMassFluxBC
variable = P
boundary = 'left'
value = 1
[../]
#[./P_in]
# type = DirichletBC
# variable = P
# boundary = 'left'
# value = 40e6
#[../]
[./P_out]
type = DirichletBC
variable = P
boundary = 'right'
value = 30e6
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '0'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1.0e-12
porosity = 0.1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-7
nl_max_its = 500
#nl_rel_tol = 1e-5
#nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
#exodus = true
[./Exodus]
type = Exodus
file_base = out
#use_problem_dimension = false
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
[Postprocessors]
[./SideFluxAverage]
type = PTMassSideFluxAverage
variable = P
boundary = 'right'
[../]
[]
############################################################
(test/tests/PT_H_permeability_function/test2.i)
############################################################
[Mesh]
file = mesh.e
[]
############################################################
[Variables]
[./P]
initial_condition = 20e6
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./k]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./k]
type = PTPermeabilityAux
variable = k
[../]
[]
############################################################
[BCs]
[./P_in]
type = DirichletBC
variable = P
boundary = '5'
value = 20e6
[../]
[./P_out]
type = DirichletBC
variable = P
boundary = '6'
value = 1e6
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '1'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability_function_type = '2'
permeability = 4.368e-11
porosity = 1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 0.0
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[./matrix]
type = PTGeothermal
block = '2'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1e-20
porosity = 0.01
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 3
gravity = 0.0
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'PJFNK' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 1000
l_tol = 1e-3
nl_max_its = 500
#nl_rel_tol = 1e-12
nl_abs_tol = 1e-9
[]
############################################################
[Outputs]
[./Exodus]
type = Exodus
file_base = out2
[../]
[./CSV]
type = CSV
file_base = out2_post
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
[Postprocessors]
[./InMassFlowRate]
type = PTMassSideFluxIntegral
variable = P
boundary = '5'
[../]
[./OutMassFlowRate]
type = PTMassSideFluxIntegral
variable = P
boundary = '6'
[../]
[]
(test/tests/PT_H_gravity/inp.i)
#### Test case: P-T based mass balance due to gravity
#### Author: Yidong Xia (INL)
############################################################
[Mesh]
file = mesh.e
[]
############################################################
[Variables]
[./P]
order = FIRST
family = LAGRANGE
initial_condition = 101325
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = FIRST
family = MONOMIAL
[../]
[./v_y]
order = FIRST
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_time_derivative]
type = PTMassTimeDerivative
variable = P
[../]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[]
############################################################
[BCs]
[./P_top]
type = DirichletBC
variable = P
boundary = '3'
value = 101325
[../]
[]
############################################################
[Materials]
[./porous]
type = PTGeothermal
block = 1
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1.0e-10 # default = 1.000e-12
porosity = 0.4 # default = 0.300e+00
compressibility = 4.0e-10 # default = 1.000e-05
density_rock = 2480 # default = 2.500e+03
density_water = 999.526 # default = 1.000e+03
viscosity_water = 0.001 # default = 1.200e-04
specific_heat_rock = 920 # default = 0.920e+03
specific_heat_water = 4180 # default = 4.186e+03
thermal_conductivity = 1.0 # default = 2.500e+00
gravity = 9.80665 # default = 9.80665
gravity_direction = '0 -1 0' # default = '0 0 -1'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Transient # Steady | Transient
num_steps = 5
dt = 1296000
solve_type = 'PJFNK' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-12
nl_max_its = 15
nl_rel_tol = 1e-5
nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
file_base = out
exodus = true
[./Console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
(test/tests/PT_TH_SGW2016_Reference_Case/inp3d_run.i)
#### The reference test case in the paper of Xia et al.
#### at 2016 Stanford Geothermal Workshop.
#### The PDF paper is uploaded separately.
#### This is the input file for the Production Run
############################################################
[Mesh]
file = out3d_ini.e
[]
############################################################
[Functions]
[./initial_T_distribution]
type = ParsedFunction
value = 493.15-0.065*z
[../]
[]
############################################################
[Variables]
[./P]
initial_from_file_var = P
initial_from_file_timestep = 2
[../]
[./T]
[./InitialCondition]
type = FunctionIC
function = initial_T_distribution
[../]
[../]
[]
############################################################
[ICs]
[]
############################################################
[Preconditioning]
[./SMP]
#type = FDP
type = SMP
off_diag_row = 'P'
off_diag_column = 'T'
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./v_z]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_time_derivative]
type = PTMassTimeDerivative
variable = P
[../]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[./T_time_derivative]
type = PTEnergyTimeDerivative
variable = T
[../]
[./T_energy_residual]
type = PTEnergyResidual
variable = T
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./vz]
type = PTFluidVelocityAux
variable = v_z
component = 2
[../]
[]
############################################################
[BCs]
[./P_in]
type = PTMassFluxBC
variable = P
boundary = '1001 1002 1003 1004 1005'
value = 0.5
[../]
[./P_out]
type = DirichletBC
variable = P
boundary = '1006 1007 1008 1009 1010' # T = 203.75 deg. C
value = 30e6
[../]
[./T_in]
type = DirichletBC
variable = T
value = 353.15
boundary = '1001 1002 1003 1004 1005 1111'
[../]
[./T_out]
type = PTEnergyOutFlowBC
variable = T
boundary = '1006 1007 1008 1009 1010'
[../]
[]
############################################################
[Materials]
[./matrix]
type = PTGeothermal
block = '4'
#pressure = P
temperature = T
reference_pressure = 30e6
fluid_property_formulation = 'wseos'
stabilizer = 'none'
permeability = 1e-20
porosity = 0.01
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 3
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[./fractures]
type = PTGeothermal
block = '3'
pressure = P
temperature = T
fluid_property_formulation = 'wseos'
stabilizer = 'supg'
permeability = 5e-12
porosity = 0.1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5 #100 work |
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Transient
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
#dt = 1e1
dtmin = 1
[./TimeStepper]
type = SolutionTimeAdaptiveDT
dt = 10
[../]
end_time = 1e9
num_steps = 2 # set to 1000000 for complete simulation
l_max_its = 30
l_tol = 1e-3 # could use 1e-3 for faster stepping
nl_max_its = 10
nl_rel_tol = 1e-2 # could use 1e-3 for faster stepping
nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
[./Exodus]
type = Exodus
file_base = out3d_run
#use_problem_dimension = false
interval = 100
[../]
[./CSV]
type = CSV
file_base = out3d_run
interval = 10
[../]
[./CONSOLE]
type = Console
output_linear = true
output_nonlinear = true
execute_postprocessors_on = 'none'
interval = 1
[../]
[]
############################################################
[Postprocessors]
[./DT]
type = TimestepSize
[../]
[./P_1001]
type = SideAverageValue
variable = P
boundary = '1001'
[../]
[./P_1002]
type = SideAverageValue
variable = P
boundary = '1002'
[../]
[./P_1003]
type = SideAverageValue
variable = P
boundary = '1003'
[../]
[./P_1004]
type = SideAverageValue
variable = P
boundary = '1004'
[../]
[./P_1005]
type = SideAverageValue
variable = P
boundary = '1005'
[../]
[./T_1006]
type = SideAverageValue
variable = T
boundary = '1006'
[../]
[./T_1007]
type = SideAverageValue
variable = T
boundary = '1007'
[../]
[./T_1008]
type = SideAverageValue
variable = T
boundary = '1008'
[../]
[./T_1009]
type = SideAverageValue
variable = T
boundary = '1009'
[../]
[./T_1010]
type = SideAverageValue
variable = T
boundary = '1010'
[../]
[./MassRate_1006]
type = PTMassSideFluxIntegral
variable = P
boundary = '1006'
[../]
[./MassRate_1007]
type = PTMassSideFluxIntegral
variable = P
boundary = '1007'
[../]
[./MassRate_1008]
type = PTMassSideFluxIntegral
variable = P
boundary = '1008'
[../]
[./MassRate_1009]
type = PTMassSideFluxIntegral
variable = P
boundary = '1009'
[../]
[./MassRate_1010]
type = PTMassSideFluxIntegral
variable = P
boundary = '1010'
[../]
[]
############################################################
(test/tests/PT_TH_elder/step2_comp.i)
#### FALCON test case: PT_TH_elder
#### Author: Yidong Xia (INL)
#### Description:
#### P-T based coupled mass & energy balance –--
#### thermally induced buoyant convection and instability
#### (using water-steam EOS)
############################################################
[Mesh]
file = out1.e
[]
############################################################
[Variables]
[./P]
initial_from_file_var = P
initial_from_file_timestep = 2
[../]
[./T]
initial_condition = 285.15
[../]
[]
############################################################
[Preconditioning]
[./SMP]
#type = FDP
type = SMP
off_diag_row = 'P'
off_diag_column = 'T'
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = FIRST
family = MONOMIAL
[../]
[./v_y]
order = FIRST
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_time_derivative]
type = PTMassTimeDerivative
variable = P
[../]
[./P_mass_residual]
type = PTMassResidual
variable = P
coupled_temperature = T
[../]
[./T_time_derivative]
type = PTEnergyTimeDerivative
variable = T
[../]
[./T_energy_residual]
type = PTEnergyResidual
variable = T
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[]
############################################################
[BCs]
#### 1 = bottom left half (20 Celsius)
#### 2 = bottom right half (12 Celsius)
#### 3 = top boundary (12 Celsius)
[./P_top]
type = DirichletBC
variable = P
boundary = '3'
value = 101325
[../]
[./T_top]
type = DirichletBC
variable = T
boundary = '3'
value = 285.15
[../]
[./T_bottom1]
type = DirichletBC
variable = T
boundary = '1'
value = 293.15
[../]
[./T_bottom2]
type = DirichletBC
variable = T
boundary = '2'
value = 285.15
[../]
[]
############################################################
[Materials]
[./porous]
type = PTGeothermal
block = 1
pressure = P
temperature = T
fluid_property_formulation = 'compressibility'
stabilizer = 'none'
permeability = 1.0e-10 # default = 1.000e-12
porosity = 0.4 # default = 0.300e+00
compressibility = 1.0e-10 # default = 1.000e-05
density_rock = 2480 # default = 2.500e+03
density_water = 999.526088 # default = 1.000e+03
viscosity_water = 0.0012389 # default = 1.200e-04
specific_heat_rock = 920 # default = 0.920e+03
specific_heat_water = 4180 # default = 4.186e+03
thermal_conductivity = 1.0 # default = 2.500e+00
gravity = 9.80665 # default = 9.80665
gravity_direction = '0 -1 0' # default = '0 0 -1'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Transient
solve_type = 'PJFNK'
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
num_steps = 1 # use 50000 for complete run
dt = 10000
end_time = 630720000
l_max_its = 30
l_tol = 1e-3 # could use 1e-3 for faster stepping
nl_max_its = 15
nl_rel_tol = 1e-3 # could use 1e-3 for faster stepping
nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
file_base = out2_comp
# use 100 for complete run
exodus = true
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
(test/tests/PT_TH_SGW2016_Reference_Case/inp3d_ini.i)
#### The reference test case in the paper of Xia et al.
#### at 2016 Stanford Geothermal Workshop.
#### The PDF paper is uploaded separately.
#### This is the input file for the Initial Run
############################################################
[Mesh]
file = mesh3d.e
[]
############################################################
[Variables]
[./P]
initial_condition = 30e6
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./v_z]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./vz]
type = PTFluidVelocityAux
variable = v_z
component = 2
[../]
[]
############################################################
[BCs]
[./P_in]
type = PTMassFluxBC
variable = P
boundary = '1001 1002 1003 1004 1005'
value = 0.5
[../]
[./P_out]
type = DirichletBC
variable = P
boundary = '1006 1007 1008 1009 1010'
value = 30e6
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '3'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 5e-12
porosity = 0.1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[./matrix]
type = PTGeothermal
block = '4'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1e-20
porosity = 0.01
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 3
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-3
nl_max_its = 500
#nl_rel_tol = 1e-5
nl_abs_tol = 1e-6
[]
############################################################
[Outputs]
#exodus = true
[./Exodus]
type = Exodus
file_base = out3d_ini
#use_problem_dimension = false
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
(test/tests/PT_H_SideFluxIntegral/inp.i)
############################################################
[Mesh]
type = GeneratedMesh
dim = 3
nx = 20
ny = 2
nz = 2
xmin = 0
xmax = 20
ymin = 0
ymax = 2
zmin = 0
zmax = 2
[]
############################################################
[Variables]
[./P]
initial_condition = 30e6
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./v_z]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./vz]
type = PTFluidVelocityAux
variable = v_z
component = 2
[../]
[]
############################################################
[BCs]
[./P_in]
type = PTMassFluxBC
variable = P
boundary = 'left'
value = 1
[../]
#[./P_in]
# type = DirichletBC
# variable = P
# boundary = 'left'
# value = 40e6
#[../]
[./P_out]
type = DirichletBC
variable = P
boundary = 'right'
value = 30e6
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '0'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1.0e-12
porosity = 0.1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 9.80665
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-7
nl_max_its = 500
#nl_rel_tol = 1e-5
#nl_abs_tol = 1e-7
[]
############################################################
[Outputs]
#exodus = true
[./Exodus]
type = Exodus
file_base = out
#use_problem_dimension = false
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
[Postprocessors]
[./SideFluxIntegral]
type = PTMassSideFluxIntegral
variable = P
boundary = 'right'
[../]
[]
############################################################
(test/tests/PT_H_permeability_function/test1.i)
############################################################
[Mesh]
file = mesh.e
[]
############################################################
[Variables]
[./P]
initial_condition = 0.0
[../]
[]
############################################################
[AuxVariables]
[./v_x]
order = CONSTANT
family = MONOMIAL
[../]
[./v_y]
order = CONSTANT
family = MONOMIAL
[../]
[./k]
order = CONSTANT
family = MONOMIAL
[../]
[]
############################################################
[Kernels]
[./P_mass_residual]
type = PTMassResidual
variable = P
[../]
[]
############################################################
[AuxKernels]
[./vx]
type = PTFluidVelocityAux
variable = v_x
component = 0
[../]
[./vy]
type = PTFluidVelocityAux
variable = v_y
component = 1
[../]
[./k]
type = PTPermeabilityAux
variable = k
[../]
[]
############################################################
[BCs]
[./P_in]
type = DirichletBC
variable = P
boundary = '5'
value = 10e6
[../]
[./P_out]
type = DirichletBC
variable = P
boundary = '6'
value = 0.0
[../]
[]
############################################################
[Materials]
[./fractures]
type = PTGeothermal
block = '1'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability_function_type = '1'
initial_fluid_pressure = 0.0
total_overburden_confining_pressure = 10e6 #1000e6
fitting_parameter = 1.0
permeability = 1e-12
porosity = 1
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 1.5
gravity = 0.0
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[./matrix]
type = PTGeothermal
block = '2'
pressure = P
#temperature = T
fluid_property_formulation = 'constant'
stabilizer = 'none'
permeability = 1e-20
porosity = 0.01
compressibility = 4.0e-10
density_rock = 2700
density_water = 1000
viscosity_water = 0.001
specific_heat_rock = 790
specific_heat_water = 4181
thermal_conductivity = 3
gravity = 0.0
gravity_direction = '0 0 0'
constant_pressure_gradient = '0 0 0'
[../]
[]
############################################################
[Executioner]
type = Steady
solve_type = 'NEWTON' # default = PJFNK | NEWTON
petsc_options_iname = '-pc_type -pc_hypre_type
-ksp_gmres_restart -snes_ls
-pc_hypre_boomeramg_strong_threshold'
petsc_options_value = 'hypre boomeramg 201 cubic 0.7'
l_max_its = 400
l_tol = 1e-3
nl_max_its = 500
#nl_rel_tol = 1e-5
nl_abs_tol = 1e-6
[]
############################################################
[Outputs]
[./Exodus]
type = Exodus
file_base = out1
[../]
[./CSV]
type = CSV
file_base = out1_post
[../]
[./console]
type = Console
output_linear = true
output_nonlinear = true
[../]
[]
############################################################
[Postprocessors]
[./OutMassFlowRate]
type = PTMassSideFluxIntegral
variable = P
boundary = '6'
[../]
[]