1D solution with 8 inlets
 

Hi everyone,

My problem is a fluid flow simulation in a pipe with 8 inlets. I am writing a 1D solution based in FVM, but I have some problems that I cannot handle.
The flow rate is increased along the pipe.
I have been calculated velocity by continuity equation, and added implicit source at energy equation. The pressure is calculated by Darcy Equation, and the density by the ideal gas model.
I don't know why a have the temperature profile as is shown in the figure uploaded.

Part of my code is:

bcs is a list that stores the control volumes with an inlet.

for energy source:
if CV in bcs:
a_P_S[CV] = flowrate at CV inlet
else:
a_P_S[CV] = 0.

for continuity equation:

if CV in bcs:
rho_U_in = flowrate at CV inlet

rho_U_w =flowrate[CV-1]
['U'][CV] = (rho_U_w - rho_U_in) / (rho[CV]*area[CV])


Can anyone help me?
Thanks in advance