Processing math: 100%
The general model formulated as first-order ODEs
m
u
″
+
β
|
u
′
|
u
′
+
s
(
u
)
=
F
(
t
)
.
Rewritten as first-order system:
u
′
=
v
,
v
′
=
1
m
(
F
(
t
)
−
β
|
v
|
v
−
s
(
u
)
)
.
To solution methods:
Forward-backward scheme
Centered scheme on a staggered mesh
