5

We can use FEM as it is to solve this problem in versions 12-12.1.1 Needs["NDSolve`FEM`"]; k = 2; Pr0 = 1; Ra = 10^4; R = Ra*Pr0; a = 0; d = 1/2; b = 1/3; Da = 1; Ha = 1; Q = 1; reg = Rectangle[{0., 0.}, {1., 1.}]; {UX, VY, PK, TK} = NDSolveValue[ {{Inactive[Div][{{-\[Mu], 0}, {0, -\[Mu]}} . Inactive[Grad][u[x, y], {x, y}], {x, y}] + ...


4

One way to deal with this would be to change the blotc function to generate assymmetric blobs in the first place. blotc2[smoothness_ : 20, points_Integer : 10, symm_ : 0.5] := With[{fun = Exp[-smoothness #.#] &, pts = Transpose[{RandomReal[1 - symm, points], RandomReal[symm, points]}]}, With[{fc = Compile[{xl, yl}, Total[fun[# - {xl, yl}] &...


4

VoronoiMesh[X, MeshCellStyle -> {{1, "Interior"} -> {Thick, Red}, {1, "Frontier"} -> {Thick, Red}}] Alternatively, VoronoiMesh[X, MeshCellStyle -> {{1, "Boundary"} -> Opacity[0], {1, All} -> Directive[Thick, Red]}] same picture


3

Following the discussion in the comment section with @TumbiSapichu, I've found a possible solution to this problem. As mentioned, instead of translating the seeds, we could simply add more points, enough so that, upon drawing a rectangle centred in this new mesh, you simply pick the first n cells which seeds intersect the rectangle, with increasing size, ...


2

ListLinePlot[data, PlotRange -> {{0, 0.5}, {0, 1}}, PlotStyle -> Arrowheads[ConstantArray[.02, 10]], Frame -> True, ImageSize -> Large, PlotRangeClipping -> False, Epilog -> {Black, PointSize[Large], Point[e0], Point[e1], Point[e2], Point[e3]}] /. Line -> Arrow


1

There are two perspectives in the question a) PlotRange -> All b) Filling -> {1 -> {2}} to a) It is necessary to put in both intervals explicitly into the options of PlotRange: PlotRange -> {{-1.*10^-15, 1.*10^-15}, {-1.0, 2.1}} will present the desired result. to b) This appears to me as taken from the example for Filling Fill between curves 1 ...


1

Instead of Labeled, use Text Clear["Global`*"]; eqns = {(1 - x)/2, -2 - 3 x}; pt = {x, eqns[[1]]} /. Solve[Equal @@ eqns, x, Reals][[1]] (* {-1, 1} *) plot = Plot[eqns, {x, -3, 1}, PlotRange -> {-1, 3}, PlotStyle -> Thick, PlotLabels -> {y == (1 - x)/2, y == -2 - 3 x}, AxesLabel -> {x, y}, ImageSize -> {Automatic, 200}];...


Only top voted, non community-wiki answers of a minimum length are eligible