Pole obliczeniowej dynamiki płynów (CFD) poświęcone jest rozwiązywaniu równań Naviera-Stokesa (lub ich uproszczeniu). Podzbiór modeli CFD, oceanicznych i atmosferycznych rozwiązuje liczbowo te same równania dla realistycznych zastosowań. Jakie są różnice i kompromisy między ogólnym podejściem do CFD a zastosowanymi realistycznymi przypadkami?
10
Odpowiedzi:
Atmosfera i ocean mają wysoce rozwarstwione przepływy, w których siła Coriolisa jest głównym źródłem dynamiki. Utrzymanie równowagi geostroficznej jest niezwykle ważne, a wiele schematów numerycznych ma być dokładnie zgodnych (przynajmniej przy braku topografii), aby uniknąć promieniowania energii w falach grawitacyjnych. Ze względu na stratyfikację ograniczenie dyfuzji numerycznej w pionie jest niezwykle ważne i do tego celu często stosuje się specjalne siatki (szczególnie w oceanie). Wiele metod to skutecznie preparaty 2,5-wymiarowe.
W przypadku symulacji klimatu w długich okresach często uważa się, że zachowanie energii i innych strumieni (takich jak sól) ma krytyczne znaczenie dla wyników istotnych statystycznie. Metody, które są mniej dokładne i mają określone artefakty numeryczne, mogą być wybrane w celu uniknięcia tłumienia dynamiki. Należy zauważyć, że długoterminowa dynamika może nie ujednolicić się w skalach kontynentalnych uśrednionych z wielu dziesięcioleci.
Przemysłowe solwery CFD są zwykle używane do przepływów, które są bardziej izotropowe (autentycznie 3D) i często pomijają Coriolisa. Często mają silniejsze siły, a tym samym mniej krytyczne wymagania dotyczące oszczędzania energii. Często występują silne wstrząsy, w których to przypadkach należy stosować nieliniowe dyskretyzacje przestrzenne, mimo że są bardziej rozpraszające.
Ponieważ eksperymenty laboratoryjne można faktycznie przeprowadzić dla większości aplikacji przemysłowych, oprogramowanie podlega większej weryfikacji. Modele pogodowe również podlegają stałej walidacji, ale modele klimatyczne są prawie niemożliwe do walidacji ze względu na stosowaną skalę czasową i nieuniknione nadmierne dopasowanie.
źródło
Jed Brown opisał tradycyjne podejście stosowane w modelach mezoskalowych i większych. W rzeczywistości w mikroskali wiele modeli atmosferycznych jest bardzo zbliżonych do tradycyjnych kodów CFD, stosuje podobne dyskretyzacje o skończonej objętości, podobne siatki 3D, w których pion jest traktowany podobnie jak poziomo i tak dalej. W zależności od rozdzielczości nawet funkcje takie jak budynki są rozwiązywane przy użyciu tych samych metod znanych z inżynieryjnego CFD, takich jak metody zanurzenia na granicy lub siatki montowane na nadwoziu.
Możesz napotkać wszystkie techniki dyskretyzacji, które znasz z inżynierskiego CFD, takie jak skończone różnice, skończone objętości, pseudo-spektralne, a nawet skończone elementy. Te same metody korekcji ciśnienia (krok ułamkowy) są często stosowane do rozwiązywania nieściśliwych równań Naviera-Stokesa (z Boussinesq lub anelastycznymi określeniami wyporu).
Oczywiście powszechnie stosuje się różne parametryzacje strumieni ciepła i pędu w pobliżu powierzchni, biorąc pod uwagę specyfikę interakcji ląd-powierzchnia, takich jak podobieństwo Monina-Obuchowa lub inne relacje półempiryczne.
Cała metoda symulacji dużych wirów (LES), obecnie bardzo popularna w inżynierii, faktycznie wywodzi się z meteorologii warstw granicznych. Powiedziałbym nawet, że wielu modelarzy atmosferycznych w tej skali nie zawahałoby się nazywać swojej pracy CFD.
W wielu (ale nie wszystkich) aplikacjach musisz również dodać siłę Coriolisa. Schematy nie muszą być dobrze wyważone, jest to tylko jedna dodatkowa siła objętości. Jeśli obliczymy również takie procesy, jak tworzenie się chmur, opady i promieniowanie, sprawy staną się bardziej skomplikowane, ale to samo dotyczy modeli inżynieryjnych, które rozwiązują kinetykę reakcji, spalanie i podobne.
Ta klasa modeli obejmuje również te uwzględniające interakcje ocean-atmosfera, o które prosiłeś, patrz na przykład https://ams.confex.com/ams/pdfpapers/172658.pdf
źródło
Różnica między oprogramowaniem do prognozowania pogody a „zwykłym solwerem CFD” polega na tym, jak prognozowanie pogody działa z przejściem wody. Woda jest traktowana jako drugi komponent, więc model staje się trójwymiarowy z 2 komponentami.
źródło