Hi Henrik, That is indeed the claim. Let us know if it works.... ;) David On Fri, 16 Oct 2015 at 12:09 Buesing, Henrik <HBuesing@eonerc.rwth-aachen.de> wrote:
Dear Firedrakers,
Since https://github.com/firedrakeproject/firedrake/pull/600 has been merged will
*bc0 = DirichletBC(W.sub(0), Expression("1.0e7"), 2,method="geometric")*
now work on a DG space on an extruded mesh (see below for actual code)?
Do I just have to update firedrake for this to work?
Thank you!
Henrik
--
Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
------------------------------------------------------
Mathieustr. 10 | Tel +49 (0)241 80 49907
52074 Aachen, Germany | Fax +49 (0)241 80 49889
------------------------------------------------------
http://www.eonerc.rwth-aachen.de/GGE
hbuesing@eonerc.rwth-aachen.de
------------------------------------------------------
*Von:* firedrake-bounces@imperial.ac.uk [mailto: firedrake-bounces@imperial.ac.uk] *Im Auftrag von *Buesing, Henrik *Gesendet:* 23 September 2015 12:06 *An:* firedrake@imperial.ac.uk *Betreff:* [firedrake] Dirichlet BC on extruded mesh with discontinuous function space
Dear Firedrakers,
I’m constructing a DG space on an extruded mesh.
*meshbase = RectangleMesh(Nx, Ny, Lx, Ly, quadrilateral=True) mesh = ExtrudedMesh(meshbase, Nz, Delta_z) horiz_elt = FiniteElement("DG", quadrilateral, 0) vert_elt = FiniteElement("DG", interval, 0) elt = OuterProductElement(horiz_elt, vert_elt)*
*DG = FunctionSpace(mesh, elt)*
*W=DG*DG*
Now I would like to impose Dirchlet BC on one of the boundaries.
I tried
*bc0 = DirichletBC(W.sub(0), Expression("1.0e7"), 2,method="geometric")*
and get a
*ValueError: Geometric boundary conditions are not yet supported on extruded meshes*
Using the standard topological way, seems to have no effect at all.
Is there a way around this? I would be fine with fixing the value on the most outer cell layer, if that is possible…
Henrik
--
Dipl.-Math. Henrik Büsing
Institute for Applied Geophysics and Geothermal Energy
E.ON Energy Research Center
RWTH Aachen University
------------------------------------------------------
Mathieustr. 10 | Tel +49 (0)241 80 49907
52074 Aachen, Germany | Fax +49 (0)241 80 49889
------------------------------------------------------
http://www.eonerc.rwth-aachen.de/GGE
hbuesing@eonerc.rwth-aachen.de
------------------------------------------------------