[cmigueis]

Hi, chpjz0391!

The addition of the collisionCloud into the coalChemistry Foam was accomplished!
You can analize the code in the reply that I've sent in december, 28, in this thread!

About the alphac, it's only required when your have a dense particulate flow - the magical number (I don't know how it was determined) is about 10 or 12% (or greater) of particle fraction (volumetric) into your domain.
I mean:

Vparticles/Vcell >= 0,12.

I hope this was helpful!

Regards from Brazil!

[chpjz0391]

Originally Posted by cmigueis

Hi, chpjz0391!

The addition of the collisionCloud into the coalChemistry Foam was accomplished!
You can analize the code in the reply that I've sent in december, 28, in this thread!

About the alphac, it's only required when your have a dense particulate flow - the magical number (I don't know how it was determined) is about 10 or 12% (or greater) of particle fraction (volumetric) into your domain.
I mean:

Vparticles/Vcell >= 0,12.

I hope this was helpful!

Regards from Brazil!

HI César,
Thanks for your reply.
I modified the lib and compiled it successfully. However, when I complied my solver I got the following error.

5): undefined reference to `Foam::ThermoParcel<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > >::debug'/home/ping/OpenFOAM/root-3.0.1/platforms/linux64GccDPInt32Opt/lib/libmycollision.so: undefined reference to `Foam::WallModel<Foam::CollidingCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > >::dictionaryConstructorTablePtr_'/home/ping/OpenFOAM/root-3.0.1/platforms/linux64GccDPInt32Opt/lib/libmycollision.so: undefined reference to `Foam::PairModel<Foam::CollidingCloud<Foam::KinematicCloud<Foam::Cloud<Foam::ReactingMultiphaseParcel<Foam::ReactingParcel<Foam::ThermoParcel<Foam::CollidingParcel<Foam::KinematicParcel<Foam::particle> > > > > > > > >::dictionaryConstructorTablePtr_'collect2: error: ld returned 1 exit status

Then I modified "makeBasicReactingMultiphaseParcelSubmodels.C" as well as "defineBasicReactingMultiphaseParcel.C"like yours.


makeBasicReactingMultiphaseParcelSubmodels.C

"ttypedef basicReactingMultiphaseCloud::collidingCloudType collidingParcelType;
makeParcelCollisionModels(collidingParcelType);
"

defineBasicReactingMultiphaseParcel.C

     typedef    Foam::ThermoParcel<CollidingParcel<KinematicParcel<particle> > > thermoParcel;
    typedef    Foam::ReactingParcel<ThermoParcel<CollidingParcel<KinematicParcel<particle> > > > reactingParcel;

    defineTemplateTypeNameAndDebug(thermoParcel, 0);
   defineTemplateTypeNameAndDebug(reactingParcel, 0);

However, this time I cannot even compile the lib .

Do you have any suggestion? I am using OpenFOAM 3.0.1 . Is the version problem? The attachment is my code.

Kind Regards,

Ping

mycollision.zip

[amin_jalalian]

Hi,
I think you must write "- alpha*rho.dimensionedInternalField()*g
==
"

[amin_jalalian]

Originally Posted by cmigueis

Hi Shuai! Did you have any progress?

I started to modify the coalChemistryFoam code to include the void fraction into the transport equations. In the coalChemistryFoam.c file I included the following code:

// Update continuous phase volume fraction field
        alpha = max(1.0 - coalCloud.theta() - limestoneCloud.theta(), 0.01);
        alpha.correctBoundaryConditions();
        alphaf = fvc::interpolate(alpha);
        alphaPhi = alphaf*phi;

that I adapted from the DPMFoam solver. The cloudName.theta() return the volumetric fraction of particles in the cell, so the alpha is the void fraction of the cell.

So, I modified the UEqn.h file into this (I'm gonna show only the modified part):

fvVectorMatrix UEqn
    (
        fvm::ddt(rho*alpha, U)
      + fvm::div(alphaPhi, U)
      + alpha*fvc::grad(p)
      + alpha*turbulence->divDevRhoReff(U)
      + alpha*rho.dimensionedInternalField()*g
     ==
        coalParcels.SU(U)
      + limestoneParcels.SU(U)
      + fvOptions(rho, U)
    );

The continuity equation looks like:

    {
    fvScalarMatrix rhoEqn
    (
        fvm::ddt(alpha,rho)
      + fvc::div(alpha,phi)
      ==
        coalParcels.Srho(rho)
    limestoneParcels.Srho(rho)
      + fvOptions(rho)
    );

    fvOptions.constrain(rhoEqn);

    rhoEqn.solve();

    fvOptions.correct(rho);
}

    );

Now I'm trying tho modify the energy equation, but I'm having a hard time trying to understand the code (I've already modified the temporal and the convective terms of the internal energy and the kinetic energy).

{
    volScalarField& he = thermo.he(); 

    fvScalarMatrix EEqn
    (
        fvm::ddt(alpha*rho, he) + mvConvection->fvmDiv(alphaPhi, he)
      + fvc::ddt(alpha*rho, K) + fvc::div(alphaPhi, K)
      + (
            he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt
        )
      - fvm::laplacian(turbulence->alphaEff(), he)
     ==
        combustion->Sh()
      + coalParcels.Sh(he)
      + limestoneParcels.Sh(he)
      + radiation->Sh(thermo)
      + fvOptions(rho, he)
    );

Seems to me that this code should represent the Energy Equation in the for of the equation 12 in the article about the energy equation in the CFD DIrect ( http://cfd.direct/openfoam/energy-equation/ )



I'm stucked in this part:

+ (
            he.name() == "e"
          ? fvc::div
            (
                fvc::absolute(phi/fvc::interpolate(rho), U),
                p,
                "div(phiv,p)"
            )
          : -dpdt
        )

Can someone explain the code above? Which part of Eq. 12 (assuming that is the form of Eq 12 indeed) it represents?

Thanks in advance to everybody and regards from Brazil!

[ Moderator note: Moved from this thread: http://www.cfd-online.com/Forums/ope...ification.html - Reason: consolidation of the ongoing development ]

Hi,
I think you must write
"- alpha*rho.dimensionedInternalField()*g
== "
in UEqn.H