GlaDS - how to modify geothermal heat flux?

Extension of Elmer in computational glaciology
Post Reply
awhite
Posts: 3
Joined: 30 May 2019, 03:07
Antispam: Yes

GlaDS - how to modify geothermal heat flux?

Post by awhite »

Hi,

I'm setting up a GlaDS model and am trying to see conditions necessary to get channelization under my glacier.
e.g. how strong high geothermal heat flux would needs to be.

How do I adjust geothermal heat flux in GlaDS? Or alternatively the m_b water input term from Equation 12 in Gagliardini '18

From [ELMERICE]/Tests/GlaDS/ i can see the input:

Code: Select all

$Source = 0.0                 ! all input in moulins
I have no moulins, and dont want to add them if possible.

Thanks very much,

Arran
gagliar
Posts: 79
Joined: 04 Sep 2009, 16:34
Location: LGGE - Grenoble
Contact:

Re: GlaDS - how to modify geothermal heat flux?

Post by gagliar »

Dear Arran,

GlaDS doesn't see the geothermal heat flux, but only melt, trough the source term: Hydraulic Potential Volume Source. If you want to include melt from geothermal heat flux (and basal friction and internal deformation) you will have to compute the temperature field and evaluate how much water is produced at the base of the glacier.

Hope it helps,
Regards,
Olivier
awhite
Posts: 3
Joined: 30 May 2019, 03:07
Antispam: Yes

Re: GlaDS - how to modify geothermal heat flux?

Post by awhite »

Dear Olivier,

Thanks for your reply,

What does 'Hydraulic Potential Volume Source' correspond to in your '18 paper with Werder?

Is it m_b (equation 12) from the paper?

Thanks very much,

arran
gagliar
Posts: 79
Joined: 04 Sep 2009, 16:34
Location: LGGE - Grenoble
Contact:

Re: GlaDS - how to modify geothermal heat flux?

Post by gagliar »

Yes, this is m_b from Eq. (12). So in fact, contrary to what I was saying in my previous post, m_b can be directly evaluated as a function of the geothermal heat flux and the frictional heating.
awhite
Posts: 3
Joined: 30 May 2019, 03:07
Antispam: Yes

Re: GlaDS - how to modify geothermal heat flux?

Post by awhite »

Okay great, thanks for that.
On [ELMERICE]/Tests/GlaDS/shmip_B5.sif it says that this source is 'moulin input' - which needs updating.

I am also trying to run a drainage model with a poroelastic sheet where sheet thickness is directly related to effective pressure, h = h(N) --- no opening or closing terms and no dependence on basal velocity.
e.g. equation 5b in hewitt '13 https://doi.org/10.1016/j.epsl.2013.04.022

Is this possible with Elmer/GlaDs?
Id just need to switch off the cavity solver and calculate h_cav as a function of N.

Cheers
arran
gagliar
Posts: 79
Joined: 04 Sep 2009, 16:34
Location: LGGE - Grenoble
Contact:

Re: GlaDS - how to modify geothermal heat flux?

Post by gagliar »

It is not possible to switch off the cavity (just the channels). I guess you will not be able to do what you want directly with the GlaDS solver but you will have to write your own solver with the appropriate h=h(N) function implemented.
Regards,
Olivier
Post Reply