Hi,
In am using FLAC3D in Thermal only mode.
I am trying to compute the thermal heat flux along a boundary were I have fixed the temperature at gridpoints.
The way I do it:
- I collect all the zones that have a face along the boundary.
- Get the zone average heat flux vector, take the component of this vector that is normal to the face and multiply it by the face area. This gives me the boundary heat flux for the given face.
Now comes my question. When zones shared only an edge with the boundary such as in the figure below ( in red zones sharing a face with the boundary and in green zones only sharing and edge) how to handle it?
From my understanding of the numerical formulation in FLAC3D these zones should be ignored as they do not contribute to the boundary heat flux.
Would someone be able to tell me if this approach is correct or if I am misunderstanding the nodal formulation of the energy equation in FLAC3D and my approach is wrong.
To give context of why I am trying to to this it is because I wish coupling FLAC3d with another model along the boundary. Flac3d receives a temperature along the boundary and I need to pass the heat flux along the boundary to the other model.
Thank you in advance!