化学種&反応

化学種&反応

DEFINE VR RATE()とDEFINE NET REACTION RATE()の使い分け

    • FAQFAQ
      Participant

      化学種輸送モデル(Species Transport)を利用する場合、UDFにより反応速度を定義することが可能です。 利用できるマクロは下記になります。 DEFINE_VR_RATEとDEFINE_NET_REACTION 使い分けは以下に示します。 (1) 化学種輸送モデル: Finite-Rate/No TCI (有限反応速度) Finite-Rate/Eddy-Dissipation (有限反応速度/渦消散モデル) Eddy-Dissipation (渦消散モデル) 反応速度を定義するマクロ: DEFINE_VR_RATE 各反応式の体積反応速度を定義します。 (The MOLAR(per reaction) reaction rate [kgmol/m3-s]) (2) 化学種輸送モデル: Finite-Rate/No TCI (Stiff Chemistry Solver) Eddy-Dissipation Concept 反応速度を定義するマクロ: DEFINE_NET_REACTION_RATE 各化学種の正味反応速度(生成、消滅速度)を定義します。(The NET (per species) reaction rate[kgmol/m3-s]) 以下に同様な反応式を利用し、各マクロを使用するサンプルを示します。 反応: CH4+H2O=CO+3H2 (1) CO+H2O=CO2+H2 (2) 関与する化学種: CH4, H2O, CO, CO2, H2 /* DEFINE_VR_RATEを利用し各反応式の反応速度を定義する */ DEFINE_VR_RATE(user_rate, c, t, r, wk, yk, rr_l, rr_t) { real rate1, rate2; /* reaction rate of reaction 1, 2 */ if (!strcmp(r->name, “reaction-1”)) { *rr_l = rate1; /* 層流反応速度 */ *rr_t = *rr_l; /* 乱流反応速度 */ } else if (!strcmp(r->name, “reaction-2”)) { *rr_l = rate2; *rr_t = *rr_l; } } /* DEFINE_NET_REACTONを利用し各化学種の正味の反応速度を定義する */ DEFINE_NET_REACTION_RATE(my_net_rr, c,t,particle, pressure, temp, yi, rr, jac) { int i; real rate1, rate2; /* reaction rate of reaction1, 2 */ rr[ch4index] = -(rate1); rr[h2oindex] = -(rate1)-(rate2); rr[coindex] = (rate1)-(rate2); rr[h2index] = 3.0*(rate1)+(rate2); rr[co2index] = (rate2); }