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

sclsVOFFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 2, 2017, 14:37
Default sclsVOFFoam
  #1
New Member
 
Masoumeh Gharaati
Join Date: Apr 2017
Posts: 19
Rep Power: 5
mgharaati is on a distinguished road
Dear all

I 'm using sclsVOFFoam solver. http://www.tfd.chalmers.se/~hani/kur...ankarMenon.pdf
I 'm modeling droplet impact into liquid layer. the liquid is water and the gas is air. When I 'm using the
Quote:
transportProperties
it works but when I change viscosity to air and gas I face to these errors.

Quote:
p { margin-bottom: 0.1in; line-height: 120%; } #0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)
Would you please help me for solving this?
mgharaati is offline   Reply With Quote

Old   June 20, 2017, 23:59
Default
  #2
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 113
Rep Power: 8
ndtrong is on a distinguished road
"Foam::symGaussSeidelSmoother::smooth(Foam::wo rd const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:? "

I think there is problem of variable name, you should keep same variable name. On the other hand, you simply change value of the variable you want to change
ndtrong is offline   Reply With Quote

Old   June 22, 2017, 12:51
Default
  #3
New Member
 
Masoumeh Gharaati
Join Date: Apr 2017
Posts: 19
Rep Power: 5
mgharaati is on a distinguished road
I just changed the value of variables. I didn't change any name.
mgharaati is offline   Reply With Quote

Old   June 22, 2017, 13:27
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,924
Rep Power: 35
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi all,

@mgharaati

Your solution diverges. Why your solution diverges? It is difficult to say from what you have posted.

This

Quote:
when I change viscosity to air and gas
is rather enigmatic.
alexeym is offline   Reply With Quote

Old   June 26, 2017, 16:49
Default
  #5
New Member
 
Masoumeh Gharaati
Join Date: Apr 2017
Posts: 19
Rep Power: 5
mgharaati is on a distinguished road
It 's related to Courant number. I check the log file and Courant number got too much. Then, I tried to keep Courant number under 0.1 but it didn't work too.
mgharaati is offline   Reply With Quote

Old   June 28, 2017, 04:05
Default
  #6
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 113
Rep Power: 8
ndtrong is on a distinguished road
I don't think it work by changing of variable's value. I mean you change the name of variable, you should look into the createFields file to see how your code read variable name.
ndtrong is offline   Reply With Quote

Old   June 28, 2017, 04:07
Default
  #7
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 113
Rep Power: 8
ndtrong is on a distinguished road
please upload your transportProperties file in constant folder
ndtrong is offline   Reply With Quote

Old   July 5, 2017, 15:17
Default
  #8
New Member
 
Masoumeh Gharaati
Join Date: Apr 2017
Posts: 19
Rep Power: 5
mgharaati is on a distinguished road
Please find the attached file.
p.zip
mgharaati is offline   Reply With Quote

Old   July 6, 2017, 06:25
Default
  #9
Senior Member
 
Nguyen Duy Trong
Join Date: Apr 2014
Posts: 113
Rep Power: 8
ndtrong is on a distinguished road
Quote:
Originally Posted by mgharaati View Post
Please find the attached file.
Attachment 57111
I am so sorry. I am wrong with my opinion.
Maybe there are another reason for your errors.

Sorry about that
ndtrong is offline   Reply With Quote

Old   September 30, 2018, 02:17
Default
  #10
Senior Member
 
krishna kant
Join Date: Feb 2016
Location: Hyderabad, India
Posts: 108
Rep Power: 6
kk415 is on a distinguished road
Hi All


I am also solving sclsVOFFoamsf and getting a similar error like shown in the original post. It is coming from the reinitialization step. Can anyone guide what could possibly went wrong?


solve the reinitialization equation


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? at sigaction.c:?
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#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 ??:?
#5 ? at ??:?
#6 __libc_start_main at ??:?
#7 ? at ??:?
Floating point exception (core dumped)


I have individually look on the output of all individual values involved like psi, grad(psi) deltaTau they all look fine. The boundary condition is cyclic for all variables.
kk415 is offline   Reply With Quote

Old   September 30, 2018, 03:32
Default
  #11
New Member
 
Masoumeh Gharaati
Join Date: Apr 2017
Posts: 19
Rep Power: 5
mgharaati is on a distinguished road
Hi

I am afraid to say that I couldn't solve the problem and I still get the same error.


Quote:
Originally Posted by kk415 View Post
Hi All


I am also solving sclsVOFFoamsf and getting a similar error like shown in the original post. It is coming from the reinitialization step. Can anyone guide what could possibly went wrong?


solve the reinitialization equation


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? at sigaction.c:?
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#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 ??:?
#5 ? at ??:?
#6 __libc_start_main at ??:?
#7 ? at ??:?
Floating point exception (core dumped)


I have individually look on the output of all individual values involved like psi, grad(psi) deltaTau they all look fine. The boundary condition is cyclic for all variables.
mgharaati is offline   Reply With Quote

Old   October 31, 2018, 22:53
Default
  #12
Senior Member
 
Huang Xianbei
Join Date: Sep 2013
Location: Yangzhou,China
Posts: 298
Rep Power: 10
huangxianbei is on a distinguished road
Quote:
Originally Posted by kk415 View Post
Hi All


I am also solving sclsVOFFoamsf and getting a similar error like shown in the original post. It is coming from the reinitialization step. Can anyone guide what could possibly went wrong?


solve the reinitialization equation


#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? at sigaction.c:?
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#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 ??:?
#5 ? at ??:?
#6 __libc_start_main at ??:?
#7 ? at ??:?
Floating point exception (core dumped)


I have individually look on the output of all individual values involved like psi, grad(psi) deltaTau they all look fine. The boundary condition is cyclic for all variables.
Hi, I got the same error. I tried several turbulence models and still got the error. I think it should be attributed to the solver rather than the turbulence models.

The problem remains in the solveLSFunction.H:
psi = psi + psi0/mag(psi0)...

psi can be zero so that the divide operator will lead to divergence.

Change:

psi = psi + psi0/max(mag(psi0),VSMALL)

This can avoid the error. However, I still can not get a converged solution for my case. The calculation crashed after a few iterations

Last edited by huangxianbei; November 1, 2018 at 22:20.
huangxianbei is offline   Reply With Quote

Reply

Tags
clsvof, viscosity

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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



All times are GMT -4. The time now is 12:33.