# Could someone explain the implementation of the convection scheme

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

 January 20, 2008, 20:18 It seems that something wrong #1 New Member   Liu Huafei Join Date: Mar 2009 Location: Shanghai, China Posts: 20 Rep Power: 10 Sponsored Links 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 Φf =ΦN+(1-fx)*ψ(rf) * (ΦP-ΦN) =(1-fx) *ψ(rf)* ΦP +[1-(1-fx) *ψ(rf)]* * Φ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 tmp limitedSurfaceInterpolationScheme::weights ( const GeometricField& phi ) const { const fvMesh& mesh = this->mesh(); tmp tWeights(this->limiter(phi)); this isψ(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?

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Romuald Skoda Main CFD Forum 10 August 3, 2010 12:23 kallipygian OpenFOAM Running, Solving & CFD 0 October 13, 2008 07:29 alberto OpenFOAM Running, Solving & CFD 0 March 13, 2008 22:18 Markus Main CFD Forum 4 January 6, 2007 02:53 Andy R FLUENT 2 December 13, 2006 06:16