MOOSE is not working on porousflow through 1D fractures in 2D domain.

[Yulongggggg]

The .msh files look fine, trying to simulate the 1D fractures with 1D element in MOOSE and put them into the BLOCK=, however the moose is not working and is having error:

*** ERROR ***
Materials/pf_temp_fracs/block: the following blocks (ids) do not exist on the mesh: fractures (65535)
Blocks names (resp. ids) that do exist: cap1, cap2, block (4243, 4244, 4245)

Here's the GMSH file and .i files. https://github.com/Yulongggggg/THM

[GiudGiud]

What does --mesh-only show when loading the mesh through a FileMeshGenerator?

[Yulongggggg]

The output .e file looks right too

[GiudGiud]

fractures is a sideset in the mesh not a subdomain.
if you need it to be a subdomain you can concert it with this mesh generator
https://mooseframework.inl.gov/source/meshgenerators/LowerDBlockFromSidesetGenerator.html

else you should consider where you are using fractures as a subdomain (for example kernels) and see if it should be used as a sideset (for a boundary condition for example)

[Yulongggggg]

i see, I am trying to model the fracture networks for geothermal, if it becomes a 1D element, I can use it as: block = fractures, right?

[Yulongggggg]

I am writing my own python script to generate the gmsh files, thank you for the prompt reply.

[Yulongggggg]

fractures is a sideset in the mesh not a subdomain. if you need it to be a subdomain you can concert it with this mesh generator https://mooseframework.inl.gov/source/meshgenerators/LowerDBlockFromSidesetGenerator.html

else you should consider where you are using fractures as a subdomain (for example kernels) and see if it should be used as a sideset (for a boundary condition for example)

Hi I don't know if I could do this:
[Mesh]
[from_file]
type = FileMeshGenerator
file = FN.msh
[]
[fracture_block]
type = LowerDBlockFromSidesetGenerator
input = from_file
sidesets = 'fractures'
new_block_id = 5000
new_block_name = 'fractures'
[]
[]