SUBPC_ERROR while using multiphase model #28937

ashishdhole · 2024-10-25T23:24:11Z

ashishdhole
Oct 25, 2024

Hello,

I am working on phase transformation in steels. I am using multiphase model at the moment because:
a. I will be using thermodynamics database for free energies of each phase and will have multiple grains per phase. Usinf lagrange's multiplier It will be easy for me to keep the sum of order parameter = 1 at every point. (I tried KKS system, but fr some reason Lagrange multiplier was not working)
b. I have to introduce elastic field later in my model.
c. incorporation of discrete nucleation model.

With all these requirements I thought multiphase model is good. I made the following input file

[Mesh]
    type = GeneratedMesh
    dim = 2
    nx = 50
    ny = 50
    xmin = 0
    xmax = 50
    ymin = 0
    ymax = 50
    elem_type = QUAD4
[]

[GlobalParams]
    op_num = 2
    var_name_base = eta
[]

[Variables]
    # carbon concentration
[c]
    order = FIRST
    family = LAGRANGE
[InitialCondition]
    type = ConstantIC
    value = 0.1
[]
[]
[w]
    order = FIRST
    family = LAGRANGE
[]
[eta0]
    order = FIRST
    family = LAGRANGE
[InitialCondition]
    type = SmoothCircleIC
    x1 = 25
    y1 = 25
    radius = 10
    invalue = 0.9
    outvalue = 0.1
    int_width = 1
[]
[]
[eta1]
    order = FIRST
    family = LAGRANGE
[InitialCondition]
    type = SmoothCircleIC
    x1 = 25
    y1 = 25
    radius = 10
    invalue = 0.1
    outvalue = 0.9
    int_width = 1
[]
[]
[PolycrystalVariables]
    order = FIRST
    family = LAGRANGE
[]

# Lagrange multiplier
[lambda]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.0
[]
[]

[AuxVariables]
[bnds]
[]
[]
[AuxKernels]
[BndsCalc]
    type = BndsCalcAux
    variable = bnds
    execute_on = 'timestep_end'
[]
[]

[BCs]
[Periodic]
[all]
    auto_direction = ''
[]
[]
[]

[Materials]
[fb]
    type = DerivativeParsedMaterial
    property_name = fb
    expression = 'c_fe:=1-c_mn-c; G := 8.3145*T*(1.0*if(c_fe > 1.0e-15,c_fe*log(c_fe),0) + 1.0*if(c_mn > 1.0e-15,c_mn*log(c_mn),0) + 1.0*if(c > 1.0e-15,c*log(c),0) + 1.0*if(c_va > 1.0e-15,c_va*log(c_va),0))/(1.0*c_fe + 1.0*c_mn + 1.0*c) + 8.3145*T*if(T < 95.0*c_mn*c_va,-3.76634024154742e-33*T^15/(c_mn*c_va + 1.05263157894737e-11)^15 - 1.23049163062049e-20*T^9/(c_mn*c_va + 1.05263157894737e-11)^9 - 2.03517994308988e-7*T^3/(c_mn*c_va + 1.05263157894737e-11)^3 + 1 - 0.86033875460538*(95.0*c_mn*c_va + 1.0e-9)/T,if(T < -285.0*c_mn*c_va,-2.62482727189424e-40*T^15/(-c_mn*c_va + 3.50877192982456e-12)^15 - 6.25154514362899e-25*T^9/(-c_mn*c_va + 3.50877192982456e-12)^9 - 7.53770349292548e-9*T^3/(-c_mn*c_va + 3.50877192982456e-12)^3 + 1 - 0.86033875460538*(-285.0*c_mn*c_va + 1.0e-9)/T,if(T > -285.0*c_mn*c_va & 285.0*c_mn*c_va < 0,-80269907448.7061*(-c_mn*c_va + 3.50877192982456e-12)^5/T^5 - 9.00929501779248e+33*(-c_mn*c_va + 3.50877192982456e-12)^15/T^15 - 6.6889616090583e+57*(-c_mn*c_va + 3.50877192982456e-12)^25/T^25,if(T > 95.0*c_mn*c_va & 285.0*c_mn*c_va > 0,-330328837.237474*(c_mn*c_va + 1.05263157894737e-11)^5/T^5 - 6.27873260157897e+26*(c_mn*c_va + 1.05263157894737e-11)^15/T^15 - 7.8945491943478e+45*(c_mn*c_va + 1.05263157894737e-11)^25/T^25,0))))*log(-0.66*c_mn*c_va*if(0.66*c_mn*c_va >= 0,-0.333333333333333,1.0) + 1)/(1.0*c_fe + 1.0*c_mn + 1.0*c) + (-36732.0*c_fe*c_mn*c - 10184.0*c_fe*c_mn*c_va - 34671.0*c_fe*c*c_va - 52204.0*c_mn*c*c_va)/(1.0*c_fe + 1.0*c_mn + 1.0*c) + (c_fe*c*(if(T < 103.0,-0.09009204*T - 2.75e-5*T^3.0 - 1049.140845,if(T < 350.0,1.76583*T*log(T) - 7.39898691*T - 0.01706952*T^2.0 - 988.25091,12000000000.0*T^(-3.0) - 264300000.0*T^(-2.0) + 2562600.0*1/T - 24.3*T*log(T) + 170.73*T - 0.0004723*T^2.0 - 17368.441)) + if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) + 80000.0) + c_fe*c_va*(if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) + 4745.0) + c_mn*c*(27.46525*T + if(T < 103.0,-0.09009204*T - 2.75e-5*T^3.0 - 1049.140845,if(T < 350.0,1.76583*T*log(T) - 7.39898691*T - 0.01706952*T^2.0 - 988.25091,12000000000.0*T^(-3.0) - 264300000.0*T^(-2.0) + 2562600.0*1/T - 24.3*T*log(T) + 170.73*T - 0.0004723*T^2.0 - 17368.441)) + if(T < 36.7,-1.0947477*T - 4.25e-6*T^4.0 - 4901.51892,if(T < 165.0,5.82736*T*log(T) - 20.019344*T - 0.12386*T^2.0 + 8.75e-5*T^3.0 - 4822.7003,if(T < 1519.0,69827.0*1/T - 23.4582*T*log(T) + 130.059*T - 0.00734768*T^2.0 - 8115.28,1.656847e+30*T^(-9.0) - 48.0*T*log(T) + 312.2648*T - 28733.41)))) + c_mn*c_va*if(T < 36.7,-1.0947477*T - 4.25e-6*T^4.0 - 4901.51892,if(T < 165.0,5.82736*T*log(T) - 20.019344*T - 0.12386*T^2.0 + 8.75e-5*T^3.0 - 4822.7003,if(T < 1519.0,69827.0*1/T - 23.4582*T*log(T) + 130.059*T - 0.00734768*T^2.0 - 8115.28,1.656847e+30*T^(-9.0) - 48.0*T*log(T) + 312.2648*T - 28733.41))))/(1.0*c_fe + 1.0*c_mn + 1.0*c); G/100000'
    coupled_variables = 'c'
    constant_names = 'T c_mn c_va'
    constant_expressions = '900 0.0001 0.01'
[]
[ff]
    type = DerivativeParsedMaterial
    property_name = ff
    expression = 'c_fe:=1-c_mn-c; G := 8.3145*T*(1.0*if(c_fe > 1.0e-15,c_fe*log(c_fe),0) + 1.0*if(c_mn > 1.0e-15,c_mn*log(c_mn),0) + 1.0*if(c > 1.0e-15,c*log(c),0) + 1.0*if(c_va > 1.0e-15,c_va*log(c_va),0))/(1.0*c_fe + 1.0*c_mn + 1.0*c) + 8.3145*T*if(T < 689.333333333333*c_fe*c_mn*c_va*(c_fe - c_mn) + 760.666666666667*c_fe*c_mn*c_va + 67.0*c_fe*c + 67.0*c_fe*c_va + 540.0*c_mn*c_va,-1.05647165480942e-46*T^15/(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^15 - 9.09581551984581e-29*T^9/(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^9 - 3.9645187279453e-10*T^3/(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^3 + 1 - 0.86033875460538*(689.333333333333*c_fe*c_mn*c_va*(c_fe - c_mn) + 760.666666666667*c_fe*c_mn*c_va + 67.0*c_fe*c + 67.0*c_fe*c_va + 540.0*c_mn*c_va + 1.0e-9)/T,if(T < -2068.0*c_fe*c_mn*c_va*(c_fe - c_mn) - 2282.0*c_fe*c_mn*c_va - 201.0*c_fe*c - 201.0*c_fe*c_va - 1620.0*c_mn*c_va,-7.36273260959473e-54*T^15/(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^15 - 4.62115303553615e-33*T^9/(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^9 - 1.46834026960937e-11*T^3/(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^3 + 1 - 0.86033875460538*(-2068.0*c_fe*c_mn*c_va*(c_fe - c_mn) - 2282.0*c_fe*c_mn*c_va - 201.0*c_fe*c - 201.0*c_fe*c_va - 1620.0*c_mn*c_va + 1.0e-9)/T,if(T > -2068.0*c_fe*c_mn*c_va*(c_fe - c_mn) - 2282.0*c_fe*c_mn*c_va - 201.0*c_fe*c - 201.0*c_fe*c_va - 1620.0*c_mn*c_va & 2068.0*c_fe*c_mn*c_va*(c_fe - c_mn) + 2282.0*c_fe*c_mn*c_va + 201.0*c_fe*c + 201.0*c_fe*c_va + 1620.0*c_mn*c_va < 0,-2.64184091899323e+15*(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^5/T^5 - 3.21182969926491e+47*(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^15/T^15 - 2.58302259505429e+80*(-0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) - c_fe*c_mn*c_va - 0.0880806310254163*c_fe*c - 0.0880806310254163*c_fe*c_va - 0.709903593339176*c_mn*c_va + 4.38212094653812e-13)^25/T^25,if(T > 689.333333333333*c_fe*c_mn*c_va*(c_fe - c_mn) + 760.666666666667*c_fe*c_mn*c_va + 67.0*c_fe*c + 67.0*c_fe*c_va + 540.0*c_mn*c_va & 2068.0*c_fe*c_mn*c_va*(c_fe - c_mn) + 2282.0*c_fe*c_mn*c_va + 201.0*c_fe*c + 201.0*c_fe*c_va + 1620.0*c_mn*c_va > 0,-10871773329190.3*(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^5/T^5 - 2.23837934085496e+40*(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^15/T^15 - 3.04857467251016e+68*(0.906222611744084*c_fe*c_mn*c_va*(c_fe - c_mn) + c_fe*c_mn*c_va + 0.0880806310254163*c_fe*c + 0.0880806310254163*c_fe*c_va + 0.709903593339176*c_mn*c_va + 1.31463628396144e-12)^25/T^25,0))))*log((-2.1*c_fe*c - 2.1*c_fe*c_va - 1.86*c_mn*c_va)*if(2.1*c_fe*c + 2.1*c_fe*c_va + 1.86*c_mn*c_va >= 0,-0.333333333333333,1.0) + 1)/(1.0*c_fe + 1.0*c_mn + 1.0*c) + (c_fe*c*(-15.877*T + if(T < 1811.0,-1.15*T*log(T) + 8.282*T + 0.00064*T^2.0 + if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) - 1462.4,4.9251e+30*T^(-9.0) + 0.940009*T + if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) - 1713.815) + if(T < 103.0,-0.09009204*T - 2.75e-5*T^3.0 - 1049.140845,if(T < 350.0,1.76583*T*log(T) - 7.39898691*T - 0.01706952*T^2.0 - 988.25091,12000000000.0*T^(-3.0) - 264300000.0*T^(-2.0) + 2562600.0*1/T - 24.3*T*log(T) + 170.73*T - 0.0004723*T^2.0 - 17368.441)) + 77207.0) + c_fe*c_va*if(T < 1811.0,-1.15*T*log(T) + 8.282*T + 0.00064*T^2.0 + if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) - 1462.4,4.9251e+30*T^(-9.0) + 0.940009*T + if(T < 50.0,-9.929623*T - 1.975e-6*T^4.0 + 4701.3908,if(T < 170.0,8.75517*T*log(T) - 41.511332*T - 0.13106*T^2.0 + 9.2615e-5*T^3.0 + 4871.6848,if(T < 1811.0,77359.0*1/T - 23.5143*T*log(T) + 124.134*T - 0.00439752*T^2.0 - 5.8927e-8*T^3.0 + 1225.7,2.29603e+31*T^(-9.0) - 46.0*T*log(T) + 299.31255*T - 25383.581))) - 1713.815) + c_mn*c*(13.659*T + if(T < 103.0,-0.09009204*T - 2.75e-5*T^3.0 - 1049.140845,if(T < 350.0,1.76583*T*log(T) - 7.39898691*T - 0.01706952*T^2.0 - 988.25091,12000000000.0*T^(-3.0) - 264300000.0*T^(-2.0) + 2562600.0*1/T - 24.3*T*log(T) + 170.73*T - 0.0004723*T^2.0 - 17368.441)) + if(T < 36.7,-1.0947477*T - 4.25e-6*T^4.0 - 4901.51892,if(T < 165.0,5.82736*T*log(T) - 20.019344*T - 0.12386*T^2.0 + 8.75e-5*T^3.0 - 4822.7003,if(T < 1519.0,69827.0*1/T - 23.4582*T*log(T) + 130.059*T - 0.00734768*T^2.0 - 8115.28,1.656847e+30*T^(-9.0) - 48.0*T*log(T) + 312.2648*T - 28733.41)))) + c_mn*c_va*if(T < 1519.0,69600.0*1/T - 24.5177*T*log(T) + 131.884*T - 0.006*T^2.0 - 3439.3,2.205113e+30*T^(-9.0) - 2.5984*T + if(T < 36.7,-1.0947477*T - 4.25e-6*T^4.0 - 4901.51892,if(T < 165.0,5.82736*T*log(T) - 20.019344*T - 0.12386*T^2.0 + 8.75e-5*T^3.0 - 4822.7003,if(T < 1519.0,69827.0*1/T - 23.4582*T*log(T) + 130.059*T - 0.00734768*T^2.0 - 8115.28,1.656847e+30*T^(-9.0) - 48.0*T*log(T) + 312.2648*T - 28733.41))) + 2663.31))/(1.0*c_fe + 1.0*c_mn + 1.0*c) + (c_fe*c_mn*c*(20082.0 - 11.6312*T) + c_fe*c_mn*c_va*(3.865*T - 7762.0) - 259.0*c_fe*c_mn*c_va*(c_fe - c_mn) - 34671.0*c_fe*c*c_va - 41333.0*c_mn*c*c_va)/(1.0*c_fe + 1.0*c_mn + 1.0*c); G/100000'
    coupled_variables = 'c'
    constant_names = 'T c_mn c_va'
    constant_expressions = '900 0.0001 0.01'
[]

[free_energy_A]
    type = DerivativeParsedMaterial
    property_name = Fa
    coupled_variables = 'c'
    expression = '(c-0.1)^2'
    derivative_order = 2
    enable_jit = true
[]
[free_energy_B]
    type = DerivativeParsedMaterial
    property_name = Fb
    coupled_variables = 'c'
    expression = '(c-0.9)^2'
    derivative_order = 2
    enable_jit = true
[]

# h(eta)
[h0]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta0
    function_name = h0
[]
[h1]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta1
    function_name = h1
[]
# g(eta)
[barrier]
    type = MultiBarrierFunctionMaterial
    etas = 'eta0 eta1'
[]

# constant properties
[constants]
    type = GenericConstantMaterial
    prop_names = 'M   L   kappa_eta kappa_c'
    prop_values = '0.1 1 1 1'
[]

[free_energy]
    type = DerivativeMultiPhaseMaterial
    property_name = F
    fi_names = 'fb ff' # free energy from thermodynamic database
    hi_names = 'h0  h1'
    etas     = 'eta0 eta1'
    coupled_variables = 'c'
    W = 1
[]
[]

[Kernels]
[c_res]
    type = SplitCHParsed
    variable = c
    f_name = F
    kappa_name = kappa_c
    w = w
    coupled_variables = 'eta0 eta1'
[]
[w_res]
    type = SplitCHWRes
    variable = w
    mob_name = M
[]
[time1]
    type = CoupledTimeDerivative
    variable = w
    v = c
[]

# Allen Cahn for eta_s
[deta0dt]
    type = TimeDerivative
    variable = eta0
[]
[ACBulk0]
    type = AllenCahn
    variable = eta0
    coupled_variables = 'c eta1'
    f_name = F
[]
[ACInterface0]
    type = ACInterface
    variable = eta0
    kappa_name = kappa_eta
[]
[lagrange0]
    type = SwitchingFunctionConstraintEta
    variable = eta0
    h_name   = h0
    lambda = lambda
[]

#
[deta1dt]
    type = TimeDerivative
    variable = eta1
[]
[ACBulk1]
    type = AllenCahn
    variable = eta1
    coupled_variables = 'c eta0'
    f_name = F
[]
[ACInterface1]
    type = ACInterface
    variable = eta1
    kappa_name = kappa_eta
[]
[lagrange1]
    type = SwitchingFunctionConstraintEta
    variable = eta1
    h_name   = h1
    lambda = lambda
[]

[lagrange]
    type = SwitchingFunctionConstraintLagrange
    variable = lambda
    etas    = 'eta0 eta1'
    h_names = 'h0   h1'
    epsilon = 0
[]
[]

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

[Executioner]
    type = Transient
    scheme = 'bdf2'

    solve_type = 'NEWTON'      
    line_search = none

    l_max_its = 15
    l_tol = 1.0e-6

    nl_max_its = 50
    nl_rel_tol = 1.0e-8
    nl_abs_tol = 1.0e-10

    start_time = 0.0
    num_steps = 1000
    dt = 0.1
    #dtmin = 0.01
[]

[Debug]
    # show_var_residual_norms = true
[]

[Outputs]
    execute_on = 'timestep_end'
    exodus = true
[]

My simulation runs only for some time. and reach an error

Time Step 31, time = 3.1, dt = 0.1
 0 Nonlinear |R| = 1.111361e+00
      0 Linear |R| = 1.111361e+00
      1 Linear |R| = 6.090967e-02
      2 Linear |R| = 3.624371e-03
      3 Linear |R| = 5.666513e-04
      4 Linear |R| = 9.757378e-05
      5 Linear |R| = 2.279333e-05
      6 Linear |R| = 6.711429e-06
      7 Linear |R| = 1.893094e-06
      8 Linear |R| = 6.346924e-07
 1 Nonlinear |R| = 9.296200e-03
      0 Linear |R| = 9.296200e-03
      1 Linear |R| = 4.807362e-04
      2 Linear |R| = 4.283650e-05
      3 Linear |R| = 6.014787e-06
      4 Linear |R| = 6.285949e-07
      5 Linear |R| = 1.301810e-07
      6 Linear |R| = 2.560505e-08
      7 Linear |R| = 6.969041e-09
 2 Nonlinear |R| = 9.270756e-05
      0 Linear |R| = 9.270756e-05
      1 Linear |R| = 5.288269e-06
      2 Linear |R| = 4.785061e-07
      3 Linear |R| = 8.712457e-08
      4 Linear |R| = 8.980206e-09
      5 Linear |R| = 1.801790e-09
      6 Linear |R| = 3.482315e-10
      7 Linear |R| = 9.738035e-11
      8 Linear |R| = 3.046745e-11
 3 Nonlinear |R| = 8.602855e-07
      0 Linear |R| = 8.602855e-07
      1 Linear |R| = 4.493172e-08
      2 Linear |R| = 4.348744e-09
      3 Linear |R| = 1.175112e-09
      4 Linear |R| = 4.701996e-10
      5 Linear |R| = 1.917130e-10
      6 Linear |R| = 5.935528e-11
      7 Linear |R| = 1.765063e-11
      8 Linear |R| = 4.577082e-12
      9 Linear |R| = 1.715771e-12
     10 Linear |R| = 3.829734e-13
 4 Nonlinear |R| = 8.633468e-09
 Solve Converged!

Outlier Variable Residual Norms:
  eta0: 6.104782e-09
  eta1: 6.104786e-09

Time Step 32, time = 3.2, dt = 0.1
 0 Nonlinear |R| = 1.110754e+00
      0 Linear |R| = 1.110754e+00
      1 Linear |R| = 6.094338e-02
      2 Linear |R| = 3.626445e-03
      3 Linear |R| = 5.685494e-04
      4 Linear |R| = 9.735900e-05
      5 Linear |R| = 2.273736e-05
      6 Linear |R| = 6.687574e-06
      7 Linear |R| = 1.884154e-06
      8 Linear |R| = 6.317986e-07
 1 Nonlinear |R| = 9.294519e-03
    Linear solve did not converge due to DIVERGED_PC_FAILED iterations 0
                   PC failed due to SUBPC_ERROR 
  Nonlinear solve did not converge due to DIVERGED_FNORM_NAN iterations 1
 Solve Did NOT Converge!
Aborting as solve did not converge

Time Step 32, time = 3.15, dt = 0.05
 0 Nonlinear |R| = 9.873180e-01

at this line the simulation stops.

Am I doing anything wrong?

Thank you
Ashish

Answered by ashishdhole

Oct 29, 2024

This issue was solved by changing the initial condition of lambda in

[lambda]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.0
[]
[]

I used initial_condition = 1e-6

and later the

[lagrange]
    type = SwitchingFunctionConstraintLagrange
    variable = lambda
    etas    = 'eta0 eta1'
    h_names = 'h0   h1'
    epsilon = 0
[]

epsilon = 1e-9

Thank you

View full answer

GiudGiud · 2024-10-29T03:17:47Z

GiudGiud
Oct 29, 2024
Collaborator

Hello

Let's see if a more robust PC can pass this point in the simulation.
Use -pc_type LU with a nonzero -pc_factor_shift_type.

Also output all the material properties. If one of them evaluates to NaN anywhere on the mesh, that would take out the subPC

1 reply

ashishdhole Oct 29, 2024
Author

This issue was solved by changing the initial condition of lambda in

[lambda]
    order = FIRST
    family = LAGRANGE
    initial_condition = 0.0
[]
[]

I used initial_condition = 1e-6

and later the

[lagrange]
    type = SwitchingFunctionConstraintLagrange
    variable = lambda
    etas    = 'eta0 eta1'
    h_names = 'h0   h1'
    epsilon = 0
[]

epsilon = 1e-9

Thank you

Answer selected by GiudGiud

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

SUBPC_ERROR while using multiphase model #28937

Uh oh!

{{title}}

Uh oh!

Replies: 1 comment 1 reply

Uh oh!

{{title}}

Uh oh!

Uh oh!

{{title}}

Uh oh!

Select a reply

Uh oh!

SUBPC_ERROR while using multiphase model #28937

Uh oh!

ashishdhole Oct 25, 2024

Replies: 1 comment · 1 reply

Uh oh!

GiudGiud Oct 29, 2024 Collaborator

Uh oh!

ashishdhole Oct 29, 2024 Author

ashishdhole
Oct 25, 2024

Replies: 1 comment 1 reply

GiudGiud
Oct 29, 2024
Collaborator

ashishdhole Oct 29, 2024
Author