24СʱÈÈÃÅ°æ¿éÅÅÐаñ    

²é¿´: 1403  |  »Ø¸´: 5
¡¾ÐüÉͽð±Ò¡¿»Ø´ð±¾ÌûÎÊÌ⣬×÷Õßkeen_xq½«ÔùËÍÄú 35 ¸ö½ð±Ò

keen_xq

Ìú³æ (СÓÐÃûÆø)

[ÇóÖú] DPMºêÃüÁîBegin_particle_cell_loop(p,c,t)µÄʹÓÃÎÊÌâÒÑÓÐ1È˲ÎÓë

¸÷λͬѧ ÓÐûÓÐʹÓÃFluentÖÐ Õâ¸ö begin_particle_cell_loop(p,c,t)ºêÄØ£¬¸Ãºê³öÏÖÔÚDEFINE_DPM_SPRAY_COLLIDEÖУ¬¿´°ïÖúÉϵĽâÊÍ£¬¸ÃºêÓ¦¸ÃÊDZéÀúµ¥ÔªÖÐËùÓеĿÅÁ£¡£ÎÒÏëÓøúêÇó³öµ¥ÔªÖÐ ¿ÅÁ£µÄÌå»ý·ÖÊý£¬ÃüÁîÈçÏ¡£ ¸ÃÃüÁîÔÚ±àÒëÖÐûÎÊÌ⣬¼ÆËãʱ±ÀÀ£¡£ÎÒµ÷ÊÔ¹ý¶à´Î£¬×îºó·¢ÏÖ begin_particle_cell_loop(pi,c,t)²»ÄÜÔÙDEFINE_DPM_DRAGÖÐʹÓã¬Ò»µ©Ê¹Óã¬Çó½â¾Íû°ì·¨½øÐУ¬ÎÒÒ²³¢ÊÔÔÚÆäËûºê ÈçDEFINE_DPM_SCALE_UPDATA ºÍDEFINE_DPM_BODY_FORCEÖе÷Óþùʧ°Ü¡£ÓÐûÓдóţʹÓùý¸Ãºê£¬ÇóÖ¸µãÃÔ½ò£¬±»ÕÛÄ¥Á˺þÃÁË¡£


DEFINE_DPM_DRAG(particle_drag_force,Re,p)
{
        real dd,Vsum=0.0;
        cell_t c=P_CELL(p);
        Thread *t =P_CELL_THREAD(p);
        Particle *pi;
        begin_particle_cell_loop(pi,c,t)
        {
                Vsum+=M_PI*pow(P_DIAM(pi),3.0)/6.0;
        }end_particle_cell_loop(pi,c,t)

        C_UDMI(c,t,0)=Vsum/C_VOLUME(c,t);  //ÇóÌåÏà»ý·ÖÊý£»//
        if(C_UDMI(c,t,0)>1)
                C_UDMI(c,t,0)=1.0;

        if(Dragsign==1)
                dd=Cdrag_force1(p,Re);
        else if(Dragsign==2)
                dd=Cdrag_force2(p,Re);
        else if(Dragsign==3)
                dd=Cdrag_force3(p,Re);
        else if(Dragsign==4)
                dd=Cdrag_force4(p,Re);
        else
                dd=Cdrag_force5(p,Re);
        return dd*pow(1-C_UDMI(c,t,0),5.0/3.0)/pow(1-C_UDMI(c,t,0),2.0);
}
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÎðÍü³õÐÄ£¬·½µÃʼÖÕ¡£
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

liukan12

Òø³æ (СÓÐÃûÆø)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ÄãµÄÀí½âÓÐÎóÇø ʵ¼ÊÉÏÕâ¸öºêÐèÒªÏÈÒª°ÑparticleÔÚcellÖм佨Á¢Ë÷ÒýÁбí ËùÒÔÐèÒªÏÈ°Ñparticle·âװһϠ¾ÍÊÇ˵ Ê×ÏÈÄãÐèÒªÒ»¸öparticle bin
·ñÔòÄãÊǻᱨ´íµÄ
Æäʵһ¸öµ¥ÔªÄڵĿÅÁ£×ÜÊýÊÇÓмòµ¥µÄÇ󷨵Ä

» ±¾ÌûÒÑ»ñµÃµÄºì»¨£¨×îÐÂ10¶ä£©

ANSYSFluentÁ÷Ìå·ÂÕæÑо¿ÈËÔ±¡£
2Â¥2019-06-15 14:43:25
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

keen_xq

Ìú³æ (СÓÐÃûÆø)

Ĥ°Ý´óÉñ¡£ÇëÎÊ´óÉñ£¬ÄãËù˵µÄparticle binÈçºÎ´´½¨£¬»òÕßÈçºÎµÃµ½£¿ÔÚfluent UDF ÊÖ²áÖйØÓÚÕâ¸öºêÎÒÖ»·¢ÏÖÔÚDEFINE_DPM_SPRAY_COLLIDEÖдæÔÚ£¬¹ØÓÚÕâ¸öºêµÄ½éÉÜÒ²ÊÇÒ»±Ê´ø¹ý¡£´óÉñ£¬ÄãÄÜÏêϸ˵Ï£¬»òÔòÍƼö¸ö×ÊÁÏÂð£¬udfÊÖ²áÖеĽéÉܺܲ»Ïêϸ£¬Óкܶ඼ÐèÒª×Ô¼º½áºÏ³ÌÐòÍƳöÊÇʲôÒâ˼¡£¡£´óÉñ£¬ÎÒudfµÄÄ¿µÄÊÇÏë±éÀúÿ¸öÁ÷Ìåµ¥ÔªÖÐËùÓеÄÆøÅÝ£¬Í¨¹ýÀÛ¼ÓÆäÌå»ýÈ»ºóÇóµÃÆøÏàµÄÌå»ý·ÖÊý¡£´óÉñ£¬ÄãÄÜÍƼöϼòµ¥µÄ·½·¨Â𣿲»Ê¤¸Ð¼¤¡£¡£»¹Óнð±Ò·îÉÏ

·¢×ÔСľ³æAndroid¿Í»§¶Ë
ÎðÍü³õÐÄ£¬·½µÃʼÖÕ¡£
3Â¥2019-06-16 00:11:52
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

keen_xq

Ìú³æ (СÓÐÃûÆø)

Ëͺ컨һ¶ä
ÒýÓûØÌû:
2Â¥: Originally posted by liukan12 at 2019-06-15 14:43:25
ÄãµÄÀí½âÓÐÎóÇø ʵ¼ÊÉÏÕâ¸öºêÐèÒªÏÈÒª°ÑparticleÔÚcellÖм佨Á¢Ë÷ÒýÁбí ËùÒÔÐèÒªÏÈ°Ñparticle·âװһϠ¾ÍÊÇ˵ Ê×ÏÈÄãÐèÒªÒ»¸öparticle bin
·ñÔòÄãÊǻᱨ´íµÄ
Æäʵһ¸öµ¥ÔªÄڵĿÅÁ£×ÜÊýÊÇÓмòµ¥µÄÇ󷨵Ä

Ĥ°Ý´óÉñ¡£ÇëÎÊ´óÉñ£¬ÄãËù˵µÄparticle binÈçºÎ´´½¨£¬»òÕßÈçºÎµÃµ½£¿ÔÚfluent UDF ÊÖ²áÖйØÓÚÕâ¸öºêÎÒÖ»·¢ÏÖÔÚDEFINE_DPM_SPRAY_COLLIDEÖдæÔÚ£¬¹ØÓÚÕâ¸öºêµÄ½éÉÜÒ²ÊÇÒ»±Ê´ø¹ý¡£´óÉñ£¬ÄãÄÜÏêϸ˵Ï£¬»òÔòÍƼö¸ö×ÊÁÏÂð£¬udfÊÖ²áÖеĽéÉܺܲ»Ïêϸ£¬Óкܶ඼ÐèÒª×Ô¼º½áºÏ³ÌÐòÍƳöÊÇʲôÒâ˼¡£¡£ÎÒudfµÄÄ¿µÄÊÇÏë±éÀúÿ¸öÁ÷Ìåµ¥ÔªÖÐËùÓеÄÆøÅÝ£¬Í¨¹ýÀÛ¼ÓÆäÌå»ýÈ»ºóÇóµÃÆøÏàµÄÌå»ý·ÖÊý¡£´óÉñ£¬ÄãÄÜÍƼöϼòµ¥µÄ·½·¨Â𣿲»Ê¤¸Ð¼¤¡£¡£»¹Óнð±Ò·îÉÏ

·¢×ÔСľ³æAndroid¿Í»§¶Ë
ÎðÍü³õÐÄ£¬·½µÃʼÖÕ¡£
4Â¥2019-06-16 00:15:38
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Ç®äì

гæ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ÒýÓûØÌû:
4Â¥: Originally posted by keen_xq at 2019-06-16 00:15:38
Ĥ°Ý´óÉñ¡£ÇëÎÊ´óÉñ£¬ÄãËù˵µÄparticle binÈçºÎ´´½¨£¬»òÕßÈçºÎµÃµ½£¿ÔÚfluent UDF ÊÖ²áÖйØÓÚÕâ¸öºêÎÒÖ»·¢ÏÖÔÚDEFINE_DPM_SPRAY_COLLIDEÖдæÔÚ£¬¹ØÓÚÕâ¸öºêµÄ½éÉÜÒ²ÊÇÒ»±Ê´ø¹ý¡£´óÉñ£¬ÄãÄÜÏêϸ˵Ï£¬»òÔòÍƼö¸ö×ÊÁÏ ...

ÇëÎÊÂ¥Ö÷½â¾öÕâ¸öÎÊÌâÁËÂð£¬ÎÒÒ²ÊÇÏ£Íûͨ¹ýUDF»ñµÃÁ÷ÌåÓòµÄÆøÒº±È£¬½á¹û£¬³ö¸÷ÖÖÎÊÌâ
5Â¥2020-05-19 16:44:38
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

Ç®äì

гæ (³õÈëÎÄ̳)

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

ÒýÓûØÌû:
5Â¥: Originally posted by Ç®äì at 2020-05-19 16:44:38
ÇëÎÊÂ¥Ö÷½â¾öÕâ¸öÎÊÌâÁËÂð£¬ÎÒÒ²ÊÇÏ£Íûͨ¹ýUDF»ñµÃÁ÷ÌåÓòµÄÆøÒº±È£¬½á¹û£¬³ö¸÷ÖÖÎÊÌâ...

There is a loop available for looping over all the particles in a cell. But it does not work if used directly. Actually the array to bin the particles in the cell is not allocated by default. It has to be done before using this loop and then it will run ok.

So, at the top of the routine, please put,

Alloc_Storage_Vars(domain, SV_DPM_PARTICLE_BIN, SV_NULL);
bin_particles_in_cells(domain);

and at the botton, free up the memory with,

Free_Storage_Vars(domain, SV_DPM_PARTICLE_BIN, SV_NULL);
I am attaching my test UDF. It worked nicely for a simple 2d case in Linux.

#include "udf.h"
#include "surf.h"

DEFINE_ON_DEMAND(cell)
{

Domain *d=Get_Domain(1);
cell_t c;
int id;
Particle *p;

Thread *c_thread=Lookup_Thread(d,2); /*In my case 2 was the id of the fluid thread. Please change this number accordingly */

Alloc_Storage_Vars(d, SV_DPM_PARTICLE_BIN, SV_NULL);
bin_particles_in_cells(d);

begin_c_loop(c, c_thread)

{
/* Message("Cell_index=%dn",c);*/
begin_particle_cell_loop(p,c,c_thread)
{
id=p->part_id;
Message(" cell index=%d part_id=%d ",c,id);
}
end_particle_cell_loop(p,c,c_thread)
}

end_c_loop(c, c_thread)

Free_Storage_Vars(d, SV_DPM_PARTICLE_BIN, SV_NULL);

}
6Â¥2020-05-19 18:00:43
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌøת ÎÒÒª¶©ÔÄÂ¥Ö÷ keen_xq µÄÖ÷Ìâ¸üÐÂ
²»Ó¦Öú È·¶¨»ØÌûÓ¦Öú (×¢Ò⣺ӦÖú²Å¿ÉÄܱ»½±Àø£¬µ«²»ÔÊÐí¹àË®£¬±ØÐëÌîд15¸ö×Ö·ûÒÔÉÏ)
×î¾ßÈËÆøÈÈÌûÍƼö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 311Çóµ÷¼Á +6 ¸Â¶àÀÁÑòÑò 2024-03-28 8/400 2024-03-28 21:38 by xiejq2010
[½Ìʦ֮¼Ò] ±ÏÒµ6Ä꣬¸Ð¾õ¿ÆÑеÄ·Խ×ßÔ½Õ­ÁË£¬²»ÖªµÀÔõô°ì +13 ÓàÒ¦ºÍг 2024-03-27 13/650 2024-03-28 21:02 by ÓãÏèdzµ×1
[¿¼ÑÐ] Çóµ÷¼Á +3 085600h 2024-03-28 3/150 2024-03-28 20:56 by Kilig0317
[²ÄÁϹ¤³Ì] Äϲý´óѧ²ÄÁϹ¤³Ì297Çóµ÷¼Á +9 ºÎij~~~ 2024-03-25 10/500 2024-03-28 18:10 by »ú³µÌÀÄ·É­
[¿¼ÑÐ] ±¾¿ÆË«·ÇÒ»±¾£¬Ò»Ö¾Ô¸ÖÐÄÏ´óѧ085600²ÄÁÏÓ뻯¹¤309·Ö£¬ÇóÄÏ·½Ñ§Ð£µ÷¼Á +10 À´Çï_ 2024-03-23 10/500 2024-03-28 17:47 by LLLOOOLLLPPP
[¿¼ÑÐ] ±¾¿ÆÌ«Ô­Àí¹¤´óѧ£¬300·ÖÇóµ÷¼Á +5 ²Ø´ºÉ½ 2024-03-27 5/250 2024-03-28 16:12 by shiyihelulu
[¿¼ÑÐ] 299Çóµ÷¼Á +7 Ò»¿ÚÄϹϱýr 2024-03-23 7/350 2024-03-28 15:13 by qlm5820
[¿¼ÑÐ] 302Çóµ÷¼Á +6 ·¶ÖÇÇǦŠ2024-03-26 6/300 2024-03-28 14:55 by qlm5820
[¿¼ÑÐ] 301Çóµ÷¼Á +6 СӨ@ 2024-03-28 6/300 2024-03-28 11:21 by Áúɽ110
[¿¼²©] ÕÐÊÕÒ»Ãû·ÖÎö»¯Ñ§»òÓлú»¯Ñ§·½ÏòµÄ²©Ê¿£¨»ªÄÏÀí¹¤´óѧ£© +3 contact_angle 2024-03-22 5/250 2024-03-28 11:02 by TME¡¢
[Óлú½»Á÷] ºÏ³ÉʵÑé½ô¼±ÇóÖú£¬¸÷λ´óÉñ°ï棡 50+4 1595515531 2024-03-26 5/250 2024-03-28 09:33 by ¸ßÊ÷×®
[¿¼ÑÐ] 323²ÄÁÏѧ˶Çóµ÷¼Á£¬Ò»Ö¾Ô¸985 +8 hzxhan 2024-03-27 10/500 2024-03-28 07:12 by EBSD
[¿¼ÑÐ] »·¾³Ñ§Ë¶334Çóµ÷¼Á +6 Ì«Ñô553 2024-03-23 8/400 2024-03-27 20:50 by AÎÒÒªÉÏѧ
[¿¼ÑÐ] 309Çóµ÷¼Á +8 С±¿ÐÜ12138 2024-03-23 8/400 2024-03-27 16:26 by tkdmxl
[¿¼ÑÐ] 321·Ö²ÄÁÏÀà ±¾¿Æ¸ß·Ö×ÓÇóµ÷¼Á +16 Tumi1218 2024-03-22 19/950 2024-03-27 16:18 by tkdmxl
[¿¼ÑÐ] Ò»Ö¾Ô¸Î人Àí¹¤²ÄÁÏѧ ѧ˶ 300·ÖÇóµ÷¼Á +4 ÏëÏòÉÏpao 2024-03-25 8/400 2024-03-27 08:49 by δÀ´¿ÉÆÚ¼òÖ±²»¶
[¿¼ÑÐ] 09ʳ¼Ó 285AÇøÇóµ÷¼Á 211±¾¿Æ +3 ÐÄÖ®ÑúÑúѽ 2024-03-26 4/200 2024-03-27 01:47 by ÐÄÖ®ÑúÑúѽ
[¿¼ÑÐ] 335Çóµ÷¼Á +6 ÔÆÔõôÍìÁô 2024-03-26 7/350 2024-03-26 18:29 by Whisper890
[¿¼ÑÐ] 284Çóµ÷¼Á +7 TwistedC 2024-03-22 7/350 2024-03-26 13:00 by helloѽ6
[¿¼ÑÐ] 327Çóµ÷¼Á +4 ÍüÈ´ÆßÃëÇé 2024-03-22 4/200 2024-03-24 14:01 by СÓê¸ç
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û