Problem about mesh

[Johnney12138]

Hello,
I was trying to solve an phase field problem. After running, it gave me an error message:
*** ERROR *** No mesh file was supplied and no generation block was provided
The complete code is attached as follows:

[Mesh]
  [gen]
    type = GeneratedMeshGenerator
    dim = 2               
    nx = 10               
    ny = 10               
    xmin = 2
    xmax = 4              
    ymin = 2
    ymax = 4
  []
[]
  
[Variables]
  [./w]       # chemical potential (concentration).
  [../]
  [./phi]     # phase
  [../]
  [./T]       # temperature
  [../]
[]
  
[ICs]
  [./w_InitialCondition]
    type = ConstantIC
    value = 0.1
    variable = w
  [../]
  [./phi_InitialCondition]
    type = SmoothCircleIC
      variable = phi
      int_width = 0.01
      x1 = 3
      y1 = 3
      radius = 0.1
      outvalue = 0
      invalue = 1
  [../]
  [./T_InitialCondition]
    type = ConstantIC
    value = 10
    variable = T
  [../]
[]
  
[Kernels]
  # Order parameter phi
  [./AC_bulk]
    type = AllenCahn      
    variable = phi        
    f_name = F            
    coupled_variables = 'w T'     
  [../]
  [./anisoACinterface1]
    type = ACInterfaceKobayashi1                 
    variable = phi                               
    mob_name = M
  [../]
  [./anisoACinterface2]
    type = ACInterfaceKobayashi2                 
    variable = phi                               
    mob_name = M
  [../]
  [./phi_dot]
    type = TimeDerivative   
    variable = phi
  [../]
  # Chemical potential (concentration). 
  [./w_dot]
    type = SusceptibilityTimeDerivative   
    variable = w
    f_name = chi                          
    coupled_variables = 'phi'
  [../]
  [./Diffusion]
    type = MatDiffusion                   
    variable = w
    diffusivity = Dchi                    
  [../]
  [./coupled_phidot]
    type = CoupledSusceptibilityTimeDerivative
    variable = w
    v = phi
    f_name = ft
    coupled_variables = 'phi'
  [../]
    # adding antitrapping current term in the model.
  [./coupled_phidot_int]
    type = AntitrappingCurrent
    variable = w
    v = phi
    f_name = rhodiff
  [../]
  # temperature equation.        
  [./coupled_tempdiff]
    type = TimeDerivative
    variable = T
  [../]
  [./coef_diffusion]
    type = CoefDiffusion
    variable = T
    coef = 10
  [../]
  [./coef_coupledtimederi]
    type = CoefCoupledTimeDerivative
    variable = T
    v = phi
    coef = -0.5
  [../]
[]
  
  [Materials]
  [./mobility]
    type = DerivativeParsedMaterial
    property_name = M
    coupled_variables = 'phi w'
    material_property_names = 'phiy:=D[phi,y] phix:=D[phi,x]'
    expression = 'M_0*(1+eps*(phix^4-6*phix^2*phiy^2+phiy^4)/(phix^2+phiy^2)^2)^2*(1/Le+NCi*(1+(1-k)*w))'
    enable_jit = true
    derivative_order = 2
  [../]
  [./buck free energy]
    type = DerivativeParsedMaterial
    property_name = F
    coupled_variables = 'phi w T'
    expression = 'M*(phi^3 - phi + lambda*(1-phi^2)^2*(NCi*w+T))'
    enable_jit = true
    derivative_order = 2
  [../]
  [./antitrapping coef]
    type = DerivativeParsedMaterial           
    args = 'w'
    f_name = rhodiff
    function = '-tc/sqrt(2)*((1-k)*w+1)'   
  [../]
  [./chi]
    type = DerivativeParsedMaterial
    args = 'phi'        
    f_name = chi
    function = '1+k-(1-k)*phi'    
  [../]
  [./Dchi]
    type = DerivativeParsedMaterial
    args = 'phi'
    f_name = Dchi
    function = 'Di*(phi-1)'        
  [../]
  [./ft]
    type = DerivativeParsedMaterial
    args = 'w'
    f_name = ft
    function = '(1-k)*w+1'
  [../]
  [./const]
    type = GenericConstantMaterial
    prop_names =  'M_0   eps    Le  NCi   k    lambda tc    Di'
    prop_values = '0.0005 0.02  1  0.15  0.15  3.1913  0.01  2'
  [../]
[]
  
[Preconditioning]
  [./SMP]
    type = SMP        
    full = true      
  [../]
[]
  
[Executioner]
  type = Transient
  scheme = bdf2           
  solve_type = PJFNK      
  petsc_options_iname = '-pc_type -pc_hypre_type -ksp_gmres_restart'
  petsc_options_value = 'hypre    boomeramg      31'
  l_tol = 1.0e-3          
  l_max_its = 100          
  nl_max_its = 100         
  nl_rel_tol = 1.0e-8     
  nl_abs_tol = 1e-8       
  # end_time = 2          
  [./TimeStepper]
    type = IterationAdaptiveDT  
    dt = 0.0005                 
    cutback_factor = 0.7        
    growth_factor = 1.2         
  [../]
[]

[Adaptivity]
initial_steps = 3                   
max_h_level = 3                     
initial_marker = err_eta           
marker = err_bnds                  
[./Markers]                          
  [./err_eta]
    type = ErrorFractionMarker     
    coarsen = 0.3
    refine = 0.95
    indicator = ind_eta
  [../]
  [./err_bnds]
    type = ErrorFractionMarker
    coarsen = 0.3
    refine = 0.95
    indicator = ind_bnds
  [../]
[../]
[./Indicators]                       
  [./ind_eta]
    type = GradientJumpIndicator     
    variable = etaa0
  [../]
  [./ind_bnds]
    type = GradientJumpIndicator    
    variable = bnds
  [../]
[../]
[]
    
[Outputs]
  interval = 10     
  exodus = true
[]

I am confused since I did had the mesh block. Is there any wrong with my input? Thanks!

[GiudGiud]

Hello

This should not happen, as you do have a Mesh block
Can you type the exact command you use here to run the code?
You may find it in bash_history, a file often at ~/.bash_history

Guillaume

[Johnney12138]

Hello @GiudGiud
Thank you for pointing out the problem with the code. If I want to declare that phi is a function of x and y, which block do I need to write such code?

Nan Liu

[GiudGiud]

x and y can be written directly into the expression
it's these dpx and dpy that need to be specified either as a function, postprocessor etc

[Johnney12138]

Hello @GiudGiud
Thank you for your reply.
I still haven't solved the problem. In order to facilitate understanding, I use mathematical expressions to describe this problem. In a phase field model, Allen-Cahn equation has an anisotropic mobility. The anisotropic mobility is given as follows:
Tau=tau_0A^2{1/Le+Mc*[1+(1-k)*U]},

where tau_0, Le, Mc, and k are constants, and U is temperature. A=1+epicos(4theta) is a function that describes the anisotropy and theta is the angle between the direction normal to the interface and the x (horizontal) axis. In Moose, I don't seem to find objects or examples of how to deal with the anisotropic mobility of the AC equation. I tried to write the form of mobility directly in the code, but there was a problem with the expression of theta. I would like to know, in moose, are there objects that can represent the mobility or theta of the AC equation?

Thanks!
Nan Liu

[GiudGiud]

I m not familiar with that.
Maybe @dschwen or @recuero

[Johnney12138]

Hello @GiudGiud
Thank you for your reply.

Nan Liu