[itsme_kit]

Hi

Any one has template in writing these equations in C++

Thanks in advance

[chriss85]

Just take a look at a solver that solves a problem similar to your task.

[Annier]

Hi,
I am replying an answer to this question:
How does the C++ library of OpenFOAM represent Navier-Stokes equation? (I request you to edit your question if it is somewhat like this.)

The OpenFOAM source library describes the Navier-Stokes equation in its applications/solvers.
For examples,
1.In OpenFOAM/OpenFOAM-2.3.x/applications/solvers/incompressible/icoFoam/icoFoam.C, the Navier-Stokes equation is represented as:

fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu, U)
        );

        solve(UEqn == -fvc::grad(p));

//--- PISO LOOP
...

2. Similarly, in OpenFOAM/OpenFOAM-2.3.x/applications/solvers/incompressible/pisoFoam/pisoFoam.C, the Navier-Stokes equation is represented as:

 // Momentum predictor

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

            UEqn.relax();

            if (momentumPredictor)
            {
                solve(UEqn == -fvc::grad(p));
            }
//--- PISO LOOP
...

3. ...
4. ...

Yours
Anil Kunwar