When using the dual mesh (vertex-centered) for finite volume methods, you end up with a cell center at the boundaries between materials. It is possible that the equations being solved in each material are different, or at the very least don't have continuous derivatives at the boundary (for example: poisson's equation for electric potential moving from gas to metal).

My question is, how do I handle these cells. Do I split the cell in half so each partial cell is entirely contained in one material or do I treat the boundaries of the cell as the material boundaries (i.e. purely gas cells see the cell boundaries as metal and metal cells see the cell boundaries as gas)? If I split the cell, do I move the cell center for each off of the new boundary? Any book or reference handling this would be greatly appreciated, since all the references I have found for dual mesh FVM limit themselves to dirichlet conditions, and don't hanlde interfaces at all.

  • $\begingroup$ Do you have same equation but different material properties ? $\endgroup$ – cfdlab Apr 8 '14 at 14:33
  • $\begingroup$ Sometimes. For poisson's equation, yes. For the energy equation, I have standard conductivity in the metal, but in the gas I have it split between the neutral gas, and electrons (so 2 equations). Others (cons. of mass) are only solved in one of the materials. $\endgroup$ – Godric Seer Apr 8 '14 at 14:47
  • $\begingroup$ I have seen something like this, but unfortunately I cannot remember where. They resolved the interface by the control-volumes and introduced a redundant variable at the control-volume faces located at the interface. $\endgroup$ – Christian Waluga Apr 9 '14 at 9:33

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.