Accurate and fast calculation of frequency-dependent leakage inductance is critical for the optimization of the high-frequency transformer with Litz wire winding. The analytical method is widely used due to its high computation speed. However, its computation accuracy is limited due to the following factors. i) One-dimensional (1D) or two-dimensional magnetic field model (2D) has unsatisfying accuracy for the transformer with a rotationally asymmetric structure. ii) The skin and proximity effect factors based on the multilayer foil winding are inaccurate for the Litz wire winding in the calculation of frequency-dependent magnetic energy. This paper proposes a double-2D frequency-dependent leakage inductance model for the transformer with Litz wire winding. It uses the image method to calculate the 2D leakage magnetic field inside and outside the core window, respectively, which considers the rotationally asymmetric structure in leakage inductance calculation. The round conductor model is proposed for the frequency-dependent magnetic energy calculation. The single strand in a Litz wire is taken as the basic modeling object of the round conductor model, which is more detailed than the multilayer foil winding model for the Litz wire winding. The comparison between the experiment, FEM, state-of-the-art models and the proposed model is conducted for eight transformers with E-core, U-core, round or rectangular Litz wire winding. The proposed model shows satisfying computation accuracy in the whole frequency band. And the computational speed of the proposed model is 100 times that of the FEM.