089 4

p

[mPq] 5-

ł [*owk]

Rys. 3.6. Przebieg ciśnienia gazu w przestrzeni sprężania p_c i rozprężania p_£

metrów roboczych silnika i jego podzespołów. Dane te są niekiedy bardzo szczegółowe, stąd ich uzyskanie musi być poprzedzone projektowaniem wstępnym oraz obliczeniami optymalizacyjnymi silnika prostszym modelem drugiego rzędu. Niektóre parametry początkowe, jak: prędkość, temperatura i masa gazu w poszczególnych objętościach kontrolnych należy przyjąć arbitralnie, gdyż nie są one znane na tym etapie obliczeń. Stąd, po zrealizowaniu pełnego cyklu obliczeń należy sprawdzić warunek poprawności otrzymanych wyników. Za odpowiednie kryterium przyjęto stan zbilansowania na zero całkowitego ciepła wymienionego między gazem a wkładem regeneratora w czasie realizacji pełnego obiegu. Na ogół uzyskanie takiego stanu wymaga wykonania kilkudziesięciu pełnych cykli obliczeń. Stwierdzonó przydatność stosowania metody wyznaczania wartości początkowych nowego cyklu opartej na średniej arytmetycznej z wartości początkowej i końcowej poprzedniego cyklu obliczeń. Wyjątek stanowi temperatura wkładu regeneratora, dla której istnieje możliwość wprowadzenia zmiennego współczynnika korekcyjnego 5, przyśpieszającego moment spełnienia końca obliczeń

T»T +5(T    - T )

(3.70)

gdzie: T - temperatura wkładu regeneratora na początku obiegu, T - tempera-

mO    ml

tura wkładu regeneratora na końcu obiegu, 5 - współczynnik korekcyjny z przedziału (0+1).

Niektóre parametry określane w każdym kroku obliczeń, zwłaszcza ciśnienia wewnątrz objętości kontrolnych, wykazują tendencje do oscylacji. Oscylacje te należy przypisać bardziej samej metodzie obliczeń niż faktycznym zjawiskom fizycznym. W celu wyeliminowania ich wpływu na obliczenia zastosowano specjalną metodę stabilizacji, polegającą na wykorzystaniu obliczonych wartości ciśnienia wg modelu do wyznaczenia wartości na granicach objętości kontrolnych i osta-

92