Funkcja Greena

Funkcja Greena jest standardowym narzędziem do rozwiązywania liniowego równania różniczkowego z wyrazem źródłowym, gdy warunki brzegowe lub początkowe są już ustalone. Mówiąc prosto, opisuje ona, jak układ reaguje na jednostkowe źródło punktowe umieszczone w punkcie $\xi$.

Dla operatora liniowego $L$ funkcję Greena $G(x,\xi)$ definiuje się przez problem ze źródłem punktowym

$$L\,G(x,\xi) = \delta(x-\xi)$$

w zmiennej $x$, wraz z tymi samymi warunkami brzegowymi lub początkowymi, które ma spełniać końcowe rozwiązanie. Tutaj $\delta(x-\xi)$ to delta Diraca, czyli dystrybucja reprezentująca jednostkowe źródło skupione w punkcie $\xi$.

Jeśli taka funkcja Greena istnieje i problem jest liniowy, to rozwiązanie równania

$$L\,u(x) = f(x)$$

często można zapisać w postaci

$$u(x) = \int G(x,\xi) f(\xi)\, d\xi.$$

To jest główna idea: raz rozwiązać problem ze źródłem punktowym, a potem złożyć te odpowiedzi punktowe, aby obsłużyć ogólne źródło $f$.

Co oznacza funkcja Greena

Źródło delta sprawia, że ta definicja jest użyteczna. Zamiast od razu rozwiązywać pełne równanie, najpierw pytasz o odpowiedź na możliwie najbardziej skupione wymuszenie.

Dlatego $G(x,\xi)$ można czytać jako „wpływ w punkcie $x$ jednostkowego impulsu przyłożonego w punkcie $\xi$”. Bardziej złożone źródło $f$ traktuje się wtedy jako ciągłą superpozycję źródeł punktowych. Ten krok działa tylko dlatego, że problem jest liniowy.

Dlaczego warunki brzegowe zmieniają odpowiedź

Funkcja Greena nie jest wyznaczona wyłącznie przez operator różniczkowy. Należy do całego problemu, łącznie z warunkami.

Na przykład operator $-u''$ na przedziale $0<x<1$ z warunkami $u(0)=u(1)=0$ ma inną funkcję Greena niż ten sam operator z warunkami brzegowymi Neumanna. Równanie może wyglądać tak samo, ale dopuszczalne rozwiązania się zmieniają, więc jądro też się zmienia.

To jeden z najczęstszych błędów. Nie istnieje jedna funkcja Greena dla „tego równania”, jeśli warunki nie są już częścią sformułowania problemu.

Przykład: problem Dirichleta dla $-u''(x)=f(x)$

Rozważ problem brzegowy

$$-u''(x) = f(x), \qquad 0 < x < 1,$$

z warunkami

$$u(0)=0, \qquad u(1)=0.$$

Dla tego konkretnego problemu funkcja Greena ma postać

$$G(x,\xi) = \begin{cases} x(1-\xi), & x \le \xi, \\ \xi(1-x), & x \ge \xi. \end{cases}$$

Jest ona ciągła względem $x$, znika dla $x=0$ i $x=1$, a jej pierwsza pochodna ma właściwy skok w punkcie $x=\xi$. To właśnie ten skok wytwarza źródło delta w równaniu $-G''(x,\xi)=\delta(x-\xi)$.

Gdy jądro jest już znane, każdy wyraz źródłowy dla tego samego problemu brzegowego korzysta z tego samego wzoru:

$$u(x)=\int_0^1 G(x,\xi) f(\xi)\, d\xi.$$

Wybierzmy prosty wyraz źródłowy:

$$f(x)=1.$$

Wtedy

$$u(x)=\int_0^1 G(x,\xi)\, d\xi = \int_0^x \xi(1-x)\, d\xi + \int_x^1 x(1-\xi)\, d\xi.$$

Obliczmy każdą część:

$$\int_0^x \xi(1-x)\, d\xi = (1-x)\frac{x^2}{2},$$ $$\int_x^1 x(1-\xi)\, d\xi = x\frac{(1-x)^2}{2}.$$

Po zsumowaniu dostajemy

$$u(x)=\frac{x(1-x)}{2}.$$

Szybkie sprawdzenie to potwierdza:

$$u'(x)=\frac{1-2x}{2}, \qquad u''(x)=-1,$$

więc

$$-u''(x)=1=f(x),$$

a wartości brzegowe wynoszą

$$u(0)=0, \qquad u(1)=0.$$

Ten przykład dobrze pokazuje korzyść z tej metody. Nie budujesz całego rozwiązania od nowa dla każdego nowego $f$; zachowujesz to samo jądro i zmieniasz tylko całkę.

Typowe błędy związane z funkcją Greena

1. Traktowanie funkcji Greena tak, jakby działały w ten sam sposób dla problemów nieliniowych. Krok superpozycji zależy od liniowości.
2. Zapominanie, że warunki brzegowe lub początkowe są częścią definicji.
3. Zakładanie, że jądro zawsze musi być symetryczne. Symetria zwykle wymaga dodatkowej struktury, na przykład odpowiedniego ustawienia samosprzężonego.
4. Mylenie funkcji Greena z rozwiązaniem fundamentalnym. Są ze sobą blisko związane, ale nie zawsze są tym samym obiektem.
5. Używanie wzoru całkowego bez sprawdzenia, czy dla danego problemu funkcja Greena w ogóle istnieje.

Gdzie używa się funkcji Greena

Funkcje Greena pojawiają się w równaniach różniczkowych zwyczajnych i cząstkowych wszędzie tam, gdzie liniowy problem jest wymuszany przez wyraz źródłowy. Są powszechne w elektrostatyce, dyfuzji, zagadnieniach falowych, mechanice kwantowej i teorii sprężystości.

W liniowych układach niezmienniczych względem czasu są też ściśle związane z odpowiedzią impulsową. Nazewnictwo zmienia się między dziedzinami, ale idea jest podobna: zrozumieć odpowiedź na jedno skupione wymuszenie, a potem zbudować odpowiedź na ogólne wymuszenie na tej podstawie.

Jak to szybko zapamiętać

Jeśli zapomnisz formalnej definicji, zapamiętaj tę wersję:

$$\text{Funkcja Greena} = \text{odpowiedź na jedno źródło punktowe}.$$

Pełne rozwiązanie powstaje wtedy przez sumowanie takich odpowiedzi punktowych w całej dziedzinie, co właśnie robi całka.

Spróbuj podobnego wyrazu źródłowego

Zachowaj te same warunki brzegowe, ale zamień $f(x)=1$ na $f(x)=x$. Jądro pozostaje takie samo, a zmienia się tylko ostatnia całka:

$$u(x)=\int_0^1 G(x,\xi)\,\xi\, d\xi.$$

To dobry następny krok, bo wyraźnie rozdziela dwa elementy: funkcja Greena należy do ustawienia problemu, natomiast wyraz źródłowy może zmieniać się w zależności od przypadku. Jeśli chcesz pójść dalej, spróbuj własnej wersji z innym źródłem i sprawdź otrzymaną całkę względem warunków brzegowych.