underworldcode / underworld2

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

A question about strata results from coupled model #689

Closed tyszwh closed 4 months ago

tyszwh commented 5 months ago

Hi Dr Julian I have some question about the coupled model, and I'm not sure it's a problem with the code. I was hoping you could give me some advice. When I use the coupled uwgeodynamics (2d) and badlands (2.5d) model, the position of the basin sequences formed by the surface process model seems to be out of sync with the position of the basin sediment on the underworld2 model, there is clearly some deviation. How is this caused? Is this prossible to fix this problem? Please refer to the attached figure for details.

Thanks

Model_Resutls

julesghub commented 5 months ago

Hi @tyszwh, how are you visualising these results? what version of the code are you running?

The mesh for badlands and the upper surface of Underworld2 are not the same and information is shared between the two meshes using a virtual interface. I believe this is the cause of the deviation you are observing.

tyszwh commented 5 months ago

These results were visualised using paraview and zoomed in 5x on the vertical. I have also used python plot strata results and the results are similar. I am using underworld version 2.16.x and badlands version 2.2.4. Is the 'np_surface' defined in the code a virtual surface? Any suggestions on how to fix this?

tyszwh commented 5 months ago

Hi,@julesghub I did some more testing and I found that even though the two models share the same velocity field, they don't seem to deform in the same way in the areas of greater deformation. Is this related to mesh accuracy? Any suggestions on fixing this issue? Thanks