Hello,
I am running successfully two-phase flow computations with elmer (such ad dam break prolems, sloshing) for moderate ratios of density.
In order to compare with reference solutions, I would like to know if it is possible to save the position of the interface (zero value of the levelet) with Elmer along a line.
The saveLine does almost the job , but I would like to know if it is possible to save only the position of the isozero along this line.
Solver 6
Equation = "SaveLine"
Procedure = File "SaveData" "SaveLine"
Filename = "line.dat"
File Append = Logical True
Variable 1 = String Levelset
Variable 2 = String Coordinate 2
Polyline Coordinates(2,2) = Real 0.0 0.0 0.0 2.0 ! (n,DIM) Real X1 Y1 X2 Y2
End
Two-Phase flow (save interface position)
-
- Posts: 9
- Joined: 23 Jan 2011, 16:45
Two-Phase flow (save interface position)
- Attachments
-
- dambreak.sif
- Case dambreak
- (4.12 KiB) Downloaded 455 times
-
- Site Admin
- Posts: 4812
- Joined: 22 Aug 2009, 11:57
- Antispam: Yes
- Location: Espoo, Finland
- Contact:
Re: Two-Phase flow (save interface position)
Hi
You could try the IsosurfaceSolver that creates a surface and let ResultOutputSolver save it. See an example I just sent:
viewtopic.php?f=18&t=1757
-Peter
You could try the IsosurfaceSolver that creates a surface and let ResultOutputSolver save it. See an example I just sent:
viewtopic.php?f=18&t=1757
-Peter
-
- Posts: 9
- Joined: 23 Jan 2011, 16:45
Re: Two-Phase flow (save interface position)
Thanks, I tried to use it but the plugin crashes it does not work.
I used gmsh as input for the mesh. Here is my geo file.
At line 808 of file Isosurface.f90
Traceback: not available
Fortran runtime error: Deallocated a bad pointer
which appear here in the fortran code :
DefElement => AllocateElement()
DefElement % TYPE => GetElementType(NewElemType)
END IF
END DO
DEALLOCATE(DefElement) !! error here
END FUNCTION CreateSurfaces
I also noticed in thi scase that elmer has some troubles with the postpro of the boundary elements when writing the boundary elements (the post file in gmsh did not contain the wright number of #elements).
In order to postprocess with gmsh I had to change the following two variables in ResultOutputSolve.src at line 1502
>>NumberOfAllElements = Model % NumberOfBulkElements
>> nsize = NumberOfAllElements
Thanks in advance,
Emilie
I used gmsh as input for the mesh. Here is my geo file.
At line 808 of file Isosurface.f90
Traceback: not available
Fortran runtime error: Deallocated a bad pointer
which appear here in the fortran code :
DefElement => AllocateElement()
DefElement % TYPE => GetElementType(NewElemType)
END IF
END DO
DEALLOCATE(DefElement) !! error here
END FUNCTION CreateSurfaces
I also noticed in thi scase that elmer has some troubles with the postpro of the boundary elements when writing the boundary elements (the post file in gmsh did not contain the wright number of #elements).
In order to postprocess with gmsh I had to change the following two variables in ResultOutputSolve.src at line 1502
>>NumberOfAllElements = Model % NumberOfBulkElements
>> nsize = NumberOfAllElements
Thanks in advance,
Emilie
- Attachments
-
- dambreak.geo
- (660 Bytes) Downloaded 421 times