  Coronal loops act as resonant cavities for low frequency fluctuations that
are transmitted from the deeper layers of the solar atmosphere and are
amplified in the corona, triggering nonlinear interactions. However trapping is
not perfect, some energy leaks down to the chromosphere, thus limiting the
turbulence development and the associated heating. We consider the combined
effects of turbulence and leakage in determining the energy level and
associated heating rate in models of coronal loops which include the
chromosphere and transition region. We use a piece-wise constant model for the
Alfven speed and a Reduced MHD - Shell model to describe the interplay between
turbulent dynamics in the direction perpendicular to the mean field and
propagation along the field. Turbulence is sustained by incoming fluctuations
which are equivalent, in the line-tied case, to forcing by the photospheric
shear flows. While varying the turbulence strength, we compare systematically
the average coronal energy level (E) and dissipation rate (D) in three models
with increasing complexity: the classical closed model, the semi-open corona
model, and the corona-chromosphere (or 3-layer) model, the latter two models
allowing energy leakage. We find that:
  (i) Leakage always plays a role (even for strong turbulence), E and D are
systematically lower than in the line-tied model. (ii) E is close to the
resonant prediction, i.e., assuming effective turbulent correlation time longer
than the Alfven coronal crossing time (Ta). (iii) D is close to the value given
by the ratio of photospheric energy divided by Ta (iv) The coronal spectra
exibits an inertial range with 5/3 spectral slope, and a large scale peak of
trapped resonant modes that inhibit nonlinear couplings. (v) In the realistic
3-layer model, the two-component spectrum leads to a damping time equal to the
Kolmogorov time reduced by a factor u_rms/Va_corona
