[palarcon]

Hello everybody

I'm working on topology optimization with OpenFOAM and I'm using and adapting the solver adjointShapeOptimizationFoam.
In order to do that, I need to compute and print the objective function used on that solver (which is not the same used on the paper which is based on).
The real objective function used in this solver is the total pressure drop, which mathematically is



I created a .H file to compute that pressure drop, but I'm not really sure that the implementation is correct, given the fact that the pressure is multiplied by mesh.Sf (I always gets confused with the normalized pressure used in simpleFOAM, which the solver is based on). Could somebody please tell me if I'm right?

scalar jDissPower (0);

for (iLoop=0; iLoop<nObjPatch; iLoop++)
{
	if (objFunction==1) {
		jDissPower =
		jDissPower - sum(
		(p.boundaryField()[objPatchList[iLoop]]*mesh.Sf().boundaryField()[objPatchList[iLoop]]
		+ 0.5*magSqr(U.boundaryField()[objPatchList[iLoop]]))
		);
	}
}

where nObjPatch is the number of patchs (usually 2, inlet and outlet) and the sum function to compute the integral over the patchs.

Thank you all in advance