I did last year some simulations of a jet in cross-flow using the version 5.1.6 :
I try this year to run the same case in 5.3.0 (the switch of version is motivated by the use of the HTLES model developed in the thesis of V. Duffal).
I expected to find the same results, but it was not the case . Using the bdiff tool, i found a maximum difference about 10 m/s (for a maximum velocity in the domain around 20 m/s) located just at the outlet of the jet (Uj on the picture).
I try to find the origin of the discrepancy, without success. The latest (and simpler) setup I have tried is :
- k-omega SST model using the Neumann condition for omega (I have modified the clptur.f90 file for the v530 to be sure that this condition was chosen by the code)
- Spatial schemes : upwind for all variables
- Linear solvers : Conjugate Gradient with Multigrid preconditionner for Pressure and Symetric Gauss-Seidel (without preconditionner) for the other variables
The two setup as well as the mesh used are available here (the coarser mesh size is 18Mo...) : mesh and cases
If you could have a look to this problem, it would be wonderful
(Please note that the incoming profile are defined into two .dat files stored in the DATA folder. The interpolation at the boundary of the domain is done using a Python3 script also in the DATA folder. You should change the system command line into cs_user_boundary_conditions.f90 to match with your installation of Python).
Thanks in advance !
Regards
Franck