Hybrydowe schematy przestrzenne dla CFD: jakaś wada mieszania kontra zamiana?

9

Czy oprócz dodatkowych kosztów obliczeniowych związanych z koniecznością obliczenia obu strumieni w określonym regionie, czy jest jakaś wada łącząca dwie oceny strumienia dla schematu hybrydowego w metodzie o skończonej objętości? Ocena strumienia wyglądałaby następująco:

faja+12)=Λja+12)faja+12)do+(1-Λja+12))faja+12)u

Przełącznik oparty jest na czujniku gradientu ciśnienia i / lub gęstości, w zależności od zastosowania. fado jest systemem centralnym (McCormack, kompaktowy, ...) i faujest schematem podmuchu wiatru, jak podział różnicy strumienia z rekonstrukcją MUSCL. Czy są jakieś problemy dotyczące liczb, konserwatywnych właściwości, jeśli mieszam dwa schematy za pomocą funkcji ciągłej dlaΛ w przeciwieństwie do zwykłego przełączania między schematami za pomocą Λ wyceniony jako 0 lub 1?

FrenchKheldar
źródło
Klasycznym przykładem błędnych przełączników jest schemat Murmana-Cole'a z 1970 r. Dla równań potencjałów transonicznych (ustalony w 1973 r.). Jeśli zmiana lub łączenie nie daje teleskopowej sumy, naruszasz ochronę. Testowanie jest proste.
tpg2114
Jak rozwiążesz równania? Implikowana dyskretyzacja (pseudo) czasu? Przełączanie może mieć negatywny wpływ na konwergencję iteracyjną.
Chris
Byłoby to przy użyciu wyraźnego schematu czasowego, predyktora-korektora lub podstawowej Runge-Kutty.
FrenchKheldar

Odpowiedzi:

4

Podejście, które zastosujesz, zachowa ochronę w obu przypadkach. Istnieją inne oczywiste podejścia, które nie są konserwatywne i mogą powodować problemy.

Możliwe jest (a nawet prawdopodobne), że stracisz porządek dokładności w regionie, w którym się przełączasz, jeśli zbadasz lokalny błąd obcięcia. Ale zazwyczaj ten błąd jest lokalizowany, więc błąd globalny nadal ma oczekiwaną kolejność. Z mojego doświadczenia wynika, że ​​zobaczysz zasadniczo to samo zachowanie, niezależnie od tego, czy używasz twardego przełącznika, czy regionu przejściowego.

Mam rękopis (mniej więcej) na ten temat: Analiza błędów jawnych partycjonowanych schematów Runge-Kutta dla praw ochrony .

Byłbym bardzo zainteresowany wiedząc, co widzisz, gdy wypróbowuję te dwa podejścia, jeśli różni się to od tego, co sugeruję.

David Ketcheson
źródło
super, miałem nadzieję, że będziesz w stanie to skomentować.
Aron Ahmadia
Czy istnieje teoretyczny sposób sprawdzenia, podobny do sum teleskopowych?
tpg2114,