[mostafa kareem]

Hello,
I was solving emprical model(suzen model) for DBD plasma flow control. the target is to calculate the body force contours which equal to

Body force = charge density*ElectricFied &
ElectricFied = grad(electric potential).

I created the solver for charge density and electric potential. the contours for both of them match well with validation result, However On refining the mesh the intgrated body force remain the same, but the mximum value keeps increasing as a result of singularity at the electrode corner. I think the problem with E = grad(electric potential) : gradient operator at sharp coner , I kept changing gradient schemes (cell based and node based ), but the problem ramains unsolvable, so any suggestion on how to solve the problem