The direct and indirect PFBEM(power flow boundary element method) for the treatment of the 3 dimensional multi-domain problems are proposed to predict the acoustic energy density in medium to high frequency ranges. In the proposed method, the equation is derived in a matrix form by considering coupled relationships of the power flow at the interface of given domains. The proposed method can successfully obtain the analytical solutions for the problems of coupled cubes and the small-scale reverberant chamber. Then the experiment is carried out to obtain STL(sound transmission loss) by using small-scale reverberant chamber and the results are compared with analysis results.