Density-based, fully-coupled, compressible solver in OpenFOAM?
I would like to discuss and to share opinions with all you Foamers out there about the lack of a density-based, fully coupled, compressible solver in OpenFOAM.
I've been studying OpenFOAM for about a month and, as far as I have understood, the only density-based solver (by "density-based" I mean that density is calculated from continuity equation) is rhoCentralFoam, which however solves the governing equations separately one from each other (in a so-called segregated way). All the other solvers are of pressure-based type. First of all: is this true?
Second: how hard do you think it would be to create a new density-based solver, which solves the governing equations simultaneously (i.e. a coupled solver), provided with a robust upwind-like scheme such as Roe's or AUSM?
Third: I have found, searching the Forum, that few years ago two Italian students developed a solver called aeroFoam
which apparently satisfies all the features that I mentioned. However, the solver has not been officially included into one of the OpenFOAM releases and the project has been apparently abandoned. Natural question: why? :D
Hope the discussion will be constructive,
In addition, the rest of OpenFOAM would be need to made of aware of these new classes, e.g. the linear solvers.
There has been some ongoing development in the -extend branch of OpenFOAM to add block coupled matrix classes. I don't believe a coupled N-S stokes solver has been built yet with it but I have seen work that uses it to couple other quantities. Dig around the last couple of workshops for material.
So as you can imagine, the work is significant and non-trivial.
thanks for your answer. I realize the strong effort needed to extend OpenFOAM towards a coupled strategy. For those interested, the subject is studied in more detail in the following slides:
On the other hand, I think that implementing an upwind-like flux scheme in rhoCentralFoam, for instance, should be much easier. Any experience/opinion?
One other effort of note is some density-based (segregated) solvers which include some more modern convective schemes, e.g. AUSM, etc, for turbomachinery in the extend repo:
I have seen very little documentation about it's use (presentatations or otherwise) but it's there for you to explore.
Very interesting, thank you very much!
|All times are GMT -4. The time now is 12:57.|