Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOp
 

Hello Dear Foamers,
I am a newbie to OpenFoam and would like to seek your help to solve an internal flow problem inside a valve.
Its a Steady state turbulent flow and I get the error after 6-8 iterations as quoted below while trying to solve with SimpleFoam.

Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#10
at simpleFoam.C:0
#11
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"

I hereby attach my fv solution, fv schemes and log files. Any help/guidance is much appreciated.

Thanks
Nagesh

Hi,

could you post:

1. checkMesh output (using CODE tag or as an attachment)
2. your boundary conditions (maybe at archive of your 0 folder)

Hi Alexey Matveichev...
Thanks for the reply.Here are my boundary conditions file. Will upload my checkMesh file in the following post.
I now realise that nut file wasn't needed for a K-epsilon model. But just want to clarify if it does any harm afterall?

Regards
Nagesh

checkMesh file
 

And here is the checkMesh file.

Thanks and regards
Nagesh

Correction with error message
 

My apologies.Due to a wrong epsilon value, I had got the error posted above.
Although I corrected it according to my problem, the error still persists and is as below.

Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#8
at simpleFoam.C:0
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"

Well, already at this point

it's more-or-less obvious, you've got problem with turbulence model. It can be due to mesh, ICs, or BCs.

How did you calculate IC and BC values for k and epsilon? You're using fixedValue BC for both, though maybe it's better to use turbulentIntensityKineticEnergyInlet and turbulentMixingLengthDissipationRateInlet with intensity and mixing length estimated from Re and hydraulic radius.

Hello Alexy,
Thanks for your insight to my problem. I calculated the values for k and epsilon using the expressions which is herewith attached file.
I shall try out your suggestion and get back if the problem still persists.

Sorry if the attached file is of any inconvenience.

Thanks and regards
Nagesh

Hello again,
As per your suggestions I made the changes accordingly.And I get quite a different error but at the same point of time. (27th iteration).
Vaguely I understand that it might be a problem with Mesh.
However I would like to know if mesh is the only problem or might there be problem with BC s and IC s again.

Here is the error. Kindly provide with your inputs.


#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::multiply<Foam::Tensor<double> >(Foam::Field<Foam::Tensor<double> >&, Foam::UList<double> const&, Foam::UList<Foam::Tensor<double> > const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#4 void Foam::multiply<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<Foam::Tensor<d ouble>, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#5 Foam::tmp<Foam::GeometricField<Foam::Tensor<double >, Foam::fvPatchField, Foam::volMesh> > Foam::operator*<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&, Foam::tmp<Foam::GeometricField<Foam::Tensor<double >, Foam::fvPatchField, Foam::volMesh> > const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libincompressibleTurbulenceModel.so"
#6 Foam::incompressible::RASModels::kEpsilon::divDevR eff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"

Hi,

for the future: it'll be more convenient, if you provide part of the log before the error.

As for the post, suggestions will be rather generic:

1. You've got non-hex cells in your mesh, use "leastSquares" instead of "Gauss linear" for gradSchemes

2. Set divSchemes to first-order upwind

3. Reduce relaxation factors for k and epsilon (let's say 0.3). Though if the problem is elsewhere this will just move FPE to later time.

4. Finally, if the problem persists, try moving from "corrected" schemes to "limited corrected 0.5".

Hi,
Firstly thank you for your valuable inputs. Although they helped me to keep my solver running, I get worst results.
I hereby would like to know if its because of mesh/ICs BCs or could there still be problem with my fv schemes.
I have attached the log file herewith,the solver is still running,however results are very bad.
I have an inlet velocity of 5 m/s and at the end of 5th iteration,it shoots upto 36 m/s and the same with pressure as well.(above 500 at the end of 5th time step).
I am having a tough time with my little knowledge to dodge this problem.
Any inputs is highly appreciated.

Well, dynamics of residuals looks quite promising, i.e. they are reducing during iterations.

Though if you plan to compare results of the simulation with experimental values (or just would like to get something meaningful), you should set convergence criterion for SIMPLE algorithm. Currently you've got none:

add residualControl dictionary to your SIMPLE dictionary in fvSolution, so it looks like:

maybe you don't need 1e-6, maybe 1e-4 will be enough.

Aha yes.I thought the same about residuals as well. I have set up the convergence criterian now.
But I feel I might have gone wrong with the mixing length.Could you please tell me how to calculate mixing length?

1. http://www.cfd-online.com/Wiki/Turbu...ary_conditions
2. http://www.cfd-online.com/Wiki/Turbulent_length_scale

Usually I use 7% of hydraulic radius.

Thank you so much for all your suggestions.But for the fact that results seem to be too unrealistic,the solver works fine. Will work on it to get a satisfying result.

Same Problem persists
 

Hello,
I had previously solved for one half of a Valve and now I am trying to analyse flow inside a Complete Valve.Unfortunately, I get the error which I had got earlier during Half-valve model.
Obviously I have changed BCs and ICs as per flow and model requirements and the other settings remaining the same as before.
It would be of great help for me if anyone could throw some light on where the problem lies.
Herewith I am attaching all the related files,which might be helpful.

Thanks
Nagesh

Attachment 35338

Attachment 35339

Attachment 35340

Attachment 35341

Attachment 35342

Attachment 35343Attachment 35344

Oops, and this is the error msg that I get

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#8
at simpleFoam.C:0
#9
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
in "/opt/openfoam222/platforms/linux64GccDPOpt/bin/simpleFoam"

Hi Tobi,
First of all thanks for the reply. And I am afraid I cant post the pictures due to confidentiality issues.However,I can look for it myself if you let me know what aspect to monitor.Was it because of snap layers( Morphing) warning?

You suggested me not to use smooth solver.May I know what I can use instead?

And sorry for the inconvenience with the error msges. I shall correct myself in the future posts.

Thanks and regards
Nagesh

Hello, Dear Nagesh Atreyas:
It seems that I came cross a problem similar to this one.
Could you share your way to solve this problem for me?
I am looking forward to your help. Thank you very much.

We need more input to solve your question. Can you please share your logfile (stdout and stderr please).

Dear Nagesh,

Please excuse for the late reply. Try the following changes in the attachments:

and

Re-run your case, I think this will work

-
Best Luck!

Hi Tushar,
Thanks for the suggestions. I have solved this problem already.I changed the grad schemes to gaussian and others to first order,increased all the relaxations apart from p to 0.7, also there were few problems with bc.

However, I am keen on understanding about your suggestion. Changing to linear schemes is clear, it is because of accuracy,but why were you hinting on changing the relaxation factor alone to 0.7??! Is it very evident from the error to do so or how?Thanks.

Regards
Nagesh

Dear Nagesh,

Many things you will come to know with experience. In short, I would say that it is more related to theory.

_
Best Luck!

Foam::error::printStack(Foam::Ostream&) at ??:?
 

Hi all,
I am getting the same error as discussed in above posts .I am pasting what error i am getting.

Reading field U

Reading/calculating face flux field phi

No finite volume options present


SIMPLE: no convergence criteria found. Calculations will run for 0.5 steps.


Starting time loop

Time = 0.005

smoothSolver: Solving for Ux, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::divide(Foam::Field<double>&, double const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::operator/<Foam::fvPatchField, Foam::volMesh>(Foam::dimensioned<double> const&, Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5 ? at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? at ??:?
Floating point exception (core dumped)

Can anybody help me with this . where is the error. Any help will be appreciated....

Regards
AJAY

Hey Ajay,
In my experience this error is often the result of incorrectly initialized values. Check with ur bcs'. It happens when the solver comes across expressions like 0/0 or something that is equally strange.
Or your mesh Quality is not good enough. Check with the tool checkMesh.
Initial and final residuals=0 is something that I havent come across until now. May be someone else can comment better on these values?!
Cheers.

Nagesh,

Hi were u able to figure out the problem ??

Hi,

I am new to OpenFOAM. Currently I'm trying to simulate flow through an S duct and I'm getting the following error when I'm trying to run simpleFoam

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libpthread.so.0"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::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)

Can someone pls explain me where I have gone wrong? Pls help me if you can :(

See your #1 error. There it says sigFpe error.

It means that you have some variable value which is making 0/0 (divison by zero) form somewhere in your solver.

Also you have to check the boundary conditions of your variables in the solver. This error also comes when there are wrong BC.

Hope this helps.

Best
Ajay

I just want to mention one more thing. In general this error means a division by zero. As Ajay said, you can check out the error messages to get the piece of code where the problem occur. But normally it is related to the initial set-up. If you get the error right after the start (without any solving), then there should be a problem with the BC (as Ajay mentioned). If you have some iterations before, then you have either a mesh problem or boundary problem. But there are also other things that can throw out this exception.

Thank you Ajay and Tobi.
But I'm not understanding where to check for the 0/0 error or the BC error.
I've given in my 0 directory the values for k and epsilon I've calculated and walls have been assigned specific wall criteria. I'm not quite sure where to edit:(

Hi,
I'm working on a 65mmX65mm inlet and exit, s shaped duct of 90 degree rotation. u=15.79m/s at inlet, pressure 0 at exit. epsilon in 0 file=24.92, k=0.35, nut=4.424e-4, nuTilda=0. In transport properties, nu=1.50934e-5.
I'm using pimplefoam for solving the system directories are as below:

controldict

application pimpleFoam;
startFrom latestTime;
startTime 0;
stopAt endTime;
endTime 1;
deltaT 0.0001;
writeControl adjustableRunTime;
writeInterval 0.001;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
adjustTimeStep yes;
maxCo 5;

fvschemes

ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,k) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,R) bounded Gauss upwind;
div(R) Gauss linear;
div(phi,nuTilda) bounded Gauss upwind;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian(rAUf,p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}
interpolationSchemes
{
default linear;
interpolate(U) linear;
}
snGradSchemes
{
default corrected;
}

fvsolution

solvers
{
p
{
solver GAMG;
tolerance 1e-7;
relTol 0.01;
smoother DICGaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
pFinal
{
$p;
relTol 0;
}
"(U|k|epsilon)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-05;
relTol 0.1;
}
"(U|k|epsilon)Final"
{
$U;
relTol 0;
}
}
PIMPLE
{
nNonOrthogonalCorrectors 0;
nCorrectors 2;
}

Please help if you find any mistake in my calculations or the system properties. I've tried all what I've understood. Not sure why the following error is coming up :

PIMPLE: iteration 1
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0083465, No Iterations 1000
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0303338, No Iterations 1000
GAMG: Solving for p, Initial residual = 1, Final residual = 1.54301e+34, No Iterations 1000
time step continuity errors : sum local = 2.68132e+31, global = 1.9573e+26, cumulative = 1.9573e+26
GAMG: Solving for p, Initial residual = 0.970831, Final residual = 2.76941e+30, No Iterations 1000
time step continuity errors : sum local = 7.67121e+61, global = 2.29341e+57, cumulative = 2.29341e+57
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libpthread.so.0"
#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 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#10 Foam::RASModels::kEpsilon<Foam::IncompressibleTurb ulenceModel<Foam::transportModel> >::correct() at ??:?
#11 ? at ??:?
#12 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13 ? at ??:?
Floating point exception (core dumped)

Hi,

first of all, please use code tags. The formating of your post is really everything else than readable-friendly. In the first lines of the solver, you can see that you have 1000 iterations for U and p and the final residuals are increasing. This give you the hint that your set-up of the BCs are probably complete wrong and the matrix system you are building and trying to solve has a lot of solutions. Check them, make sure that it is physical. In addition a mesh problem, complete wrong initialization parameters and other stuff can cause this. But having the final residuals going to 1e30 should be a good indicator that there is an fatal mistake (just check the continuity error). Good luck.

Hi,

Since your sigFpe error is coming in initial iteration there is a high probability that you have some wrong BC in your case which is leading to this error.

Also as Tobi mentioned this error also comes when your meshing is wrong.

I have encountered this error many times. So, by changing BC, Analyzing your solver code to see if any variable value making 0/0 form there and checking the mesh length and number of divisons made in x and y directions in your blockMeshDict file helps a lot to debug this error.

Hope this helps

Best
Ajay

Hi,
I simulate the water-gas two phase using interFoam(openfoam3.0.0). But i got a error message. I change the boundary conditions, and try it again and again. It still does not work. There are always the same mistakes(Foam::error...). Could you please help me understand this error?

Why my deltaT is so small? I couldn't understand. It has been bothering me for a long time. I really hope someone here could help me. Thank you.