Foam::lduMatrix::faceH

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 September 22, 2009, 05:21 Foam::lduMatrix::faceH #1 New Member   Yijin Li Join Date: Mar 2009 Posts: 13 Rep Power: 9 Based on the pressure solution, assemble conservatie face flux F F=s_f*H(u)-a_p_n(p_n-p_p) Correct face flux phi -= pEqn.flux(); I found that in the code faceHpsi[face] = Upper[face]*psi[u[face]] - Lower[face]*psi[l[face]]; why use this expression to correct face flux? //////////////////////////////////////////////////////////////// Foam::lduMatrix::faceH template Foam::tmp > Foam::lduMatrix::faceH(const Field& psi) const { tmp > tfaceHpsi(new Field (lduAddr().lowerAddr().size())); Field& faceHpsi = tfaceHpsi(); if (lowerPtr_ || upperPtr_) { const scalarField& Lower = const_cast(*this).lower(); const scalarField& Upper = const_cast(*this).upper(); // Take refereces to addressing const unallocLabelList& l = lduAddr().lowerAddr(); const unallocLabelList& u = lduAddr().upperAddr(); for (register label face=0; face::zero; } return tfaceHpsi; } //////////////////////////////////////////////////////////////// fumiya likes this.

 September 23, 2009, 03:41 #2 Senior Member   Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,783 Rep Power: 22 Because of consistency: since the pressure equation IS the continuity equation, fuxes must be built in exactly the same way as the pressure matrix coefficients and sources. Work out the maths by hand and you will see that what you need is the off-diagonal pressure coefficient (+ explicit corrections, eg. non-orthogonality, which are treated separately). Hope this helps, Hrv fumiya likes this. __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

 September 23, 2009, 08:00 #3 New Member   Yijin Li Join Date: Mar 2009 Posts: 13 Rep Power: 9 Thank you for your help. I see.

 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

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

 Contact Us - CFD Online - Top