Processing math: 100%
Extension to quadratic damping
mu″+β|u′|u′+s(u)=F(t).
Centered scheme: |u′|u′ gives rise to a nonlinearity.
Linearization trick: use a geometric mean,
[|u′|u′]n≈|[u′]n−12|[u′]n+12.
Scheme:
[mDtDtu]n+β|[Dtu]n−12|[Dtu]n+12+s(un)=Fn.