How Lnorms are used to study stability of a numerical scheme?
Greetings,
Can someone tell me in what way L1, L2 or infinity norm is used to study stability analysis of a numerical scheme? I am trying to develop a code to solve compressibleNavierStokes equation using compact scheme for spatial discretization and rungekutta scheme for temporal discretization. What other methods can be used for stability analysis of this scheme? Thank you. 
Quote:

Quote:
Thanks for your reply! Now, in addition to studying the norms for vanishing step sizes, I would like to know how can the norms be used in determining stability of a numerical scheme. My next question is about that. I applied the scheme mentioned in my previous post to one dimensional linear wave equation (u_t + a * u_xx = 0) with Gaussian function as initial condition. I used cfl of 0.5. I have attached a plot in this post which shows the L1, L2 and infinity norm of the error. The error is defined as the difference between exact solution and analytical solution. The plots show that the Lnorms are uniformly increasing with number of time iterations. Now, from this plot, what can I say about the stability of this scheme? This is the link to the plot. http://dl.dropbox.com/u/56389861/norms%20of%20error.png 
a) The linear wave equation is u_t + a*u_x = 0, is a first order equation, not a second order.
b) the stability of a numerical scheme can be studied by the von Neumann analysis, the magnitude of the amplification factor must be less the unity. c) The Lax theorem implies that a consistent and stable linear scheme must converge to the exact solution Therefore, the stability in your case can be studied analitically, you have no reason to do a numerical case. You could use a norm on the error to verify the error slope and the accuracy order. 
Quote:

Thank you for your replies!!!
Sorry for the error. I did mean the equation u_t + a*u_x = 0 and the plot is for this equation. I used this equation just as a simple example to ask my question about the use of Lnorms in determining the stability of a numerical scheme. Now, vonNeumann stability analysis cannot be used for compressibleNavierStokes equations. So, can the concept of linear growth of Lnorms of error be used as a stability criterion? What other methods are used to study the stability of a numerical scheme meant for compressibleNavierStokes equations? Kind regards. 
Dear Ravindra,
i suggest you to read: LeVeque: Finite Volume Methods for Hyperbolic Problems  Chapter 8 Hirsch: Numerical Computation of Internal and External Flows  Chapter 7, 8, 9 where the concepts of norms, stability and convergence are fully clarified. However, as you said, classical Von Neumann analysis (which, by the way, is stability in L2 norm) is not suitable for non linear problems. In that case (but this really is not my field) i think you have to move to more general concepts like Total Variation Bounding where, still, some specific norm is applied to some specific quantity. However, i don't think that there are general stability results concerning systems of nonlinear equations like the compressible NSE. 
However, as said by cfdnewbie, as long as the error norms grow linearly (and not exponentially) you can conclude that there is no instability (it is just the buildup of error), no matter what system you are considering. The converse, of course, could not be true.

All times are GMT 4. The time now is 06:51. 