|
[Sponsors] |
dispersion model with lagragian particle tracking model for incompressible flows |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#41 |
Member
India
Join Date: Oct 2012
Posts: 84
Rep Power: 13 ![]() |
Hi Eelcov,
I am using your solver to simulate particle in a single bifurcating tube with one inlet and 2 outlets of different diameters. I have generated particle position file and corresponding velocity file for particles. How can I specify initial velocity of each particle using your solver? or alternatively how can I make initial particle velocity equal to fluid velocity profile at inlet? Regards, Mayank. |
|
![]() |
![]() |
![]() |
![]() |
#42 |
Member
India
Join Date: Oct 2012
Posts: 84
Rep Power: 13 ![]() |
Hi Eelcov ,
Your transient particle solver is two way coupled. In case I want to make it one way coupled solver, what changes do I need to make in UEqn.H, pEqn and pimpleCoupledKinematicParcelFoam.C? Regards, Mayank. |
|
![]() |
![]() |
![]() |
![]() |
#43 |
New Member
Join Date: May 2012
Posts: 3
Rep Power: 14 ![]() |
Mayank,
Regarding the coupling there is a switch at the very top of the kinematicParticleProperties file within the /constant/ directory set it to false/off... Regarding the velocities of the particles you would either have to change the velocity file, within the lagrangian directory of your first time step directory, manually (line by line, ex. using a macro in VI) . Or you would have to specify one of the existing injectors where you can prescribe a starting velocity...see the lagrangian tutorials for more details on how the injectors work. Hope this helps Koulis |
|
![]() |
![]() |
![]() |
![]() |
#44 |
New Member
Azadeh Saeedi
Join Date: Mar 2014
Location: Canada
Posts: 23
Rep Power: 12 ![]() |
Hello is there any one here could help me a bit about how to do calculate velocity gradient on tracked particles? I want to calculate the shear, so I need it.
Thanks, ![]() |
|
![]() |
![]() |
![]() |
![]() |
#45 | |
Senior Member
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13 ![]() |
Quote:
Dear, jiejie, I encountered the same problem as you i got the outputs stays zero when i use eelcovv's solver. May be the way i compile the new class "parcelStatistics" in a wrong way. I copy the intermediate to my $FOAM_RUN/....../intermediate and the put the parcelStatistics in the CloudFunctions folder and changed the corresponding lines in the makParcelCloudSubmodels.H in the parcels folder. Then i use the DPMFoam as the test case, yes ,the new class is used but firestly, "Duplicate entry dual in runtime selection table AveragingMethod etc. many as these lines, but they say the same thing" errors was happening when the DPMFoam is exected. Secondly, the output was staying zero as i have described above. Please help me, thank you ! regards, wen |
||
![]() |
![]() |
![]() |
![]() |
#46 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#47 | |
Senior Member
Freedom
Join Date: May 2014
Posts: 209
Rep Power: 13 ![]() |
Quote:
Thank you, master!!! I have solved the problem as soon as you give me the link. But that's ok when i run the DPMFoam , however when i changed to the reactingParcelFoam and use the same method things gone wrong have the same point of "duplicate......." Make/options: Code:
EXE_INC = \ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I${LIB_SRC}/meshTools/lnInclude \ -I$(LIB_SRC)/turbulenceModels/compressible/turbulenceModel \ -I$(LIB_SRC)/lagrangian/basic/lnInclude \ -I/root/OpenFOAM/root-2.3.0/src/myIntermediate/lnInclude \ -I$(LIB_SRC)/lagrangian/coalCombustion/lnInclude \ ................... EXE_LIBS = \ -lfiniteVolume \ -lmeshTools \ -lcompressibleTurbulenceModel \ -lcompressibleRASModels \ -lcompressibleLESModels \ -llagrangian \ -L$(FOAM_USER_LIBBIN) -lmyintermediate \ -llagrangianTurbulence \ Code:
//- Post-move hook virtual void postMove ( const typename CloudType::parcelType& p, const label cellI, const scalar dt, const point& position0, bool& keepParticle ); Code:
//- Post-move hook virtual void postMove ( const parcelType& p, const label cellI, const scalar dt, bool& keepParticle ); thank you again! you are really a expert!!! best regards, Wen Last edited by wyldckat; November 17, 2014 at 15:50. Reason: Added [CODE][/CODE] |
||
![]() |
![]() |
![]() |
![]() |
#48 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Hi Wen,
I finally managed to look into this. It seems you've spotted a bug that has already been fixed several months ago (on the 18th of February 2014) in OpenFOAM 2.3.x: https://github.com/OpenFOAM/OpenFOAM...35a5209f11d31a Since you depend on this, perhaps you should either install OpenFOAM 2.3.x or at least get the "intermediate" library source code from it. Best regards, Bruno
__________________
|
|
![]() |
![]() |
![]() |
![]() |
#49 |
New Member
Amjad
Join Date: May 2012
Posts: 21
Rep Power: 14 ![]() |
Hello Eelc,
could you tell me pls, how u estimated the massflowrate in your case ? regards, Amjad |
|
![]() |
![]() |
![]() |
![]() |
#50 |
Member
Ping Chang
Join Date: Feb 2016
Location: Perth
Posts: 93
Rep Power: 10 ![]() |
Hey Mayank,
could you post ur solver for incompressible flow particles tracking solver for steady state. Actually I am looking for such a solver, however, I am a newbie of OF, I am not able to compile the solver by myself. Thank you |
|
![]() |
![]() |
![]() |
![]() |
#51 | |
Member
Ping Chang
Join Date: Feb 2016
Location: Perth
Posts: 93
Rep Power: 10 ![]() |
Quote:
Can I get the particles concentration(kg/m3) by this way?If I can could you tell me how Ping |
||
![]() |
![]() |
![]() |
![]() |
#52 | |
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 11 ![]() |
Quote:
I am quite interested in Eelcos work and everybodys enlightment on the topic of this thread. I would very much like to dig into the solver and the elegant code. But before that I think I need some help about the physical description about the attached test case and also the theoretical background. So I am going to restart this thread. Could anybody give me some guidance on the physical and mathematical interpretation of the pimpleLPTBubbleFoam or the bubcol3d case? I work with OF-2.3.0 and know the pimple solver, but I am quite new to the lagrangian world of OF. Thank you in advance! Best, Karekle |
||
![]() |
![]() |
![]() |
![]() |
#53 |
Member
Yong Wang
Join Date: Apr 2009
Posts: 34
Rep Power: 17 ![]() |
Hi Eelco,
Your solver looks great. Before trying, could you please let me know whether it can consider heat transfer or not? I want to use it for particle tracking in a room where the turbulent flow is driven by temperature difference. Thanks. Best, Yong |
|
![]() |
![]() |
![]() |
![]() |
#54 |
Senior Member
Eelco van Vliet
Join Date: Mar 2009
Location: The Netherlands
Posts: 124
Rep Power: 19 ![]() |
Dear All
Thanks for the likes for this solver I posted a couple of years ago. Could not answer all the question as I am not online on the forum that much anymore Anyhow, I needed to use it myself lately and found that the OpenFOAM 5.x branch has an improved particle tracking algorithm, see here I have updated the particle tracking solver for simulating bubbles in an incompressible flow for the OpenFOAM 5.x version (from openfoam.org) The good news is that you do not have to compile any libraries anymore. Just the solver is enough. I have included a small example of usage as well. Good luck with it! Regards Eelco |
|
![]() |
![]() |
![]() |
![]() |
#55 |
Member
vishal
Join Date: Mar 2013
Posts: 73
Rep Power: 13 ![]() |
Hello,
Did anyone solved the void fraction issue in CollidingParcel case ? can anyone please share their file or suggest me where and how to incorporate the line " td.cloud().functions().postMove(p, cellI, dt);" in CollidingParcel.C. Moreover, dt is not mention in the CollidingParcel.H file. so it will give an error while compiling I am using OF-2.2.2. |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Particle tracking results are nonsense | Paul | FLUENT | 16 | October 17, 2019 14:12 |
Water subcooled boiling | Attesz | CFX | 7 | January 5, 2013 03:32 |
Check particle impaction with User Fortran | Julian K. | CFX | 3 | January 12, 2012 09:46 |
massless particle tracking problem | Renold | FLUENT | 0 | January 26, 2011 14:23 |
A model is used to model particle interaction with the geometry | dungsidietquy | FLUENT | 0 | November 4, 2009 16:05 |