[firedrake] building a function space for slice vector fields

Cotter, Colin J colin.cotter at imperial.ac.uk
Mon Jun 15 18:07:43 BST 2015


I'm not quite sure what I need to try, exactly?

--cjc
________________________________
From: firedrake-bounces at imperial.ac.uk [firedrake-bounces at imperial.ac.uk] on behalf of David Ham [David.Ham at imperial.ac.uk]
Sent: 15 June 2015 18:04
To: firedrake
Subject: Re: [firedrake] building a function space for slice vector fields

Have you tried it? I think it might Just Work. The first thing which might break is the HDiv operator, but you never know....
On Mon, 15 Jun 2015 at 11:58 Cotter, Colin J <colin.cotter at imperial.ac.uk<mailto:colin.cotter at imperial.ac.uk>> wrote:
Dear all,
  A lot of what we do in the vertical slice model would be a lot easier if we could combine together our u-w components of our velocity field (stored as an RT field in the slice) with our v component (stored as a DG field in the slice) so that the basis functions are 3D vectors on a 2D x-z plane. In particular, then I could assemble the inverse of the combined u-w-v velocity mass matrix combined with Coriolis, and make things a lot more implicit in the Coriolis term, leading to bigger timesteps when we sent the Coriolis parameter to infinity (something we have to do in our test case unfortunately).

Is this supported by current code, if not, how much work would it take for us to have this?

all the best
--cjc
-------------- next part --------------
A non-text attachment was scrubbed...
Name: winmail.dat
Type: application/ms-tnef
Size: 5330 bytes
Desc: not available
URL: <http://mailman.ic.ac.uk/pipermail/firedrake/attachments/20150615/161a279f/attachment.bin>


More information about the firedrake mailing list