Hi there,(adsbygoogle = window.adsbygoogle || []).push({});

I am working on a project in a subject called "Multiphase Flow in Pipes." We are to make a "simulator", describing the pressure (step-wise) in a pipe inserted vertically into water. Air is injected at the bottom of the pipe with a given mass rate, m, and for simplicity, ideal gas, no slip conditions, and a given QL (liquid flow rate) are assumed. The pressure-losses taken into account are hydrostatic and friction only. (I can explain further on request, but my problem is, I hope, code-related, and not physics :P)

I started writing what I thought was a fairly simple piece of "code", and decided to test my program with one given m and QL. However, based on my graph, it looks as though the pressure-drop upwards in the pipe is constant, and does not change, no matter what m and QL I try.

I think I may have done something wrong in my "for" loop, but since I am very inexperienced with Matlab, I thought I would ask if anyone might give me a small pointer..

Code follows:

------------------------------

clf

% First, we set up the main constants

g = 9.81;

h = 50;

D = 0.15;

rhoL = 1000;

P(1) = 1*10^5 + (rhoL*h*g);

Ptop = 1.013*10^5;

Pref = 1.013*10^5;

nu = 1.5*10^-5;

viscL = 0.001;

rhoref = 1.2;

% Then we start entering the equations necessary to iterate.

dH = 0.5;

H = 0:dH:h;

NH = length(H);

m = 1;

ql = 0.001;

for i = 2:NH % Loop, upwards in the well.

% Calculate gas density and friction pressure drop

rhoG = rhoref * (P(i-1)/Pref); % density of gas

qg = (m/rhoG); % volumetric flow rate of gas

EG = (qg)/(qg + ql); % gas fraction

EL = (ql)/(qg + ql); % liquid fraction

rhom = (EL * rhoL) + (EG * rhoG); % density of gas/liquid mixture

viscm = (EL * viscL) + (EG * rhoG * nu); % viscosity of gas/liquid mixture

U = (ql+qg)/((pi/4)*D^2); % flow speed

phyd = rhom * g; % hydrostatic pressure loss

pfrict = (4/D) * 0.046 * ((rhom*U*D)/(viscm))^(-0.2); % frictional pressure loss

ptot = phyd + pfrict; % total pressure loss

P(i) = P(i-1) - (ptot*dH)

end

figure(1)

plot(H,P)

---------------------------------------

Thank you very much for your time.

Yours Sincerely,

André Røsbak

**Physics Forums | Science Articles, Homework Help, Discussion**

Dismiss Notice

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Matlab iteration issue (I think :P)

**Physics Forums | Science Articles, Homework Help, Discussion**