[manu305]

Dear Foamers,


Greetings,


Currently, I am trying to implement thermophoretic force acting on a particle in lagrangian particle tracking of OpenFOAM. Hope someone might have already implemented it since it is an important aspect to consider in particle tracking.



I came across an old implementation in the issue tracking of OpenFOAM foundation. https://bugs.openfoam.org/view.php?id=2549


With little modifications to the code (Thermophoretic_2.zip in above link), the code was succesfully compiled. I thought it would be appropriate to consider thermophoretic force as NonCoupled force instead of Coupled force. So I changed it inside the code. However when I am running the simulation with reactingParcelFoam, I am ending with an error. It would be great if someone can suggest a way out from here.


---------------------------------------------------------------------------

Solving 3-D cloud reactingCloud1
#0 Foam::error:rintStack(Foam::
Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2 ? in "/lib64/libc.so.6"
#3 Foam::interpolationCell<double>::interpolate(Foam: :Vector<double> const&, int, int) const at ??:?
#4 Foam::Thermophoretic<Foam::KinematicCloud<Foam::Cl oud<Foam::ReactingMultiphaseParcel<Foam::ReactingP arcel<Foam::ThermoParcel<Foam::KinematicParcel<Foa m:article> > > > > > >::calcNonCoupled(Foam::ReactingMultiphaseParcel<F oam::ReactingParcel<Foam::ThermoParcel<Foam::Kinem aticParcel<Foam:article> > > > const&, double, double, double, double) const at ??:?
#5 ? at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 ? at ??:?
#9 ? at ??:?
#10 ? at ??:?
#11 ? at ??:?
#12 ? at ??:?
#13 ? at ??:?
#14 __libc_start_main in "/lib64/libc.so.6"
#15 ? at /home/abuild/rpmbuild/BUILD/glibc-2.22/csu/../sysdeps/x86_64/start.S:121
Segmentation fault (core dumped)



------------


I tried to analyse the reason by using the constant thermal condcutivity of fluid as additional input parameter (Thermophoretic.zip in above link), it works fine. So I came to the conclusion that constructing the kappa field and interpolating to particle position on fly is causing the problem. In my case, it is pretty obvious that thermal conductivity of fluid changes in the boundary layer due to large temperature gradients. It can significantly affect the magnitude of thermophoretic force.


Instead of constructing the kappa field on the fly, I thouht of creating it as an additional field (in createFields.H) and read it inside the Thermophoretic.C and interpolate to the posiition of particle. However, I ended up without any clue how to do it ? All the other particle forces, need to read only one field, and in this I need to read two them at a time. It would be great help, if anyone can help me with reading of the new kappa field from Thermophoretic.C ?




Thanks in advance,

Manohar

[jines0708]

Were you able to solve this ?