[tgvosk]

Hello,

I am trying to scale a volScalarField (myField) by cell volume. Based on other threads on this topic, I have the following code (which works)

    volScalarField cellVolume
    (
        IOobject
        (
            "cellVolume",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        mesh,
        dimensionedScalar("zero", dimensionSet(0, 3, 0, 0, 0, 0, 0), 0.0)
    );

    cellVolume.internalField() = mesh.V();

    myField = myField * cellVolume;

What I would like to know is whether there is a way of doing this directly with field operations (without creating an intermediate volScalarField to store the cell volume or looping through all cells and multiplying cell-by-cell).

Any suggestions?

[stevenvanharen]

Use this:

T.internalField() *=  mesh.V();

[tgvosk]

Steven,

Thanks, that's what I was looking for! Both of the following work as expected.

 myField = mag(fvc::laplacian(T));
 myField.internalField() *= mesh.V();

myField.internalField() = mag(fvc::laplacian(T)) * mesh.V();