Convergence Issues #29644

julianseb · 2025-01-05T15:57:59Z

julianseb
Jan 5, 2025

Hello everybody,

im trying to implement a 2-phase flow in my model. Right now im running into some convergence issues. Later on i'm trying to implement a open water body to the right and simulate the thermal influence of burried electrical cables on groundwater and open water bodies.

I have a groundwaterflow from the left to the right that i implement over a FunctionDirichletBC with the porepressure. The gas pressure i set fixed to atmospheric pressure (101325 Pa). My gradient is 0,1 m/m and my model does have a width of 20 m. So the height difference of the groundwater from left to right border is 2 m.

For the groundwater i use the Water97FluidProperties and for the gas phase i planed to use the NitrogenFluidProperties but then the code takes forever to compute the residual, so i use the IdealGasFluidProperties. In hopes to simplify my code and better the convergence i also used the SimpleFluidProperties but they gave worse convergence than the Water97Properties. I also reduced my timestep in order to improve convergence but that didn't have a great impact.

This is my first big project and i'm kinda lost on what to do so i hope somebody could help me or point me in the right direction as to what is causing the bad convergence. I added my current state.

current_state.zip

This is my current output:

Time Step 0, time = 0

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

 0 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 3.086649e-10
 1 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 2.826211e-10
 2 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 2.955420e-10
 3 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.976411e-10
 4 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.876559e-10
 5 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 3.006491e-10
 6 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.638493e-10
...

44 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.807103e-10
45 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.695321e-10

This is my mesh:

This is my current code:

[Mesh]
  [mesh]
    type = FileMeshGenerator
    file = einfach.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'
  []
[]

[Functions]
  [gradient]
    type = ParsedFunction
    vars = 'rho g gradient left_head'
    vals = '1000 9.8065 0.1 10'
    expression = '101325 + (left_head-y-(x*gradient))*(g*rho)'
  []
[]

[BCs]
    [flow]
      type = FunctionDirichletBC
      function = gradient      
      variable = pp
      boundary = 'rechts links unten'
    []
    [flow_gp]
      type = PorousFlowOutflowBC 
      flux_type = fluid
      variable = gp
      boundary = 'rechts'
      PorousFlowDictator = dictator_boden
    []
    [top]
      type = DirichletBC
      boundary = 'oben'
      value = 101325
      variable = 'gp'
    []
[]

[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
  []
  [pressure_dt_gp]
    type = PorousFlowMassTimeDerivative
    variable = gp
    PorousFlowDictator = dictator_boden
    block = 'boden'
  []
  [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
  []
  [gp_advective]
    type = PorousFlowAdvectiveFlux
    variable = gp
    PorousFlowDictator = dictator_boden
    block = 'boden'
  []
[]

[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
  []
  [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
  []
  [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
  []
  [density_boden]
    type = PorousFlowPropertyAux
    variable = density
    property = density
    block = boden
    execute_on = 'LINEAR TIMESTEP_END'
    PorousFlowDictator = dictator_boden
  []
[]

[FluidProperties]
  [wasser]
    type = Water97FluidProperties
#    tolerance = 1e-8
#    type = SimpleFluidProperties
  []
  [dampf]
#    type = NitrogenFluidProperties
    type = IdealGasFluidProperties
  []
[]



[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'
      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 = 'pp gp'
    []
    [diff_boden]
      type = PorousFlowDiffusivityConst
      diffusion_coeff = '1E-6 1E-5 1E-6 1E-5'
      tortuosity = '0.5 0.5'
      block = 'boden'
      PorousFlowDictator = dictator_boden
    []
[]

[UserObjects]
  [dictator_boden]
    type = PorousFlowDictator
    number_fluid_components = 2
    number_fluid_phases = 2
    porous_flow_vars = 'pp temperature gp'
  []
  [capillary_pressure_boden]
    type = PorousFlowCapillaryPressureVG
    alpha = 1E-4
    m = 0.5
    block = 'boden'
#    s_scale = 0.9
  []
[]

# [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 -pc_factor_mat_solver_package'
    petsc_options_value = ' lu       mumps'
  []
[]

[Executioner]
  type = Transient
  solve_type = NEWTON
  automatic_scaling = true
  compute_scaling_once = true
  end_time = 1E8
  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
  []
[]

[Outputs]
    [exodus]
    type = Exodus
    []
    # [csv]
    #   type = CSV
    # []
[]

Answered by cpgr

Jan 8, 2025

I agree with @lynnmunday that it doesn't look big enough to run out of memory.

You could try a different preconditioned option instead of mumps. Another common one used for porous flow is

petsc_options_iname = '-pc_type -sub_pc_type -sub_pc_factor_shift_type'
petsc_options_value = ' asm      lu           NONZERO'

One more thing - your PorousFlowOutflowBC might work better if you specify the parameter

mass_fraction_component = 1

so it doesn't try and remove non-existant water in the gas phase.

View full answer

GiudGiud · 2025-01-07T00:48:26Z

GiudGiud
Jan 7, 2025
Collaborator

Hello

For the groundwater i use the Water97FluidProperties and for the gas phase i planed to use the NitrogenFluidProperties but then the code takes forever to compute the residual, so i use the IdealGasFluidProperties. In hopes to simplify my code and better the convergence i also used the SimpleFluidProperties but they gave worse convergence than the Water97Properties. I also reduced my timestep in order to improve convergence but that didn't have a great impact.

This is good feedback. I think these other properties have been implemented for correctness rather than performance. The trick here to speed it up is either do what you did, or use a TabulatedFluidProperties. This can be used to generate a tabulation (using a different fp object) then load the tabulation and only use that (no need to specify fp anymore)

With regards to convergence issues, this page should help:
https://mooseframework.inl.gov/moose/application_usage/failed_solves.html
notably it would be good to know (on a small mesh) if the problem is well posed and does have any singular value

5 replies

cpgr Jan 7, 2025
Collaborator

Hey @julianseb - have you managed to get a simpler model working (so maybe one with constant temperature or even single phase flow or even a mesh without a hole)? I always find it easiest to get something simple working and then build on that.

I can't see anything that is jumping out with your input file that is obviously wrong.

cpgr Jan 7, 2025
Collaborator

Also, interesting that the convergence is worse with Ideal gas and simple fluid properties, as they definitely have the correct derivatives! So it must be something else

julianseb Jan 7, 2025
Author

Thanks for your replies.

I've added the Debug and Outputoptions from the failed solves site. Seems like my gaspressure (gp) is causing the high residual:

Time Step 0, time = 0

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999169
                gp:          2.70791
 0 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 3.086649e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999167
                gp:          2.70791
 1 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 2.826211e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999175
                gp:          2.70791
 2 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 2.955420e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999183
                gp:          2.70791
 3 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.976411e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999185
                gp:          2.70791
 4 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.876559e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999189
                gp:          2.70791
 5 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 3.006491e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999189
                gp:          2.70791

...

49 Nonlinear |R| = 2.709754e+00
      0 Linear |R| = 2.709754e+00
      1 Linear |R| = 2.981997e-10
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999282
                gp:          2.70791
50 Nonlinear |R| = 2.709754e+00
  Nonlinear solve did not converge due to DIVERGED_MAX_IT iterations 50
 Solve Did NOT Converge!

@GiudGiud I don't care that much about the speed of my calculation. Would the TabulatedFluidProperties be less precise than the Water97FluidProperties? And for the check with a small mesh. Do you mean dimensionwise like 1 x 1 m with nothing in the mesh apart from the block?

@cpgr Yes i got the same mesh working with PorousFlow1PhaseP and a temperature at the cable.

Isn't the model i sent with constant temperature? Because i only set an IC for the Temperature so the FluidProperties could be choosen.

To the difference with the IdealGasFluidProperties to the NitrogenFluidProperties. They don't really make a difference it just takes longer. So the only difference is with the SimpleFluidProperties. The top output is with IdealGasFluidProperties and the following is with the NitrogenFluidProperties:

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    Finished Computing Residual                                                          [  6.30 s] [  186 MB]
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999172
                gp:          2.70778
 0 Nonlinear |R| = 2.709625e+00
      0 Linear |R| = 2.709625e+00
      1 Linear |R| = 3.121295e-10
    Finished Computing Residual                                                          [  6.24 s] [  204 MB]
    Finished Computing Residual                                                          [  6.17 s] [  204 MB]
    Finished Computing Residual                                                          [  6.18 s] [  204 MB]
    Finished Computing Residual                                                          [  6.17 s] [  204 MB]
    Finished Computing Residual                                                          [  6.35 s] [  204 MB]
    Finished Computing Residual                                                          [  6.20 s] [  204 MB]
    Finished Computing Residual                                                          [  6.15 s] [  204 MB]
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999171
                gp:          2.70778
 1 Nonlinear |R| = 2.709624e+00
      0 Linear |R| = 2.709624e+00
      1 Linear |R| = 2.919320e-10
    Finished Computing Residual                                                          [  6.24 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.20 s] [  206 MB]
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999179
                gp:          2.70778
 2 Nonlinear |R| = 2.709624e+00
      0 Linear |R| = 2.709624e+00
      1 Linear |R| = 2.768339e-10
    Finished Computing Residual                                                          [  6.37 s] [  206 MB]
    Finished Computing Residual                                                          [  6.18 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.16 s] [  206 MB]
    Finished Computing Residual                                                          [  6.17 s] [  206 MB]
    Finished Computing Residual                                                          [  6.16 s] [  206 MB]
    Finished Computing Residual                                                          [  6.16 s] [  206 MB]
    Finished Computing Residual                                                          [  6.16 s] [  206 MB]
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999186
                gp:          2.70778

On the Troubleshooting page is for failing nonlinear solve at the end this sentence: "If they are more than an order of magnitude off, then use the scaling parameter in the variables block to scale the smaller variable up." I already use the automatic scaling, but the values are several orders of magnitude apart. Would it be better if i do it manually or deactivate the "compute_scaling_once = true" ?

Thanks again to you

lynnmunday Jan 7, 2025

Why doesn't the nonlinear residual change between nonlinear iterations? Does that mean one of the kernels is missing a jacobian term? Maybe using the PJFNK solve type would help if it is missing a jacobian term.

julianseb Jan 8, 2025
Author

Thanks for your message.
I tested my inputfile with and without the suggestions from @cpgr and the PJFNK solver.
Without the suggestions this is the result:

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0999169
                gp:          2.70791
 0 Nonlinear |R| = 2.709755e+00
      0 Linear |R| = 2.709755e+00
      1 Linear |R| = 2.709749e+00
      2 Linear |R| = 2.709747e+00
      3 Linear |R| = 2.709746e+00
      4 Linear |R| = 2.709741e+00
      5 Linear |R| = 2.709738e+00
      6 Linear |R| = 2.709734e+00
      7 Linear |R| = 2.709734e+00

and it just keeps running. An with the suggestions it fails the same way as the Newton solver, but now it is saying, that there is a naN-issue for the temperature of the fluid:

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 0 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 

Temperature -nan 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 -nan 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_LINE_SEARCH iterations 0
 Solve Did NOT Converge!
Aborting as solve did not converge

cpgr · 2025-01-08T00:13:17Z

cpgr
Jan 8, 2025
Collaborator

I think I have found your main problem. You have

[mass_frac_boden]
  type = PorousFlowMassFraction
  PorousFlowDictator = dictator_boden
  block = 'boden'
  mass_fraction_vars = 'pp gp'
[]

The mass_fraction_vars parameter expects the mass fraction of species 1 in fluid phase 0 etc, but you are giving it the pressure variables. These are not in the range [0,1] so it is struggling to solve.

I think that is you have two fluid phases with a single fluid component in each (so only water in liquid phase, only gas in the gas phase) that you can just use mass_fraction_vars = '1 0' (so mass fraction of water in liquid is 1, mass fraction of water in gas is 0, and the others are computed from them).

10 replies

julianseb Jan 8, 2025
Author

Thanks for your message.
So i made the changes you suggested:

    [mass_frac_boden]
      type = PorousFlowMassFraction
      PorousFlowDictator = dictator_boden
      block = 'boden'
      mass_fraction_vars = '1 0'
    []

and added the fluid compontent for the Kernels:

[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]
    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
  []
[]

But now the Residual is even higher and i get an error message.

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 0 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 1 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 2 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 3 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44
 4 Nonlinear |R| = 1.218438e+03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_OUTMEMORY 
    |residual|_2 of individual variables:
                temperature: 0.000152083
                pp:          0.0917713
                gp:          1218.44

Do i understand the error message right, in that the preconditioner is running out of memory so it shuts the linear solve down?

Thanks again :)

lynnmunday Jan 8, 2025

Your problem doesn't look big enough to be running out of memory with LU, sometimes it reports this if the solve fails: #25820
But LU should work if you have under a million dofs.
Can you simplify the problem by removing the temperatre variable and kernels? Does it work with a single phase? Did you create this input file from a prexisting input file that runs. Can you try a uniform mesh with a square hole so all elements are the same size. I don't know where the problem is but creating a simpler problem will help you narrow it down.

julianseb Jan 8, 2025
Author

Thanks for your message.
Yes, i got the same mesh from the current_state.zip to work with a single phase. But there i used the SimpleFluidProperties, because initially i didn't plan on modeling two phase flow. My starting point for this file was the one that worked for the single phase and i extended it for two phases.
The Inputfile for the single phase is in the zip folder further down.

Unfortunately i can't remove the temperature variable. After removing the variable from the variable block and the dictator, the Kernels (PorousFlowEnergyTimeDerivative, PorousFlowHeatAdvection, PorousFlowHeatConduction) and the PorousFlowTemperature material i got a missing material property error. So i also excluded the PorousFlowMatrixInternalEnergy material. But still got the missing material property error. Did i miss anything involving the temperature?

For the uniform mesh. I made this one. Is it to coarse?

But then i get this error:

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000131056
                pp:          0.5283
                gp:          649.651
 0 Nonlinear |R| = 6.496510e+02
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to FACTOR_NUMERIC_ZEROPIVOT 
    |residual|_2 of individual variables:
                temperature: 0.000131056
                pp:          0.5283
                gp:          649.651
 1 Nonlinear |R| = 6.496510e+02
      0 Linear |R| = 6.496510e+02
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Error in external library
[0]PETSC ERROR: MUMPS error in solve: INFOG(1)=-3 (see users manual https://mumps-solver.org/index.php?page=doc "Error and warning diagnostics")
[0]PETSC ERROR: WARNING! There are unused option(s) set! Could be the program crashed before usage or a spelling mistake, etc!
[0]PETSC ERROR:   Option left: name:-i value: dbg_simple.i source: command line
[0]PETSC ERROR:   Option left: name:-snes_converged_reason value: ::failed source: code
[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.21.4, unknown 
[0]PETSC ERROR: ./konzept-opt on a  named f6-ing-ubuntuvm-004 by julian.munz Wed Jan  8 18:49:04 2025
[0]PETSC ERROR: Configure options --with-64-bit-indices --with-cxx-dialect=C++17 --with-debugging=no --with-fortran-bindings=0 --with-mpi=1 --with-openmp=1 --with-strict-petscerrorcode=1 --with-shared-libraries=1 --with-sowing=0 --download-fblaslapack=1 --download-hpddm=1 --download-hypre=1 --download-metis=1 --download-mumps=1 --download-ptscotch=1 --download-parmetis=1 --download-scalapack=1 --download-slepc=1 --download-strumpack=1 --download-superlu_dist=1 --with-hdf5-dir=/home/[email protected]/miniforge/envs/moose --with-make-np=16 --COPTFLAGS=-O3 --CXXOPTFLAGS=-O3 --FOPTFLAGS=-O3 --with-x=0 --with-ssl=0 --with-mpi-dir=/home/[email protected]/miniforge/envs/moose AR=/home/[email protected]/miniforge/envs/moose/bin/x86_64-conda-linux-gnu-ar RANLIB=/home/[email protected]/miniforge/envs/moose/bin/x86_64-conda-linux-gnu-ranlib CFLAGS="-march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/[email protected]/miniforge/envs/moose/include -fdebug-prefix-map=/data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work=/usr/local/src/conda/moose-petsc-3.21.4 -fdebug-prefix-map=/home/[email protected]/miniforge/envs/moose=/usr/local/src/conda-prefix" CXXFLAGS="-fvisibility-inlines-hidden -fmessage-length=0 -march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/[email protected]/miniforge/envs/moose/include " CPPFLAGS="-DNDEBUG -D_FORTIFY_SOURCE=2 -O2 -isystem /home/[email protected]/miniforge/envs/moose/include" FFLAGS="-march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/[email protected]/miniforge/envs/moose/include -fdebug-prefix-map=/data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work=/usr/local/src/conda/moose-petsc-3.21.4 -fdebug-prefix-map=/home/[email protected]/miniforge/envs/moose=/usr/local/src/conda-prefix" FCFLAGS="-march=nocona -mtune=haswell -ftree-vectorize -fPIC -fstack-protector-strong -fno-plt -O2 -ffunction-sections -pipe -isystem /home/[email protected]/miniforge/envs/moose/include -fdebug-prefix-map=/data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work=/usr/local/src/conda/moose-petsc-3.21.4 -fdebug-prefix-map=/home/[email protected]/miniforge/envs/moose=/usr/local/src/conda-prefix" LDFLAGS="-Wl,-O2 -Wl,--sort-common -Wl,--as-needed -Wl,-z,relro -Wl,-z,now -Wl,--disable-new-dtags -Wl,--gc-sections -Wl,--allow-shlib-undefined -Wl,-rpath,/home/[email protected]/miniforge/envs/moose/lib -Wl,-rpath-link,/home/[email protected]/miniforge/envs/moose/lib -L/home/[email protected]/miniforge/envs/moose/lib" --prefix=/home/[email protected]/miniforge/envs/moose/petsc
[0]PETSC ERROR: #1 MatSolve_MUMPS() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/mat/impls/aij/mpi/mumps/mumps.c:1495
[0]PETSC ERROR: #2 MatSolve() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/mat/interface/matrix.c:3633
[0]PETSC ERROR: #3 PCApply_LU() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/pc/impls/factor/lu/lu.c:169
[0]PETSC ERROR: #4 PCApply() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/pc/interface/precon.c:497
[0]PETSC ERROR: #5 PCApplyBAorAB() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/pc/interface/precon.c:767
[0]PETSC ERROR: #6 KSP_PCApplyBAorAB() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/include/petsc/private/kspimpl.h:469
[0]PETSC ERROR: #7 KSPGMRESCycle() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/ksp/impls/gmres/gmres.c:146
[0]PETSC ERROR: #8 KSPSolve_GMRES() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/ksp/impls/gmres/gmres.c:227
[0]PETSC ERROR: #9 KSPSolve_Private() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/ksp/interface/itfunc.c:905
[0]PETSC ERROR: #10 KSPSolve() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/ksp/ksp/interface/itfunc.c:1078
[0]PETSC ERROR: #11 SNESSolve_NEWTONLS() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/snes/impls/ls/ls.c:220
[0]PETSC ERROR: #12 SNESSolve() at /data/civet1/build/conda/conda-bld/moose-petsc_1723759504406/work/src/snes/interface/snes.c:4753
libMesh terminating:
MUMPS error in solve: INFOG(1)=-3 (see users manual https://mumps-solver.org/index.php?page=doc "Error and warning diagnostics")
Abort(1) on node 0 (rank 0 in comm 0): application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

I've added my current inputfile and the uniform mesh in this zip:
latest.zip

cpgr Jan 8, 2025
Collaborator

I agree with @lynnmunday that it doesn't look big enough to run out of memory.

You could try a different preconditioned option instead of mumps. Another common one used for porous flow is

petsc_options_iname = '-pc_type -sub_pc_type -sub_pc_factor_shift_type'
petsc_options_value = ' asm      lu           NONZERO'

One more thing - your PorousFlowOutflowBC might work better if you specify the parameter

mass_fraction_component = 1

so it doesn't try and remove non-existant water in the gas phase.

Answer selected by julianseb

julianseb Jan 9, 2025
Author

Thanks for the input.
That does seem to work for the beginning, but at higher timesteps it seems to fail again.

For example i let it run for the uniform mesh:

Time Step 1, time = 1, dt = 1

Performing automatic scaling calculation

    |residual|_2 of individual variables:
                temperature: 0.000131056
                pp:          0.5283
                gp:          0.000235155
 0 Nonlinear |R| = 5.283001e-01
      0 Linear |R| = 5.283001e-01
      1 Linear |R| = 4.155742e-12
    |residual|_2 of individual variables:
                temperature: 1.50261e-07
                pp:          0.00104698
                gp:          1.49071e-09
 1 Nonlinear |R| = 1.046983e-03
      0 Linear |R| = 1.046983e-03
      1 Linear |R| = 5.446954e-15
    |residual|_2 of individual variables:
                temperature: 1.89305e-13
                pp:          7.07728e-11
                gp:          1.47879e-10
 2 Nonlinear |R| = 1.639424e-10
 Solve Converged!

Outlier Variable Residual Norms:
  gp: 1.478792e-10

Time Step 2, time = 3, dt = 2
    |residual|_2 of individual variables:
                temperature: 0.000101855
                pp:          0.0430427
                gp:          1.47879e-10
 0 Nonlinear |R| = 4.304286e-02
      0 Linear |R| = 4.304286e-02
      1 Linear |R| = 3.387551e-13
    |residual|_2 of individual variables:
                temperature: 2.93435e-10
                pp:          2.90018e-06
                gp:          8.97002e-14
 1 Nonlinear |R| = 2.900179e-06
      0 Linear |R| = 2.900179e-06
      1 Linear |R| = 4.826495e-19
    |residual|_2 of individual variables:
                temperature: 7.61002e-14
                pp:          2.61994e-11
                gp:          0
 2 Nonlinear |R| = 2.619953e-11
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 2.619942e-11

Time Step 3, time = 7, dt = 4
    |residual|_2 of individual variables:
                temperature: 0.000101047
                pp:          0.0412823
                gp:          0
 0 Nonlinear |R| = 4.128239e-02
      0 Linear |R| = 4.128239e-02
      1 Linear |R| = 7.882306e-15
    |residual|_2 of individual variables:
                temperature: 5.79485e-10
                pp:          5.73685e-06
                gp:          0
 1 Nonlinear |R| = 5.736851e-06
      0 Linear |R| = 5.736851e-06
      1 Linear |R| = 5.541206e-19
    |residual|_2 of individual variables:
                temperature: 3.97652e-14
                pp:          1.05583e-11
                gp:          0
 2 Nonlinear |R| = 1.055841e-11
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 1.055834e-11

Time Step 4, time = 15, dt = 8
    |residual|_2 of individual variables:
                temperature: 0.000101025
                pp:          0.0410618
                gp:          0
 0 Nonlinear |R| = 4.106189e-02
      0 Linear |R| = 4.106189e-02
      1 Linear |R| = 6.580356e-15
    |residual|_2 of individual variables:
                temperature: 1.13279e-09
                pp:          1.12142e-05
                gp:          0
 1 Nonlinear |R| = 1.121419e-05
      0 Linear |R| = 1.121419e-05
      1 Linear |R| = 2.132151e-18
    |residual|_2 of individual variables:
                temperature: 1.90351e-14
                pp:          6.11194e-12
                gp:          0
 2 Nonlinear |R| = 6.111972e-12
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 6.111943e-12

Time Step 5, time = 31, dt = 16
    |residual|_2 of individual variables:
                temperature: 0.000101011
                pp:          0.0405702
                gp:          0
 0 Nonlinear |R| = 4.057037e-02
      0 Linear |R| = 4.057037e-02
      1 Linear |R| = 1.285182e-14
    |residual|_2 of individual variables:
                temperature: 2.16505e-09
                pp:          2.14305e-05
                gp:          0
 1 Nonlinear |R| = 2.143047e-05
      0 Linear |R| = 2.143047e-05
      1 Linear |R| = 8.095591e-18
    |residual|_2 of individual variables:
                temperature: 1.01259e-14
                pp:          9.57052e-12
                gp:          0
 2 Nonlinear |R| = 9.570522e-12
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 9.570517e-12

Time Step 6, time = 63, dt = 32
    |residual|_2 of individual variables:
                temperature: 0.000100985
                pp:          0.0396091
                gp:          0
 0 Nonlinear |R| = 3.960924e-02
      0 Linear |R| = 3.960924e-02
      1 Linear |R| = 2.460333e-14
    |residual|_2 of individual variables:
                temperature: 3.95813e-09
                pp:          3.91682e-05
                gp:          0
 1 Nonlinear |R| = 3.916816e-05
      0 Linear |R| = 3.916816e-05
      1 Linear |R| = 2.912404e-17
    |residual|_2 of individual variables:
                temperature: 7.25249e-15
                pp:          4.27596e-11
                gp:          0
 2 Nonlinear |R| = 4.275961e-11
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 4.275961e-11

...

Time Step 23, time = 2.62144e+06, dt = 1.04858e+06
    |residual|_2 of individual variables:
                temperature: 2.18385e-05
                pp:          0.000213111
                gp:          1.1951e-17
 0 Nonlinear |R| = 2.142274e-04
      0 Linear |R| = 2.142274e-04
      1 Linear |R| = 1.578726e-12
    |residual|_2 of individual variables:
                temperature: 1.96539e-05
                pp:          0.000179868
                gp:          3.45051e-12
 1 Nonlinear |R| = 1.809389e-04
      0 Linear |R| = 1.809389e-04
      1 Linear |R| = 1.407470e-12
    |residual|_2 of individual variables:
                temperature: 1.76878e-05
                pp:          0.000153433
                gp:          5.92844e-12
 2 Nonlinear |R| = 1.544496e-04
      0 Linear |R| = 1.544496e-04
      1 Linear |R| = 1.257432e-12
    |residual|_2 of individual variables:
                temperature: 1.59184e-05
                pp:          0.00013233
                gp:          7.64295e-12
 3 Nonlinear |R| = 1.332836e-04
      0 Linear |R| = 1.332836e-04
      1 Linear |R| = 1.126859e-12
    |residual|_2 of individual variables:
                temperature: 1.4326e-05
                pp:          0.000115356
                gp:          8.7629e-12
...
(still Timestep 23)
18 Nonlinear |R| = 2.258367e-05
      0 Linear |R| = 2.258367e-05
      1 Linear |R| = 1.988794e-13
    |residual|_2 of individual variables:
                temperature: 1.52874e-06
                pp:          1.65406e-05
                gp:          3.46199e-12
19 Nonlinear |R| = 1.661108e-05
      0 Linear |R| = 1.661108e-05
      1 Linear |R| = 1.287722e-13
    |residual|_2 of individual variables:
                temperature: 8.9859e-09
                pp:          1.12735e-05
                gp:          1.86524e-12
20 Nonlinear |R| = 1.127355e-05
      0 Linear |R| = 1.127355e-05
      1 Linear |R| = 3.068889e-16
    |residual|_2 of individual variables:
                temperature: 1.47373e-13
                pp:          2.4968e-10
                gp:          5.79438e-18
21 Nonlinear |R| = 2.496796e-10
 Solve Converged!

Outlier Variable Residual Norms:
  pp: 2.496796e-10

Time Step 24, time = 4.71859e+06, dt = 2.09715e+06
    |residual|_2 of individual variables:
                temperature: 1.15035e-05
                pp:          0.000132808
                gp:          5.79438e-18
 0 Nonlinear |R| = 1.333055e-04
      0 Linear |R| = 1.333055e-04
      1 Linear |R| = 1.317020e-12
    |residual|_2 of individual variables:
                temperature: 1.45014e-05
                pp:          0.000107648
                gp:          1.39441e-12
 1 Nonlinear |R| = 1.086199e-04
      0 Linear |R| = 1.086199e-04
      1 Linear |R| = 2.085017e-12
    |residual|_2 of individual variables:
                temperature: 1.3073e-05
                pp:          8.09952e-05
                gp:          2.13668e-12
 2 Nonlinear |R| = 8.204340e-05
      0 Linear |R| = 8.204340e-05
      1 Linear |R| = 2.055298e-12
    |residual|_2 of individual variables:
                temperature: 1.17854e-05
                pp:          7.55753e-05
                gp:          2.63689e-12
 3 Nonlinear |R| = 7.648866e-05
      0 Linear |R| = 7.648866e-05
      1 Linear |R| = 2.016316e-12
    |residual|_2 of individual variables:
                temperature: 1.12005e-05
                pp:          7.45407e-05
                gp:          2.64867e-12
 4 Nonlinear |R| = 7.537747e-05
      0 Linear |R| = 7.537747e-05
      1 Linear |R| = 1.946223e-12
    |residual|_2 of individual variables:
                temperature: 1.06447e-05
                pp:          7.39642e-05
                gp:          2.64562e-12
 5 Nonlinear |R| = 7.472626e-05
      0 Linear |R| = 7.472626e-05
      1 Linear |R| = 1.873821e-12
    |residual|_2 of individual variables:
                temperature: 1.05915e-05
                pp:          7.36266e-05
                gp:          2.63356e-12
...
45 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806229e-12
    |residual|_2 of individual variables:
                temperature: 1.01896e-05
                pp:          7.10504e-05
                gp:          2.5413e-12
46 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806225e-12
    |residual|_2 of individual variables:
                temperature: 1.01896e-05
                pp:          7.10504e-05
                gp:          2.5413e-12
47 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806224e-12
    |residual|_2 of individual variables:
                temperature: 1.01896e-05
                pp:          7.10504e-05
                gp:          2.5413e-12
48 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806232e-12
    |residual|_2 of individual variables:
                temperature: 1.01896e-05
                pp:          7.10504e-05
                gp:          2.5413e-12
49 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806227e-12
    |residual|_2 of individual variables:
                temperature: 1.01896e-05
                pp:          7.10504e-05
                gp:          2.5413e-12
50 Nonlinear |R| = 7.177738e-05
      0 Linear |R| = 7.177738e-05
      1 Linear |R| = 1.806223e-12
  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 50
 Solve Did NOT Converge!
Aborting as solve did not converge

Solve failed, cutting timestep.

I let it run to the end an this were the last iterations:

Time Step 40, time = 2.72429e+06, dt = 3.63798e-12
    |residual|_2 of individual variables:
                temperature: 1.06688e-05
                pp:          0.000131931
                gp:          4.00379e-18
 0 Nonlinear |R| = 1.323615e-04
      0 Linear |R| = 1.323615e-04
      1 Linear |R| = 1.891523e-20
  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 0
 Solve Did NOT Converge!
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 40, time = 2.72429e+06, dt = 1.81899e-12
    |residual|_2 of individual variables:
                temperature: 1.06688e-05
                pp:          0.000131931
                gp:          4.00379e-18
 0 Nonlinear |R| = 1.323615e-04
      0 Linear |R| = 1.323615e-04
      1 Linear |R| = 2.229316e-20
  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 0
 Solve Did NOT Converge!
Aborting as solve did not converge

Solve failed, cutting timestep.

Time Step 40, time = 2.72429e+06, dt = 1e-12
    |residual|_2 of individual variables:
                temperature: 1.06688e-05
                pp:          0.000131931
                gp:          4.00379e-18
 0 Nonlinear |R| = 1.323615e-04
      0 Linear |R| = 1.323615e-04
      1 Linear |R| = 1.677338e-20
  Nonlinear solve did not converge due to DIVERGED_LINE_SEARCH iterations 0
 Solve Did NOT Converge!
Aborting as solve did not converge


*** ERROR ***
/home/[email protected]/projects/konzept/dbg_simple.i:413.3:
The following error occurred in the TimeStepper 'IterationAdaptiveDT' of type IterationAdaptiveDT.

Solve failed and timestep already at dtmin, cannot continue!

cpgr Jan 9, 2025
Collaborator

Is that for the really simple square mesh? It might have reached a steady state and can't solve any further? Do the results look like what you expect? Does the mesh with a circle work with those changes?

julianseb Jan 9, 2025
Author

Yes, thats for the simple uniform mesh.
I would say it looks like expected. The artifacts are probably because of the coarse grid. And the temperature increase could be handled by a OutflowBC.

For the complex grid with the circle it shows a simmilar behaviour. For the first time steps it is running fine and at some point ist starts stagnating. It just shows a corious behaviour for the gas pressure:

And following on that also for the saturations of the two phases.

Would you say it would be the easiest to just adapt the IC and BC for the gaspressure? I would have thougt that the line of saturation change would be higher, closer to where i set the pp with the gradient. So going from 10 on the left to 8 on the right.

But otherwise i would say the convergence issues are handeld. Thank you !!!

lynnmunday Jan 9, 2025

@cpgr is an expert in this so definitely follow all of his suggestions before trying mine.

Anyway, I think you can get rid of temperature as a variable by creating a PorousFlowTemperature material called temperature:

[Materials]
  [temperature]
    type = PorousFlowTemperature
    temperature = T
  []
[]

You can initialize temperature with this:

[AuxVariables]
  [T]
    initial_condition = 293
    order = CONSTANT
    family = MONOMIAL
  []
[]

This should then work with [FluidProperties] type = TabulatedBicubicFluidProperties

I think your smaller mesh should work. With a small mesh you can output the jacobian and see if it looks bad:

petsc_options = '-pc_svd_monitor'
petsc_options_iname = '-pc_type'
petsc_options_value = 'svd'

cpgr · 2025-01-08T00:14:13Z

cpgr
Jan 8, 2025
Collaborator

Also just noticed that you don't specify fluid_component = 1 for every kernel for the gp variable which you should do!

0 replies

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Convergence Issues #29644

{{title}}

Replies: 3 comments 15 replies

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{editor}}'s edit

{{editor}}'s edit

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

{{title}}

Select a reply

Convergence Issues #29644

julianseb Jan 5, 2025

Replies: 3 comments · 15 replies

GiudGiud Jan 7, 2025 Collaborator

cpgr Jan 7, 2025 Collaborator

cpgr Jan 7, 2025 Collaborator

julianseb Jan 7, 2025 Author

lynnmunday Jan 7, 2025

julianseb Jan 8, 2025 Author

cpgr Jan 8, 2025 Collaborator

julianseb Jan 8, 2025 Author

lynnmunday Jan 8, 2025

julianseb Jan 8, 2025 Author

cpgr Jan 8, 2025 Collaborator

julianseb Jan 9, 2025 Author

cpgr Jan 9, 2025 Collaborator

julianseb Jan 9, 2025 Author

lynnmunday Jan 9, 2025

cpgr Jan 8, 2025 Collaborator

julianseb
Jan 5, 2025

Replies: 3 comments 15 replies

GiudGiud
Jan 7, 2025
Collaborator

cpgr Jan 7, 2025
Collaborator

cpgr Jan 7, 2025
Collaborator

julianseb Jan 7, 2025
Author

julianseb Jan 8, 2025
Author

cpgr
Jan 8, 2025
Collaborator

julianseb Jan 8, 2025
Author

julianseb Jan 8, 2025
Author

cpgr Jan 8, 2025
Collaborator

julianseb Jan 9, 2025
Author

cpgr Jan 9, 2025
Collaborator

julianseb Jan 9, 2025
Author

cpgr
Jan 8, 2025
Collaborator