{"id":311896,"date":"2023-10-17T17:56:04","date_gmt":"2023-10-17T17:56:04","guid":{"rendered":"\/forum\/forums\/topic\/udf-for-interstitial-heat-transfer-coefficient-in-non-equilibrium-porous-media\/"},"modified":"2023-10-17T17:56:04","modified_gmt":"2023-10-17T17:56:04","slug":"udf-for-interstitial-heat-transfer-coefficient-in-non-equilibrium-porous-media","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/udf-for-interstitial-heat-transfer-coefficient-in-non-equilibrium-porous-media\/","title":{"rendered":"UDF for Interstitial heat transfer coefficient in non equilibrium porous media"},"content":{"rendered":"

Hello, I am trying to create a UDF to fit a numerical model which must be applied as UDFs in ANSYS Fluent to model a phase change material in a porous medium. I have attached the mathematical model I am trying to follow as well as my current code. I am currently stuck on getting the interstitial heat transfer coefficient to come up as a UDF, that can be selected in ANSYS Fluent, for the non-equilibrium model of porous medium. Any help would be greatly appreciated. Should it be define_profile or define_property<\/p>\n

#include <\/span>“udf.h”<\/span><\/h1>\n

#include <\/span>“mem.h”<\/span><\/p>\n

#define<\/span> Tsolidus <\/span>300.0<\/span><\/p>\n

#define<\/span> Tliquidus <\/span>304.0<\/span><\/p>\n

#define<\/span> L_f <\/span>250<\/span> \/\/Latent heat of fusion<\/span><\/p>\n

 <\/p>\n

#define<\/span> UDM_RHO_F_PREV <\/span>0<\/span>  <\/span>\/\/ UDM index for storing previous rho_f<\/span><\/p>\n

#define<\/span> UDM_LAMBDA_PREV <\/span>1<\/span>  <\/span>\/\/ UDM index for storing previous lambda<\/span><\/p>\n

 <\/p>\n

DEFINE_PROFILE<\/span>(h_sfq, c, t)<\/span><\/p>\n

{<\/span><\/p>\n

    <\/span>real dl, epsilon, mu, nu, pr, rho_pcm, u, v, mu_pcm, k_pcm, Re_d, h_sfq, local_h_sfq;<\/span><\/p>\n

    <\/span><\/span>cell_t<\/span> cell;<\/span><\/p>\n

 <\/p>\n

    <\/span>dl <\/span>= <\/span>0.0194<\/span>;<\/span> \/\/Diameter Matlab<\/span><\/p>\n

    <\/span>epsilon <\/span>= <\/span>0.85<\/span>;<\/span> \/\/Porosity Matlab<\/span><\/p>\n

    <\/span>mu <\/span>= <\/span>C_MU_L<\/span>(c, t);<\/span>        <\/span>\/\/ Laminar viscosity<\/span><\/p>\n

    <\/span>rho_pcm <\/span>= <\/span>C_R<\/span>(c, t);<\/span>          <\/span>\/\/ Density<\/span><\/p>\n

    <\/span>nu <\/span>=<\/span> mu<\/span>\/<\/span>rho_pcm;<\/span>             <\/span>\/\/ Kinematic Viscosity<\/span><\/p>\n

    <\/span>pr <\/span>=<\/span> nu<\/span>\/<\/span>rng_alpha<\/span>(<\/span>1.<\/span>, mu, mu);<\/span> \/\/ Prandtl Number<\/span><\/p>\n

 <\/p>\n

    <\/span>u <\/span>= <\/span>C_U<\/span>(c, t);<\/span><\/p>\n

    <\/span>v <\/span>= <\/span>C_V<\/span>(c, t);<\/span><\/p>\n

    <\/span>mu_pcm <\/span>= <\/span>C_MU_L<\/span>(c ,t);<\/span><\/p>\n

    <\/span>k_pcm <\/span>= <\/span>C_K_L<\/span>(c, t);<\/span><\/p>\n

 <\/p>\n

    <\/span>Re_d <\/span>=<\/span> (rho_pcm <\/span>* <\/span>sqrt<\/span>(u<\/span>*<\/span>u <\/span>+<\/span> v<\/span>*<\/span>v) <\/span>*<\/span>dl) <\/span>\/<\/span> (epsilon <\/span>*<\/span> mu_pcm);<\/span><\/p>\n

 <\/p>\n

    <\/span><\/span>if<\/span> (Re_d <\/span>> <\/span>0 <\/span>&&<\/span> Re_d <\/span><= <\/span>40<\/span>)<\/span><\/p>\n

        <\/span>local_h_sfq <\/span>=<\/span> (<\/span>0.76 <\/span>* <\/span>pow<\/span>(Re_d, <\/span>0.4<\/span>) <\/span>* <\/span>pow<\/span>(pr, <\/span>0.37<\/span>) <\/span>*<\/span> k_pcm)<\/span>\/<\/span>dl;<\/span><\/p>\n

    <\/span><\/span>else if<\/span> (Re_d <\/span>> <\/span>40 <\/span>&&<\/span> Re_d <\/span><= <\/span>1000<\/span>)<\/span><\/p>\n

        <\/span>local_h_sfq <\/span>=<\/span> (<\/span>0.52 <\/span>* <\/span>pow<\/span>(Re_d, <\/span>0.5<\/span>) <\/span>* <\/span>pow<\/span>(pr, <\/span>0.37<\/span>) <\/span>*<\/span> k_pcm)<\/span>\/<\/span>dl;<\/span><\/p>\n

    <\/span><\/span>else if<\/span> (Re_d <\/span>> <\/span>1000 <\/span>&&<\/span> Re_d <\/span><= <\/span>20000<\/span>)<\/span><\/p>\n

        <\/span>local_h_sfq <\/span>=<\/span> (<\/span>0.26 <\/span>* <\/span>pow<\/span>(Re_d, <\/span>0.6<\/span>) <\/span>* <\/span>pow<\/span>(pr, <\/span>0.37<\/span>) <\/span>*<\/span> k_pcm)<\/span>\/<\/span>dl;<\/span><\/p>\n

    <\/span><\/span>else <\/span><\/p>\n

    <\/span><\/span><\/p>\n

    <\/span><\/span>return<\/span> local_h_sfq;<\/span><\/p>\n

}<\/span><\/p>\n

<\/p>\n","protected":false},"template":"","class_list":["post-311896","topic","type-topic","status-closed","hentry","topic-tag-fluent","topic-tag-pcm-charging-discharging-2","topic-tag-porosity-1","topic-tag-porous","topic-tag-udf"],"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\n\t\t\n\t\t\n\t\t