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!