[rsaintlege]

Hello,

I'm trying to find a function to multiply a field by another but I can't find what I want.
However I found the functions 'sum' and 'subtract' so I'm assuming that it is possible to multiply fields as well with functions...

I also found the function 'scale' but it allows to multiply a field by a scalar only.

Any ideas?

[HPE]

Hi,

- I assume you really meant `function object` by `function`.
- To my knowledge, there is no compiled FO for field multiplications, like there is `add`
- You can use `codedFunctionObject`s.
- You can modify an existing function object for your requirements.

Hope this helps for a start.

[rsaintlege]

Thank you for your answer it helped me a lot. I tried to use codedFunctionObject and I wrote this (I'm not sure about it):

    compute_Power
    {
        functionObjectLibs ( "libutilityFunctionObjects.so" );
        type coded;
        enabled true;
        redirectType Power;
	executeControl  writeTime;
	writeControl    writeTime;
        code
        #{
            const volScalarField& T = mesh().lookupObject<volScalarField>("T");
            const volVectorField& U = mesh().lookupObject<volVectorField>("U");
            volScalarField Power
            (
                IOobject
                (
                    "Power",
                    mesh().time().timeName(),
                    Power.mesh(),
                    IOobject::NO_READ,
                    IOobject::AUTO_WRITE
                ),
                Power= 4185*1000*mag(U)*(T-293)
            );
	    Power.write();
        #};
    }

However, when I run the code, it seems to work fine but when I check the results in paraview, it appears that the field that I call "power" doesn't show up. Did I write something wrong or forget something in the code above?

[rsaintlege]

This is finally the best I managed to do:

compute_Power
	{
	    libs        ("libutilityFunctionObjects.so");

	    type coded;
	    name compute_Power;
	    outputControl  writeTime;

	    codeWrite
	    #{
		    const volScalarField& h = mesh().lookupObject<volScalarField>("h");
		    const volScalarField& h_0 = mesh().lookupObject<volScalarField>("h_0");
		    const volVectorField& U = mesh().lookupObject<volVectorField>("U");
		    volScalarField Power = 1000*mag(U)*(h-h_0);
		    Power.write();
	    #};
	}

That seems to work but the generated files don't have a real name, it's something like "(1000*mag(U)*(h-h_0))". Paraview cannnot therefore open those files. How can I name those files?

[JNSN]

Do you need this field just for post proccesing? Then you could also do this just in paraview.

[rsaintlege]

Originally Posted by JNSN

Do you need this field just for post proccesing? Then you could also do this just in paraview.

Thank you for your answer.

Yes I just need it for post processing actually. How would you do to get this field with paraview? I previously did it by exporting the needed datas as csv files in Excel but I can't do it 100 times

I'm new on OpenFOAM so I don't really know what are the possibilities yet.

I'm still interested about making the previous function object work though. I think it would be more convenient, if you have an idea.

*EDIT*: I finally managed to do it with the 'calculator' option in Paraview

[thomas_mhn]

Originally Posted by rsaintlege

Thank you for your answer.

Yes I just need it for post processing actually. How would you do to get this field with paraview? I previously did it by exporting the needed datas as csv files in Excel but I can't do it 100 times

I'm new on OpenFOAM so I don't really know what are the possibilities yet.

I'm still interested about making the previous function object work though. I think it would be more convenient, if you have an idea.

*EDIT*: I finally managed to do it with the 'calculator' option in Paraview

Hi I'm also interested in making this kind of function.

Did you manage to change the name of the output file? I also have something like ( field1 * field2 ) as the name of the file.

[little_fish]

Hello,



i m still interested to use this postprocess function object, because i don t need to compute this field for each iteration, and i also need the field for calculation (so paraview is not the best option...)



here my coded function :

field_UthetaT
    {
        libs        ("libutilityFunctionObjects.so");

        type coded;
        name UthetaT;
        // executeControl  writeTime;
        // writeControl    writeTime;
        // region         region1;
        // patches         (cyclic1);

        codeWrite
        #{
            const volScalarField& Utheta = mesh().lookupObject<volScalarField>("Utheta.gz");
            const volScalarField& T = mesh().lookupObject<volScalarField>("T.gz");
            volScalarField UthetaT = Utheta * T;
            UthetaT.write();
        #};
        }

and here the error i get when i use the function.

Create time

Create mesh for time = 20000

Using dynamicCode for functionObject field_UthetaT at line 3 in "functions.field_UthetaT"
Time = 20000

Reading fields:

Executing functionObjects


--> FOAM FATAL ERROR: (openfoam-2112 patch=220610)

    failed lookup of Utheta.gz (objectRegistry region0)
    available objects of type volScalarField:
0()


    From const Type& Foam::objectRegistry::lookupObject(const Foam::word&, bool) const [with Type = Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>]
    in file /usr/lib/openfoam/openfoam2112/src/OpenFOAM/lnInclude/objectRegistryTemplates.C at line 463.

FOAM exiting

If someone could help me, it would be very cool.



Regards

[little_fish]

ps: i did not write the function in the control dict because i have already some datas i need to postprocess, so i m using the following command :



postProcess -func field_UthetaT -latestTime



cordially

[cnsidero]

Checkout multiply Class Reference

You can find an example of it in the openfoam.com tutorials:

$FOAM_TUTORIALS/incompressible/pisoFoam/RAS/cavity/system/FOs/FOmultiply