Some confusion in the adjustment of free_surface.prm

Hello, experts and scholars.I recently made some changes based on the free_surface.prm file in the manual, and I want to run a model to explore a hot mass of material in the mantle rising up from the mantle and escaping to the surface.
I encountered some confusion in the modification of the model:

  1. I want to add a tidal heating model to my model, but I don’t seem to see the option of tidal heating in the heating model in the manual. If I need to add a tidal heating model, what should I do?

  2. When I increase the radius of the sphere in the free_surface-prm source file, the shape of the sphere becomes less regular. Is there any way to change this?

Hello ddarkerlu,

Thanks for the post! For your first point, you are right that there is currently no tidal heating model in ASPECT, so this would have to be something that you add in a plugin. You can find some general documentation for writing a plugin here [How to write a plugin — ASPECT 2.6.0-pre], and more specific documentation on heating_model plugins in the manual here [Heating models — ASPECT 2.6.0-pre], but don’t hesitate to ask more questions on the forum if you need help with this.

For your second point, from your screenshot it seems to me like this might just be a resolution issue, if you increase the resolution does the sphere become more regular? If increasing the resolution doesn’t help, would you mind posting the changes you made to free_surface.prm here so we can take a closer look?

Cheers,
Daniel

Thank you very much for your reply. After adjusting the resolution, the shape has been greatly improved. Now I will try to write a plug-in to add an addition to the tidal heating model。

1 Like

Dear daniel,
I’m sorry to bother you again. Because I am not very familiar with C++, I wrote a plug-in about tidal heating with GPT, and the following error occurred when compiling, so what should I do to modify the error in tidal_heating.cc file to make it compile properly?

Thank you very much!
Lu


I suspect you are missing a statement like

#include <aspect/material_model/interface.h>

at the top of the file.

Best
W.

Hi Lu,

Yes I also suspect the same thing as Wolfgang. Also, this is something that other ASPECT users might want to use in the future at some point, so it would be great if you opened a pull request to try to get this tidal heating module merged into the main branch (see this page in the manual for some documentation on this process), and as before don’t hesitate to ask more questions on the forum!

Cheers,
Daniel

Dear Wolfgang,
Unfortunately, after adding this, the error remains unchanged, so I wonder if there are some steps or files missing when I wrote this plug-in.
As for the writing of this plugin, at present I have a tidal_heaing.cc file and a CMakeLists.txt modified according to the CMakeLists.txt of the session-2 plugin. Is there anything else I need to do?

Best
Lu

Dear daniel,
Thank you very much for your reply. At present, I am trying to modify my plug-in to make it run normally, but unfortunately, it still has some problems. If in the future it perfectly solves my problem, I’d love to contribute to the ASPECT community!

Best
Lu