CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Stability and Amplification Factor

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   January 28, 2010, 18:05
Default Stability and Amplification Factor
  #1
hyd
New Member
 
Anonymous
Join Date: Jan 2010
Posts: 1
Rep Power: 0
hyd is on a distinguished road
I have a set of equations where the equations can be written in the following format; dU/dt+A*dU/dx=0 (these are partial differential)
If I discretize this with explicit FTCS
(Ui^n+1)-(Ui^n)=-(A*Delta_t/(2*Delta_x))*((Ui+1^n)-(Ui-1^n)) and try to find out the stability limit, I employ the Von Neumann stability analysis.

For single equations all there is that needs to be done is find G and then find the stability limit from |G|^2<1 (smaller than or equal to), but in this case there is a set of equations which complicates the matter. Let's say
A is a 2x2 matrix
A=[ a b
c a]
Now what needs to be done is finding the eigenvalues and what I end up with is; (I=imaginary)
eigen1=(1-(Delta_t/Delta_x)*((b*c)^0.5)*sintheta)-((Delta_t/Delta_x)*I*a*sintheta)
eigen2=(1+(Delta_t/Delta_x)*((b*c)^0.5)*sintheta)-((Delta_t/Delta_x)*I*a*sintheta)

The stability limit is |eigen*conjugateeigen|<1 (smaller than or equal to 1)
eigen1*conjugateeigen1=(1-(Delta_t/Delta_x)*((b*c)^0.5)*sintheta)^2+((Delta_t/Delta_x)*a*sintheta)^2

The answer is unconditionally unstable for sintheta=-1, but what about sintheta=1. The result should still be unconditionally unstable but it doesn't seem like it.
With sintheta=1
eigen1*conjugateeigen1=(1-(Delta_t/Delta_x)*((b*c)^0.5))^2+((Delta_t/Delta_x)*a)^2 where this can be smaller than 1.
Where am I going wrong?
Thank you.
hyd is offline   Reply With Quote

Reply

Tags
amplification factor, stability

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On



All times are GMT -4. The time now is 16:22.