Flow characteristics for a wide range of Reynolds number up to turbulent gas flow regime, including flow choking were numerically investigated with a microtube discharged into the atmosphere. The numerical methodology is based on the Arbitrary-Lagrangian-Eulerian (ALE) method. The LB1 turbulence model was used in the turbulent flow case. Axis-symmetric compressible momentum and energy equations of an ideal gas are solved to obtain the flow characteristics. In order to calculate the underexpanded (choked) flow at the microtube outlet, the computational domain is extended to the downstream region of the hemisphere from the microtube outlet. The back pressure was given to the outside of the downstream region. The computations were performed for adiabatic microtubes whose diameter ranges from 10 to 500 mu m and whose aspect ratio is 100 or 200. The stagnation pressure range is chosen in such a way that the flow becomes a fully underexpanded flow at the microtube outlet. The results in the wide range of Reynolds number and Mach number were obtained including the choked flow. With increasing the stagnation pressure, the flow at the microtube outlet is underexpanded and choked. Although the velocity is limited, the mass flow rate (Reynolds number) increases. In order to further validate the present numerical model, an experiment was also performed for nitrogen gas through a glass microtube with 397 mu m in diameter and 120 mm in length. Three pressure tap holes were drilled on the glass microtube wall. The local pressures were measured to determine local values of Mach numbers and friction factors. Local friction factors were numerically and experimentally obtained and were compared with empirical correlations in the literature on Moody's chart. The numerical results are also in excellent agreement with the experimental ones. (C) 2018 Elsevier Ltd. All rights reserved.