#include?“udf.h“
#include?“edem_coupling.h“
#include?“compute_particle_heat_flux.h“

#define?P_PROJECTED_AREA（P）（0.25*M_PI*SQR（P_DIAM（P）））

#if?!RP_HOST
static?void
compute_next_velocity（real?dt?real?a[3]?real?b?real?V0[3]?real?Vn[3]）
{
??real?beta_fact?=?1./（1.?+?b*dt）;
??int?i;

??for?（i=0;?i<3;?++i）
????Vn[i]?=?（V0[i]?+?dt?*a[i]）?*?beta_fact;
}
#endif

void?compute_particle_forces_using_Fluent（Tracked_Particle?*tp?Injection?*I?cxboolean?fluid_source_terms）
{
#if?!RP_HOST
??Particle?*pp;

??loop?（pp?I->p）
????{
??????real?N3V_VEC（acc）;?????????/*?explicit?part?of?particle?acceleration?*/
??????real?N3V_VEC（dvdt）?=?{0.};?/*?velocity?derivative?not?used?*/
??????real?beta;?????????????????/*?drag?factor?*/
??????real?vmf_inv;??????????????/*?virtual?mass?factor?*/
??????real?next_vel[6];??????????/*?first?3?components?contain?new?velocity?last?3?components?carry?old?velocity?*/
??????cxboolean?update_randoms?=?FALSE;?/*?do?not?update?particle?random?numbers?for?Brownian?motion?inside?ParticleAcceleration?*/
??????cxboolean?cphase_interaction?=?TRUE;?/*?fill?relevant?fluid?source?terms?inside?HeatMassUpdate.?*/
??????real?strength?=?P_N（pp）/solver_par.flow_time_step;???????/*?number?of?particles?per?second?typically?1?per?time?step*/
??????particle_state_t?p_old;????/*?needed?to?restore?unintended?changes?in?HeatMassUpdate?*/

??????/*?Fill?tracked?particle?with?information?in?pp?and?interpolate?cphase_state?*/
??????init_tracked_particle（tp?pp?dpm_par.unsteady_tracking?FALSE?FALSE）;
??????P_DT（tp）?=?pp->next_time_step;
??????P_TIME（tp）?+=?P_DT（tp）;
??????
??????if?（DEM_DRAG_FLUENT_P（edem_coupling））
????????{
??????????/*?Now?let‘s?do?a?virtual?step?into?the?particle‘s?direction?for?the?given?flow?time?step?*/
??????????ParticleAcceleration（tp?acc?dvdt?&beta?&vmf_inv?update_randoms）;
??????????compute_next_velocity（P_DT（tp）?acc??beta?P_VEL（tp）?next_vel）;

??????????/*?acc?is?carrying?explicit?part?of?particle?acceleration?
???????????*?let?us?add?implicit?part?and?subtract?gravity?included?in?body?forces
???????????*?to?provide?fluid?related?acceleration?to?EDEM?as?a?force?*/
??????????N3V_VS（acc?-=?P_VEL（pp）?*?beta）;
/*???????????N3V_V?（acc?-=?tp->source.bf_acc）;?*/
??????????N3V_V?（acc?-=?solver_par.G）;

??????????/*?Now?let‘s?compute?the?particle?force?from?the?acceleration?*/
??????????ND_VS（P_DEM_FORCE_X（pp）P_DEM_FORCE_Y（pp）P_DEM_FORCE_Z（pp）=?acc?*?P_MASS（pp））;

??????????/*?Update?velocity?of?tracked?particle?to?be?considered?in?AddSources?*/
??????????N3V_V（P_VEL（tp）?=?next_vel）;
????????}

??????/*?need?to?call?this?always?since?some?DPM?source?terms?are?treated?inside?*/
??????/*?save?state?before?calling?HeatMassUpdate?*/
??????memcpy（（char?*）?&（p_old）?（char?*）?&（tp->state）?sizeof（particle_state_t））;
??????HeatMassUpdate（tp?cphase_interaction）;

??????/*?Consider?Fluent?heat?transfer?here?since?we?have?tp?already?set?*/
??????if?（DEM_HEAT_COUPLED_P（edem_coupling）?&&?DEM_HEAT_TRANSFER_FLUENT_P（ed

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????8540??2016-06-04?19:43??AdaptorInterface.h
?????文件????????6710??2016-06-04?19:43??compute_particle_forces.c
?????文件?????????206??2016-06-04?19:43??compute_particle_forces.h
?????文件????????2891??2016-06-04?19:43??compute_particle_heat_flux.c
?????文件?????????321??2016-06-04?19:43??compute_particle_heat_flux.h
?????目錄???????????0??2019-08-30?15:45??documentation\
?????文件?????1657357??2016-06-04?19:43??documentation\EDEM?Fluent?Coupling?Compilation?Windows.pdf
?????文件?????1077199??2016-06-04?19:43??documentation\EDEM?Fluent?Coupling?Usage.pdf
?????目錄???????????0??2019-08-30?15:45??edem_cfd_adaptor_build\
?????文件?????????515??2016-06-04?19:43??edem_cfd_adaptor_build\README.Linux
?????文件?????????520??2016-06-04?19:43??edem_cfd_adaptor_build\README.Windows
?????文件?????????352??2016-06-04?19:43??edem_cfd_adaptor_build\SConstruct
?????目錄???????????0??2019-08-30?15:45??edem_cfd_adaptor_build\src\
?????文件???????13098??2019-08-01?11:11??edem_cfd_adaptor_build\src\AdaptorInterface.cpp
?????文件????????8794??2016-06-04?19:43??edem_cfd_adaptor_build\src\AdaptorInterface.h
?????文件????????3456??2019-08-01?11:11??edem_cfd_adaptor_build\src\CAdaptorQuaternion.cpp
?????文件????????3008??2016-06-04?19:43??edem_cfd_adaptor_build\src\CAdaptorQuaternion.h
?????文件????????1185??2016-06-04?19:43??edem_cfd_adaptor_build\src\CCellAndThreadData.cpp
?????文件????????2200??2016-06-04?19:43??edem_cfd_adaptor_build\src\CCellAndThreadData.h
?????文件????????5355??2016-06-04?19:43??edem_cfd_adaptor_build\src\CFluentParticleData.cpp
?????文件????????3536??2016-06-04?19:43??edem_cfd_adaptor_build\src\CFluentParticleData.h
?????文件????????6408??2016-06-04?19:43??edem_cfd_adaptor_build\src\CParticleData.cpp
?????文件????????5311??2016-06-04?19:43??edem_cfd_adaptor_build\src\CParticleData.h
?????文件?????????834??2016-06-04?19:43??edem_cfd_adaptor_build\src\SConscript
?????文件????????3414??2016-06-04?19:43??edem_cfd_adaptor_build\src\udfParticle.h
?????目錄???????????0??2019-08-30?15:45??edem_cfd_adaptor_build\win64\
?????文件???????34387??2019-08-01?11:28??edem_cfd_adaptor_build\win64\AdaptorInterface.obj
?????文件???????11298??2019-08-01?11:28??edem_cfd_adaptor_build\win64\CAdaptorQuaternion.obj
?????文件?????????709??2019-08-01?11:28??edem_cfd_adaptor_build\win64\CCellAndThreadData.obj
?????文件??????192158??2019-08-01?11:28??edem_cfd_adaptor_build\win64\CFluentParticleData.obj
?????文件??????163782??2019-08-01?11:28??edem_cfd_adaptor_build\win64\CParticleData.obj
............此處省略73個文件信息