[firedrake] Second derivatives
Justin Chang
jychang48 at gmail.com
Tue Jul 21 20:17:39 BST 2015
I am looking at equal order mixed formulations. Two discretizations in
particular:
1) Our group recently released a paper on enforcing local mass balance and
discrete maximum principles through the least-squares finite element
method. The paper can found here:
http://arxiv.org/pdf/1506.06099v1.pdf
Equation 4.13 is what I want to solve using firedrake. Everything in that
paper was written in MATLAB because the authors could not get FEniCS to do
what they want. I am guessing it's because there were no Q1 elements
available within FEniCS (double derivative of P1 elements makes no sense).
2) The other discretization, based on the variational multi-scale
formulation, for Darcy-Brinkmann equations is described in this paper:
http://onlinelibrary.wiley.com/doi/10.1002/fld.2544/abstract
Though I am currently looking at the steady-state version of this. Appendix
A1 describes how they discretize the \delta v and \delta w terms.
I could perhaps figure these out through trial and error, but I got my
hands full with other things ATM :)
Thanks,
Justin
On Tue, Jul 21, 2015 at 4:33 AM, David Ham <David.Ham at imperial.ac.uk> wrote:
> Hi Justin,
>
> Can you provide a little more information about the sort of discretisation
> you're talking about: if we can see what you're talking about then we'd be
> in a better position to tell you whether Firedrake can do that.
>
> Cheers,
>
> David
>
> On Fri, 17 Jul 2015 at 07:39 McRae, Andrew <a.mcrae12 at imperial.ac.uk>
> wrote:
>
>> It's possible to represent second derivatives in a form, such as
>> assemble(div(grad(f))*dx). I assume this is true in FEniCS as well as
>> Firedrake, because this is just UFL/FFC/FIAT functionality. This would
>> produce the sum (over cells) of the Laplacian of f on each cell.
>>
>> However, like FEniCS, none of our function spaces have more than C^0
>> continuity. That is, the functions are at most continuous, but won't have
>> continuous derivatives. It's likely that you won't want to use second
>> derivatives of C^0 functions in a practical discretisation (though I'm sure
>> there are methods that *do* do this).
>>
>> Andrew
>>
>> On 17 July 2015 at 01:50, Justin Chang <jychang48 at gmail.com> wrote:
>>
>>> Hi everyone,
>>>
>>> Is it possible to do second derivatives? What I mean by that is things
>>> like div[grad[u]] and grad[grad[u]]. I haven't tried this out yet, but we
>>> use these discretizations a lot for our research, and FEniC's inability to
>>> do this made us sad (although our more major qualms had to do with its
>>> inability to support quads).
>>>
>>> Thanks,
>>> Justin
>>>
>>>
>>
> _______________________________________________
> firedrake mailing list
> firedrake at imperial.ac.uk
> https://mailman.ic.ac.uk/mailman/listinfo/firedrake
>
>
-------------- next part --------------
HTML attachment scrubbed and removed
More information about the firedrake
mailing list