[blackk100]

I'm setting up a custom modular solver to simulate electrochemical cells in OpenFOAM 11. I've created a derived class from multicomponentFluid that compiles and runs.

I'm currently attempting to add in another scalar field variable for the electric potential, but I cannot figure out how to do that exactly (OF10 & OF2312 are straightforward with this).

This is what I attempted (does not compile):
customFluid.C:

Foam::solvers::customFluid::customFluid(fvMesh& mesh)
:
multicomponentFluid(mesh),
E_(
    IOobject("E", runTime.name(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE),
    mesh
),
E(E_)
{}

customFluid.H:

...
class customFluid: public multicomponentFluid
{
protected:
    volScalarField &E_;
private:
public:
    const volScalarField &E;

   TypeName("customFluid");
...
}

I did notice the scalarTransport function, but it does not cover all of my needs (I need to add terms dependent on E to the U and Y equations).

Would appreciate any help on this!