{"id":371074,"date":"2024-07-09T08:32:44","date_gmt":"2024-07-09T08:32:44","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/source-item-udf\/"},"modified":"2024-07-09T08:32:44","modified_gmt":"2024-07-09T08:32:44","slug":"source-item-udf","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/source-item-udf\/","title":{"rendered":"Source item udf"},"content":{"rendered":"

I am doing research on foaming simulation. The udfs of the water component conversion rate Xw and hydroxyl component conversion rate Xoh source terms are as follows. Why does the rate of change become very small as the calculation proceeds and can almost be ignored? Yes, please give me some advice<\/p>\n

DEFINE_SOURCE(uds_source_W, cell, sec_th, dS, eqn)<\/div>\n
{<\/div>\n
real P = C_P(cell, sec_th)+101325;<\/div>\n
real T = C_T(cell, sec_th);<\/div>\n
real Xw = C_UDSI(cell, sec_th, 0);<\/div>\n
real Xoh = C_UDSI(cell, sec_th, 1);<\/div>\n
real P_bar = P \/ 100000.0;<\/div>\n
 <\/div>\n
real source_W;<\/div>\n
Thread*mix_th;<\/div>\n
Thread*pri_th;<\/div>\n
mix_th = THREAD_SUPER_THREAD(sec_th);<\/div>\n
pri_th = THREAD_SUB_THREAD(mix_th, 0);<\/div>\n
sec_th = THREAD_SUB_THREAD(mix_th,1);<\/div>\n
source_W =K_X*C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aw * exp(-Ew \/ (R * T )) * (1 – Xw);<\/div>\n
dS[eqn] =K_X*C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aw * exp(-Ew \/ (R * T ));<\/div>\n
return source_W;<\/div>\n
}<\/div>\n
 <\/div>\n
DEFINE_SOURCE(uds_source_OH, cell, sec_th, dS, eqn)<\/div>\n
{<\/div>\n
real P = C_P(cell,sec_th)+101325;<\/div>\n
real T = C_T(cell, sec_th);<\/div>\n
real Xw = C_UDSI(cell, sec_th, 0);<\/div>\n
real Xoh = C_UDSI(cell, sec_th, 1);<\/div>\n
real P_bar = P \/ 100000.0;<\/div>\n
real source_OH;<\/div>\n
Thread*mix_th;<\/div>\n
Thread* pri_th;<\/div>\n
    mix_th = THREAD_SUPER_THREAD(sec_th);<\/div>\n
pri_th = THREAD_SUB_THREAD(mix_th, 0);<\/div>\n
sec_th = THREAD_SUB_THREAD(mix_th,1);<\/div>\n
source_OH = K_X * C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aoh * exp(-Eoh \/ (R * T )) * coh0 * (1 – Xoh) * (cnco0 \/ coh0 – 2 * (cw0 \/ coh0) * Xw – Xoh);<\/div>\n
dS[eqn] = K_X * C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aoh * exp(-Eoh \/ (R * T )) * coh0 * ((-2 * cw0 \/ coh0) * (1 – Xw) * cnco0 \/ coh0 – 2 * Xw * cw0 \/ coh0 – Xoh);<\/div>\n
return source_OH;<\/div>\n
}<\/div>\n
 <\/div>\n
DEFINE_SOURCE(uds_source_Qe, cell, thread, dS, eqn)<\/div>\n
{<\/div>\n
real P = C_P(cell, thread)+101325 ;<\/div>\n
real T = C_T(cell, thread);<\/div>\n
real Xw = C_UDSI(cell, thread, 0);<\/div>\n
real Xoh = C_UDSI(cell, thread, 1);<\/div>\n
real P_bar = P \/ 100000.0;<\/div>\n
real Hoh = 7.705e4;<\/div>\n
real Hw = 8.6e4;<\/div>\n
real dXoh_dt = Aoh * exp(-Eoh \/ (R * T)) * coh0 * (1 – Xoh) * (cnco0 \/ coh0 – 2 * (cw0 \/ coh0) * Xw – Xoh) ;<\/div>\n
real dXw_dt = Aw * exp(-Ew \/ (R * T )) * (1 – Xw) ;<\/div>\n
real source_Qe;<\/div>\n
Thread* pri_th;<\/div>\n
Thread* sec_th;<\/div>\n
pri_th = THREAD_SUB_THREAD(thread, 0);<\/div>\n
sec_th = THREAD_SUB_THREAD(thread,1);<\/div>\n
source_Qe = K_X * C_R(cell,  sec_th) * C_VOF(cell, sec_th) * (Hoh * coh0 * dXoh_dt + Hw * cw0 * dXw_dt) \/ Rpu;<\/div>\n
dS[eqn] = 0;<\/div>\n
return source_Qe;<\/div>\n
}<\/div>\n

 <\/p>\n","protected":false},"template":"","class_list":["post-371074","topic","type-topic","status-closed","hentry"],"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