# averaging in les

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

 February 5, 2010, 04:21 averaging in les #1 Member   Kapa Lilla Join Date: Mar 2009 Location: Bruxelles, Belgium Posts: 45 Rep Power: 8 Dear All, I made a compressible LES simulation and used the averaging function to get the mean and r.m.s. pressure. A bit disturbing, but the r.m.s. pressure contains negative values in some cells. How can this be? Thanks, Lilla

 February 9, 2010, 05:02 #2 Senior Member   Eugene de Villiers Join Date: Mar 2009 Posts: 725 Rep Power: 12 I don't think the Prime2Mean values are being calculated in the way you think, i.e Field_rms = sqrt(Mean(sqr(Field-FieldMean))) Remember, this is a running average, so the initial value for FieldMean is not accurate. As a result the mean squares are calculated as: Field_ms_new = alpha * Field_ms_old + beta*sqr(Field) - sqr(FieldMean_old) + sqr(FieldMean_new) where alpha = (N-1)/N and beta = 1/N, which can lead to negative values for small N. It could also happen, that you restart the calculation with a FieldMean present, but with no Field_ms. This could lead to large initial negative values, as the initial Field_ms is calculated as: Field_ms_init = sqr(Field) - sqr(FieldMean) A better initialisation would probably be: Field_ms_init = sqr(Field - FieldMean)

 February 10, 2010, 03:29 #3 Member   Kapa Lilla Join Date: Mar 2009 Location: Bruxelles, Belgium Posts: 45 Rep Power: 8 Dear Eugene, thanks for your answer. Indeed, I checked the implementation as well, and the problem can arise from my test case. I have a pulsating source around a constant flowfield. The source term I guess makes harder to compute the mean, so I decided to fix the mean flow by myself and than calculate p^2 for noise computation. Best regards, Lilla

 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 Sungho Yoon CFX 6 July 23, 2008 10:34 iko FLUENT 3 July 3, 2008 05:55 fab Main CFD Forum 2 June 16, 2005 04:53 Bin Li Main CFD Forum 2 February 20, 2004 10:58 Kim Main CFD Forum 0 November 21, 2001 08:18

All times are GMT -4. The time now is 12:33.