Thermodynamic Fundamentals of Cellular Automata Model of the Process of Solidification of Metals and Alloys Considering the Phase Transition Tatyana Selivorstovaa , Vadim Selivorstova , Anton Gudaa and Katerina Ostrovskaa a The National Metallurgical Academy of Ukraine, Gagarina avenue 4, Dnipro, Ukraine, 49005 Abstract We considered the mathematical apparatus that can be used to calculate the thermal problem associated with the features of solidification of metals and alloys. Thermodynamic laws are presented, that made it possible to correctly take into account the phase transition, a feature of which is to release latent heat. A one-dimensional problem and its generalization to a cylindrical coordinate system are considered. Keywords Cellular automata, thermodynamic model, temperature, amount of liquid phase, solidification, phase transition 1. Introduction In industrial process control systems, including foundries, embedded systems based on microcon- trollers are widely used. Unfortunately, direct simulation on microcontrollers is impossible due to limited computing resources. Using cellular automata approach allows simplified calculations on embedded systems. Description of real metallurgical processes associated with solidification [1, 2] leads to signif- icant differences from the classical boundary value problem due to many factors [3, 4], the most important of which are due to: heat release during the phase transition [5], which is described by a complex diagram states [6]; heat transfer due to convective flows in the liquid phase [7, 8]; taking into account the real boundary conditions of heat transfer [9]; inhomogeneity of alloy properties [10, 11]; the dependence of thermophysical parameters on temperature [12]; the complexity of the geometry of the product [13, 14]. The analysis of methods for modelling metallurgical processes showed that the existing analytical solutions of the solidification problem provide a high accuracy of the solution [15, 16], but can be used only in some simple cases and are rather methodological in nature [17, 18]. The description of the solidification problem in the form of partial differential equations leads to the use of numerical solution methods [19, 20], for which a wide range of modifications has been developed [19]. They are united by the presence of artificial methods, such as catching ICT&ES-2020: Information-Communication Technologies & Embedded Systems, November 12, 2020, Mykolaiv, Ukraine " tatyanamikhaylovskaya@gmail.com (T. Selivorstova); seliverstovvy@gmail.com (V. Selivorstov); atu.guda@gmail.com (A. Guda); kuostrovskaya@gmail.com (K. Ostrovska)  0000-0002-2470-6986 (T. Selivorstova); 0000-0002-1916-625X (V. Selivorstov); 0000-0003-1139-1580 (A. Guda); 0000-0002-9375-4121 (K. Ostrovska) Β© 2020 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings http://ceur-ws.org ISSN 1613-0073 CEUR Workshop Proceedings (CEUR-WS.org) the front at the grid node, straightening the fronts, successively approximating the position of the transition phase, smoothing the phase transition, introducing the effective heat capacity [21] and then determining the boundaries of the two-phase zone by interpolation, which lead to a distortion of the original physical setting of the considered problem [22, 23, 24]. The use of thermodynamic models [25] presupposes the existence of a large-scale database of experimental data [26] on the redistribution of components between phases, the dependence of the composition and number of precipitated phases on temperature changes under conditions of equilibrium and nonequilibrium crystallization [27]. Thus, the development of mathematical models of metallurgical processes with phase tran- sitions is an urgent task, which is caused by the necessity of improvements of the quality of products of metallurgical production [28, 29, 30] and their cost reduction [31, 32]. The use of cellular automata in the development of mathematical models of solidification [33, 34, 35] makes it possible to describe various nonlinear processes immediately in a discrete language [36, 37] and has a number of advantages associated with the possibility of organizing high-speed parallel computations [38, 39, 40], the obviousness of algorithms, the possibility of using them to describe processes that are difficult or even impossible to describe by partial differential equations. 2. Cellular automata solidification model Solidification of metals and alloys is accompanied by complex physical and chemical processes in the melt. The most important of them are due to: heat release during the phase transition; heat transfer due to convective flows in the liquid phase; the complexity of the real boundary conditions of heat transfer; heterogeneity of alloy properties; dependence of thermophysical parameters on temperature; the complexity of the geometry of the product [41]. Thermodynamic phenomena accompanying the process of solidification of metals is described by the Fourier equation [42] π‘˜π‘† β–³ 𝑇 △𝑄 = 𝜏, (1) β„Ž where △𝑄 [𝐽 ] – the amount of heat passed through the surface 𝑆 [π‘š2 ] for the time 𝜏 [𝑠] with temperature differences △𝑇 [𝐾 ], π‘˜ [π‘Š /(π‘šπΎ )] – – coefficient of thermal conductivity, β„Ž [π‘š] – layer thickness. Equilibrium Specific Energy π‘ž [𝐽 /π‘˜π‘”] is uniquely related to the temperature and phase state of the unit volume 𝑄 π‘ž= , (2) πœŒβ„Žπ‘† where 𝜌 [π‘˜π‘”/π‘π‘š3 ] – density. 2.1. One-dimensional cellular automata thermodynamic model The cellular automaton mathematical model assumes the description of nonlinear processes in a discrete language at once and allows avoiding the description of thermodynamic processes in the form of partial differential equations [43]. To construct a one-dimensional thermodynamic model of solidification, taking into account the phase transition, it is necessary to determine the rules of heat transfer between cells. Consider a cell with a Neumann neighbourhood for the one-dimensional case (Fig. 1). Figure 1: Schematic representation of a group of cells and symbols for the one-dimensional case Fig. 1 notation accepted: π‘‡π‘–βˆ’1 , 𝑇𝑖 , 𝑇𝑖+1 – temperatures in 𝑖 βˆ’ 1, 𝑖, 𝑖 + 1 adjacent cells, respectively, β„Ž – distance between cell centers, π‘„π‘–βˆ’1,𝑖 , 𝑄𝑖+1,𝑖 – energy entering the cell from 𝑖 βˆ’ 1 and 𝑖 + 1 respectively. Energy entering the cell 𝑖 from 𝑖 βˆ’ 1 cell, is determined based on the replacement in (1) of the differential by the difference π‘˜π‘–,π‘–βˆ’1 𝑆(π‘‡π‘–βˆ’1 βˆ’ 𝑇𝑖 )𝜏 π‘„π‘–βˆ’1,𝑖 = , (3) β„Ž where π‘˜π‘–,π‘–βˆ’1 – effective coefficient of thermal conductivity, determined from the equality of heat flows at the border of neighboring cells. A schematic representation of the heat flow motion for the one-dimensional case is shown on Fig.2. Figure 2: Diagram of heat flow movement between cells Heat flow 𝑃 [𝐽 /𝑠] (fig. 2), corresponding to the amount of energy transmitted through the isothermal surface per unit of time is determined by the ratio π‘˜π‘†π‘‡ 𝑃= . (4) β„Ž Heat flow between 𝑖 and 𝑖 βˆ’ 1 cells is equal to π‘˜π‘–βˆ’1,𝑖 𝑆(𝑇𝑖 βˆ’ π‘‡π‘–βˆ’1 ) π‘ƒπ‘–βˆ’1,𝑖 = . β„Ž Heat flow to the border that separates cells: β€’ from 𝑖 βˆ’ 1 side π‘˜π‘–βˆ’1 𝑆(𝑇 βˆ’ π‘‡π‘–βˆ’1 ) π‘ƒπ‘–βˆ’1 = , β„Ž/2 β€’ from 𝑖 cell side π‘˜π‘– 𝑆(𝑇𝑖 βˆ’ 𝑇 ) 𝑃𝑖 = , β„Ž/2 where 𝑇 – the current temperature at the border of the considered cells. Based on the equality of heat flows, at the border of the cells π‘ƒπ‘–βˆ’1 = 𝑃𝑖 , we can find 𝑇 – temperature on the border between them: π‘˜π‘–βˆ’1 π‘‡π‘–βˆ’1 + π‘˜π‘– 𝑇𝑖 𝑇 = . π‘˜π‘–βˆ’1 + π‘˜π‘– Equating heat flows π‘ƒπ‘–βˆ’1 = 𝑃𝑖 , and we obtain the value of the effective thermal conductivity coefficient 2π‘˜π‘– βˆ’ 1π‘˜π‘– π‘˜π‘–βˆ’1,𝑖 = . (5) π‘˜π‘– βˆ’ 1 + π‘˜π‘– So the change in energy △𝑄 for 𝑖 cells, based on (3), is determined by the relation πœπ‘† △𝑄 = (π‘˜π‘–βˆ’1,𝑖 (π‘‡π‘–βˆ’1 βˆ’ 𝑇𝑖 ) + π‘˜π‘–,𝑖+1 (𝑇𝑖+1 βˆ’ 𝑇𝑖 )). (6) β„Ž The above approach for calculating the change in cell energy can be extended to the case of two-dimensional and three-dimensional problems in Cartesian coordinates. 2.2. Three-dimensional cylindrical cellular automata thermodynamic model A significant part of technical objects has a cylindrical shape, that allows using the property of axial symmetry to reduce the dimension of the model and reduce the amount of computation [44]. Let’s consider the motion of heat flows for a group of cells with a Neumann neighbourhood in the presence of axial symmetry of the simulation object. Fig. 3 schematically shows a cell in a cylindrical coordinate system. Before determining the change in cell energy, it is necessary to determine the areas of the lateral surfaces of the cell. Figure 3: Three-dimensional image of a cell Let’s introduce the notation: 𝑖 – cell number in the radial direction, 𝑖 = 1...𝑖𝑛 ; 𝑗 – sectoral cell number, 𝑗 = 1...𝑗𝑛 ; π‘˜ – cell number along the axis of symmetry, π‘˜ = 1...π‘˜π‘› ; β„Žπ‘– – radial cell length; β„Žπ‘— – the sectoral width of the cell depends on the position of the cell in the radial direction; β„Žπ‘˜ – cell height along the axis of symmetry; 𝑛𝑖 – number of cells in the radial direction; 𝑛𝑗 – sectoral number of cells; π‘›π‘˜ – number of cells along the axis of symmetry. Fig. 4, 5, 6 show a schematic representation of a cell (𝑖, 𝑗, π‘˜) and notations of its lateral sides. Figure 4: Schematic representation of a cell and notations of lateral sides The sectoral width of the cell depends on the position of the cell in the radial direction, i.e. from 𝑖: 2πœ‹β„Žπ‘– 𝑖 β„Žπ‘— (𝑖) = . 𝑛𝑗 The lateral sides of a cell are expressed in terms of the cell number, cell length and the number of cells along the axes: 𝑆𝑖 = β„Žπ‘– β„Žπ‘˜ , 2πœ‹β„Žπ‘– 𝑖 𝑆𝑗 (𝑖) = β„Žπ‘˜ , 𝑛𝑗 πœ‹β„Ž2𝑖 (2𝑖 βˆ’ 1) π‘†π‘˜ (𝑖) = . 𝑛𝑗 Figure 5: Schematic representation of cells in a cylindrical coordinate system (top view) Figure 6: Schematic representation of cells in a cylindrical coordinate system (side view) The change in cell energy for a cylindrical coordinate system consists of the sums of energy changes in the radial, sectoral directions and along the axis of symmetry of the body. Thus, the change in energy for 𝑖, 𝑗, π‘˜β€“th cell β€’ in radial direction: 𝑗,π‘˜ 𝑗,π‘˜ π‘˜π‘–βˆ’1,𝑖 𝑆𝑗 (𝑖) [π‘‡π‘–βˆ’1 βˆ’ 𝑇𝑖 ] 𝜏 π‘„π‘–βˆ’1,𝑖 = , β„Žπ‘– 𝑗,π‘˜ 𝑗,π‘˜ π‘˜π‘–+1,𝑖 𝑆𝑗 (𝑖 + 1) [𝑇𝑖+1 βˆ’ 𝑇𝑖 ] 𝜏 𝑄𝑖+1,𝑖 = ; β„Žπ‘– β€’ in the sectoral direction: 𝑖,π‘˜ 𝑖,π‘˜ π‘˜π‘—βˆ’1,𝑗 𝑆𝑗 (𝑖) [π‘‡π‘—βˆ’1 βˆ’ 𝑇𝑗 ] 𝜏 π‘„π‘—βˆ’1,𝑗 = , β„Žπ‘— (𝑖) 𝑖,π‘˜ 𝑖,π‘˜ π‘˜π‘—+1,𝑗 𝑆𝑗 (𝑖 + 1) [𝑇𝑖+1 βˆ’ 𝑇𝑖 ] 𝜏 𝑄𝑗+1,𝑗 = ; β„Žπ‘— (𝑖) β€’ along the axis of symmetry: 𝑖,𝑗 𝑖,𝑗 π‘˜π‘˜βˆ’1,π‘˜ π‘†π‘˜ (𝑖) [π‘‡π‘˜βˆ’1 βˆ’ π‘‡π‘˜ ] 𝜏 π‘„π‘˜βˆ’1,π‘˜ = , β„Žπ‘˜ 𝑖,𝑗 𝑖,𝑗 π‘˜π‘˜+1,π‘˜ π‘†π‘˜ (𝑖 + 1) [π‘‡π‘˜+1 βˆ’ π‘‡π‘˜ ] 𝜏 π‘„π‘˜+1,π‘˜ = . β„Žπ‘˜ In the above formulas, superscripts indicate the same coordinates for neighbouring cells, and subscripts indicate the direction of the considered interaction. A cell energy change is the sum of energy changes in all directions: βˆ‘ 𝑄𝑖,𝑗,π‘˜ = (7) In the above formulas, the unknown parameter is the effective thermal conductivity, that is calculated from the equality of heat flows between cells for a cylindrical coordinate system (Fig. 7). Figure 7: Heat flows between cells for radial heat distribution Heat flow between cells 𝑖 and 𝑖 + 1 is determined by ratio π‘˜π‘–,𝑖+1 (𝑇𝑖 βˆ’ 𝑇𝑖+1 ) 2πœ‹β„Žπ‘– β„Žπ‘˜ 𝑖 𝑃𝑖,𝑖+1 = . β„Žπ‘– 𝑛𝑗 Heat flow to the border that separates cells: β€’ from 𝑖 cell’s side π‘˜π‘– (𝑇𝑖 βˆ’ 𝑇 ) 2πœ‹β„Žπ‘– β„Žπ‘˜ (𝑖 βˆ’ 0.25) 𝑃𝑖 = , β„Žπ‘– /2 𝑛𝑗 β€’ from 𝑖 + 1 cell’s side π‘˜π‘–+1 (𝑇 βˆ’ 𝑇𝑖+1 ) 2πœ‹β„Žπ‘– β„Žπ‘˜ (𝑖 βˆ’ 0.25) 𝑃𝑖+1 = . β„Žπ‘– /2 𝑛𝑗 Based on their equality of heat flows at the cell border 𝑃𝑖 = 𝑃𝑖+1 we can find 𝑇 – cell border temperature: π‘˜π‘– 𝑇𝑖 (𝑖 βˆ’ 0.25) + π‘˜π‘–+1 𝑇𝑖+1 (𝑖 + 0.25) 𝑇 = . (8) π‘˜π‘– (𝑖 βˆ’ 0.25) + π‘˜π‘–+1 (𝑖 + 0.25) Equate 𝑃𝑖,𝑖+1 = 𝑃𝑖 and we obtain the value of the effective coefficient of thermal conductivity for radial heat propagation 2π‘˜π‘– (𝑖 βˆ’ 0.25)π‘˜π‘–+1 (𝑖 + 0.25) π‘˜π‘–βˆ’1,𝑖 = . (9) 𝑖(π‘˜π‘– (𝑖 βˆ’ 0.25) + π‘˜π‘–+1 (𝑖 + 0.25)) When calculating the change in energy along the axis of symmetry and along the circum- ference, the effective coefficient of thermal conductivity is calculated similarly to the one- dimensional case using the formula (5). 3. Conclusions The presented mathematical model can be integrated into an automated production control system, the components of which communicate in order to optimise the parameters of the production cycle. We have given a description of a cellular-automaton thermodynamic model of solidification, taking into account the phase transition for the Cartesian and cylindrical coordinate systems, which differs from the existing ones by correctly considering the thermodynamic features of the solidification process. References [1] T. Selivyorstova, A. Mikhalyov, Mathematical model of the two-phase zone supply of solidified metal castings under the influence of adjustable gas pressure, in: 2019 9th Inter- national Conference on Advanced Computer Information Technologies (ACIT), IEEE, 2019. URL: https://doi.org/10.1109/acitt.2019.8779914. doi:10.1109/acitt.2019.8779914. [2] T. Selivyorstova, A. Mikhalyov, Analysis of prediction mathematical models of shrinkage defects in castings, in: 2018 IEEE First International Conference on System Analysis & Intelligent Computing (SAIC), IEEE, 2018. URL: https://doi.org/10.1109/saic.2018.8516811. doi:10.1109/saic.2018.8516811. [3] V. Kalinin, V. Khrychikov, V. Krivosheyev, V. Y. Seliverstov, Y. V. Dotsenko, A. Kondrat, Advanced technologies of cast iron complex alloying and inoculation for mining and smelting equipment parts casting, Metallurgical and Mining Industry 2 (2010) 13. [4] K. Borodianskiy, V. Selivorstov, Y. Dotsenko, M. Zinigrad, Effect of additions of ceramic nanoparticles and gas-dynamic treatment on al casting alloys, Metals 5 (2015) 2277–2288. URL: https://doi.org/10.3390/met5042277. doi:10.3390/met5042277. [5] L. Du, Y. Cao, Z. Gao, H. Du, Phase-field modelling and simulation of morphology transition induced by thermoelectric effects during solidification in pure metals, Materials Research Express 5 (2018). URL: https://doi.org/10.1088/2053-1591/aad48c. doi:10.1088/ 2053-1591/aad48c. [6] C. Li, Z. Yu, J. Gao, J. Zhao, X. Han, Numerical simulation and experimental study of cladding fe60 on an ASTM 1045 substrate by laser cladding, Surface and Coatings Technology 357 (2019) 965–977. URL: https://doi.org/10.1016/j.surfcoat.2018.10.099. doi:10. 1016/j.surfcoat.2018.10.099. [7] Y. Jaluria, K. E. Torrance, Computational Heat Transfer, Routledge, 2017. URL: https: //doi.org/10.1201/9781315140018. doi:10.1201/9781315140018. [8] S. Chakraborty, P. Dutta, Three-dimensional double-diffusive convection and macrosegre- gation during non-equilibrium solidification of binary mixtures, International Journal of Heat and Mass Transfer 46 (2003) 2115–2134. URL: https://doi.org/10.1016/s0017-9310(02) 00530-6. doi:10.1016/s0017-9310(02)00530-6. [9] J. D. Miller, T. M. Pollock, The effect of processing conditions on heat transfer during directional solidification via the bridgman and liquid metal cooling processes, Metal- lurgical and Materials Transactions A 45 (2013) 411–425. URL: https://doi.org/10.1007/ s11661-013-1986-2. doi:10.1007/s11661-013-1986-2. [10] V. Selivorstov, Y. Dotsenko, K. Borodianskiy, Influence of low-frequency vibration and modification on solidification and mechanical properties of al-si casting alloy, Materials 10 (2017) 560. URL: https://doi.org/10.3390/ma10050560. doi:10.3390/ma10050560. [11] Y. V. Dotsenko, V. Y. Selivorstov, T. Selivorstova, N. Dotsenko, Influence of heterogeneous crystallization conditions of aluminum alloy on its plastic properties., Scientific Bulletin of National Mining University (2015). [12] W. hua SUN, L. jun ZHANG, M. WEI, Y. DU, B. yun HUANG, Effect of liquid diffusion coefficients on microstructure evolution during solidification of al356.1 alloy, Transactions of Nonferrous Metals Society of China 23 (2013) 3722–3728. URL: https://doi.org/10.1016/ s1003-6326(13)62922-2. doi:10.1016/s1003-6326(13)62922-2. [13] F. Nedopekin, V. Borodin, Mathematical modeling of processes in continuously cast ingot, Russ. Metall. (1985) 62–68. [14] R. Tuttle, S. Ramrattan, L. Wells, As-cast surface characterization for steel using disk- shaped chemically bonded sand specimens, International Journal of Metalcasting (2020). URL: https://doi.org/10.1007/s40962-020-00520-4. doi:10.1007/s40962-020-00520-4. [15] E. L. Kuznetsova, I. A. Selin, V. F. Formalev, ANALYTICAL STUDY OF STEFAN-TYPE PROB- LEMS IN COMPOSITES WITH AN ARBITRARY NUMBER OF MOVING BOUNDARIES OF PHASE TRANSITIONS, Composites: Mechanics, Computations, Applications, An In- ternational Journal 1 (2010) 25–35. URL: https://doi.org/10.1615/compmechcomputapplintj. v1.i1.20. doi:10.1615/compmechcomputapplintj.v1.i1.20. [16] W. Wang, S. Luo, M. Zhu, Performance optimization and evaluation of a 3d CA-FVM model for dendritic growth of fe-c alloy, in: The Minerals, Metals & Materials Series, Springer In- ternational Publishing, 2016, pp. 101–108. URL: https://doi.org/10.1007/978-3-319-65133-0_ 13. doi:10.1007/978-3-319-65133-0_13. [17] C. Su, Research on the early temperature field characteristics of double-block track, Journal of Railway Engineering Society 36 (2019) 33–40. [18] T. Mikhaylovskaya, V. Selivorstov, Computer program sbheat - engineering calculation of the temperature field of a semi-bounded casting, 2007. Certificate of registration of copyright to the work β„– 22135. [19] I. Smyrnova, T. Selivyorstova, S. Liulchak, I. Sezonova, R. Yuriy, V. Liashenko, The results of simulation of the process of occurrence of damages to the semiconductor elements of radio-electronic equipment under the influence of multi-frequency signals of short duration, International Journal of Advanced Trends in Computer Science and Engineering 9 (2020) 3053–3056. URL: https://doi.org/10.30534/ijatcse/2020/86932020. doi:10.30534/ ijatcse/2020/86932020. [20] I. Smyrnova, V. Horbenko, A. Lutsyshyn, V. Kaminskyi, Z. Sasiuk, T. Selivyorstova, I. Ienina, The method of determining the probability of affection of the semiconductor elements under the influence of the multifrequency space-time signals, International Journal of Emerging Trends in Engineering Research 8 (2020) 1776–1779. URL: https://doi.org/10. 30534/ijeter/2020/46852020. doi:10.30534/ijeter/2020/46852020. [21] K. Venkateshwar, S. Tasnim, H. Simha, S. Mahmud, Effect of spatially varying morphologies of metal foams on phase change process, Thermal Science and Engineering Progress 19 (2020) 100667. URL: https://doi.org/10.1016/j.tsep.2020.100667. doi:10.1016/j.tsep. 2020.100667. [22] J. A. Lyon, A. J. Magana, A review of mathematical modeling in engineering education, The International journal of engineering education 36 (2020) 101–116. [23] N. I. Nikitenko, Y. F. Snezhkin, N. N. Sorokovaya, Development of a theory and methods for calculating the heat and mass transfer in drying a porous body with multicomponent vapor and liquid phases, Journal of Engineering Physics and Thermophysics 81 (2008) 1153–1167. URL: https://doi.org/10.1007/s10891-009-0132-x. doi:10.1007/s10891-009-0132-x. [24] L. S. Vasil’ev, Possible mechanisms for the nucleation of primary fracture zones dur- ing deformation-induced phase transformations in solids: IV. kinetic mechanisms of submicrocrack nucleation, Russian Metallurgy (Metally) 2008 (2008) 71–75. URL: https: //doi.org/10.1134/s003602950801014x. doi:10.1134/s003602950801014x. [25] M. Kawamoto, T. Nagira, Flowering of continuous casting process for steel in japan and new fundamental seeds to the future, Tetsu-to-Hagane 100 (2014) 472–484. URL: https://doi.org/10.2355/tetsutohagane.100.472. doi:10.2355/tetsutohagane.100.472. [26] T. Iida, R. Guthrie, Predictive models for the velocity of sound in liquid metallic elements at their melting-point temperatures, Metallurgical and Materials Transac- tions B 40 (2009) 949–958. URL: https://doi.org/10.1007/s11663-009-9261-z. doi:10.1007/ s11663-009-9261-z. [27] A. Romanov, Effect of diffusion on thermal processes of alloy crystallization. ii. application of heat flow method to analysis of experimental alloy cooling curves., Rasplavy (2004) 14–20. [28] D. You, S. K. Michelic, P. Presoly, J. Liu, C. Bernhard, Modeling inclusion formation during solidification of steel: A review, Metals 7 (2017) 460. URL: https://doi.org/10.3390/ met7110460. doi:10.3390/met7110460. [29] L. Boriwal, R. M. Sarviya, M. M. Mahapatra, Modeling the resistance spot welding of galvanized steel sheets using neuro-fuzzy method, in: Advances in Intelligent Sys- tems and Computing, Springer Singapore, 2017, pp. 37–50. URL: https://doi.org/10.1007/ 978-981-10-5272-9_4. doi:10.1007/978-981-10-5272-9_4. [30] R. von Hanxleden, E. A. Lee, C. Motika, H. Fuhrmann, Multi-view modeling and prag- matics in 2020, in: Large-Scale Complex IT Systems. Development, Operation and Man- agement, Springer Berlin Heidelberg, 2012, pp. 209–223. URL: https://doi.org/10.1007/ 978-3-642-34059-8_11. doi:10.1007/978-3-642-34059-8_11. [31] H. Nishihara, T. Kyotani, Zeolite-templated carbons – three-dimensional microporous graphene frameworks, Chemical Communications 54 (2018) 5648–5673. URL: https://doi. org/10.1039/c8cc01932k. doi:10.1039/c8cc01932k. [32] V. V. Mazalov, Method of cumulative sums and problems of optimal control by technological processes, International Journal of Software Engineering and Knowl- edge Engineering 07 (1997) 219–229. URL: https://doi.org/10.1142/s0218194097000126. doi:10.1142/s0218194097000126. [33] M. M. Sarsengeldin, S. Guvercin, S. Kassabek, Solution of one phase inverse stefan problem by IEF method, Journal of Physics: Conference Series 495 (2014) 012021. URL: https: //doi.org/10.1088/1742-6596/495/1/012021. doi:10.1088/1742-6596/495/1/012021. [34] D. S. Svyetlichnyy, J. Nowak, Ł. Łach, Modeling of recrystallization with recovery by frontal cellular automata, in: International Conference on Cellular Automata, Springer, 2012, pp. 494–503. [35] Ł. Łach, J. Nowak, D. Svyetlichnyy, The evolution of the microstructure in AISI 304l stain- less steel during the flat rolling – modeling by frontal cellular automata and verification, Journal of Materials Processing Technology 255 (2018) 488–499. URL: https://doi.org/10. 1016/j.jmatprotec.2017.12.001. doi:10.1016/j.jmatprotec.2017.12.001. [36] S. P. Gill, M. G. Cornforth, A. C. Cocks, Modelling microstructure evolution in engineering materials, International Journal of Plasticity 17 (2001) 669–690. URL: https://doi.org/10. 1016/s0749-6419(00)00066-8. doi:10.1016/s0749-6419(00)00066-8. [37] M. B. Cortie, Simulation of metal solidification using a cellular automaton, Metal- lurgical Transactions B 24 (1993) 1045–1053. URL: https://doi.org/10.1007/bf02660996. doi:10.1007/bf02660996. [38] A. N. Chichko, Y. V. Yatskevich, Three-dimensional computer modeling of the casting cooling and selection of the gating system runner on basis of heat conduction equation, Inzhenerno-Fizicheskii Zhurnal 72 (1999) 797–802. [39] D. Griffeath, C. Moore, New constructions in cellular automata, Oxford University Press on demand, 2003. [40] N. Konno, Quantum walks and quantum cellular automata, in: Lecture Notes in Com- puter Science, Springer Berlin Heidelberg, 2008, pp. 12–21. URL: https://doi.org/10.1007/ 978-3-540-79992-4_3. doi:10.1007/978-3-540-79992-4_3. [41] V. M. Golod, L. Y. Dobosh, Computation material science of structural-phase transformation in casting aluminium alloys, IOP Conference Series: Materials Science and Engineering 192 (2017) 012027. URL: https://doi.org/10.1088/1757-899x/192/1/012027. doi:10.1088/ 1757-899x/192/1/012027. [42] F. He, L. Zhou, Formation mechanism and influence factors of the sticker between solidified shell and mold in continuous casting of steel, High Temperature Materials and Processes 38 (2019) 192–198. URL: https://doi.org/10.1515/htmp-2018-0005. doi:10. 1515/htmp-2018-0005. [43] T. Toffoli, What automata can provide a medium for life?, in: Cellular Automata and Discrete Complex Systems, Springer International Publishing, 2016, pp. 10–25. URL: https: //doi.org/10.1007/978-3-319-39300-1_2. doi:10.1007/978-3-319-39300-1_2. [44] T. Mikhailovskaya, A. Mikhalev, Three-dimensional cellular automata model of the solidification process of metals and alloys, Modern problems of metallurgy. Scientific bulletin 9 (2006) 102–112.