3 недели назад
Решил визуализировать несколько основных сейсмических волн. Все волны ниже — решения одного уравнения эластики (уравнение движения упругой среды): ρ ∂²u/∂t² = (λ + 2μ) ∇(∇·u) − μ ∇×∇×u. где ρ — плотность, λ и μ — параметры Ламе, ν — коэффициент Пуассона. Подстановка u = ∇φ + ∇×ψ (разложение Гельмгольца) разделяет его на два волновых уравнения: ∂²φ/∂t² = α² ∇²φ (P-волна), ∂²ψ/∂t² = β² ∇²ψ (S-волна). Солвер решает плоскую деформацию (компоненты ux, uz, поляризация P–SV), поэтому всё ниже, кроме волны Лява (но ее я визуализирую в 3D, где есть SH), в нём считается...
