How to apply cell volum

Hello,

I have a term cv^(2/3) in my form, where cv is the cell volumn. But my element is hexahedron, so I get the error message

 warnings.warn("Only know how to compute the cell volume of an affine cell.")

Does anyone know what I could do in this case? Thanks!

Please provide a minimum working example so that someone who wants to help can reproduce your error.

As stated by the error message, a hexahedron is not in general affine, and thus one cannot compute it just based on its vertices.

You could compute the volume through a projection or assembly of the rhs in a DG space, see Volume of an element - #4 by dokken

1 Like