Hi all,
I am working on computing the melt fraction and total melt production for no-melt material models like the Steinberger model. Specifically, I aim to compute the melt fraction at each time step at each evaluation point based on P and T, and then integrate the melt volume over time. My model assumes that melts are not extracted and do not affect the model’s behavior, similar to the approach described inLi et al. (2016)
My current thought is to creat an additional generic compositional field called “melt_fractions” that could be updated at each time step by the material model. This way, I could use the existing melt statistic post-processor to compute the total melt volume.
Does this approach sound reasonable? If so, I would appreciate any guidance on which source files I should modify or add to implement this. Any thoughts or suggestions would be greatly appreciated!
Best,
Lucy
Hi Lucy,
Sorry for the late reply.
Just to make sure I understand this correctly:
When you say you want to integrate the melt volume over time, what you want is to compute–at each point in time–is the maximum melt fraction each parcel of rock has experienced at any point in time before, and then integrate that over the domain? In other words, if a part of the model is partially molten in one step, and then it’s still molten in the next, you do not want to count it twice? It’s like we pretend the melt is being extracted (we just do not let that affect the model) and then want to know the integrated “depletion” in each step? So in the end, we get the total amount of melt that has been produced in the model?
In this case, I think the melt statistics postprocessor is not useful, because it gives you the melt fraction at each point in time, so you can not track if any point in the model has experienced melting before. Instead, I think what you suggested with the additional compositional field makes sense. Let’s say we call the field “depletion”. You can implement a reaction term that compares the current melt fraction to whatever is stored in the depletion field, and if the current melt fraction is larger, add the difference to the depletion field. That is basically what the melt material models already do, so you can have a look at them. You can get the melt fraction from the “katz2003_model” reaction plugin as in melt_simple (see here), or compute your own melt fraction.
This means you only need to modify the material model you are using. You can then simply use the composition statistics postprocessor to get the total depletion (= total melt production).
Does that help?
-Juliane