
[Sponsors] 
November 14, 2011, 02:42 
rhoCentralFoam linearUpwind

#1 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 8 
Dear foamers.
I want to try the schemes linearUpwind in interpolationSchemes for rhoCentralFoam solver. In the begining my fvSchemes file looks so: Code:
fluxScheme Tadmor; ddtSchemes { default Euler; } gradSchemes { default cellLimited leastSquares 1; } divSchemes { default none; div(tauMC) Gauss linearUpwind; div(phi,k) Gauss linearUpwind; div(phi,omega) Gauss linearUpwind; } laplacianSchemes { default Gauss linear corrected; } interpolationSchemes { default linear; reconstruct(rho) upwind; reconstruct(U) upwind; reconstruct(T) upwind; } snGradSchemes { default corrected; } Code:
interpolationSchemes { default linear; reconstruct(rho) linearUpwind; reconstruct(U) linearUpwind; reconstruct(T) linearUpwind; } Code:
[3] > FOAM FATAL IO ERROR: [3] Grad scheme not specified Valid grad schemes are : 8 ( Gauss cellLimited cellMDLimited extendedLeastSquares faceLimited faceMDLimited fourth leastSquares ) Code:
interpolationSchemes { default linear; reconstruct(rho) Gauss linearUpwind; reconstruct(U) Gauss linearUpwind; reconstruct(T) Gauss linearUpwind; } Code:
[17] > FOAM FATAL IO ERROR: [17] Unknown discretisation scheme Gauss Valid schemes are : 51 ( Gamma Gamma01 MUSCL MUSCL01 Minmod OSPRE QUICK SFCD SuperBee UMIST biLinearFit blended clippedLinear cubic cubicUpwindFit downwind filteredLinear filteredLinear2 filteredLinear3 fixedBlended harmonic limitWith limitedCubic limitedCubic01 limitedGamma limitedLimitedCubic limitedLimitedLinear limitedLinear limitedLinear01 limitedMUSCL limitedVanLeer linear linearFit linearPureUpwindFit linearUpwind localBlended localMax localMin midPoint outletStabilised quadraticFit quadraticLinearFit quadraticLinearUpwindFit quadraticUpwindFit reverseLinear skewCorrected upwind vanAlbada vanLeer vanLeer01 weighted ) Thanks. 

November 14, 2011, 04:46 

#2 
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 12 
To use the linearUpwind scheme you have to specify also a gradient scheme. The correct syntax in OF2.0.x should be (for instance):
div(phi,k) Gauss linearUpwind grad(k); where grad(k) is the scheme specified for k in the gradSchemes section (in your case, the default cellLimited leastSquares 1 scheme will be selected). In OF1.7.x this syntax will not be recognized, because you have to specify explicitly the gradient scheme, like in the following example: div(phi,k) Gauss linearUpwind cellLimited leastSquares 1; Finally, just a note about div(tauMC): it is an explicitly calculated divergence term, not a convective flux term, thus the only scheme that makes sense is Gauss linear. Hope this helps V. 

November 14, 2011, 11:11 

#3 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 8 
Thanks for your fast reply,
You tip has really helped me, but i have same questions: If i use Code:
interpolationSchemes { default linear; reconstruct(rho) upwind cellLimited leastSquares 1; reconstruct(U) upwind cellLimited leastSquares 1; reconstruct(T) upwind cellLimited leastSquares 1; } 

November 14, 2011, 11:30 

#4  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 12 
Quote:
Code:
interpolationSchemes { default linear; reconstruct(rho) Gamma 1; reconstruct(U) GammaV 1; reconstruct(T) Gamma 1; } V. 

November 14, 2011, 11:43 

#5 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 8 
I already tried to use this:
Code:
interpolationSchemes { default linear; reconstruct(rho) Gamma 1; reconstruct(U) GammaV 1; reconstruct(T) Gamma 1; } my case has the decision only Code:
interpolationSchemes { default linear; reconstruct(rho) upwind; reconstruct(U) upwind; reconstruct(T) upwind; } Thanks. 

November 14, 2011, 11:51 

#6  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 12 
Quote:
V. 

November 14, 2011, 12:02 

#7 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 8 
I solve a problem of the expiration stream (М~25) in vacuum what solver can you recommend? I used sonicFoam, but he required too much time for the decision.


November 14, 2011, 12:27 

#8  
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 12 
Quote:
V. 

November 14, 2011, 12:34 

#9 
New Member
Join Date: Nov 2011
Posts: 28
Rep Power: 8 
Thanks a lot, I will try select schemes.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
A few simple questions about linearUpwind and limitedLinear  chegdan  OpenFOAM  16  October 13, 2015 03:59 
rhoCentralFoam with totalPressure/totalTemperature at inlet of subsonic channel  deepblue17  OpenFOAM Running, Solving & CFD  5  February 11, 2013 03:42 
localEuler in rhoCentralFoam  praveen  OpenFOAM  4  October 21, 2012 02:19 
linearUpwind scheme in OpenFOAM 2.0.1 ???  cabul  OpenFOAM  8  November 9, 2011 07:57 
Always crash when solve a CD nozzle flow field using rhoCentralFoam  hawklion  OpenFOAM Running, Solving & CFD  0  March 9, 2011 07:13 