Hi everyone,
I come back to you because I would like to extract from a mesh, made
with GMSH and imported via the
GmshImporter3D() function, a specific region. My goal is to use this
information to impose a specific
BC (convective!) to this area.
I looked at the FiPy documentation and mails about imposing BCs and
I'm a little bit confusing with
the use of a mask with my mesh which is 3D and not cubic i.e. you
don't have a front face, a back face and so on !
So I have two questions:
1. In GMSH you can group some mesh elements and name them (usually
with a code number).
Is it possible to get this code number within FiPy to be able to work
on these elements specifically.
For instance, a method like mesh.getFaceColor() ?
If not, I guess that the unique solution is to used face coordinates
as in many examples you gave?
2. Once my faces were extracted and stored in a variable, for example,
maskmyfaces, I would like to impose
a convective flux (h(T-Tair)) like in this example which works for a
cubic like geometry
maskmyfaces= mesh.getFacesTop()
outCoeffTop = (maskmyfaces* [[0], [h]]).getDivergence()
TopConvec = ImplicitSourceTerm(coeff=outCoeffTop) - outCoeffTop * Tair
My question is about the variable outCoeffTop which is, in this
particular case, related to the y direction. But
how to define the same variable when the flux direction is imposed by
the normal direction of faces.
Thanks in advance for your help
Best regards
Marc
