Hello all,
attached to this post there are scripts I made internal heat transfer in multispherical particles the changes involve:
a) thermal capacity and thermal conductivity are now property/atom
b) new property/atom thermal conductivity tangential added for anisotropic exchange and property/global contactsurface added
c) heat transfer in anisotropic solid, you can manage this with particle types: if particle A is of type 1 and particle 2 is of type 2 conductivity tangential will be used, if instead particle A is of the same type as particle B normal conductivity will be used
d) every particle A in the multisphere assembly is counted as a discretized piece of the multisphere particle, who exchanges heat with the others, if B is a particle in contact: (mean conductivity A, B)*sqrt(contactSurfaceA*contactSurfaceB)
to install these new scripts all you have to do overwrite the existing in the src directory.
attached there is also an exemple, if you need any information or I did a mistake I will gladly answer in private message, however all this need to be tested
sorry if some comments are in italian :D
hope this could be usefull,
| Attachment | Size |
|---|---|
 fix_heat_gran.cpp_.txt | 6.77 KB |
| fix_heat_gran.h.txt | 2.02 KB |
| fix_heat_gran_conduction.cpp_.txt | 14.82 KB |
| fix_heat_gran_conduction.h.txt | 2 KB |
| fix_scalar_transport_equation.cpp_.txt | 10.15 KB |
| fix_scalar_transport_equation.h.txt | 1.96 KB |
| in.piroplacn.txt | 3.65 KB |
| piro.data_.txt | 74.28 KB |
ckloss | Thu, 09/25/2014 - 08:40
Hi Riccardo,
Hi Riccardo,
thanks! If that's ok for you, I'll have a look at it before we have the next release. We might include some of your code in the PUBLIC version
Cheers
Christoph
Riccardo Maione | Mon, 09/29/2014 - 13:39
Ok, let me know if something
Ok, let me know if something isn't right I will gladly help
ckloss | Thu, 10/30/2014 - 14:34
Hi Riccardo,
Hi Riccardo,
thanks again for sharing your model! But actually I think the current way how things are implemented in LIGGGHTS with respect to conduction should still be the default behavior, and your model the way it is implemented would change this default. So I'd implement it as a derived class FixHeatGranConductionNonUniform, which derives from FixHeatGranConduction.
I will write this on my list for the next release - this release I had just limited time to implement new stuff.
Thanks and best wishes
Christoph
Riccardo Maione | Tue, 11/04/2014 - 09:00
I know, this was the way I
I know, this was the way I wanted it to work for my problem so I didn't care too much, hope it is still interesting however,
Best regards,
Riccardo
EmersonSouza | Wed, 12/21/2022 - 09:03
Files are not available
Hello Riccardo,
Could you please upload the files that you about multisphere heat transfer again?
I new in liggghts and I would like to start a simulation using this type of approach, and as the current version is restricted for multisphere heat transfer, your code would be very helpful.
Thank you for your attention.