Groundwater flow - error at the hydraulic-head distribution

[Leandra Drewes]

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

[Thomas Kalbacher]

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

--
You received this message because you are subscribed to the Google Groups "ogs-devs" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ogs-devs+u...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

-- 

---------------------------------------------------------------------
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
thomas.k...@ufz.de / http://www.ufz.de

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: http://www.opengeosys.org/
e-mail: in...@opengeosys.org
--------------------------------------------------------------------- 

[Leandra Drewes]

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