[vonboett]

I would use salome to create a tet mesh and convert to OpenFOAM (some nice tutorials on youtube).
However, I use hex mesh and snappyHexMesh so I can't help you here.

Good Luck!

[vishwesh]

Originally Posted by vonboett

I would use salome to create a tet mesh and convert to OpenFOAM (some nice tutorials on youtube).
However, I use hex mesh and snappyHexMesh so I can't help you here.

Good Luck!

Thanks for the help! That's what I am planning to do

[vishwesh]

Hi!

I was able to make a tet mesh using SALOME. I have another doubt regarding the solver technique proposed by Aurelia in her thesis.

The inject function that we are using, can it be used multiple times to inject multiple particles at different time instances? Also what happens to the particles that have been injected earlier...

[vishwesh]

Hi!

Thanks to the help, I was able to make changes to the solver and now it's working fine. I have one problem where I am stuck, following Aurelia's work and adding rhoInterp, UInterp, etc. in solidParticleCloud.C file (move function), when I run the solver, it gets stuck as soon as it reaches the interpolation part. I don't know why that is happening. I left the solver working for a whole day on 32 processors but it didn't proceed further than that.

Here is the relevant solidParticleCloud.C section:

    Info<<"Function called"<<endl;
    const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho");
    Info<<"rho: "<<rho<<endl;
    const volVectorField& U = mesh_.lookupObject<const volVectorField>("U");
    const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu");
    const volScalarField& alphaW = mesh_.lookupObject<const volScalarField>("alpha.water");
    const volScalarField& p = mesh_.lookupObject<const volScalarField>("p");

    Info<<"lookup object worked"<<endl;

    Info<<"rho interpolation"<<endl;
    interpolationCellPoint<scalar> rhoInterp(rho);
    Info<<"U interpolation"<<endl;
    interpolationCellPoint<vector> UInterp(U);
    Info<<"nu interpolation"<<endl;
    interpolationCellPoint<scalar> nuInterp(nu);
    Info<<"alpha.water interpolation"<<endl;
    interpolationCell<scalar> alphaWInterp(alphaW);
    Info<<"p interpolation"<<endl;
    interpolationCellPoint<scalar> pInterp(p);

    Info<<"Interpolation complete"<<endl;

And the results:

Injecting particle with diameter: 0.0301191m, center position: (0.153843 0.0928162 0.102841)m, velocity: (0 -2 0)m/s
Function called
rho: dimensions      [1 -3 0 0 0 0 0];

internalField   uniform 1;

boundaryField
{
    inlet
    {
        type            symmetry;
    }
    outlet
    {
        type            symmetry;
    }
    walls
    {
        type            symmetry;
    }
    procBoundary0to1
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to2
    {
        type            processor;
        value           uniform 1;
    }
    procBoundary0to3
    {
        type            processor;
        value           uniform 1;
    }
}

lookup object worked
rho interpolation

[yuanlee2011]

Originally Posted by alexlupo

I am not sure if I understood your question properly, but the solidParticle class should not give any problems relating phase migration. The issues you can experience (and by the way that I am experiencing) are related to the particles going out of the boundaries, thus showing unexpected behavior, or simulations crashing in parallel due to the injection of particles in more than one processor boundary.
If your question is related to particles moving from one phase to another, I guess that could eventually happen, but from my point of view that could also have a physical explanation, taking into account the difference of density between both particles and fluids.
I guess if both phases are clearly separated, you could just create a boundary in between so that you would not allow the particles to move from phase to the other.
Let me know if this was what you intend to do and if not, could you explain it further?
Cheers,

The simulation could perform smoothly when run in serial， while crashed in parallel，what‘s the reasons？

[soroush.amd]

Hello all,

I am using OpenFoam 4.1 and have coupled solidParticle class with icoFoam for particle tracking. It seems to work fine except for two issues that I just came across.
1. I have noticed that sometimes the interpolation looks unreasonable and trapped in a very long loop. this only happens when my particle is on a face. something like below:

Time = 3.028 trackTime = 0.001 tEnd = 0.001 steptFraction() = 0
in solidParticle, position = (0.0381269 0.095 0.000327212), dt = 0.001, U = (0.685765 0.0165633 0)
U_ = (0.685765 0.0165633 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999879 steptFraction() = 0.000120856
in solidParticle, position = (0.038127 0.095 0.000327212), dt = 0.000999879, U = (0.685766 0.0165632 0)
U_ = (0 -0.0165632 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999879 steptFraction() = 0.000120856
in solidParticle, position = (0.038127 0.095 0.000327281), dt = 0.000999879, U = (1 0 0)
U_ = (1 0 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999879 steptFraction() = 0.000120856
in solidParticle, position = (0.038127 0.095 0.000327315), dt = 0.000999879, U = (0.685766 0.0165632 0)
U_ = (0.685766 0.0165632 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999124 steptFraction() = 0.000875538
in solidParticle, position = (0.0381275 0.095 0.000327315), dt = 0.000999124, U = (0.685768 0.0165623 0)
U_ = (0 -0.0165623 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999124 steptFraction() = 0.000875538
in solidParticle, position = (0.0381275 0.095 0.000327383), dt = 0.000999124, U = (1 0 0)
U_ = (1 0 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.000999124 steptFraction() = 0.000875538
in solidParticle, position = (0.0381275 0.095 0.000327417), dt = 0.000999124, U = (0.685768 0.0165623 0)
U_ = (0.685768 0.0165623 0)
1moving in solidParticle.C
Time = 3.028 trackTime = 0.001 tEnd = 0.00099837 steptFraction() = 0.00163026
in solidParticle, position = (0.038128 0.095 0.000327417), dt = 0.00099837, U = (0.685771 0.0165615 0)
........ and it goes on for much longer

It is a cavity case with structured grid, and (0.0381269 0.095 0.000327212) falls on a face. One thing to note is that the interpolated velocity does not make sense at all as it fluctuates between (0.685768 0.0165623 0) and (1 0 0).

2. The second issue that I have noticed is that although my case is 2D and there is no z-velocity component, my final particle position depends on its initial z-coordinate which is not physical. I have noticed that the code itself change the particle z-coordinate sometimes but not sure why this happens, its definitely not due to z-velocity component as its always zero and I have confirmed it. So the solidParticle library is doing something under the hood which I do not know and whatever it is it does not seem correct to me.

I was wondering if anyone has faced such issues? I have heard in other threads that solidParticle library has bugs and should not be used, should I use something else such as intermediate library instead?

Thank you for your help,

[vonboett]

have you tried in 3D?
I recommend DPMFoam as a good Lagrangian coupled solver, or IcoUncoupledKinematicParcelFoam if you only need particles as trackers.
If your particle is larger then the cell size at the face, your particle will stick to the wall, it will "sink into it" and cause pressure oscillatios.

[soroush.amd]

Thank you for your reply. Yes I have also tried it in 3D, although I did not try to reproduce this error because I am interested in 2D problems only. I found solidParticle library somehow immature because in intermediate library there are some functionally specifically for 2D to restrain the particle within the mesh and its interpolated velocity, also to make sure particle are injected via one processor only which is a problem in solidParticle library if your injection occurs on a processor patch. So I have moved to intermediate library as its more mature and reliable.

Originally Posted by vonboett

have you tried in 3D?
I recommend DPMFoam as a good Lagrangian coupled solver, or IcoUncoupledKinematicParcelFoam if you only need particles as trackers.
If your particle is larger then the cell size at the face, your particle will stick to the wall, it will "sink into it" and cause pressure oscillatios.

[neilk]

Originally Posted by houwy

Hi! I still have some problem with it.
code:

void Foam::solidParticleCloud::inject(solidParticle::trackingData &td)
{
  label cellI = 1;
  label tetFaceI = 1;
  label tetPtI = 1;
  mesh_.findCellFacePt(td.cloud().posP1_, cellI, tetFaceI, tetPtI);
  solidParticle* ptr1=new solidParticle(*this,td.cloud().posP1_,cellI,tetFaceI, tetPtI,
  td.cloud().dP1_,td.cloud().UP1_);
  Cloud<solidParticle>::addParticle(ptr1);
  mesh_.findCellFacePt(td.cloud().posP1_, cellI, tetFaceI, tetPtI);
  solidParticle* ptr2=new solidParticle(*this,td.cloud().posP2_,cellI,tetFaceI, tetPtI,
  td.cloud().dP2_,td.cloud().UP2_);
  Cloud<solidParticle>::addParticle(ptr2);
}

when I compile it , There are still some errors.

solidParticleCloud.C: In member function ‘void Foam::solidParticleCloud::inject(Foam::solidParticle::trackingData&)’:
solidParticleCloud.C:81:34: error: no matching function for call to ‘Foam::solidParticle::solidParticle(Foam::solidParticleCloud&, Foam::vector&, Foam::label&, Foam::label&, Foam::label&, Foam::scalar&, Foam::vector&)’
solidParticleCloud.C:81:34: note: candidates are:
solidParticle.H:135:9: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, Foam::Istream&, bool)
solidParticle.H:135:9: note:   candidate expects 3 arguments, 7 provided
solidParticleI.H:47:8: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, const vector&, Foam::label, Foam::label, Foam::label, Foam::scalar, const vector&)
solidParticleI.H:47:8: note:   no known conversion for argument 1 from ‘Foam::solidParticleCloud’ to ‘const Foam::polyMesh&’
solidParticle.H:58:7: note: Foam::solidParticle::solidParticle(const Foam::solidParticle&)
solidParticle.H:58:7: note:   candidate expects 1 argument, 7 provided
solidParticleCloud.C:85:34: error: no matching function for call to ‘Foam::solidParticle::solidParticle(Foam::solidParticleCloud&, Foam::vector&, Foam::label&, Foam::label&, Foam::label&, Foam::scalar&, Foam::vector&)’
solidParticleCloud.C:85:34: note: candidates are:
solidParticle.H:135:9: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, Foam::Istream&, bool)
solidParticle.H:135:9: note:   candidate expects 3 arguments, 7 provided
solidParticleI.H:47:8: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, const vector&, Foam::label, Foam::label, Foam::label, Foam::scalar, const vector&)
solidParticleI.H:47:8: note:   no known conversion for argument 1 from ‘Foam::solidParticleCloud’ to ‘const Foam::polyMesh&’
solidParticle.H:58:7: note: Foam::solidParticle::solidParticle(const Foam::solidParticle&)
solidParticle.H:58:7: note:   candidate expects 1 argument, 7 provided 

Hello,

Were you able to solve this problem? I am trying to compile the same solver and it's giving me the same error.

Thank you

Nilay

[neilk]

Originally Posted by houwy

Hi! I still have some problem with it.
code:

void Foam::solidParticleCloud::inject(solidParticle::trackingData &td)
{
  label cellI = 1;
  label tetFaceI = 1;
  label tetPtI = 1;
  mesh_.findCellFacePt(td.cloud().posP1_, cellI, tetFaceI, tetPtI);
  solidParticle* ptr1=new solidParticle(*this,td.cloud().posP1_,cellI,tetFaceI, tetPtI,
  td.cloud().dP1_,td.cloud().UP1_);
  Cloud<solidParticle>::addParticle(ptr1);
  mesh_.findCellFacePt(td.cloud().posP1_, cellI, tetFaceI, tetPtI);
  solidParticle* ptr2=new solidParticle(*this,td.cloud().posP2_,cellI,tetFaceI, tetPtI,
  td.cloud().dP2_,td.cloud().UP2_);
  Cloud<solidParticle>::addParticle(ptr2);
}

when I compile it , There are still some errors.

solidParticleCloud.C: In member function ‘void Foam::solidParticleCloud::inject(Foam::solidParticle::trackingData&)’:
solidParticleCloud.C:81:34: error: no matching function for call to ‘Foam::solidParticle::solidParticle(Foam::solidParticleCloud&, Foam::vector&, Foam::label&, Foam::label&, Foam::label&, Foam::scalar&, Foam::vector&)’
solidParticleCloud.C:81:34: note: candidates are:
solidParticle.H:135:9: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, Foam::Istream&, bool)
solidParticle.H:135:9: note:   candidate expects 3 arguments, 7 provided
solidParticleI.H:47:8: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, const vector&, Foam::label, Foam::label, Foam::label, Foam::scalar, const vector&)
solidParticleI.H:47:8: note:   no known conversion for argument 1 from ‘Foam::solidParticleCloud’ to ‘const Foam::polyMesh&’
solidParticle.H:58:7: note: Foam::solidParticle::solidParticle(const Foam::solidParticle&)
solidParticle.H:58:7: note:   candidate expects 1 argument, 7 provided
solidParticleCloud.C:85:34: error: no matching function for call to ‘Foam::solidParticle::solidParticle(Foam::solidParticleCloud&, Foam::vector&, Foam::label&, Foam::label&, Foam::label&, Foam::scalar&, Foam::vector&)’
solidParticleCloud.C:85:34: note: candidates are:
solidParticle.H:135:9: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, Foam::Istream&, bool)
solidParticle.H:135:9: note:   candidate expects 3 arguments, 7 provided
solidParticleI.H:47:8: note: Foam::solidParticle::solidParticle(const Foam::polyMesh&, const vector&, Foam::label, Foam::label, Foam::label, Foam::scalar, const vector&)
solidParticleI.H:47:8: note:   no known conversion for argument 1 from ‘Foam::solidParticleCloud’ to ‘const Foam::polyMesh&’
solidParticle.H:58:7: note: Foam::solidParticle::solidParticle(const Foam::solidParticle&)
solidParticle.H:58:7: note:   candidate expects 1 argument, 7 provided 

Hello,
Were you able to solve this problem? I am trying to compile the same solver and I am stuck at this point. Could you please help me out?

Thank you

Nilay

[thintn222]

Thank you so much for this post !

[kk415]

hello All


I am also working in coupling of interFoam with lagrangian particle. I still don't understand what is the meaning of tetFacei and tetPtI. If I need to define a particle I need its position, diameter and velocity. what is importance of tetFaceI and tetPtI. If I want to create a particle at position X which lies inside cellI then physically what is tetFaceI and tetPtI ???

[wangziyang]

Originally Posted by alexlupo

Hi there!
I havent been able to validate because I am working now on the same kind of solver but with pimpleFoam.
Anyway, as a starting point, there is a very good tutorial on how to couple both interFoam and the solidParticle class:
http://www.tfd.chalmers.se/~hani/kur...LPT_120911.pdf
The tutorial is for a previous version of OpenFoam so you will have to upgrade it in case you are using 2.2.x.
Take a look and just let me know if you need some help. I will be glad to help you.
Cheers,

Hi Formers!

Referring to the previous paper and guide, i want to insert the solidParticle into the pisoFoam.

The operiations are followed by the Ds. Vallier guide, i used the of3.0.x

when i complie the new solver pisoLPTFoam, it appear an error:

pisoLPTFoam.C: In function ‘int main(int, char**)’:
pisoLPTFoam.C:87:18: error: ‘g’ was not declared in this scope
   particles.move(g);

When i complie the particle into the InterFoam, it works well, but fail for the pisoFoam. The error says the 'g' in the move function is not declaration.

How can i modify the code?

Thanks in advances!

ziyang

[avinashv20]

Hello foamers!


While coupling the solidParticle library with and interFoam based solver I am getting following error:


solidParticleCloud.C:75:27: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
mesh_.findCellFacePt(td.cloud().posP_, cellI, tetFaceI, tetPtI);
^~~~~
solidParticleCloud.C:76:50: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
solidParticle* ptr=new solidParticle(mesh_, td.cloud().posP_, cellI, tetFaceI, tetPtI, td.cloud().dP_, td.cloud().UP_);
^~~~~
solidParticleCloud.C:76:93: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
solidParticle(mesh_, td.cloud().posP_, cellI, tetFaceI, tetPtI, td.cloud().dP_, td.cloud().UP_);
^~~~~
solidParticleCloud.C:76:109: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
esh_, td.cloud().posP_, cellI, tetFaceI, tetPtI, td.cloud().dP_, td.cloud().UP_);
^~~~~
solidParticleCloud.C: In member function ‘void Foam::solidParticleCloud::move(const dimensionedVector&)’:
solidParticleCloud.C:95:30: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
if(mesh_.time().value()> td.cloud().tInjStart_&&
^~~~~
solidParticleCloud.C:96:28: error: ‘class Foam::solidParticle::trackingData’ has no member named ‘cloud’
mesh_.time().value()< td.cloud().tInjEnd_)





my solidParticleCloud.C file is as follows:


void Foam::solidParticleCloud::inject(solidParticle::tr ackingData& td)
{
label cellI = 1;
label tetFaceI = 1;
label tetPtI = 1;
mesh_.findCellFacePt(td.cloud().posP_, cellI, tetFaceI, tetPtI);
solidParticle* ptr=new solidParticle(mesh_, td.cloud().posP_, cellI, tetFaceI, tetPtI, td.cloud().dP_, td.cloud().UP_);
Cloud<solidParticle>::addParticle(ptr);

}

void Foam::solidParticleCloud::move(const dimensionedVector& g)
{
const volScalarField& rho = mesh_.lookupObject<const volScalarField>("rho");
const volVectorField& U = mesh_.lookupObject<const volVectorField>("U");
const volScalarField& nu = mesh_.lookupObject<const volScalarField>("nu");

interpolationCellPoint<scalar> rhoInterp(rho);
interpolationCellPoint<vector> UInterp(U);
interpolationCellPoint<scalar> nuInterp(nu);

solidParticle::trackingData
td(*this, rhoInterp, UInterp, nuInterp, g.value());

Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue());
if(mesh_.time().value()> td.cloud().tInjStart_&&
mesh_.time().value()< td.cloud().tInjEnd_)
{
this->inject(td);
}
}


Please help me resolve this.


I am using OpenFOAM-in-Box-18.06v1

[wind_]

Originally Posted by wangziyang

Hi Formers!

Referring to the previous paper and guide, i want to insert the solidParticle into the pisoFoam.

The operiations are followed by the Ds. Vallier guide, i used the of3.0.x

when i complie the new solver pisoLPTFoam, it appear an error:

pisoLPTFoam.C: In function ‘int main(int, char**)’:
pisoLPTFoam.C:87:18: error: ‘g’ was not declared in this scope
   particles.move(g);

When i complie the particle into the InterFoam, it works well, but fail for the pisoFoam. The error says the 'g' in the move function is not declaration.

How can i modify the code?

Thanks in advances!

ziyang

Hi, ziyang
I'm facing the same problem now, just replaced pisoFoam with pimpleFoam
Thanks!