CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Solving ODEs in OpenFOAM

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

LinkBack Thread Tools Display Modes
Old   April 19, 2010, 07:50
Default Solving ODEs in OpenFOAM
Wolfgang W.
Join Date: Nov 2009
Location: Switzerland
Posts: 57
Rep Power: 9
WiWo is on a distinguished road
Hi Foamers,

I'm having a quite unexpected difficulty at the moment. For the equations in my current solver I need to determine a single scalar parameter from a non-linear 2nd order ordinary differential equation (initial value problem) at every time step. This parameter does not belong to any field - it's just a number (on per time step) I need to plug into the solver.
Now it doesn't seem to be straight forward to solve an ODE (in my case I want to get the radius of a micro bubble from the Rayleigh-Plesset equation) for just this one parameter in OpenFOAM. Does OpenFOAM offer the capabilities of solving such ODEs?

I'd appreciate any help or advice.

WiWo is offline   Reply With Quote

Old   April 20, 2010, 08:49
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,803
Rep Power: 24
hjasak will become famous soon enough
Yes, absolutely. We have the ODE library and ODE solvers. You can write your ODE as a derived class and use an ODE solver with it.

There is a relatively simple example in 1.5-dev under silvers/basic.

Good hunting,

Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote

Old   April 27, 2010, 07:54
Wolfgang W.
Join Date: Nov 2009
Location: Switzerland
Posts: 57
Rep Power: 9
WiWo is on a distinguished road
Dear Prof. Jasak,

The hunt was very successful - thank you very much for referencing me to the sixDOFSolver!

I was able to adapt the ODE solver and its dedicated classes to fit my needs. It was then pretty straight forward to include this package into my overall solver - all done, verified and working :-)

I recognized that there is quite a significant difference between the ODESolver class of OF-1.5-dev and OF-1.6. My first attempt was actually to work with the ODETest case (applications/test) in OF-1.5-dev, but it was evidently designed to work with the 1.6(.x)-version of ODESolver. What is the reason for this discrepancy?

Thank you for your help!

Best regards,
WiWo is offline   Reply With Quote

Old   August 4, 2013, 02:59
Default about the solution of rayleigh plesset equation
Anil Kunwar
Join Date: Jun 2013
Posts: 55
Rep Power: 5
Annier is an unknown quantity at this point
dear Wolfgang,
would you please guide me with the steps of solving the rayleigh plesset equation in OpenFOAM?
thanking you in advance.
Annier 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16
Unknown error sivakumar OpenFOAM Pre-Processing 9 September 9, 2008 12:53
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

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