eqn = Q[z] D[u[x, y, z, t], t] + Q[z] U[z] D[u[x, y, z, t], x] +
Q[z] w[x, y, z, t] D[U[z], z] == -D[p[x, y, z, t], x] +
M[z]/3 (D[u[x, y, z, t], {x, 2}] +
D[u[x, y, z, t], {x, 1}, {y, 1}] +
D[u[x, y, z, t], {x, 1}, {z, 1}]) +
D[M[z], z] (D[w[x, y, z, t], x] + D[u[x, y, z, t], z])
方程我是故意简化了的,本来是一个非常复杂的方程组,为了方便,我取其中一个方程,并且省略了很多项。为的是专注于方法,而非某个解决具体问题。
Q[z] w[x, y, z, t] D[U[z], z] == -D[p[x, y, z, t], x] +
M[z]/3 (D[u[x, y, z, t], {x, 2}] +
D[u[x, y, z, t], {x, 1}, {y, 1}] +
D[u[x, y, z, t], {x, 1}, {z, 1}]) +
D[M[z], z] (D[w[x, y, z, t], x] + D[u[x, y, z, t], z])
方程我是故意简化了的,本来是一个非常复杂的方程组,为了方便,我取其中一个方程,并且省略了很多项。为的是专注于方法,而非某个解决具体问题。