
[Sponsors] 
February 7, 2007, 13:19 
power law code

#1 
Guest
Posts: n/a

hi,how can take the fluent code for power law?(i know that fluent simulate the power law fluids)


February 7, 2007, 15:21 
Re: power law code

#2 
Guest
Posts: n/a

yes, you can configure it in the Define>Materials panel. You can select "nonnewtonian powerlaw" instead of "constant" for viscosity properties. When you select Edit, you can enter the consistency, the powerlaw index, the maximum and minimum value for viscosity. Note that Fluent enables the energy equation. If you are not interested in the solution of the energy equation, you can deselect "Energy" in the Solve>Controls>Solution panel and Fluent won't solve this equation. I hope I answered to your question.
Antoine 

February 8, 2007, 11:08 
Re: power law code

#3 
Guest
Posts: n/a

thanks Antoine, power law in fluent defined as : k exp(T0/T),but i need to change the power law define to : k exp(T0/(TT1)),then i need the know how to change the source code of power law fluid in fluent,if can you help me PLEASE help me.
majid reza 

February 12, 2007, 20:04 
Re: power law code

#4 
Guest
Posts: n/a

Sorry to give an answer so late; I'm interested in customizing laws for viscosity too, and I tried with your case without succeeding very well.
So, you can modify the law for viscosity in Fluent using an UDF DEFINE_PROFILE. This UDF is called in the "Material panel" in Fluent selecting "userdefined" for viscosity. In this UDF, you will define your variable for viscosity (ex: mu_app) and define a law for this variable. Then with the line "return mu_app", Fluent will consider the calculated value in the UDF as a value for the viscosity. When you define the law for viscosity, you have to access the shear rate variable and the temperature variable calculated by Fluent. This can be done using C_STRAIN_RATE_MAG(c,t) and C_T(c,t). Now I will give you the UDF I've done. It works but the problem is that the userdefined viscosity is not initialized correctly: when you select in Fluent the classic non newtonian power law, apparent viscosity is set to the value of k at initialisation. But when you initialize Fluent with your UDF (defining mu_app=k), apparent viscosity is set to the max value (because shear_rate=0), even if you hook a DEFINE_INIT macro. So results calculated with this method are not correct. I'm sorry to have written so much but if you find the way to initialize correctly this function, please tell me. Here is the UDF: #include "udf.h" real mu_app; real k; DEFINE_INIT(init_mu_app,d) { cell_t c; Thread *t; thread_loop_c(t,d) { begin_c_loop_all(c,t) { k=1; mu_app=k; } end_c_loop_all(c,t) } } DEFINE_PROPERTY(user_power_law,c,t) { real visc_dyn_min=0.0001; real visc_dyn_max=10000; real index_n=0.9; real T0=0; real T1=0; real temp=C_T(c,t); real shear_rate=C_STRAIN_RATE_MAG(c,t); real mu_temp; mu_temp=(k)*pow(shear_rate,index_n1)*exp(T0/(tempT1)); /*!!1/pow()is correct only if index_n < 1; this is written like this only to set mu_app to visc_dyn_max when shear_rate is equal to 0*/ if ((mu_temp>=visc_dyn_min)&&(mu_temp<=visc_dyn_max)) { mu_app=mu_temp; } else if (mu_temp>=visc_dyn_max) { mu_app=visc_dyn_max; } else if (mu_temp<=visc_dyn_min) { mu_app=visc_dyn_min; } return mu_app; } 

February 13, 2007, 17:28 
Re: power law code

#5 
Guest
Posts: n/a

Hi Antonie, Thank you very very much, i try to find it and then tell you again,thank you
mahjoob 

February 17, 2007, 19:42 
Re: power law code

#6 
Guest
Posts: n/a

I've finally found the causes of the problem. To initialize correctly the apparent viscosity, I used the current iteration number macro N_ITER, so when it is < 1 (computation is not launched) apparent viscosity is set to a fixed value. Also, in my previous udf, I defined k because I'm used to see a minus sign before k in the definition of the shear stress. But the apparent viscosity is positive so the minus sign is a source of error because mu_app always took a negative value < visc_dyn_min... I've also modified some lines of the last version of the udf because they finally weren't necessary. Now the udf works fine, I tested the case of a simple powerlaw and I obtained errors below 0.5% compared to the fluent integrated powerlaw. Here is the tested last version of the user defined powerlaw. Even if I don't find problems anymore, you should test it before launching important computations.
#include "udf.h" DEFINE_PROPERTY(user_power_law,c,t) { real visc_dyn_min=0.0001; /*Define the minimum limit for the value of viscosity*/ real visc_dyn_max=10000; /*Define the maximum limit for the value of viscosity*/ real index_n=0.9; /*Define the power law index*/ real k=1; /*Define the consistency index*/ real T0=0; /*Define the reference temperature T0*/ real T1=0; /*Define the fixed temperature T1*/ real temp=C_T(c,t); /*Define the access to the temperature calculated by Fluent through the temp variable*/ real mu_app; /*Define the variable for the apparent viscosity which will be returned to the solver*/ real mu_temp; /*Define a temporary variable to stock the user defined apparent viscosity*/ mu_temp=k*pow(C_STRAIN_RATE_MAG(c,t),index_n1)*exp(T0/(tempT1)); if (N_ITER<1) { mu_app=k; /*Initialize mu_app to the value of k*/ } else if ((N_ITER>=1)&& (mu_temp>visc_dyn_min)&& (mu_temp< visc_dyn_max)) } mu_app=mu_temp; /*Replace the value of mu_app by the calculated userdefined temporary viscosity when the computation has been launched*/ } else if ((N_ITER>=1)&&(mu_temp>=visc_dyn_max)) { mu_app=visc_dyn_max; } else if ((N_ITER>=1)&&(mu_temp<=visc_dyn_min)) { mu_app=visc_dyn_min; } return mu_app; /*Set the value of the userdefined viscosity in Fluent*/ } Best regards Antoine 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
NonNewtonian Power Law for Viscosity  John  FLUENT  6  November 26, 2014 09:44 
DPMSource Code for Standard Vaporization Law  Woo Meng Wai  FLUENT  6  August 15, 2013 08:10 
How to implement power law or hybrid schemes??  kamkari  OpenFOAM  10  May 9, 2011 10:29 
power law inlet condition  asharma  OpenFOAM  8  December 3, 2010 07:32 
new parameter for non newtonian power law model!  Antoine  FLUENT  0  July 21, 2008 04:30 