Ting,
I am afraid ther is no such formula. That is why simulation is the normal method to estimate the performance. Such simulation usually takes hours to complete.
For very challenging application like this one, the high frequency resistance Rs must be made very small, and coil diameter large enough. Multi-strand wire can reduce Rs and can be used to wind the coil.
Best regards,
Don