Acessibilidade / Reportar erro

Application of self-organizing maps to evaluate the influence and behavior of the film formed during salting of Prato cheese

Abstract

Problems related to high blood pressure have led consumers to choose foods with low levels of sodium chloride, and an alternative to reduce the salt content of this salt is its partial replacement with KCl. In order to not affect the sensory properties of cheese, salting was performed using a solution containing 70% of NaCl and 30% of KCl in a static and dynamic system. During the salting process by immersion, the mass transfer is affected by the formation of a film on the cheese surface. The analysis of the diffusion of salts in the film can be performed using self-organizing map (SOM) combined with the Finite Element Method (FEM). Through these tools, a greater influence of the film was observed when the static system was applied and that the diffusion is different according to the studied position. On the biosolid sides the diffusion was more pronounced than in the center, indicating a decrease in the film thickness towards the edges. The SOM combined with 3D modeling showed to be an efficient tool to investigate the formation, influence and behavior of the film during the diffusion of Na+ and K+ on salting of Prato cheese.

Keywords:
neural networks; diffusion; finite element methods; simulation

1 Introduction

Currently, a direct relation between sodium content in food and its influence on blood pressure has led consumers to choose healthier products. Thus, the production of foods with low sodium content is necessary, but sodium chloride plays a fundamental role in the final quality of a product. Therefore, its total substitution is not recommended (Bordin et al., 2019Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
). KCl can replace NaCl by up to 30% without affecting the acceptability of the final product. It is possible to replace part of the sodium chloride with potassium chloride, as this salt does not alter the sensory, physicochemical or microbiological characteristics of the product. In addition, potassium intake increases kidney excretion of sodium, resulting in an antihypertensive effect and also exhibits antimicrobial activity (Bordin et al., 2019Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
; Borsato et al., 2012Borsato, D., Moreira, M. B., Moreira, I., Pina, M. V. R., Silva, R. S., & Bona, E. (2012). Saline distribution during multicomponent salting in pre-cooked quail eggs. Food Science and Technology, 32(2), 281-288. http://dx.doi.org/10.1590/S0101-20612012005000060.
http://dx.doi.org/10.1590/S0101-20612012...
).

Cheese is a product widely used by Brazilians. The most consumed types are: Prato, Mozzarella, Minas and Parmesan. Among them, Prato cheese stands out for its high nutritional value. A homogeneous distribution of NaCl is essential for its final quality (Bona et al., 2006Bona, E., Silva, L. H. M., Borsato, D., Silva, R. S. F., Fidelis, D. A. de S., & Araújo, A. (2006). Optimization of space and time discretization during the finite element method application to multicomponent diffusion simulation. Acta Scientiarum. Technology, 28, 141-150., 2010Bona, E., Silva, R. S., Borsato, D., Silva, L. H. M., & Fidelis, D. A. S. (2010). Multicomponent diffusion during Prato cheese ripening: mathematical modeling using the finite element method. Food Science and Technology, 30(4), 955-963. http://dx.doi.org/10.1590/S0101-20612010000400018.
http://dx.doi.org/10.1590/S0101-20612010...
; Silva et al., 2017Silva, H. L. A., Balthazar, C. F., Esmerino, E. A., Vieira, A. H., Cappato, L. P., Cucinelli, R. P. No., Verruck, S., Cavalcanti, R. N., Portela, J. B., Andrade, M. M., Moraes, J., Franco, R. M., Tavares, M. I. B., Prudencio, E. S., Freitas, M. Q., Nascimento, J. S., Silva, M. C., Raices, R. S. L., & Cruz, A. G. (2017). Effect of sodium reduction and flavor enhancer addition on probiotic Prato cheese processing. Food Research International, 99(Pt 1), 247-255. http://dx.doi.org/10.1016/j.foodres.2017.05.018. PMid:28784481.
http://dx.doi.org/10.1016/j.foodres.2017...
). Therefore, salting is an important step in cheesemaking, consisting of its immersion in brine, static or dynamic, in which the salt spreads into the solid by mass diffusion mechanisms (Albarracín et al., 2011Albarracín, W., Sánchez, I. C., Grau, R., & Barat, J. M. (2011). Salt in food processing; usage and reduction: a review. International Journal of Food Science & Technology, 46(7), 1329-1336. http://dx.doi.org/10.1111/j.1365-2621.2010.02492.x.
http://dx.doi.org/10.1111/j.1365-2621.20...
; Guinee, 2004Guinee, T. P. (2004). Salting and the role of salt in cheese. International Journal of Dairy Technology, 57(2‐3), 99-109. http://dx.doi.org/10.1111/j.1471-0307.2004.00145.x.
http://dx.doi.org/10.1111/j.1471-0307.20...
).

Studies on the physicochemical properties and characteristics of various cheeses have recently been conducted. Matera et al. (2018)Matera, J., Luna, A. S., Batista, D. B., Pimentel, T. C., Moraes, J., Kamimura, B. A., Ferreira, M. V. S., Silva, H. L. A., Mathias, S. P., Esmerino, E. A., Freitas, M. Q., Raices, R. S. L., Quitério, S. L., Sant’Ana, A. S., Silva, M. C., & Cruz, A. G. (2018). Brazilian cheeses: a survey covering physicochemical characteristics, mineral content, fatty acid profile and volatile compounds. Food Research International, 108, 18-26. http://dx.doi.org/10.1016/j.foodres.2018.03.014. PMid:29735047.
http://dx.doi.org/10.1016/j.foodres.2018...
studied the physicochemical characteristics of Brazilian cheeses Salum et al. (2019)Salum, P., Erbay, Z., & Selli, S. (2019). The compositional properties, proteolytic–lipolytic maturation parameters and volatile compositions of commercial enzyme‐modified cheeses with different cheese flavours. International Journal of Dairy Technology, 72(3), 416-426. http://dx.doi.org/10.1111/1471-0307.12591.
http://dx.doi.org/10.1111/1471-0307.1259...
the compositional properties of enzyme-modified commercial cheeses, Punoo et al. (2018)Punoo, H. A., Patil, G., & Bijoy, R. R. (2018). Textural and microstructural properties of Kradi cheese (an indigenous cheese of Jammu and Kashmir, India). International Journal of Dairy Technology, 71(2), 372-381. http://dx.doi.org/10.1111/1471-0307.12459.
http://dx.doi.org/10.1111/1471-0307.1245...
the textural and microstructural properties of Kradi cheese and Rafiq et al. (2018)Rafiq, S., Huma, N., Gulzar, N., Murtaza, M. A., & Hussain, I. (2018). Effect of cheddar cheese peptide extracts on growth inhibition, cell cycle arrest and apoptosis induction in human lung cancer (H‐1299) cell line. International Journal of Dairy Technology, 71(4), 975-980. http://dx.doi.org/10.1111/1471-0307.12533.
http://dx.doi.org/10.1111/1471-0307.1253...
cheddar cheese peptide extracts in anticancer activity. Researches related to Prato cheese include the partial replacement of NaCl by KCl and the addition of flavor enhancers to the probiotic Prato cheese (Silva et al., 2018Silva, H. L. A., Balthazar, C. F., Esmerino, E. A., Cucinelli, R. P. No., Rocha, R. S., Moraes, J., & Santos, J. S. (2018). Partial substitution of NaCl by KCl and addition of flavor enhancers on probiotic Prato cheese: a study covering manufacturing, ripening and storage time. Food Chemistry, 248, 192-200. http://dx.doi.org/10.1016/j.foodchem.2017.12.064.
http://dx.doi.org/10.1016/j.foodchem.201...
), the change in volatile compound concentration during the maturation of Prato cheese by reducing fat and adding whey protein (Domingos et al., 2019Domingos, L. D., Souza, H. A. L., Mariutti, L. R. B., Benassi, M. T., Bragagnolo, N., & Viotto, W. H. (2019). Fat reduction and whey protein concentrate addition alter the concentration of volatile compounds during Prato cheese ripening. Food Research International, 119, 793-804. http://dx.doi.org/10.1016/j.foodres.2018.10.062. PMid:30884718.
http://dx.doi.org/10.1016/j.foodres.2018...
), and attenuation of cigarette smoke-induced lesions in mice by the probiotic Prato cheese (Vasconcelos et al., 2019Vasconcelos, F. M., Silva, H. L. A., Poso, S. M. V., Barroso, M. V., Lanzetti, M., Rocha, R. S., Graça, J. S., Esmerino, E. A., Freitas, M. Q., Silva, M. C., Raices, R. S. L., Granato, D., Pimentel, T. C., Sant’Ana, A. S., Cruz, A. G., & Valença, S. S. (2019). Probiotic Prato cheese attenuates cigarette smoke-induced injuries in mice. Food Research International, 123, 697-703. http://dx.doi.org/10.1016/j.foodres.2019.06.001. PMid:31285019.
http://dx.doi.org/10.1016/j.foodres.2019...
).

According to Schwartzberg & Chao (1982)Schwartzberg, H. G., & Chao, R. Y. (1982). Solute diffusivities in leaching processes. Food Technology, 36(2), 73-86., when a fluid is in contact with a solid, a film is formed on its surface. If there is a mass transfer between the surface and the fluid, the current must pass through the stationary layer, which acts as a resistive barrier, and therefore there is a specific diffusion coefficient that describes the mass transfer in this film. The question is that this coefficient cannot be determined only experimental, and so to evaluate the influence of the formed film it is necessary to take into account the ratio between the resistance to internal and external mass diffusivity that can be quantified by the Biot mass number. According to Bona et al. (2007)Bona, E., Carneiro, R. L., Borsato, D., Silva, R. S., Fidelis, D. A. S., & Silva, L. H. M. (2007). Simulation of NaCl and KCl mass transfer during salting of Prato cheese in brine with agitation: a numerical solution. Brazilian Journal of Chemical Engineering, 24(3), 337-349. http://dx.doi.org/10.1590/S0104-66322007000300004.
http://dx.doi.org/10.1590/S0104-66322007...
, high values of the Biot number indicate that the internal resistance is limiting, and during the diffusion process as the Biot number decreases, the resistance of the film formed at the solid/solution interface increases (Rakotondramasy-Rabesiaka et al., 2010Rakotondramasy-Rabesiaka, L., Havet, J.-L., Porte, C., & Fauduet, H. (2010). Estimation of effective diffusion and transfer rate during the protopine extraction process from Fumaria officinalis L. Separation and Purification Technology, 76(2), 126-131. http://dx.doi.org/10.1016/j.seppur.2010.09.030.
http://dx.doi.org/10.1016/j.seppur.2010....
; Schwartzberg & Chao, 1982Schwartzberg, H. G., & Chao, R. Y. (1982). Solute diffusivities in leaching processes. Food Technology, 36(2), 73-86.). Thus, the Biot number is related to the film coefficient, since the higher its value, the lower is the influence of the formed film.

When it comes to food biosolids the diffusion on surface as well as film formation may be dependent on the geometry and position where the diffusion occurs. Thus, to better understand the effects of these factors during the process, some data analysis tools can be used through unconventional statistics (Cremasco et al., 2016Cremasco, H., Borsato, D., Angilelli, K. G., Galão, O. F., Bona, E., & Valle, M. E. (2016). Application of self‐organising maps towards segmentation of soybean samples by determination of inorganic compounds content. Journal of the Science of Food and Agriculture, 96(1), 306-310. http://dx.doi.org/10.1002/jsfa.7094. PMid:25641560.
http://dx.doi.org/10.1002/jsfa.7094...
, 2019Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917.
http://dx.doi.org/10.1590/fst.40917...
; Haykin, 2001Haykin, S. (2001). Neural networks: a comprehensive foundation. Englewood Cliffs: Preditice Hall.).

One tool is the SOM-type ANN. It has been applied to solve various types of problems in different areas of science such as engineering, chemistry, medicine, bioenergy, among others (Bona et al., 2012Bona, E., Silva, R. S. S. F., Borsato, D., & Bassoli, D. G. (2012). Self-organizing maps as a chemometric tool for aromatic pattern recognition of soluble coffee. Acta Scientiarum. Technology, 34(1), 111-119. http://dx.doi.org/10.4025/actascitechnol.v34i1.10892.
http://dx.doi.org/10.4025/actascitechnol...
; Borsato et al., 2011Borsato, D., Pina, M. V. R., Spacino, K. R., Scholz, M. B. S., & Androcioli, A. Fo (2011). Application of artificial neural networks in the geographical identification of coffee samples. European Food Research and Technology, 233(3), 533-543. http://dx.doi.org/10.1007/s00217-011-1548-z.
http://dx.doi.org/10.1007/s00217-011-154...
; Cremasco et al., 2016Cremasco, H., Borsato, D., Angilelli, K. G., Galão, O. F., Bona, E., & Valle, M. E. (2016). Application of self‐organising maps towards segmentation of soybean samples by determination of inorganic compounds content. Journal of the Science of Food and Agriculture, 96(1), 306-310. http://dx.doi.org/10.1002/jsfa.7094. PMid:25641560.
http://dx.doi.org/10.1002/jsfa.7094...
; Lemes & Dal Pino, 2008Lemes, M. R., & Dal Pino, A. Jr. (2008). Periodic table of the elements in the perspective of artificial neural networks. Química Nova, 31(5), 1141-1144. http://dx.doi.org/10.1590/S0100-40422008000500040.
http://dx.doi.org/10.1590/S0100-40422008...
; Lindsey et al., 2018Lindsey, R., Daluiski, A., Chopra, S., Lachapelle, A., Mozer, M., Sicular, S., Hanel, D., Gardner, M., Gupta, A., Hotchkiss, R., & Potter, H. (2018). Deep neural network improves fracture detection by clinicians. Proceedings of the National Academy of Sciences of the United States of America, 115(45), 11591-11596. http://dx.doi.org/10.1073/pnas.1806905115.
http://dx.doi.org/10.1073/pnas.180690511...
; Walkoff et al., 2017Walkoff, A. R., Antunes, S. R. M., Arrúa, M. E. P., Silva, L. R. C., Borsato, D., & Rodrigues, P. R. P. (2017). Self-organizing maps neural networks applied to the classification of ethanol samples according to the region of commercialization. Orbital: The Electronic Journal of Chemistry, 9(4), 248-255.). The purpose of this tool is to transform an input pattern of arbitrary dimension into a discrete one- or two-dimensional map, and perform this transformation adaptively in a topologically ordered manner. The interesting thing about it is that the result can be observed through a topological map whose function is to facilitate the interpretation of the results (Kohonen & Maps, 1995Kohonen, T., & Maps, S.-O. (1995). Self-organizing maps (Springer Series in Information Sciences, No. 30). Berlin: Springer. http://dx.doi.org/10.1007/978-3-642-97610-0.
http://dx.doi.org/10.1007/978-3-642-9761...
).

The objective of this work was to investigate the formation, influence and behavior of the film during multicomponent diffusion of NaCl and KCl, in static and dynamic brine, applying ANN type SOM with 3D computational modeling using the Finite Element Method (FEM).

2 Materials and methods

2.1 Prato cheese

It was used 3 kg of Prato cheese with rectangular geometry, provided by Laticinios Campina Alta (Manoel Ribas-PR). The cheese was divided into standardized geometry samples of 0.04 m × 0.04 m × 0.02 m (Figure 1).

Figure 1
Solid automatically generated by the software and the positions used in the simulation associated with a set of coordinates x, y, z.

2.2 Brine preparation and sampling

It was used 15 L of brine with 5% (w/v) salt concentration, containing 30% of potassium chloride (KCl) and 70% of sodium chloride (NaCl). For the system with agitation was used a pump with 520 L h-1 of circulation flow, with a constant temperature of 20 (± 1 °C). Prato cheese samples were completely immersed in the static and dynamic brine being collected, and their dimensions measured using a digital caliper at times of 0, 0.25, 0.75, 1.0, 2.5, 6.5, 7.5, 9.5, 10.5, 11.5, 17.5, 19.0, 29.0, 36.0, 40.0, 50.0 and 62.0 hours.

2.3 Determination of sodium and potassium

NaCl and KCl concentrations in Prato cheese samples were determined according to the methodology described by Bordin et al. (2019)Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
with modifications, using the Micronal photometer, model B-462, with air pressure of 0.8 kgf cm-2 and 1.5 kgf cm-2 air pump pressure using butane gas.

2.4 Three-dimensional modeling

The modeling was performed using the finite element method considering a rectangular domain, three-dimensional mass transfer, the generalized equations of Fick second law and Onsager (1945Onsager, L. (1945). Theories and problems of liquid diffusion. Annals of the New York Academy of Sciences, 46(5), 241-265. http://dx.doi.org/10.1111/j.1749-6632.1945.tb36170.x. PMid:21024247.
http://dx.doi.org/10.1111/j.1749-6632.19...
) equations. The diffusion coefficient was considered constant, regardless of the position and the immersion time in the brine. The solute diffusion occurred under isothermal conditions (20 °C) and the cheese contraction during salting was not considered. To evaluate the influence of the film formed on the cheese surface, in the two salting systems used, the Cauchy boundary condition was considered (Bona et al., 2007Bona, E., Carneiro, R. L., Borsato, D., Silva, R. S., Fidelis, D. A. S., & Silva, L. H. M. (2007). Simulation of NaCl and KCl mass transfer during salting of Prato cheese in brine with agitation: a numerical solution. Brazilian Journal of Chemical Engineering, 24(3), 337-349. http://dx.doi.org/10.1590/S0104-66322007000300004.
http://dx.doi.org/10.1590/S0104-66322007...
). The finite element formulation followed the procedures was established by Cremasco et al. (2019)Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917.
http://dx.doi.org/10.1590/fst.40917...
.

2.5 Finite element method simulation

The simulation was performed using COMSOL Multiphysics® software and the standard “Transport of Diluted Species (tds)” physical interface. The parameters used in the simulation were the main and cross diffusion coefficients, Biot mass number and the relationship between the mass transfer coefficient and the mass conductivity (hmm), that were adjusted using the super modified simplex method (Bona et al., 2000Bona, E., Borsato, D., Silva, R. S. S. F., & Herrera, R. P. (2000). Software for optimization using the sequential simplex method. Acta Scientiarum, 22, 1201-1206.). Figure 1 shows the solid generated automatically by the software, showing part of the extremely fine mesh used and the equidistant points chosen for the study of Na+ and K+ diffusion on the surface of the film formed in the Prato cheese.

2.6 Artificial Neural Networks (ANN)

For the analysis of the influence of the film on the Prato cheese surface, a Kohonen self-organizing map (SOM) was used. The applied network presented a 6x6 hexagonal topology with 8000 training epochs, with learning rate training starting at 0.2 and decaying exponentially with the training epochs to 6.71x10-5, and with the initial neighborhood relation of 3.5 decaying to 0.045.

2.7 Computational processing

For the simulation it was used the software COMSOL Multiphysics® version 5.2 (COMSOL, Inc., Burlington, MA) based on the finite element method. The neural network routine developed by our research group was used according to the algorithm described in Haykin (2001)Haykin, S. (2001). Neural networks: a comprehensive foundation. Englewood Cliffs: Preditice Hall. and processed by Matlab® R2007b software to evaluate the influence and behavior of NaCl and KCl concentrations in the Prato cheese surface.

3 Results and discussion

The main and cross diffusion coefficients, the film coefficient, the ratio between the mass transfer coefficient and the mass conductivity, and the Biot number of the Prato cheese, subjected to salting in the static and dynamic process, were determined by simulation using the finite element method with the application of a simplex optimization coupled with the desirability functions (Bordin et al., 2019Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
). All simulations were performed using a 3D geometry modeling, with time-dependent solver Backward Differentiation Formula (BDF) with extremely fine tetrahedral finite element mesh, composed of 130,382 elements with 370,052 degrees of freedom. Since the method is iterative and constrained, to avoid possible oscillations and an excessive number of vertices in the simplex procedure (Bona et al., 2000Bona, E., Borsato, D., Silva, R. S. S. F., & Herrera, R. P. (2000). Software for optimization using the sequential simplex method. Acta Scientiarum, 22, 1201-1206.; Bordin et al., 2019Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
), the lower and upper limits of each variable were defined through previous tests. The results obtained were within the range of the independent variables chosen and convergence to the optimal was achieved faster (Angilelli et al., 2015Angilelli, K. G., Orives, J. R., Silva, H. C., Coppo, R. L., Moreira, I., & Borsato, D. (2015). Multicomponent diffusion during osmotic dehydration process in melon pieces: influence of film coefficient. Journal of Food Processing and Preservation, 39(4), 329-337. http://dx.doi.org/10.1111/jfpp.12236.
http://dx.doi.org/10.1111/jfpp.12236...
).

The values of the main (Dii) and cross (Dij) diffusion coefficients, the film coefficient (hm), the ratio between the mass transfer and the mass conductivity coefficients (hmm), the percentage errors obtained from the difference between experimental and simulated data, and the Biot number, using the COMSOL Multiphysics® software, during the diffusion process, without and with agitation, are presented in Table 1.

Table 1
Parameters obtained by simplex optimization during the sodium and potassium chloride diffusion processes in static and dynamic salting.

The optimized values of Biot numbers and film coefficients (Table 1) show a difference in mass transfer between the static and the dynamic process. The lower value for the Biot number and the higher value for hm/ λm show that the influence of the barrier on the solution/biosolid interface is greater in the static system, indicating ions resistance to transfer from the solution to the Prato cheese (Bordin et al., 2019Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188. PMid:30292382.
http://dx.doi.org/10.1016/j.foodchem.201...
).

According to Table 1, it is observed that the main and cross coefficients in the dynamic and the static brine are the same, since these parameters are related to the mass transfer inside the biosolid, being not dependent on the external disturbance of the system. The main diffusion coefficient of Na+ (D11) is 1.7 times higher than the main diffusion coefficient of K+ (D22). Also, the cross coefficients presented smaller values than the main ones, showing that the diffusion of the solutes in their own flow is more important than the interference between them.

According to Cremasco et al. (2019)Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917.
http://dx.doi.org/10.1590/fst.40917...
, in a study of the diffusion in the solution/quail egg interface, the influence of the film formed during the diffusion process can be studied by analyzing the solute concentrations at different positions over time. The diffusion profile in static and dynamic systems, comparing experimental and simulated concentrations in g salt/100g (solution) is shown in Figure 2. Therefore, there must be a physical barrier on the outer surface of Prato cheese as it was observed that the mass transfer was faster when the stirring system was applied (Figure 2b). In addition, the Biot number for the dynamic system was 84 and for the static system 20.52 (Table 1), indicating the greatest influence of the film formed on the cheese surface when the salting process is not agitated (Figure 2a).

Figure 2
Diffusion profile of experimental and simulated salt concentrations during the static brine (a) and dynamic brine (b) experiments. NaCl experimental data are represented by (▲) and KCl by (■), NaCl and KCl simulated data are represented by (-) and (…).

The NaCl and KCl concentrations used in brine were 597.83 and 201.48 mol.m-3, respectively. These should be the concentrations obtained on the surface of the cheese immediately after immersion if the Dirichlet boundary conditions were applied (Chung, 1978Chung, T. J. (1978). Finite element analysis in fluid dynamics (NASA STI/Recon Technical Report A, No. 78). New York: McGraw Hill.). However, after 15 minutes of salting, the simulated average values found in the static system were 322.76 and 100.40 mol m-3 for NaCl and KCl, respectively, and in the dynamic brine were 507.07 and 162.30 mol m- 3 for NaCl and KCl, respectively. These values prove the influence of the film formed on the surface and that it was not eliminated even with agitation, selecting the Cauchy conditions (Bona et al., 2007Bona, E., Carneiro, R. L., Borsato, D., Silva, R. S., Fidelis, D. A. S., & Silva, L. H. M. (2007). Simulation of NaCl and KCl mass transfer during salting of Prato cheese in brine with agitation: a numerical solution. Brazilian Journal of Chemical Engineering, 24(3), 337-349. http://dx.doi.org/10.1590/S0104-66322007000300004.
http://dx.doi.org/10.1590/S0104-66322007...
). Therefore, as salting time increases, these values tend to approach the initial concentration of the brine used.

To analyze the influence of the film formed on the surface of Prato cheese the concentration at each chosen point during 45 hours was simulated through FEM. The values obtained were analyzed through the application of self-organizing map (SOM)-type neural networks. In a preliminary study, 8000 epochs were used, where it was possible to verify that the stabilization of the mean quantization error occurred after 6000 epochs, so this value was used in the network training.

Figure 3 shows the topological maps with concentration distribution at the points located in the upper and side surface of the cheese (Figure 1), during the Na+ and K+ multicomponent diffusion process. Figures 33b represent the maps of the upper points on the x-axis and Figures 33d the side points on the z-axis for both dynamic (D) and static systems (S). According to the Figure 3 it is noted that the influence of the formed film is not the same for Na+ and K+, since the topological map presented changes when analyzed in the same axis. Thus, is possible to say that not only the intrinsic properties of the two ions, such as their ionic radius, charge density, among others, that determine their mobility on the cheese surface, but also the position and type of ion considered.

Figure 3
Distribution of the analyzed points according to the neuron for sodium (a) and potassium (b) on the x-axis and sodium (c) and potassium (d) on the z-axis, for the static (S) and dynamic (D) systems.

Figure 3a shows some positions that are different from the topological map shown in Figure 3b, i.e. the neuron where the 1D and 13D positions are located, for the K+ on the x-axis has similarity to the neurons that are located in the positions 3D, 4D, 11D, 2D and 12D. However, this behavior is not repeated in the topological map shown in Figure 3a. On the other hand, the same similarity relation between neurons observed in Figure 3c occurs in Figure 3d. In both axes the neural network used showed that on the surface of the cheese there is a clear separation between the dynamic and static processes, because in the topological map, the more distant the points the lower the similarity between them (Cremasco et al., 2019Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917.
http://dx.doi.org/10.1590/fst.40917...
).

As the analyzed positions move away from the center point (position 7 of Figure 1) towards the cheese edge, topological maps show that their similarity decreases, indicating that the thickness of the formed film varies along the x-axis. Thus, at position 7 the resistance to ion transfer is greater, decreasing as it moves away from the center of the cheese surface. The same behavior is observed when analyzing the positions located on the z-axis.

Figures 4 and 5 show the weight maps at the different times evaluated in the simulation (FEM) for Na+ (Figure 4) and K+ (Figure 5) in the static and dynamic systems, for x and z-axis. In these maps, the positions in the darker regions have greater weight, and as the color intensity decreases a lower weight is attributed.

Figure 4
Weight map of NaCl diffusion on the biosolid surface for the x-axis and z-axis in static (S) and dynamic (D) systems.
Figure 5
Weight maps of KCl diffusion on the biosolid surface for the x-axis and z-axis in the static (S) and dynamic (D) systems.

It is possible to observe that there are differences in ion concentration values in relation to the diffusion time in different positions of the biosolid. Note that there was a separation of the data obtained between the static and the dynamic system, in the positions studied during the simulation of the diffusion in Prato cheese. Each position defined for the static system presented lower concentration values when compared with the concentrations of the same positions in the dynamic system. This observation indicates that the influence of the film is more pronounced in the static system since the concentrations of species are lower than those obtained for the dynamic, because in this system the brine agitation reduces the film thickness, and consequently the mass transfer is more effective (Cremasco et al., 2019Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917.
http://dx.doi.org/10.1590/fst.40917...
).

Among the studied positions, both in the x and z-axes, as well as in the systems and ions, it was observed that the positions 1, 13 and 23 presented higher concentrations. The biosolid is a rectangular block and the diffusion happens faster in the positions located at the edges of the cheese. Therefore, it is suggested that the film formed is thinner at the edge and thicker at the central external part. In these positions, the time corresponding to 0.25h of salting is sufficient to reach the brine concentration in the dynamic system and only after 28h of salting the static system reaches this equilibrium.

In Figure 4, on average, a difference in concentration is observed as diffusion time increases. The Na+ concentration in the x-axis reached 85.71% and 50.00% of brine concentration after 0.25h of salting in the dynamic and static system, respectively. In the z-axis it reached 85.71% and 64.30% in the dynamic and static system, respectively. The Na+ reached 97.80% of brine concentration in both axes after 28h and 97.14% after 6h of brine in the static and dynamic systems, respectively.

According to Figure 5, in the x-axis, 96.00% of the K+, reached the brine concentration after 28 and 15 hours of salting in the static and dynamic system, respectively. On the z-axis, 97.00% of the K+ ions reached brine concentration after 24 hours in the static system and 18 hours in the dynamic system.

4 Conclusions

The optimization of the main, cross and hmm coefficients was performed by associating the simplex optimization with the finite element method (FEM). The comparison between the static and dynamic systems showed that the diffusion process is influenced by the film formed at the biosolid/solution interface, being minimized with the brine agitation.

The combination of SOM-type artificial neural networks with FEM simulation proved to be an appropriate tool, as a new application of Prato cheese salting, to evaluate the formation of the film, showing that this phenomenon influences the mass transfer and diffusion time.

  • Practical Application: Simulation of situations that will be very useful for industries working with diffusive processes.
  • Erratum

    Talita Fogaça de OLIVEIRA1, Marco Aurélio Jeanegitz CLEMENTE1, Diego GALVAN1, Gustavo FIX1, Ana Carolina Gomes MANTOVANI2, Dionísio BORSATO1* https://orcid.org/0000-0003-2281-0242
    1Departamento de Química, Universidade Estadual de Londrina – UEL, Londrina, PR, Brasil
    2Departamento de Física, Universidade Estadual de Londrina – UEL, Londrina, PR, Brasil
    Due to authors honest mistake the article “Application of self-organizing maps to evaluate the influence and behavior of the film formed during salting of Prato cheese” (DOI https://doi.org/10.1590/fst.33219), published in Food Science and Technology, Campinas, 40(Suppl. 2): 482-488, Dec 2020, was published without the acknowledgements.
    On page 487, after the Conclusions, it should be read:
    Acknowledgements
    National Council for Scientific and Technological Development – CNPq.
    The authors apologize for the mistake.

References

  • Albarracín, W., Sánchez, I. C., Grau, R., & Barat, J. M. (2011). Salt in food processing; usage and reduction: a review. International Journal of Food Science & Technology, 46(7), 1329-1336. http://dx.doi.org/10.1111/j.1365-2621.2010.02492.x
    » http://dx.doi.org/10.1111/j.1365-2621.2010.02492.x
  • Angilelli, K. G., Orives, J. R., Silva, H. C., Coppo, R. L., Moreira, I., & Borsato, D. (2015). Multicomponent diffusion during osmotic dehydration process in melon pieces: influence of film coefficient. Journal of Food Processing and Preservation, 39(4), 329-337. http://dx.doi.org/10.1111/jfpp.12236
    » http://dx.doi.org/10.1111/jfpp.12236
  • Bona, E., Borsato, D., Silva, R. S. S. F., & Herrera, R. P. (2000). Software for optimization using the sequential simplex method. Acta Scientiarum, 22, 1201-1206.
  • Bona, E., Carneiro, R. L., Borsato, D., Silva, R. S., Fidelis, D. A. S., & Silva, L. H. M. (2007). Simulation of NaCl and KCl mass transfer during salting of Prato cheese in brine with agitation: a numerical solution. Brazilian Journal of Chemical Engineering, 24(3), 337-349. http://dx.doi.org/10.1590/S0104-66322007000300004
    » http://dx.doi.org/10.1590/S0104-66322007000300004
  • Bona, E., Silva, L. H. M., Borsato, D., Silva, R. S. F., Fidelis, D. A. de S., & Araújo, A. (2006). Optimization of space and time discretization during the finite element method application to multicomponent diffusion simulation. Acta Scientiarum. Technology, 28, 141-150.
  • Bona, E., Silva, R. S. S. F., Borsato, D., & Bassoli, D. G. (2012). Self-organizing maps as a chemometric tool for aromatic pattern recognition of soluble coffee. Acta Scientiarum. Technology, 34(1), 111-119. http://dx.doi.org/10.4025/actascitechnol.v34i1.10892
    » http://dx.doi.org/10.4025/actascitechnol.v34i1.10892
  • Bona, E., Silva, R. S., Borsato, D., Silva, L. H. M., & Fidelis, D. A. S. (2010). Multicomponent diffusion during Prato cheese ripening: mathematical modeling using the finite element method. Food Science and Technology, 30(4), 955-963. http://dx.doi.org/10.1590/S0101-20612010000400018
    » http://dx.doi.org/10.1590/S0101-20612010000400018
  • Bordin, M. S. P., Borsato, D., Cremasco, H., Galvan, D., Silva, L. R. C., Romagnoli, É. S., & Angilelli, K. G. (2019). Mathematical modeling of multicomponent NaCl and KCl diffusion process during the salting of pre-cooked champignon mushrooms. Food Chemistry, 273, 99-105. http://dx.doi.org/10.1016/j.foodchem.2018.01.188 PMid:30292382.
    » http://dx.doi.org/10.1016/j.foodchem.2018.01.188
  • Borsato, D., Moreira, M. B., Moreira, I., Pina, M. V. R., Silva, R. S., & Bona, E. (2012). Saline distribution during multicomponent salting in pre-cooked quail eggs. Food Science and Technology, 32(2), 281-288. http://dx.doi.org/10.1590/S0101-20612012005000060
    » http://dx.doi.org/10.1590/S0101-20612012005000060
  • Borsato, D., Pina, M. V. R., Spacino, K. R., Scholz, M. B. S., & Androcioli, A. Fo (2011). Application of artificial neural networks in the geographical identification of coffee samples. European Food Research and Technology, 233(3), 533-543. http://dx.doi.org/10.1007/s00217-011-1548-z
    » http://dx.doi.org/10.1007/s00217-011-1548-z
  • Chung, T. J. (1978). Finite element analysis in fluid dynamics (NASA STI/Recon Technical Report A, No. 78). New York: McGraw Hill.
  • Cremasco, H., Borsato, D., Angilelli, K. G., Galão, O. F., Bona, E., & Valle, M. E. (2016). Application of self‐organising maps towards segmentation of soybean samples by determination of inorganic compounds content. Journal of the Science of Food and Agriculture, 96(1), 306-310. http://dx.doi.org/10.1002/jsfa.7094 PMid:25641560.
    » http://dx.doi.org/10.1002/jsfa.7094
  • Cremasco, H., Galvan, D., Angilelli, K. G., Borsato, D., & Oliveira, A. G. (2019). Influence of film coefficient during multicomponent diffusion–KCl/NaCl in biosolid for static and agitated system using 3D computational simulation. Food Science and Technology, 39(Suppl. 1), 173-181. http://dx.doi.org/10.1590/fst.40917
    » http://dx.doi.org/10.1590/fst.40917
  • Domingos, L. D., Souza, H. A. L., Mariutti, L. R. B., Benassi, M. T., Bragagnolo, N., & Viotto, W. H. (2019). Fat reduction and whey protein concentrate addition alter the concentration of volatile compounds during Prato cheese ripening. Food Research International, 119, 793-804. http://dx.doi.org/10.1016/j.foodres.2018.10.062 PMid:30884718.
    » http://dx.doi.org/10.1016/j.foodres.2018.10.062
  • Guinee, T. P. (2004). Salting and the role of salt in cheese. International Journal of Dairy Technology, 57(2‐3), 99-109. http://dx.doi.org/10.1111/j.1471-0307.2004.00145.x
    » http://dx.doi.org/10.1111/j.1471-0307.2004.00145.x
  • Haykin, S. (2001). Neural networks: a comprehensive foundation Englewood Cliffs: Preditice Hall.
  • Kohonen, T., & Maps, S.-O. (1995). Self-organizing maps (Springer Series in Information Sciences, No. 30). Berlin: Springer. http://dx.doi.org/10.1007/978-3-642-97610-0
    » http://dx.doi.org/10.1007/978-3-642-97610-0
  • Lemes, M. R., & Dal Pino, A. Jr. (2008). Periodic table of the elements in the perspective of artificial neural networks. Química Nova, 31(5), 1141-1144. http://dx.doi.org/10.1590/S0100-40422008000500040
    » http://dx.doi.org/10.1590/S0100-40422008000500040
  • Lindsey, R., Daluiski, A., Chopra, S., Lachapelle, A., Mozer, M., Sicular, S., Hanel, D., Gardner, M., Gupta, A., Hotchkiss, R., & Potter, H. (2018). Deep neural network improves fracture detection by clinicians. Proceedings of the National Academy of Sciences of the United States of America, 115(45), 11591-11596. http://dx.doi.org/10.1073/pnas.1806905115
    » http://dx.doi.org/10.1073/pnas.1806905115
  • Matera, J., Luna, A. S., Batista, D. B., Pimentel, T. C., Moraes, J., Kamimura, B. A., Ferreira, M. V. S., Silva, H. L. A., Mathias, S. P., Esmerino, E. A., Freitas, M. Q., Raices, R. S. L., Quitério, S. L., Sant’Ana, A. S., Silva, M. C., & Cruz, A. G. (2018). Brazilian cheeses: a survey covering physicochemical characteristics, mineral content, fatty acid profile and volatile compounds. Food Research International, 108, 18-26. http://dx.doi.org/10.1016/j.foodres.2018.03.014 PMid:29735047.
    » http://dx.doi.org/10.1016/j.foodres.2018.03.014
  • Onsager, L. (1945). Theories and problems of liquid diffusion. Annals of the New York Academy of Sciences, 46(5), 241-265. http://dx.doi.org/10.1111/j.1749-6632.1945.tb36170.x PMid:21024247.
    » http://dx.doi.org/10.1111/j.1749-6632.1945.tb36170.x
  • Punoo, H. A., Patil, G., & Bijoy, R. R. (2018). Textural and microstructural properties of Kradi cheese (an indigenous cheese of Jammu and Kashmir, India). International Journal of Dairy Technology, 71(2), 372-381. http://dx.doi.org/10.1111/1471-0307.12459
    » http://dx.doi.org/10.1111/1471-0307.12459
  • Rafiq, S., Huma, N., Gulzar, N., Murtaza, M. A., & Hussain, I. (2018). Effect of cheddar cheese peptide extracts on growth inhibition, cell cycle arrest and apoptosis induction in human lung cancer (H‐1299) cell line. International Journal of Dairy Technology, 71(4), 975-980. http://dx.doi.org/10.1111/1471-0307.12533
    » http://dx.doi.org/10.1111/1471-0307.12533
  • Rakotondramasy-Rabesiaka, L., Havet, J.-L., Porte, C., & Fauduet, H. (2010). Estimation of effective diffusion and transfer rate during the protopine extraction process from Fumaria officinalis L. Separation and Purification Technology, 76(2), 126-131. http://dx.doi.org/10.1016/j.seppur.2010.09.030
    » http://dx.doi.org/10.1016/j.seppur.2010.09.030
  • Salum, P., Erbay, Z., & Selli, S. (2019). The compositional properties, proteolytic–lipolytic maturation parameters and volatile compositions of commercial enzyme‐modified cheeses with different cheese flavours. International Journal of Dairy Technology, 72(3), 416-426. http://dx.doi.org/10.1111/1471-0307.12591
    » http://dx.doi.org/10.1111/1471-0307.12591
  • Schwartzberg, H. G., & Chao, R. Y. (1982). Solute diffusivities in leaching processes. Food Technology, 36(2), 73-86.
  • Silva, H. L. A., Balthazar, C. F., Esmerino, E. A., Vieira, A. H., Cappato, L. P., Cucinelli, R. P. No., Verruck, S., Cavalcanti, R. N., Portela, J. B., Andrade, M. M., Moraes, J., Franco, R. M., Tavares, M. I. B., Prudencio, E. S., Freitas, M. Q., Nascimento, J. S., Silva, M. C., Raices, R. S. L., & Cruz, A. G. (2017). Effect of sodium reduction and flavor enhancer addition on probiotic Prato cheese processing. Food Research International, 99(Pt 1), 247-255. http://dx.doi.org/10.1016/j.foodres.2017.05.018 PMid:28784481.
    » http://dx.doi.org/10.1016/j.foodres.2017.05.018
  • Silva, H. L. A., Balthazar, C. F., Esmerino, E. A., Cucinelli, R. P. No., Rocha, R. S., Moraes, J., & Santos, J. S. (2018). Partial substitution of NaCl by KCl and addition of flavor enhancers on probiotic Prato cheese: a study covering manufacturing, ripening and storage time. Food Chemistry, 248, 192-200. http://dx.doi.org/10.1016/j.foodchem.2017.12.064
    » http://dx.doi.org/10.1016/j.foodchem.2017.12.064
  • Vasconcelos, F. M., Silva, H. L. A., Poso, S. M. V., Barroso, M. V., Lanzetti, M., Rocha, R. S., Graça, J. S., Esmerino, E. A., Freitas, M. Q., Silva, M. C., Raices, R. S. L., Granato, D., Pimentel, T. C., Sant’Ana, A. S., Cruz, A. G., & Valença, S. S. (2019). Probiotic Prato cheese attenuates cigarette smoke-induced injuries in mice. Food Research International, 123, 697-703. http://dx.doi.org/10.1016/j.foodres.2019.06.001 PMid:31285019.
    » http://dx.doi.org/10.1016/j.foodres.2019.06.001
  • Walkoff, A. R., Antunes, S. R. M., Arrúa, M. E. P., Silva, L. R. C., Borsato, D., & Rodrigues, P. R. P. (2017). Self-organizing maps neural networks applied to the classification of ethanol samples according to the region of commercialization. Orbital: The Electronic Journal of Chemistry, 9(4), 248-255.

Publication Dates

  • Publication in this collection
    01 June 2020
  • Date of issue
    Dec 2020

History

  • Received
    19 Nov 2019
  • Accepted
    27 Jan 2020
Sociedade Brasileira de Ciência e Tecnologia de Alimentos Av. Brasil, 2880, Caixa Postal 271, 13001-970 Campinas SP - Brazil, Tel.: +55 19 3241.5793, Tel./Fax.: +55 19 3241.0527 - Campinas - SP - Brazil
E-mail: revista@sbcta.org.br