
[Sponsors] 
Initial and Final Residual of omega 0 by calculation with kw turbulence model 

LinkBack  Thread Tools  Search this Thread  Display Modes 
August 9, 2019, 17:03 
Initial and Final Residual of omega 0 by calculation with kw turbulence model

#1 
New Member
Join Date: May 2015
Posts: 16
Rep Power: 6 
Hi,
i'm trying to simulate an airfoil of an linear cascade. The mesh contains one airfoil and has periodic boundaries near the pressure and suction side. It is a q3dcase with only one layer, with a slip BC on the one sidewall and a symmetry BC on the other. The thermophysical properties are Perfect gas with constant thermodynamics and Sutherland for transport. Check mesh only complains about about high aspect ratio cells. But these are necessary for y+ values less then 1. I'm using the kw turbulence model and the Smooth Solver for p, U, h, k and omega. The divergence scheme is Lust. As for the boundary conditions i choosed k = Low Re Wall Function and omega = Standard Wall Function (nu_t = Low Re Wall Function and alpha_T Standard) on the walls. The simulation runs fine and the residuals went beyond 1e5 for Ux, Uy, Uz, p and k (h get stuck by ~1e4). But here is the problem: The solver won't solve omega from the first timestep. Code:
velocityDampingConstraint velocityDamping damped 0 (0%) of cells smoothSolver: Solving for Ux, Initial residual = 3.597191e007, Final residual = 6.192554e008, No Iterations 2 smoothSolver: Solving for Uy, Initial residual = 2.276461e007, Final residual = 4.008005e008, No Iterations 2 smoothSolver: Solving for Uz, Initial residual = 1.496093e005, Final residual = 5.666532e008, No Iterations 8 smoothSolver: Solving for h, Initial residual = 8.749293e005, Final residual = 1.187253e008, No Iterations 8 smoothSolver: Solving for p, Initial residual = 8.486576e007, Final residual = 5.21028e007, No Iterations 102 smoothSolver: Solving for p, Initial residual = 6.523312e007, Final residual = 4.989128e007, No Iterations 102 smoothSolver: Solving for p, Initial residual = 5.678153e007, Final residual = 4.820913e007, No Iterations 102 time step continuity errors : sum local = 7.096876e012, global = 6.886245e012, cumulative = 0.0002020726 rho max/min : 1.209845 1.162966 smoothSolver: Solving for omega, Initial residual = 0, Final residual = 0, No Iterations 0 smoothSolver: Solving for k, Initial residual = 4.188657e007, Final residual = 2.28007e008, No Iterations 2 ExecutionTime = 934.098 s ClockTime = 934 s HTML Code:
https://www.cfdonline.com/Forums/openfoamsolving/81527problemsstmodelstrangebehaviour.html Does anyone had to deal with the same error and found a solution? And what does it mean for the results when one field variable isn't solved? Greetings, Stuntmanbob Edit: I'm working on Windows10 with the free Version of SimFlow 3.1, which uses the OpenFoam 16.12. I also tried the same simulation with OpenFoam 7 (19.06) on ubuntu installed via the windows subsystem for linux. But the result is the same Last edited by Stuntmanbob; August 9, 2019 at 17:10. Reason: Additional information 

August 11, 2019, 09:10 

#2  
Member
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 82
Rep Power: 5 
Quote:
I've never used komega, but used kepsilon quite a lot, which I think is pretty similar to komega in terms of how things are solved. 

August 16, 2019, 15:30 

#3 
New Member
Join Date: May 2015
Posts: 16
Rep Power: 6 
Hi,
thanks for the reply ich checked the omega and k distribution (see pics in the attachment). The omega values near the surface of the blade and in the wake are lower than in the main flow. The same with the k distribution. That seems fine to me. Another point is, that the min, max and average values of omega change during the simulation: Code:
volAverage() of omega = 2645.706 min(omega) = 1.153726 at location (0.2775581 0.01820756 0.0694632) max(omega) = 3.789602e+008 at location (0.1311357 0.2398932 0.0694632) Greetings, Stuntmanbob 

Today, 05:02 

#4  
Member
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 82
Rep Power: 5 
Okay, so omega does have a nonuniform distribution, which would mean that the omega equation is actually solved. Not sure what's going on. One thing that is similar to this is when you use the compressible solvers, for density, the log file always read
Quote:
In the kepsilon model, epsilon is the turbulence dissipation, so its maximum value could be very high. I've had something in the order of e+006, so I guess your max(omega) should be reasonable? Need some kw experts here! 

Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Segmentation fault when using reactingFOAM for Fluids  Tommy Floessner  OpenFOAM Running, Solving & CFD  4  April 22, 2018 12:30 
chtMultiRegionSimpleFoam turbulent case  Aditya Patil  OpenFOAM Running, Solving & CFD  6  April 24, 2017 22:13 
simpleFoam error  "Floating point exception"  mbcx4jc2  OpenFOAM Running, Solving & CFD  12  August 4, 2015 02:20 
pimpleFoam: turbulence>correct(); is not executed when using residualControl  hfs  OpenFOAM Running, Solving & CFD  3  October 29, 2013 08:35 
calculation stops after few time steps  sivakumar  OpenFOAM Running, Solving & CFD  7  March 17, 2013 06:37 