Reg. generation of a mesh file with heterogeneous permeability

1

Dear Sir,

I want to assign a random distribution of material parameters to the mesh (e.g., porosity and permeability).

I have seen such a mesh in Tests\Data\Parabolic\ComponentTransport\heterogeneous\ogs5_H_2D.

How can I generate such a mesh file?

Thanks
Pavan.

1 Like
2

I am using meshio import meshio, then you can use a lines like

mesh = meshio.read(inputfile) 
mesh.cell_data['vm_permeability'] = [np.array(k_cell)]

of course behind the array there must be some logic (even random numbers, since I assume some depth-trend), evaluations of coordinates, interpolation etc.

3

Dear Dominik,

In case if I need to give a heterogeneous permeability tensor for 2D approximation, do I need to generate a 3x3 tensor or 4x1 vector for simplification?

Because if I need to give MeshElement type permeability, it might result in an error if any tensorial element (xx,yy,zz,xy) is not provided.

Let me know.

Thanks
Pavan

4

Never mind, I was able to develop the code in Python. Thank you.

5

Dear Pavan,

nice that you found a solution. Is it possible to share your solution? I think that would be helpful also for others that want to do similar things.

Best regards,
Tom

6

Dear Thomas,

Sure, why not!

But, the Python code was developed for quad elements. However, you can modify the code for other mesh element types. The code basically generates a random log-normal distribution of permeability between maximum and minimum values. You can perform the needful changes.

Thanks
Pavan.

rand_lognormal.zip (907 Bytes)