UDF Energy Source for DPM

ArikGuyz · May 28, 2020, 07:18

Hello Everyone,

I am trying to add an energy source to a moving particle using DEFINE_DPM_SOURCE. I did go over the UDF manual but could not find a suitable example for my case. Depending on the particle position amount of heat source added to that particle will vary. Can anyone help me out?

Also I don't fully understand how DEFINE_DPM_SOURCE works, when I tried to add constant heat source to the particle, entire domain gets heat up instead of particle.

Thanks in advance!

vinerm · May 28, 2020, 07:26

The name DPM_SOURCE is actually a misnomer; there are no source terms in Lagrangian equation. Source terms are meant for Eulerian form. If you want to add energy to particles, then you have to calculate the expected temperature change for a particle and update the particle temperature. There is no macro or function to add energy to a particle; it's not required because particle has no temperature variation.

ArikGuyz · May 28, 2020, 07:47

Thank you for the fast response. Well, I know how much a particle must generate heat. So there is no way for me to heat up the particles with a known heat generation rate?

vinerm · May 28, 2020, 07:56

There is no direct option to setup heat generation in particle. So, you can either use a particle scalar, in case the formulation for the determination of the heat generated is rather involved, or any DEFINE_DPM_ function that does not otherwise affect the particles. Within this function, you can calculate what should be the particle temperature due to this heat generation and its current temperature and then set that temperature as the particle temperature using TP_T (or using P_T, if you are using older versions).

ArikGuyz · September 18, 2020, 04:32

Hello Again,

First thank you for the advice I did manage to change particle temperature using UDF. However I have a question about it. I am using DPM_LAW and LAW1 custom inert heating law to change the temperature of the particle. I am applying additional heat source as a dT to particles current temperature. And I am modelling DPM with two way coupling. However, when I add additional temperature for ex. 200K, I don't see any change in flow temperature or velocity or even in particle velocity. The only difference is particle temperature. With two way coupling enabled shouldn't particle properties influence the continuous phase?

May 28, 2020, 07:18	UDF Energy Source for DPM	#1
ArikGuyz New Member Join Date: Oct 2019 Posts: 12 Rep Power: 7	Hello Everyone, I am trying to add an energy source to a moving particle using DEFINE_DPM_SOURCE. I did go over the UDF manual but could not find a suitable example for my case. Depending on the particle position amount of heat source added to that particle will vary. Can anyone help me out? Also I don't fully understand how DEFINE_DPM_SOURCE works, when I tried to add constant heat source to the particle, entire domain gets heat up instead of particle. Thanks in advance! qinchangliang likes this.

May 28, 2020, 07:26	Particle Source	#2
vinerm Senior Member Vinerm Join Date: Jun 2009 Location: Nederland Posts: 2,946 Blog Entries: 1 Rep Power: 36	The name DPM_SOURCE is actually a misnomer; there are no source terms in Lagrangian equation. Source terms are meant for Eulerian form. If you want to add energy to particles, then you have to calculate the expected temperature change for a particle and update the particle temperature. There is no macro or function to add energy to a particle; it's not required because particle has no temperature variation. ArikGuyz and qinchangliang like this. __________________ Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.

May 28, 2020, 07:56	Heat Generation	#4
vinerm Senior Member Vinerm Join Date: Jun 2009 Location: Nederland Posts: 2,946 Blog Entries: 1 Rep Power: 36	There is no direct option to setup heat generation in particle. So, you can either use a particle scalar, in case the formulation for the determination of the heat generated is rather involved, or any DEFINE_DPM_ function that does not otherwise affect the particles. Within this function, you can calculate what should be the particle temperature due to this heat generation and its current temperature and then set that temperature as the particle temperature using TP_T (or using P_T, if you are using older versions). ArikGuyz likes this. __________________ Regards, Vinerm PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[OpenFOAM.com] swak4foam compiling issues on a cluster	saj216	OpenFOAM Installation	5	January 17, 2023 17:05
Problem compiling a custom Lagrangian library	brbbhatti	OpenFOAM Programming & Development	2	July 7, 2014 12:32
[swak4Foam] Swak4FOAM 0.2.3 / OF2.2.x installation error	FerdiFuchs	OpenFOAM Community Contributions	27	April 16, 2014 16:14
[swak4Foam] build problem swak4Foam OF 2.2.0	mcathela	OpenFOAM Community Contributions	14	April 23, 2013 14:59
OpenFOAM on MinGW crosscompiler hosted on Linux	allenzhao	OpenFOAM Installation	127	January 30, 2009 20:08

May 28, 2020, 07:47		#3
ArikGuyz New Member Join Date: Oct 2019 Posts: 12 Rep Power: 7	Thank you for the fast response. Well, I know how much a particle must generate heat. So there is no way for me to heat up the particles with a known heat generation rate?

September 18, 2020, 04:32		#5
ArikGuyz New Member Join Date: Oct 2019 Posts: 12 Rep Power: 7	Hello Again, First thank you for the advice I did manage to change particle temperature using UDF. However I have a question about it. I am using DPM_LAW and LAW1 custom inert heating law to change the temperature of the particle. I am applying additional heat source as a dT to particles current temperature. And I am modelling DPM with two way coupling. However, when I add additional temperature for ex. 200K, I don't see any change in flow temperature or velocity or even in particle velocity. The only difference is particle temperature. With two way coupling enabled shouldn't particle properties influence the continuous phase?