[mgharaati]

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

transportProperties

it works but when I change viscosity to air and gas I face to these errors.

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?

[ndtrong]

"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

[mgharaati]

I just changed the value of variables. I didn't change any name.

[alexeym]

Hi all,

@mgharaati

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

This

when I change viscosity to air and gas

is rather enigmatic.

[mgharaati]

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.

[ndtrong]

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]

please upload your transportProperties file in constant folder

[mgharaati]

Please find the attached file.
p.zip

[ndtrong]

Originally Posted by mgharaati

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

[kk415]

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]

Hi

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

Originally Posted by kk415

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.

[huangxianbei]

Originally Posted by kk415

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