Re: [firedrake] upwind example remark; extension to CN in time
Onno did attach rectmesh.msh... Onno, is there a reason why you didn't use one of the built-in meshes? cheers --cjc On 25 June 2016 at 16:18, Francis Poulin <fpoulin@uwaterloo.ca> wrote:
Hello Onno,
I just tried running your code but see that I can't because I don't have the rectmesh.msh. Might you be able to send that so we can see exactly what you see?
Francis
P.S. I would also be curious to see the gmsh file that generated the msh file if possible.
------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo
email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637
------------------------------ *From:* firedrake-bounces@imperial.ac.uk [firedrake-bounces@imperial.ac.uk] on behalf of Onno Bokhove [O.Bokhove@leeds.ac.uk] *Sent:* Saturday, June 25, 2016 11:05 AM *To:* firedrake *Subject:* [firedrake] upwind example remark; extension to CN in time
Hi FDs,
(i) The upwind example is odd in that the real boundaries are treated
in a special way rather than setting the outside value in
the same numerical flux as anywhere else;
my question is how this can be fixed?
So no special terms at the in- and outflow, just specify
the outside velocity used. This is imo the proper DG way not was is done
now.
(ii) Attached my attempt to do a time dependent advection
with Crank Nicolson, which works but does not yet
give the right answer, in part due to (i).
[I don't understand why I have to make
my theta variable that a constant but other wise it did not work.]
Any suggestions, corrections, also on other silly things done by me?
Next step is to extend to system with 3 equations; hence the few extra but unused variables.
Thanks, from midnight Seoul in old style Korean guesthouse.
Onno ------------------------------
-- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916
Thanks Colin for pointing that out. Sorry that I didn't see the file but for some strange reason Outlook will not tell me there's a .msh file. I found it through another means though. I just ran it and will try and better understand it. Cheers, Francis ------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo email: fpoulin@uwaterloo.ca Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637 ________________________________ From: firedrake-bounces@imperial.ac.uk [firedrake-bounces@imperial.ac.uk] on behalf of Colin Cotter [colin.cotter@imperial.ac.uk] Sent: Saturday, June 25, 2016 4:22 PM To: firedrake Subject: Re: [firedrake] upwind example remark; extension to CN in time Onno did attach rectmesh.msh... Onno, is there a reason why you didn't use one of the built-in meshes? cheers --cjc On 25 June 2016 at 16:18, Francis Poulin <fpoulin@uwaterloo.ca<mailto:fpoulin@uwaterloo.ca>> wrote: Hello Onno, I just tried running your code but see that I can't because I don't have the rectmesh.msh. Might you be able to send that so we can see exactly what you see? Francis P.S. I would also be curious to see the gmsh file that generated the msh file if possible. ------------------ Francis Poulin Associate Professor Department of Applied Mathematics University of Waterloo email: fpoulin@uwaterloo.ca<mailto:fpoulin@uwaterloo.ca> Web: https://uwaterloo.ca/poulin-research-group/ Telephone: +1 519 888 4567 x32637<tel:%2B1%20519%20888%204567%20x32637> ________________________________ From: firedrake-bounces@imperial.ac.uk<mailto:firedrake-bounces@imperial.ac.uk> [firedrake-bounces@imperial.ac.uk<mailto:firedrake-bounces@imperial.ac.uk>] on behalf of Onno Bokhove [O.Bokhove@leeds.ac.uk<mailto:O.Bokhove@leeds.ac.uk>] Sent: Saturday, June 25, 2016 11:05 AM To: firedrake Subject: [firedrake] upwind example remark; extension to CN in time Hi FDs, (i) The upwind example is odd in that the real boundaries are treated in a special way rather than setting the outside value in the same numerical flux as anywhere else; my question is how this can be fixed? So no special terms at the in- and outflow, just specify the outside velocity used. This is imo the proper DG way not was is done now. (ii) Attached my attempt to do a time dependent advection with Crank Nicolson, which works but does not yet give the right answer, in part due to (i). [I don't understand why I have to make my theta variable that a constant but other wise it did not work.] Any suggestions, corrections, also on other silly things done by me? Next step is to extend to system with 3 equations; hence the few extra but unused variables. Thanks, from midnight Seoul in old style Korean guesthouse. Onno ________________________________ -- http://www.imperial.ac.uk/people/colin.cotter www.cambridge.org/9781107663916<http://www.cambridge.org/9781107663916> [http://assets.cambridge.org/97811076/63916/cover/9781107663916.jpg]
participants (2)
-
Colin Cotter
-
Francis Poulin