Please run the attached code. The density plot has some artifacts on it (The light blue parts around the circle) Is there a way to resolve this issue. The density plot is showing velocity of a liquid passing through the channel. Regards.
[CapitalOmega] = ImplicitRegion[ 0 <= x <= 2.0 && 0 <= y <= 0.48208 && ! (x >= 1.0 && y <= 0.19) && ! (x >= 1.0 && y >= 0.31) && ! ((x - 0.76)^2 + (y - 0.24104)^2 <= 0.238125^2), {x, y}]; RegionPlot[[CapitalOmega], AspectRatio -> Automatic, ImageSize -> 300] stokesFlowOperator = {Div[({{-1, 0}, {0, -1}}.Grad[ u[x, y], {x, y}]), {x, y}] + \!(*SuperscriptBox[(w), TagBox[ RowBox[{"(", RowBox[{"1", ",", "0"}], ")"}], Derivative], MultilineFunction->None])[x, y], Div[({{-1, 0}, {0, -1}}.Grad[v[x, y], {x, y}]), {x, y}] + \!(*SuperscriptBox[(w), TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None])[x, y],
\!(*SuperscriptBox[(v), TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None])[x, y] + \!(*SuperscriptBox[(u), TagBox[ RowBox[{"(", RowBox[{"1", ",", "0"}], ")"}], Derivative], MultilineFunction->None])[x, y]}; Subscript[[CapitalGamma], D] = {DirichletCondition[{u[x, y] == -1, v[x, y] == 0}, x == 0.], DirichletCondition[{u[x, y] == 1., v[x, y] == 0}, 0 < x < 2.0], DirichletCondition[w[x, y] == 0.1, x == 2.0]}; {xVel, yVel, pressure} = NDSolveValue[{stokesFlowOperator == {0, 0, 0}, Subscript[[CapitalGamma], D]}, {u, v, w}, {x, y} [Element] [CapitalOmega], Method -> {"FiniteElement", "InterpolationOrder" -> {u -> 2, v -> 2, w -> 1}}]; rmf = RegionMember[[CapitalOmega]]; Show[BoundaryDiscretizeRegion[[CapitalOmega]], DensityPlot[{xVel[x, y], yVel[x, y]}, {x, 0, 2.0}, {y, 0, 0.48208}, RegionFunction -> Function[{x, y}, rmf[{x, y}]], AspectRatio -> Automatic, PlotLegends -> Automatic, ColorFunction -> ColorData["BlueGreenYellow"], PlotPoints -> 20], ImageSize -> 1000]
Attachments:
