Hi Francis,
Thank you for your advice. It also lets me solve my problem. BR. |
Hi,
I am using a Cht K epsion model using Openfoam and I keep getting this error. I have tried changing the maxiter to 10000 and still get the problem too. Did anyone get to the bottom of this? I completely remade my mesh and still get this same error although I have a few overly skewed faces (64 failed skewness faces) and 900 faces failing nonOrthag check, however due to the geometry even if I used SHM these errors would still occur. I tried Francis's method and this has't solved it. I have noticed a sharp (impossible) rise in density in certain regions of the fluid, this was so sharp that it shouted alarm. However I am unfamiliar with this type of modelling so thought I would ask for some advice. I previously attended the Advanced Openfoam training course and have emailed the guy who ran it (as has my work mate) regarding problems which he encouraged. Unfortunately he does not respond to emails, which has left me a bit discouraged. Cheers |
Hey everyone,
I recently ran into the same issue using the sprayFoam solver with version OF2.4 and an implemented ignition. The error especially ocurrs when I ignite more than one cell and I think it must be due to the Newton solver used to calculate the temperature. Eventually it ocurrs with only one ingiting cell as well though. Has anyone solved this issue so far? @sahul could you provide the implented brent solver? Thanks in advance Moritz |
Hi Moritz,
Have you had any luck on this one? I have noticed that my temperatures are diverging too for no reason. I riginally solved the issue by switching to the plane wall 2d solver, however I have got the error again now. |
Maximum number of iterations exceeded
I've been having a lot of trouble tracing the source of this error. I know it happens in the thermoI.H file, but I can't trace what is calling that, or figure out any reason why this is happening. I've had no luck trying to fix this on my own for a while now.
The case i am running is a simple wedge (meshed with snappy hex mesh) inside a blockmesh region using rhoCentral foam. any help would be really appreciated this is the error log file Code:
/*---------------------------------------------------------------------------*\ p: Code:
dimensions [1 -1 -2 0 0 0 0]; Code:
dimensions [0 0 0 1 0 0 0]; Code:
dimensions [0 1 -1 0 0 0 0]; |
Hi Henry,
Most probably your T value is going to negative or extremely high. To see in into T.Eqn file add Info << T.min() << " " << T.min() <<endl; then you can understand reason. hope it helps you. Baris |
It looks like that was the problem, but I'm lost as to how to go about solving it. What can I do to prevent T from going out of control like that?
|
Hey together,
could anyone please explain me how I can compile the thermo files if I want for example some output debugging in the thermoI.H file??? I somehow don't understand where these files are included in any Make/files... Any help would be sincerely appreciated Thanks and cheers moritz |
go to the directory that includes the thermoI.H file, and once you edit whatever debugging code you want into that file, type "wclean" and then "wcompile", which should recompile all the files in that directory for you
|
Hi henry,
I didnt want to do it as root but I guess it is easier. But still it's not giving me the info in my logfile I acutally wrote in the code. Why? Code:
template<class Thermo, template<class> class Type> |
Ok, sorry for the stupid question. You have to compile the whole thermophysical folder if anyone faces the same problem at some point.
Thanks anyway guys! |
chtMultiRegionSimpleFoam - Solution suggestion
Hello Everybody,
Increasing the Mesh quality by increasing the number of cells solved my divergence problem and no error massage appeared anymore. - My mesh is created by blockMesh - chtMultiRegionSimpleFoam has been used as a test - I left the maximum iterations by 100 (did not followed the suggestion above) - in my test I gave: runTimeModifiable false I hope it can helps Regards Peter |
Dear foamer,
I solved the problem. In my case, the error occurred because of different values for T_common in the therm.dat file. Make sure that all reactants involved in the chemical reactions have the same value for T_common. But I think in the current version this bug is already fixed. Cheers |
Quote:
|
Resolved
1) Reduce delta T to a order 10000 times lower than the current value.
2) Set Max Courant Number to 1 3) Set ´adjustTimeStep´ to yes try running the solver again, if it crashes yet again, reduce the delta T even lower then try running again. This worked for me. The problem was with higher courant number and that leading to a diverged solution. |
Quote:
We are also struggling with this issue of temperature update using N-R method. Can you share your brent method code implemented in thermoI.H file? We also feel that by using Brent method, this problem can be solved. Hrushilesh |
hi
Did you find a solution to your problem? I'm having the same problem with the compressible solver rhoSimpleFoam. Thanks in advance! |
Hi
I'm having the same ṕroblem you mentioned. Did you manage to solve it? what did you do? Thanks in advance Andrea |
Quote:
Hi srahul I have some questions, How can make the thermoI.H file? Can you send brent solver text in openFoam? thanks |
My mistake was...
Hello Everybody!
chtMultiRegionSimpleFoam has been used here! I faced the problem also like many of you. I tried almost everything has been written here without a success... I increased the number of cells to a very high number (1 000 000) cells for 100 x 100 x 100 mm region(s). That worked, but the results where very strange. And I found the mistake I made! I made the folowing: - 3D Design FreeCAD then export as stl - Used blockMesh to make my initial mesh - In openFoam using snappyMesh to build my final mesh The dimensions in my stl file were in mm. That means FreeCAD exports the 100... snappyMesh read this 100 as 100 m instead of the 100 mm exported My blockMesh mesh is also too big made (factor 1000)... Reducing the dimensions in factor 1000 bei: convertToMeters 0.001 in blockMeshDict AND in my exported stl file (FreeCAD) Solved my problem! Check up the dimensions (in your blockMeshDict as example) is my suggestion here. Regards Peter |
All times are GMT -4. The time now is 08:38. |