underworldcode / underworld2

underworld2: A parallel, particle-in-cell, finite element code for Geodynamics.
http://www.underworldcode.org/
Other
168 stars 58 forks source link

May I ask how to use Underworld2 modeling the melt-related process? #694

Open lsli0116 opened 5 months ago

lsli0116 commented 5 months ago

Hello everyone,

I'm currently exploring the utilization of Underworld2 for simulating melting processes, encompassing melt generation and migration. While I've come across a melting example in UWGeodynamics documentation, which influences the viscosity and density of the swarm, I haven't been able to locate any resources implementing a similar Darcy melt model, akin to the two-phase system described by McKenzie (1984). Additionally, I've noticed that Sim et al. (2020) discuss calculations that can be implemented in Underworld2, but it based on the TerraFERMA implementation here, raising uncertainties regarding compatibility and feasibility within Underworld2.

I'd appreciate any insights or resources regarding existing examples or code for implementing a Darcy melt model in Underworld2. Furthermore, I'm keen to understand the potential avenues and considerations for undertaking such an implementation within the Underworld2 framework.

Thank you for your attention and assistance.

Reference: McKenzie, D., 1984. The Generation and Compaction of Partially Molten Rock. Journal of Petrology 25, 713–765. https://doi.org/10.1093/petrology/25.3.713

Sim, S.J., Spiegelman, M., Stegman, D.R., Wilson, C., 2020. The influence of spreading rate and permeability on melt focusing beneath mid-ocean ridges. Physics of the Earth and Planetary Interiors 304, 106486. https://doi.org/10.1016/j.pepi.2020.106486