nucleated grains not captured by grain tracker

[ashishdhole]

Hello,

I am working on phase transformation with nucleation in a two phase system. I am using multiphase model for the phase transformation. I want to nucleate new grains with particular phase in the system of poly crystals.
following is the input file I made

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

[GlobalParams]
    op_num = 4
    var_name_base = eta
[]

[Variables]
    # carbon concentration
[c]
    order = FIRST
    family = LAGRANGE
[InitialCondition]
    type = ConstantIC
    value = 0.1
[]
[]
[w]
    order = FIRST
    family = LAGRANGE
[]
[PolycrystalVariables]
    order = FIRST
    family = LAGRANGE
[]

# Lagrange multiplier
[lambda]
    order = FIRST
    family = LAGRANGE
    initial_condition = 1e-5
[]
[]

[ICs]
[PolycrystalICs]
[PolycrystalColoringIC]
    polycrystal_ic_uo = voronoi
[]
[]
[]
[AuxVariables]
[bnds]
[]
[unique_grains]
    order = CONSTANT
    family = MONOMIAL
[]
[]
[AuxKernels]
[BndsCalc]
    type = BndsCalcAux
    variable = bnds
    execute_on = 'timestep_end'
[]
[unique_grains]
    type = FeatureFloodCountAux
    variable = unique_grains
    field_display = UNIQUE_REGION
    execute_on = 'initial timestep_end'
    flood_counter = grain_tracker
[]
[]
[UserObjects]
[grain_tracker]
    type = GrainTracker
    verbosity_level = 0
    compute_var_to_feature_map = true
    reserve_op = 1
[]
[voronoi]
    type = PolycrystalVoronoi
    grain_num = 3 # Number of grains
    rand_seed = 123
    coloring_algorithm = bt
[]
[inserter]
    type = DiscreteNucleationFromFile
    hold_time = 100
    file = nuclei2.csv
    radius = 4
[]
[map]
    type = DiscreteNucleationMap
    int_width = 0.5
    periodic = eta3
    inserter = inserter
[]
[]
[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 = '1100 0.0001 0.001'
[]
[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 = '1100 0.0001 0.001'
[]

[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
[]
[nucleation]
    # The nucleation material is configured to insert nuclei into the free energy
    # tht force the concentration to go to 0.95, and holds this enforcement for 500
    # time units.
    type = DiscreteNucleation
    property_name = Fn
    op_names  = eta3
    op_values = 0.90
    map = map
    outputs = exodus
[]
[free_energy_nuc]
    # add the chemical and nucleation free energy contributions together
    type = DerivativeSumMaterial
    derivative_order = 2
    coupled_variables = 'eta3 c'
    sum_materials = 'Fa Fn'
    property_name = Fan
[]

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

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

[free_energy]
    type = DerivativeMultiPhaseMaterial
    property_name = F
    fi_names = 'Fa Fb Fb Fan' # free energy from thermodynamic database
    hi_names = 'h0  h1 h2 h3'
    etas     = 'eta0 eta1 eta2 eta3'
    coupled_variables = 'c'
    W = 1
[]
[]

[Kernels]
[eta_force]
    type = DiscreteNucleationForce
    variable = eta3
    map = map
    no_nucleus_value = 0
    nucleus_value = 0.9
[]
[eta_react]
    type = Reaction
    variable = eta3
[]
[c_res]
    type = SplitCHParsed
    variable = c
    f_name = F
    kappa_name = kappa_c
    w = w
    coupled_variables = 'eta0 eta1 eta2 eta3'
[]
[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 eta2 eta3'
    f_name = F
[]
[ACInterface0]
    type = ACInterface
    variable = eta0
    coupled_variables = 'eta1 eta2 eta3'
    kappa_name = kappa_eta
[]
[lagrange0]
    type = SwitchingFunctionConstraintEta
    variable = eta0
    h_name   = h0
    lambda = lambda
    coupled_variables = 'eta1 eta2 eta3'
[]

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

#
[deta2dt]
    type = TimeDerivative
    variable = eta2
[]
[ACBulk2]
    type = AllenCahn
    variable = eta2
    coupled_variables = 'c eta0 eta1 eta3'
    f_name = F
[]
[ACInterface2]
    type = ACInterface
    variable = eta2
    kappa_name = kappa_eta
    coupled_variables = 'eta0 eta1 eta3'
[]
[lagrange2]
    type = SwitchingFunctionConstraintEta
    variable = eta2
    h_name   = h2
    lambda = lambda
    coupled_variables = 'eta1 eta0 eta3'
[]

#
[deta3dt]
    type = TimeDerivative
    variable = eta3
[]
[ACBulk3]
    type = AllenCahn
    variable = eta3
    coupled_variables = 'c eta0 eta1 eta2'
    f_name = F
[]
[ACInterface3]
    type = ACInterface
    variable = eta3
    kappa_name = kappa_eta
    coupled_variables = 'eta0 eta1 eta2'
[]
[lagrange3]
    type = SwitchingFunctionConstraintEta
    variable = eta3
    h_name   = h3
    lambda = lambda
    coupled_variables = 'eta1 eta2 eta0'
[]

[lagrange]
    type = SwitchingFunctionConstraintLagrange
    variable = lambda
    etas    = 'eta0 eta1 eta2 eta3'
    h_names = 'h0   h1 h2 h3'
    epsilon = 1e-6
[]
[]

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

I was able to get the phase transformation effect without nucleation. But I think I am missing soething in the nucleation blocks that is preventing it to capture the newly formed nuclei.

at t=0

at t=10

but nucleation shows no effect

Thank you
Ashish

[GiudGiud]

@laagesen

[ashishdhole]

Hello,

Can you explain what is happening with my nucleus?

This is my input file

[Mesh]
    type = GeneratedMesh
    dim = 2
    nx = 150
    ny = 150
    xmin = 0
    xmax = 100
    ymin = 0
    ymax = 100
    elem_type = QUAD4
    parallel_type = DISTRIBUTED
[]

[GlobalParams]
    op_num = 5
    var_name_base = eta
[]

[Variables]
    # carbon concentration
[c]
    order = FIRST
    family = LAGRANGE
[InitialCondition]
    type = ConstantIC
    value = 0.1900000000
[]
[]
[w]
    order = FIRST
    family = LAGRANGE
[]
[PolycrystalVariables]
    order = FIRST
    family = LAGRANGE
[]

# Lagrange multiplier
[lambda]
    order = FIRST
    family = LAGRANGE
    initial_condition = 1e-5
[]
[]
#[Modules]
#[PhaseField]
#[EulerAngles2RGB]
#    crystal_structure = cubic
#    euler_angle_provider = euler_angle_file
#    grain_tracker = grain_tracker
#[]
#[]
#[]
[ICs]
[PolycrystalICs]
[PolycrystalColoringIC]
    polycrystal_ic_uo = voronoi
[]
[]
[]
[AuxVariables]
[bnds]
[]
[unique_grains]
    order = CONSTANT
    family = MONOMIAL
[]
[]
[AuxKernels]
[BndsCalc]
    type = BndsCalcAux
    variable = bnds
    execute_on = 'timestep_end'
[]
[unique_grains]
    type = FeatureFloodCountAux
    variable = unique_grains
    field_display = UNIQUE_REGION
    execute_on = 'initial timestep_end'
    flood_counter = grain_tracker
[]
[]
[UserObjects]
[grain_tracker]
    type = GrainTracker
    verbosity_level = 0
    compute_var_to_feature_map = false
    reserve_op = 1
    #connecting_threshold = 0.1
    remap_grains = true
[]
[voronoi]
    type = PolycrystalVoronoi
    grain_num = 4 # Number of grains
    rand_seed = 8
    coloring_algorithm = bt
[]
#[inserter1]
#    type = DiscreteNucleationInserter
#    hold_time = 100
#    probability = P1
#    seed = 12346
#    radius = 4
#[]
#[inserter2]
#    type = DiscreteNucleationInserter
#    hold_time = 100
#    probability = P2
#    seed = 123
#    radius = 4
#[]
[inserter1]
    type = DiscreteNucleationFromFile
    hold_time = 100
    file = nuclei1.csv
    radius = 3
[]
#[inserter2]
#    type = DiscreteNucleationFromFile
#    hold_time = 100
#    file = nuclei2.csv
#    radius = 3
#[]
[map1]
    type = DiscreteNucleationMap
    int_width = 0.2
    periodic = eta4
    inserter = inserter1
[]
#[map2]
#    type = DiscreteNucleationMap
#    int_width = 0.1
#    #periodic = eta4
#    inserter = inserter2
#[]
#[euler_angle_file]
#    type = RandomEulerAngleProvider
#    grain_tracker_object = grain_tracker
#    execute_on = 'initial timestep_begin'
#[]
[]
[BCs]
[Periodic]
[all]
    auto_direction = 'x y'
[]
[]
[]

[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 = '1058 0.0001 0.001'
[]
[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 = '1058 0.0001 0.001'
[]

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

#[probability1]
#    # This is a made up toy nucleation rate it should be replaced by
#    # classical nucleation theory in a real simulation.
#    type = ParsedMaterial
#    property_name = P1
#    coupled_variables = eta0
#    expression = 'if(eta0 < 0.6, if(eta0 > 0.4, 0.001, 0), 0)'
#    outputs = exodus
#[]
#[probability2]
#    # This is a made up toy nucleation rate it should be replaced by
#    # classical nucleation theory in a real simulation.
#    type = ParsedMaterial
#    property_name = P2
#    coupled_variables = eta1
#    expression = 0#'if(eta1 < 0.6, if(eta1 > 0.4, 0.0000, 0), 0)'
#    outputs = exodus
#[]
# h(eta)
[h0]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta0
    function_name = h0
[]
[h1]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta1
    function_name = h1
[]
[h2]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta2
    function_name = h2
[]
[h3]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta3
    function_name = h3
[]
[h4]
    type = SwitchingFunctionMaterial
    h_order = HIGH
    eta = eta4
    function_name = h4
[]
# g(eta)
[barrier]
    type = MultiBarrierFunctionMaterial
    etas = 'eta0 eta1 eta2 eta3 eta4'
[]

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

[free_energy]
    type = DerivativeMultiPhaseMaterial
    property_name = F
    fi_names = 'Fb Fb Fb Fb Fa' # free energy from thermodynamic database
    hi_names = 'h0  h1 h2 h3 h4'
    etas     = 'eta0 eta1 eta2 eta3 eta4'
    coupled_variables = 'c'
    W = 1
[]

[material_anisotropic]
    type = InterfaceOrientationMaterial
    op = eta4
    reference_angle = 45
    anisotropy_strength = 0.045
    eps_bar = 0.4
    mode_number = 1
[]
[]

[Kernels]
    #[c_force1]
    #    
    #    type = DiscreteNucleationForce
    #    variable = eta3
    #    map = map1
    #       
    #    no_nucleus_value = 0.1
    #    nucleus_value = 0.9
    #  
    #[]
    #[eta_react1]
    #    type = Reaction
    #    variable = eta3
    #[]
[c_force2]
    
    type = DiscreteNucleationForce
    variable = eta4
    map = map1
       
    no_nucleus_value = 0
    nucleus_value = 0.9
  
[]
[eta_react2]
    type = Reaction
    variable = eta4
[]
[c_res]
    type = SplitCHParsed
    variable = c
    f_name = F
    kappa_name = kappa_c
    w = w
    coupled_variables = 'eta0 eta1 eta2 eta3 eta4'
[]
[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 eta2 eta3 eta4'
    f_name = F
    mob_name = L
[]
[ACInterface0]
    type = ACInterface
    variable = eta0
    coupled_variables = 'eta1 eta2 eta3 eta4'
    kappa_name = kappa_eta
[]
[lagrange0]
    type = SwitchingFunctionConstraintEta
    variable = eta0
    h_name   = h0
    lambda = lambda
    coupled_variables = 'eta1 eta2 eta3 eta4'
[]

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

#
[deta2dt]
    type = TimeDerivative
    variable = eta2
[]
[ACBulk2]
    type = AllenCahn
    variable = eta2
    coupled_variables = 'c eta0 eta1 eta3 eta4'
    f_name = F
    mob_name = L
[]
[ACInterface2]
    type = ACInterface
    variable = eta2
    kappa_name = kappa_eta
    coupled_variables = 'eta0 eta1 eta3 eta4'
[]
[lagrange2]
    type = SwitchingFunctionConstraintEta
    variable = eta2
    h_name   = h2
    lambda = lambda
    coupled_variables = 'eta1 eta0 eta3 eta4'
[]

#
[deta3dt]
    type = TimeDerivative
    variable = eta3
[]
[ACBulk3]
    type = AllenCahn
    variable = eta3
    coupled_variables = 'c eta0 eta1 eta2 eta4'
    f_name = F
    mob_name = L
[]
[ACInterface3]
    type = ACInterface
    variable = eta3
    kappa_name = kappa_eta
    coupled_variables = 'eta0 eta1 eta2 eta4'
[]
[lagrange3]
    type = SwitchingFunctionConstraintEta
    variable = eta3
    h_name   = h3
    lambda = lambda
    coupled_variables = 'eta1 eta2 eta0 eta4'
[]

#
[deta4dt]
    type = TimeDerivative
    variable = eta4
[]
[ACBulk4]
    type = AllenCahn
    variable = eta4
    coupled_variables = 'c eta0 eta1 eta2 eta3'
    f_name = F
    mob_name = L
[]
[ACInterface4]
    type = ACInterface
    variable = eta4
    kappa_name = kappa_eta
    coupled_variables = 'eta0 eta1 eta2 eta3'
[]
#[anisoACinterface1]
#    type = ACInterfaceKobayashi1
#    variable = eta4
#    mob_name = L
#    coupled_variables = 'eta0 eta1 eta2 eta3'
#[]
#[anisoACinterface2]
#    type = ACInterfaceKobayashi2
#    variable = eta4
#    mob_name = L
#    coupled_variables = 'eta0 eta1 eta2 eta3'
#[]
[lagrange4]
    type = SwitchingFunctionConstraintEta
    variable = eta4
    h_name   = h4
    lambda = lambda
    coupled_variables = 'eta1 eta2 eta0 eta3'
[]

[lagrange]
    type = SwitchingFunctionConstraintLagrange
    variable = lambda
    etas    = 'eta0 eta1 eta2 eta3 eta4'
    h_names = 'h0 h1 h2 h3 h4'
    epsilon = 1e-6
[]
[]

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

[Executioner]
    type = Transient
    scheme = 'bdf2'

    solve_type = 'newton'

    l_max_its = 30
    l_tol = 1.0e-3

    nl_max_its = 30
    nl_rel_tol = 1.0e-3
    nl_abs_tol = 1.0e-3

    start_time = 0.0
    num_steps = 100000
[TimeStepper]
    type = IterationAdaptiveDT
    cutback_factor = 0.1
    growth_factor = 1.1
    dt = 0.05
    optimal_iterations = 8
[]
[]

[Debug]
    # show_var_residual_norms = true
[]

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

output.mov

Thank you

[GiudGiud]

Do you mean the oscillations? Are they physical or numerical?

[ashishdhole]

Yes. The newly formed grain is supposed to stay and grow based on the free energy I have supplied. But it seems multiple nucleation are taking place at the same place. I think it is physical because the free energy is generic single well function.

[GiudGiud]

I'd be concerned about numerical instabilities. How converged is this? Does it still appear with a tight convergence criterion?

@dschwen is the "multiple nucleation" scenario a possibility?

[ashishdhole]

So in the meantime I did some changes they are

remap_grains = false

in GrainTracker

and commented #periodic = eta4
from Map1 in the input file.

This is what I am getting now

output.mov

But if I look at the grain boundaries only, I am still not sure if this is correct. you can see the grey region inside newly formed grain.

output1.mov