Dear Patrick, Thanks for your help! It essentially worked, however, for some reason the nonlinear solver does not allow me to use direct methods as the inner solver, i.e., ('ksp_type': 'preonly', 'pc_type': 'lu', 'pc_factor_mat_solver_type': 'mumps') produces some sort of error, possibly related to local_to_global_map, I do not understand the error actually. Thanks anyway! On Tue, Apr 7, 2020 at 8:00 AM <firedrake-request@imperial.ac.uk> wrote:

Send firedrake mailing list submissions to firedrake@imperial.ac.uk

To subscribe or unsubscribe via the World Wide Web, visit https://mailman.ic.ac.uk/mailman/listinfo/firedrake or, via email, send a message with subject or body 'help' to firedrake-request@imperial.ac.uk

You can reach the person managing the list at firedrake-owner@imperial.ac.uk

When replying, please edit your Subject line so it is more specific than "Re: Contents of firedrake digest..."

Today's Topics:

1. BC's - 1D Hermite elements (Roberto Federico Ausas) 2. Re: BC's - 1D Hermite elements (Patrick Farrell) 3. Problem installing Firedrake (Smail Merabet)

----------------------------------------------------------------------

Message: 1 Date: Mon, 6 Apr 2020 16:49:25 -0300 From: Roberto Federico Ausas <rfausas@icmc.usp.br> To: firedrake@imperial.ac.uk Subject: [firedrake] BC's - 1D Hermite elements Message-ID: <CAH6Ec4nWnhmM=HeMtd8uU9VQ9NF0_hVdorMFdK= p8Xc3Gus4xQ@mail.gmail.com> Content-Type: text/plain; charset="utf-8"

Hello,

My name is Roberto. I'm currently a professor in the Applied Mathematics Department of Univ. of S?o Paulo at S?o Carlos.

I'm working on some fluid-structure interaction problems in which the solid will be modeled as a 1D beam. Our reference integration domain will be the interval [0,L].

We are planning to use firedrake and Hermite elements. The problem is nonlinear.

I already have a working example. Values degrees of freedom are easy to impose by using DirichletBC(), however I'm not figuring out how to enforce derivatives degrees of freedom as boundary conditions (let say at x=0), which I need for a clamped beam.

Any help or trick would be appreciated Regards

-- Roberto F. Ausas Instituto de Ci?ncias Matem?ticas e de Computa??o - ICMC/USP S?o Carlos, SP - Brasil *http://www.lmacc.icmc.usp.br/~ausas* <http://www.lmacc.icmc.usp.br/~ausas

...