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

where does this term come from

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 12, 2013, 04:43
Default where does this term come from
  #1
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 153
Rep Power: 8
Jonathan is on a distinguished road
Hi everyone,

I am solving a turbine case in OF2.1.1 using SRFSimpleFoam using kOmegaSST. I am having some bounding and convergence issues but am confused about a particular term in the divSchemes subdictionary in fvSchemes:

Does anyone know which equation:
Code:
div((nuEff*dev(T(grad(Urel))))) Gauss upwind; //linear;
comes from and why i cannot using Gauss upwind to discretize it - I get fatal read beyond EOF error if i dont use linear with it.

thanks very much in advance
jonathan
Jonathan is offline   Reply With Quote

Old   June 12, 2013, 06:00
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,205
Rep Power: 17
immortality is on a distinguished road
Hi Jonathan
I had a problem like this with tauMC in rhoCentralFoam and Bruno said me that its a tensor field not a scalar or vector field that can use upwind schemes.
if you can get more information or others its nice to share.
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   June 12, 2013, 08:51
Default from kOmegaSST turbulence model
  #3
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 153
Rep Power: 8
Jonathan is on a distinguished road
hi,

so the term comes from the turbulence model (kOmegaSST):

From kOmegaSST.H
Code:
        //- Return the source term for the momentum equation
        virtual tmp<fvVectorMatrix> divDevReff(volVectorField& U) const;
From kOmegaSST.C
Code:
tmp<fvVectorMatrix> kOmegaSST::divDevReff(volVectorField& U) const
{
    return
    (
      - fvm::laplacian(nuEff(), U)
      -  fvc::div(nuEff()*dev(T(fvc::grad(U))))
    );
}
As immortality says, this might be a tensor field, and so the standard divSchemes as per the UG do not apply, but i am still not 100% sure.
Jonathan is offline   Reply With Quote

Old   June 16, 2013, 11:36
Default
  #4
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,511
Blog Entries: 34
Rep Power: 86
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

@Jonathan: The post Ehsan was talking about is this one: error on linearUpwind - post #4.
And yes, this is basically the same pattern of mathematics going on in this case.

If you run the following command, you'll see that all tutorials in OpenFOAM use "Gauss linear" for this particular kind of div:
Code:
find $FOAM_TUTORIALS -name fvSchemes | xargs grep "T(grad(U"
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   June 17, 2013, 07:52
Default thanks
  #5
Senior Member
 
Join Date: Mar 2010
Location: Cape Town, SA
Posts: 153
Rep Power: 8
Jonathan is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Greetings to all!

@Jonathan: The post Ehsan was talking about is this one: error on linearUpwind - post #4.
And yes, this is basically the same pattern of mathematics going on in this case.

If you run the following command, you'll see that all tutorials in OpenFOAM use "Gauss linear" for this particular kind of div:
Code:
find $FOAM_TUTORIALS -name fvSchemes | xargs grep "T(grad(U"
Best regards,
Bruno
Hi Bruno,

Thanks very much - your help is much appreciated.

Best regards
Jon
Jonathan is offline   Reply With Quote

Reply

Tags
divschemes, equation, fvschemes

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
momentum source term zwdi FLUENT 13 December 5, 2013 18:58
How does fluent handles the momentum source term? QBeast FLUENT 0 April 22, 2013 14:12
Source Term Question Blackhawks84 Fluent UDF and Scheme Programming 1 July 31, 2012 15:44
Large source term in species equation MACFD FLUENT 4 January 4, 2011 15:16
bouyancy term in epsilon equation Michael Main CFD Forum 1 June 25, 1999 10:20


All times are GMT -4. The time now is 21:10.