[chen112p]

Hello, I am wondering if anyone has help me on this issue.

I am trying to use scalarCodedSource to inject some particles in the domain and transport these particles in a transient simulation. The issue I'm experiencing is the selectionMode setting (cellSet) in scalarCodedSourceCoeffs does not allow me to only apply the source to selected cell set. When using mesh.C() or mesh.V(), the entire mesh got selected.

nox_01
{
    type            scalarTransport;
    libs            ("libsolverFunctionObjects.so");
    field           nox;
	D 				1e-5;
    log no;
    writeControl     writeTime; 
    fvOptions
    {
        noxSource
    {
    type            scalarCodedSource;
	name            noxSource;
    scalarCodedSourceCoeffs
    {
        selectionMode   cellSet;
		cellSet			nox;
        fields          (nox);
        name            noxSource;

        codeInclude
        #{
			#include "OFstream.H"
			#include "fvMesh.H"
			#include "fvMatrices.H"
			#include "DimensionedField.H"
			#include "fvm.H"
        #};

        codeCorrect
        #{
            //Pout<< "**codeCorrect**" << endl;
        #};

        codeAddSup
        #{
	        const Time& time = mesh().time();
            const scalarField& V = mesh_.V();
            const vectorField& C = mesh_.C();
            scalarField& heSource = eqn.source();
            Pout<<V.size();
            forAll(C, i)
            {
                heSource[i] += 10;
            };
        #};
        codeConstrain
        #{
            //Pout<< "**codeConstrain**" << endl;
        #};
    }
}
    }
}

Im also uploading my test case here if anyone wants to do a test run.

[chen112p]

I guess the only way to select cells is still using cellZone like this...

const label zoneID = mesh_.cellZones().findZoneID("nox");
const cellZone& zone = mesh_.cellZones()[zoneID];
const cellZoneMesh& zoneMesh = zone.zoneMesh();
const labelList& cellsZone = zoneMesh[zoneID]; 
volScalarField Flg_
(	
IOobject
  (
	"Flg_",
        mesh_.time().timeName(),
	mesh_
  ),
mesh_,
dimensionedScalar("Flg_", dimensionSet(0,0,-1,0,0,0,0), 0.0)
);
forAll(cellsZone, cI)
{
	Flg_[cellsZone[cI]] += 1.0 * 10;
}

eqn += Flg_;