I'm encountering some strange convergence issues in the temperature solver in a coupled NS/Heat equation + mesh update spinup run.

Although my simulation is transient, I'm attempting to solve the steady state temperature equation by specifying:

Code: Select all

`Steady State Condition = Real 1.0 `

Interestingly, this seems to work (though it takes a while to converge) if I switch off the Mesh Update solver, but not if it is on. The strange part is, if I specify in the body force:

Code: Select all

` Mesh Update = Real 0.0 `

This is really confusing me. It seems as if MeshSolve.F90 is doing something, aside from displacing the mesh, which is messing up the temperature solution on the following timestep.

A final piece of the puzzle: If I turn off:

Code: Select all

```
Temp Load = Variable Velocity 1
Real Procedure "ElmerIceUSF" "getFrictionLoads"
```

If anyone has any ideas what might be causing this, I'd be very grateful!