Fluent 2019R2-EDEM 2018 耦合文件基于coupling2.1文件

1、ANSYS 2019R2和EDEM 2018安装在默认位置（C:\Program Files\ANSYS Inc 或 C:\Program Files\DEM Solutions） 2、食用过程中有任何问题请参考流沙大神的耦合文件编译贴或CSDN搜索到的编译指导。

#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个文件信息