{"id":90,"date":"2022-09-07T09:34:20","date_gmt":"2022-09-07T00:34:20","guid":{"rendered":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en2\/?post_type=case&p=90"},"modified":"2022-09-07T11:02:39","modified_gmt":"2022-09-07T02:02:39","slug":"periodic_boundary_surface","status":"publish","type":"case","link":"https:\/\/www.ssil.co.jp\/product\/EMSolution\/en\/case\/periodic_boundary_surface\/","title":{"rendered":"Current Field Analysis on Periodic Boundaries"},"content":{"rendered":"

Summary<\/h3>\n

When defining "PHICOIL"<\/font><\/a>, "SUFCUR"<\/font><\/a> and "DCCURR"<\/font><\/a> as coils and bulk conductors across a periodic boundary surface, it was previously necessary to define gap elements on the coil or conductor cross section inside the mesh. This was no longer possible when defining a current inflow surface at the periodic boundary with a gap element. We are pleased to report that it is now possible to define the current inflow surface on the periodic boundary without gap elements. The analysis can also be performed for steady-state current field analysis. <\/p>\n

Explanation<\/h3>\n

This function is a type of periodic boundary condition, also called an equidifferential periodic boundary condition (Ref. [1]). To provide a potential difference between periodic boundaries, a gap element was used to provide a potential difference as in the 360-degree model, but now it can be solved by providing a potential difference between the current inflow and outflow surfaces. This allows the "PHICOIL at periodic boundaries"<\/font><\/a> and "PHICOIL in rotational z-anti-mirroring symmetry"<\/font><\/a> models to be analyzed without gap elements at the periodic boundaries. <\/p>\n

\u25a0 PHICOIL at periodic boundaries (without gap element)<\/h4>\n
\n

As a validation calculation, we use the model of "PHICOIL at the periodic boundaries"<\/font><\/a> to define a PHICOIL current input surface at the periodic boundaries.
\nFig. 1 shows the calculation results together with the results defined using gap elements inside the mesh. It can be seen that the results are in agreement, as expected. Although the gap element must be defined one layer beyond the coil cross section, it is now possible to define the current input surface only on the coil cross section as well as the mirror symmetry surface on the periodic boundary plane, which simplifies the setup.<\/p>\n

\n
\n
\n \"\"<\/a>
\n<\/p>\n

(a) Gap Element Definition<\/p>\n<\/p><\/div>\n

\n
\n \"\"<\/a>
\n<\/p>\n

(b) Periodic boundary surface definition<\/p>\n<\/p><\/div>\n

Fig.1 Current density vector distribution $[A\/m^2]$<\/p>\n<\/div>\n

\u25a0 Periodic boundary conditions in steady-state current field analysis<\/h4>\n
\n

Consider the model shown in Fig. 2, in which the conductor with holes is modeled with periodic boundaries, and a steady-state current field analysis[1]<\/a><\/sup> is performed. Assume that the periodic boundary condition is satisfied in the x-direction and that current continuity is satisfied. In this case, the use of periodic boundary conditions means that a part of the analytical model is taken out, so the voltage applied to both ends of the model is input. If the model is a 1\/10 model, the analysis can also be done by inputting the voltage across the entire model and setting REGION_FACTOR=10. <\/p>\n

Fig. 3 shows the current distribution obtained for the periodic boundary condition and the mirror symmetry condition. In the case of the mirror symmetry condition, the current flows in and out perpendicularly at the boundary, whereas in the periodic boundary condition, the current density has oblique components on both sides. In this model, there is only a slight difference, but depending on the problem, the difference may not be negligible, so we hope you will use this function. <\/p>\n

\n
\n
\n \"\"<\/a>\n <\/div>\n

Fig.2 Model of conductor with holes<\/p>\n<\/div>\n

\n
\n
\n \"\"<\/a>
\n<\/p>\n

(a) Mirror symmetry condition<\/p>\n<\/p><\/div>\n

\n
\n \"\"<\/a>
\n<\/p>\n

(b) Isoperiodic boundary condition<\/p>\n<\/p><\/div>\n

Fig.3 Current density vector distribution $[A\/m^2]$<\/p>\n<\/div>\n

Briefly, we have introduced a function that allows the analysis of periodic boundaries using PHICOIL and steady-state current field analysis without gap elements. As mentioned earlier, this function can be used with "PHICOIL"<\/font><\/a>, "SUFCUR"<\/font><\/a>, "DCCURR"<\/font><\/a> and "Steady Current Field Analysis Function"<\/font><\/a>. As with the mirror symmetry surface (Bn=0 surface), a current input surface can be defined as the periodic boundary surface.<\/p>\n

References<\/h3>\n

[1]: Tokumasu, Fujita: "Remaining issues in two-dimensional analysis," Joint Research Institute of Electrostatic Stator and Rotating Machinery, SA-05-73, RM-05-80 (2005)<\/small><\/p>\n

<\/p>\n

How to use<\/h3>\n

\u25a0 PHICOIL at periodic boundaries (without gap element)<\/h4>\n

Define PHICOIL (Handbook "17.4 Potential Current Sources"). Define a periodic boundary current inflow surface for SMAT_ID_IN.<\/p>\n

* PHICOIL * SERIES_NO * NO_PARTS *
\n PHICOIL 1 1
\n* MAT_ID * S_MAT_ID_IN * CURRENT * SIGMA *
\n 1 102 1.0 5.0e7\n<\/p>\n

\u25a0 Periodic boundary condition in steady-state current field analysis (isoperiodic boundary condition)<\/h4>\n

Define periodic boundary conditions (Handbook "14. Periodic Boundary Conditions, Sliding Surfaces"). <\/p>\n

* CYCLIC *<\/font> SYMMETRY * SLIDE_NDIV * SLIDE_MAT * FITNESS * LINE_INPUT * N_CYCLIC *
\n 1<\/font> 0 0 0 0.00000 0 1
\n* DX * DY * DZ *<\/font>
\n 2.5000000e+000 0.0000000e+000 0.0000000e+000<\/font>\n<\/p>\n

Download<\/h3>\n

\u25a0 PHICOIL at periodic boundaries\u3000\u3000