[ravi]

Hi,

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.

Ravi

[mattijs]

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

[ravi]

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.

Ravi

[hjasak]

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:

system/fvSchemes

gradSchemes
{
default leastSquares;
}

Enjoy,

Hrv

[ravi]

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

Ravi