| Author |
 Topic  |
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
 Posted - 02/11/2010 : 17:56:06
|
Hello,
I am working with Hydrus 3D. My work is to elaborate a 3D model for an aquifer and study the recharge process in this aquifer.
My first step in this project is to edit the geometry of this aquifer. In fact, I have the DEM of the hole basin. I have a problem in editing the 3D geometry in which we will simulate the water flow.
I have read that we should define thikness vectors after editing the base surface. I can't see how to integrate the DEM (which represents the top surface) in this geometry.
thanks for helping me |
|
|
Jirka
Moderator
USA
1788 Posts |
Posted - 02/11/2010 : 18:27:31
|
Mirek will probably answer this later and better than I can. But anyway, in the meantime look at the tutorial 2.07 (3D Domain, Solid 6 - import of surface points from a GIS (in text format)). This will probably answer most of your questions.
Jirka
http://www.pc-progress.com/en/Default.aspx?h3d-tutorials
|
 |
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 02/11/2010 : 18:44:19
|
thanks for replying, I will focus in the tutorial and I will wait also for Mirek's reply
Thanks a lot |
|
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 02/11/2010 : 19:27:12
|
Hi,
If I had undrestand the tutorial 2.07, the file imported from the GIS contains the hole geometry of the aquifer studied, I mean that it includes the thikness of the aquifer. However in my case, the DEM contains the coordinates (x,y,z) of the surface, it doesn't includes the thikness of the aquifer. I have to complete these informations by editing the 3D geometry. am I right???
Regards |
|
|
|
Mirek
Moderator
Czech Republic
194 Posts |
Posted - 02/11/2010 : 21:17:29
|
Hi,
The simplest way to create a domain from your 3D points is described below.
Regards Mirek
1/ Read HYDRUS Help: Graphical User Interface -> Menu -> File -> Import and Export ->Import Geometry from a Text File, KEY_WORD: THICKNESS_ARR3Z_NLAYERS
2/ I assume that you have an array of 3D points representing the top surface of your domain and that the bottom surface is a plane with z-coordinate = Z0
3/ Create a text file like this:
OBJECT=THICKNESS_ARR3Z_NLAYERS
4.500000e+000 4.000000e+000 0.000000e+000 0.000000e+000 1.300000e+000
5.500000e+000 2.500000e+000 0.000000e+000 0.000000e+000 1.300000e+000
7.000000e+000 4.500000e+000 0.000000e+000 0.000000e+000 1.400000e+000
7.000000e+000 2.500000e+000 0.000000e+000 0.000000e+000 1.200000e+000
6.500000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.600000e+000
5.000000e+000 2.000000e+000 0.000000e+000 0.000000e+000 1.200000e+000
4.500000e+000 2.000000e+000 0.000000e+000 0.000000e+000 1.400000e+000
5.000000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.500000e+000
5.500000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.100000e+000
;
where the number of rows (coordinates) = number of your 3D points and where the columns are:
Column1 = X–coordinates of your points
Column2 = Y–coordinates of your points
Column3 = Z0
Column4 = Z0
Column5 = Z–coordinates of your points
4/ Run “File -> Import and Export ->Import Geometry from a Text File” and import data from the text file. This will import thickness vectors that will define your 3D domain.
5/ Create the 3D domain similarly as shown in tutorial 2.07. |
|
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 02/12/2010 : 10:54:02
|
Hi,
Thank you Mirek, I will focus on your recommendations.
Regards |
|
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 02/16/2010 : 09:29:42
|
Hi,
I had imported my 3D geometry successfully. But, while creating the 3D domain similarly as shown in tutorial 2.07, I received a warninig ("Points 28411 and 28498 ovelap")when trying to create surface via boundary, I can't undrestand the problem.
PS: my DEM contains about 9500 points with their correspondant thikness points. It's too hard to join the points of the surface base boundary manually. Can we have other alternative to import or create this boundary???
Regards. |
|
|
|
Mirek
Moderator
Czech Republic
194 Posts |
Posted - 02/16/2010 : 10:40:16
|
Hi,
The message "Points 28411 and 28498 overlap" means that these points have identical coordinates, which is not allowed. I think that you probably created a duplicated point when you were inserting boundary curves. You can find the point using the „Find“ function (menu -> Edit -> Find) and fix the problem manually or you can fix all duplicated points automatically using the “Repair geometry” function (menu Tools -> Repair Geometry).
You can import all objects (including boundary lines and surfaces) from the text file. You will find all needed commands (KEY_WORDS) in HYDRUS Help: Graphical User Interface -> Menu -> File -> Import and Export ->Import Geometry from a Text File.
However, I believe that graphical definition of the domain should be possible even if the number of imported points is relatively large. Could you send me the file with your coordinates (mireks@pc-progress.cz)? I will try to work with this project and perhaps I will be able to give you some hints.
Mirek |
|
|
|
Mirek
Moderator
Czech Republic
194 Posts |
Posted - 02/18/2010 : 11:30:38
|
Hello Maroua,
Thank you for the file with coordinates. The main problem is that your points are given in some global coordinates (e.g. x=587012.594 m , y=397585.9776 m), which is not good. The diameter of your domain is about 500 m (which is OK) but the domain is located far away from the origin [0,0,0]. Its center is at point x=586853 m, y=397802 m. HYDRUS was able to import this domain but the incorrect position of the domain may cause many problems and it is necessary to move the domain close to the origin [0,0,0]. I modified X,Y coordinates in the input file so that the domain was moved to [0,0,0] and then it was relatively easy to import thickness vectors, create boundary lines and the solid. You can see the result on attached figures. Then I generated FE-mesh with 30 horizontal layers (refined close to the top surface). The mesh is rather large (it has over 500 000 elements) and reducing the number of elements is almost impossible because of the very large number of imported geometrical points (all these points must be included in the FE-mesh). To work with the model effectively, you either need a very fast PC or you could reduce the number of geometrical points in the input file, which would allow you reducing the FE-mesh size.
I was able to prepare the model and generate the FE-mesh within one hour. I will send you the project (without the FE-mesh but ready for mesh generation). However, I can imagine that for some less-experienced HYDRUS users this can take a longer time. I could prepare a video tutorial showing several “tricks” that may be useful when working with domains like this one. Let me know if I may use your domain for this tutorial. We will not publish the coordinates (which may be protected), just the video and images displaying the domain.
Regards
Mirek
|
Edited by - Mirek on 02/18/2010 11:31:22 |
|
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 03/04/2010 : 20:18:46
|
Hi Mirek,
I still have one more question about the solid geometry. Can we consider a variable base surface. If we are working on modelling aquifers, we needs to consider a base surface similar to the topography of the top surface. I think that this option is not integrated in the 1.10 version of Hydrus. am I right?
Regards
|
|
|
|
Mirek
Moderator
Czech Republic
194 Posts |
Posted - 03/04/2010 : 20:48:08
|
Hello Maroua,
The base surface must be planar but bottom points of thickness vectors need not lie in the base surface (see Hydrus help – >Using Hydrus -> Domain Geometry -> Geometric Objects -> Thicknessses -> Thickness Vector). This means that the bottom surface of your domain can be variable as well as the top surface (which is what you need). Z-coordinates of the bottom surface should be defined in the fourth column of the THICKNESS_ARR3Z_NLAYERS block:
OBJECT=THICKNESS_ARR3Z_NLAYERS
4.500000e+000 4.000000e+000 0.000000e+000 0.000000e+000 1.300000e+000
5.500000e+000 2.500000e+000 0.000000e+000 0.000000e+000 1.300000e+000
7.000000e+000 4.500000e+000 0.000000e+000 0.000000e+000 1.400000e+000
7.000000e+000 2.500000e+000 0.000000e+000 0.000000e+000 1.200000e+000
6.500000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.600000e+000
5.000000e+000 2.000000e+000 0.000000e+000 0.000000e+000 1.200000e+000
4.500000e+000 2.000000e+000 0.000000e+000 0.000000e+000 1.400000e+000
5.000000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.500000e+000
5.500000e+000 1.000000e+000 0.000000e+000 0.000000e+000 1.100000e+000
;
Regards
Mirek
Btw: sorry for the delay of the promised video tutorial. Now I’m very busy but I’m going to do it later. |
|
|
|
MOD_LAB
Average Member
Tunisia
14 Posts |
Posted - 03/05/2010 : 20:47:38
|
Hello Mirek,
Thanks a lot for replying.
ok, no problem about the delay, take your time.
Regards
Maroua |
|
|
| |
Topic |
|
Geometry editing in Hydrus 3D
