CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

rhoCentralFoam Runge Kutta

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

Like Tree20Likes

LinkBack Thread Tools Search this Thread Display Modes
Old   May 12, 2015, 15:58
New Member
Join Date: Mar 2012
Posts: 1
Rep Power: 0
philem is on a distinguished road
Henning, thanks! I am new to OpenFOAM source code. Could you please explain why there are two sets of RK4values (RK4values and RK4values2)?


Plz ignore this question. I've found the classical 4th RK in Jameson's paper.

Last edited by philem; May 15, 2015 at 11:04.
philem is offline   Reply With Quote

Old   June 18, 2015, 12:12
Join Date: Nov 2011
Posts: 32
Rep Power: 14
Andy_bm is on a distinguished road
Thank you for your code.

I have been working with rhoCentralFoam solver several years. And I have some problem. Does anyone has the follow problem: rhocentralFoam is good working in 3d cases, BUT only with UPWIND reconstruct scheme? If I change to vanleer scheme or other(gamma, superbee) the solver crashed. Also I tried different grad scheme but it has no result. This problem is present in all my tasks: external flow at M~2-4 or jet stream from nozzle.

Andy_bm is offline   Reply With Quote

Old   July 20, 2016, 02:58
New Member
Marko Josic
Join Date: Dec 2010
Posts: 20
Rep Power: 15
emjay is on a distinguished road

i have/had the same problem. Does anybody know the discretization settings for a second order in time and space?

with "backward" the solver crashes, and Euler is only first order and for my unsteady problem to diffusive.

and for spatial dissertation is the same problem, only "upwind" works, any other higher order discretization with limiter function produces overshoots, when i compare the numerical with the analytical solution.
emjay is offline   Reply With Quote

Old   November 3, 2016, 09:01
New Member
Wenkun Zhao
Join Date: Mar 2015
Location: Nanjing, China
Posts: 14
Rep Power: 11
Eric Brant is on a distinguished road
Hi Henning,

I am quite interested in your rhoCentralFoamRK4 and how about the results compared to rhoCentralFoam? Have you ever calculated some cases and compared with other results obtained by FLUNET etc.

Eric Brant is offline   Reply With Quote

Old   December 13, 2016, 10:12
Default RK4 - rhoCentralFoam
Senior Member
Join Date: Sep 2015
Location: Singapore
Posts: 102
Rep Power: 10
usv001 is on a distinguished road
Hello everyone,

I am trying to implement RK4 for rhoCentralFoam as well. I want to clarify a couple of issues pertaining to Henning86's implementation first.
  1. Can RK4 only be applied in a fully explicit equation? I ask this because in the uploaded solvers, the 'fvm::laplacian()' operators have been changed to 'fvc::laplacian()'.
  2. In a completely explict equation, there is no need to have separate predictor and corrector steps which is how this current solver is coded I believe, However, wouldn't this negatively affect accuracy?

On a separate note, I would like to know whether it is feasible to implement RK4 as a standalone ddtScheme like EulerDdtScheme because, so far, I don't think it is possible.

It would be extremely helpful if anyone could enlighten me on these issues.

Many thanks,
usv001 is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Main advantage of using Runge Kutta of higher order? jakubstary Main CFD Forum 14 August 20, 2019 16:15
1D Burgers euqation with 4th Runge Kutta dokeun Main CFD Forum 3 August 8, 2011 06:34
Runge Kutta Discontinuous Galerkin procedure? Ameya J Main CFD Forum 6 July 29, 2011 05:56
Runge Kutta Optimization vasanth Main CFD Forum 6 December 2, 2005 13:07
Diagonally Dominate Runge Kutta Method Anthony Iannetti Main CFD Forum 0 January 23, 2001 21:27

All times are GMT -4. The time now is 22:47.