To calculate NORMAl=normal*are, that is always outward, on can call following function from kratos/utilities
BodyNormalCalculationUtils().CalculateBodyNormals((fluid_model_part).Elements(),domain_size)
