[Nicholas_Langbein]

Hello guys,

i am new to openfoam, i have a variable volScalarField f_sl which is calculated from strain rate multiplied with some constants and then dimensionedScalar (with dimension length) -- now volScalarField is m/s velocity

i want to make this volVectorField with f_sl as Z Component to create fvc::flux() using this volVectorField to use for convection.

Please can someone give an idea how to make volScalarField to volVectorField

i saw volScalarField to volVectorField

but Uvec.internalField().replace(vector::X,avgUx.inter nalField()); dosent work

[nhj]

Hi Nicholas,


off the top of my head, you could try iterating over your volVectorField and assigning the component from f_sl directly to each field element, sort of along the lines of:

// Define your vector field

volVectorField f;


// This loops over each cell of your mesh
forAll(f_sl, celli){
    //Access the z component of your vector field and assign the cell value of the strain rate 
    f[celli].z() = f_sl[celli]; 

}

// ...

Let me know if you figure out a more efficient way of doing this though.

cheers

[dlahaye]

Not sure.

Using the .flux member function of the volScalarField?

[Tobermory]

Originally Posted by Nicholas_Langbein

Hello guys,
i am new to openfoam, i have a variable volScalarField f_sl which is calculated from strain rate multiplied with some constants and then dimensionedScalar (with dimension length) -- now volScalarField is m/s velocity

i want to make this volVectorField with f_sl as Z Component

If I understand right, you want to create a volVectorField where the z component is the value of the scalar, and the x/y components are zero? If so, then that's really easy - try something like:

        vector v0(0,0,1);
        volVectorField v1
                (
                   IOobject ("v1", mesh().time().name(), mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ),
                   mesh(),
                   v0
                );
        v1 *= f_sl;