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
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.
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
Never mind, I was able to develop the code in Python. Thank you.