We have been numerically investigating the thermal-fluid characteristic of sCO2 laminar flowing in a 0.5 mm diameter circular microtube in upward direction. Depending on sharp density change along the tube, this case can be admitted as a compressible flow. When applied pressure-based solver option, the pressure firstly decreases until a certain point along the tube, then increases towards the outlet. So, this situation is impossible. Can this situation result from pressure-based solver? Can it be fixed by using density-based solver? I would like to be grateful if you could share your recommendations. (Boundary conditions: mass flow inlet, pressure outlet (75 bar), constant wall heat flux)