Analiza porównawcza schematów rozwiązywaczy dla przepływów ściśliwych w aerodynamice dżetów niedorozprężonych ze weryfikacją Schlieren

Dlaczego szybkie dżety gazowe mają znaczenie

Od dysz rakietowych po przemysłowe wyloty i nawet erupcje wulkaniczne — szybkie dżety gazowe kształtują rozkład ciągu, hałasu i ciepła w otaczającym powietrzu. Dokładne przewidzenie zachowania takich dżetów jest trudne, ponieważ fale sprężonego i rozprężonego gazu tworzą złożone wzory, które standardowe narzędzia komputerowe mają trudność śledzić. W tym badaniu stawiamy praktyczne pytanie: czy tańszy typ rozwiązującego równania, zwykle przeznaczony do wolniejszych przepływów, poradzi sobie z tymi wymagającymi zadaniami bez utraty dokładności?

Dwa różne sposoby rozwiązania tego samego przepływu

Inżynierowie często polegają na obliczeniowej dynamice płynów (CFD), aby badać projekty przed ich wykonaniem. Dla szybkich przepływów gazu tradycyjnie stosowano solver „oparty na gęstości”, ponieważ jest on dostosowany do obsługi gwałtownych zmian ciśnienia i gęstości wokół fal uderzeniowych. Solver „oparty na ciśnieniu”, przeciwnie, jest powszechnie używany dla przepływów wolnozmiennych, bliskich nieściśliwości, i jest znany z większej oszczędności. W ostatnich latach niektórzy badacze zgłaszali, że solvery oparte na ciśnieniu mogą też działać dobrze przy wyższych prędkościach, ale ich wydajność dla silnie niedorozprężonych dżetów — gdzie gaz opuszcza dyszę przy znacznie wyższym ciśnieniu niż powietrze otoczenia — pozostawała niepewna.

Budowa dżetu w laboratorium i fotografowanie niewidocznych fal

Aby sprawiedliwie przetestować oba podejścia, autorzy skonstruowali prosty, lecz starannie kontrolowany eksperyment. Sprężone powietrze wypuszczano ze zbiornika przez małą dyszę skupiającą do nieruchomego powietrza w pomieszczeniu. Poprzez regulację ciśnienia zasilania ustawiono cztery stosunki ciśnień między zbiornikiem a atmosferą, od nieco powyżej warunku „dławienia”, gdy gaz opuszcza dyszę z prędkością dźwięku, aż po wyraźnie niedorozprężone dżety z silnymi wzorami fal. Aby zobaczyć te w przeciwnym razie niewidoczne struktury, wykorzystano układ optyczny Schlieren: punktowe źródło światła, zwierciadło wklęsłe, żyletkę i aparat cyfrowy tak ustawione, że drobne zmiany gęstości powietrza pojawiają się jako jasne i ciemne pasy. Te obrazy dostarczyły wizualnej mapy komórek falowych tworzących się wzdłuż strugi za dyszą.

Postawienie modeli komputerowych na próbę

Równolegle te same warunki dyszy i przepływu odtworzono w kodzie CFD ANSYS Fluent. Badacze użyli identycznych siatek, warunków brzegowych i modeli turbulencji dla obu solverów, zmieniając jedynie zasadniczy schemat rozwiązania. Śledzili ciśnienie, prędkość i parametry turbulencji wzdłuż i wokół dżetu oraz sprawdzali, jak szybko każdy solver zbiega do stabilnego wyniku. Oba narzędzia odtworzyły kluczowe cechy: powstanie rdzenia potencjalnego przy warunku dławienia, pojawienie się diamentowych komórek falowych gdy dżet stał się niedorozprężony, oraz sposób, w jaki komórki te wydłużały się i wzmacniały wraz ze wzrostem stosunku ciśnień. Pierwsza główna fala uderzeniowa, znana jako tarcza Macha, pojawiła się niemal w tych samych lokalizacjach w obu symulacjach i w opublikowanych danych eksperymentalnych, z różnicami w szczytowych prędkościach na poziomie tylko kilku procent.

Gdzie solvery się zgadzają, a gdzie się różnią

Wzdłuż linii środkowej dżetu oba solvery wygenerowały niemal identyczne oscylacje ciśnienia i liczby Macha, odzwierciedlające powtarzalny wzór sprężania i rozprężania. Różnice pojawiały się dalej w dół strugi, gdzie dżet zwolnił do prędkości poddźwiękowych i dominowała turbulencja. Tam solver oparty na gęstości miał skłonność do przewidywania wyższej kinetycznej energii turbulencji niż solver oparty na ciśnieniu, zgodnie ze znanymi ograniczeniami schematów opartych na gęstości przy niskich prędkościach. Mimo to oba solvery dały prawie taki sam współczynnik wydajności przepływu (discharge coefficient), miarę efektywności przejścia masy przez dyszę, z różnicami poniżej 0,2 procent. Pod względem nakładu obliczeniowego solver oparty na ciśnieniu osiągnął ten sam rygorystyczny cel zbieżności w około 22 procent krótszym czasie CPU.

Co to oznacza dla zastosowań praktycznych

Dla inżynierów, którzy muszą symulować szybkiego dżety z relatywnie prostych dysz, badanie to przynosi uspokajającą wiadomość. Solver oparty na ciśnieniu, choć pierwotnie zaprojektowany do wolniejszych, nieściśliwych przepływów, potrafi uchwycić główne wzory falowe i metryki wydajności dżetów nadźwiękowych niedorozprężonych z dokładnością porównywalną do wyspecjalizowanego solvera opartego na gęstości, przy jednoczesnym mniejszym zużyciu czasu obliczeniowego. Autorzy zastrzegają, że schematy oparte na gęstości wciąż mają przewagi w niektórych reżimach, zwłaszcza tam, gdzie występują bardzo silne uderzenia i złożona chemia. Jednak dla wielu praktycznych problemów dysz ta tańsza metoda oparta na ciśnieniu może być wiarygodnym wyborem, pozwalając projektantom badać więcej opcji bez wygórowanych kosztów symulacji.

Cytowanie: Alsaedi, S.S., Al-Sadawi, L.A., Al-Haddad, L.A. et al. Comparative analysis of compressible solver schemes for underexpanded jet aerodynamics with Schlieren validation. Sci Rep 16, 15724 (2026). https://doi.org/10.1038/s41598-026-44651-w

Słowa kluczowe: dżet niedorozprężony, dysza nadźwiękowa, obrazowanie Schlieren, solver CFD, fale uderzeniowe