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