Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- eqns = {D[c[x, t], t] == (2 c[x, t] + 1)*D[v[x, t], x] -
- D[c[x, t], x]*v[x, t] - D[v[x, t], x] c[x, t],
- 2 D[v[x, t], t] == (c[x, t] + 1)*(8 c[x, t]*D[c[c, t], x] + 4 D[c[x, t], x])*
- D[v[x, t], x] + (c[x, t] + 1)*(4 c[x, t]^2 + 4 c[x, t] + 1)*D[v[x, t], x, x] +
- 4 (c[x, t] + 1)*D[c[x, t], x]*mu[x, t] + (2 c[x, t] + 1)*(c[x, t] + 1)*
- D[mu[x, t], x] + (c[x, t] + 1)*D[v[x, t], x, x] 11/3 - 2 v[x, t]*D[v[x, t], x],
- mu[x, t] - log[c[x, t]] - 1 - 4 c[x, t] + 3 D[c[x, t], x, x] == 0};
- bcs = {v[0, t] == 0, v[1, t] = 0, Derivative[1, 0][c][0, t] == 0, Derivative[1, 0][c][1, t] == 0};
- ics = {c[x, 0] == 0.3, v[x, 0] == 0, mu[x, 0] == 0.6960};
- sol = NDSolveValue[{eqns, bcs, ics}, {c, v,mu}, {x, 0, 1}, {t, 0, 3},
- Method -> {"IndexReduction" -> Automatic, "EquationSimplification" -> "Residual",
- "PDEDiscretization" -> {"MethodOfLines",
- "SpatialDiscretization" -> {"TensorProductGrid", "MinPoints" -> 200, "MaxPoints" -> 200}}}]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement