[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