[aliyah.]

Hello

I am trying to use codestream to create several oval shapes in my alpha domains. These ovals are placed throughout the domain*at various x coordinates, each of which is determined by a for loop with j label, and here's my 0/alpha file.

dimensions      [0 0 0 0 0 0 0];


internalField    #codeStream
{
	codeInclude
	#{
		#include "fvCFD.H"
	#};
        
	codeOptions
	#{
		-I$(LIB_SRC)/finiteVolume/lnInclude \
		-I$(LIB_SRC)/meshTools/lnInclude
	#};

	codeLibs
	#{
		-lmeshTools \
		-lfiniteVolume
	#};
	code
	#{
    
    const IOdictionary& d = static_cast<const IOdictionary&>(dict);
	const fvMesh& mesh = refCast<const fvMesh>(d.db());
		
   scalarField alpha(mesh.nCells(), 0);

    for (label j=0; j <= 10; j++) 
    {

    scalar lambda =  2.29922953e-3;
    scalar he = j*lambda;
    scalar ke = 0.5*lambda;
    scalar ae = 0.2*lambda;
    scalar be = 0.1*lambda;

		forAll(alpha, i)
		{
			const scalar x = mesh.C()[i][0];
			const scalar y = mesh.C()[i][1];

			if (  pow(y-ke,2) <= ((1 - pow(x-he,2)/pow(ae,2) )*pow(be,2)) )
			{
				alpha[i] = 1;
			}
		}
		alpha.writeEntry("", os);
    }
	#}
    ;
};


boundaryField
{
    up
    {
        type            inletOutlet;
        inletValue      uniform 1;
        value           uniform 0;
    }
    "wall.*"
    {
        type            cyclic;
    }
    down
    {
        type            fixedValue;
        value           uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
}
}

However, only one oval is created (the one with j equal to 0), and the remaining ovals are not formed.

What is the problem here and what I should do? Thank you!

[aliyah.]

The problem with this codestream was the location of writeEntry. It should be right before #}.

So it will be like

internalField    #codeStream
{
	codeInclude
	#{
		#include "fvCFD.H"
	#};
        
	codeOptions
	#{
		-I$(LIB_SRC)/finiteVolume/lnInclude \
		-I$(LIB_SRC)/meshTools/lnInclude
	#};

	codeLibs
	#{
		-lmeshTools \
		-lfiniteVolume
	#};
	code
	#{
    
    const IOdictionary& d = static_cast<const IOdictionary&>(dict);
	const fvMesh& mesh = refCast<const fvMesh>(d.db());
		
   scalarField alpha(mesh.nCells(), 0);

    for (label j=0; j <= 10; j++) 
    {

    scalar lambda =  2.29922953e-3;
    scalar he = j*lambda;
    scalar ke = 0.5*lambda;
    scalar ae = 0.2*lambda;
    scalar be = 0.1*lambda;

		forAll(alpha, i)
		{
			const scalar x = mesh.C()[i][0];
			const scalar y = mesh.C()[i][1];

			if (  pow(y-ke,2) <= ((1 - pow(x-he,2)/pow(ae,2) )*pow(be,2)) )
			{
				alpha[i] = 1;
			}
		}
    }
		alpha.writeEntry("", os);
	#}
    ;
};