Hot answers tagged

4

The explanation in the book does not use von Neumann analysis at all but the absolute stability regions and the eigenvalues of the discrete Laplacian operator. For the result you specifically mentioned we use the fact that the maximum eigenvalues is $$ \lambda_m \approx -\frac{4}{h^2} $$ from the expression given. We then want this eigenvalue to lie inside ...


3

I believe Von Neumann's stability analysis would give you the answer here. Consider the heat transfer equation: $$\frac{\partial \mathcal{T}}{\partial t} = \alpha \frac{\partial^{2} \mathcal{T}}{\partial x^{2}}$$ By using Forward Euler time integration and central difference in space discretization: $$\mathcal{T}^{t+\Delta t}_{x} = \mathcal{T}^{t}_{x} + \...


2

Your nonlinear equation is a viscous Burgers' equation and could be linearized easily by using Cole-Hopf transformation. If I rewrite your equation as: $$u_{t} = 2 c \Big (\frac{\alpha}{2 c} u_{x} - \frac{u^{2}}{2}\Big)_{x}$$ And take: $$\mu = \frac{\alpha}{2c}$$ $$u = -2 \mu \frac{\phi_{x}}{\phi}$$ I would have: $$u_{x} = \frac{2 \mu}{\phi^{2}} (\phi_{x}^{2}...


2

You could refine your discretization and then compare the logarithm of the error ($\log |e|$) with the logarithm of the size ($\log h$) of your elements. Using a linear regression you could obtain an approximation of the order of convergence. Keep in mind that this order of convergence is asymptotic when $h \rightarrow 0$, so, for "large" $h$ you ...


Only top voted, non community-wiki answers of a minimum length are eligible