16.920J/SMA 5212 Numerical Methods for PDEs E△E Complementary Particular(steady-state) (transient )solution solution where (cy)=[ The stability analysis of the space discretization, keeping time continuous, is based on the eigenvalue structure of a The exact solution of the system ofequations is determined by the eigenvalues and eigenvectors ofa STABILITY ANALYSIS Coupled ODEs to Uncoupled ODEs We can think of the solution to the semi-discretized problem u=Elce-EA-E-lb as a superposition of eigenmodes of the matrix operator A Each mode j contributes a( transient )time behaviour of the form d to the time-dependent part of the solution Since the transient solution must decay with time Real()≤0 for all j This is the criterion for stability of the space discretization(of a parabolic PDE) keeping time continuous Slide 5
16.920J/SMA 5212 Numerical Methods for PDEs 11 Evaluating, ( ) t 1 1 u EU E ce E E b λ − − = = − Λ ✂ ✂ ✁✁✂ ✂ ( ) 1 2 1 w 1 2 1 here j N T t t t t t j N ce c e c e c e c e λ λ λ λ λ − − ✄ ☎ = ✆ ✝ ✞✁✞✁✞✟✞✠ The stability analysis of the space discretization, keeping time continuous, is based on the eigenvalue structure of A. The exact solution of the system of equations is determined by the eigenvalues and eigenvectors of A. Slide 14 STABILITY ANALYSIS Coupled ODEs to Uncoupled ODEs We can think of the solution to the semi-discretized problem as a superposition of eigenmodes of the matrix operator A. Each mode contributes a (transient) time behaviour of the form to the time-dependent part of the solution. j t j e λ Since the transient solution must decay with time, Real ( ) 0 λ j ≤ for all j This is the criterion for stability of the space discretization (of a parabolic PDE) keeping time continuous. Slide 15 Complementary (transient) solution Particular (steady-state) solution ( ) t 1 1 u E ce E E b λ − − = − Λ ☛ ✡✁✡✁✡✟✡☛ ☛
16.920J/SMA 5212 Numerical Methods for PDEs STABILITY ANALYSIS Use of Modal (scalar) Equation It may be noted that since the solution iyTs expressed as a contribution from all the modes of the initial solution which have propagated or(and) diffused with the eigenvalue A. and a contribution from the source term b. all the properties of the time integration (and their stability properties) can be analysed separately for each mode with he scalar equation dU dt Slide 16 STABILITY ANALYSIS Use of Modal ( scalar) Equation The spatial operator A is replaced by an eigenvalue 1, and the above modal equation will serve as the basic equation of the stability of a time-integ (yet to be introduced) as a function of the eigenvalues a of the space-discretization operators This analysis provides a general technique for the determination of time integration methods which lead to stable algorithms for a given space discretization Slide 7 EXAMPLE 1 Continuous Time Operator Consider a set of coupled ODEs (2 equations only) a1l1+a12l2 dt 141 +all du lide 18
16.920J/SMA 5212 Numerical Methods for PDEs 12 STABILITY ANALYSIS Use of Modal (Scalar) Equation It may be noted that since the solution is expressed as a contribution from all the modes of the initial solution, which have propagated or (and) diffused with the eigenvalue , and a contribution fr j u λ om the source term , all the properties of the time integration (and their stability properties) can be analysed separately for each mode with the scalar equation j b Slide 16 STABILITY ANALYSIS Use of Modal (Scalar) Equation The spatial operator A is replaced by an eigenvalue λ, and the above modal equation will serve as the basic equation for analysis of the stability of a time-integration scheme (yet to be introduced) as a function of the eigenvalues λ of the space-discretization operators. This analysis provides a general technique for the determination of time integration methods which lead to stable algorithms for a given space discretization. Slide 17 EXAMPLE 1 Continuous Time Operator Consider a set of coupled ODEs (2 equations only): 1 11 1 12 2 2 21 1 22 2 du a u a u dt du a u a u dt = + = + 1 11 12 2 21 22 Let , u a a du u A Au u a a dt ✁✄✂ ✁ ✂ = = ☎ = ✆ ✝ ✆ ✝ ✞✄✟ ✞ ✟ ✠ ✠ ✠ Slide 18 j dU U F dt λ ✡ ☛ = + ☞ ✌ ✍ ✎