Initiation Issue

[julianseb]

This is in continuation of #29644. So maybe @cpgr or @lynnmunday could help me out again. Any other help is also very appreciated.

After the other model worked out fine, i tried to implement the open water body i was talking about in the previous discussion. I set up another PorousFlowDictator (dictator_wasser) with only one component and one phase, because i see the open water body as fully saturated. I set the porosity to 0.9999.

I set up some OutflowBCs for temperature and the porepressure at the bottom of the water body (boundary: grund) in each direction with the complimentary dictators. I set the porepressure at the water surface (boundary: wasserspiegel) the same as the gas pressure in the ground at atmospheric pressure (101325 Pa).

I've managed to get a 1-phase model of it running, but there i only used one PorousFlowDictator.

But when i try to run the inputfile it initializes, tabulates the fluidproperties for the gas phase but then shows this error:

Execution Information:
  Executioner:             Transient
  TimeStepper:             IterationAdaptiveDT
  TimeIntegrator:          ImplicitEuler
  Solver Mode:             NEWTON
  PETSc Preconditioner:    asm 
  MOOSE Preconditioner:    SMP

dampf_tabulated: Generating (p, T) tabulated data


*** ERROR ***
The following error occurred in the Problem 'MOOSE Problem' of type FEProblem.

The following parallel-communicated exception was detected during INITIAL evaluation:
Pressure 0 is outside the range of tabulated pressure (12523, 3e+06).
Because this did not occur during residual evaluation, there is no way to handle this, so the solution is aborting.

So there must be a problem with the gas pressure, because the gas phase is the one with the tabulated properties. But i didn't do any changes to that in comparrison to the aformentioned case that was running.

Does anybody have an idea why the pressure apperently drops to 0 Pa? Is this because of a shared boundary between the dictators?

This is my mesh:

Here are my inputfile and the mesh:
Current_state.zip

I also added my inputfile here:

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = Gewaesser_nah_klein.msh
  []
[]

[GlobalParams]
  gravity = '0 -9.8065 0'
[]

[Variables]
  [temperature]
    initial_condition = 293.15
  []
  [pp]
  []
  [gp]
  []
[]

[ICs]
  [InitialConditionBoden]
    type = FunctionIC
    variable = pp
    function = gradient
    block = 'boden'
  []
  [ICDampf]
    type = ConstantIC
    value = 101325
    variable = gp
    block = 'boden'
  []
  [InitialConditionWasser]
    type = FunctionIC
    variable = pp
    function = wasser
    block = 'wasser'
  []
[]

[Functions]
  [gradient]
    type = ParsedFunction
    symbol_names = 'rho g gradient left_head'
    symbol_values = '1000 9.8065 0.1 10'
    expression = '101325 + (left_head-y-(x*gradient))*(g*rho)'
  []
  [t_kabel]
      type = ParsedFunction
      expression = 'if(t <= 14400, (293.15 + (363.15 - 293.15) * t / 14400), 363.15)'
  []
  [wasser]
    type = ParsedFunction
    symbol_names = 'rho g left_head'
    symbol_values = '1000 9.8065 10'
    expression = '101325 + (left_head-y)*(g*rho)'
  []
[]

[BCs]
    [flow]
      type = FunctionDirichletBC
      function = gradient      
      variable = 'pp'
      boundary = 'rechts links unten oben'
    []
    [outflow_gp]
      type = PorousFlowOutflowBC 
      flux_type = fluid
      variable = 'gp'
      boundary = 'rechts'
      PorousFlowDictator = dictator_boden
      mass_fraction_component = 1
    []
    [top]
      type = DirichletBC
      boundary = 'oben'
      value = 101325
      variable = 'gp'
    []    
    [top_wasserspiegel]
      type = FunctionDirichletBC
      variable = pp
      boundary = wasserspiegel
      function = wasser
    []
    # [temp_kabel]
    #   type = FunctionDirichletBC
    #   boundary = 'kabel'
    #   function = t_kabel
    #   variable = temperature
    # []
    [temp_out]
      type = PorousFlowOutflowBC
      PorousFlowDictator = dictator_boden
      boundary = 'rechts'
      variable = temperature
      flux_type = heat
    []
    [temp_in]
      type = DirichletBC
      boundary = links
      value = 293.15
      variable = temperature
    []
    [temp_unten]
      type = ADConvectiveHeatFluxBC
      T_infinity = 293.15
      boundary = 'unten rechts'
      heat_transfer_coefficient = 20
      variable = temperature
    []
    [temp_atm]
      type = ADConvectiveHeatFluxBC
      T_infinity = 293.15
      boundary = 'oben wasserspiegel'
      heat_transfer_coefficient = 5
      variable = temperature
    []
    [outflow_pp_wasser]
      type = PorousFlowOutflowBC
      variable = pp
      flux_type = FLUID
      boundary = 'grund'
      PorousFlowDictator = dictator_wasser
    []
    [outflow_pp_boden]
      type = PorousFlowOutflowBC
      variable = pp
      flux_type = FLUID
      boundary = 'grund'
      PorousFlowDictator = dictator_boden
    []

    [t_grund_boden]
      type = PorousFlowOutflowBC
      PorousFlowDictator = dictator_boden
      boundary = grund
      variable = temperature
      flux_type = HEAT
    []
    [t_grund_wasser]
      type = PorousFlowOutflowBC
      variable = temperature
      boundary = grund
      PorousFlowDictator = dictator_wasser
      flux_type = HEAT
    []
[]

[Kernels]
  [temperature_dt_boden]
    type = PorousFlowEnergyTimeDerivative
    variable = temperature
    PorousFlowDictator = dictator_boden
    block = 'boden'
  []
  [advectionT_boden]
    type = PorousFlowHeatAdvection
    variable = temperature
    block = 'boden'
    PorousFlowDictator = dictator_boden
  []
  [conductionT_boden]
    type = PorousFlowHeatConduction
    variable = temperature
    block = 'boden'
    PorousFlowDictator = dictator_boden
  []
  [pressure_dt_pp_boden]
    type = PorousFlowMassTimeDerivative
    variable = pp
    block = 'boden'
    PorousFlowDictator = dictator_boden
    fluid_component = 0
  []
  [pressure_dt_gp]
    type = PorousFlowMassTimeDerivative
    variable = gp
    PorousFlowDictator = dictator_boden
    block = 'boden'
    fluid_component = 1
  []
  [pp_disp_boden]
    type = PorousFlowDispersiveFlux
    disp_long = '1E-4 1E-4'
    disp_trans = '1E-5 1E-5'
    fluid_component = 0
    block = 'boden'
    variable = pp
    PorousFlowDictator = dictator_boden
  []
  [gp_disp_boden]
    type = PorousFlowDispersiveFlux
    disp_long = '5E-4 5E-4'
    disp_trans = '5E-5 5E-5'
    fluid_component = 1
    block = 'boden'
    variable = gp
    PorousFlowDictator = dictator_boden
  []
  [pp_advective_boden]
    type = PorousFlowAdvectiveFlux
    variable = pp
    block = 'boden'
    PorousFlowDictator = dictator_boden
    fluid_component = 0
  []
  [gp_advective]
    type = PorousFlowAdvectiveFlux
    variable = gp
    PorousFlowDictator = dictator_boden
    block = 'boden'
    fluid_component = 1
  []
  [temperature_dt_wasser]
    type = PorousFlowEnergyTimeDerivative
    variable = temperature
    PorousFlowDictator = dictator_wasser
    block = 'wasser'
  []
  [advectionT_wasser]
    type = PorousFlowHeatAdvection
    variable = temperature
    block = 'wasser'
    PorousFlowDictator = dictator_wasser
  []
  [conductionT_wasser]
    type = PorousFlowHeatConduction
    variable = temperature
    block = 'wasser'
    PorousFlowDictator = dictator_wasser
  []
  [pressure_dt_pp_wasser]
    type = PorousFlowMassTimeDerivative
    variable = pp
    block = 'wasser'
    PorousFlowDictator = dictator_wasser
    fluid_component = 0
  []
  [pp_disp_wasser]
    type = PorousFlowDispersiveFlux
    disp_long = '1E-6'
    disp_trans = '1E-7'
    block = 'wasser'
    variable = pp
    fluid_component = 0
    PorousFlowDictator = dictator_wasser
  []
  [pp_advective_wasser]
    type = PorousFlowAdvectiveFlux
    variable = pp
    block = 'wasser'
    fluid_component = 0
    PorousFlowDictator = dictator_wasser
  []
[]

[AuxVariables]
  [darcy_x_p0]
    family = MONOMIAL
    order = CONSTANT
  []
  [darcy_y_p0]
    family = MONOMIAL
    order = CONSTANT
  []
  [darcy_x_p1]
    family = MONOMIAL
    order = CONSTANT
  []
  [darcy_y_p1]
    family = MONOMIAL
    order = CONSTANT
  []
  [sat_p0]
    family = MONOMIAL
    order = CONSTANT
  []
  [sat_p1]
    family = MONOMIAL
    order = CONSTANT
  []
  [cp]
    family = MONOMIAL
    order = CONSTANT
  []
  [vis]
    family = MONOMIAL
    order = CONSTANT
  []
  [density]
    family = MONOMIAL
    order = CONSTANT
  []
[]

[AuxKernels]
  [d_v_x_boden_p0]
    type = PorousFlowDarcyVelocityComponent
    component = x
    variable = darcy_x_p0
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    fluid_phase = 0
  []
  [d_v_y_boden_p0]
    type = PorousFlowDarcyVelocityComponent
    component = y
    variable = darcy_y_p0
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    fluid_phase = 0
  []
  [d_v_x_boden_p1]
    type = PorousFlowDarcyVelocityComponent
    component = x
    variable = darcy_x_p1
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    fluid_phase = 1
  []
  [d_v_y_boden_p1]
    type = PorousFlowDarcyVelocityComponent
    component = y
    variable = darcy_y_p1
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    fluid_phase = 1
  []
  [d_v_x_wasser_p0]
    type = PorousFlowDarcyVelocityComponent
    component = x
    variable = darcy_x_p0
    block = 'wasser'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_wasser
    fluid_phase = 0
  []
  [d_v_y_wasser_p0]
    type = PorousFlowDarcyVelocityComponent
    component = y
    variable = darcy_y_p0
    block = 'wasser'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_wasser
    fluid_phase = 0
  []
  [saturation_boden_p0]
    type = PorousFlowPropertyAux
    variable = sat_p0
    property = saturation
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    phase = 0
  []
  [saturation_boden_p1]
    type = PorousFlowPropertyAux
    variable = sat_p1
    property = saturation
    block = 'boden'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    phase = 1
  []
  [saturation_wasser]
    type = PorousFlowPropertyAux
    variable = sat_p0
    property = saturation
    block = 'wasser'
    execute_on = 'INITIAL TIMESTEP_END'
    PorousFlowDictator = dictator_wasser
    phase = 0
  []
  [capillary_pressure]
    type = PorousFlowPropertyAux
    property = capillary_pressure
    variable = cp
    block = 'boden'
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_boden
    liquid_phase = 0
  []
  [viscosity_boden]
    type = PorousFlowPropertyAux
    variable = vis
    property = viscosity
    block = boden
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_boden
  []
  [viscosity_wasser]
    type = PorousFlowPropertyAux
    variable = vis
    property = viscosity
    block = wasser
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_wasser
  []
  [density_boden]
    type = PorousFlowPropertyAux
    variable = density
    property = density
    block = boden
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_boden
  []
  [density_wasser]
    type = PorousFlowPropertyAux
    variable = density
    property = density
    block = wasser
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_wasser
  []
[]

[FluidProperties]
  [wasser]
    type = Water97FluidProperties
  []
  [dampf]
    type = NitrogenFluidProperties
  []
  [dampf_tabulated]
    type = TabulatedBicubicFluidProperties
    fp = dampf
    interpolated_properties = 'density enthalpy internal_energy viscosity k c cv cp'
    construct_pT_from_ve = false
    construct_pT_from_vh = false
    temperature_min = 274
    temperature_max = 400
    pressure_min = 12523
    pressure_max = 3e6
    tolerance = 1e-8
    T_initial_guess = 293.15
    p_initial_guess = 1e5
  []
[]



[Materials]
    [internal_energy_boden]
      type = PorousFlowMatrixInternalEnergy
      density = 2600
      specific_heat_capacity = 1000
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [thermal_conductivity_boden]
      type = PorousFlowThermalConductivityIdeal
      dry_thermal_conductivity = '0.4 0 0  0 0.4 0 0 0 0.4'
      wet_thermal_conductivity = '1 0 0    0 1 0   0 0 1'
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [porosity_boden]
      type = PorousFlowPorosityConst
      porosity = 0.3
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [permeablity_boden]
      type = PorousFlowPermeabilityConst
      block = 'boden'
      permeability = '1E-12 0 0    0 1E-12 0    0 0 1E-12'
      PorousFlowDictator = dictator_boden
    []
    [fluid_boden]
      type = PorousFlowSingleComponentFluid
      fp = 'wasser'
      phase = 0
      PorousFlowDictator = dictator_boden
      block = 'boden'
    []
    [dampf]
      type = PorousFlowSingleComponentFluid
      fp = 'dampf_tabulated'
      phase = 1
      PorousFlowDictator = dictator_boden
      block = 'boden'
    []
    [temperature_qp_boden]
      type = PorousFlowTemperature
      temperature = temperature
      PorousFlowDictator = dictator_boden
      block = 'boden'
    []
    [saturation_calculator_boden]
      type = PorousFlow2PhasePP
      PorousFlowDictator = dictator_boden
      phase0_porepressure = pp
      phase1_porepressure = gp
      capillary_pressure = capillary_pressure_boden
      block = 'boden'
    []
    [relperm_boden_dampf]
      type = PorousFlowRelativePermeabilityCorey
      n = 2
      phase = 1
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [relperm_boden_fluid]
      type = PorousFlowRelativePermeabilityCorey
      n = 2
      s_res = 0.1
      sum_s_res = 0.1
      phase = 0
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [mass_frac_boden]
      type = PorousFlowMassFraction
      PorousFlowDictator = dictator_boden
      block = 'boden'
      mass_fraction_vars = '1 0'
    []
    [diff_boden]
      type = PorousFlowDiffusivityConst
      diffusion_coeff = '1E-6 1E-5 1E-6 1E-5'
      tortuosity = '0.5 0.5'
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
    [internal_energy_wasser]
      type = PorousFlowMatrixInternalEnergy
      density = 1000
      specific_heat_capacity = 4200
      block = 'wasser'
      PorousFlowDictator = dictator_wasser
    []
    [thermal_conductivity_wasser]
      type = PorousFlowThermalConductivityIdeal
      dry_thermal_conductivity = '0.6 0 0    0 0.6 0   0 0 0.6'
      wet_thermal_conductivity = '0.6 0 0    0 0.6 0   0 0 0.6'
      block = 'wasser'
      PorousFlowDictator = dictator_wasser
    []
    [porosity_wasser]
      type = PorousFlowPorosityConst
      porosity = 0.9999
      block = 'wasser'
      PorousFlowDictator = dictator_wasser
    []
    [permeablity_wasser]
      type = PorousFlowPermeabilityConst
      block = 'wasser'
      permeability = '1E-4 0 0    0 1E-4 0    0 0 1E-4'
      PorousFlowDictator = dictator_wasser
    []
    [fluid_wasser]
      type = PorousFlowSingleComponentFluid
      fp = 'wasser'
      phase = 0
      PorousFlowDictator = dictator_wasser
      block = 'wasser'
    []
    [temperature_qp_wasser]
      type = PorousFlowTemperature
      temperature = temperature
      PorousFlowDictator = dictator_wasser
      block = 'wasser'
    []
    [saturation_calculator_wasser]
      type = PorousFlow1PhaseP
      capillary_pressure = capillary_pressure_wasser
      PorousFlowDictator = dictator_wasser
      porepressure = pp
      block = 'wasser'
    []
    [relperm_wasser]
      type = PorousFlowRelativePermeabilityConst
      PorousFlowDictator = dictator_wasser
      phase = 0
      kr = 1
      block = 'wasser'
    []
    [mass_frac_wasser]
      type = PorousFlowMassFraction
      PorousFlowDictator = dictator_wasser
      block = 'wasser'
    []
    [diff_wasser]
      type = PorousFlowDiffusivityConst
      diffusion_coeff = '1.4E-7'
      tortuosity = '1'
      block = 'wasser'
      PorousFlowDictator = dictator_wasser
    []
[]

[UserObjects]
  [dictator_boden]
    type = PorousFlowDictator
    number_fluid_components = 2
    number_fluid_phases = 2
    porous_flow_vars = 'pp temperature gp'
  []
  [dictator_wasser]
    type = PorousFlowDictator
    number_fluid_components = 1
    number_fluid_phases = 1
    porous_flow_vars = 'pp temperature'
  []
  [capillary_pressure_boden]
    type = PorousFlowCapillaryPressureVG
    alpha = 1E-4
    m = 0.5
    block = 'boden'
  []
  [capillary_pressure_wasser]
    type = PorousFlowCapillaryPressureConst
    block = 'wasser'
    pc = 0
  []
[]

# [Postprocessors]
#   [heat]
#     type = PorousFlowHeatEnergy
#     outputs = csv
#     PorousFlowDictator = dictator_boden
#   []
#   [Fluid]
#     type = PorousFlowFluidMass
#     PorousFlowDictator = dictator_boden
#     outputs = csv
#     fluid_component = 0
#   []
#   [Gas]
#     type = PorousFlowFluidMass
#     PorousFlowDictator = dictator_boden
#     outputs = csv
#     fluid_component = 1
#   []
# []


[Preconditioning]
  [basic]
    type = SMP
    full = true
    petsc_options_iname = '-pc_type -sub_pc_type -sub_pc_factor_shift_type'
    petsc_options_value = ' asm      lu           NONZERO'
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  automatic_scaling = true
  compute_scaling_once = true
  end_time = 1E9
  nl_abs_tol = 1E-6
  nl_rel_tol = 1E-6
  nl_max_its = 20
  
  [TimeStepper]
    type = IterationAdaptiveDT
    dt = 1
    growth_factor = 2
    cutback_factor = 0.5
  []
[]

[Debug]
  show_var_residual_norms = true
[]

[Outputs]
  exodus = true
  print_linear_residuals = true
  # [csv]
  #   type = CSV
  # []
[]

[lynnmunday]

I've never seen a simulation with more than one PorousFlowDictator objects. I didn't think you could do that. I looked through the test suite and I didn't find an example that used more than one. Did you find any tests that do something similar with only one PorousFlowDictator? Can you use one?

[julianseb]

No i did not find tests doing that. But i thought since the PorousFlowDictator is a UserObject that is called upon by the Kernels, Materials, etc. i could use a arbitrary number of Dictators as long as they are used by different Objects. At least it was worth a try.

I guess i will have to use one PorousFlowDictator then or maybe look at an multiapp system (which seems a bit overwhelming at the moment).

But anyway. Thanks a lot :)

[lynnmunday]

@WilkAndy or @cpgr can multiple porousFlowDictator's be used in a single input file?

[julianseb]

So i rewrote my inputfile so it only contains one PorousFlowDictator.

If i don't have any BC on the border between the ground and the open water body it runs as expected like #29644 but with a no flow boundary in its way.

But i want to have a coupling of the to blocks so that heat and porepressure are exchanged. I tried to implement this with the OutflowBCs. For only the porepressure it seems fine. When i only engage the temperature it also runs, but i get artifacts on the border from the waterbody to the unsaturated zone and there seem to be some cooling i can not explain.

Sometimes i get an error message like this that would explain that there is cooling but unfortunately not why it is happening:

Temperature 270.67 is out of range in wasser: inRegion()
To recover, the solution will fail and then be re-attempted with a reduced time step.

A MooseException was raised during FEProblemBase::computeResidualTags
Temperature 270.67 is out of range in wasser: inRegion()
To recover, the solution will fail and then be re-attempted with a reduced time step.

  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 8
 Solve Did NOT Converge!

And when i engage the temperature and the porepressure OutflowBCs at the same time the convergence is really bad and it cuts the timestep constantly.

...
Time Step 1, time = 0.03125, dt = 0.03125
    |residual|_2 of individual variables:
                temperature: 0.0676454
                pp:          0.00275448
                gp:          1.55676
 0 Nonlinear |R| = 1.558232e+00
      0 Linear |R| = 1.558232e+00
      1 Linear |R| = 3.733900e-09
    |residual|_2 of individual variables:
                temperature: 0.180709
                pp:          0.00227147
                gp:          0.000387099
 1 Nonlinear |R| = 1.807236e-01
      0 Linear |R| = 1.807236e-01
      1 Linear |R| = 2.343165e-10
    |residual|_2 of individual variables:
                temperature: 0.0450336
                pp:          0.000574951
                gp:          1.89311e-09
 2 Nonlinear |R| = 4.503732e-02
      0 Linear |R| = 4.503732e-02
      1 Linear |R| = 1.773325e-10
    |residual|_2 of individual variables:
                temperature: 0.0236775
                pp:          0.000298875
                gp:          7.56092e-10
 3 Nonlinear |R| = 2.367936e-02
      0 Linear |R| = 2.367936e-02
      1 Linear |R| = 2.251880e-10
    |residual|_2 of individual variables:
                temperature: 0.0171082
                pp:          0.000216034
                gp:          5.8726e-10
 4 Nonlinear |R| = 1.710954e-02
      0 Linear |R| = 1.710954e-02
      1 Linear |R| = 5.675462e-11
    |residual|_2 of individual variables:
                temperature: 0.0142076
                pp:          0.000178531
                gp:          3.70942e-10
 5 Nonlinear |R| = 1.420868e-02
      0 Linear |R| = 1.420868e-02
      1 Linear |R| = 3.797841e-11
    |residual|_2 of individual variables:
                temperature: 0.0106693
                pp:          0.00013447
                gp:          2.0394e-10
 6 Nonlinear |R| = 1.067016e-02
      0 Linear |R| = 1.067016e-02
      1 Linear |R| = 1.301987e-10
    |residual|_2 of individual variables:
                temperature: 0.00889562
                pp:          0.00011211
                gp:          2.42013e-10
 7 Nonlinear |R| = 8.896326e-03
      0 Linear |R| = 8.896326e-03
      1 Linear |R| = 6.176494e-11
    |residual|_2 of individual variables:
                temperature: 0.00763035
                pp:          9.6169e-05
                gp:          2.65802e-10
 8 Nonlinear |R| = 7.630958e-03
      0 Linear |R| = 7.630958e-03
      1 Linear |R| = 4.222395e-11
    |residual|_2 of individual variables:
                temperature: 0.00747023
                pp:          9.40895e-05
                gp:          2.80426e-10
 9 Nonlinear |R| = 7.470826e-03
      0 Linear |R| = 7.470826e-03
      1 Linear |R| = 2.566226e-11
    |residual|_2 of individual variables:
                temperature: 0.00707274
                pp:          8.90493e-05
                gp:          2.89188e-10
10 Nonlinear |R| = 7.073304e-03
  Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 10
  Finished Solving                                                                       [103.61 s] [  343 MB]
 Solve Did NOT Converge!
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 0.015625, dt = 0.015625
    |residual|_2 of individual variables:
                temperature: 0.0676454
                pp:          0.00275448
                gp:          1.55676
 0 Nonlinear |R| = 1.558232e+00
      0 Linear |R| = 1.558232e+00
      1 Linear |R| = 3.131965e-09
    |residual|_2 of individual variables:
                temperature: 0.21916
                pp:          0.00273185
                gp:          0.000386187
 1 Nonlinear |R| = 2.191774e-01
      0 Linear |R| = 2.191774e-01
      1 Linear |R| = 2.852329e-10
    |residual|_2 of individual variables:
                temperature: 0.0555337
                pp:          0.00070943
                gp:          4.09194e-09
 2 Nonlinear |R| = 5.553823e-02
      0 Linear |R| = 5.553823e-02
      1 Linear |R| = 3.365325e-10
    |residual|_2 of individual variables:
                temperature: 0.0301049
                pp:          0.000374472
                gp:          5.37534e-09
 3 Nonlinear |R| = 3.010720e-02
      0 Linear |R| = 3.010720e-02
      1 Linear |R| = 2.476050e-10
    |residual|_2 of individual variables:
                temperature: 0.0252734
                pp:          0.000314313
                gp:          4.5801e-09
 4 Nonlinear |R| = 2.527531e-02
      0 Linear |R| = 2.527531e-02
      1 Linear |R| = 1.504118e-10
    |residual|_2 of individual variables:
                temperature: 0.0202637
                pp:          0.00025191
                gp:          3.51368e-09
 5 Nonlinear |R| = 2.026522e-02
      0 Linear |R| = 2.026522e-02
      1 Linear |R| = 5.070768e-11
    |residual|_2 of individual variables:
                temperature: 0.0169667
                pp:          0.000212282
                gp:          2.28632e-10
 6 Nonlinear |R| = 1.696800e-02
      0 Linear |R| = 1.696800e-02
      1 Linear |R| = 1.885023e-10
    |residual|_2 of individual variables:
                temperature: 0.0151194
                pp:          0.000189108
                gp:          2.69138e-10
 7 Nonlinear |R| = 1.512060e-02
      0 Linear |R| = 1.512060e-02
      1 Linear |R| = 1.330392e-10
    |residual|_2 of individual variables:
                temperature: 0.0135125
                pp:          0.000168974
                gp:          2.85915e-10
 8 Nonlinear |R| = 1.351356e-02
      0 Linear |R| = 1.351356e-02
      1 Linear |R| = 1.037414e-10
    |residual|_2 of individual variables:
                temperature: 0.0121035
                pp:          0.000151336
                gp:          3.11829e-10
 9 Nonlinear |R| = 1.210447e-02
      0 Linear |R| = 1.210447e-02
      1 Linear |R| = 8.473254e-11
    |residual|_2 of individual variables:
                temperature: 0.011714
                pp:          0.000146466
                gp:          3.34518e-10
10 Nonlinear |R| = 1.171495e-02
  Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 10
 Solve Did NOT Converge!
  Finished Solving                                                                       [103.06 s] [  343 MB]
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 1, time = 0.0078125, dt = 0.0078125
    |residual|_2 of individual variables:
                temperature: 0.0676454
                pp:          0.00275448
                gp:          1.55676
 0 Nonlinear |R| = 1.558232e+00
...

Also adding the gaspressure would be nice, but isn't really necessary.

Is ther a better way to couple two blocks like there is no boundary than with OutflowBCs?

This is my current BC block. I also added my current inputfile and mesh:
Current_state.zip

[BCs]
    [flow]
      type = FunctionDirichletBC
      function = gradient      
      variable = 'pp'
      boundary = 'rechts links unten oben'
    []
    [outflow_gp]
      type = PorousFlowOutflowBC 
      flux_type = fluid
      variable = 'gp'
      boundary = 'rechts'
      PorousFlowDictator = dictator_boden
      mass_fraction_component = 1
    []
    [top]
      type = DirichletBC
      boundary = 'oben wasserspiegel'
      value = 101325
      variable = 'gp'
    []
    [top_wasserspiegel]
      type = FunctionDirichletBC
      variable = pp
      boundary = wasserspiegel
      function = wasser
    []
    # [temp_kabel]
    #   type = FunctionDirichletBC
    #   boundary = 'kabel'
    #   function = t_kabel
    #   variable = temperature
    # []
    [temp_out]
      type = PorousFlowOutflowBC
      PorousFlowDictator = dictator_boden
      boundary = 'rechts oben'
      variable = temperature
      flux_type = heat
    []
    [temp_in]
      type = DirichletBC
      boundary = links
      value = 293.15
      variable = temperature
    []
    [temp_unten]
      type = ADConvectiveHeatFluxBC
      T_infinity = 293.15
      boundary = 'unten rechts'
      heat_transfer_coefficient = 20
      variable = temperature
    []
    [temp_atm]
      type = DirichletBC
      boundary = 'oben wasserspiegel'
      variable = temperature
      value = 293.15
    []
    [pp_grund]
      type = PorousFlowOutflowBC
      variable = pp
      flux_type = FLUID
      boundary = 'grund'
      PorousFlowDictator = dictator_boden
      mass_fraction_component = 0
    []
    [t_grund]
      type = PorousFlowOutflowBC
      PorousFlowDictator = dictator_boden
      boundary = grund
      variable = temperature
      flux_type = HEAT
    []
    # [gp_grund]
    #   type = PorousFlowOutflowBC
    #   variable = gp
    #   flux_type = FLUID
    #   boundary = 'grund'
    #   PorousFlowDictator = dictator_boden
    #   mass_fraction_component = 1
    # []
[]