Hi Jur,
The difference between your -aa 0.08 simulations could arise from having -ar too large. One option would be to set a smaller -ar with a larger -aa to prevent the run time from increasing too much.
Do you have measurements to compare your simulation results to? Irradiance caching is a biased algorithm, so I would not immediately assume that the Radiance result is more accurate. Accelerad's built-in overture pass generally increases its reliability compared to Radiance.
For your -aa 0 results, it looks like you are using the default -lw value of 2e-3. This is probably too large, which would cause you to calculate a lower daylight factor. I recommend setting -lw to less than 1/ad, so no more than 4e-4 in your case.
Regarding triangulation, Accelerad does require that each polygon has a consistent winding direction. In Radiance, some programs also require consistent winding direction, but most do not. This mainly affects polygons that contain internal holes. Most geometry export tools wind polygons correctly, but there are a few such as su2rad which do not. If you are creating your .rad files with an export tool that does not wind polygons correctly, then the solution is to mesh them before exporting.
Radiance and Accelerad also both have limits on the size and closeness of objects due to their use of floating point precision, though their limits are different. This should not affect you unless your model is located very far from the coordinate system origin.
Nathaniel