|
[Sponsors] |
September 12, 2013, 18:36 |
Implementation of 2nd order upwind scheme
|
#1 |
New Member
Jason Pearl
Join Date: Jul 2013
Location: Burlington
Posts: 10
Rep Power: 12 |
Hi FOAMers,
I want to implement a 2nd order upwind scheme. However i get the error message below whenever I try to use a scheme like QUICK or linearUpwind for my InterpolationSchemes, gradSchemes or my LaplacianSchemes. I can only implement a second order upwind scheme in the divergence section. My fvSchemes file is pasted below Code:
attempt to read beyond EOF file: /home/owner/OpenFOAM/owner-2.2.0/run/Thesis/sonicFoam/Sept12/system/fvSchemes.interpolationSchemes.interpolate(U_0) at line 59. From function ITstream::read(token&) in file db/IOstreams/Tstreams/ITstream.C at line 83. FOAM exiting Code:
ddtSchemes { default Euler; } gradSchemes { default none; grad(U) Gauss upwind phi; grad(p) Gauss upwind phi; grad(e) Gauss upwind phi; } divSchemes { default none; div(phi,U) Gauss upwind phi; div(phid,p) Gauss upwind phi; div(phi,e) Gauss upwind phi; div(phi,K) Gauss upwind phi; div(phiv,p) Gauss upwind phi; div((muEff*dev2(T(grad(U))))) Gauss linear 1; } laplacianSchemes { default none; laplacian(muEff,U) Gauss linear corrected; laplacian(alphaEff,e) Gauss linear corrected; laplacian(Dp,p) Gauss linear corrected; } interpolationSchemes { default none; interpolate(rho) upwind phi; interpolate(rho_0) upwind phi; interpolate((((1|A(U))*rho_0)*U_0)) upwind phi; interpolate(HbyA) upwind phi; interpolate(((1|A(U))*rho_0)) upwind phi; interpolate(thermo:psi) upwind phi; interpolate(U_0) QUICK; } snGradSchemes { default corrected; } fluxRequired { default no; p ; } what I want to do is substitute linearUpwind for all upwind terms in the code as to make it second order and therefore have to deal with less diffusion type inaccuracies. any suggestions?? Am I going about implementing a second order upwind scheme incorrectly? I am using sonicFoam as a solver but I had the same problems with rhoCentralFoam any help would be greatly appreciated. |
|
Tags |
fvschemes, rhocentralfoam, second order accuracy, second order upwind, upwind |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
High Resolution (CFX) vs 2nd Order Upwind (Fluent) | gravis | ANSYS | 3 | March 24, 2011 02:43 |
2nd order conservative schemes | taw | Main CFD Forum | 1 | September 16, 2008 07:05 |
2nd order interpolation for NS solver | Quarkz | Main CFD Forum | 6 | July 15, 2005 21:50 |
1st order temporal & 2nd order spatial | Prateep Chatterjee | FLUENT | 0 | January 19, 2003 00:31 |
second order FD upwind scheme | Heinz Wilkening | Main CFD Forum | 2 | November 3, 1998 14:33 |