About how-to set surface tension in Level Set Method

Numerical methods and mathematical models of Elmer
Post Reply
SagittariusChiron
Posts: 3
Joined: 11 Jan 2021, 05:36
Antispam: Yes

About how-to set surface tension in Level Set Method

Post by SagittariusChiron »

Hellow.

I am a beginer of the Elmer in Japan and started the Elmer only one month ago.

I am interested in a behavior of a free surface, which is in the condition of steady flow and has surface tension coefficient.

I have found two relevant discussion in this post:

viewtopic.php?f=3&t=3303&p=10533&hilit= ... ion#p10533

viewtopic.php?f=3&t=808&p=10690&hilit=s ... ion#p10690

I read those,and I undestand that the "surface tension coefficient" means the "curvature coefficient " in levelsetcurvature solver.

Now, I already done some tutorials.
"Axisymmetric coating prosess" as 10th tutorial of Elmer non-GUI Tutorials (August 3,2020)is useful for me. But, this tutorial is included in only free-suface-reduced and not included in level-set-method and surface-tension-coefficient.


Could you please teach me below,

Where is the tutorial about level set method and surface tension coefficient ?
or
Where is another tutorials from Elmer GUI Tutorials and Elmer non-GUI Tutorials ?
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: About how-to set surface tension in Level Set Method

Post by raback »

Hi

Tutorials exist only for a small fraction of all the capabilities. Levelset has not really been in the focus lately. Probably the only test case with surface tension is the test case "LevelsetFallingDrop". Additionally I would refer to the Models Manual. There are also some other levelset cases among tests but no other includes surface tension.

-Peter
SagittariusChiron
Posts: 3
Joined: 11 Jan 2021, 05:36
Antispam: Yes

Re: About how-to set surface tension in Level Set Method

Post by SagittariusChiron »

Thank for your advice.

I got the test case "LevelsetFallingDrop",and executed blow comand in Windows10:
- ElmerGrid 1 2 square.grd
- elmerf90 DropStuff.f90 -o DropStuff.dll
- ElmerSolver


ELMER SOLVER (v 9.0) STARTED AT: 2021/01/12 22:36:26
ParCommInit: Initialize #PEs: 1
MAIN:
MAIN: =============================================================
MAIN: ElmerSolver finite element software, Welcome!
MAIN: This program is free software licensed under (L)GPL
MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
MAIN: Version: 9.0 (Rev: Release, Compiled: 2020-11-11)
MAIN: Running one task without MPI parallelization.
MAIN: Running with just one thread per task.
MAIN: Lua interpreted linked in.
MAIN: =============================================================
ERROR:: ElmerSolver: Unable to find ELMERSOLVER_STARTINFO, can not execute.
STOP 1


There is a file "ELMERSOLVER_STARTINFO.txt" in a folder "ELMERSOLVER_STARTINFO".
Could you please teach me how to run the test case.
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: About how-to set surface tension in Level Set Method

Post by kevinarden »

ELMERSOLVER_STARTINFO is a text file that contains the name of the sif file to be run

It isn't required to be used you can do

ElmerSolver mycase.sif

to specify the sif file

although the test case looks to be correctly set up
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: About how-to set surface tension in Level Set Method

Post by raback »

The "ELMERSOLVER_STARTINFO" is a bare file. You probably introduced the .txt suffix yourself when copying the file into Windows.
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: About how-to set surface tension in Level Set Method

Post by kevinarden »

"ELMERSOLVER_STARTINFO.txt" windows wants to add a suffix to every file, but depending on options does not always show the suffix. This file name has no suffix, but it is a text file. If your suffix display option is off you will not see that it has a suffix, which it should no have one.
SagittariusChiron
Posts: 3
Joined: 11 Jan 2021, 05:36
Antispam: Yes

Re: About how-to set surface tension in Level Set Method

Post by SagittariusChiron »

Thanks !
I succeeded to run the test case "LevelsetFallingDrop".

Windows add a suffix "txt" to a file ELMERSOLVER_STARTINFO, but I don't want it.

By the way
this "LevelsetFallingDrop" case means the VOF method, unlike that "Axisymmetric coating process" case means the ALE(Arbitrary Lagrangian-Eulerian) method. And, I did not find "surface tension coefficient" in the "drop.sif" file.

I want to add the surface tension effect to free surface which transformed with ALE method.

Is there some hint for it ?
Attachments
temp.JPG
temp.JPG (43.83 KiB) Viewed 2365 times
raback
Site Admin
Posts: 4812
Joined: 22 Aug 2009, 11:57
Antispam: Yes
Location: Espoo, Finland
Contact:

Re: About how-to set surface tension in Level Set Method

Post by raback »

Hi

There is no "surface tension coefficient" since you don't have a surface in levelset. Instead the surface is smeared down and the surface tension becomes a force of type c\grad{d}. This is communicated to the Navier-Stokes equation by the following way:

Code: Select all

  Potential Force = Logical True
  Potential Field = Equals Surface
  Potential Coefficient = Equals Curvature
Here c is "curvature" and d is "surface". These are computed by the LevelSet solvers. I think equations 32.9-32.11 in Elmer Models Manual (page 190) are needed to close the problem.

No keyword is named "surface tension coefficient" because the levelset solver does not really know anything about the physics. It just moves a signed distance function and can compute some quantities on it. In the Navier-Stokes solver the force of this type can also be some electrostatic force, for example. Hence it would be limiting to narrow the keywords to the singular use type related to surface tension. So this is a side product of favoring atomic solvers that carry out their work independently.

-Peter
Post Reply