CFD Online Discussion Forums - RhoCentralFoam detail

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - RhoCentralFoam detail (https://www.cfd-online.com/Forums/openfoam-solving/58569-rhocentralfoam-detail.html)

Dear Luca

I ran a supersonic flat plate flow with turbulence, but when I change different limiter of:
reconstruct(rho) Minmod;
reconstruct(U) Minmod;
reconstruct(T) Minmod;
such as Minmod, vanLeer, and superBee, the skin friction is quite differnt from each other, the minmod gave the highest value, while the superBee gave the lowest.
As for a supersonic flat plate, there is no shock in the boundary layer, i do not know why there is so much difference, when I look into the U file alone y direction, near the wall the limiters influence it. Do you noticed that when you did some simulation?
regards,
Jason

Hi, foamers,
when I search the "neg", the officical document show this:

Quote:

neg
Foam::dimensionSet::neg()
Foam::neg(const dimensionedScalar &ds)
Foam::neg(const dimensionSet &)
Foam::neg(const label s)
Foam::neg(const Scalar s)
negate
negIndexFoam::BiIndirectList
negListFoam::BiIndirectList
negPart
negSumDiag

where is the function using field?

Quote:

Originally Posted by david_h (Post 197407)

Alberto,

I agree with the above statements, if the "face area vectors" point outward from cell i,j.

If all the "face area vectors" are pointed either toward the right or to the top, then:

on face the between i-1 and i (i-1/2),
"pos" would return rho[i-1]
"neg" would return rho[i]
on face the between i and i+1 (i+1/2),
"pos" would return rho[i]
"neg" would return rho[i+1]

The above is for zero-order reconstruction. For "cell-based" higher order reconstruction (e.g. limited-linear), the face values would take the form:
rho_pos[i+1/2] = rho[i] + drho[i]*(x[i+1/2] - x[i])
rho_neg[i+1/2] = rho[i] + drho[i]*(x[i+1/2] - x[i+1])

where drho[i] is limited gradient of rho at cell "i"

Dave

Good morning,
sorry to reopen this old thread, but I'm studying rhoCentraFoam and I would like to understand better how this solver works.

I can't understand when David writes:

Quote:

rho_neg[i+1/2] = rho[i] + drho[i]*(x[i+1/2] - x[i+1])

I think "negative interpolation" interpolates from x[i+1], so that:
rho_neg[i+1/2] = rho[i+1] + drho[i+1]*(x[i+1/2] - x[i+1])
Is it corret, or am I missing something?

In general, refers to the article:
"Implementation of semi-discrete, non-staggered central schemes in a colocated, polyhedral, finite volume framework, for high-speed viscous flows" by C.J. Greenshields et Al.
it is explained how to compute positive interpolation of a general quantity:
$\Psi_{f+}=(1-g_{f+})*\Psi_{P}+g_{f+}*\Psi_{N}$

where:
$g_{f+}=\beta*(1-w_f)$
However I can't understand how to compute $\Psi_{f-}$ .
I hope someone can help me.
Thank you