Have error:"*** Warning *** Missing coupled variables {temperature} (add them to coupled_variables parameter of TM_all0) *** Warning *** Missing coupled variables {temperature} (add them to coupled_variables parameter of TM_all1)"

[Yulongggggg]

I am simulating the injection and production in geothermal systems, keep having the error:"*** Warning ***
Missing coupled variables {temperature} (add them to coupled_variables parameter of TM_all0)

*** Warning ***
Missing coupled variables {temperature} (add them to coupled_variables parameter of TM_all1)"

even though I have set like this in porousflow:

"[PorousFlowFullySaturated]
fluid_properties_type = PorousFlowSingleComponentFluid
add_darcy_aux = true
coupling_type = ThermoHydroMechanical
gravity = '0 -9.8 0'
porepressure = porepressure
temperature = temperature
fp = FluidProperties
use_displaced_mesh = False
biot_coefficient = 1
multiply_by_density = true
eigenstrain_names = 'thermal_contribution ini_strain'
block = 'cap1 cap2 block'
[]"

And I was also wondering if it's ok I setup the "eigenstrain_names = 'thermal_contribution ini_strain'" in both porous flow and solid mechanics? Will they have some discrepancies or it's ok? And i have also setted up two two boundary conditions at the injection well, one is
' [T_inlet]
type = DirichletBC
boundary = injection
variable = temperature
value = 300
[]' to make sure the injection temperature and the other

' [q_in]
type = PorousFlowSink
boundary = injection
variable = temperature
flux_function = inj_flux
fluid_phase = 0
use_enthalpy = true
[]' to simulate the injected heat with the injected water, is this ok or can i have them both?

Here's the link of my code 'https://github.com/Yulongggggg/THM/blob/main/THM.i'

[GiudGiud]

Hello

I think it's from this block

[Physics]
  [SolidMechanics]
    [QuasiStatic]
      [all]
        strain = SMALL
        incremental = false
        add_variables = true
        generate_output = 'vonmises_stress'
        material_output_family = 'MONOMIAL'
        material_output_order = 'FIRST'
        block = 'cap1 cap2 block'
        planar_formulation = 'PLANE_STRAIN'
        eigenstrain_names = 'ini_strain thermal_contribution'
      []
    []
  []
[]

maybe you need to specify the temperature there? is there a coupled_variable parameter there?

[Yulongggggg]

Thank you so much, I think that's the reason. Do you know the answer of my other questions? thank you :"And I was also wondering if it's ok I setup the "eigenstrain_names = 'thermal_contribution ini_strain'" in both porous flow and solid mechanics? Will they have some discrepancies or it's ok? And i have also setted up two two boundary conditions at the injection well, one is
' [T_inlet]
type = DirichletBC
boundary = injection
variable = temperature
value = 300
[]' to make sure the injection temperature and the other

' [q_in]
type = PorousFlowSink
boundary = injection
variable = temperature
flux_function = inj_flux
fluid_phase = 0
use_enthalpy = true
[]' to simulate the injected heat with the injected water, is this ok or can i have them both?

Here's the link of my code 'https://github.com/Yulongggggg/THM/blob/main/THM.i'"

[GiudGiud]

if it's ok I setup the "eigenstrain_names = 'thermal_contribution ini_strain'" in both porous flow and solid mechanics?

No right now you should likely delete the mechanics block, because the porous flow hydromechanics is already creating the mechanics.

Use Debug/show_actions=true and/or Problem/verbose_setup=true to verify that