DG0 elements in 3D mesh ; TensorProductCell error
Hi all, Could somebody help me with this? I'm having the following error msg when trying to assemble the bilinear form bellow. a = jump(u)*jump(v)*dS petsc_a = assemble(a).M.handle ValueError: interior_facet integral cannot be used with a
TensorProductCell; need to distinguish between vertical and horizontal contributions.
I'm using a 3D mesh built as an extruded quadrilateral mesh. u and v are in DG0 function space - FunctionSpace(mesh, "DG", 0). Ultimately, I'd like to solve a heterogeneous diffusion equation in 3D. Thanks in advance -- Rodrigo
On Mon, May 18, 2020 at 4:56 PM Rodrigo Piccinini <rbpiccinini@gmail.com> wrote:
Hi all, Could somebody help me with this?
I'm having the following error msg when trying to assemble the bilinear form bellow.
a = jump(u)*jump(v)*dS petsc_a = assemble(a).M.handle
ValueError: interior_facet integral cannot be used with a
TensorProductCell; need to distinguish between vertical and horizontal contributions.
The error message is saying that dS is not alright because integrating over a vertical face is different than over a horizontal face. There must be different measures for them, like dH and dV, but I don't know what they are. Thanks, Matt
I'm using a 3D mesh built as an extruded quadrilateral mesh. u and v are in DG0 function space - FunctionSpace(mesh, "DG", 0).
Ultimately, I'd like to solve a heterogeneous diffusion equation in 3D.
Thanks in advance -- Rodrigo _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
Thanks, Matthew. I just found a doc page about this issue. https://www.firedrakeproject.org/extruded-meshes.html Em seg, 18 de mai de 2020 18:01, Matthew Knepley <knepley@gmail.com> escreveu:
On Mon, May 18, 2020 at 4:56 PM Rodrigo Piccinini <rbpiccinini@gmail.com> wrote:
Hi all, Could somebody help me with this?
I'm having the following error msg when trying to assemble the bilinear form bellow.
a = jump(u)*jump(v)*dS petsc_a = assemble(a).M.handle
ValueError: interior_facet integral cannot be used with a
TensorProductCell; need to distinguish between vertical and horizontal contributions.
The error message is saying that dS is not alright because integrating over a vertical face is different than over a horizontal face. There must be different measures for them, like dH and dV, but I don't know what they are.
Thanks,
Matt
I'm using a 3D mesh built as an extruded quadrilateral mesh. u and v are in DG0 function space - FunctionSpace(mesh, "DG", 0).
Ultimately, I'd like to solve a heterogeneous diffusion equation in 3D.
Thanks in advance -- Rodrigo _______________________________________________ firedrake mailing list firedrake@imperial.ac.uk https://mailman.ic.ac.uk/mailman/listinfo/firedrake
-- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
participants (2)
-
Matthew Knepley
-
Rodrigo Piccinini