Warunek Couranta-Friedrichsa-Lewy’ego

Z Wikipedii, wolnej encyklopedii
Stabilność rozwiązania równania z warunkiem początkowym i brzegowym dla różnych kroków czasowych.

Warunek Couranta-Friedrichsa-Lewy’ego (warunek CFL) – matematyczny warunek zbieżności numerycznych metod rozwiązywania pewnych równań różniczkowych cząstkowych (zwłaszcza równań hiperbolicznych). Pojawia się przy analizie stabilności jawnych metod numerycznych dla zagadnień zależnych od czasu (lub równoważnych im).

Warunek CFL głosi, że długość kroku czasowego używanego w przybliżeniu numerycznym równania różniczkowego nie może przekroczyć pewnej wielkości granicznej, gdyż w przeciwnym wypadku metoda numeryczna straci stabilność, a uzyskane za jej pomocą rozwiązanie będzie diametralnie odbiegać od rozwiązania rzeczywistego.

Nazwa warunku pochodzi od nazwisk trzech niemieckich matematyków: Richarda Couranta, Kurta Friedrichsa i Hansa Lewy’ego, którzy wyprowadzili go w 1928 r.[1][2]

Definicja[edytuj | edytuj kod]

Warunek CFL zwykle dotyczy równań różniczkowych z członem adwekcyjnym opisującym propagację fal, np. równania adwekcji, równania Burgersa, równania Naviera-Stokesa. W przypadku jednowymiarowym, warunek CFL wyraża się wzorem

gdzie:

– (maksymalna) prędkość fali (L/T),
– numeryczny krok czasowy (T),
– wartość stałej sieci w modelu numerycznym (L),

a wartość stałej zależy od postaci rozwiązywanego równania oraz zastosowanego modelu numerycznego, nie zależy zaś od ani W bardzo wielu zastosowaniach praktycznych

W przypadku dwuwymiarowym warunek CFL przyjmuje postać

Znaczenie fizyczne[edytuj | edytuj kod]

Fizyczne znaczenie warunku CFL można przedstawić następująco. Jeżeli oryginalne równanie różniczkowe opisuje propagację fali, to w modelu numerycznym, w którym przestrzeń ciągłą przybliżono dyskretną siatką punktów, fala musi przechodzić pomiędzy sąsiednimi punktami siatki w czasie nie dłuższym niż czas potrzebny rzeczywistej fali na pokonanie tej samej odległości:

Stąd wniosek, że wraz ze zmniejszaniem odległości między punktami siatki obliczeniowej redukcji ulega też maksymalna wartość kroku czasowego używanego w symulacji. Warunek CFL oznacza, że tzw. obszar zależności rozwiązania numerycznego musi zawierać cały obszar zależności rozwiązania analitycznego, gdyż tylko wtedy algorytm numeryczny ma dostęp do wszystkich danych niezbędnych do wyznaczenia poprawnego rozwiązania zagadnienia analitycznego. Z warunku CFL wynika, że każdej próbie zmniejszenia błędu obcięcia modelu numerycznego poprzez zagęszczenie siatki obliczeniowej powinno towarzyszyć odpowiednie zmniejszenie kroku czasowego, co dodatkowo zwiększa czasochłonność obliczeń.

Liczba Couranta[edytuj | edytuj kod]

Bezwymiarowy parametr

zwany jest liczbą Couranta.

W przypadku trójwymiarowym liczba Couranta zdefiniowana jest wzorem

Przypisy[edytuj | edytuj kod]

  1. R. Courant, K. Friedrichs, H. Lewy, Über die partiellen Differenzengleichungen der mathematischen Physik, Mathematische Annalen, vol. 100, nr. 1, s. 32–74, (1928).
  2. R. Courant, K. Friedrichs, H. Lewy, On the partial difference equations of mathematical physics, IBM Journal, Marzec 1967, s. 215–234. Jest to angielskie tłumaczenie oryginalnej pracy z 1928 r., wersja elektroniczna: tutaj.

Bibliografia[edytuj | edytuj kod]

  • John C. Tannehill, Dale A. Anderson i Richard H. Pletcher, Computational Fluid Mechanics and Heat Transfer (Second Edition), Francis & Taylor, Philadelphia, 1997, ISBN 1-56032-046-X.