
[Sponsors] 
May 30, 2011, 14:37 

#181 
Senior Member
n/a
Join Date: Sep 2009
Posts: 198
Rep Power: 9 
Meant u'T', v'T' .


June 1, 2011, 14:49 

#182 
New Member
Join Date: Apr 2011
Posts: 20
Rep Power: 7 
hi
i guess flow isnt , fully devesloped , and a length needs to be fully developed turbulent flow . 

July 4, 2011, 11:21 

#183  
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 8 
Quote:
hi gaby, i'm confronted with the same problem you had. how did you manage it? i hope you're still visiting the forum from time to time best regards grandgo 

July 5, 2011, 03:32 

#184 
Member
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 80
Rep Power: 9 
Hi,
depends how you want to evaluate the energy spectrum (at a specific point, or over the whole domain). In the dnsFoam the energy spectrum is calculated, so you can simply code your own post processing tool by copying the routines. Code:
#include "fvCFD.H" #include "Kmesh.H" #include "UOprocess.H" #include "fft.H" #include "calcEk.H" #include "graph.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // int main(int argc, char *argv[]) { #include "setRootCase.H" #include "createTime.H" #include "createMeshNoClear.H" #include "createFields.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Info<< nl << "Starting time loop" << endl; instantList timeDirs = timeSelector::select0(runTime, args); Kmesh K(mesh); forAll(timeDirs, timeI) { runTime.setTime(timeDirs[timeI], timeI); Info << "Time =" << runTime.timeName() << endl; Info<< "Reading field U\n" << endl; volVectorField U ( IOobject ( "U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh ); calcEk(U, K).write(runTime.timePath()/"Ek", runTime.graphFormat()); Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s" << " ClockTime = " << runTime.elapsedClockTime() << " s" << nl << endl; } Info<< "End\n" << endl; return 0; } gregor Last edited by gregor; July 11, 2011 at 04:37. Reason: corrected code to create U for every time dir 

July 5, 2011, 04:55 

#185  
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 8 
Quote:
hi gregor, thanks for your post. i'm simulating a pipe flow with LES and i want to measure the energy spectrum "one dimensional" over the radius of the pipe. in dnsFoam the whole domain ist calculated if i'm right (?). so i could sample the E(k) data with a function in the controlDict file, right? best regards, grandgo 

July 5, 2011, 05:28 

#186  
Member
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 80
Rep Power: 9 
Quote:
Quote:
gregor 

July 5, 2011, 06:01 

#187  
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 8 
Quote:
one more: in what way do i have to change the code from above, so that the utility evaluates the energy spectrum of one specific point? EDIT: i used your code to compile a utility but i get this message: Code:
> FOAM FATAL ERROR: calculated number of cells is incorrect From function Kmesh::Kmesh(const fvMesh& mesh) in file Kmesh/Kmesh.C at line 87. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) in "/sw/OpenFOAM/OpenFOAM1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/sw/OpenFOAM/OpenFOAM1.7.x/lib/linux64GccDPOpt/libOpenFOAM.so" #2 Foam::Kmesh::Kmesh(Foam::fvMesh const&) in "/sw/OpenFOAM/OpenFOAM1.7.x/lib/linux64GccDPOpt/librandomProcesses.so" #3 in "/home/stss8/OpenFOAM/stss81.7.x/applications/bin/linux64GccDPOpt/energySpec" #4 __libc_start_main in "/lib64/libc.so.6" #5 at /usr/src/packages/BUILD/glibc2.11.2/csu/../sysdeps/x86_64/elf/start.S:116 any idea? Last edited by grandgo; July 5, 2011 at 09:13. 

July 6, 2011, 03:55 

#188  
Member
Gregor Olenik
Join Date: Jun 2009
Location: http://greole.github.io/
Posts: 80
Rep Power: 9 
Quote:
EDIT: The error message is because your domain is not a box. If look at the KMesh.C file it is counting all cells along every direction. Then it multiplies the cell numbers and if they do not match the total number of cells it raises an error. gregor Last edited by gregor; July 6, 2011 at 04:21. 

July 6, 2011, 04:29 

#189  
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 8 
Quote:
yes, it was my worry, that it has something do to with my domain not being cubeshaped. thanks anyway. best regards grandgo 

July 8, 2011, 14:03 

#190 
Member
Join Date: Oct 2010
Location: Stuttgart
Posts: 35
Rep Power: 8 
hi foamers,
i need your opinions. i have a LES case with a domain of about 6.5e6 cells, reynolds number is 23000, the courant number is about 1, deltaT 1e6 sec . i'm using a modified pisoFoam parallel with 16 cores, 2,4ghz each. the calculation time for one time step is about 200 seconds. i'm wondering if the calculation time is acceptable or not. what do you think? within reason or do i have to change anything? best regards grandgo 

July 9, 2011, 02:55 

#191  
Member
Bernhard Grieser
Join Date: Mar 2010
Location: Zurich, Switzerland
Posts: 30
Rep Power: 9 
Quote:
The PISO algorithm requires lower CFL numbers, so if I were you I'd reduce the time step significantly until maxCourant is below 0.1, maybe even less. You're modeling a pipe flow, don't you? For that given Reynolds number you can go coarser on the grid; maybe reducing the total cell count by a factor of two would be sufficient. Also, if possible, I recommend you to double your proc number to 32. Which sgsModel are you using, btw? 

July 11, 2011, 17:14 
Regarding u'^2, T'^2 and Yi'^2

#192 
Member
N. A.
Join Date: May 2010
Posts: 64
Rep Power: 8 
Hello Folks,
I think my question is suited for this post and I am hoping I get to learn more about some of the postprocessing of variables in openFOAM. Following are my questions regarding u'^2, T'^2 andYi'^2. 1) Is u'^2 that we obtained from prime2Mean same as the u_sgs^2? Similarly is T'^2 that we get from prime2Mean also same as T_sgs^2? 2) If not, then how are Tprime2Mean related to T_sgs? 3) I am trying to compare <T_measured> (time averaged mean values) with simulations and heance I compare it with T_mean as obtained from mean suing fieldAverage calculation in OpenFOAM. I am wondering if this is an appletoapple comparison, bececause T_mean from openFOAM is mean of the resolved scale values and not the mean of instantaneous values as calulated for measurements. Please share your thoughts and I think this is an important point for folks using LES and comparison with measurements. Thanks, Nir Last edited by N. A.; July 12, 2011 at 11:06. 

July 21, 2011, 22:58 
hi

#193 
New Member
Join Date: Apr 2011
Posts: 20
Rep Power: 7 
hi fomers
i want to simulate , turbulent flow in a square channel with oodles solver , i use LES model with smagorinsky , i want to know what i set for boundary condition for nu Sgs in inlet , outlet , wall , and symmetrry plane ? is zerogradient suitable ? 

August 25, 2011, 09:21 

#194  
Member

Quote:
really thanks for your post. I used LES and I want to draw E(k)K, but my domian is not a box (a cylinder). how can I use the code in dnsfoam? is it possible? I would appreciate any help from you. best regards. 

March 13, 2012, 07:11 
perturbU_compilation

#195  
Member
supercommandodhruv
Join Date: Sep 2011
Posts: 57
Rep Power: 7 
Hi Alberto,
I downloaded your version of perturbU since I was having the same problem as you described below. However, now it throws up a similar error saying Quote:
1) What is the mathematical constant? 2) How can I implement it in the code and compile it. 3) What is the significance of the constant? Thanks in advance, Dhruv. Quote:


June 12, 2012, 12:14 
LES OneEqEddy in OF 2.1.1

#196 
Member
achinta
Join Date: May 2010
Location: Sydney
Posts: 66
Rep Power: 8 
Hello everyone,
I am using OF2.1.1 OneEqEddy LES model for my case having 6.5 million cells. I imported the mesh from fluent. checkMesh said the mesh is OK but it has non orthogonality:  Checking geometry... Overall domain bounding box (1.47 0.17 0.19) (0.384 0.17 0.19) Mesh (nonempty, nonwedge) directions (1 1 1) Mesh (nonempty) directions (1 1 1) Boundary openness (9.64958e17 2.34786e15 6.13025e16) OK. Max cell openness = 4.58304e15 OK. Max aspect ratio = 299.485 OK. Minumum face area = 3.41806e11. Maximum face area = 0.000249108. Face area magnitudes OK. Min volume = 1.22828e14. Max volume = 4.16394e07. Total volume = 0.20005. Cell volumes OK. Mesh nonorthogonality Max: 87.0638 average: 19.2898 *Number of severely nonorthogonal faces: 18931. Nonorthogonality check OK. <<Writing 18931 nonorthogonal faces to set nonOrthoFaces Face pyramids OK. Max skewness = 2.88371 OK. Coupled point location match (average 0) OK. Mesh OK.  I gave appropriate boundary conditions. I mention some of them which could be wrong(please correct me). 1) There are 4 mass flow inlets where i gave the following in /0/U file:  type flowRateInletVelocity; flowRate constant <value>; value uniform (0 0 0);  There is a velocity inlet for which i gave conditions as usual. type fixedValue; value uniform (<value> 0 0);  2) pressure: All patches/walls except outlet(p=0) have 0 gradient. 3) B Internal field, velocity and mass flow inlets have 0 value. exit and walls have 0 gradient. 4) k (turbulent kinetic energy) Internal field has k=0.8 (same as velocity inlet) mass flow inlets(assumed to have negligible turbulence) have k=1e2 (arbitrary value) Velocity inlet has k=0.8 exit: 0 gradient walls: {type kqRWallFunction; value uniform 0.8;} 5)nuSgs internalfield: 1e7 (arbitrary value) Inlets and exit: 0 gradient Walls: { type nuSgsUSpaldingWallFunction; value uniform 1e7; }  For delta, i am using 'cubeRootVol' fvSchemes and fvSolution are the same as the one used in TUTORIALS "incompressible/pisoFoam/les/pitzDaily". I am running it on 4 nodes. After few time steps it diverges. k goes unbound and residuals of p increase:  Time = 1e09 Courant Number mean: 3.82832e07 max: 0.00098636 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 8.89837e10, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.24966e10, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 1.81297e09, No Iterations 1 DICPCG: Solving for p, Initial residual = 1, Final residual = 0.0492102, No Iterations 841 time step continuity errors : sum local = 3.97898e10, global = 7.55946e12, cumulative = 7.55946e12 DICPCG: Solving for p, Initial residual = 0.156543, Final residual = 1.67279e06, No Iterations 1001 time step continuity errors : sum local = 7.68918e13, global = 7.51877e15, cumulative = 7.56698e12 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 1.83321e06, No Iterations 1 bounding k, min: 3.30181e05 max: 8.65 average: 1.64454e06 ExecutionTime = 1372 s ClockTime = 1383 s .. .. ... ... Time = 5e09 Courant Number mean: 1.71033e06 max: 0.089907 DILUPBiCG: Solving for Ux, Initial residual = 0.00110495, Final residual = 7.5732e09, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.0157435, Final residual = 2.62783e07, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.0154086, Final residual = 3.92554e07, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.529496, Final residual = 0.0254955, No Iterations 14 time step continuity errors : sum local = 9.93601e10, global = 2.7549e12, cumulative = 1.30333e11 DICPCG: Solving for p, Initial residual = 0.437915, Final residual = 2.59646e06, No Iterations 1001 time step continuity errors : sum local = 9.63605e14, global = 7.5182e16, cumulative = 1.30341e11 DILUPBiCG: Solving for k, Initial residual = 0.00816952, Final residual = 1.52233e06, No Iterations 2 bounding k, min: 0.0225029 max: 319157 average: 0.485327 ExecutionTime = 5485.57 s ClockTime = 5603 s Time = 6e09 Courant Number mean: 1.71111e06 max: 0.126074 DILUPBiCG: Solving for Ux, Initial residual = 0.000781675, Final residual = 3.83783e08, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.0120933, Final residual = 5.34931e07, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.0111096, Final residual = 9.24812e07, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.69218, Final residual = 0.0339392, No Iterations 13 time step continuity errors : sum local = 1.60353e09, global = 7.22851e13, cumulative = 1.3757e11 DICPCG: Solving for p, Initial residual = 0.684253, Final residual = 4.81795e06, No Iterations 1001 time step continuity errors : sum local = 2.54543e13, global = 1.56233e15, cumulative = 1.37585e11 DILUPBiCG: Solving for k, Initial residual = 0.0184975, Final residual = 2.59942e07, No Iterations 3 bounding k, min: 6.54063 max: 1.7305e+07 average: 23.2609 ExecutionTime = 6482.16 s ClockTime = 6622 s Courant Number mean: 1.72122e06 max: 0.459119 DILUPBiCG: Solving for Ux, Initial residual = 0.000715073, Final residual = 1.65796e07, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 0.0119322, Final residual = 2.4439e06, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.00885946, Final residual = 2.70834e06, No Iterations 1 DICPCG: Solving for p, Initial residual = 0.853146, Final residual = 1.58563e05, No Iterations 1001 time step continuity errors : sum local = 1.40116e12, global = 2.70491e15, cumulative = 1.37612e11 DICPCG: Solving for p, Initial residual = 0.842176, Final residual = 2.86118e05, No Iterations 1001 time step continuity errors : sum local = 3.19604e12, global = 1.44523e14, cumulative = 1.37757e11 DILUPBiCG: Solving for k, Initial residual = 0.0295402, Final residual = 6.14694e06, No Iterations 4 bounding k, min: 1010.36 max: 1.49557e+09 average: 1435.33 ExecutionTime = 7972.54 s ClockTime = 8197 s  THIS IS THE ERROR MESSAGE  #0 Foam::error:rintStack(Foam::Ostream&) in "/gt/home/h112812/OpenFOAM/ OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #1 Foam::sigFpe::sigHandler(int) in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1 .1/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #2 __restore_rt at sigaction.c:0 [3] #3 Foam:ILUPreconditioner::calcReciprocalD(Foam::Fi eld<double>&, Foam::ld uMatrix const&) in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64Gc cDPOpt/lib/libOpenFOAM.so" [3] #4 Foam:ILUPreconditioner:ILUPreconditioner(Foam: :lduMatrix::solver con st&, Foam::dictionary const&) in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platf orms/linux64GccDPOpt/lib/libOpenFOAM.so" [3] #5 Foam::lduMatrix:reconditioner::addasymMatrixCons tructorToTable<Foam: ILUPreconditioner>::New(Foam::lduMatrix::solver const&, Foam::dictionary const&) in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libO penFOAM.so" [3] #6 Foam::lduMatrix:reconditioner::New(Foam::lduMatr ix::solver const&, Foa m::dictionary const&) in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/lin ux64GccDPOpt/lib/libOpenFOAM.so" [3] #7 Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, uns igned char) const in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64 GccDPOpt/lib/libOpenFOAM.so" [3] #8 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/gt/home/h112 812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [3] #9 Foam::fvMatrix<double>::solve() in "/gt/home/h112812/OpenFOAM/OpenFOAM2 .1.1/platforms/linux64GccDPOpt/bin/pisoFoam" [3] #10 Foam::incompressible::LESModels:neEqEddy::correc t(Foam::tmp<Foam::Geo metricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) i n "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libinc ompressibleLESModels.so" [3] #11 Foam::incompressible::LESModel::correct() in "/gt/home/h112812/OpenFOAM /OpenFOAM2.1.1/platforms/linux64GccDPOpt/lib/libincompressibleLESModels.so" [3] #12 main in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64GccD POpt/bin/pisoFoam" [3] #13 __libc_start_main in "/lib64/tls/libc.so.6" [3] #14 _start in "/gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux64Gc cDPOpt/bin/pisoFoam" [cubad10269:07980] *** Process received signal *** [cubad10269:07980] Signal: Floating point exception (8) [cubad10269:07980] Signal code: (6) [cubad10269:07980] Failing at address: 0x11c6c00001f2c [cubad10269:07980] [ 0] /lib64/tls/libc.so.6 [0x384322e2f0] [cubad10269:07980] [ 1] /lib64/tls/libc.so.6(gsignal+0x3d) [0x384322e25d] [cubad10269:07980] [ 2] /lib64/tls/libc.so.6 [0x384322e2f0] [cubad10269:07980] [ 3] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditioner15calcR eciprocalDERNS_ 5FieldIdEERKNS_9lduMatrixE+0xb5) [0x2a97b8eff5] [cubad10269:07980] [ 4] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam18DILUPreconditionerC1ERKNS _9lduMatrix6sol verERKNS_10dictionaryE+0x159) [0x2a97b8f819] [cubad10269:07980] [ 5] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam9lduMatrix14preconditioner3 1addasymMatrixC onstructorToTableINS_18DILUPreconditionerEE3NewERK NS0_6solverERKNS_10dictionaryE +0x3c) [0x2a97b8f92c] [cubad10269:07980] [ 6] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libOpenFOAM.so(_ZN4Foam9lduMatrix14preconditioner3 NewERKNS0_6solv erERKNS_10dictionaryE+0x7ef) [0x2a97b7ebaf] [cubad10269:07980] [ 7] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam5PBiCG5solveERNS_5FieldIdE ERKS2_h+0x6e3) [0x2a97b83853] [cubad10269:07980] [ 8] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE5solveERKNS _10dictionaryE+ 0x157) [0x2a969c3df7] [cubad10269:07980] [ 9] pisoFoam(_ZN4Foam8fvMatrixIdE5solveEv+0xd3) [0x421a13] [cubad10269:07980] [10] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libincompressibleLESModels.so(_ZN4Foam14incompress ible9LESModels9 oneEqEddy7correctERKNS_3tmpINS_14GeometricFieldINS _6TensorIdEENS_12fvPatchFieldE NS_7volMeshEEEEE+0x743) [0x2a95d9c443] [cubad10269:07980] [11] /gt/home/h112812/OpenFOAM/OpenFOAM2.1.1/platforms/linux 64GccDPOpt/lib/libincompressibleLESModels.so(_ZN4Foam14incompress ible8LESModel7c orrectEv+0x35) [0x2a95d10275] [cubad10269:07980] [12] pisoFoam [0x418ed4] [cubad10269:07980] [13] /lib64/tls/libc.so.6(__libc_start_main+0xdb) [0x384321c3 fb] [cubad10269:07980] [14] pisoFoam [0x41548a] [cubad10269:07980] *** End of error message ***  Could some please tell me where i went wrong? 1) Is it because of nonorthogonality of the mesh? 2) Domain size is 1.8*0.4*0.4 cubic meters. and it has 6.5 million cells. Is it too refined to use wall functions near the wall? 3) Should i change discretization schemes?  ddtSchemes { default backward; } gradSchemes { default Gauss linear; grad(p) Gauss linear; grad(U) Gauss linear; } divSchemes { default none; div(phi,U) Gauss filteredLinear; div(phi,k) Gauss limitedLinear 1; div(phi,B) Gauss limitedLinear 1; div(phi,nuTilda) Gauss limitedLinear 1; div(B) Gauss linear; div((nuEff*dev(T(grad(U))))) Gauss linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss linear corrected; laplacian((1A(U)),p) Gauss linear corrected; laplacian(DkEff,k) Gauss linear corrected; laplacian(DBEff,B) Gauss linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss linear corrected; } interpolationSchemes { default linear; interpolate(U) linear; } snGradSchemes { default corrected; } fluxRequired { default no; p ; }  4) Should i change fvSolution? I am using DICPCG for p and DILUPBiCG for U and k with  PISO { nCorrectors 2; nNonOrthogonalCorrectors 0; pRefCell 0; pRefValue 0; }  Please let me know if you need more details to figure out the problem. Kind regards, Achinta Last edited by achinta; June 19, 2012 at 09:22. Reason: improvement 

June 18, 2012, 17:15 

#197 
Member
Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 7 
Well, two hints that subsequetly come to my mind when I read your post.
Firstly, try to set nNonOrthogonalCorrectors higher. Your value of 0 won't correct nonOrthogonality... Try to go with 2 since your mesh with about 6.5 million cells will take a lot more with increasing the number of correctors. Secondly, you should initialise k with a higher value (the internalFields). That might help. Use GAMG for p. (There are a lot of tutorials on how to change the fvSolution in order to use GAMG. 

June 19, 2012, 08:13 

#198 
Member
achinta
Join Date: May 2010
Location: Sydney
Posts: 66
Rep Power: 8 
Hi Rob,
Thanks for the reply. 1) I tried different value for k (updated in the previous thread). 2) I used nNonOrthogonalCorrectors=1 since nNonOrthogonalCorrectors=2 takes lot of time. 3) I have read in some posts that GAMG is good for serial computing and PCG is good for parallel computing So i continued with PCG only.(please check http://www.cfdonline.com/Forums/ope...ureissue.html ) Please tell me your opinion. Is there a great difference between GAMG and PCG results? With the above corrections, the solution is diverging again. Do you think nNonOrthogonalCorrectors=2/GAMG can stabilize the results? I kindly request other users opinion also . Kind regards, Achinta Last edited by achinta; June 19, 2012 at 09:26. 

June 19, 2012, 11:38 

#199 
Member
Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 7 
I doubt that GAMG will stabilise your solution.
Maybe you should send someone your case so he or she can check some possibilities... Or you should at least provide some more information By the way, I usually use "Euler" instead of "backward" for ddtSchemes. 

June 20, 2012, 11:57 

#200 
New Member
João Duarte Miranda
Join Date: Jan 2012
Posts: 5
Rep Power: 7 
Does anyone have the compressible LES lowReOneEqEddy reference? I can't find it anywhere?
Thanks a lot for your help. João. 

Thread Tools  
Display Modes  

