[MensaFoam]

Hello all,

We (two students) try to combine SolidParticle with a simple incompressible solver. We also inserted an injection method by following the tutorial of Anton Persson (Chalmers) http://www.tfd.chalmers.se/~hani/kur...tonPersson.pdf
Complilation and simulation seems to be successful (no errors). Although we had some problems to define the arguments for creating the new particle.

Our Problem now:
Its seems to simulate successfully (increasing calculationtime for each time step, caused by more particles in the domain), but we can't visualize the injected particles, only particles which are predefined in "0"-Folder can be seen.

Our code:
part of injectSolidParticleCloud class:

void Foam::injectSolidParticleCloud::move(const dimensionedVector& g)
{
    // Injector 1
    //Set injection position (z=0 if 2d)
    scalar posy=9e-2;
    scalar posz=0.05;
    scalar posx=1e-2;
    vector pos = vector(posx,posy,posz);
    //Set initial velocity vector
    vector vel=vector(0,0,0);
    //Particle diameter
    scalar d = 2.0e-3;
    // Find cell at specified injection position and add particle here
    label cellI=mesh_.findCell(pos);
    mesh_.findCellFacePt(pos, cellI, 1, 1);
    solidParticle* ptr= new solidParticle
            (
                mesh_,        //const polyMesh& mesh,
                pos,        //const vector& position,
                cellI,        //const label celli,
                1,        //const label tetFacei,
                1,        //const label tetPtI,
                d,        //const scalar d,
                vel        //const vector& U
            );
    Cloud<injectorSolidParticle>::addParticle(ptr);


    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(td, mesh_.time().deltaTValue());
}

Thanks a lot in advance!

[vishwesh]

Originally Posted by MensaFoam

Hello all,

We (two students) try to combine SolidParticle with a simple incompressible solver. We also inserted an injection method by following the tutorial of Anton Persson (Chalmers) http://www.tfd.chalmers.se/~hani/kur...tonPersson.pdf
Complilation and simulation seems to be successful (no errors). Although we had some problems to define the arguments for creating the new particle.

Our Problem now:
Its seems to simulate successfully (increasing calculationtime for each time step, caused by more particles in the domain), but we can't visualize the injected particles, only particles which are predefined in "0"-Folder can be seen.

Our code:
part of injectSolidParticleCloud class:

void Foam::injectSolidParticleCloud::move(const dimensionedVector& g)
{
    // Injector 1
    //Set injection position (z=0 if 2d)
    scalar posy=9e-2;
    scalar posz=0.05;
    scalar posx=1e-2;
    vector pos = vector(posx,posy,posz);
    //Set initial velocity vector
    vector vel=vector(0,0,0);
    //Particle diameter
    scalar d = 2.0e-3;
    // Find cell at specified injection position and add particle here
    label cellI=mesh_.findCell(pos);
    mesh_.findCellFacePt(pos, cellI, 1, 1);
    solidParticle* ptr= new solidParticle
            (
                mesh_,        //const polyMesh& mesh,
                pos,        //const vector& position,
                cellI,        //const label celli,
                1,        //const label tetFacei,
                1,        //const label tetPtI,
                d,        //const scalar d,
                vel        //const vector& U
            );
    Cloud<injectorSolidParticle>::addParticle(ptr);


    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(td, mesh_.time().deltaTValue());
}

Thanks a lot in advance!

Hi!
I am also having similar problem. Did you find any solution yet?