Fig. 1. Frozen water cylinder 2D cross-section

Fig. 2. Distribution of PH electric field under y-polarized x-axis incident plane wave illumination

Fig. 3. When$\mathrm{o}\mathrm{p}\mathrm{e}\mathrm{n}\mathrm{i}\mathrm{n}\mathrm{g}\mathrm{a}\mathrm{s}\mathrm{y}\mathrm{m}\mathrm{m}\mathrm{e}\mathrm{t}\mathrm{r}\mathrm{y}b=\mathrm{1.56},\mathrm{o}\mathrm{p}\mathrm{e}\mathrm{n}\mathrm{i}\mathrm{n}\mathrm{g}\mathrm{s}\mathrm{y}\mathrm{m}\mathrm{m}\mathrm{e}\mathrm{t}\mathrm{r}\mathrm{y}b=a$，the relationship between the characteristic parameters of PH and M coordinates

Fig. 4. Distribution of PH electric field under y-polarized x-axis incident plane wave illumination

Fig. 5. The relationship between the characteristic parameters of PH and the depth of the depression in seven cases，$b=a=\mathrm{0.26},b=a=\mathrm{0.52},b=a=\mathrm{0.78},b=a=\mathrm{1.04},b=a=\mathrm{1.30},b=a=\mathrm{1.56},b=a=\mathrm{1.82}$

Fig. 6. Adjust coefficient A、 B、 and C to control the structure of frozen water cylinder

Fig. 7. Distribution of electric field for different coefficients A，B，C and D