[IC01]

Hi
I have the problem with vectorCodedSource; in fvOption file
I'm modifying incompressible/simpleFoam/rotorDisk tutorial. The target is to have either polynomial or tabular velocity in the RotatingZone. Please let me know if there is a simpler way
Right now I'm trying to put velocity in the y-direction that corresponds to the y coordinate of cell. It gives me not feasible values

Thank you

{
 

    type            vectorCodedSource;
    active          true;
    selectionMode   cellZone;
    cellZone        rotatingZone;

    fields          (U);
    name            codedSource;
    codeAddSup
    #{
        const Time& time = mesh().time();
	const scalarField& V = mesh_.V();
        vectorField& USource = eqn.source();
       // const scalarField& cellx = mesh_.C().component(0) ; 
        const scalarField& celly = mesh_.C().component(1) ;
      //  const scalarField& cellz = mesh_.C().component(2) ;
      

        // Start time
        const scalar startTime = 20.0;  

 	if (time.value() > startTime)

    
        forAll(celly, i)
        {
            scalar v1 = 0 ;
	    scalar v2 = celly[i] ;
	    scalar v3 = 0 ;


            USource[i][0] +=celly[i]; //celly[i];
        }
    #};

}

[Mahmoud Abbaszadeh]

Originally Posted by IC01

Hi
I have the problem with vectorCodedSource; in fvOption file
I'm modifying incompressible/simpleFoam/rotorDisk tutorial. The target is to have either polynomial or tabular velocity in the RotatingZone. Please let me know if there is a simpler way
Right now I'm trying to put velocity in the y-direction that corresponds to the y coordinate of cell. It gives me not feasible values

Thank you

{
 

    type            vectorCodedSource;
    active          true;
    selectionMode   cellZone;
    cellZone        rotatingZone;

    fields          (U);
    name            codedSource;
    codeAddSup
    #{
        const Time& time = mesh().time();
	const scalarField& V = mesh_.V();
        vectorField& USource = eqn.source();
       // const scalarField& cellx = mesh_.C().component(0) ; 
        const scalarField& celly = mesh_.C().component(1) ;
      //  const scalarField& cellz = mesh_.C().component(2) ;
      

        // Start time
        const scalar startTime = 20.0;  

 	if (time.value() > startTime)

    
        forAll(celly, i)
        {
            scalar v1 = 0 ;
	    scalar v2 = celly[i] ;
	    scalar v3 = 0 ;


            USource[i][0] +=celly[i]; //celly[i];
        }
    #};

}

Hi,

Did you manage to resolve this issue?

[Haitham Osman CFD]

I hope your issue was solved.

If NOT, I think you should multiply your equation by the cell volume. So, the equation should be like

USource[i][0] +=celly[i] *V[i]

Best wishes