Geometrical Platform of Big Database Computing for Modeling of Complex Physical Phenomena in Electric Current Treatment of Liquid Metals

: According to the principles of multiphysical, multiscale simulations of phenomena and processes which take place during the electric current treatment of liquid metals, the need to create an adjustable and concise geometrical platform for the big database computing of mathematical models and simulation is justiﬁed. In this article, a geometrical platform was developed based on approximation of boundary contours using arcs for application of the integral equations method and matrix transformations. This method achieves regular procedures using multidimensional scale matrices for big data transfer and computing. The e ﬃ ciency of this method was veriﬁed by computer simulation and used for di ﬀ erent model contours, which are parts of real contours. The obtained results showed that the numerical algorithm was highly accurate based on the presented geometrical platform of big database computing and that it possesses a potential ability for use in the organization of computational processes regarding the modeling and simulation of electromagnetic, thermal, hydrodynamic, wave, and mechanical ﬁelds (as a practical case in metal melts treated by electric current). The e ﬃ ciency of this developed approach for big data matrices computing and equation system formation was displayed, as the number of numerical procedures, as well as the time taken to perform them, were much smaller when compared to the ﬁnite element method used for the same model contours.


Introduction
The study of complex processes and phenomena by mathematical modeling is a natural way to improve and enhance the efficiency of technical methods and technologies.Today, efficiency upgrades in every industry and scientific field involve data processing and mathematical modeling.In some areas of human activity, information technology for data processing is an integral part of production process management.Mathematical models are also often used alongside data obtained from nature experiments, and physical model construction is firmly established in the practice of engineering design.
There are also technical spheres where the experimental study of the internal content of current processes and real measurements of local physical parameters of reagents or interacting media are not currently possible.In these cases, the only way to identify the influence of setting exposures on controlled parameters, intermediate data, and final results is by using a computational method that utilizes generalized data of empirical dependencies or numerical simulations based on existing theories of the corresponding phenomena.In this article, special attention is paid to the geometric platform of big database computing to model complex physical phenomena in the electric current treatment of liquid metals.

Problem Statement
In light of the noted specificity of the problems of numerical simulation of CECT processes, a promising approach to solving a set of mathematical physics equations seems to be the application of a more concise and compact method instead of (or together with) FEM, based on an integral representation of the characteristics of physical fields in the area where they interact with the medium.This especially concerns the situation when the medium, besides all, changes its properties.
In general, this approach is called "the method of integral equations (IEs)".Its advantages become apparent in situations where the spatial interaction region (i.e., the set where a process or phenomenon is localized) can be represented equivalently on a set of smaller dimensionality: from the volume of the body (space) to its surface area; and in two-dimensional problems, from the figure area to the boundary contour [24].
In fact, such a transition is completely analogous to the representation of solutions of the Dirichlet and Neumann problems for the differential Laplace and Poisson equations by Green's integral formulas, in which the main properties of physical fields are displayed by simple and double-layer potentials [25].These potentials are expressed by integral operators (IOs) with the kernels K 1 (P,M) and K 2 (P,M) of the form: where r PM = (x P − x M ) 2 + (y P − y M ) 2 + (z P − z M ) 2 is the radius vector drawn from the integration point M on the surface (or contour) to some observation point P.That is, kernels are functions of two points' coordinates; the kernel K 2 is actually a partial derivative of K 1 in the direction normal to the surface of integration.The derivative above is calculated by the following formulae.
For the 3D model: and for the 2D model: where cos(r PM , n P ) = cos(r PM , → e x ) • cos(n P , → e x ) + sin(r PM , → e y ) • sin(n P , → e y ) is the cosine of the angle between radius-vector r PM and normal n P at point P ( → e x is the unit vector of axis x, and → e y is the unit vector of axis y).
In most cases, the numerical solution of an IE with kernels of the specified type is carried out by replacing the integrals with finite sums having a sufficiently high degree of discretization of the integration surface (contour) and reducing them to systems of linear algebraic equations (SLAEs).However, it should be noted that, depending on the boundary conditions and parameters connected with the characteristics of the media, the SLAEs may be ill-conditioned-the task becomes incorrect [26], and its solution requires the application of special methods such as singular value decomposition (SVD).The effectiveness of this procedure has been repeatedly noted [27,28] (particularly in [29]), but its implementation is connected with multiple matrix transformations, the computation of which requires the creation of a sufficiently large database related to the configuration of the integration surfaces.
From the above, it is clear that in the transition from IE to SLAE and the success of their solution, the primary point of methodological importance is the procedure of the rational presentation of the geometrical parameters of the corresponding sets, accounting for the need to use them explicitly in the algorithmic structure of mathematical operations of differentiation and integration on curvilinear surfaces, as well as calculating two-point operators in integral equations.Thus, in order to ensure the computational process, first of all, for the formation of matrices of coefficients of the SLAE, which approximates the set of IEs of the complex of physical processes of CECT, it is necessary to create a ramified scheme for calculating the geometric parameters of the elements of curvilinear surfaces and contour lines, their curvatures, tangents and normal vectors to them, and so on-that is, a kind of geometrical platform for the prompt update of large databases of numerical models.
The basis for the implementation of the above approach should be a rational method of setting the configuration of computational areas (bodies) and the optimal procedures for simulating the transduction of their influence into local zones of field.Medium interactions should serve as the basis for the implementation of the stated approach.
Despite the existence of a number of rather good algorithmic procedures [30], the composition of adequate geometric models of complex technical objects for the numerical analysis of various physical processes is still an urgent problem, since the well-known and accessible software from this arsenal is not always suitable due to its profile specialization for some classes of tasks [31].
Particularly, traditional combinational or generational approaches impose definite limitations: firstly, they practically do not allow any variations; secondly, being a locally deterministic procedure, they do not guarantee the convergence of the grid to a given area on a complicated surface [32].
It is precisely this situation that we encounter when considering a complex of interconnected physical phenomena with a high degree of uncertainty of the external and internal boundaries of their localization areas (patches, bubbles) in the CECT process.
The purpose of this work was to create a geometrical platform for the operative formation and processing of large databases of numerical models of physical processes and phenomena in the technological operations of a foundry, which are inaccessible to direct observation.

Method for Specifying Geometric Objects of Complicated Shape
The main idea of the procedure for forming a geometrical platform which serves as the basis for the computational operations of the outlined problems is similar to that presented in [2].In this paper, this concept has undergone further development and acquired a complete form.Unlike the three conceptually different approaches specified in [30], it is synthetic, uniting the elements of object modeling according to drawings of their projections, the parametric description of surfaces, and the composition of basic topological primitives.Its essence is as follows.
Virtually any complicated detail or construction can be described (depicted) by a set of its projections and sections by parallel (or non-parallel) planes, so that the borders of the computational areas will represent closed contours formed by curves or pieces of straight lines.Examples include wing profiles, ship hull frames, magnetohydrodynamic (MHD) flow channels, melt ladle profiles, etc.A set of such two-dimensional (flat) contours allows a 3D-model of objects to be obtained from them rather simply from slices.In [32], this technique was called the "method of frame sections".A similar idea was implemented with the help of topological primitives and R-functions in [33], but this approach seems rather cumbersome for tasks with multiple runs of varying modes of CECT processes.
So, the problem comes down to the task of assigning a set of closed curvilinear contours with different curvature segments in the simplest way.It is almost always possible to approximate such a contour with sufficient accuracy by a sequence of joined arcs of circles of suitable radius.The agnetoarsenal of approximation methods is unlimited: from the visual splitting of the contour into segments and the assignment of coordinates of the docking points, to the fully automated conversion of scanned drawings into digital arrays.
Thus, in our algorithm an arc of a circle is taken as a basic primitive; its limiting forms are the corner point with a zero (conditionally) radius, and a straight line segment with an infinitely large radius.
For a representation of the preparatory stage of the procedure for forming the geometrical platform, let us consider the configuration of the ladle section with the molten metal, into which a pulsing current is fed through the electrodes (Figure 1).For simplicity, we consider that the ladle is extended along the Z axis (beyond the drawing plane), and by following literally all the bends of the border surfaces, we obtain a shape whose contour needs to be approximated by the mentioned primitives.According to Figure 1, we visually select 11 segments of straight lines and 2 arcs of the contour on its curvilinear part.Thus, we define the total number of approximating segments (arcs): M = 13.
counterclockwise according to the standard positive direction at planimetric integration.The initial (IP) and the final (EP) points of each arc are designated by the corresponding indices, so that the end points of the previous arcs coincide with the initial points of the subsequent arcs: EPk ≡ IPk + 1 (Figure 2).
In the next step we specify the coordinates (Xip,j, Yip,j) of the first and every subsequent initial points of arcs, which coincide with the coordinates of the end points of previous arcs (Xep,j−1,Yep,j−1), where j is the arc number from the common quantity of arcs M.
After that we define the curvature radii of the arcs Rc,j, but as for their value there exists the restriction that Rc,j must be not less than half of the distance la,j from the beginning to the end of the arc (i.e., of the chord).If Rc,j = 0.5 la,j, it means the arc is a semicircle; if Rc,j > 100 la,j, then the arc becomes a piece of the straight line.Then, we perform the operation of specifying the spatial position and key geometric parameters of the arcs.
Firstly, all arcs and contour segments, starting at any point, are numbered one by one counterclockwise according to the standard positive direction at planimetric integration.The initial (IP) and the final (EP) points of each arc are designated by the corresponding indices, so that the end points of the previous arcs coincide with the initial points of the subsequent arcs: EP k ≡ IP k + 1 (Figure 2).
primitives.According to Figure 1, we visually select 11 segments of straight lines and 2 arcs of the contour on its curvilinear part.Thus, we define the total number of approximating segments (arcs): Then, we perform the operation of specifying the spatial position and key geometric parameters of the arcs.
Firstly, all arcs and contour segments, starting at any point, are numbered one by one counterclockwise according to the standard positive direction at planimetric integration.The initial (IP) and the final (EP) points of each arc are designated by the corresponding indices, so that the end points of the previous arcs coincide with the initial points of the subsequent arcs: EPk ≡ IPk + 1 (Figure 2).
In the next step we specify the coordinates (Xip,j, Yip,j) of the first and every subsequent initial points of arcs, which coincide with the coordinates of the end points of previous arcs (Xep,j−1,Yep,j−1), where j is the arc number from the common quantity of arcs M.
After that we define the curvature radii of the arcs Rc,j, but as for their value there exists the restriction that Rc,j must be not less than half of the distance la,j from the beginning to the end of the arc (i.e., of the chord).If Rc,j = 0.5 la,j, it means the arc is a semicircle; if Rc,j > 100 la,j, then the arc becomes a piece of the straight line.In the next step we specify the coordinates (X ip,j , Y ip,j ) of the first and every subsequent initial points of arcs, which coincide with the coordinates of the end points of previous arcs (X ep,j−1 ,Y ep,j−1 ), where j is the arc number from the common quantity of arcs M.
After that we define the curvature radii of the arcs R c,j , but as for their value there exists the restriction that R c,j must be not less than half of the distance l a,j from the beginning to the end of the arc (i.e., of the chord).If R c,j = 0.5 l a,j , it means the arc is a semicircle; if R c,j > 100 l a,j , then the arc becomes a piece of the straight line.
Moreover, we assign to the radius the sign: "+" or "−".A positive value of the R c,j indicates that the center of arc curvature is situated on the right side of chord direction, passing the arc from the initial point to the end point.Accordingly, a negative value of Rc,j indicates that the center of the arc curvature is situated on the left side of the chord direction.
The position of the center of arc curvature relative to the chord, which is oriented from the beginning to the end along the contour, characterizes the bending of the arc-that is, its concavity or convexity.If the radius is positive, the arc will be concave, will pass to the left of the chord, penetrating into the area covered by the contour, and reducing its area; with a negative radius, the arc will be convex, located to the right of the chord, respectively, expanding the area.
Data 2019, 4, 136 6 of 18 Therefore, we have a set of parameters for each approximating arc: X, Y coordinates for initial and end points, as well as the radius of curvature with the sign ±R.For convenience of use and the further computer processing of these data by the software of widespread versions of MATLAB or other languages, it is expedient to present them as a compact file in the form of a text file.Into this file, we enter in the usual way a two-dimensional matrix with the size of 5 × M, where M is the total number of approximating arcs.The structure of the matrix is as follows: 1st line: Xip,j-coordinates of the initial points of the arcs; 2nd line: Yip,j-coordinates of the initial points of the arcs; An example listing (listing 1) of such a file is given below (see Table 1); the listing data represent a contour in the form of a circle of unit radius R c = −1 with the center at the beginning of the coordinates (0,0, 0), which is composed of six arcs of 60 degrees each with the countdown from the abscissa axis.

Calculation of the Geometrical Parameters of Approximating Arcs
To perform the numerical integration of a vector function along a curvilinear contour, it is necessary, to divide it into elementary sections (i.e., arcs), which involves calculating the coordinates of the nodal points (i.e., the points of the conjugation of elementary arcs and their lengths, as well as the direction cosines of the tangents to them).In addition, IO kernels of the form (1) require the operation of the form (2), (3) by means of which from the scalar product of vectors r PM and n P the cosine of the angle between them is found.
The accuracy of integration is ensured by setting the smallest possible size of elementary arcs; however, their total number is limited by the capabilities of computing means.Now, on the basis of item 3 data, it is possible to calculate every geometrical parameter necessary for the computation of the kernels of integral operators, which is executed in a few stages.
First of all, we need to calculate the coordinates of arcs' curvature centers, sector angles, directions of chords, and normals to the corresponding arcs.
To get the formulae for calculating the coordinates (X cc,j ,Y cc,j ) of the arcs' curvature centers and the sector angles γ j that correspond to each arc, let us consider Figure 3, which shows the coordinate plane with one arc as a separate part of the entire contour.The main designations in Figure 3 are IP-initial point, EP-end point, CC-center of the circle, MP-middle point of the arc's chord, α-angle of the arc chord inclination to the vertical, β-angle of the arc chord inclination to the horizontal.
First, we perform some elementary auxiliary counting: -The coordinates of the CP point are calculated as the half-sum of corresponding coordinates of the initial and end points of a given arc half of the length of IP-EP arc chord denoted by l c : -The height h c of an isosceles triangle, formed by the chord and two radii R c,l and R c,r (left and right), drawn from the common curvature center of the arc to its ends IP and EP, which is the bisectrix of the sector angle γ as well: -The length of each of the arcs and the total length of the contour: For the convergence of subsequent calculations, it is very important to identify correctly the position of the chord and the direction of the normal with respect to the axes of coordinates by means of the angles α and β (Figure 3).In [2], this situation was resolved by the combinatorial method-by searching all possible variants of the arc, chord, and normal in the quadrants of the coordinate system, taking into account the convexity or concavity of the arc.Although the result thus obtained is reliable, the logic of the computational algorithm was rather cumbersome.
Note that it is possible to simplify the procedure for determining these angles: the normal to the arc is the prolongation of the radius drawn from the center to the midpoint of the arc, and its direction coincides with the direction of this radius.In this case, as previously indicated, the chord perpendicular to this radius has an orientation from the beginning of the arc to its end.That is, it can be Data 2019, 4, 136 8 of 18 considered as a vector which is rotated with respect to the normal at an angle of 90 • counterclockwise.If such a vector is presented in complex form by the number then its module will give the length of the chord, where lc is a half-chord according to (4), and the argument will determine the required angle: Operations ( 8)-( 10) are directly performed by MATLAB functions, and the result is transmitted without any additional manipulations to the corresponding executable block of the calculation program.We give a concise listing of these calculations in Table 2, where the number z ch is given the name "comHord".The next parameter related to the angle β is the angle α, which is the angle of inclination of the normal to the arc at its midpoint to the horizontal, coinciding with the angle of inclination of the chord to the vertical.Thus, the normal in the complex plane is 90 • behind the chord, and therefore the angle α is determined quite simply: Finally, we determine the coordinates of the point CC (i.e., the center of arc curvature): The last item of the set of approximating arcs parameters that will be required for further calculations are the angles α IP and α EP , which determine the slope of the normals at the initial and end points of a given arc, which serve as reference points:

Procedure of Approximating Arcs Partition into Elementary Sections and Determination of Their Geometric Parameters
To obtain a numerical solution of integral equations (IEs) which describes the physical phenomena in the CECT process, it is necessary to reduce from their continuous (integral) formulation to a discrete one as a sum of values of operators like (1) on elementary parts of the integration contour breakdown.At the same time, the number of such parts, that is, arc finite elements (AFEs), should be large enough to ensure the necessary accuracy of the finite element approximation of contour integrals.
Experience with similar calculations has shown that acceptable (engineering) accuracy can be achieved by setting the size (length) of the SLAE to at least less than 1% of the length of the integration contour or, in other words, more than 100 AFEs per contour.This means that the square matrix of the Data 2019, 4, 136 9 of 18 SLAE coefficients, which corresponds to the IEs, will contain more than 10 4 elements.If we take for example the contour in the form of a full circle of 360 • , represented by Table 1, and divide it into AFEs by 1 • , then we obtain a matrix of 130,000 elements.
Therefore, the first thing required to define the set of AFEs is to appoint the rational number of partitions of each arc into elementary sections and total the number of them.
The number of elementary sections m j into which every arc is parted is preset in such a way that the ratio ∆l j = l j /m j ≈ Σl j /Σ m j = L/M is approximately maintained.Here j is the number of the approximating arc along the contour.
The values of m j may be set initially together with the array entArray when entering the input data of arcs approximating the contour, through the program interface (if provided), or manually using the command input ('arcNDel ='), where arcNDel is a vector-string of input values.
Then, we calculate: -The sector angle of each elementary arc ∆γ j = γ j /m j ; - The corresponding lengths of each AFE as ∆l j = l j /m j .
On the basis of these data, using the values of angles ( 13), we calculate the coordinates of the nodal points, that is, points of AFEs linking: where the index i denotes the number of a nodal point on the j-th approximating arc.Now we can find the most important parameters of each AFE-the angle β j,i of inclination of its chord to the horizontal, and then the angle α j,i of inclination of the normal n j,i at the midpoint of the i-th elementary arc to the x-axis.To do this, using the coordinates of the nodal points (x j,i , y j,i ), we again turn to the complex form of representing the vectors of chords and normals and apply the same procedure to each elementary arc, which is described by Formulas (8), (10), and (11): Finally, we form a system of collocation points (ξ j,i , η j,i ) at which the values of integral functionals will be calculated and solutions of IEs will be built.Additionally, the convergence of the calculation process and accuracy of the obtained results will be checked at these points.It is convenient to place these points in the middle between the node points, and then their coordinates are determined similarly to Formula (13): and with a sufficiently large number of the contour breakdowns on AFEs, it may be counted even simpler: ξ j,i = (x j,i + x j,i+1 )/2 and η j,i = (y j,i + y j,i+1 )/2 .(17)

Organizing the Computation of Integral Operator Kernels
The calculation of the kernels of integral operators is performed by Formulae (2) and (3).All of the variable values necessary for their use are determined according to procedures outlined in the previous section.Here we point out some features of their applications and a form of the database in which calculation results are arranged.
All variables in the given expressions contain two indices: the number of the approximating arc j and the number of a point on this arc i.Therefore, in the process of computing, the values of each variable are grouped into M linear matrices (row vectors) generally of different lengths (with a different number of elements).
At the same time, the result of calculating Formulae (2) and ( 3) is a value related to any two points P and M selected on the contour, which means that it is naturally supplied with two indexes-the numbers of these points kP and kM according to numbering on the contour.Therefore, the indexation of variable values must be transformed from two-digit to single-digit-that is, one vector must be built from M vectors by combining them in the required sequence.This procedure is called concatenation, and in the MATLAB environment it is performed as a usual matrix operation [34].

Pattern of the Procedure for the Formation of One-Dimensional (Linear) Arrays of Nodal Points Coordinates
To clarify the approach to the implementation of this procedure, we give a fragment of the program code with arrays named XNODEJ and YNODEJ (see Table 3):  Recall that "arcNumber" is the total number of approximating arcs (M), and "arcNDel(in)" is the number of partitions of the arc number "in" on AFE; the calculations in the internal cycle correspond to Formula (14).
Linear matrices of the coordinates of collocation points (XCOLLJ and YCOLLJ-Formula (16)) and complex values of the normals (UNITNORM-Formula (15)) are formed in a completely analogous way.
To ensure the flexibility (adaptability) of integral operators' computing to different classes of physical phenomena, in addition to the specified arrays we also enter the matrices RCURVE and DELARC to the database of the geometrical platform, which contain the values of the curvature radius and the length of the elementary arc (AFE) for each point according to numbering along the contour.
The size of the matrices corresponds to the total number of node points on the contour N.

Pattern of Procedure for the Formation of the Database for Solving Integral Equations
The last stage of forming of the database for the solution of integral equations is computing a square matrix of SLAE coefficients, by which the IEs are replaced.It is most convenient to carry out this stage of preparation of the simulation of complex physical phenomena in CECT by a separate procedure, which in the MATLAB environment is made out in the form of the user function.Its listing is given in Table 4 below.
This function has two input arguments: the numbers kP and kM of points P and M on the contour, and one output argument "Kernel", which transmits the calculated kernel value for these points to the main program.The filling of the matrix of SLAE coefficients of N × N size is performed by calculating the function "yadro" in a double cycle by the numbers of points P and M.

Examples and Estimates of Efficiency
Testing of the developed procedures for the formation of the geometrical platform and corresponding databases is done on a control contour in the form of a circumference with radius R = 1, divided into six even arcs, which is presented in Table 1.The number of arc partitions was varied from 5 to 50.Graphically, the results of the calculation are shown in Figure 4a (i.e., the distribution of nodal and collocation points (120)) and Figure 4b (i.e., the distribution of normal vectors).In all variants of partitioning, the sum of AFEs was equal to circumference length (2π) with accuracy not less than 10 characters.This function has two input arguments: the numbers kP and kM of points P and M on the contour, and one output argument "Kernel", which transmits the calculated kernel value for these points to the main program.The filling of the matrix of SLAE coefficients of N × N size is performed by calculating the function "yadro" in a double cycle by the numbers of points P and M.

Examples and Estimates of Efficiency
Testing of the developed procedures for the formation of the geometrical platform and corresponding databases is done on a control contour in the form of a circumference with radius R = 1, divided into six even arcs, which is presented in Table 1.The number of arc partitions was varied from 5 to 50.Graphically, the results of the calculation are shown in Figure 4a (i.e., the distribution of nodal and collocation points (120)) and 4b (i.e., the distribution of normal vectors).In all variants of partitioning, the sum of AFEs was equal to circumference length (2π) with accuracy not less than 10 characters.The functionality of the submitted geometrical platform was also evaluated on more difficult configurations: a hexagonal figure with arches of different curvature and camber (Figure 5), as well as a trigonous figure (Figure 6).These results indicate full suitability of the developed tool for the solution of tasks with a multilink boundary configuration.
For the circumference represented by six equal arcs, all components of the geometrical platform specified in the text were calculated, including the coordinates of the points, the directing cosines, and the arc finite elements (AFEs).The figures show the absolute (with the 0th error) The functionality of the submitted geometrical platform was also evaluated on more difficult configurations: a hexagonal figure with arches of different curvature and camber (Figure 5), as well as a trigonous figure (Figure 6).These results indicate full suitability of the developed tool for the solution of tasks with a multilink boundary configuration.As another criterion for the quality of approximation of different figures by means of the geometrical platform, authors accepted the extent of coincidence of the sum of all AFEs with the theoretically exact length of a contour.In particular, for a circumference with different amounts of AFEs (from 30 to 300) even at the roughest discretization (30 AFEs), the values coincided up to the tenth digit inclusively.Table 5 demonstrates the same results for other forms of contour.In fact, this testing was carried out for a dozen different configurations, close to the outlines of real objects, but are not presented here for the sake of brevity.In all cases, the accuracy of reproducing objects' geometrical parameters was more than sufficient.The functional validity of the geometrical platform (GP) for execution of the main taskcalculating the kernel of integral operators-was tested on a number of contours, ranging from a As another criterion for the quality of approximation of different figures by means of the geometrical platform, authors accepted the extent of coincidence of the sum of all AFEs with the theoretically exact length of a contour.In particular, for a circumference with different amounts of AFEs (from 30 to 300) even at the roughest discretization (30 AFEs), the values coincided up to the tenth digit inclusively.Table 5 demonstrates the same results for other forms of contour.In fact, this testing was carried out for a dozen different configurations, close to the outlines of real objects, but are not presented here for the sake of brevity.In all cases, the accuracy of reproducing objects' geometrical parameters was more than sufficient.The functional validity of the geometrical platform (GP) for execution of the main taskcalculating the kernel of integral operators-was tested on a number of contours, ranging from a For the circumference represented by six equal arcs, all components of the geometrical platform specified in the text were calculated, including the coordinates of the points, the directing cosines, and the arc finite elements (AFEs).The figures show the absolute (with the 0th error) coincidence of the normals to the arcs with the directions of the radii of arcs' curvature for all contours.

Kind of contour
As another criterion for the quality of approximation of different figures by means of the geometrical platform, authors accepted the extent of coincidence of the sum of all AFEs with the theoretically exact length of a contour.In particular, for a circumference with different amounts of AFEs (from 30 to 300) even at the roughest discretization (30 AFEs), the values coincided up to the tenth digit inclusively.Table 5 demonstrates the same results for other forms of contour.In fact, this testing was carried out for a dozen different configurations, close to the outlines of real objects, but are not presented here for the sake of brevity.In all cases, the accuracy of reproducing objects' geometrical parameters was more than sufficient.
The functional validity of the geometrical platform (GP) for execution of the main task-calculating the kernel of integral operators-was tested on a number of contours, ranging from a circle to curved acute-angled triangles.In particular, we obtained data characterizing the accuracy of kernel value calculation and the duration of computation for a circle, with respect to which the exact theoretical value is known: K 2 (P, M) = 1/2R c .Accordingly, for R c = 1 the kernel K 2 (P, M) = 0.5.The results of calculating kernel values, accuracy, and duration of computing for different numbers of AFEs are presented in Table 6.Now, it is appropriate to evaluate the general options of the obtained database and to compare them with those that characterize the base for solving similar problems by means of differential equations (DEs) using FEM.As a reference point for comparison, the parameters of test contour may be taken from six even arcs, which is presented in Listing 1-a circle of radius R = 1, broken down into 360 AFEs (through 1 • ).Then, the number of sought-for values of variables at the collocation points will be N = 360, and the SLAE matrix will have 360 × 360 = 129,600 elements.
To confirm and illustrate the full functional validity of GP processing for kernel calculation, a visual representation of a kernel values array is shown in Figure 7 for two contours (i.e., Figures 5 and 6).
Data 2019, 4, 136 13 of 18 circle to curved acute-angled triangles.In particular, we obtained data characterizing the accuracy of kernel value calculation and the duration of computation for a circle, with respect to which the exact theoretical value is known: K2 (P, M) = 1/2Rc.Accordingly, for Rc = 1 the kernel K2(P, M) = 0.5.The results of calculating kernel values, accuracy, and duration of computing for different numbers of AFEs are presented in Table 6.Now, it is appropriate to evaluate the general options of the obtained database and to compare them with those that characterize the base for solving similar problems by means of differential equations (DEs) using FEM.As a reference point for comparison, the parameters of test contour may be taken from six even arcs, which is presented in Listing 1-a circle of radius R = 1, broken down into 360 AFEs (through 1°).Then, the number of sought-for values of variables at the collocation points will be N = 360, and the SLAE matrix will have 360 × 360 = 129,600 elements.
To confirm and illustrate the full functional validity of GP processing for kernel calculation, a visual representation of a kernel values array is shown in Figure 7 for two contours (i.e., Figures 5  and 6).
Therefore, the obtained results give the authors all grounds to claim that the submitted geometrical platform is a reliable software tool of procedures that are entirely suitable for solving integral equations of the considered types.If for a task with differential equations we consider the computational domain in the form of a square equivalent in area to the mentioned circle with a side 1772 π = , breaking its perimeter into the same AFEs as the circle, we get 406 contour elements, and the points of collocations where it is necessary to find the variables' values (i.e., the order of the equations system) will be approximately (406/2) 2 ≈ 41,000.The order of the coefficient matrix for such a system will reach more than 1 billion elements.Certainly, such systems are not solved directly, but are reduced to sparse matrices.However, the counting duration and the accumulation of errors (residuals) of their solution increase very significantly.Computing such systems of equations using a computer with a 2-4-core processor can take more than 5-6 h.Therefore, the obtained results give the authors all grounds to claim that the submitted geometrical platform is a reliable software tool of procedures that are entirely suitable for solving integral equations of the considered types.
If for a task with differential equations we consider the computational domain in the form of a square equivalent in area to the mentioned circle with a side √ π = 1772, breaking its perimeter into the same AFEs as the circle, we get 406 contour elements, and the points of collocations where it is necessary to find the variables' values (i.e., the order of the equations system) will be approximately (406/2) 2 ≈ 41,000.The order of the coefficient matrix for such a system will reach more than 1 billion elements.Certainly, such systems are not solved directly, but are reduced to sparse matrices.However, the counting duration and the accumulation of errors (residuals) of their solution increase very significantly.Computing such systems of equations using a computer with a 2-4-core processor can take more than 5-6 h.
To estimate the temporal characteristics of the presented geometric platform, a timer was built into the program code, which captures the computing time.As a result of the test, the above-mentioned matrix of SLAE coefficients of order 360 was filled for 5 s, and its 300-fold repetition lasted 4 min and 49 s.Based on the quadratic trend of time growth from the order of the matrix (i.e., SLAE), it can be expected that computing a system of order of about 1000 would take 50-60 s.
Moreover, the authors executed an extended test of the entire package of GP program blocks, applying it directly to solve the integral equations on the contours, the outlines of which are close to actual cross-sections of metallurgical ladles.In particular, two examples of field problem solving by IEM on the basis of submitted GP were considered: section of type "keg" (8a) and section of type "jug" (9a), placed in an external uniform exciting field.The solutions of equations were found as the distribution of the normal component of the resulting field along the contours of the conditional ladle section.In Figures 8 and 9, the plots of mentioned field distribution are shown in comparison with the exciting field.Authors must remark that the presented examples are only illustrative, but there are no principal complications to expanding this approach for other CECT tasks.To estimate the temporal characteristics of the presented geometric platform, a timer was built into the program code, which captures the computing time.As a result of the test, the abovementioned matrix of SLAE coefficients of order 360 was filled for 5 s, and its 300-fold repetition lasted 4 min and 49 s.Based on the quadratic trend of time growth from the order of the matrix (i.e., SLAE), it can be expected that computing a system of order of about 1000 would take 50-60 s.
Moreover, the authors executed an extended test of the entire package of GP program blocks, applying it directly to solve the integral equations on the contours, the outlines of which are close to actual cross-sections of metallurgical ladles.In particular, two examples of field problem solving by IEM on the basis of submitted GP were considered: section of type "keg" (8a) and section of type "jug" (9a), placed in an external uniform exciting field.The solutions of equations were found as the distribution of the normal component of the resulting field along the contours of the conditional ladle section.In Figures 8 and 9, the plots of mentioned field distribution are shown in comparison with the exciting field.Authors must remark that the presented examples are only illustrative, but there are no principal complications to expanding this approach for other CECT tasks.Based on the results, the authors may declare with full responsibility that the submitted geometric platform as an ordered set of program codes is completely ready to be used for solving manifold tasks connected with CECT.
Therefore, this approach will be applied in future research dealing with the practical solution, using the method of integral equations (MIE), of real problems during the CECT of liquid metal in genuine ladles based on the submitted geometric platform for big data computing.To estimate the temporal characteristics of the presented geometric platform, a timer was built into the program code, which captures the computing time.As a result of the test, the abovementioned matrix of SLAE coefficients of order 360 was filled for 5 s, and its 300-fold repetition lasted 4 min and 49 s.Based on the quadratic trend of time growth from the order of the matrix (i.e., SLAE), it can be expected that computing a system of order of about 1000 would take 50-60 s.
Moreover, the authors executed an extended test of the entire package of GP program blocks, applying it directly to solve the integral equations on the contours, the outlines of which are close to actual cross-sections of metallurgical ladles.In particular, two examples of field problem solving by IEM on the basis of submitted GP were considered: section of type "keg" (8a) and section of type "jug" (9a), placed in an external uniform exciting field.The solutions of equations were found as the distribution of the normal component of the resulting field along the contours of the conditional ladle section.In Figures 8 and 9, the plots of mentioned field distribution are shown in comparison with the exciting field.Authors must remark that the presented examples are only illustrative, but there are no principal complications to expanding this approach for other CECT tasks.Based on the results, the authors may declare with full responsibility that the submitted geometric platform as an ordered set of program codes is completely ready to be used for solving manifold tasks connected with CECT.
Therefore, this approach will be applied in future research dealing with the practical solution, using the method of integral equations (MIE), of real problems during the CECT of liquid metal in genuine ladles based on the submitted geometric platform for big data computing.Based on the results, the authors may declare with full responsibility that the submitted geometric platform as an ordered set of program codes is completely ready to be used for solving manifold tasks connected with CECT.
Therefore, this approach will be applied in future research dealing with the practical solution, using the method of integral equations (MIE), of real problems during the CECT of liquid metal in genuine ladles based on the submitted geometric platform for big data computing.

Conclusions
The management of casting quality in modern foundry production has paramount importance for which, in particular, the electric current treatment of liquid metal is used.In this process, the interaction of electric current with the melt is accompanied by the emergence of a complex of interrelated physical fields of different nature-electromagnetic, thermal, hydrodynamic, wave and mechanical fields, which eventually predetermine the quality of castings.However, the totality of physical phenomena in this process is largely unavailable for direct observation and measurement.
Therefore, a key role in the choice of the most appropriate treatment parameters for the directed quality control of castings is played by the mathematical modeling and simulation of specified complexes of physical fields-primarily electromagnetic.The theoretical basis for studying the interaction with the medium (liquid and solidifying metal) of the listed fields are the known differential equations in various interpretations.At present, traditional approaches to their solution in the context of CECT problems are based mainly on use of the finite element method and its modifications.
However, the complicated configuration of computational areas and the uncertainty of boundary conditions in boundary value problems often lead to the fact that a very significant degree of discretization of the DE solution area is required in order to ensure the reliability of results.The number of FEs and, accordingly, the number of sought-for variable values is a very large, and may be more than one hundred thousand.At the same time, the database for preparing and executing convergent decision-making procedures according to FEM must have at least a quadratic trend relative to dimensionality of the problem (i.e., it will exceed 1 billion components).In another case, the algorithmic optimization of the database dimension with the use of disperse matrices causes an increase in the volume of computational operations and in the duration of counting in geometrical progression.A standard computer with a 4-core processor spends 5-6 h to solve such tasks.
The main contribution of this work is the development of an alternative methodological approach to formulating problems of mathematical modeling and complex physical phenomenon simulation during the CECT of liquid metals.The approach developed by the authors is guided by the use of integral formulations of fundamental solutions of field theory problems in the form of Green's formulae and connected relationships, which are generalized in the method of integral equations (MIE).
The solution of tasks in this case is built on the boundary set, and as a result the spatial dimensionality of the original (primary) problem is reduced by at least one order of magnitude and the number of sought unknowns is reduced by almost by two orders.
For transition of the task formulation from DE to IE and the equivalent SLAE, as well as the success of their solution, the procedure of rationally representing the geometric parameters of computational area boundaries-accounting for the need of their explicit use in the algorithmic structure of the mathematical operations of differentiation and integration on curvilinear surfaces as well as the calculation of two-point operators in the integral equations-has great methodological importance.
In other words, in these circumstances, necessity appeared to form an adaptive and concise geometric platform for computing big databases in problems of studying complex physical phenomena during the CECT of liquid metals.The authors obtained the fundamental results in solving this problem by: (a) developing a mathematical model for representing intricate configurations of computational areas by using fairly simple procedures, (b) developing a scheme for the organization of computational processes, and (c) implementing the geometrical platform in the form of program code that may provide big-database computing for the whole variety of problems arising from CECT tasks.
Finally, the computational efficiency of the created platform can be characterized by the following comparison: for the boundary contour representing a ladle section with melt shown in Figure 2, the order of the equations system (i.e., the number of unknown values) will be (depending on the required accuracy) 500-600, compared to 100,000 in the solution by FE method.Naturally, the computing time is reduced to several minutes instead of 5 h.
The functionality of the proposed geometrical platform was evaluated on complicated configurations: a hexagonal figure with arches of different curvature and camber, and a trigonous configuration.These results show the full suitability of the developed tool for the solution of tasks with a multilink configuration of boundaries.
Future research will deal with the use of MIE for the practical solution of real problems during the CECT of liquid metal in genuine ladles based on the submitted geometric platform for big data computing.

Figure 1 .
Figure 1.Principal constructive scheme of the conductive electric current treatment (CECT).

Figure 2 .
Figure 2. Ladle contour approximated by straight lines and arcs.

Figure 1 .
Figure 1.Principal constructive scheme of the conductive electric current treatment (CECT).

Figure 1 .
Figure 1.Principal constructive scheme of the conductive electric current treatment (CECT).

Figure 2 .
Figure 2. Ladle contour approximated by straight lines and arcs.

Figure 2 .
Figure 2. Ladle contour approximated by straight lines and arcs.

Figure 3 .Figure 3 .
Figure 3.For determination of the main geometrical parameters of a separate arc.

Figure 4 .
Figure 4. Distribution of points and normals on a circumference: (a) distribution of nodal and collocation points; (b) distribution of normal vectors.

Figure 4 .
Figure 4. Distribution of points and normals on a circumference: (a) distribution of nodal and collocation points; (b) distribution of normal vectors.

Figure 5 .Figure 6 .
Figure 5. Distribution of points and normals on a hexagonal contour: (a) distribution of nodal and collocation points on a contour; (b) distribution of normal vectors on a contour.

Figure 5 .Figure 5 .Figure 6 .
Figure 5. Distribution of points and normals on a hexagonal contour: (a) distribution of nodal and collocation points on a contour; (b) distribution of normal vectors on a contour.

Figure 6 .
Figure 6.Distribution of points and normals on a trigonous contour: (a) distribution of nodal and collocation points on a contour, (b) distribution of normal vectors on a contour.

Figure 7 .
Figure 7. Plots of kernel values array: (a) of the contour in Figure 5 with 30 AFEs; (b) of the contour in Figure 6 with 15 AFEs.

Figure 7 .
Figure 7. Plots of kernel values array: (a) of the contour in Figure 5 with 30 AFEs; (b) of the contour in Figure 6 with 15 AFEs.

Table 1 .
The contour input data matrix.

Table 2 .
Calculation of the chord inclination angle.

Table 3 .
Formation of arrays of nodal points coordinates with names XNODEJ and YNODEJ.

Table 4 .
Computation of a square matrix of system of linear algebraic equations (SLAE) coefficients.

Table 5 .
Comparison of calculated contour length with theoretical value.

Table 5 .
Comparison of calculated contour length with theoretical value.

Table 5 .
Comparison of calculated contour length with theoretical value.

Table 6 .
Kernel values and duration of computing.

Table 6 .
Kernel values and duration of computing.