Linking subglacial hydrology to Heat Equation

Extension of Elmer in computational glaciology

Linking subglacial hydrology to Heat Equation

Postby joeatodd » 20 May 2013, 14:28

Hi all,

I'm seeking some advice on linking the residual of the limited heat equation to subglacial hydrology at the bed.

My model setup has a Dirichlet temperature condition at the bed (Pressure Melting Point) and, if I understand correctly, the heat sinks and sources imposed in satisfying this Dirichlet condition should be directly related to the melting/refreezing at the bed by the value of the latent heat of fusion.

From reading the documentation, I understand that

Code: Select all
Calculate Loads = Logical True

can be set, which makes available the "Temperature Loads" variable.

My plan would be then to have the source term set such that:

Code: Select all

$LHF = latentheatoffusion

Water Flux = Variable Temperature Loads
   Real MATC "tx/$LHF"

Can anyone advise me if this will achieve the desired result or if I've missed something?


Posts: 34
Joined: 02 Feb 2012, 18:49

Re: Linking subglacial hydrology to Heat Equation

Postby raback » 20 May 2013, 14:43

Hi Joe

The nodal loads are in Units of power, e.g. Watts. So it would be natual to use the them directly as loads in the coupled equation as well. So this "Water Flux" should then be directly a nodal quantity and you would add it directly to the r.h.s. vector. This is actually automatically available as a library function. Assuming that you would apply it to the N-S equation you could then set the source term for continuity equation as

Flow Solution 4 Load = Variable Temperature Loads
Real MATC "tx/Cmelt"

Or something like that...

In order to couple equations on a discrete level the mesh must of course be the same. If the method is applicable then it is often the preferred method.

Site Admin
Posts: 3110
Joined: 22 Aug 2009, 11:57
Location: Espoo, Finland

Return to Elmer/Ice

Who is online

Users browsing this forum: No registered users and 1 guest