We demonstrate experimentally an efficient control of light intensity distribution inside a random scattering system. The adaptive wavefront shaping technique is applied to a silicon waveguide containing scattering nanostructures, and the on-chip coupling scheme enables access to all input spatial modes. By selectively coupling the incident light to open or closed channels of the disordered system, we not only vary the total energy stored inside the system by 7.4 times, but also change the energy density distribution from an exponential decay to a linear decay and to a profile peaked near the center. This work provides an on-chip platform for controlling light-matter interactions in turbid media.PACS numbers: 42.25. Bs, 42.25.Dd, It has long been known that in disordered media there are many fascinating counter-intuitive effects resulting from interferences of multiply scattered waves [1,2]. One of them is the creation of transmission eigenchannels which can be broadly classified as open and closed [3,4]. Existence of high-transmission (open) channels allows an optimally prepared coherent input beam transmitting through a lossless diffusive medium with order unity efficiency. Opposite to that, waves injected to lowtransmission (closed) channels can barely penetrate the medium and are mostly reflected instead. In general, the penetration depth and energy density distribution of multiply scattered waves inside a disordered medium are determined by the spatial profiles of the transmission eigenchannels that are excited by the incident light. The distinct spatial profiles of open and closed channels suggest that selective coupling of incident light to these channels enables an effective control of total transmission and energy distribution inside the random medium [5,6]. Since the energy density determines the light-matter interactions inside a scattering system, manipulating its spatial distribution opens the door to tailoring optical excitations as well as linear and nonlinear optical processes such as absorption, emission, amplification, and frequency mixing inside turbid media. The potential applications range from photovoltaics [7,8] In recent years there have been numerous theoretical and experimental studies on transmission eigenchannels [5,[13][14][15][16][17]. While they can be deduced from the measured transmission matrix [18][19][20][21], it is difficult to directly probe their spatial profiles inside three-dimensional (3D) random media. So far, the open and closed channels are observed only with acoustic wave inside a two-dimensional (2D) disordered waveguide [22], but controlling the energy density distribution has not been realized due to lack of an efficient wavefront modulator for acoustic wave or microwave. The advantage for optical wave is the availability of spatial light modulator (SLM) with many degrees of freedom, however, the commonly used samples in the optics experiment have an open slab geometry, making it impossible to control all input modes due to finite numerical aperture of the imaging...