Changes between Version 13 and Version 14 of u/erica/2D_Godunov
- Timestamp:
- 08/07/13 15:08:03 (11 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/2D_Godunov
v13 v14 11 11 Here are some plots of the initial condition. You will see a circle of high density/pressure in the center of the mesh, with zero velocity in both x and y. 12 12 13 Density :13 Density (pressure looks similar): 14 14 15 15 [[Image(rhoInit3d.png, 25%)]] 16 16 17 Velocity:17 X-velocity (velocity in y is also zero): 18 18 19 19 [[Image(uInit.png, 25%)]] … … 27 27 [[latex($\vec{u}^n = \{ \vec{u}^n_{i,j} \} ~\forall ~i,~\forall~j$)]] 28 28 29 where u is the vector of primitive variables. To solve the Euler equations for this given initial condition, we splitthe 2D Euler equations,29 where u is the vector of primitive variables. To solve the Euler equations for this given initial condition, we "split" the 2D Euler equations, 30 30 31 31 [[latex($\frac{\partial }{\partial t} \begin{pmatrix} \rho \\ \rho u \\ E \\ \rho v \end{pmatrix} + \frac{\partial }{\partial x} \begin{pmatrix} \rho u \\ \rho u^2 + p \\ u(E+p) \\ \rho u v \end{pmatrix} + \frac{\partial }{\partial y} \begin{pmatrix} \rho v \\ \rho v^2 + p \\ v(E+p) \\ \rho u v \end{pmatrix} = 0$)]]