[Franko]

Hi!
Could anybody tell me, is there any difference between using these two solver's codes? Can I use the first code instead of the second code and vice versa? (Personally I think they have the same result):

So, the first code:

    surfaceScalarField phi ("phi", A);

    surfaceScalarField phiflux = fvc::flux(phi, B, "div(phi,B)");

     fvScalarMatrix sEqn
        (
              fvm::ddt(s) + fvc::div(phiflux)
        );    
     sEqn.solve();

and the second code:

     surfaceScalarField phi ("phi", A);

     fvScalarMatrix sEqn
      (
        fvm::ddt(s) + fvc::div(phi, B, "div(phi,B)")
      );
     sEqn.solve();

[sharonyue]

The reason to use fvc::flux() MIGHT be to use flux limiter scheme. such as vanleer~

http://www.cfd-online.com/Forums/ope...interfoam.html