PL 
EnglishDomyślny
Українська〰️ Propagacja fal w ośrodku lepkosprężystym
Fale w lepkosprężystym ośrodku Maxwella. Poniżej częstotliwości przejściowej ω_c = G/η dyfundują, powyżej — propagują się. Obserwuj, jak krzywe dyspersji i tłumienie zmieniają się z czasem relaksacji.
Jak to działa
Na lewej granicy wymuszane jest sinusoidalne przemieszczenie. Symulacja propaguje jednowymiarową falę ścinającą przez lepkosprężysty ośrodek Maxwella za pomocą przesuniętego schematu różnic skończonych prędkość-naprężenie (Virieux 1986): prędkość znajduje się w węzłach siatki, a naprężenie w węzłach połówkowych, więc obie pochodne przestrzenne są wycentrowane nad jedną komórką, a szachownicowy tryb nieparzysto-parzysty schematu skolokowanego nie może się pojawić. Podkrok spełnia warunek CFL c·Δt/Δx ≤ 1, gdzie c = √(G/ρ), a człon relaksacji σ/τ jest całkowany dokładnie, więc schemat pozostaje stabilny nawet gdy τ jest znacznie krótsze niż krok czasowy. Model Maxwella daje zależny od częstotliwości zespolony moduł ścinania G*(ω) = G·iωτ/(1+iωτ), który wyznacza zespoloną liczbę falową k(ω) = ω·√(ρ/G*(ω)).
Przy niskich częstotliwościach (ω << ωc = 1/τ) dominuje część urojona k — zaburzenia zanikają w obrębie ułamka długości fali (reżim dyfuzyjny). Przy wysokich częstotliwościach (ω >> ωc) dominuje część rzeczywista k — spójne fale rozchodzą się z prędkością fazową ω/Re(k), zbliżającą się do c = √(G/ρ) (reżim sprężysty). Obszar ma długość 20 jednostek, co przy domyślnych ustawieniach odpowiada mniej więcej trzem długościom fali i 4,3 długościom tłumienia, więc fala widocznie zanika podczas przechodzenia przez ośrodek, zamiast wygasać przy źródle. Przerywany czerwony znacznik na panelu przemieszczenia znajduje się w 1/Im(k), gdzie obwiednia spadła do 1/e. Dolny panel pokazuje krzywą dyspersji na żywo.
ω_c = 1/τ = G/η [częstotliwość przejścia]
G*(ω) = G · iωτ / (1 + iωτ)
k = ω √(ρ / G*(ω))
Najczęściej zadawane pytania
Czym jest materiał lepkosprężysty?
Materiał lepkosprężysty wykazuje zarówno zachowanie lepkie (podobne do cieczy), jak i sprężyste (podobne do ciała stałego), w zależności od skali czasowej odkształcenia. Przykłady obejmują stopione polimery, żele biologiczne i „głupi glut”.
Czym jest lepkosprężysty model Maxwella?
Model Maxwella przedstawia lepkosprężystość jako sprężynę (moduł sprężystości G) i tłumik (lepkość η) połączone szeregowo. Naprężenie relaksuje wykładniczo z czasem relaksacji τ = η/G.
Czym jest czas relaksacji w lepkosprężystości?
Czas relaksacji τ = η/G to czas potrzebny, by naprężenie zanikło do 1/e swojej wartości początkowej po narzuconym skokowym odkształceniu. Oddziela on reżim sprężysty (szybki) od reżimu lepkiego (wolny).
Dlaczego fale dyfundują przy niskiej częstotliwości w cieczy Maxwella?
Przy częstotliwościach poniżej ωc = G/η = 1/τ materiał ma czas na relaksację i zachowuje się jak ciecz lepka. Zaburzenia zanikają dyfuzyjnie, zamiast propagować się jako spójne fale.
Dlaczego fale propagują się przy wysokiej częstotliwości w cieczy Maxwella?
Przy częstotliwościach powyżej ωc materiał nie zdąży relaksować, zanim przejdzie fala, i zachowuje się sprężyście, wspierając propagujące fale ścinające z prędkością c = √(G/ρ).
Jaka jest relacja dyspersji dla fal lepkosprężystych?
Dla cieczy Maxwella: k² = ρω²/G·(1 + i/(ωτ)). Część rzeczywista daje prędkość fazową, część urojona daje tłumienie przestrzenne.
Czym jest tłumienie w propagacji fal?
Tłumienie to zanik amplitudy fali z odległością: amplituda maleje jako exp(−αx), gdzie α = Im(k). W cieczy Maxwella zwiększenie lepkości η wydłuża czas relaksacji τ = η/G, co czyni ośrodek bardziej sprężystym, więc α maleje, a długość tłumienia 1/Im(k) → 2τ·√(G/ρ) rośnie. (To model Kelvina-Voigta, z tłumikiem połączonym równolegle, ma tłumienie rosnące wraz z lepkością.) Odczyt Dł. tłum. 1/Im(k) pokazuje to bezpośrednio.
Dlaczego prędkość fazowa nie jest po prostu równa √(G/ρ)?
√(G/ρ) to prędkość fali ścinającej dla nierelaksowanego (nieskończonej częstotliwości) ośrodka, a ośrodek Maxwella osiąga ją tylko przy ωτ → ∞. Przy skończonym ωτ ośrodek jest miększy niż G, więc rzeczywista prędkość fazowa to ω/Re(k), z k wyznaczonym z relacji dyspersji Maxwella, i spada poniżej √(G/ρ) — zbliżając się do √(2ωη/ρ) w granicy dyfuzyjnej. Odczyt Prędkość fazowa podaje właśnie ω/Re(k), tę samą wielkość wykreślaną w panelu dyspersji.
Jakie są praktyczne przykłady efektów fal lepkosprężystych?
Zjawiska fal lepkosprężystych są istotne w tłumieniu sejsmicznym w skałach osadowych, ultradźwiękach w tkance biologicznej, pochłanianiu wstrząsów w piankach polimerowych oraz tłumieniu drgań w strukturach kompozytowych.
Czym jest liczba Deborah?
Liczba Deborah De = τ/tobs porównuje czas relaksacji materiału ze skalą czasową obserwacji. De >> 1 oznacza zachowanie sprężyste; De << 1 oznacza zachowanie lepkie.
Czym model Kelvina-Voigta różni się od modelu Maxwella?
Model Kelvina-Voigta umieszcza sprężynę i tłumik równolegle, dając odpowiedź pełzania, ale bez relaksacji naprężenia. Model Maxwella umieszcza je szeregowo, dając relaksację naprężenia, ale nieograniczone pełzanie.
O tej symulacji
Solwer różnic skończonych typu przesunięta prędkość-naprężenie krokuje naprężenie na komórkach połówkowych i prędkość na węzłach 400-węzłowej siatki 1D, wykorzystując prawo konstytutywne Maxwella dσ/dt = G(dv/dx) − σ/τ i wymuszany sinusoidalnym przemieszczeniem na lewej krawędzi. Ponieważ czas relaksacji τ = η/G wyznacza częstotliwość przejścia ωc = G/η, ten sam ośrodek może wyglądać jak dyfundująca ciecz albo dzwoniące ciało sprężyste, zależnie wyłącznie od tego, jak szybko się nim potrząsa. Obszar jest dobrany względem własnej długości tłumienia ośrodka 1/Im(k), więc przy typowych ustawieniach fala przechodzi przez kilka długości fali, jednocześnie zanikając — można zobaczyć tłumienie, zamiast jedynie je wywnioskować.
🔬 Co pokazuje
Wymuszaną falę przemieszczenia w górnym panelu, z przerywanym znacznikiem przy 1/Im(k), gdzie obwiednia spadła do 1/e, oraz diagram dyspersji Re(k)/Im(k) na żywo poniżej, z zaznaczoną częstotliwością przejścia ωc, oddzielającą reżim dyfuzyjny od reżimu fali propagującej.
🎮 Jak korzystać
Dostosuj Moduł ścinania G, Lepkość η, Częstotliwość wymuszenia ω i Amplitudę A suwakami. Iloczyn ωτ pokazany na płótnie jest parametrem kontrolnym: przesuń go daleko poniżej 1 (niskie η, niskie ω) dla reżimu dyfuzyjnego, w którym fala umiera w ułamku obszaru, i daleko powyżej 1 dla reżimu sprężystego, w którym fala dzwoni przez cały ośrodek. Naciśnij P/R, by wstrzymać/zresetować.
💡 Czy wiesz, że…
„Głupi glut” to klasyczna sztuczka imprezowa demonstrująca tę fizykę: ciągnięty powoli płynie jak ciecz (poniżej ωc), ale uderzony gwałtownie pęka jak kruche ciało sprężyste (powyżej ωc).
Najczęściej zadawane pytania
Dlaczego odczyt Reżimu przełącza się między Sprężystym a Lepkim?
Kod porównuje częstotliwość wymuszenia omega z częstotliwością przejścia wc = G/eta, obliczaną w każdej klatce; gdy omega*tau przekracza 1, oznacza ośrodek jako Sprężysty (fala), w przeciwnym razie jako Lepki (dyfuzyjny), zgodnie z teorią modelu Maxwella.
Co pokazuje pomarańczowa linia nałożona na wykres przemieszczenia?
Wykreśla ona wewnętrzne pole naprężenia, które równanie konstytutywne Maxwella ewoluuje niezależnie od przemieszczenia, ukazując, jak naprężenie i odkształcenie mogą być przesunięte w fazie w ośrodku lepkosprężystym. Jest przechowywane na komórkach połówkowych przesuniętej siatki, w połowie drogi między węzłami prędkości.
Jak obliczana jest krzywa dyspersji w dolnym panelu?
Dla każdej częstotliwości na osi x kod oblicza zespoloną liczbę falową z k^2 = rho*omega^2*(1 - i/(omega*tau))/G, wykreślając Re(k) na niebiesko (propagacja) i Im(k) na czerwono (tłumienie) na wspólnej skali.
Dlaczego zwiększenie lepkości zmienia częstotliwość przejścia?
Czas relaksacji tau = eta/G rośnie wraz z lepkością, a ponieważ wc = 1/tau, wyższa wartość eta obniża częstotliwość przejścia, sprawiając, że ośrodek zachowuje się sprężyście w szerszym zakresie częstotliwości — a ponieważ długość tłumienia 1/Im(k) dąży do 2*tau*sqrt(G/rho), fala może przebyć większą odległość, zanim zaniknie.
Dlaczego siatka jest przesunięta?
Prędkość znajduje się na węzłach, a naprężenie na węzłach połówkowych, więc każda pochodna przestrzenna jest wycentrowana dokładnie nad jedną komórką. Schemat skolokowany, różniczkujący na dwóch komórkach, pozostawia węzły parzyste i nieparzyste nierozłączone i dopuszcza tryb piłokształtny (szachownicowy), który rośnie bez ograniczeń wraz ze zmniejszaniem się rozmiaru komórki. Przesunięcie eliminuje ten tryb, co pozwala na zmniejszenie komórek na tyle, by rozdzielić falę.
Skąd wiadomo, że solwer jest stabilny?
Jawny schemat leapfrog jest stabilny, dopóki liczba Couranta c*dt/dx pozostaje poniżej 1, gdzie c = sqrt(G/rho) to nierelaksowana prędkość fali; kod przelicza podkrok na podstawie bieżącego G w każdej klatce i utrzymuje liczbę Couranta na poziomie 0,4. Człon relaksacji naprężenia jest całkowany dokładnie, a nie jawnym krokiem Eulera, więc bardzo krótki czas relaksacji również nie może zdestabilizować schematu.
Co dzieje się na prawej krawędzi siatki przemieszczenia?
Stosowana jest tam charakterystyczna (jednokierunkowa) warunek pochłaniający: fala biegnąca w prawo w tym ośrodku spełnia sigma = -rho*c*v, więc ostatni węzeł prędkości jest ustalany na podstawie ostatniej wartości naprężenia. Fala opuszcza obszar, zamiast odbijać się z powrotem.