{"id":462273,"date":"2026-08-05T14:57:38","date_gmt":"2026-08-05T14:57:38","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/manual-implementation-pressure-strain-model\/"},"modified":"2026-08-05T14:57:38","modified_gmt":"2026-08-05T14:57:38","slug":"manual-implementation-pressure-strain-model","status":"publish","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/manual-implementation-pressure-strain-model\/","title":{"rendered":"Manual implementation pressure-strain model"},"content":{"rendered":"

Hello,<\/p>\n

When a Reynolds stress turbulence model is selected, Fluent allows to specify which scale variable (epsilon or omega) and which formulation for the pressure-strain source term (LRR or SSG) the user wants to use. <\/p>\n

The LRR model is made of four terms which depend on the coefficients C_1 and C_2. The first affects the slow term and the second affects the remaining 3 terms (rapid terms). See Theory Guide, equation 4.247. No wall reflection is needed when BSL-omega is the scale-trasporting equation (like in my case). <\/p>\n

I’ve been trying to reimplement the rapid terms manually through UDFs:<\/p>\n

    \n
  1. the slow term remains untouched.<\/li>\n
  2. the coefficient C_2 is changed so that the term with gamma_0 goes to 0.<\/li>\n
  3. the remaining terms (with alpha_0 and beta_0) are canceled out by a UDF containing six DEFINE_SOURCE functions, one for each of the components of the turbulent stress tensor. This UDF works as a sink.<\/li>\n
  4. a second UDF, containing six DEFINE_SOURCE functions, reinjects the rapid terms. <\/li>\n<\/ol>\n

    The idea is that simulations with this model built through UDFs should give the same result as the simulation using the standard built-in LRR term. I found this to be somehow true after a few thousands of iterations. After that, if the simulation keeps running, the solution begins to deteriorate and the residuals steadily rise. The domain is 2D, axisymmetric swirl. The deterioration seems to begin in the axis region and then propagates to the rest of the domain.<\/p>\n

    My question is whether Fluent imposes limiters and\/or realizability checks that are not reported in the guides.<\/p>\n

    Thank you.<\/p>\n","protected":false},"template":"","class_list":["post-462273","topic","type-topic","status-publish","hentry","topic-tag-udf-2","topic-tag-reynolds-stress-model","topic-tag-turbulence-models-1"],"aioseo_notices":[],"aioseo_head":"\n\t\t\n\t\n\t\n\t\n\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t\n\t\t