Even though the existence of meridional currents in stars has been known for quite a long time (Eddington 1925, Vogt 1925), its exact structure as well as its influence on stellar evolution is still unclear. Some authors concentrated on finding the exact shape of meridional circulation in a rotating star, while others tried to model its effect on the chemical distribution in the interior. In all studies performed so far however, meridional circulation is considered in an asymptotic regime in which the advection of entropy by the meridional currents is supposed to balance exactly the source term of the non-zero radiative flux divergence. Other terms could however be added to that asymptotic regime which could turn out to dominate the transport of chemicals. We wish to present here preliminary results of 3D numerical simulation attempted to tackle this problem.