none3We propose a model to describe lava tube formation, considering lava as a Newtonian fluid moving downslope in a rectangular channel. We obtain flow velocity using an analytical steady-state solution of the Navier-Stokes equation. Shear stress is also calculated from velocity for a Newtonian incompressible and isotropic fluid. A 2-D model with heat flux assigned at the upper surface is introduced to describe lava cooling by radiation into the atmosphere and to obtain the flow temperature. Lava crust is considered as a plastic body, and its rheology is described through the introduction of the yield strength as a function of temperature. It describes the capacity of crustal structure to prevent shear deformation for lower shear stress v...