[ashish.svm]

I want to include frictional pressure drop calculation in rhoPisoFoam.

frictional pressure drop can be calculated by 0.5*ff*rho*U^2
where ff is fanning friction factor, which can be calculated by using standard correlations based on reynolds number.

I have modified UEqn as following.

fvVectorMatrix UEqn
(
    fvm::ddt(rho, U)
  + fvm::div(phi, U)
);

if(includeFrictionFactor)
{
  UEqn += fvm::div(phiFrictional, U);
}

phiFrictional is calculated as following.

phiFrictional =
    fvc::interpolate(fanningFrictionFactor)*
    (
        (fvc::interpolate(U) & mesh.Sf())
      + fvc::ddtPhiCorr(rUA, rho, U, phiFrictional)
    );

and added the following term in pressure equation

if(includeFrictionFactor)
    {
      pEqn += fvc::div(phiFrictional);
    }

but using this formulation leads to incorrect results.

Please let me know what is wrong in my formulation.

Or

Please suggest me some other solver which has the pressure drop calculation using fanning friction factor.

[Santiago]

Originally Posted by ashish.svm

I want to include frictional pressure drop calculation in rhoPisoFoam.

frictional pressure drop can be calculated by 0.5*ff*rho*U^2
where ff is fanning friction factor, which can be calculated by using standard correlations based on reynolds number.

I have modified UEqn as following.

fvVectorMatrix UEqn
(
    fvm::ddt(rho, U)
  + fvm::div(phi, U)
);

if(includeFrictionFactor)
{
  UEqn += fvm::div(phiFrictional, U);
}

phiFrictional is calculated as following.

phiFrictional =
    fvc::interpolate(fanningFrictionFactor)*
    (
        (fvc::interpolate(U) & mesh.Sf())
      + fvc::ddtPhiCorr(rUA, rho, U, phiFrictional)
    );

and added the following term in pressure equation

if(includeFrictionFactor)
    {
      pEqn += fvc::div(phiFrictional);
    }

but using this formulation leads to incorrect results.

Please let me know what is wrong in my formulation.

Or

Please suggest me some other solver which has the pressure drop calculation using fanning friction factor.

Its wrong, the velocity in your formula is the square of the magnitude of the velocity, or inner product. Besides, why you take the divergence of such force?

[ashish.svm]

Originally Posted by Santiago

Its wrong, the velocity in your formula is the square of the magnitude of the velocity, or inner product. Besides, why you take the divergence of such force?

Sorry, my understanding of the friction factor was wrong. I found out that this force is applied to the wetted surface.

But I am unable to understand how to modify the momentum equation to include fanning friction term.