July 2021 - firedrake - imperial.ac.uk

installing older version of firedrake
by Miroslav Kuchta 28 Jul '21

28 Jul '21

Hi, is there a recommended way of installing older versions of firedrake? For example, I have previously (some months ago) ran the firedrake-install script on a different system and have the downloaded firedrake folder backed up. Now with a newer version of ubuntu can I somehow reuse the old script/source code downloaded? To provide some context for the question, I am trying to solve Biot equations with stabilized BDM_1-DG_0 elements and observe a different behavior between the versions. However, it's not only firedrake that is different but also the underlying system so I'd like to see where the error comes from. Thanks, Miro

4 5

Problem with Dirichlet BC on edges in a 3D mesh
by Venketeswaran, Abhishek (FN) 19 Jul '21

19 Jul '21

Hi, I created a a couple of physical groups within GMSH with appropriate labels/tags. I then tried to use the same labels after importing the .msh file in firedrake. I got the following error: Please find below the relevant section in the .py file that generates the gmsh mesh and saves it in .msh file gmsh.option.setNumber("Mesh.MshFileVersion", 2) # Physical curves defined with appropriate tags starting from 1000 gmsh.model.addPhysicalGroup(1, b_o_c_arcs,1001) gmsh.model.addPhysicalGroup(1, b_i_c_arcs,1002) gmsh.model.addPhysicalGroup(1, list(t_o_c_arcs),1003) gmsh.model.addPhysicalGroup(1, list(t_i_c_arcs),1004) # Physical surfaces defined with appropriate tags starting from 2000 gmsh.model.addPhysicalGroup(2, circ_pl_surf,2001) gmsh.model.addPhysicalGroup(2, extsurfs,2002) # Physical volumes defined with appropriate tags starting from 3000 gmsh.model.addPhysicalGroup(3, voltags,3001) gmsh.model.mesh.generate(3) gmsh.write("Pipe3D.msh") Here is the snippet of the .py file executed within firedrake: mesh=Mesh('Pipe3D.msh') bc1=DirichletBC(V,zero,1003) When I execute the solve command I get this error: Error: {1003} are not a valid markers (not in {2001,2002}) I think I have a workaround that is not yet giving me errors wherein I define a class derived from DirichletBC and redefine the @utils.cached_property nodes function. However, It would be great if I could do it using the existing DirichletBC class. I would be very grateful if you could help me out here. Thanks -- Abhishek Venketeswaran NETL Research Associate - ORISE National Energy Technology Laboratory Department of Energy abhishek.venketeswaran(a)netl.doe.gov Work: 412-386-4833 Mobile: 716-507-7890 www.netl.doe.gov<http://www.netl.doe.gov> [cid:0311bcdf-830c-456c-a671-ff3e827d18f2]

3 4

Firedrake tutorial on ARCHER2, Online, August 2021
by Betteridge, Jack D 15 Jul '21

15 Jul '21

The Firedrake team will be delivering an online training session on 23rd August 2021, 13:00 - 17:00 BST (12:00 - 16:00 UTC). This course will introduce users to Firedrake and will also be an opportunity to get hands on experience using ARCHER2, the UK national supercomputer. Participants will gain experience writing Firedrake code to solve PDEs, configure solvers and preconditioners for HPC and run the resultant code on ARCHER2. Our course is aimed at researchers who have little or no previous experience of using Firedrake and have a desire to run on HPC. For further details and registration see: https://www.archer2.ac.uk/training/courses/210823-firedrake/

1 0

Building FireDrake + RedHat (Centos,Fedora)
by Régio Pires 14 Jul '21

14 Jul '21

Dear Support, Please, can you help me... FireDrake is compatible with RedHat? Are there a form to build them on this Operation System? (sorry my english) Best Regards -- Régio da Silva Pires Laboratório Nacional de Computação Científica - LNCC regiopires(a)lncc.br 24.2233-6285

2 1

Re: [firedrake] [EXTERNAL] Reg spectral (hp) tetrahedral element
by Lawrence Mitchell 09 Jul '21

09 Jul '21

I realised I forget to leave the mailing list in CC, reintroducing it. > On 9 Jul 2021, at 18:21, Venketeswaran, Abhishek (FN) <Abhishek.Venketeswaran(a)netl.doe.gov> wrote: > > Hi Lawrence, > Thanks for getting back to me. I am currently working on simulation of ultrasonic guided waves on the surface of cylindrical pipes. I am currently using hexahedral elements with GLL quadrature within firedrake to perform this simulation. > Next, I am looking to introduce notch and hole defects on the pipe wall and rerun the simulation. I am not sure on how to create such defects on the pipe wall with the current limitation of defining hexahedral elements using only extrusion of the 2D geometry. I thought, I could circumnavigate this problem by meshing the 3D geometry using tetrahedral elements and then leverage the benefits of spectral element method using the collapsed tetrahedral elements. Kindly let me know how I can solve this problem in Firedrake. Ah, I see. I think that while you can do the 3D meshing with tets, the spectral element stuff is a lot of implementation away from working. If you want to give it a go we can try and provide advice. We're hopeful that unstructured hexahedral meshes will be supported by the end of the year, so that might be swifter! Lawrence

2 3

Reg spectral (hp) tetrahedral element
by Venketeswaran, Abhishek (FN) 08 Jul '21

08 Jul '21

Hi, I would like to use the spectral (hp) tetrahedral elements proposed in the paper below: "A NEW TRIANGULAR AND TETRAHEDRAL BASIS FOR HIGH-ORDER (hp) FINITE ELEMENT METHODS, SPENCER J. SHERWIN AND GEORGE EM KARNIADAKIS, INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 38, 3775-3802 (1995) I believe that I most probably saw in the UFL documentation that this element formulation was already implemented. However, I am not able to track that documentation down now. Can you kindly help me with how I may implement this element in Firedrake? Thanks -- Abhishek Venketeswaran NETL Research Associate - ORISE National Energy Technology Laboratory Department of Energy abhishek.venketeswaran(a)netl.doe.gov Work: 412-386-4833 Mobile: 716-507-7890 www.netl.doe.gov<http://www.netl.doe.gov> [cid:9f512d72-18a5-491c-b5bd-71b2da00af09]

1 0

reading group
by Cotter, Colin J 02 Jul '21

02 Jul '21

Dear Firedrakers, Just to remind people that our Firedrake reading group continues on Mondays at 4pm UK time. It's quite a general and very informal reading group on all scientific computing/computational mathematics/numerical analysis topics, just a way to keep us talking together when we can't meet at conferences etc. PhD researchers are especially welcome (and indeed Masters or undergraduate students)! For more details please join #readinggroup on Slack. all the best --cjc Professor Colin Cotter (he/him) Department of Mathematics 755, Huxley Building Imperial College London South Kensington Campus United Kingdom of Great Britain and Northern Ireland +44 2075943468

1 0