CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

Having difficulties with coupling interFoam with Level Set method

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 21, 2017, 07:41
Default Having difficulties with coupling interFoam with Level Set method
  #1
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 124
Rep Power: 12
ndtrong is on a distinguished road
Hi everyone,

I am coupling interFoam vs Level Set (LS) method. The coupling manner is followed by tutorial of Sankar Menon at:
http://www.tfd.chalmers.se/~hani/kurser/OS_CFD/

for the first step, if I only care about the surface tension (sclsVOFsf solver), it can done quite well. However, when I take into account correction of density and viscosity (sclsVOF solver) it requires declaration of Heaviside function H in fvOption. I guess H and alpha.water are similar then I declared they are same solver. However, it crashed as in the attached file.

It would be appreciate if someone could help me solve this problem.

Thank you in advance.
Attached Files
File Type: txt log.myCLSVOF.txt (22.3 KB, 33 views)
File Type: txt fvSolution.txt (3.3 KB, 28 views)
ndtrong is offline   Reply With Quote

Old   October 4, 2017, 04:58
Default error in Couple Level Set with VOF interFoam
  #2
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 124
Rep Power: 12
ndtrong is on a distinguished road
Dear Foamers,

I am trying to recompile the sclsVOFFoamsf of Dr. T. Yamamoto for OpenFoam v5. Compilation seems good and solver can run for few step, however, I got the errors as follows:
[4] #0 [5] [3] #0 Foam::error:rintStack(Foam::Ostream&)#0 Foam::error:rintStack(Foam::Ostream&)Foam::error :rintStack(Foam::Ostream&) at ??:?
at ??:?
[3] #1 [5] #1 Foam::sigFpe::sigHandler(int)Foam::sigFpe::sigHand ler(int) at ??:?
[4] #1 Foam::sigFpe::sigHandler(int) at ??:?
[4] #2 ? at ??:?
[5] #2 ? at ??:?
[3] #2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[4] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[3] #4 at ??:?
[4] #5 at ??:?
[5] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&)?Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam:perator/<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
[3] #5 at ??:?
[5] #5 at ??:?
[4] #6 __libc_start_main?? in "/lib/x86_64-linux-gnu/libc.so.6"
[4] #7 at ??:?
[5] #6 __libc_start_main at ??:?
[3] #6 __libc_start_main? in "/lib/x86_64-linux-gnu/libc.so.6"
[5] #7 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #7 ? at ??:?

It would be appreciate if someone could help me to fix this error. The log file of solver is as in the attached file.

Thank you
Attached Files
File Type: txt log.sclsVOFFoam.txt (10.4 KB, 13 views)
ndtrong is offline   Reply With Quote

Old   October 7, 2017, 08:32
Default
  #3
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10
saddy is on a distinguished road
plz give more information regarding where did you get the code of sclsvofsf and what changes have you made in order to upgrade it to compile it on openfoam 5.
i encountered similar problems while recompiling evapvofhardt solver.
it happens so that compiling successfully and running successfully are two different things
saddy is offline   Reply With Quote

Old   October 9, 2017, 13:09
Default
  #4
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 124
Rep Power: 12
ndtrong is on a distinguished road
Quote:
Originally Posted by saddy View Post
plz give more information regarding where did you get the code of sclsvofsf and what changes have you made in order to upgrade it to compile it on openfoam 5.
i encountered similar problems while recompiling evapvofhardt solver.
it happens so that compiling successfully and running successfully are two different things
Dear Sandy

I changed values of cAlpha and maxDeltaT and it can run well.
Thank you so much for your comment
ndtrong is offline   Reply With Quote

Old   October 13, 2017, 04:44
Default
  #5
New Member
 
Sumeet Sambhus
Join Date: Aug 2017
Posts: 2
Rep Power: 0
sumeet_s is on a distinguished road
Dear Nguyen,

I am new to Openfoam and want to solve problems using CLSVOF solver.
Could you please provide your solver ?

Thanks in advance!

regards,
Sumeet.
sumeet_s is offline   Reply With Quote

Old   October 18, 2017, 04:39
Default
  #6
Senior Member
 
sandy
Join Date: Feb 2016
Location: .
Posts: 117
Rep Power: 10
saddy is on a distinguished road
here's the link to solver and report + presentation
www.tfd.chalmers.se/~hani/kurser/OS_CFD/

ctr+f: sankar menon
saddy is offline   Reply With Quote

Old   October 18, 2017, 06:36
Default
  #7
New Member
 
Sumeet Sambhus
Join Date: Aug 2017
Posts: 2
Rep Power: 0
sumeet_s is on a distinguished road
Thanks sandy.
sumeet_s is offline   Reply With Quote

Old   January 21, 2018, 18:45
Default
  #8
cyw
Member
 
Yuanwei Cao
Join Date: Jun 2017
Location: Munich,Germany
Posts: 54
Rep Power: 9
cyw is on a distinguished road
Hi, have you tested the sclsvofFoam solver? Such as dambreak, rum in parallel?
cyw is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
setting of level set method as initial condition sooroo3 FLUENT 4 July 31, 2018 05:20
Level set method---reinitiallization and boundary conditions wenxu OpenFOAM Running, Solving & CFD 2 February 28, 2018 12:43
Level Set Method - OpenFOAM-1.5-dev Ingenierias2003 OpenFOAM 2 August 7, 2012 05:36
calculating Normal vector in level set method amir2920 Main CFD Forum 1 July 21, 2009 08:25
Level set method for detonation Amir Main CFD Forum 1 July 2, 2008 16:25


All times are GMT -4. The time now is 21:35.