±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 2254  |  »Ø¸´: 0

ÈÕÔÂÐÇ¿Õ

Ìú³æ (³õÈëÎÄ̳)

[ÇóÖú] fluentÔËÐйý³ÌÖгöÏÖ´íÎó£ºDivergence detected in AMG solver

ÎÒ×öµÄÊÇÃÜÏàÊäÔË´²ÌáÉý¹ÜÄÚÃºÆø»¯µÄÊýֵģÄâ¡£·¢ÉúµÄ¶à¸ö»¯Ñ§·´Ó¦ÐèҪдUDF¼ÓÈë½øÈ¥£¬Ã»µ¼ÈëUDF֮ǰ·´Ó¦»¹ÄܽøÐУ¬¿ÉÊÇÒ»µ¼ÈëUDF¾Í²»ÐÐÁË¡£
ËùÒÔÎÒÅжÏÓ¦¸ÃÊÇUDFдµÄÓÐÎÊÌ⣬¿Éµ¼ÈëµÄʱºòÏÔʾÓ﷨ʲôµÄû´íÎ󣬾ßÌåÄÄÀï²»ºÏÀí£¬ÎÒ̽¾¿Á˺ܾÃҲûÕÒ³ö벡¡£»¹Ï£Íû´óÉñÃÇÄܹ»¸øÓèÖ¸µ¼¡£
fluentÌáʾµÄ´íÎóÈçÏ£º
# Divergence detected in AMG solver: mp-x-momentum -> Decreasing coarsening group size!
# Divergence detected in AMG solver: mp-x-momentum -> Increasing relaxation sweeps!
# Divergence detected in AMG solver: pressure correction -> Turning off correction scaling!
# Divergence detected in AMG solver: pressure correction -> Increasing relaxation sweeps!
Errorivergence detected in AMG solver:pressure correction

Errorivergence detected in AMG solver:pressure correction.
ÆäÖÐÒ»¸öUDFÈçÏ£º

#include "udf.h"
#include "stdio.h"
#include "time.h"
DEFINE_HET_RXN_RATE(c_arr,c,t,hr,mw,yi,rr,rr_t)
{
Domain **domain_reactant = hr->domain_reactant;
real *stoich_reactant = hr->stoich_reactant;
int *reactant = hr->reactant;
int i;
int sp_id;
int dindex;
Thread *t_reactant;
real pp;
real T,C;
real p1,p3,p4,p5,mi,msum,m0,m1,m2,m3,m4,m5,m6,m7,m8,m9,ci;
/* instead of compute rr directly, compute log(rr) and then take exp */
*rr = 0;


i=0;
sp_id = reactant; /* species ID to access mw and yi */
if (sp_id == -1) sp_id = 0; /* if phase does not have species,
mw, etc. will be stored at index 0 */
dindex = DOMAIN_INDEX(domain_reactant);
/* domain index to access mw & yi */
t_reactant = THREAD_SUB_THREAD(t,dindex);
/* get conc. */
pp = 6*C_P(c,t_reactant)/101325;
m0 = yi[dindex][0]/mw[dindex][0];
m1 = yi[dindex][1]/mw[dindex][1];
m2 = yi[dindex][2]/mw[dindex][2];
m3 = yi[dindex][3]/mw[dindex][3];
m4 = yi[dindex][4]/mw[dindex][4];
m5 = yi[dindex][5]/mw[dindex][5];
m6 = yi[dindex][6]/mw[dindex][6];
m7 = yi[dindex][7]/mw[dindex][7];
m8 = yi[dindex][8]/mw[dindex][8];
m9 = yi[dindex][9]/mw[dindex][9];
msum=m0+m1+m2+m3+m4+m5+m6+m7+m8+m9;
p3=pp*yi[0][3]/mw[0][3]/msum;
p4=pp*yi[0][4]/mw[0][4]/msum;
p5=pp*yi[0][5]/mw[0][5]/msum;
ci = MAX(p3,0);
*rr+= log(ci);
ci = MAX(p4,0);
*rr+= log(ci);
ci = MAX(p5,0);
*rr+= log(ci);


i=1;
sp_id = reactant; /* species ID to access mw and yi */
if (sp_id == -1) sp_id = 0; /* if phase does not have species,
mw, etc. will be stored at index 0 */
dindex = DOMAIN_INDEX(domain_reactant);
/* domain index to access mw & yi */
t_reactant = THREAD_SUB_THREAD(t,dindex);
/* get conc. */
ci = yi[1][5]*C_R(c,t_reactant)/mw[1][5];/*»ñµÃµÚ¶þÏàÖз´Ó¦Îï̼µÄÖÊÁ¿·ÖÊý*/
ci = MAX(ci,0);
*rr+= log(ci);
T =C_T(c,t_reactant);
if (T>1117) T=1117;
*rr+= log(2250)-42000/(8.314*T)-(17.29-16326/T);
/* 1.e-40 < rr < 1.e40 */
*rr = MAX(*rr,-40);
*rr = MIN(*rr,40);
*rr = exp(*rr);

}

Ï£Íû´óÉñ¿´µ½ºóÄܹ»¸øÓèÖ¸µ¼£¬¶àлÁË¡£
»Ø¸´´ËÂ¥
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ ÈÕÔÂÐÇ¿Õ µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] µ÷¼Á +5 ÐܶþÏëÉϰ¶ 2026-04-04 5/250 2026-04-04 13:51 by ССÊ÷2024
[¿¼ÑÐ] ¿¼Ñе÷¼Á +4 zybz³å³å³å 2026-04-03 6/300 2026-04-04 13:08 by zybz³å³å³å
[¿¼ÑÐ] 301Çóµ÷¼Á +16 ÂæÍÕÄÐÈË 2026-04-02 16/800 2026-04-04 11:59 by asdhh1991
[¿¼ÑÐ] 265Çóµ÷¼Á +20 ÁºÁºÐ£Ð£ 2026-04-01 21/1050 2026-04-04 00:38 by userper
[¿¼ÑÐ] 283·ÖÇóµ÷¼Á +6 СÄô°®Ñ§Ï° 2026-04-03 6/300 2026-04-03 20:46 by zhq0425
[¿¼ÑÐ] 303Çóµ÷¼Á +9 DLkz1314. 2026-03-30 9/450 2026-04-03 18:34 by lsÁõ˧
[¿¼ÑÐ] ѧ˶»úе¹¤³Ì303Çóµ÷¼Á +6 ÎÞÃûËùÒÔ½ÐÎâÃ÷ 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[¿¼ÑÐ] ²ÄÁÏר˶µ÷¼Á +18 Ò¬Ò¬¡£ 2026-03-29 18/900 2026-04-03 16:45 by ÁáÁá0606
[¿¼ÑÐ] 11408£¬284·Ö£¬¶þÕ½Õæ³ÏÇóµ÷¼Á +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[¿¼ÑÐ] Ò»Ö¾Ô¸±±¾©½»Í¨´óѧ²ÄÁϹ¤³Ì×Ü·Ö358 +4 cs0106 2026-04-03 4/200 2026-04-03 13:41 by °ÙÁéͯ888
[¿¼ÑÐ] ÊýÒ»Ó¢Ò»285Çóµ÷¼Á +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[¿¼ÑÐ] 0703»¯Ñ§ +7 goldtt 2026-04-02 9/450 2026-04-03 09:36 by À¶ÔÆË¼Óê
[¿¼ÑÐ] ²ÄÁϵ÷¼Á +7 Ò»ÑùYWY 2026-04-02 7/350 2026-04-02 21:49 by dongzh2009
[¿¼ÑÐ] 366Çóµ÷¼ÁÒ»Ö¾Ô¸¶«±±´óѧ +8 ÔËÆøÀ´µÃÈôÓÐËÆÎ 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[¿¼ÑÐ] 362Çóµ÷¼Á +14 Î÷ÄϽ»²ÄÁÏר˶3 2026-03-31 14/700 2026-04-02 17:50 by yunlongyang
[¿¼ÑÐ] 413Çóµ÷¼Á +3 ¿Âijij 2026-03-31 3/150 2026-04-02 16:59 by zzsw+
[¿¼ÑÐ] 318Çóµ÷¼Á +3 óÆÐÐÖÂÔ¶. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[¿¼ÑÐ] ÄÜÔ´¶¯Á¦ µ÷¼Á +3 ²»ÆÆ²»Á¢0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[¿¼ÑÐ] 354Çóµ÷¼Á +4 lxb598 2026-03-31 5/250 2026-04-02 09:55 by Jaylen.
[¿¼ÑÐ] Ò»Ö¾Ô¸Öк£Ñó320»¯Ñ§¹¤³ÌÓë¼¼Êõѧ˶Çóµ÷¼Á +8 ÅûÐÇºÓ 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û