[srinath]

Hello

I would like to write a Riemann solver for Euler's equations, as i have all the usual functions for calculating the intercell fluxes.
But i am unsure about how to implement it in openFoam.

I need to do the following
1)Loop over all internal faces
2)Find the 2 cells sharing a particular face
3)Extract values such as u,rho,p etc from these 2 cells and assign a flux to that face(after solving a Riemann problem)
But how would i do this in OpenFoam?

Just a couple of code snippets and some info on where to look in the source code should suffice.
None of the existing solvers, seem to do anything like this

Thanks
Srinath

[pbohorquez]

Did you have a look into g3dFoam? I think this can help you.

g3dFoam: Implementing third order compressible flow solver for hexahedral meshes in OpenFoam, by Martin Olausson.

http://www.tfd.chalmers.se/~hani/kurser/OF_phD_2007/

[srinath]

Thanks Patricio

I will take a look.

Srinath