Due to layered structure, the electronic conduction in high-T_c cuprate superconductors is highly anisotropic. Therefore, in the present work, we have made an attempt to study the anisotropy in supercurrent density in layered cuprate superconductors likeYBa₂Cu₃O_6+x. We have considered tight binding bilayered Hubbard Hamiltonian. The situation considered here for the bilayered cuprates is similar to a SIS junction. Using Green function technique, we have got the expressions for superconducting order parameters, carrier density and supercurrent densities.The numerical analysis shows that the anisotropy in supercurrent density depends on the temperature and various other microscopic parameters of the model Hamiltonian. The anisotropy in supercurrent density increases with temperature and diverges as system reaches towards the superconducting transition temperature. On increasing the carrier density within the plane, the anisotropy in supercurrent density decreases because around optimal doping region, bilayer system(YBaCuO) behaves like almost isotropic three dimensional superconductors. On increasing the interlayer attractive interaction, the anisotropy in supercurrent density decreases. It is also shown that when interlayer attractive interaction (W) is just equal to the intralayer pairing interaction (U), the anisotropy in supercurrent density vanishes i.e. the bilayer Quasi-2D system becomes almost isotropic(three dimensional) in supercurrent flow in the superconducting state. We have shown, using tight binding bilayer Hubbard model within BCS meanfield formalism, that the anisotropy in supercurrent density depend on various microscopic intra and interlayer interactions as well as on the thermodynamics of superconducting order parameter in an essential way. Finally we have compared our theoretical results of anisotropy in supercurrent density with the existing experimental findings in bilayered cuprate superconductors.

PACS: 74.50.+r, 74.72.Bk, 74.81.Fa.