CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

dispersion model with lagragian particle tracking model for incompressible flows

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree28Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   September 10, 2013, 00:02
Default
  #41
Member
 
India
Join Date: Oct 2012
Posts: 84
Rep Power: 5
mayank.dce2k7 is on a distinguished road
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.
mayank.dce2k7 is offline   Reply With Quote

Old   March 28, 2014, 14:20
Default
  #42
Member
 
India
Join Date: Oct 2012
Posts: 84
Rep Power: 5
mayank.dce2k7 is on a distinguished road
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.
mayank.dce2k7 is offline   Reply With Quote

Old   April 2, 2014, 07:20
Default
  #43
New Member
 
Join Date: May 2012
Posts: 3
Rep Power: 6
koulis is on a distinguished road
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
koulis is offline   Reply With Quote

Old   July 16, 2014, 15:16
Arrow How to calculate velocity gradient in particle tracking
  #44
Azy
New Member
 
Azadeh Saeedi
Join Date: Mar 2014
Location: Canada
Posts: 23
Rep Power: 4
Azy is on a distinguished road
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,
Azy is offline   Reply With Quote

Old   November 2, 2014, 05:06
Default
  #45
Senior Member
 
Wen Xu
Join Date: May 2014
Posts: 189
Rep Power: 4
wenxu is on a distinguished road
Quote:
Originally Posted by jiejie View Post
Hi eelcovv

Could I ask how come my kinematicCloudTheta is like the following after I add the voidFraction in the kinematicCloudProperties with the cavityMassPartUncoupledSD_clear case?

Code:
dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    movingWall
    {
        type            calculated;
        value           uniform 0;
    }
    walls
    {
        type            calculated;
        value           uniform 0;
    }
}
Thank you very much.

Jie

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
wenxu is offline   Reply With Quote

Old   November 2, 2014, 08:45
Default
  #46
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,523
Blog Entries: 36
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quote:
Originally Posted by wenxu View Post
"Duplicate entry dual in runtime selection table AveragingMethod etc. many as these lines, but they say the same thing"
Quick answer - Study this thread: Weird segmentation fault!
wyldckat is offline   Reply With Quote

Old   November 3, 2014, 02:37
Default
  #47
Senior Member
 
Wen Xu
Join Date: May 2014
Posts: 189
Rep Power: 4
wenxu is on a distinguished road
Quote:
Originally Posted by wyldckat View Post
Quick answer - Study this thread: Weird segmentation fault!
Dear Bruno,
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 \
PS: i also find a bug(maybe).. the member function in the CloudFunctionObject.H (base class of the CloudFunctionObjects)
Code:
//- Post-move hook
            virtual void postMove
            (
                const typename CloudType::parcelType& p,
                const label cellI,
                const scalar dt,
                const point& position0,
                bool& keepParticle
            );
but in the derived class,say, VoidFraction, the function is:
Code:
//- Post-move hook
            virtual void postMove
            (
                const parcelType& p,
                const label cellI,
                const scalar dt,
                bool& keepParticle
            );
SO when we use the function in the CloudFunctionObject, wen can not use the derived class function with the same name( polymorphism ! )

thank you again! you are really a expert!!!
best regards,
Wen

Last edited by wyldckat; November 17, 2014 at 16:50. Reason: Added [CODE][/CODE]
wenxu is offline   Reply With Quote

Old   November 17, 2014, 17:03
Default
  #48
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,523
Blog Entries: 36
Rep Power: 97
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
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
__________________
___
I'll be at OFW11 in Portugal
wyldckat is offline   Reply With Quote

Old   January 31, 2015, 11:10
Default
  #49
New Member
 
Amjad
Join Date: May 2012
Posts: 21
Rep Power: 5
Amjad Asad is on a distinguished road
Hello Eelc,
could you tell me pls, how u estimated the massflowrate in your case ?
regards,
Amjad
Amjad Asad is offline   Reply With Quote

Old   May 13, 2016, 03:02
Default
  #50
Member
 
Ping Chang
Join Date: Feb 2016
Posts: 45
Rep Power: 2
chpjz0391 is on a distinguished road
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
chpjz0391 is offline   Reply With Quote

Old   Today, 01:13
Default
  #51
Member
 
Ping Chang
Join Date: Feb 2016
Posts: 45
Rep Power: 2
chpjz0391 is on a distinguished road
Quote:
Originally Posted by eelcovv View Post
Hi wellDone,

Thanks for the comments. Yes, probably you need to tweak a bit to get the library compiled, I just compiled it too the latest 2.1.x version, but I don't keep up the versions unless I need to my self. I will have a look to the new wall function you mention.

I just got a request to share the particle statistics cloud function object. I am only attaching this because of the limited upload space (only 100kb). You should add this too the CloudFunctionObjects directory and add the following in the intermediate/parcles/include/cloudfumakeParcelCloudFunctionObjects.H file

#include "ParticleStatistics.H"

and lower

makeCloudFunctionObjectType(ParticleStatistics, CloudType);

Now you can add to the bubbleCloudProperties the following lines.

Code:
cloudFunctions
{
   voidFraction;

   particleStatistics
   {
      resetFields    true;
      prime2Mean     true;
      writeBubblePDF true;
      timeStart      100;
      timeEnd        1e10;
   }
}
The voidFraction is cloudfunction only give you the void fraction of the bubbles, but the particlestatistics also give you the particlevelocity and its rms vales, as will as the PDF of the particle distribution. Very handy, I can assure you (it is the only way to obtain mean profiles of the gas velocity)

Cheers
Eelco
Hey Eelco,

Can I get the particles concentration(kg/m3) by this way?If I can could you tell me how

Ping
chpjz0391 is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Particle tracking results are nonsense Paul FLUENT 14 October 29, 2015 15:55
Water subcooled boiling Attesz CFX 7 January 5, 2013 04:32
Check particle impaction with User Fortran Julian K. CFX 3 January 12, 2012 10:46
massless particle tracking problem Renold FLUENT 0 January 26, 2011 15:23
A model is used to model particle interaction with the geometry dungsidietquy FLUENT 0 November 4, 2009 17:05


All times are GMT -4. The time now is 19:48.