CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM

pitzDailyMapped killing fluctuations

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

Like Tree1Likes
  • 1 Post By juho

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 30, 2016, 14:12
Default pitzDailyMapped killing fluctuations
  #1
Member
 
Join Date: Apr 2016
Posts: 34
Rep Power: 7
Eman. is on a distinguished road
HI,

I was going through the pitzDailyMapped folder and it is supposed to create a turbulent inflow condition by mapping the velocity at a plane located in the downsteam of the flow back to the inlet plane. I noticed it uses averaging of velocity to maintain the volumetric flux. However, the average vector is (10,0,0), which means it kills the eddies and therefore the fluctuations over time by mapping the Uy and Uz to zero. Isn't it the whole idea of LES to take the fluctuations larger than the filter width into account? I assume the inlet profile would approach that of a laminar flow after enough time steps and if I'm right the whole simulation case of "pitzDailyMapped" makes no sense. Am I missing something here? any idea?

Many thanks in advance.
Eman. is offline   Reply With Quote

Old   January 2, 2017, 08:06
Default
  #2
Member
 
Juho Peltola
Join Date: Mar 2009
Location: Finland
Posts: 89
Rep Power: 14
juho is on a distinguished road
Here you can see how the setAverage works:

https://github.com/OpenFOAM/OpenFOAM...tchFieldBase.C

Code:
   if (setAverage_)
    {
        Type averagePsi =
            gSum(patchField_.patch().magSf()*newValues)
           /gSum(patchField_.patch().magSf());

        if (mag(averagePsi)/mag(average_) > 0.5)
        {
            newValues *= mag(average_)/mag(averagePsi);
        }
        else
        {
            newValues += (average_ - averagePsi);
        }
    }
Eman. likes this.
juho is offline   Reply With Quote

Old   January 2, 2017, 17:32
Default
  #3
Member
 
Join Date: Apr 2016
Posts: 34
Rep Power: 7
Eman. is on a distinguished road
Quote:
Originally Posted by juho View Post
Here you can see how the setAverage works:

https://github.com/OpenFOAM/OpenFOAM...tchFieldBase.C

Code:
   if (setAverage_)
    {
        Type averagePsi =
            gSum(patchField_.patch().magSf()*newValues)
           /gSum(patchField_.patch().magSf());

        if (mag(averagePsi)/mag(average_) > 0.5)
        {
            newValues *= mag(average_)/mag(averagePsi);
        }
        else
        {
            newValues += (average_ - averagePsi);
        }
    }
Thanks a lot for the reply. So, it rescales the velocity to maintain the flow rate. However, I tried to use mapped boundary condition on different geometries (2D and 3D) and everytime turbulence decays and I end up with a laminar velocity profile.

Also here on page 15 you can see:
"Care must be taken not to make the control algorithms too efficient; if the velocity control is too good then it can have the effect of completely destroying fluctuations in the velocity, thus unintentionally relaminarising the flow."
Eman. is offline   Reply With Quote

Old   January 3, 2017, 11:44
Default
  #4
Senior Member
 
chegdan's Avatar
 
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 621
Rep Power: 26
chegdan will become famous soon enoughchegdan will become famous soon enough
the way that channelFoam use to work for this type of application was to use a cyclic BC with a pressure gradient to drive the flow in a user-specified direction. Without the pressure gradient -> the flow decays. There is an fvoption in openFOAM-plus for pressure gradient source term to replicate that of channelFoam used like:

Code:
    airDeflection
    {
        type            directionalPressureGradientExplicitSource;
        active          true;

        directionalPressureGradientExplicitSourceCoeffs
        {
            selectionMode   cellZone;
            cellZone        cZone;

            fieldNames  (U);            // Name of the field
            flowDir     (1 1 0);        // Desired flow direction
            faceZone    f0Zone;         // Face zone upstream cell zone
            relaxationFactor    0.3;    // Relaxation factor for flow
                                        // deflection (default 0.3)

            //Pressure drop model [Pa]
            model       volumetricFlowRateTable;//constant;//DarcyForchheimer;

            //DarcyForchheimer model
            // deltaP = (D*mu + 0.5*rho*magUn)*magUn*length_

            D           5e7;
            I           0;
            length      1e-3;

            //constant model
            pressureDrop    40;

            //volumetricFlowRateTable model
            outOfBounds     clamp;
            fileName        "volFlowRateTable";
        }
    }
__________________
Dan

Find me on
Twitter @dancombest
LinkedIn
WordPress
chegdan is offline   Reply With Quote

Reply

Tags
large eddy simulation, mapped boundary condition, mapped condition, pitzdaily

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
How does CFX initialize velocity fluctuations in LES? dtmith CFX 3 June 4, 2017 04:09
Inlet fluctuations issue in LES Ivan Main CFD Forum 0 October 15, 2012 13:57
numerical prediction of Pressure fluctuations in hydraulic jump masoud197 CFX 1 January 22, 2012 18:58
Obtaining RMS of turbulent velocity fluctuations rks171 Main CFD Forum 0 January 19, 2012 09:45
about turbulence random fluctuations ben Main CFD Forum 2 December 23, 2004 21:23


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