[saeed jamshidi]

Hi all,
I wanna utilized the codestream for my simulation. To this end, I have prepared "S" dictionary, which calculates the curl of velocity in computaional domain.

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2112                                    |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      S;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField      #codeStream
{
    libs    ("utilityFunctionObjects.so");
	
    codeInclude
    #{
        #include "fvCFD.H"
    #};

	codeOptions
	#{
		-I$(LIB_SRC)/finiteVolume/lnInclude \
		-I$(LIB_SRC)/meshTools/lnInclude
	#};

	codeLibs
	#{
		-lmeshTools  \
		-lfiniteVolume
	#};

	code
	#{
		volScalarField S
		(
			IOobject
			(
				"S",
				mesh().time().timeName(),
				mesh(),
				IOobject::READ_IF_PRESENT,
				IOobject::AUTO_WRITE
			),
			mesh(),
			dimensionedScalar("S", dimensionSet(0, 0, 0, 0, 0, 0, 0), Zero)
		);

		const objectRegistry& db();
		const volVectorField& U = db().lookupObject<volVectorField>("U");
		const fvMesh & mesh = U.mesh();
		
		// Create a scalar field that contains the volume of each cell within the mesh
		const scalarField& v = mesh().V();

		// Calculate the curl of the flow and then take the magnitude of it
		const volScalarField& magCurlU = mag(fvc::curl(U));
		
		forAll(S, celli) 
		{
			// For curl of velocity, use the line below
			S[celli] = pow(v[celli], 1.0/3.0) * magCurlU[celli];
		}

		S.writeEntry("", os);
		
	#};

}

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    
    outlet
    {
        type            zeroGradient;
    }

    top
    {
        type            zeroGradient;
    }
    
    Bottom
    {
        type            zeroGradient;
    }

    cylinder
    {
        type            zeroGradient;
    }

    frontAndBack
    {   
        type            empty;
    }
}

// ************************************************************************* //

but, when i run it in openfoam i see some errors:

/home/myfoam/OpenFOAM/myfoam-v2112/run/tutorials/incompressible/pimpleFoam/LES/cylinder/0/S.#codeStream:47:5: error: ‘mesh’ was not declared in this scope
/home/myfoam/OpenFOAM/myfoam-v2112/run/tutorials/incompressible/pimpleFoam/LES/cylinder/0/S.#codeStream:47:5: note: suggested alternative: ‘cosh’
/usr/lib/openfoam/openfoam2112/wmake/rules/General/transform:34: recipe for target 'Make/linux64GccDPInt32Opt/codeStreamTemplate.o' failed
make: *** [Make/linux64GccDPInt32Opt/codeStreamTemplate.o] Error 1


--> FOAM FATAL IO ERROR: (openfoam-2112 patch=220610)
Failed wmake "dynamicCode/_52a2dbf780a8697296fac2199c6f254fa84ab80e/platforms/linux64GccDPInt32Opt/lib/libcodeStream_52a2dbf780a8697296fac2199c6f254fa84ab80e.so"


file: 0/S at line 17.

    From static void (* Foam::functionEntries::codeStream::getFunction(const Foam::dictionary&, const Foam::dictionary&))(Foam::Ostream&, const Foam::dictionary&)
    in file db/dictionary/functionEntries/codeStream/codeStream.C at line 192.

FOAM exiting

I would appreciate it if you helpe me to solve this problem.
thanks.

[panda007]

Hi Saeed,

I'm facing the same issue when I was trying to use #codestream for customized functions. I came-across your post, I wonder if you got any solution for this issue. Would you mind sharing your solution or any suggestion how I can fix this issue? Thank you so much!!

[geth03]

internalField  #codeStream
{
    code
    #{
        const IOdictionary& d = static_cast<const IOdictionary&>(dict);
        const fvMesh& mesh = refCast<const fvMesh>(d.db());
        scalarField fld(mesh.nCells(), 12.34);
        fld.writeEntry("", os);
    #};

    //! Optional:
    codeInclude
    #{
        #include "fvCFD.H"
    #};

    //! Optional:
    codeOptions
    #{
        -I$(LIB_SRC)/finiteVolume/lnInclude
    #};
};

post #1 does not access the mesh like in the code from the documentation?
https://www.openfoam.com/documentati...odeStream.html