A problem in Dirichlet boundary condition for diffusion

[WeiyiVirtonomy]

Please check line 178 in diffusion_dynamics.hpp:
this->all_diffusions_[m]->getInterParticleDiffusionCoff(particle_i, particle_j, e_ij);

particle_j is the index of a contact particle and cannot be used in function getInterParticleDiffusionCoff.

[WeiyiVirtonomy]

@Xiangyu-Hu Prof. Hu, could you please add my account to collaborators? Thanks a lot : )

[Chi Zhang]

Please check line 178 in diffusion_dynamics.hpp: this->all_diffusions_[m]->getInterParticleDiffusionCoff(particle_i, particle_j, e_ij);

particle_j is the index of a contact particle and cannot be used in function getInterParticleDiffusionCoff.

Would do you please specify the branch you are working on?

[WeiyiVirtonomy]

Please check line 178 in diffusion_dynamics.hpp: this->all_diffusions_[m]->getInterParticleDiffusionCoff(particle_i, particle_j, e_ij);
particle_j is the index of a contact particle and cannot be used in function getInterParticleDiffusionCoff.

Would do you please specify the branch you are working on?

That's from the master branch

[ChenxiZhaoTUM]

Hi Weiyi,
I have changed the getInterParticleDiffusionCoff() function of the Dirichlet wall boundary. (see the linked branch)

The diffusion coefficient of the diffusion process interacting with walls is unrelated to the wall boundary, so the old one getInterParticleDiffusionCoff(particle_i, particle_j, e_ij) tends to confuse the SPHinXsys users.
The new one getInterParticleDiffusionCoff(particle_i, particle_i, Vecd::Zero()) is more clear that the diffusion coefficient in this process (diffusion body contacted with Dirichlet boundary) has nothing to do with the contact walls.

For the diffusion part, we just have done the diffusion body interacting with the wall boundaries. Later, we will do the diffusion process referring to the interaction between two diffusion bodies.

[WeiyiVirtonomy]

Hi Weiyi, I have changed the getInterParticleDiffusionCoff() function of the Dirichlet wall boundary. (see the linked branch)

The diffusion coefficient of the diffusion process interacting with walls is unrelated to the wall boundary, so the old one getInterParticleDiffusionCoff(particle_i, particle_j, e_ij) tends to confuse the SPHinXsys users. The new one getInterParticleDiffusionCoff(particle_i, particle_i, Vecd::Zero()) is more clear that the diffusion coefficient in this process (diffusion body contacted with Dirichlet boundary) has nothing to do with the contact walls.

For the diffusion part, we just have done the diffusion body interacting with the wall boundaries. Later, we will do the diffusion process referring to the interaction between two diffusion bodies.

Thanks

[Chi Zhang]

Hi Weiyi, I have changed the getInterParticleDiffusionCoff() function of the Dirichlet wall boundary. (see the linked branch)

The diffusion coefficient of the diffusion process interacting with walls is unrelated to the wall boundary, so the old one getInterParticleDiffusionCoff(particle_i, particle_j, e_ij) tends to confuse the SPHinXsys users. The new one getInterParticleDiffusionCoff(particle_i, particle_i, Vecd::Zero()) is more clear that the diffusion coefficient in this process (diffusion body contacted with Dirichlet boundary) has nothing to do with the contact walls.

For the diffusion part, we just have done the diffusion body interacting with the wall boundaries. Later, we will do the diffusion process referring to the interaction between two diffusion bodies.

@ChenxiZhaoTUM So, now the diffusion process between two bodies can not be handled?