|
[Sponsors] |
July 20, 2017, 11:28 |
CMU value for k-epsilon model
|
#1 |
Member
Vedamt Chittlangia
Join Date: Feb 2016
Posts: 64
Rep Power: 9 |
Hello,
I want to define C as a local variable in FLUENT for low-Reynolds Number k-epsilon model. I am using Abe-Kondoh-Nagano model. To proceed further with full implementation I am testing only two values. I am able to get the exact results from my UDF from C = 0.09. But when I change the value of C to 0.05 in the UDF then these results don't match to results that are obtained by changing C = 0.05 from the GUI panel. Can anyone please help me resolve this issue? I am attaching my UDF below. Code:
#include "udf.h" #include "mem.h" #include "math.h" #define C_MU 0.05 #define SIG_K 1.4 #define SIG_D 1.4 /* user defined Turbulent viscosity */ DEFINE_TURBULENT_VISCOSITY(turb_vis,c,t) { real c_k, c_d, ReT, ReD, muT, fmu; c_k = C_K(c,t); c_d = C_D(c,t); /*flooring the values */ if (c_k < 1.e-10) c_k = 1.e-10; if (c_d < 1.e-10) c_d = 1.e-10; ReT = C_R(c,t)*SQR(c_k)/(C_MU_L(c,t)*c_d); ReD = C_WALL_DIST(c,t)*pow(C_MU_L(c,t)*c_d/C_R(c,t),0.25)*C_R(c,t)/C_MU_L(c,t); fmu = SQR( 1. - exp( -ReD/14.0 )) * ( 1. + 5.0*exp( -SQR(ReT/200.0 ))/(pow(ReT,0.75)) ); /* damping function */ muT = C_R(c,t)*C_MU*fmu*SQR(c_k)/c_d; if (muT > 1.e5) muT = 1.e5; return muT; } DEFINE_PRANDTL_K(tke_PR,c,t) { return SIG_K; } DEFINE_PRANDTL_D(tdr_PR,c,t) { return SIG_D; } /* AKN turbulence model: viscosity function */ |
|
July 21, 2017, 17:50 |
|
#2 |
Senior Member
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,752
Rep Power: 66 |
Is this maybe a wall function or wall modelling problem? The Cmu in the GUI panel I think only affects the Cmu for the transport equation but not the near wall stuff.
|
|
July 21, 2017, 18:09 |
|
#3 |
Member
Vedamt Chittlangia
Join Date: Feb 2016
Posts: 64
Rep Power: 9 |
I am using low-Reynolds model by Abe-Kondoh-Nagano. There are no wall functions for this model.
|
|
December 26, 2020, 12:15 |
|
#4 |
New Member
Ashhar
Join Date: Dec 2020
Posts: 1
Rep Power: 0 |
I want to change the C2_Epsilon Model constant of low Reynolds number Abe-Kondoh-Nagano model to a function using UDF functions in Fluent. I am relatively new to udf funcitons. Can anyone guide me how to change a constant model constant to a function in ANSYS Fluent.
|
|
Tags |
k-epsilon model constants, turbulence, udf and programming, viscosity |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wrong multiphase flow at rotating interface | Sanyo | CFX | 14 | February 7, 2017 18:19 |
[swak4Foam] groovyBC issue - k and epsilon | sagnikmazumdar | OpenFOAM Community Contributions | 24 | March 1, 2015 08:16 |
Jump in epsilon values near the wall : low re k-epsilon model | malaboss | OpenFOAM Verification & Validation | 1 | February 1, 2013 17:36 |
Help regarding K Epsilon model | go2saqi | Main CFD Forum | 2 | January 23, 2012 12:59 |
Questions about model constant of Cmu in Fluent | Pingzhi Fang | FLUENT | 0 | January 21, 2008 08:24 |