## Services on Demand

## Article

## Indicators

## Related links

## Share

## Brazilian Journal of Chemical Engineering

##
*Print version* ISSN 0104-6632

*On-line version* ISSN 1678-4383

### Braz. J. Chem. Eng. vol.24 no.3 São Paulo July/Sept. 2007

#### http://dx.doi.org/10.1590/S0104-66322007000300004

**FLUID DYNAMICS; HEAT AND MASS TRANSFER; AND OTHER TOPICS**

**Simulation of NaCl and KCl mass transfer during salting of Prato cheese in brine with agitation: a numerical solution**

**E. Bona ^{I,} ^{*}; R. L. Carneiro^{I}; D. Borsato^{I}; R. S. S. F. Silva^{II}; D. A. S. Fidelis^{II}; L. H. Monken e Silva^{III}**

^{I}Department of Chemistry, State University of Londrina, P.O. Box. 6001, CEP 86051-990, Londrina - PR, Brazil. E-mail: bona@uel.br

^{II}Department of Food Science and Technology, State University of Londrina, P.O. Box. 6001, CEP 86051-990, Londrina - PR, Brazil

^{III}School of Apucarana, Osvaldo de Oliveira St. 600, Apucarana - PR, Brazil

**ABSTRACT**

The association of dietary NaCl with arterial hypertension has led to a reduction in the levels of this salt in cheeses. For salting, KCl has been used as a partial substitute for NaCl, which cannot be completely substituted without affecting product acceptability. In this study a sensorially adequate saline solution (NaCl/KCl) was simultaneously diffused during salting of Prato cheese in brine with agitation. The simultaneous multicomponent diffusion during the process was modeled with Ficks second generalized law. The system of partial differential equations formed was solved by the finite element method (FEM). In the experimental data concentration the deviation for NaCl was of 7.3% and for KCl of 5.4%, both of which were considered acceptable. The simulation of salt diffusion will allow control and modulation of salt content in Prato cheese, permitting the prediction of final content from initial conditions.

**Keywords:** Multicomponent diffusion; Finite element method; Biot mass-exchange number.

**INTRODUCTION**

Sodium chloride is an important component in cheese production because it gives flavor, produces physical alterations in cheese proteins that influence texture, and also helps to control microbial growth and the activity of several enzymes in the cheese (Guinne and Fox, 1984). However, a portion of the population avoids cheese consumption to reduce salt intake due to problems related to arterial hypertension (Jameson, 1986; He and MacGregor, 2001).

An alternative proposed by Rapacci (1989) is to substitute some of the NaCl by KCl in the cheese formulation. According to the same author, the substitution of up to 30% of the NaCl by KCl does not produce significant differences in sensorial, physical, and chemical properties (p> 0.05) from traditional cheese salted with sodium chloride.

Salt transport in cheese can be simply described by Ficks second law, but many cheese and brine properties can affect the process. When several solutes diffuse simultaneously, Ficks second generalized law can be assumed (Onsager, 1945; Cussler, 1976), and in this case, in addition to the self-diffusion coefficient, cross diffusion coefficients are included to take into account the influence of a solute on the flow of the other one. Generally, for highly dissociable solutes the cross terms are smaller than the main ones (Gerla and Rubiolo, 2003). It seems to be well established that an appropriate salt concentration and its uniform distribution are important factors in cheese ripening (Guinne and Fox, 1984). Consequently, models that estimate saline concentrations as a result of processing could be useful when more than one component is considered in the salting.

Nowadays with the development of high-performance computers, it is possible to simulate a process close to reality using numerical techniques such as the finite element method, or FEM (Puri and Anantheswaran, 1993; Wang and Sun, 2003). Mass and heat transfer and solid and fluid mechanics are among the typical applications of this method (Bickford, 1990). According to Puri and Anantheswaran (1993), the main advantages of FEM are as follow: the spatial variation in material properties can be easily manipulated, irregular regions can be modeled with great precision, it is the best for nonlinear problems, the dimensions of the elements can be easily altered, the spatial interpretation is very realistic, and problems with the most diverse boundary conditions can be easily studied. However, the reliability of the solution obtained using FEM is directly related to the discretization of the model (Lyra, 1993). Silva et al. (1999) showed that due to the order of magnitude of the diffusion coefficient in cheese salting, the oscillation in the results could be dampened (Britz et al., 2002), optimizing the spatial and temporal discretization (Bona, 2004). This method has been applied successfully in simulating NaCl diffusion in cheese using Ficks second law as the theoretical support for this phenomenon (Silva et al., 1998, 1999) and in simulating simultaneous NaCl and KCl diffusion in static brine (Bona, 2004).

The objective of this study was to present an analysis of multicomponent mass transfer during salting with agitation for the NaCl-KCl-water system in Prato cheese using FEM.

**MATERIAL AND METHODS**

**Salting**

A piece of cheese measuring 24cm × 12cm × 9cm (Queijo Prato Lanche DI CARLO from Campina Alta dairy plant, located in Manoel Ribas – PR, Brazil) was used for the study. After pressing, the cheese was taken to the laboratory in the same mold where the curd had been pressed. The cheese was stored under refrigeration (10±1ºC) for approximately 12 hours before the experimental procedure. After this period it was removed from the mold and cut into smaller pieces (4cm × 9cm × 12cm) for the salting process. One of these samples was used to determine the fat, moisture, and initial sodium and potassium concentrations.

The cheeses were kept in stirred brine at 10±1ºC using an adapted domestic refrigerator (CONTINENTAL 560L). A pump (Motor Pump SARLO Better S520) with a 520L/h circulation flow was used for agitation purposes. Twenty liters of brine were prepared in accordance with Furtado (1990), who suggested a saline concentration of about 20% (w/w) and an average pH of 5.5 adjusted with HCl. To prevent Ca^{2+} loss from the cheeses and the consequent loss of rigidity, 0.5% (w/w) of CaCl_{2} was added. The brine volume used was approximately five times greater than the cheese volume to ensure that the brine concentration remained constant throughout the salting period (Gerla and Rubiolo, 2003). Before salting, a sample of brine was removed to quantify the salt contents, and another sample was removed at the end of the salting period. For studying multicomponent diffusion, approximately 70% NaCl and 30% KCl was used, as suggested by Rapacci (1989).

Because the cheese density was lower than the brine density, the cheese was immobilized by a network of cotton thread within the salting recipient (34cm × 49cm × 15cm), (Figure 1). This network guaranteed that none of the samples touched and all the surfaces were in contact with the brine, permitting uniform salting.

**Sampling**

The sampling procedure consisted of removing two cylindrical samples, each 1.5cm in diameter and 4cm in height, located as illustrated in Figure 2.

After removal, the cylinders were cut half to produce four cylindrical samples, as shown in Figure 3.

The cylinders belonging to the same face were analyzed together. The samples were removed for analysis after 1, 2, 3, 4, 5, and 7 hours of salting.

The dimensions shown in Figure 2, which were 4.5cm for the Y axis and 4.0cm for the Z axis, were sufficient to guarantee a one-dimensional approach to the problem (Gerla and Rubiolo, 2003; Luna and Bressan, 1986) because maximal NaCl penetration was estimated to be about 1.5cm for the salting time considered (Silva et al., 1998).

Two samples were removed to determine whether there was variability in the salt concentration during the extraction and quantification processes.

**Moisture and Fat**

Moisture was determined by drying in an oven at 105°C to constant weight (Richardson, 1985; AOAC, 1984) and fat, by the Soxhlet method (Cecchi, 1999).

**Sodium and Potassium**

Sodium and potassium content was determined by atomic emission (flame photometry), using a Celm FC-280 flame photometer, after incinerating the samples and extracting the salts with 0.5 mol/L HCl.

**MATHEMATICAL MODEL**

**Finite Element Formulation**

This study considered the simultaneous three-dimensional mass transfer of two solutes in a cheese that occupied a volume of W Ì R^{3}, where º [-R_{1}, R_{1}] × [-R_{2}, R_{2}] × [-R_{3}, R_{3}], associated with a system of Cartesian coordinates x, y, z with the origin located at the geometric center of the cheese. It was assumed that there were no chemical reactions; diffusion was the predominant mechanism, so there was no convective mass transfer into the cheese pores, only diffusion; the effective or apparent diffusion coefficients were constant (did not depend on concentration or coordinates); and the sample contraction was negligible. Under these conditions, multicomponent diffusion during salting in brine with agitation at a constant temperature of 10°C was modeled. With well-stirred brine we can assume that the inter-phase concentration equilibrium of the brine/cheese is instantaneously established after immersion of the cheese (Luna and Bresan, 1986).

The classical theory of transport phenomena is difficult to apply to biological materials because of their peculiar characteristics, structure, properties, etc. (Welti-Chanes et al., 2005). The effect of food microstructure on the diffusion coefficient is discussed in Aguilera and Stanley (1999). Actually, even considering cheese as a food matrix with saline solution occluded in the pores, parameters such as porosity, tortuosity, and phase ratios are not sufficient to describe the extraction/impregnation process accurately (Crossley and Aguilera, 2001). It is very probable that the quantification of food microstructures using image analysis will result in a more fundamental approach to studying mass transfer phenomena in foods and other biological materials (Aguilera, 2005). Meanwhile, many research publications on mass transfer in food engineering are oriented to solving Ficks law or other derivations of this law for some specific process (Welti-Chanes et al., 2005).

The C_{1} (x, y,z, t) and C_{2} (x, y, z, t) concentrations of the solutes NaCl and KCl, respectively, defined at point P º (x, y, z) Î W and at instant t, can be described by the Onsager equations (1945) for the solute concentrations:

where Ñ^{2}(.) = Ñ . Ñ (.) is the Laplacian operator and D_{ii} are the main diffusion coefficients and D_{ij} are the crossed diffusion coefficients that combine the flows.

The boundary conditions used for the salting process in well-stirred brine are, in mathematical terms, essential and are called Dirichlet boundary conditions and are given by

The initial conditions used were

where *C _{1,S }*and

*C*are respectively the NaCl and KCl concentrations in brine and ¶W is the set of points on the surface of the cheese, conditions (2).

_{2,S}*C*and

_{1,0 }*C*are respectively the initial NaCl and KCl concentrations in the cheese, conditions (3).

_{2,0}The partial differential equation system formed by equation (1) and by conditions (2) and (3) is difficult to solve analytically, even under the assumptions made regarding the process. One way of solving it is by using a transformation obtained with eigenvalues and eigenvectors in a matrix of differential operators that appears on the right-hand side of equation (1). The disadvantages of this method are mainly due to restrictions in system geometry. In addition, the transformed variables have no physical meaning (Bona et al., 2005).

Another alternative that does not have these restrictions and can also be applied in this study is the numerical solution of the system by a finite element method. A generalized Galerkin formulation was used (Buchanan, 1994; Huebner et al., 1995) to obtain the equation to be discretized by finite elements. The main steps are synthesized as follows:

1. The equation for the weighted residues formed for the system (1) is

where f and Y are weighting functions that are null in ¶W;

2. The right side of equation (4) is integrated by parts to obtain the following:

The discretization continues following the next steps:

3. A partition of the domain in to some subdomains called finite elements is adopted;

4. N nodes are attributed to each element, and each node has two degrees of freedom (one for each solute). The variables of the problem are interpolated locally by the interpolation equations associated with the nodes. Therefore, (5) equation can be rewritten in matrix form (local stiffness matrix) for each element as

where

5. The global stiffness matrix is the assemblage of all local stiffness matrices and it results in a system of ordinary differential equations;

6. The solution of this system is obtained by an implicit process of finite differences supplying the values of the concentrations at the element nodes. Local interpolation can determine the concentration profile at any point in the cheese at any time during salting.

**Spatial Discretization**

The spatial domain was represented by a set of hexahedral elements of the C^{0 }type, shown in Figure 4. Each element has 20 nodes with two degrees of freedom (one for each solute) on the intersecting lines and vertices, totaling forty degrees of freedom (Silva et al., 1998). Figure 4a shows the nodes and the local enumeration. The spatial orientation of the elements and enumeration of the surfaces used in the computer procedure are shown in Figures 4b and 4c. The element is defined in standard Cartesian coordinates, that is -1 __<__ x, h, z __<__ 1 (Chung, 1978).

The interpolation functions, represented by equations 7 to 10, preserve the continuity among the elements (Chung, 1978; Akin, 1982). The interpolation functions used for the nodal points of the vertices, i = 1, 2, 3, ..., 8 were

at midside nodes, i = 9, 10, 11, ..., 20.

where f_{i} is the isoparametric function (Brebbia and Ferrante, 1975; Chung, 1978; Cook, 1981; Bickford, 1990) and x, h, and z are local coordinates of the standard element. In the formulation proposed Y_{i} = f_{i} was used.

**Procedure for Integration in Time**

The Crank-Nicholson implicit algorithm, or the central difference scheme was used for time discretization. For application of the implicit scheme in the domain, and **C _{2N}** were discretized from equation (6) as follows:

where Dt is a discrete increase in time; is the solute concentration vector, which is already known for each node of the element before adding time; and represents the concentration vector that will be calculated after adding time (Dt). Organizing the terms of the equation (11) we have

Equation (12) represents the system of equations, written in matrix form, for each element. By overlapping the elements, a global system of linear equations of the **AX** = **B** type is obtained. The boundary conditions are imposed next, and the system can be solved by any technique compatible with the size and characteristics of matrix **A**.

**Solving the System of Linear Equations**

When choosing a method for this system resolution, care should be taken with respect to two problems (Sperandio et al., 2003): the first concerns the propagation of round-off errors, that is, the numerical stability of the method compared to its sensitivity in accumulating round-off errors; the second concerns the form of storage of matrix **A**, which should be in line with the algorithm used.

The LU decomposition method (Sperandio et al., 2003) was used to solve the system of linear equations . For a linear system of equations, LU factorization produces **A** = **LU** where **L** and **U** are the lower and upper triangular matrices, respectively. Therefore **LUX** = **B**, and assigning **UX** = **Y,** another linear system is formed as

The **L** and **U** matrices were generated using the Doolittle method (Sperandio et al., 2003). The computational implementation of the method was carried out by a set of subroutines (DLSLRB and DLFCRB) based on the LINPACK package (Cline et al., 1979; Dongarra et al*.*, 1979). The DLSLRB subroutine solves linear equation systems stored in a band without using interactive refinement. This subroutine uses another subroutine (DLFCRB), which estimates the condition number of the matrix and performs the LU decomposition, and the DLFSRB routine solves the system formed by the **L** and **U** matrices.

**Mesh**

A mesh generating subroutine was applied because the practical application of the finite element method requires a lot of input data. The global enumeration and the coordinates of the nodes were obtained automatically for rectangular domains. This procedure simplified node position identification, allowing sampling at any point or region of the mesh, making the computer program more versatile (Bona, 2004).

**Computational Programming **

Simul 3.0 software (Bona, 2004) in Fortran language, configured to apply the Dirichlet boundary condition, was used to simulate multicomponent diffusion during the salting with agitation. The minimal requirement for its use is a Pentium III processor and 256Mb of ram memory.

**Integration and Interpolation**

The experimental values obtained represented the mean NaCl and KCl concentrations for the cylindrical sample (Figure 3). As FEM gives point concentrations, it was necessary to interpolate the simulated points on the imaginary axis located at the center of the sample removed. The average was calculated by the integration of these points using the following equation:

The integral was calculated numerically by a Gauss quadrature formula (Sperandio et al., 2003) and the abscissas and weights were obtained from Davis and Polonsky (1965).

**Statistical Test**

The estimated and experimental values of the NaCl and KCl concentrations in the brine were compared by percentage deviation (Heldman, 1974), represented by equation , in order to assess the goodness of fit.

where _{est} is the average concentration estimated by the numerical solution, _{exp} is the experimental average concentration, and N is the number of observations considered.

**RESULTS AND DISCUSSION**

The initial brine and cheese characteristics are shown in Tables 1 and 2, respectively.

The parameters used (three steps) were obtained by optimization of the spatial temporal discretization applied to the simulation of multicomponent diffusion (Bona, 2004): (i) the simulation by finite elements was performed with a mesh of 504 elements; (ii) two 0.05h, three 0.30h, and six 1.00h intervals were used for the time discrete partitions because shorter time intervals at the start are essential (Lyra, 1993) to dampen the oscillations that damage FEM refinement; and (iii) the effective diffusion coefficients (Table 3) used in the simulation were obtained from Bona (2004) for brine at rest. Slightly different values for the diffusion coefficients were found through a one-dimensional analytical solution (Bona et al., 2005). However the deviations among the respective main coefficients were comparable with those found in published data (Zorrilla and Rubiolo, 1994). The fact that the KCl main coefficient was larger than the NaCl coefficient was due to its greater mobility in an aqueous solution (Schwartzberg and Chao, 1982). The main coefficients were greater than the cross coefficients, indicating that diffusion in relation to its own gradient was more important than diffusion due to interaction among solutes (Zorrilla and Rubiolo, 1994).

Figure 5 shows the experimental and simulated saline distribution profiles of NaCl and KCl diffusion in Prato cheese using brine with agitation. The average values for both faces (Figure 3) were used as experimental data.

Salting in brine at rest is different from salting with agitation due to the existence of an external resistance that hinders diffusion of the salts in the cheese. This resistance is inversely proportional to an adimensional number called the mass transfer Biot number (Loncin and Merson, 1979; Fryer et al., 1997). It can be expressed as

where h_{m} (g/cm^{2 }h) is the film mass transfer coefficient, l_{m} (g/cm h) is the mass conductivity, and *R _{i}* is the characteristic length (cm). For brine with agitation this resistance can be neglected (Luna and Bressan, 1986).

Simul 3.0 software was also used to assess the efficiency of the agitation, assuming that the salting was static. Results were simulated with Biots of 33.154 (Bona, 2004), 50, and 100. Table 4 shows the percentage of deviation for different Biot numbers using a Cauchy boundary condition, equation (17), and for the absence of external resistance using the Dirichlet boundary condition in the salting simulation.

The simulation with no resistance compared to that with a Biot number of 100 (Table 4) indicated that the agitation applied to the brine in this case was sufficient to reduce external resistance in a way that did not interfere with saline diffusion in the cheese. These statements are in line with Schwartzberg and Chao (1982), who reported that Biot numbers over 100 mean that the diffusion process is little influenced by external resistance.

The theoretical profiles for the well-stirred brine (no external resistance) fall rapidly toward the initial concentration at values less than 1.2cm (Figure 6). This fact needs to be taken into account when sampling cheese to obtain experimental data for comparison with theoretical data.

The curve shapes in Figure 6 are in line with other experiments on brining cheeses such as Emmental (Pajonk et al., 2003), Cuartirolo (Luna and Bressan, 1986), and Romano (Guinne and Fox, 1984).

The oscillation verified up to the fifth hour is normally expected when FEM is applied to diffusion problems dealing with food solids (Silva et al., 1999) because mass diffusivity is very low, on the order of 10^{-9} m^{2}/s (Loncin and Merson, 1974). Furthermore, the concentration gradient is very large during all brining processes. In spite of this, the oscillation was reduced applying the Crank-Nicolson scheme and spatial-time refinements (Bona, 2004).

The average concentration of the salts in whole cheeses can be calculated with equation (14). This procedure permits estimation of the ideal salting time for a desired salt content. Figure 7 shows a simulation of the total saline profile (NaCl and KCl) as a function of time, obtained from the data generated by the Simul 3.0 software. According to Furtado (1990), the salt content of Prato cheese consumed in Brazil is normally 2.0%. Based on this information, salting without agitation would take about 5.8 hours, which is twice the time necessary to salt the same Prato cheese when the brine is agitated (Figure 7).

**CONCLUSIONS**

By FEM it was possible to simulate multicomponent diffusion during the salting of Prato cheese in well-stirred brine. The data simulated were in line with experimental results, validating the application of FEM in this case. The deviation of experimental concentration data was 7.3% for NaCl and 5.4% for KCl, both of which were considered acceptable. For sodium and potassium chloride it was evident that an increase in the Biot number reduced the deviation from the experimental data. The system developed to simulate salt diffusion will allow for the control and modulation of salt content in Prato cheese (widely consumed in Brazil), predicting the final content from the initial conditions. Thus, a better quality cheese with proper salt concentrations can be produced by solely controlling the salting time. The possible reduction in salting time will allow an increase in productivity without additional investments in the process.

**ACKNOWLEDGEMENTS**

The authors thank Fundação Araucaria for its financial support and CAPES for the scholarships.

NOMENCLATURE | ||

Bi | Biot mass-exchange number | (-) |

average concentration | C (g/100g _{NaCl} + KCl + water) | |

_{est} | average concentration estimated by FEM | (g/100g _{NaCl }+ KCl + water) |

_{exp} | experimental average concentration | (g/100g _{NaCl} + KCl + water) |

C_{i,0} | initial concentration of component i | (g/100g _{NaCl} + KCl + water) |

C_{i,s} | component i concentration in brine | (g/100g _{NaCl} + KCl + water) |

C_{i} (x,y,z,t) | specific mass concentration | (g/100g _{NaCl} + KCl + water) |

D_{ii} | effective main diffusion coefficient | (cm^{2}/day) |

D_{ij} | effective cross-diffusion coefficient | (cm^{2}/day) |

h_{m} | film mass transfer coefficient | (g/cm^{2} h) |

N | number of observations used to calculate %deviation | (-) |

p | statistical significance | (p-level) |

±R_{1}, ±R_{2}, ±R_{3} | characteristic length | (cm) |

t | time | (hours) |

x, y, z | spatial coordinates | (cm) |

D^{2}(.) | Laplacian operator | (-) |

¶/¶h | normal derivative operator | (-) |

¶W | set of surface points of the cheese | (-) |

f(W)dW integration on the domain | (-) | |

f(V)dV integration on the volume | (-) | |

%deviation percentage of deviation between predicted and experimental data | (-) |

*Greek Symbols*

f, Y | weighting or interpolation functions | (-) |

l_{m} | mass conductivity | (g/cm h) |

x, h, z | standard Cartesian coordinates | (-) |

x_{i}, h_{i}, z_{i} | nodal standard cartesian coordinates | (-) |

W | cheese volume or domain | (cm^{3}) |

*Sub Indices*

i | component i |

j | component j |

1 | NaCl |

2 | KCl |

**REFERENCES**

Aguilera, J.M., Why food Microstructure? Journal of Food Engineering, vol. 67, pp. 3-11 (2005). [ Links ]

Aguilera, J.M. and Stanley, D.W., Microstructural Principles of Food Processing and Engineering. 2^{nd} ed. Gaitherburg, MD: Aspen Publishers, Inc (1999). [ Links ]

Akin, J.E., Application Implementation of Finite Element Methods. London: Academic Press, 372p. (1982). [ Links ]

AOAC – Association of Official Agricultural Chemists, Official Methods of Analysis, 14^{th} ed., 1141p. (1984). [ Links ]

Bickford, W.B., A First Course in the Finite Element Method. Boston: Irwin, 649p. (1990). [ Links ]

Bona, E., Salga Mista com NaCl e KCl em Queijo Prato – Modelagem da Difusão Multicomponente com o Emprego do Método de Elementos Finitos. Masterss thesis, Universidade Estadual de Londrina, UEL, Brazil, 120p. (2004). [ Links ]

Bona, E., Borsato, D., Silva, R.S.S.F., and Silva, L.H.M. Difusão multicomponente durante a salga mista de queijo prato. Ciência e Tecnologia de Alimentos, vol. 25, no. 2, pp. 394-400 (2005). [ Links ]

Brebbia, C.A. and Ferrante, A.J., The Finite Element Technique: An Introduction for Engineers. Porto Alegre: Ed. UFRGS, 410p. (1975). [ Links ]

Britz, D., Østerby, O., and Strutwolf, J., Damping of Crank-Nicolson error oscillations. Computational Biology and Chemistry, vol. 27, pp. 253-263 (2002). [ Links ]

Buchanan, G.R., Schaums Outline of Theory and Problems of Finite Element Analysis. New York: McGraw-Hill, 264p. (1994). [ Links ]

Cecchi, H.M., Fundamentos Teóricos e Práticos em Análise de Alimentos. Campinas: Editora da Unicamp, 212p. (1999). [ Links ]

Chung, T.J., Finite Element Analysis in Fluid Dynamics. New York: McGraw Hill, 378p. (1978). [ Links ]

Cline, A.K., Moler, C.B., Stewart, G.W., and Wilkinson, J.H., An estimate for the condition number of a matrix. SIAM Journal of Numerical Analysis, vol. 16, pp. 368-375 (1979). [ Links ]

Cook R.D., Analysis, Concepts and Application of Finite Element. New York: John Wiley & Sons, 2^{nd} ed., 535p. (1981). [ Links ]

Crossley, J. and Aguilera, J.M., Modelling the effect of microstructure on food extraction. Journal of Food Process Engineering, vol. 24, pp. 161-177 (2001). [ Links ]

Cussler, E.L., Multicomponent Diffusion. Amsterdam: Elsevier, 176p. (1976). [ Links ]

Davis, P.J. and Polonsky, I., Numerical interpolation, differentiation and integration. In: Abramowitz, M. and Stegun, I. (eds.), Handbook of Mathematical Functions. New York: Dover Publications, chap. 25, pp. 875-974 (1965). [ Links ]

Dongarra, J.J., Bunch, J.R., Moler, C.B., and Stewart, G.W., LINPACK Users Guide. Philadelphia: SIAM (1979). [ Links ]

Fryer, P.J., Pyle, D.L., and Rielly, C.D., Chemical Engineering for the Food Industry. London: Blackie Academic & Professional, 462p. (1997). [ Links ]

Furtado, M.M., A Arte e a Ciência do Queijo. São Paulo: Editora Globo S.A., 297p. (1990). [ Links ]

Gerla, P.E. and Rubiolo, A.C., A model for determination of multicomponent diffusion coefficients in foods. Journal of Food Engineering, vol. 56, pp. 401-410 (2003). [ Links ]

Guinne, T.P. and Fox, P.F., Salt in cheese: Physical, chemical and biological aspects. In: Fox, P.F. (ed.), Cheese: Chemistry, Physics and Microbiology. New York: Elsevier Applied Publ. Science, vol. 1, pp. 31-35 (1984). [ Links ]

He, F.J. and MacGregor, G.A., Beneficial effects of potassium. British Medical Journal, vol. 323, pp. 497-501 (2001). [ Links ]

Heldman, D.R., Predicting the relationships between unfrozen water fraction and temperature during food freezing using freezing point depression. Transactions ASAE, vol. 17, p. 63 (1974). [ Links ]

Huebner, K.H., Thorton, E.A., and Byron, T.G., The Finite Element Method for Engineers. New York: John Wiley & Sons, 3^{rd} ed. (1995). [ Links ]

Jameson, G.W., Dietary cheese: low fat, low salt. CSIRO Food Research Journal, vol. 46, p. 64 (1986). [ Links ]

Loncin, M. and Merson, R.L., Food Engineering Principles and Selected Applications. New York: Academic Press, 494p. (1979). [ Links ]

Luna, J.A. and Bressan, J.A., Mass transfer during brining of Cuartirolo Argentino cheese. Journal of Food Science, vol. 51, no. 3, pp. 829-831 (1986). [ Links ]

Lyra, P.R.M., Finite elements analysis of parabolic problems: combined influence of adaptive mesh refinement and automatic time step control. Revista Brasileira de Ciências Mecânicas, vol. 15, no. 2, pp. 172-198 (1993). [ Links ]

Onsager, L., Theories and problems of liquid diffusion. Ann. N. Y. Acad. Sci., vol. 46, p. 241 (1945). [ Links ]

Pajonk, A.S., Saurel, R., and Andrieu, J., Experimental study and modeling of effective NaCl diffusion: coefficients values during Emmental cheese brining. Journal of Food Engineering, vol. 60, no. 4, pp. 307-313 (2003). [ Links ]

Puri, V.M. and Anantheswaran, R.C., The finite-element method in food processing: A review. Journal of Food Engineering, vol. 19, pp. 247-274 (1993). [ Links ]

Rapacci, M., Efeito da substituição do cloreto de sódio por cloreto de potássio nas características do queijo prato. Masters thesis, Universidade Estadual de Londrina, Brazil, 106p. (1989). [ Links ]

Richardson, G.H., ed. Standard Methods for Examination of Dairy Products. 5^{th} ed. IS.ed. Washington: American Public Health Association, 412p. (1985). [ Links ]

Schwartzberg, H.G. and Chao, R.Y., Solute diffusivities in leaching processes. Food Technology, vol. 36, no. 2, pp. 73-86 (1982). [ Links ]

Silva, R.S.S.F., Borsato, D., and Silva, L. H. M. Difusão do cloreto de sódio no processo de salga de queijos: modelagem matemática com o emprego do método de elementos finitos. Ciência e Tecnologia de Alimentos, vol. 18, no. 1, pp. 87-92 (1998). [ Links ]

Silva, L.H.M., Silva, R.S.S.F., and Borsato, D. Finite element analysis of diffusion problem during cheese salting: combined influence of space and time discretization. Acta Scientiarum, vol. 21, no. 4, pp. 873-879 (1999). [ Links ]

Sperandio, D., Mendes, J.T., and Silva, L.H.M. Cálculo Numérico: Características Matemáticas e Computacionais dos Métodos Numéricos. São Paulo: Prentice Hall, 354p. (2003). [ Links ]

Wang, L. and Sun, D., Recent developments in numerical modelling of heating or cooling processes in the food industry – A review. Trends in Food Science & Technology, vol. 14, pp. 408-423 (2003). [ Links ]

Welti-Chanes, J., Vergara-Balderas, F., and Bermúdez-Aguirre, D., Transport phenomena in food engineering: Basic concepts and advances. Journal of Food Engineering, vol. 67, pp. 113-128 (2005). [ Links ]

Zorrilla, S.E. and Rubiolo, A.C., A model for using the diffusion cell in the determination of multicomponent diffusion coefficients in gels or food. Chemical Engineering Science, vol. 49, no. 13, pp. 2123-2128 (1994). [ Links ]

(Received: October 15, 2005 ; Accepted: April 18, 2007)

* To whom correspondence should be addressed