Problem with Backward Euler Results

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

 December 22, 2011, 11:51 Problem with Backward Euler Results #1 Senior Member   Join Date: Jun 2010 Posts: 111 Rep Power: 8 Hi, I have tried to solve a system of equations in the form: dy1/dt = y1 +dt*f(y1,y2) dy2/dt = y2 +dt*f(y1,y2) using different schemes such as forward Euler and backward (implicit) Euler and Runge Kutta order 4. My results seem to be OK for the forward Euler scheme and Runge Kutta above a certain number of time steps (A). But for implicit Euler using A time steps, I get an error after a certain number of time steps (say A-100). I was wondering if anyone had any ideas of what the problem might be. I can't think of anything at all. Thanks!

 December 23, 2011, 17:44 #2 Member   private Join Date: Mar 2009 Posts: 74 Rep Power: 9 [QUOTE=Hooman;336756]Hi, I have tried to solve a system of equations in the form: dy1/dt = y1 +dt*f(y1,y2) dy2/dt = y2 +dt*f(y1,y2) Do you perhaps mean f( ... ) for the first equation and g( ... ) for the second equation?

 December 25, 2011, 09:37 #3 Senior Member   Join Date: Jun 2010 Posts: 111 Rep Power: 8 yes, that's what I meant. Sorry!

 December 29, 2011, 12:53 #4 Senior Member   Join Date: Jun 2010 Posts: 111 Rep Power: 8 Any ideas?

 January 1, 2012, 11:59 #5 Senior Member   Join Date: Jun 2010 Posts: 111 Rep Power: 8 I just realized that if I use really small time steps (much smaller than that of explicit euler) my reuslts become OK. I am now thinking it could be because I used forward euler to work out the initial RHS's at the new time step, to calculate the y1 and y2 at the new time steps.

 January 1, 2012, 12:02 #6 Senior Member   Join Date: Jun 2010 Posts: 111 Rep Power: 8 ooh I just realize that I have written the equations wrong here. They should be: dy1/dt = f(y1,y2) dy2/dt = g(y1,y2)

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ismael.s CFX 13 February 27, 2012 09:00 nasir CFX 2 November 4, 2011 22:10 idrama OpenFOAM 3 March 4, 2010 06:07 Rams CFX 2 June 13, 2006 03:30 Question Man. Main CFD Forum 2 March 29, 2004 13:19

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