〰️ Поширення Хвиль у В'язкопружному Середовищі

Хвилі у в'язкопружному середовищі Максвела. Нижче ωc = G/η вони дифундують; вище — поширюються. Дивіться, як змінюються криві дисперсії та загасання.

ФізикаІнтерактивний
Верхня панель: поле переміщень · Нижня: дисперсія · P пауза · R скидання

Як це працює

Синусоїдальне переміщення збуджується на лівій межі. Симуляція поширює 1D зсувну хвилю через в'язкопружне середовище Максвела методом скінченних різниць у часовій області. Модель Максвела дає частотно-залежний комплексний модуль зсуву G*(ω) = G·iωτ/(1+iωτ), що визначає комплексний хвильовий вектор k(ω).

При низьких частотах (ω << ωc = 1/τ) уявна частина k переважає — збурення загасають експоненційно (дифузний режим). При високих частотах (ω >> ωc) дійсна частина k переважає — хвилі когерентно поширюються зі швидкістю c = sqrt(G/ρ) (пружний режим). Нижня панель показує живу дисперсійну криву.

τ = η / G [час релаксації]
ω_c = 1/τ = G/η [частота переходу]
G*(ω) = G · iωτ / (1 + iωτ)
k = ω √(ρ / G*(ω))

Часті запитання

Що таке в'язкопружний матеріал?

В'язкопружний матеріал виявляє як в'язку, так і пружну поведінку залежно від часового масштабу деформації. Приклади: розплави полімерів, біологічні гелі, пластилін.

Що таке модель Максвела для в'язкопружності?

Модель Максвела представляє в'язкопружність як пружину (модуль G) і демпфер (в'язкість η) послідовно. Напруження релаксується з часом релаксації τ = η/G.

Що таке час релаксації у в'язкопружності?

Час релаксації τ = η/G — час, за який напруження спадає до 1/e від початкового значення. Він розмежовує пружний і в'язкий режими.

Чому хвилі дифундують при низьких частотах?

При частотах нижче ωc = G/η матеріал встигає релаксувати і поводиться як в'язка рідина. Збурення загасають дифузно.

Чому хвилі поширюються при високих частотах?

При частотах вище ωc матеріал не встигає релаксувати і поводиться пружно, підтримуючи зсувні хвилі зі швидкістю c = sqrt(G/ρ).

Яке дисперсійне співвідношення для в'язкопружних хвиль?

Для рідини Максвела: k² = ρω²/G·(1 + i/(ωτ)). Дійсна частина k — фазова швидкість, уявна — просторове загасання.

Що таке загасання при поширенні хвиль?

Загасання — спад амплітуди хвилі з відстанню. У в'язкопружних середовищах виникає через в'язкісну дисипацію. Коефіцієнт загасання α = Im(k) зростає з в'язкістю.

Де важливі ефекти в'язкопружних хвиль?

В'язкопружні хвилі важливі для сейсмічного загасання в осадових породах, ультразвуку в тканинах, поглинання ударів у полімерних пінах та демпфування коливань у композитах.

Що таке число Дебори?

Число Дебори De = τ/tспост порівнює час релаксації з часом спостереження. De >> 1 — пружна поведінка; De << 1 — в'язка.

Чим модель Кельвіна-Войґта відрізняється від моделі Максвела?

Модель Кельвіна-Войґта розміщує пружину і демпфер паралельно, даючи повзучість, але не релаксацію напружень. Модель Максвела — послідовно, даючи релаксацію, але необмежену повзучість.

Про цю симуляцію

Скінченнорізницевий розв'язувач крокує напругою та швидкістю на 1D сітці з 300 точок, використовуючи визначальний закон Максвелла dσ/dt = G(dv/dx) − σ/τ, з приводом синусоїдального зміщення на лівому краю. Оскільки час релаксації τ = η/G задає частоту переходу ωc = G/η, те саме середовище може виглядати як текуча рідина або дзвінке пружне тіло — залежно виключно від того, наскільки швидко ним трясти.

🔬 Що показано

Керована хвиля зміщення у верхній панелі та живий діаграма дисперсії Re(k)/Im(k) знизу, з позначеною частотою переходу ωc, що розділяє дифузійний режим і режим поширення хвиль.

🎮 Як користуватись

Регулюйте Модуль зсуву G, В'язкість η, Частоту приводу ω та Амплітуду A повзунками; спостерігайте, як показник Regime перемикається між Elastic і Viscous, або натискайте P/R для паузи/скидання.

💡 Чи знали ви?

Play-Doh і подібна «розумна замазка» — класичний фокус для цієї фізики: тягніть повільно, і вона тече як рідина (нижче ωc), а різко смикніть — і вона трісне як крихке пружне тіло (вище ωc).

Часті запитання

Чому показник Regime перемикається між Elastic і Viscous?

Код порівнює частоту приводу omega з частотою переходу wc = G/eta, обчисленою щокадру; коли omega перевищує wc, він позначає середовище як Elastic (хвиля), інакше — Viscous (дифузія), відповідно до теорії моделі Максвелла.

Що показує помаранчева лінія поверх графіка зміщення?

Вона малює внутрішнє поле напруги sigma[i], яке рівняння Максвелла обчислює окремо від зміщення, показуючи, як напруга й деформація можуть бути не в фазі у в'язкопружному середовищі.

Як обчислюється крива дисперсії у нижній панелі?

Для кожної частоти на осі x код обчислює комплексне хвильове число з k² = ρω²(1+i/(ωτ))/G, малюючи Re(k) синім (поширення) та Im(k) червоним (затухання).

Чому збільшення в'язкості змінює частоту переходу?

Час релаксації tau = eta/G зростає з в'язкістю, а оскільки wc = 1/tau, вище значення eta знижує частоту переходу, змушуючи середовище поводитись в'язко у ширшому діапазоні частот.

Що відбувається на правому краю сітки зміщення?

Код копіює передостаннє значення сітки в останнє на кожному кроці, реалізуючи просту поглинальну межу, щоб хвилі виходили, не відбиваючись назад у домен.