Temperature blowup for with 6-DOF falling object

Andrew Wick

I am modeling an object falling through air using the dynamic mesh function and the 6-DOF solver. It works perfectly for an incompressible flow, but when I attempt to turn on the energy equation, the temperature diverges between 0.3 and 1 second of flow time. I am using Piso with Non-Iterative time advancement.

There is one grid point on the body surface with a temperature of 1 and another with a very large temperature. This causes the solution to diverge. I've tried changing the timestep, the under-relaxation factors, the grid, the convective schemes, and this always pops up. The object's BC's are adiabatic walls. When I switch to isothermal, this still pops up, only its later in time and a few grid points away from the body.

Anyone run into this problem before or can offer some insight?