[alan_fergus]

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

[gschaider]

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

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

[gschaider]

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)

[francesco_b]

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]

I've succeed in fixing the problem, sorry for my stupid question

[maxe]

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

[gschaider]

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

[sripplinger]

Bernhard,

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

[gschaider]

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

[sripplinger]

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.

[gschaider]

Hi Scott!

I sent you an EMail

Bernhard

[Ramsky]

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"

[haakon]

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!