版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者QiaoPengyu將贈送您 5 個金幣

當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖

QiaoPengyu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 29
帖子: 10
在線: 2.9小時
蟲號: 28904850
注冊: 2022-03-02
專業(yè): 工程熱物理與能源利用

[求助] fluent加載udf后初始化閃退

fluent不加載udf的時候可以計算，一加入udf就會報錯閃退（udf編譯和加載都可以成功）。udf如下大佬們看看udf是不是存在問題。
#include "udf.h"
#include "flow.h"
#include"unsteady.h"
#include"stdio.h"
#define M0 0.032
#define M1 0.028
#define k10 7.87e-9
#define k11 9.86e-9
#define k20 1541.211
#define k21 2010.908
#define k30 6.79e-10
#define k31 1.67e-9
#define k40 1968.24
#define k41 2250
#define K0 62
#define K1 19.7
#define D 0.08
#define d 1.6e-3
#define h0 12000
#define h1 18000
#define L 1035
#define hf 0.2
#define hw 3.84e-3
#define Tw 288
#define e 0.4
#define cp_s 1010    /*specific heat of solid                      [cp_s]=j/(kg.K)*/
/* Define which user-defined scalars to use.*/
enum
{
  q0,q1,Q0,Q1,OM,TM,TE,T_S
};
DEFINE_ADJUST(adjust_fcn, domain)
{
Thread* t;
cell_t c;
/*solve the q0.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(q0)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(q0)))
      {
         begin_c_loop(c, t)
         {
            real p_0, p_1, x0, x1, b0, b1;
            p_0 = 7 * C_YI(c, t, 0) * ABS_P(C_P(c, t), op_pres) / (8 - C_YI(c, t, O));
            p_1 = 8 * C_YI(c, t, 1) * ABS_P(C_P(c, t), op_pres) / (7 + C_YI(c, t, 1));
            x0 = k10 * exp(k20 / C_T(c, t));
            b0 = k30 * exp(k40 / C_T(c, t));
            b1 = k31 * exp(k41 / C_T(c, t));
            C_UDSI(c, t, q0) = x0 * p_0 / (1 + b0 * p_0 + b1 * p_1);
         }
         end_c_loop(c, t)
      }
}
/* solve the q1.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(q1)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(q1)))
      {
         begin_c_loop(c, t)
         {
            real p_0, p_1, x0, x1, b0, b1;
            x1 = k11 * exp(k21 / C_T(c, t));
            b0 = k30 * exp(k40 / C_T(c, t));
            b1 = k31 * exp(k41 / C_T(c, t));
            C_UDSI(c, t, q1) = x1 * p_1 / (1 + b0 * p_0 + b1 * p_1);
         }
         end_c_loop(c, t)
      }
}
/* solve the Q0.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(Q0)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(Q0)))
      {
         begin_c_loop(c, t)
         {
            real physical_dt;
            physical_dt = RP_Get_Real("physical-time-step"

;
            C_UDSI(c, t, Q0) = C_UDSI_M1(c, t, Q0) + K0 * physical_dt * (C_UDSI(c, t, q0) - C_UDSI_M1(c, t, Q0));
         }
         end_c_loop(c, t)
      }
}
/* solve the Q1.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(Q1)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(Q1)))
      {
         begin_c_loop(c, t)
         {
            real physical_dt;
            physical_dt = RP_Get_Real("physical-time-step"

;
            C_UDSI(c, t, Q1) = C_UDSI_M1(c, t, Q1) + K1 * physical_dt * (C_UDSI(c, t, q1) - C_UDSI_M1(c, t, Q1));
         }
         end_c_loop(c, t)
      }
}
/*  solve the OM.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(OM)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(OM)))
      {
         begin_c_loop(c, t)
         {
            C_UDSI(c, t, OM) = -(1.0 - e) * M0 * L * K0 * (C_UDSI(c, t, q0) - C_UDSI_M1(c, t, Q0));
         }
         end_c_loop(c, t)
      }
}
/* solve the TM.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(TM)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(TM)))
      {
         begin_c_loop(c, t)
         {
            C_UDSI(c, t, TM) = C_UDSI(c, t, OM) - (1.0 - e) * M1 * L * K1 * (C_UDSI(c, t, q1) - C_UDSI_M1(c, t, Q1));
         }
         end_c_loop(c, t)
      }
}

/* solve the TE.*/
thread_loop_c(t, domain)
{
      if (NULL != THREAD_STORAGE(t, SV_UDS_I(TE)) && NULL != T_STORAGE_R_NV(t, SV_UDSI_M1(TE)))
      {
         begin_c_loop(c, t)
         {
            C_UDSI(c, t, TE) = (1.0 - e) * L * (K0 * h0 * (C_UDSI(c, t, q0) - C_UDSI_M1(c, t, Q0)) + K1 * h1 * (C_UDSI(c, t, q1) - C_UDSI_M1(c, t, Q1)));
         }
         end_c_loop(c, t)
      }
}
}
      DEFINE_UDS_UNSTEADY(uns_time, c, t, i, apu, su)
      {
         real  physical_dt, vol, rho, phi_old;
         physical_dt = RP_Get_Real("physical-time-step"

;
         vol = C_VOLUME(c, t);
         rho = L;
         *apu = -rho * cp_s * vol / physical_dt;             /*implicit  part*/
         phi_old = C_STORAGE_R(c, t, SV_UDSI_M1(T_S));
         *su = rho * cp_s * vol * phi_old / physical_dt;    /*explicit  part*/
      }
      DEFINE_SOURCE(O2mass_source, c, t, dS, eqn)
      {
         real source;
         source = C_UDSI(c, t, OM);
         dS[eqn] = 0;
         return source;
      }
DEFINE_SOURCE(total_mass_source, c, t, dS, eqn)
{
real mass_source;
mass_source = C_UDSI(c, t, TM);
dS[eqn] = 0;
return  mass_source;
}
DEFINE_SOURCE(energy_source_fluid, c, t, dS, eqn)
{
real  e_source;
e_source = C_UDSI(c, t, TE);
dS[eqn] = 0;
return e_source;
}
DEFINE_SOURCE(energy_source_solid, c, t, dS, eqn)
{
real e_source;
e_source = 1.0 * C_UDSI(c, t, TE) / (1.0 - e) + 6 * hf * (C_T(c, t) - C_UDSI(c, t, T_S)) / d;
dS[eqn] = 0;
return  e_source;
}
DEFINE_DIFFUSIVITY(m_diff, c, t, i)
{
real diff;
  diff = 6.324e-5 * pow(C_T(c, t), 1.83) / ABS_P(C_P(c, t), op_pres);
return  diff;
}@wuming524

回復此樓

» 猜你喜歡

1樓 2023-11-06 09:33:49

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

QiaoPengyu

新蟲 (初入文壇)

應助: 0 (幼兒園)
金幣: 29
帖子: 10
在線: 2.9小時
蟲號: 28904850
注冊: 2022-03-02
專業(yè): 工程熱物理與能源利用

引用回帖:

2樓: Originally posted by amandahan93 at 2023-11-06 10:11:15
你的udf是用來算什么的呀？

計算氣體吸附的

回復此樓

3樓2023-11-06 10:47:51

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個回答

sallywhe1

金蟲 (小有名氣)

應助: 1 (幼兒園)
金幣: 876.4
散金: 5
帖子: 132
在線: 70.1小時
蟲號: 7342866
注冊: 2017-10-22
專業(yè): 工程熱物理相關交叉領域

你的udf是用來算什么的呀？

贊一下

回復此樓

2樓2023-11-06 10:11:15

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 3 個回答

不應助 確定回帖應助 (注意：應助才可能被獎勵，但不允許灌水，必須填寫15個字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 348求調劑 +4	wyx79618 2026-03-06	8/400	2026-03-10 14:49 by ruiyingmiao
[考研] 320求調劑 +4	魏zy 2026-03-08	4/200	2026-03-09 16:14 by ruiyingmiao
[考研] 0832食品科學與工程293調劑 +3	東東不東 2026-03-07	3/150	2026-03-09 15:58 by ls_13456
[考研] 考研調劑，一志愿山東大學材料與化工，328分，政治51 +5	關你西紅柿929 2026-03-08	6/300	2026-03-09 13:50 by 新篇章DFSS
[考研] 313求調劑 +4	Yyt楊1 2026-03-07	5/250	2026-03-09 13:48 by macy2011
[考研] 337求調劑 +3	睡醒，。 2026-03-09	3/150	2026-03-09 10:02 by 求調劑zz
[考研] 293一志愿華東理工 0817化學工程與技術調劑 +5	fjj0912 2026-03-07	5/250	2026-03-09 09:13 by 30660438
[考博] 26年博士申請 +4	科研狗111 2026-03-07	4/200	2026-03-08 21:56 by 0611517sll
[考研] 334求調劑 +8	Trying] 2026-03-06	8/400	2026-03-08 14:51 by hmn_wj
[考研] 337一志愿華南理工材料求調劑 +4	mysdl 2026-03-07	4/200	2026-03-08 11:00 by 無際的草原
[考研] 一志愿211 化學305分求調劑 +3	0703楊悅305分 2026-03-05	3/150	2026-03-07 21:15 by yinhuanshun
[考研] 一志愿南京大學資源與環(huán)境專業(yè)283分求調劑 +8	追云?? 2026-03-04	8/400	2026-03-07 09:23 by 斬魂滴兔子！
[考研] 085600材料與化工（高分子）290分求調劑 +9	wengyujian 2026-03-04	9/450	2026-03-06 15:59 by 聰明的大松鼠
[考研] 085602高分子方向求調劑 +7	tlgudy 2026-03-04	7/350	2026-03-06 14:17 by guosr9609
[考研] 316求調劑 +3	林小星發(fā)大財 2026-03-05	5/250	2026-03-06 09:05 by Iveryant
[考研] 085602 293分求調劑 +3	SivanNano. 2026-03-05	3/150	2026-03-05 19:41 by caszguilin
[考研] 化學 0703求調劑總分293 一志愿211 +9	土土小蟲 2026-03-03	11/550	2026-03-05 09:26 by vfdbghn
[考研] 316求調劑 +3	林小星發(fā)大財 2026-03-04	3/150	2026-03-05 07:49 by bxbo
[考研] 322分 085600求調劑，有互聯(lián)網+國金及主持省級大創(chuàng)經歷 +6	熊境喆 2026-03-04	6/300	2026-03-04 20:32 by kakakapanpan
[考研] 322,求調劑 +3	菜菜愛玩 2026-03-04	3/150	2026-03-04 12:15 by xiongyaxuan

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時熱門版塊排行榜

QiaoPengyu

[求助] fluent加載udf后初始化閃退

» 猜你喜歡

QiaoPengyu

sallywhe1

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频