The boundary conditions are the temperatures at the bottom of the ice layer and on top of the CO2. I expected similar values of heat flow on the surface of the CO2 layer, but results indicate some points with large variations. This is a screenshot of the solution.
The geometry and the mesh look fine, but I think that the problem could be on the thin layer. I have tried to refine this part with the ‘region refinement’ and the ‘local element size’ tools, but results don’t improve and I run out of memory easily. I would appreciate any help with this problem.
Hi Isabel! The mesh here is really coarse, I would recommend trying another setup, and make sure the thin layer has at least two elements across the z-direction For this, you can try a region refinement applied on the specific part.
I have used the region refinement tool with the default values, but I think that the mesh is still coarse and the simulation is taken too much. Is this normal?
I have tried to change the parameters to improve the mesh manually, but I run out of memory easily. I have also tried with the local element size option and I have the same problem. Maybe the extension of my model is too large?
Yes with your model you have too much domain to be simulated easily with any level of detail.
From your initial post you have 25000 m x 25000 m x 1700 m volume and you are trying to evaluate a thin sheet that consists of 15m of that 1700m. I would highly suggest trimming the area down from 25km x 25km to something maybe 1/10th of the size. You may also find you don’t even need the 1700m of height on your model and can reduce that, but I will leave that for you to discover! Best of luck. This may take a few iterations of what should be the optimal size but you will be able to do generate a more refined mesh without running out of memory.
Let me know what you think.