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

Runge Kutta 4th Order Source Code

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

LinkBack Thread Tools Display Modes
Old   February 4, 2002, 05:50
Default Runge Kutta 4th Order Source Code
Posts: n/a
Hi all. I need the source code for solving ODE using 4th order R.K method. I am not good at programming, but my research needs me to simulate fluid flow inside a micro channel. Is there anyone who can help me to develop a small programme? I am quite desperate now. Below is the rough equation (ca, w, b, and sin(theta) are constants). I need the graph for s-x. Any volunteers?

ds/dx = Ca/w [s**3/2(1/s**3 + 1/(1-s)**3)-Bsin(theta)]
  Reply With Quote

Old   February 4, 2002, 06:35
Default Re: Runge Kutta 4th Order Source Code
Christian Holm
Posts: n/a

Have you looked in the book : "Numerical Recipes in ***", where *** can be Pascal, C, C++ and Fortran (maybe there are more programming languages available). The book is from Cambridge University Press and is written by Press et al.

It is VERY good. I have it for Pascal (unfortunately they claim that the market for Delphi is too small), and I have used it a lot.

Cheers Christian

  Reply With Quote

Old   February 4, 2002, 23:46
Default Re: Runge Kutta 4th Order Source Code
Junseok Kim
Posts: n/a
I can help you, would you send me an email with more informations, related papers or the values of constants, initial contitions, etc.

And what kind of formats of programming do you want? Matlab, C, C++, f77, f90, Mathematica, etc. Actually, Matlab has Runge Kutta 4th Order routinue.

  Reply With Quote

Old   October 25, 2012, 20:53
New Member
Praveen Gambhir
Join Date: Oct 2012
Posts: 1
Rep Power: 0
gvpraveen16 is on a distinguished road
Hey guys!! I have a question regarding solution of a 3rd order ODE using 4th order Runge Kutta technique.
The function is

f= 0.5(3y-y^3)+k(3y^3-2y-y^7)

and should plot graphs between y vs f and y vs df/dy for different values of the constant k. The initial and final values of y are -1 and 1 and delta(y)=0.005. Please do help as all examples for this method show two variables and time dependency whereas the above function is time independent and has only one variable. Awaiting your reply
gvpraveen16 is offline   Reply With Quote

Old   October 26, 2012, 03:15
Senior Member
Join Date: Mar 2010
Posts: 557
Rep Power: 11
cfdnewbie is on a distinguished road
So what exactly is your question? You can program it up and plot it, if you have matlab, that would be a good choice!
cfdnewbie is offline   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
Fourth Order Runge Kutta time integration lr103476 OpenFOAM Running, Solving & CFD 22 August 21, 2015 11:09
4th order Runge-Kutte & uncoupled method Navier Stokes equations misabel Main CFD Forum 0 February 10, 2010 07:06
Understanding k-omega SST model source code tmhonka OpenFOAM Programming & Development 1 September 8, 2009 07:33
Runge Kutta Optimization vasanth Main CFD Forum 6 December 2, 2005 14:07
public CFD Code development Heinz Wilkening Main CFD Forum 38 March 5, 1999 12:44

All times are GMT -4. The time now is 13:14.