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

bounded Gauss and normal Gauss question

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By Tobi

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 11, 2013, 17:22
Default bounded Gauss and normal Gauss question
  #1
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

I modifyed the flamelet solver (libOpenSmoke) and used the nice feature in the OF 2.2.x for the equations for

  • mixture fraction Z
  • variance of mixture fraction Z
  • enthalpy H


In the old code I implemented the Sp term in the code:


Code:
        (
            fvm::div(phi, Z)
                    - fvm::Sp(fvc::div(phi), Z)        
                        - fvm::laplacian(turbulence->muEff()/sigmat, Z)
        )
In the new version I used it without:

Code:
        (
            fvm::div(phi, Z)    
                        - fvm::laplacian(turbulence->muEff()/sigmat, Z)
        )
but therefore I added bounded schemes (fvSchemes):
Code:
    div(phi,H)          bounded Gauss limitedLinear 1;
    div(phi,Z)          bounded Gauss limitedLimitedLinear 1 0 1;
    div(phi,Zvar)       bounded Gauss limitedLimitedLinear 1 0 0.25;
This settings give an identical simulation and the validation was 1:1.

OF was warning like:
Code:
--> FOAM Warning : 
    From function gaussConvectionScheme
    in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 123
    Reading "/home/shor-ty/OpenFOAM/shor-ty-2.2.x/flameletModel-2.2.x/tutorials/Sandia_COH2N2/system/fvSchemes.divSchemes.div(phi,U)" at line 30
    Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
    To remove this warning switch off 'boundedGauss' in "/home/shor-ty/OpenFOAM/OpenFOAM-2.2.x/etc/controlDict"
smoothSolver:  Solving for Ux, Initial residual = 0.00265791, Final residual = 1.16588e-05, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 3.76559e-05, Final residual = 1.86524e-07, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 0.0056559, Final residual = 4.24363e-05, No Iterations 2
GAMG:  Solving for p, Initial residual = 0.00227516, Final residual = 7.45788e-05, No Iterations 2
time step continuity errors : sum local = 0.0252954, global = 9.16694e-05, cumulative = -2.81028
rho max/min : 1.21084 0.176765
smoothSolver:  Solving for Z, Initial residual = 5.02383e-05, Final residual = 6.15323e-07, No Iterations 2
smoothSolver:  Solving for Zvar, Initial residual = 4.07254e-05, Final residual = 4.37422e-08, No Iterations 2
smoothSolver:  Solving for H, Initial residual = 5.61876e-05, Final residual = 7.14374e-07, No Iterations 2
--> FOAM Warning : 
    From function gaussConvectionScheme
    in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 123
    Reading "/home/shor-ty/OpenFOAM/shor-ty-2.2.x/flameletModel-2.2.x/tutorials/Sandia_COH2N2/system/fvSchemes.divSchemes.div(phi,epsilon)" at line 31
    Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
    To remove this warning switch off 'boundedGauss' in "/home/shor-ty/OpenFOAM/OpenFOAM-2.2.x/etc/controlDict"
smoothSolver:  Solving for epsilon, Initial residual = 0.00011511, Final residual = 5.69949e-07, No Iterations 2
--> FOAM Warning : 
    From function gaussConvectionScheme
    in file finiteVolume/convectionSchemes/gaussConvectionScheme/gaussConvectionScheme.H at line 123
    Reading "/home/shor-ty/OpenFOAM/shor-ty-2.2.x/flameletModel-2.2.x/tutorials/Sandia_COH2N2/system/fvSchemes.divSchemes.div(phi,k)" at line 32
    Unbounded 'Gauss' div scheme used in steady-state solver, use 'bounded Gauss' to ensure boundedness.
    To remove this warning switch off 'boundedGauss' in "/home/shor-ty/OpenFOAM/OpenFOAM-2.2.x/etc/controlDict"
smoothSolver:  Solving for k, Initial residual = 0.000112566, Final residual = 8.72544e-07, No Iterations 2
ExecutionTime = 651.17 s  ClockTime = 654 s
and so I used bounded schemes for all div equations.
If I do that, I get a strange solution (picture).

can someone tell why the values are so different?

picture bound: all schemes are bounded Gauss
picture unbound: only Z,Zvar,H are bounded

Thanks in advance,
Tobi
Attached Images
File Type: jpg bound.jpg (8.9 KB, 131 views)
File Type: jpg unbound.jpg (8.9 KB, 121 views)
the_ichthyologist likes this.

Last edited by Tobi; February 21, 2022 at 07:06.
Tobi is offline   Reply With Quote

Old   June 12, 2013, 04:42
Default
  #2
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi all,

maybe I made a mistake due to too less time and very tierd. I ll have a look at that the next days and Keep you posted!
Tobi is offline   Reply With Quote

Old   June 18, 2013, 01:35
Default
  #3
Senior Member
 
T. Chourushi
Join Date: Jul 2009
Posts: 321
Blog Entries: 1
Rep Power: 18
Tushar@cfd is on a distinguished road
Code:
    div(phi,H)          bounded Gauss limitedLinear 1;     
    div(phi,Z)          bounded Gauss limitedLimitedLinear 1 0 1;     
    div(phi,Zvar)       bounded Gauss limitedLimitedLinear 1 0 0.25;
Hello,

I have a doubt, I am not sure whether I am right or wrong.

I think when you are using the keyword "bounded" and "limited" together the warning might have displayed in the terminal. May be with any other bounded scheme (i.e., instead of limited) this will go off.
Tushar@cfd is offline   Reply With Quote

Reply


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
Explanation bounding k and bounding epsilon idrama OpenFOAM 42 July 13, 2017 05:05
SonicFoam crashes without any ominous signs callumso OpenFOAM Running, Solving & CFD 1 April 25, 2013 06:45


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