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

IcoLagrangianFoam and particle forces

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 10, 2007, 10:22
Default Hi what i'm trying to do is
  #21
New Member
 
alan fergus
Join Date: Mar 2009
Location: galway, ireland
Posts: 23
Rep Power: 8
alan_fergus is on a distinguished road
Hi

what i'm trying to do is include resuspension into the model. to this end when a particle is deposited i can take note of it diameter, its mass , position it was deposited at , there i can work out what energy is needed in the fluid to resuspened (re entrain) it into the fluid. at the moment i have this in a struct (see below)

struct pt{
int dep_num;
int dep_numres;
scalar tsdep;
scalar Ead;
vector posr;
scalar dp;
scalar mass;
scalar Eres;
scalar dep_res; //0-empty 1-deposit 2-resuspended
} pdr[SIZE];

what i would like to be able to do is access this information in IncompressibleCloud.C. so i can just use the inject to place back particle into cloud.

any advice would be most welcome, would like to start my test cases soon so i can be finished by deal line and this problem is something that i can't find away around.

Thanks in Advance
Alan Fergus
alan_fergus is offline   Reply With Quote

Old   May 10, 2007, 14:58
Default Hi Alan! Do I understand yo
  #22
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Alan!

Do I understand you correctly? Theses properties are to be stored per particle. Then they should not go into the Cloud-Subclass (IncompressibleCloud in IcoLagrangianFoam) but into the Particle-subclass (HardBallParticle in IcoLagrangianFoam). Look for the handling of mass_ and d_ in HardBallParticle as an example. That way it is easier to add or remover particles. Storing the particle properties in an array with a fixed size (that is how I interpret your pdr-array) is so FORTRAN. Don't get me wrong: some of the cleverest people I know program FORTRAN, it's just not the way things are done in OpenFOAM

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   July 12, 2007, 10:07
Default Hello all, I am trying to c
  #23
mayank
Guest
 
Posts: n/a
Hello all,

I am trying to compile icoLagrangianFoam on Openfoam 1.4 and get the following error:

/home/mgo/OpenFOAM/OpenFOAM-1.4/src/finiteVolume/lnInclude/volPointInterpolation .H: In constructor 'Foam::IncompressibleCloud::IncompressibleCloud(co nst Foam::volPointInterpolation&, const Foam::volVectorField&)':
/home/mgo/OpenFOAM/OpenFOAM-1.4/src/finiteVolume/lnInclude/volPointInterpolation .H:70: error: 'Foam::volPointInterpolation::volPointInterpolatio n(const Foam::volPointInterpolation&)' is private
IncompressibleCloud.C:68: error: within this context
IncompressibleCloud.C: In member function 'void Foam::IncompressibleCloud::track()':
IncompressibleCloud.C:154: error: 'track' is not a member of 'Foam::Cloud<foam::hardballparticle>'
make: *** [Make/linux64Gcc4DPOpt/IncompressibleCloud.o] Error 1

Need help to fix it.
Thanks.
mayank
  Reply With Quote

Old   July 12, 2007, 11:05
Default That is because the sources on
  #24
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
That is because the sources on the Wiki are for version 1.3 and there have been changes in the particle-codes for version 1.4. I have an adapted version for 1.4 on my hard-disk that I didn't publish yet because in my opinion it still had a problem with the second test case (particles didn't stop bouncing). I have uploaded it to a (not yet fully operational) Subversion repository. You can get it with the command (if you have Subversion installed on your machine)

svn checkout https://openfoam-extend.svn.sourceforge.net/svnroot/openfoam-extend/trunk/Breede r/solvers/other/IcoLagrangianFoam

Be aware, that this is a very untested version (as opposed to the 1.3-version which was little-tested)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   January 28, 2008, 10:19
Default Hi Bernhard, I've got icoL
  #25
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 8
francesco_b is on a distinguished road
Hi Bernhard,

I've got icoLagrangianFoam from the repository, when I compiled it I got some errors, the same as those Aurelia obtained in the Thread "Changes in OF 1.4 and compilation issues". I've fixed the error due to track, but not the other, the one at line 68.

Can you please tell me how can I solve it?

Thanks in advance

Francesco
francesco_b is offline   Reply With Quote

Old   January 28, 2008, 10:39
Default I've succeed in fixing the pro
  #26
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 8
francesco_b is on a distinguished road
I've succeed in fixing the problem, sorry for my stupid question
francesco_b is offline   Reply With Quote

Old   September 14, 2008, 12:02
Default Hi all, I am new to OpenFOA
  #27
New Member
 
Max Eggersdorfer
Join Date: Mar 2009
Posts: 1
Rep Power: 0
maxe is on a distinguished road
Hi all,

I am new to OpenFOAM and want to simulate particles in a fluid. I use OF 1.5 and have difficulties compiling icoLagrangianFoam:

Making dependency list for source file icoLagrangianFoam.C
Making dependency list for source file HardBallParticle.C
Making dependency list for source file IncompressibleCloud.C
SOURCE=icoLagrangianFoam.C ; g++ -m64 -Dlinux64 -DDP -Wall -Wno-strict-aliasing -Wextra -Wno-unused-parameter -Wold-style-cast -march=opteron -O3 -DNoRepository -ftemplate-depth-40 -I/scratch/OpenFOAM/OpenFOAM-1.5/src/cfdTools/incompressible -I/scratch/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/lnInclude -I/scratch/OpenFOAM/OpenFOAM-1.5/src/lagrangian/basic/lnInclude -I/scratch/OpenFOAM/OpenFOAM-1.5/src/cfdTools/general/lnInclude -I/scratch/OpenFOAM/OpenFOAM-1.5/src/finiteVolume/lnInclude -IlnInclude -I. -I/scratch/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/lnInclude -I/scratch/OpenFOAM/OpenFOAM-1.5/src/OSspecific/Unix/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/icoLagrangianFoam.o
In file included from icoLagrangianFoam.C:35:
HardBallParticle.H:56: error: expected template-name before '<' token
HardBallParticle.H:56: error: expected `{' before '<' token
HardBallParticle.H:56: error: expected unqualified-id before '<' token
icoLagrangianFoam.C:118: error: expected `}' at end of input
make: *** [Make/linux64GccDPOpt/icoLagrangianFoam.o] Error 1

Has anyone an updated version for 1.5 or may anyone help me to fix the problems.

And I have a more general question:
Are particles in this solver treated like pointwise particles? Is there a possibility with OpenFOAM to simulate particles with a finite Volume?

Thanks in advance
max
maxe is offline   Reply With Quote

Old   September 15, 2008, 04:27
Default Hi Max! The current version
  #28
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Max!

The current version of the icoLagrangianFoam was NOT tested with 1.5 (and I'd be surprised if it worked). A more recent example of a lagrangian solver in 1.5 can be found in $FOAM_TUTORIALS/rhoTurbTwinParcelFoam.

Yes the particles are treated pointwise. If you're looking for a continous approach have a look at $FOAM_SOLVERS/multiphase. Especially twoPhaseEulerFoam

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   October 30, 2008, 11:40
Default Bernhard, Can you send me a
  #29
Member
 
Scott Ripplinger
Join Date: Mar 2009
Posts: 30
Rep Power: 8
sripplinger is on a distinguished road
Bernhard,

Can you send me a copy of that injection through patch code? My email is sripplinger@earthnet.us
sripplinger is offline   Reply With Quote

Old   October 30, 2008, 18:13
Default Hi Scott! You are referring
  #30
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Scott!

You are referring to the posting from august 2006, are you? I can send it to you, but it would be an exercise in archeology because since then the Lagrangian-particle system went through some serious evolutions (which changed things for the better, I must add) especially in terms of the run-time-selectable injectors, I'd be surprised if the code compiles/works with recent versions of the particle system (and I currently haven't got the time to modify/test it), but if you insist, I can send it to you

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   October 31, 2008, 11:28
Default Bernhard, I am planning on
  #31
Member
 
Scott Ripplinger
Join Date: Mar 2009
Posts: 30
Rep Power: 8
sripplinger is on a distinguished road
Bernhard,

I am planning on making use of the Solid Particle class in OF v1.5. Right now, however, I already have a solver that Jeff Allen put together based on your icoLagrangianFoam that runs with v1.3. So I could still use the pertinent pieces of code for patch injection, if you can find them. If it is an inconvenience I will make do with something else.
sripplinger is offline   Reply With Quote

Old   November 3, 2008, 08:09
Default Hi Scott! I sent you an EMa
  #32
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Hi Scott!

I sent you an EMail

Bernhard
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   November 29, 2013, 15:50
Default volPointInterpolation error
  #33
New Member
 
Join Date: Jul 2013
Posts: 2
Rep Power: 0
Ramsky is on a distinguished road
Hello all,

I am using OF 2.1.x to simulate spray using sprayFoam. The fuel injected is a new fuel that has been added to the library.

The solver encounters "volPointInterpolation" error (see below) after about 5ms of flow time. However, when an existing fuel from the liquidProperties library is injected, there are no errors.

The error also refers to GeometricField and volMesh. But I'm certain that there are no issues with the snappy mesh (checkMesh -allGeometry shows no errors and face tets are OK).

Thank you !



[8] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
[8] #1 Foam::sigSegv::sigHandler(int) in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libOpenFOAM.so"
[8] #2
[8] at sigaction.c:0
[8] #3 void Foam::volPointInterpolation::interpolateInternalFi eld<Foam::Vector<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam:ointPatchField, Foam:ointMesh>&) const in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libfiniteVolume.so"
[8] #4 Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam:ointPatchField, Foam:ointMesh> > Foam::volPointInterpolation::interpolate<Foam::Vec tor<double> >(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) const in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libfiniteVolume.so"
[8] #5 Foam::interpolationCellPoint<Foam::Vector<double> >::interpolationCellPoint(Foam::GeometricField<Foa m::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libfiniteVolume.so"
[8] #6 Foam::interpolation<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::interpola tionCellPoint<Foam::Vector<double> > >::New(Foam::GeometricField<Foam::Vector<double> , Foam::fvPatchField, Foam::volMesh> const&) in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/lib/libfiniteVolume.so"
[8] #7 Foam::interpolation<Foam::Vector<double> >::New(Foam::word const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> const&) in "/home/cfduser/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64Gcc46DPOpt/bin/sprayFoam"
Ramsky is offline   Reply With Quote

Old   December 2, 2013, 03:39
Default
  #34
Senior Member
 
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 11
haakon will become famous soon enough
I think much of the particle functionality discussed in this post is out of date. A lot have happened in OpenFOAM on this field since then, especially from the 2.0-version until now. I suggest that you create a new thread for your problem, and try to describe your problem in more details. If you can, please also present examples and code snipplets. And you should always describe what you have done to fix the issue!
haakon 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
Normal particle velocity in particle deposition! Prashanth Fluent Multiphase 1 May 29, 2013 22:30
Add icoLagrangianFoam to FoamX dyxu OpenFOAM Running, Solving & CFD 10 June 11, 2008 09:53
IcoLagrangianFoam problem in contiuation run amp field reading from input stream gschaider OpenFOAM Running, Solving & CFD 2 May 27, 2008 03:45
Add icoLagrangianFoam to FoamX dyxu OpenFOAM Pre-Processing 1 August 13, 2007 04:21
Forces ? Michal Marcinkowski Phoenics 0 January 7, 2001 08:12


All times are GMT -4. The time now is 06:33.