Mesh issue: higher order U using a Cubit mesh in THM coupling

[Khaksarkn]

Hi, I am trying to use a second-order displacement with first-order pressure and temperature in my mesh. This is the part of the input file I use. It is OK with meshes generated by MOOSE's generated mesh, but it gives me this weird shape when I do this with Cubit's generated mesh. It assumes zero for values where displacement exists but p and T don't. Can you help me specify this to avoid this issue? Thanks.

just doing second_order = true results in this issue. even without order = SECOND.

Is it something I should do with Cubit? or some other specification I need to do in the input file?

I believe that it is not a visualization issue. MOOSE is confusing the file mesh. if it was just a visualization issue, the displacements should been correct but displacement results are also meaningless. So, I guess there is an issue with how I specify the input file.

[Mesh]
  type = FileMesh
  file = fiveBores.e

  second_order = true 
[]

[Variables]
  [temperature]
    initial_condition = 274.15
    family = LAGRANGE
    order = FIRST
  []
  [porePressure]
    family = LAGRANGE
    order = FIRST
  []

  [disp_x]
    family = LAGRANGE
    order = SECOND
  []
  [disp_y]
    family = LAGRANGE
    order = SECOND
  []

But this way of giving the mesh is completely OK in visualization:

[Mesh]
  type = GeneratedMesh
  dim = 2
  nx = 5
  ny = 5
  xmin = -1
  xmax = 1
  ymin = -2
  ymax = 0
  second_order = true 
[]

[Variables]
  [temperature]
    initial_condition = 274.15
    family = LAGRANGE
    order = FIRST
    scaling = 1e3
  []
  [porePressure]
    family = LAGRANGE
    order = FIRST
    scaling = 1e9
  []

  [disp_x]
    family = LAGRANGE
    order = SECOND
  []
  [disp_y]
    family = LAGRANGE
    order = SECOND
  []