sclsVOFFoam
 

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 it works but when I change viscosity to air and gas I face to these errors.

Would you please help me for solving this?

"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

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

Hi all,

@mgharaati

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

This

is rather enigmatic.

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.

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.

please upload your transportProperties file in constant folder

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

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::printStack(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::operator/<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 am afraid to say that I couldn't solve the problem and I still get the same error.

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