CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Could someone explain the implementation of the convection scheme

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 20, 2008, 20:18
Default It seems that something wrong
  #1
New Member
 
Liu Huafei
Join Date: Mar 2009
Location: Shanghai, China
Posts: 20
Rep Power: 8
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
Or
Φ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

so
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
template<class>
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

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Implementation of QUICK scheme Romuald Skoda Main CFD Forum 10 August 3, 2010 12:23
Please explain the implementation of species transport Eqn in reactingFoam kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07: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 06:10.