Enthalpy Thermomechanical Coupling of "Test Glacier"
Posted: 02 Dec 2020, 19:08
Hello,
I've just wrapped up the beginners course over the last two weeks. Thanks again for putting the course on, it was great!
As a jump start on some research objectives I’ve been trying to rewrite the thermomechanically coupled diagnostic run from the testglacier_flowline module with an enthalpy formulation. I’ve been working from the Stokes_diagnostic_thermomech.sif file from the course folder and have changed Solver 4 to solve the “Enthalpy Equation” instead of “Homologous Temperature Equation.” Following the enthalpy solver documentation on the Wiki, the enthalpy test on the github repo, and some slides on the Wiki I’ve:
I was hoping someone could provide some guidance as to where I am might be going wrong. I've attached my most up to date .sif file for those that are interested.
Thanks in advance,
Andrew
I've just wrapped up the beginners course over the last two weeks. Thanks again for putting the course on, it was great!
As a jump start on some research objectives I’ve been trying to rewrite the thermomechanically coupled diagnostic run from the testglacier_flowline module with an enthalpy formulation. I’ve been working from the Stokes_diagnostic_thermomech.sif file from the course folder and have changed Solver 4 to solve the “Enthalpy Equation” instead of “Homologous Temperature Equation.” Following the enthalpy solver documentation on the Wiki, the enthalpy test on the github repo, and some slides on the Wiki I’ve:
- added the necessary enthalpy related constants to the Constants section
- changed the material section to a “power law” Viscosity Model from “Glen”. Not totally sure about this, but all the examples I’ve come across have done this. Is this necessary because the built in “Glen” Viscosity Model is formulated as a function of temperature and not enthalpy?
- added enthalpy related parameters to the Material section (i.e. Enthalpy Density, Enthalpy Heat Diffusivity , and Enthalpy Water Diffusivity ).
Changed the basal boundary condition to be a Enthalpy Heat Flux, all of which are in consitent units (MPa - m - a, at least I think). - Changed the surface boundary condition to be a prescribed Enthalpy_h which I’ve assumed is the same prescribed temperature as before now converted to enthalpy. I got the formula from integrating the top condition of Equation 14 from Gilber et al 2014.
I was hoping someone could provide some guidance as to where I am might be going wrong. I've attached my most up to date .sif file for those that are interested.
Thanks in advance,
Andrew