|
[Sponsors] |
February 3, 2012, 09:02 |
Most common ODE solver
|
#1 |
New Member
Nereus
Join Date: Sep 2011
Posts: 18
Rep Power: 14 |
It's my understanding that the method of lines is quite commonly used in CFD. What is the most common ODE solver used with this approach? For example, LSODA, VODE, or simply 4th order Runge-Kutta?
|
|
February 3, 2012, 09:12 |
|
#2 |
Senior Member
cfdnewbie
Join Date: Mar 2010
Posts: 557
Rep Power: 20 |
||
February 3, 2012, 14:11 |
|
#3 |
Member
Join Date: Nov 2011
Location: Czech Republic
Posts: 97
Rep Power: 14 |
It also depends on solved equations. If a problem is stiff, implicit numerical methods are more suitable for that task. For non-stiff problems Dorman-Prince method is the first shot integrator in many numerical solvers.
|
|
February 3, 2012, 17:56 |
|
#4 |
Senior Member
Join Date: Aug 2011
Posts: 272
Rep Power: 15 |
Yes I agree with cfdnewbie. RK4 is known to be the most accurate method (if you exclude RK5, RK6 etc....)
But if you plan to use it to solve Navier-Stokes equations in an explicit way, RK4 is quite expensive in CPU time and storage. A special RK3 method has been specifically designed to lower that cost. |
|
February 3, 2012, 18:05 |
|
#5 |
Senior Member
Join Date: Aug 2011
Posts: 272
Rep Power: 15 |
Dormand-Prince method is linked to RK method familly and seems to be more expensive than RK4. check this link for a gpl subroutine using this method.
http://www.unige.ch/~hairer/prog/nonstiff/dopri5.f RK4 can be coded with 20 lines !!! |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
ODE solver coefficients | mrangitschdowcom | OpenFOAM Running, Solving & CFD | 2 | August 18, 2013 22:16 |
Working directory via command line | Luiz | CFX | 4 | March 6, 2011 20:02 |
Problem implementing CVODE ODE solver | markusrehm | OpenFOAM | 20 | October 13, 2010 17:02 |
Creating New Solver: For particle-laden compressible jets | sankarv | OpenFOAM | 0 | April 4, 2010 18:06 |
why the solver reject it? Anyone with experience? | bearcat | CFX | 6 | April 28, 2008 14:08 |