Groundwater flow - error at the hydraulic-head distribution

Mailing List Archives ogs-devs
1

Dear all,
I have built a groundwater flow model with OGS. The model itself works, but there is an error displaying it with ParaView. From x = 0 km to x = 2000 km the pressure in the upper layers rises only slightly and then decreases again (see ‘Integrate Variables’). At this point, the vectors (see Glyph) turn round again (see negative values by ‘NODAL_VELOCITY’) and the flux is negative. Actually this should not happen, all arrows are supposed to point from left to right in the positive flux direction.
Can someone tell me the problem, why the pressure is increasing, although it should sink continuously?
Another problem is that I get NAN-messages at x = 0 km, although a fixed pressure head is given as boundary condition at the polyline (x,z = 0).
All files can be open under: https://drive.google.com/drive/folders/0B4AN3Hp76J69NTYzelBqUDl5N28?usp=sharing
Thank you in advance!
Best regards,
Leandra

2

Dear Leandra,

  could you please prepare a mesh with less elements. There is no need for almost 500.000 elements, just for a 2D groundwater model. 10.000 should do the job and it makes visualization and analysis much easier and less time consuming.

Thanks

Thomas

···

https://drive.google.com/drive/folders/0B4AN3Hp76J69NTYzelBqUDl5N28?usp=sharing


-- --------------------------------------------------------------------- Dr. Thomas Kalbacher Department of Environmental Informatics Helmholtz Centre for Environmental Research GmbH - UFZ Helmholtz-Zentrum für Umweltforschung GmbH - UFZ Permoserstraße 15, 04318 Leipzig, Germanyphone +49 (0) 0341 235 1093 / fax +49 (0) 0341 235 1939 / Sitz der Gesellschaft: Leipzig /Registergericht: Amtsgericht Leipzig, Handelsregister Nr. B 4703 / Vorsitzender des Aufsichtsrats: MinDirig Wilfried Kraus / Wissenschaftlicher Geschäftsführer: Prof. Dr. Dr. h.c. Georg Teutsch / Administrativer Geschäftsführer: Prof. Dr. Heike Graßmann --------------------------------------------------------------------- Chairman of the OpenGeoSys Scientific Community: e-mail: ---------------------------------------------------------------------

thomas.kalbacher@ufz.dehttp://www.ufz.dehttp://www.opengeosys.org/info@opengeosys.org

3

Dear Thomas,
thank you for your answer. Unfortunately, a mesh with less elements did not help solving the problem. I upload the current mesh again in the appendix.
Are there any other solutions for my problem?

Best regards,
Leandra

Modell_10_3.msh (1.29 MB)