

July 9, 2024 at 8:32 amæ´‹ äºŽSubscriber
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
DEFINE_SOURCE(uds_source_W, cell, sec_th, dS, eqn){real P = C_P(cell, sec_th)+101325;real T = C_T(cell, sec_th);real Xw = C_UDSI(cell, sec_th, 0);real Xoh = C_UDSI(cell, sec_th, 1);real P_bar = P / 100000.0;Âreal source_W;Thread*mix_th;Thread*pri_th;mix_th = THREAD_SUPER_THREAD(sec_th);pri_th = THREAD_SUB_THREAD(mix_th, 0);sec_th = THREAD_SUB_THREAD(mix_th,1);source_W =K_X*C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aw * exp(Ew / (R * T )) * (1  Xw);dS[eqn] =K_X*C_R(cell, sec_th) * C_VOF(cell, sec_th) * Aw * exp(Ew / (R * T ));return source_W;}ÂDEFINE_SOURCE(uds_source_OH, cell, sec_th, dS, eqn){real P = C_P(cell,sec_th)+101325;real T = C_T(cell, sec_th);real Xw = C_UDSI(cell, sec_th, 0);real Xoh = C_UDSI(cell, sec_th, 1);real P_bar = P / 100000.0;real source_OH;Thread*mix_th;Thread* pri_th;Â Â mix_th = THREAD_SUPER_THREAD(sec_th);pri_th = THREAD_SUB_THREAD(mix_th, 0);sec_th = THREAD_SUB_THREAD(mix_th,1);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);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);return source_OH;}ÂDEFINE_SOURCE(uds_source_Qe, cell, thread, dS, eqn){real P = C_P(cell, thread)+101325 ;real T = C_T(cell, thread);real Xw = C_UDSI(cell, thread, 0);real Xoh = C_UDSI(cell, thread, 1);real P_bar = P / 100000.0;real Hoh = 7.705e4;real Hw = 8.6e4;real dXoh_dt = Aoh * exp(Eoh / (R * T)) * coh0 * (1  Xoh) * (cnco0 / coh0  2 * (cw0 / coh0) * Xw  Xoh) ;real dXw_dt = Aw * exp(Ew / (R * T )) * (1  Xw) ;real source_Qe;Thread* pri_th;Thread* sec_th;pri_th = THREAD_SUB_THREAD(thread, 0);sec_th = THREAD_SUB_THREAD(thread,1);source_Qe = K_X * C_R(cell,Â sec_th) * C_VOF(cell, sec_th) * (Hoh * coh0 * dXoh_dt + Hw * cw0 * dXw_dt) / Rpu;dS[eqn] = 0;return source_Qe;}Â

July 9, 2024 at 10:46 amRobForum Moderator
Check all the variables are defined. Then look at the values you're using and see what's going on in the model. Note, from the source names are these the values of the source, or a UDS?

Viewing 1 reply thread
 You must be logged in to reply to this topic.
Ansys Innovation Space
Trending discussions
 Fluent crash on starting
 Transient Porous Media Combustion
 VOF Free Surface Deformation
 Melting and solidification process for different boundary condition& time period
 Kinetic Mechanisms and Arrhenius Equation Parameters in ANSYS Fluent
 How do you approach this?
 Fluent 20.1 – bug in writing a pathlines file by GUI
 Connection Stiffness Ansys AQWA
 Convective Augmentation Factor in DEFINE_HEAT_FLUX
 Reg: Submission script to submit job on HPC
Top Contributors

341

160

130

117

99
Top Rated Tags
Â© 2024 Copyright ANSYS, Inc. All rights reserved.
Ansys does not support the usage of unauthorized Ansys software. Please visit www.ansys.com to obtain an official distribution.