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

Could someone explain the implementation of the convection scheme

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

LinkBack Thread Tools Search this Thread Display Modes
Old   January 20, 2008, 20:18
Default It seems that something wrong
New Member
Liu Huafei
Join Date: Mar 2009
Location: Shanghai, China
Posts: 20
Rep Power: 16
liuhuafei is on a distinguished road
It seems that something wrong with the implementation of the convection scheme. I guess the face value is based on the following formula

Φf =ΦC +(xf – xC) *ψ(rf) * (dΦ/ dx)f
Φf =ΦC +(xf – xC) / (xD– xC) *ψ(rf) * (ΦD-ΦC)

f ,the face
C,central node
D downwind node
ψ(rf) the flux limiter

rf is based on the following formula
rf=(ΦC-ΦU)/( ΦD-ΦC) = 2*grad(ΦC)*(xD-xC) /( ΦD-ΦC) -1

for the faceflux mf
mf>0 Φf =ΦP+fx*ψ(rf) * (ΦN-ΦP)
=(1-fx*ψ(rf)) ΦP +f x*ψ(rf)* ΦN
mf<0 &Phi;f =&Phi;N+(1-fx)*&psi;(rf) * (&Phi;P-&Phi;N)
=(1-fx) *&psi;(rf)* &Phi;P +[1-(1-fx) *&psi;(rf)]* * &Phi;N

fx is the geometric interpolation factor, defined by fx=(xN-xf)/(xN-xP)
P is owner cell
N is neighbour cell

but in the limitedSurfaceInterpolationScheme.H and limitedSurfaceInterpolationScheme.c
tmp<surfacescalarfield> limitedSurfaceInterpolationScheme<type>::weights
const GeometricField<type,>& phi
) const
const fvMesh& mesh = this->mesh();
tmp<surfacescalarfield> tWeights(this->limiter(phi)); this is&psi;(rf)!!!!!!!
surfaceScalarField& Weights = tWeights();
const surfaceScalarField& CDweights = mesh.surfaceInterpolation::weights(); fx !!!! scalarField& pWeights = Weights.internalField();
forAll(pWeights, face){
pWeights[face] = pWeights[face]*CDweights[face]
+ (1.0 - pWeights[face])*pos(faceFlux_[face]); It is right??????
Maybe openfoam is based on other formula. Could you kindly explain clearly it?
liuhuafei 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
Implementation of QUICK scheme Romuald Skoda Main CFD Forum 11 November 6, 2017 22:20
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 08:29
Deciphering interpolation scheme implementation alberto OpenFOAM Running, Solving & CFD 0 March 13, 2008 22:18
FEM Implementation of pressure-correction scheme Markus Main CFD Forum 4 January 6, 2007 02:53
Scheme / TUI implementation Andy R FLUENT 2 December 13, 2006 06:16

All times are GMT -4. The time now is 02:27.