Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 426
    • Issues 426
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 8
    • Merge requests 8
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #1009
Closed
Open
Issue created Sep 17, 2018 by Admin@OpenFOAM-adminMaintainer

kOmegaSSTDES model converges to unphysical result

I was using kOmegaSSTDES turbulence model to run a simulation at deltaT = 0.000002. The max Courant number is less than 0.01, but the result is clearly unphysical. I then decreased deltaT to 0.000001. The result is much better but the time cost is increased. I was looking for the reason so I changed the turbulence model to dynamicKEqn with deltaT = 0.000002. the converged result is also physical. It seems that there is a problem with the kOmegaSSTDES turbulence model, which will result in unphysical result unless we limit the deltaT less than 0.000001. I am wondering if there is something wrong with my case. My case files have been attached.SST_2D3.zip

There are three sets of result figures in my file. Result 1 corresponds to kOmegaSSTDES with deltaT = 0.000002, which is the unphysical result. result1_presult1_U Result 2 is kOmegaSSTDES with deltaT = 0.000001. result2_presult2_U Result 3 is dynamicKEqn with deltaT = 0.000002. result3_presult3_U

To reproduce result 1, you can use the following command decomposePar mpirun -np 8 pimpleFoam -parallel paraFoam -builtin

To reproduce result 2, please change the deltaT in controlDict to 0.000001, then decomposePar -force mpirun -np 8 pimpleFoam -parallel paraFoam -builtin

To reproduce result 3, please change the deltaT in controlDict to 0.000002, and change the turbulenceProperties to dynamicKEqn, then decomposePar -force mpirun -np 8 pimpleFoam -parallel paraFoam -builtin

Assignee
Assign to
Time tracking