Define a non-cartesian mesh in fipy

  fipy, python-3.x

I am trying to simulate the elementary unit of a 2D system that have a P6mm symmetry in fipy and I would like to define a non-cartesian mesh that describes the system described below. Yet,

mesh = fipy.Grid2D(nx = 10, ny = 10, dx = 1., dy = 1.)

only returns uniform meshes. I thought of changing the FaceVariable, but it seems it only accepts Boolean variables. I could also simulate a cartesian system equivalent to this one, but there would be redondant data. Would someone have a better approach ?

System that I want to modelise, with is a mesoporous silica. The different colours indicate different domains, that are regions with different diffusion coefficients.

Source: Python-3x Questions

LEAVE A COMMENT