CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Stability Analysis

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

LinkBack Thread Tools Display Modes
Old   December 15, 1998, 07:47
Default Stability Analysis
Posts: n/a
While trying to solve an equation of the form A * y'' - B * Y + c = 0 using finite elements in a 3d hemishpere problem I find the surface contours plot becomes unstable for higher values of Coeff B. Is there anyway I can control this parameter and get decent results?.
  Reply With Quote

Old   December 18, 1998, 13:12
Default Re: Stability Analysis
John C. Chien
Posts: n/a
The equation can be written as Y'' =(B/A) * Y - (C/A). If you use finite-difference for Y'' , you get (Y(i+1) -2Y(i) +Y(i-1)) / ((x(i+1)-x(i-1))/2)**2. Also Y can be written as Y(i). If we let dx2=((x(i+1)-x(i-1))/2)**2, then the final equation for Y(i) is Y(i)=(Y(i+1)+Y(i-1))/(2 + (B/A)*dx2) + (dx2*(C/A))/(2 + (B/A)*dx2 ). When the mesh size is small, dx2 is small. Y(i) is approximately equal to (Y(i+1)+Y(i-1))/2. So, as long as Y(i+1) and Y(i-1) are finite, the solution Y(i) ( the value between x(i+1) and x(i-1)) will be finite and bounded by Y(i+1) and Y(i-1). In this case is the average value of Y(i+1) and Y(i-1). For finite value of mesh size, dx2 will be finite ( not a very small value ). And for large value of B, the denominator (2+(B/A)*dx2) can become a large number. For example, for finite Y(i+1) and Y(i-1), the solution Y(i) can become very small or approach zero. In this case, Y(i-1) > Y(i) < Y(i+1). So a V-shaped profile is created. This result is artificial because the real solution should always close to the average value of the Y(i+1) and Y(i-1) by making mesh size very small, and thus very small dx2. So you should try to keep (B/A)*dx2 << 2, by making mesh size small. Or dx2 << (2/ (B/A)). I think, the finite element works in a similar way( with additional assumptions).
  Reply With Quote

Old   January 22, 1999, 07:27
Default Re: Stability Analysis
Posts: n/a
As this problem is in 3-d with varying material shape and thickness I could only reach an optimum mesh size. I would like to know if I can use some alogrithm to get over the wiggles created in my solution as I unable to further refine my grid ? natteri sudharsan
  Reply With Quote


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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Finite difference stability analysis Frank Main CFD Forum 2 August 31, 2006 14:30
Stability Analysis skg Main CFD Forum 3 June 15, 2006 08:41
numerical stability analysis Wen Long Main CFD Forum 1 May 22, 2002 00:14
Stability Analysis W.M.Leung Main CFD Forum 0 January 7, 2000 13:39
Stability analysis on a non-uniform grid PRAVEEN CHANDRASHEKAR Main CFD Forum 7 July 2, 1999 01:46

All times are GMT -4. The time now is 02:42.