Closed JHopeCollins closed 1 year ago
JHopeCollins
commented
1 year ago With the current resolution (12x6x12 mesh elements) a window size of 4 takes up around 45GB, so should be runable on a workstation.
I think the large footprint is due to the AssembledPC on the blocks.
JHopeCollins
commented
1 year ago I started out with a 2D bump (just the 2D example with multiple spanwise elements, but then the solution is also 2D. I changed it to 3D to check that it was actually doing something reasonable in the extra dimension.
I could add a switch for 2D / 3D geometry?
colinjcotter
commented
1 year ago Yes - I didn't describe it very well. I meant a bump that depends on x and y. Don't need it for this PR though.
From: JHopeCollins @.***> Sent: 24 February 2023 11:51 To: firedrakeproject/asQ Cc: Cotter, Colin J; Review requested Subject: Re: [firedrakeproject/asQ] Compressible 3D example (PR #108)
I started out with a 2D bump (just the 2D example with multiple spanwise elements, but then the solution is also 2D. I changed it to 3D to check that it was actually doing something reasonable in the extra dimension.
— Reply to this email directly, view it on GitHubhttps://github.com/firedrakeproject/asQ/pull/108#issuecomment-1443577200, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABOSV4UAD4GYQZ55J4W3Z43WZCOFZANCNFSM6AAAAAAVEN76LQ. You are receiving this because your review was requested.Message ID: @.***>
JHopeCollins
commented
1 year ago The bump is Gaussian in both x and y, is that what you meant?
# making a mountain out of a molehill
a = 2000.
xc = Lx/2.
yc = Ly/2.
x, y, z = fd.SpatialCoordinate(mesh)
hm = 1.
r2 = ((x - xc)/a)**2 + ((y - yc)/(2*a))**2
zs = hm*fd.exp(-r2)ah so it is! I was looking at the - not the +!
From: JHopeCollins @.***> Sent: 24 February 2023 11:55 To: firedrakeproject/asQ Cc: Cotter, Colin J; Review requested Subject: Re: [firedrakeproject/asQ] Compressible 3D example (PR #108)
The bump is Gaussian in both x and y, is that what you meant?
a = 2000. xc = Lx/2. yc = Ly/2. x, y, z = fd.SpatialCoordinate(mesh) hm = 1. r2 = ((x - xc)/a)*2 + ((y - yc)/(2a))*2 zs = hmfd.exp(-r2)
— Reply to this email directly, view it on GitHubhttps://github.com/firedrakeproject/asQ/pull/108#issuecomment-1443580849, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABOSV4X6WMAU7UDB53AWATLWZCOSTANCNFSM6AAAAAAVEN76LQ. You are receiving this because your review was requested.Message ID: @.***>
Example script for the 3D compressible Euler equations with flow over a gaussian hill in a doubly periodic domain.. Adapted from the 2D vertical slice test case of flow over a mountain.