[NewFoamer]

Hello,

My current solver is based on the incompressible pimpleFoam, in which I've added a solver for solute concentration C:

solve
    (
    fvm::ddt(C)
    + fvm::div(phi, C)
    - fvm::laplacian(Diff, C)
    );

Now, I would like to have the fluid density depend on the solute concentration according to:

rho = base_value + rho_factor * C

Not having much experience with OpenFOAM, I'm a bit in the dark about how to go about implementing this. I can't claim to fully understand the compressible solvers I've looked at yet, but I'm not sure I have to go there either? Already having the C field in place, it shouldn't be neccesary to create a rho field as well? given the simple relation that is.

Currently, the momentum equation is solved as follows:

tmp<fvVectorMatrix> UEqn
(
    fvm::ddt(U)
  + fvm::div(phi, U)
  + turbulence->divDevReff(U)
);

Basically, the question is, how do I go about implementing this? any hints or help would be appreciated!