# Prediction of gradient at the cell center for cell-centered finite volume method

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

 March 23, 2014, 06:58 Prediction of gradient at the cell center for cell-centered finite volume method #1 Senior Member   Join Date: Jan 2013 Posts: 144 Rep Power: 3 Dear CFDer, As you know, when I implement cell-centered finite volume discretization, I need to calculate the gradient at each cell center, which will be used in the numerical schemes for convection and diffusion terms. I found that this will slow down my simulations because before I start the face looping of numerical fluxes for the convection and diffusion terms, the gradients should be known in advance. I use Green-Gauss method to predict the cell center gradients, which means that the complete face looping is also needed. So, in this sense, I need to do two cycle of face sweepings: one for gradient (this should be performed first) and another for numerical fluxes. This obviously slow down my code. Since I guess this is classic problem in FV CFD, does anybody have some suggestions about the appraoches to efficiently calculate the cell ceter gradients without loss of computation speed? Thank you so much, OFFO

 March 23, 2014, 09:37 #2 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 423 Blog Entries: 13 Rep Power: 14 The standard Green-Gauss based gradient can clearly be computed trough a face loop, you just need to... do it. For each face in the loop you just need to update all the relevant gradient terms in the two cells sharing the face. Obviously, you will still need to compute the gradient in a preliminary face loop before using it. The same should be possible even for the least-squares approach (probably more cumbersone)

March 23, 2014, 09:43
#3
Senior Member

Join Date: Jan 2013
Posts: 144
Rep Power: 3
Thank you!

I did loop all the face first, which need computational time. This is my concern. Because I need loop all faces first and then do the looping again to calculate the numerical fluxes.

Thank you.

Quote:
 Originally Posted by sbaffini The standard Green-Gauss based gradient can clearly be computed trough a face loop, you just need to... do it. For each face in the loop you just need to update all the relevant gradient terms in the two cells sharing the face. The same should be possible even for the least-squares approach (probably more cumbersone)

 March 23, 2014, 15:59 #4 Senior Member     Paolo Lampitella Join Date: Mar 2009 Location: Italy Posts: 423 Blog Entries: 13 Rep Power: 14 I'm sorry, i don't know of any other obvious (or not) way to do it. Consider that, in general, you should do it more than twice to keep with the iterations. How much slow down are we talking about?

 March 23, 2014, 16:07 #5 Senior Member   Join Date: Jan 2013 Posts: 144 Rep Power: 3 Yes, you are right. Thank you.

 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 zhengjg Main CFD Forum 10 November 11, 2012 23:13 anand_30 OpenFOAM Meshing & Mesh Conversion 12 December 12, 2011 04:16 cheng1988sjtu OpenFOAM Bugs 14 September 18, 2011 13:46 WatchFirefly Main CFD Forum 12 September 12, 2010 20:58 Paul Main CFD Forum 5 June 21, 2003 12:55

All times are GMT -4. The time now is 06:43.

 Contact Us - CFD Online - Top