Документация

Проблемное определение

В этом разделе описываются проблему, которая будет решена. Это также обеспечивает необходимые уравнения и известные значения.

Поток строки Fanno является моделированием совершенного потока газа через канал постоянной области, который не изменяется со временем и который является адиабатой. Стенное трение является основным механизмом для изменения переменных потока. Этот пример смотрит на поток Fanno воздуха, вводящего канал 3 сантиметра диаметром, который 45 сантиметров длиной в Числе Маха 0,6. Условия во входном отверстии, также названном станцией 1, являются статическим давлением 150 килопаскалей и статической температурой 300 кельвинов. Канал принят, чтобы иметь коэффициент трения 0,02. Вычислите Число Маха, статическое давление и статическую температуру в выходе канала или станции 2.

ductPicture = astfrictionduct;

Данная информация в проблеме:

ductLength       = 0.45;    % Length of the duct [m]
diameter         = 0.03;    % Diameter of the duct [m]
inletMach        = 0.6;     % Mach number at the duct inlet [dimensionless]
inletPressure    = 150;     % Static pressure at the duct inlet [kPa]
inletTemperature = 300;     % Static temperature at the duct input [K]
frictionCoeff    = 0.02;    % Duct friction coefficient [dimensionless]

Жидкость является воздухом, который имеет следующее отношение удельной теплоемкости.

k = 1.4;    % Specific heat ratio [dimensionless]

Понимание параметра Fanno

Параметр Fanno является безразмерным количеством, которое указывает, сколько будет иметь трение влияния, в то время как жидкость течет через канал. Для данного канала параметр Fanno задан как

Где

Для кругового канала примите, что гидравлический диаметр является внутренним диаметром канала. Коэффициент трения, f, дан следующим выражением:

где

Отметьте, этот пример использует соглашение, в котором фактор четыре не видим в параметре Fanno. Это соглашение задает коэффициент трения как четыре раза трение кожи по динамическому давлению. Трение будет иметь больший эффект на поток в длинном канале, чем в коротком канале, потому что потоку препятствует больше стенного трения поверхности. Кроме того, трение является более доминирующим, когда канал является узким. Это вызвано тем, что пограничный слой влияет на больший фрагмент потока вдоль стен чем тогда, когда ширина канала является большой.

Трение является энергетической потерей, которая генерирует энтропию (необратимость) в системе. Увеличение энтропии заставляет поток склоняться к дросселируемому условию (Мах = 1). Дросселируемое условие происходит, если длина канала достаточно длинна. Для данного Числа Маха и отношения удельной теплоемкости, ссылочный параметр Fanno для дросселируемого потока

где

Используя параметр Fanno и FLOWFANNO, чтобы решить для свойств потока во входном отверстии

Этот пример обеспечивает длину канала, диаметр канала и коэффициент трения. Поэтому фактический параметр Fanno канала может быть вычислен можно следующим образом:

fannoParameter = frictionCoeff * ductLength / diameter;

Этот пример также обеспечивает отношение удельной теплоемкости и Число Маха. Это позволяет вам вычислить ссылочный параметр Fanno для входного условия, входного отношения температуры и вставить отношение давления. Используйте функцию flowfanno от Aerospace Toolbox:

[~, inletTempRatio, inletPresRatio, ~, ~, ~, inletFannoRef] = flowfanno(k, inletMach);

где:

Длина входной ссылки, параметр Fanno, также названный, вставил макс. длину, является длиной, что канал должен дросселировать поток для данного входного условия. Если фактическая длина является меньше, чем входное отверстие макс. длина, расширение канала необходимо для дросселируемого потока. Этот дросселируемый поток соответствует ссылочному параметру Fanno для выхода. Поскольку диаметр и коэффициент трения даны в проблемном операторе, только длины отличаются по следующему уравнению для ссылки выхода параметр Fanno:

outletFannoRef = inletFannoRef - fannoParameter;

Вычисление свойств потока при выходе с функцией FLOWFANNO

Затем, используйте flowfanno, чтобы вычислить отношения потока в станции выхода. Третий вход, 'fannosub', указывает, что второй вход, outletFannoRef, является дозвуковым входом параметра Fanno.

[outletMach, outletTempRatio, outletPresRatio] = flowfanno(k, outletFannoRef, 'fannosub');

Используйте температурные отношения, найденные во входном отверстии и выходе, чтобы вычислить температуру и давление при выходе. Ссылочные условия являются тем же самым в обеих станциях, потому что канал изолируется. Кроме того, примите, что эффекты трения действуют на обе станции таким же образом. В результате мы имеем

Поэтому температура при выходе и давлении при выходе

outletTemperature = inletTemperature / inletTempRatio * outletTempRatio;

outletPressure = inletPressure / inletPresRatio * outletPresRatio;

Значения, которые мы хотели вычислить,

outletMach          % [dimensionless]
outletTemperature   % [K]
outletPressure      % [kPa]

outletMach =

    0.7093


outletTemperature =

  292.2018


outletPressure =

  125.2332

Поскольку строка Fanno течет, где входной поток является дозвуковым, температура и давление всегда уменьшаются через канал. Для всего Fanno случается поток строки, Число Маха придвигается поближе к одному.

close(ductPicture)

Ссылка

[1] Джеймс, J. E. A. "газовая динамика, второй выпуск", Allyn and Bacon, Inc, Бостон, 1984.

%#ok<*NOPTS>