CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   StressedFoam (

ravi November 22, 2005 00:58

Hi, I am interested to use

I am interested to use 'stressedfoam' solver to calculate thermal stress. For learning, I tried to solve tutorial example 'cavity' with thermal stress 'on'.

The displacement field 'U' seems to fine but stressed field shows some numerical problem. At the block interface, it shows some kind of numerical error in stress field. I think, I am making some silly mistake ?

I would be glad, if some help is provided to solve this problem.


mattijs November 22, 2005 05:03

Does the plateHole tutorial wo
Does the plateHole tutorial work with thermal stresses? What block interface is this? cavity is only one block?

ravi November 22, 2005 06:08

Thanks for prompt response.
Thanks for prompt response.

I was trying to add thermal stress in plateHole problem by changing boundary condition for hole surface temperature (T= 400). I had changed parameter 'thermalStress' as 'yes' in file 'constant/thermalProperties'. I suppose by doing so I can solve thermal stresses.

About block interface:
In file 'blockMeshDict' we use command 'blocks' to create multiple block (in this example there are 5 blocks). We define boundary patches but there are internal patches which connect this blocks. I am getting some spurious solution near the internal connecting patches.


hjasak November 22, 2005 09:00

You need a better gradient cal
You need a better gradient calculation algorithm to deal with the local mesh distortion.

Change the default gradient scheme to least squares and all will be well. If you can still see problems, please shout and I'll have a look at it.

Something like this:


default leastSquares;



ravi November 28, 2005 04:02

Thank you very much. It works
Thank you very much. It works properly with "least squares" gradient schemes.


All times are GMT -4. The time now is 00:03.