Two-Phase flow (save interface position)

Post processing utility for Elmer
Post Reply
emarchandise
Posts: 9
Joined: 23 Jan 2011, 16:45

Two-Phase flow (save interface position)

Post by emarchandise »

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
Attachments
dambreak.sif
Case dambreak
(4.12 KiB) Downloaded 455 times
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: Two-Phase flow (save interface position)

Post by raback »

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
emarchandise
Posts: 9
Joined: 23 Jan 2011, 16:45

Re: Two-Phase flow (save interface position)

Post by emarchandise »

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
Attachments
dambreak.geo
(660 Bytes) Downloaded 421 times
Post Reply