Dear Smart $,
I tried to use different method, that is using complex number operation under Ia = 1.41 * Ia rms * Sin (ωt + ψa). Adding phase current under complex number gets 53.65A (In = -53.642 - j 0.94). This is different to 58.86A calculated.
Where am I wrong?
How am I supposed to determine where you went wrong if you only provide the format of your vector math and your wrong answer. I would need to see the entire calculation. Otherwise all I can tell you is you went wrong somewhere in the middle
If you don't want to make your calculation public, PM it to me...
Another question, can we calculate neutral current of 3 phase 4-wire system when only rms currents of phases are known? (power factor is not known).
Calculating with peak values (moduli) vs. rms values does not matter, as long as you are consistent throughout the calculations with the same type of values. As for calculating currents without knowing the power factor, you could assume all power factors are the same if all loads are of the same design. If not, it would be best to use approximate, typical power factors. As exemplified by this thread, accuracy of the resulting value can vary widely without knowing the power factor (or its equivalent).
Is the neutral current cause waste of energy under dissipated heat?
Yes there will be the typical I?R losses associated with any neutral current. But if it is cold outside and the neutral conductor is inside a facility that relies on some type of space heating to keep it warmer than outside, then the energy is not wasted until it is warmer outside than it is inside.