Prediction of drying of dairy products using a modified balance-based desorption method

Similar documents
A new simple recursive algorithm for finding prime numbers using Rosser s theorem

Water Vapour Effects in Mass Measurement

Methylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma.

Smart Bolometer: Toward Monolithic Bolometer with Smart Functions

Case report on the article Water nanoelectrolysis: A simple model, Journal of Applied Physics (2017) 122,

On size, radius and minimum degree

Can we reduce health inequalities? An analysis of the English strategy ( )

Evolution of the cooperation and consequences of a decrease in plant diversity on the root symbiont diversity

Vibro-acoustic simulation of a car window

Easter bracelets for years

On Symmetric Norm Inequalities And Hermitian Block-Matrices

Solubility prediction of weak electrolyte mixtures

Soundness of the System of Semantic Trees for Classical Logic based on Fitting and Smullyan

Quantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry

Thomas Lugand. To cite this version: HAL Id: tel

Multiple sensor fault detection in heat exchanger system

Lorentz force velocimetry using small-size permanent magnet systems and a multi-degree-of-freedom force/torque sensor

Ultra low frequency pressure transducer calibration

Completeness of the Tree System for Propositional Classical Logic

Dispersion relation results for VCS at JLab

From Unstructured 3D Point Clouds to Structured Knowledge - A Semantics Approach

Impulse response measurement of ultrasonic transducers

Interactions of an eddy current sensor and a multilayered structure

IMPROVEMENTS OF THE VARIABLE THERMAL RESISTANCE

Passerelle entre les arts : la sculpture sonore

SOLAR RADIATION ESTIMATION AND PREDICTION USING MEASURED AND PREDICTED AEROSOL OPTICAL DEPTH

Full-order observers for linear systems with unknown inputs

L institution sportive : rêve et illusion

On the longest path in a recursively partitionable graph

b-chromatic number of cacti

AC Transport Losses Calculation in a Bi-2223 Current Lead Using Thermal Coupling With an Analytical Formula

Entropies and fractal dimensions

Solving the neutron slowing down equation

DEM modeling of penetration test in static and dynamic conditions

LAWS OF CRYSTAL-FIELD DISORDERNESS OF Ln3+ IONS IN INSULATING LASER CRYSTALS

A remark on a theorem of A. E. Ingham.

Comment on: Sadi Carnot on Carnot s theorem.

Mirage detection for electrochromic materials characterization. Application to iridium oxide films

Natural convection of magnetic fluid inside a cubical enclosure under magnetic gravity compensation

Numerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral. Equation Approach

Best linear unbiased prediction when error vector is correlated with other random vectors in the model

Theoretical calculation of the power of wind turbine or tidal turbine

The influence of the global atmospheric properties on the detection of UHECR by EUSO on board of the ISS

Sound intensity as a function of sound insulation partition

A new approach of the concept of prime number

The generation of the Biot s slow wave at a fluid-porous solid interface. The influence of impedance mismatch

On production costs in vertical differentiation models

Analysis of Boyer and Moore s MJRTY algorithm

Exogenous input estimation in Electronic Power Steering (EPS) systems

New Basis Points of Geodetic Stations for Landslide Monitoring

RHEOLOGICAL INTERPRETATION OF RAYLEIGH DAMPING

Diurnal variation of tropospheric temperature at a tropical station

On The Exact Solution of Newell-Whitehead-Segel Equation Using the Homotopy Perturbation Method

FIM OBSERVATION OF MONOLAYER Pd ADSORBED ON W AND Mo SURFACES

The FLRW cosmological model revisited: relation of the local time with th e local curvature and consequences on the Heisenberg uncertainty principle

Improvement of YBCO Superconductor Magnetic Shielding by Using Multiple Bulks

On Newton-Raphson iteration for multiplicative inverses modulo prime powers

Antipodal radiation pattern of a patch antenna combined with superstrate using transformation electromagnetics

STATISTICAL ENERGY ANALYSIS: CORRELATION BETWEEN DIFFUSE FIELD AND ENERGY EQUIPARTITION

ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS

Trench IGBT failure mechanisms evolution with temperature and gate resistance under various short-circuit conditions

HIGH RESOLUTION ION KINETIC ENERGY ANALYSIS OF FIELD EMITTED IONS

On Symmetric Norm Inequalities And Hermitian Block-Matrices

A Simple Proof of P versus NP

A simple kinetic equation of swarm formation: blow up and global existence

Unfolding the Skorohod reflection of a semimartingale

The sound power output of a monopole source in a cylindrical pipe containing area discontinuities

Fast Computation of Moore-Penrose Inverse Matrices

Territorial Intelligence and Innovation for the Socio-Ecological Transition

One Use of Artificial Neural Network (ANN) Modeling in Simulation of Physical Processes: A real-world case study

Application of computational fluid dynamics to spray drying

Structural study of a rare earth-rich aluminoborosilicate glass containing various alkali and alkaline-earth modifier cations

Exact Comparison of Quadratic Irrationals

A novel method for estimating the flicker level generated by a wave energy farm composed of devices operated in variable speed mode

Numerical modeling of diffusion within composite media

Some approaches to modeling of the effective properties for thermoelastic composites

Determination of absorption characteristic of materials on basis of sound intensity measurement

Spatial representativeness of an air quality monitoring station. Application to NO2 in urban areas

On Poincare-Wirtinger inequalities in spaces of functions of bounded variation

Stickelberger s congruences for absolute norms of relative discriminants

There are infinitely many twin primes 30n+11 and 30n+13, 30n+17 and 30n+19, 30n+29 and 30n+31

Characterization of the local Electrical Properties of Electrical Machine Parts with non-trivial Geometry

Basic concepts and models in continuum damage mechanics

A numerical analysis of chaos in the double pendulum

A CONDITION-BASED MAINTENANCE MODEL FOR AVAILABILITY OPTIMIZATION FOR STOCHASTIC DEGRADING SYSTEMS

Particle-in-cell simulations of high energy electron production by intense laser pulses in underdense plasmas

A Study of the Regular Pentagon with a Classic Geometric Approach

Towards an active anechoic room

AN INTERNAL FRICTION PEAK DUE TO HYDROGEN-DISLOCATION INTERACTION IN NICKEL

Question order experimental constraints on quantum-like models of judgement

On the link between finite differences and derivatives of polynomials

Stochastic invariances and Lamperti transformations for Stochastic Processes

On path partitions of the divisor graph

Comparison of Harmonic, Geometric and Arithmetic means for change detection in SAR time series

Nonlocal computational methods applied to composites structures

Impairment of Shooting Performance by Background Complexity and Motion

On Solving Aircraft Conflict Avoidance Using Deterministic Global Optimization (sbb) Codes

Visible laser emission of Pr3+ in various hosts

NEGATIVE ION IMAGING IN FIELD ION MICROSCOPY

Some explanations about the IWLS algorithm to fit generalized linear models

Transcription:

Prediction of drying of dairy products using a modified balance-based desorption method Peng Zhu, Serge Méjean, Eric Blanchard, Romain Jeantet, Pierre Schuck To cite this version: Peng Zhu, Serge Méjean, Eric Blanchard, Romain Jeantet, Pierre Schuck. Prediction of drying of dairy products using a modified balance-based desorption method. Dairy Science & Technology, EDP sciences/springer, 213, 93 (4), pp.347-355. <1.17/s13594-12-99-9>. <hal-121418> HAL Id: hal-121418 https://hal.archives-ouvertes.fr/hal-121418 Submitted on 17 Sep 215 HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers. L archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Dairy Sci. & Technol. (213) 93:347 355 DOI 1.17/s13594-12-99-9 ORIGINAL PAPER Prediction of drying of dairy products using a modified balance-based desorption method Peng Zhu & Serge Méjean & Eric Blanchard & Romain Jeantet & Pierre Schuck Received: 17 September 212 / Revised: 12 November 212 / Accepted: 22 November 212 / Published online: 1 January 213 # INRA and Springer-Verlag France 212 Abstract A desorption method based on a thermodynamic approach has recently been developed to evaluate the behaviour of dairy concentrates during drying. Involving overall heat and mass balance throughout the dryer, this approach can determine several key gas-feed parameters for industrial spray-drying processes. Spray-drying software (SD 2 P ) was then designed following this approach to predict the optimal inlet drying air temperatures with acceptable accuracy (95 99% accuracy) for spray-drying of dairy products. However, the mass change of the sample was indirectly determined from the change in the relative humidity of the air during desorption, which could be a source of error below the detection threshold of the thermo-hygrometric sensor. In order to measure directly the mass change during drying by desorption, a modified drying-by-desorption method was investigated in this study. The novel method used a precise microbalance and a modified desorption cell which permitted measurement of the mass change of the sample and the relative humidity of the air at the same time. Different materials (water, skim milk, infant formulae, etc.) were tested using this new method. The results obtained with direct (microbalance) and indirect (thermo-hygrometer) measurements were found to be highly consistent (coefficient of determination 1). This confirmed that the mass change estimation by the original desorption method was correct overall under current conditions. Moreover, the proposed new desorption method makes it possible to monitor water transfer with constant accuracy over the entire desorption process, P. Zhu : E. Blanchard Laiterie de Montaigu, 856 Montaigu, France E. Blanchard LISA.RD, 856 Montaigu, France S. Méjean : R. Jeantet : P. Schuck (*) INRA, UMR1253, 35 Rennes, France e-mail: Pierre.Schuck@rennes.inra.fr S. Méjean: R. Jeantet : P. Schuck Agrocampus Ouest, UMR1253, 35 Rennes, France

348 P. Zhu et al. thus permitting reliable study of the mass transfer phenomena throughout experiments. Keywords Desorption. Spray-drying. Water transfer. Dairy product 1 Introduction Spray-drying is a widely used technique to remove water from dairy products. The principle of this technique is to spray a concentrated liquid in individual small droplets and then mix them with hot air (temperature ranging from 1 to 3 C). Evaporation of water from the droplets is facilitated by the temperature and partial pressure gradients between the droplet and hot air. The temperature at the surface of the droplets is believed to be close to the wet bulb temperature, i.e., in the range of 3 to 5 C, and the total duration of drying is only a few seconds (Schuck et al. 1998; 29). This technique can therefore provide reasonable preservation of the physiochemical and nutritional quality of the products. Due to the increasing need for stringent product quality and strict function requirements, different dairy concentrates with a variety of constituents are currently spraydried. For example, the formulation of an infant formula may involve the mixing of over 4 different constituents before the evaporation and spray-drying processing stages. The differences in product composition may result in different drying behaviours, and the drying of such products may require adjustment of the spray-drying parameters. An inappropriate drying operation may lead to non-compliant product quality (physiochemical: water activity; nutritional: Maillard reaction; physical: cakiness and stickiness) and substantial economic losses. There are two main methods to obtain optimal drying parameters for spray-drying of food products: the first is through empirical pilot experiments, but this method is usually expensive, time consuming and sometimes unreliable. The other method to determine the optimal drying parameters is through spray-drying modelling based on drying kinetics (Patel et al. 21). A drying-by-desorption method and Spray-drying parameter simulation and determination software (SD2P software) were developed in 25 at the Institut National de la Recherche Agronomique in France (Schuck et al. 29). This desorption method using an air-tight stainless-steel cylindrical cell (Fig. 1) to dry a milk concentrate sample (16±1 mg) contained in a small plastic cup at 45 C. The desorption cell was filled with 12±1 ml zeolite particles (Zhu et al. 211a). Water then transferred from the concentrate to the zeolites by means of the pressure gradient between the concentrate and the zeolites. Changes in the relative humidity of the air inside the cell during desorption were continuously monitored using a relative humidity sensor placed close to the surface of the milk sample. The change in relative humidity of the air between the concentrate and the zeolites can be measured by a thermo-hygrometric sensor connected to a computer, and the drying rate of the concentrate can be calculated from the desorption curve using the SD2P software. By integrating several parameters of the dryer, the concentrate and the final product requirements (e.g. evaporation capacity, air flow rates and humidity, total solid contents and temperature of the

Modified balance-based desorption method 349 Relative humidity sensor Cup Zeolites Fig. 1 Diagram of the original desorption experimental setup concentrate, water content of the powder in relation to water activity, energy cost, etc.), the main operational parameters of spray-drying (e.g. inlet and outlet air temperatures, concentrate and powder flow rates, specific energy consumption, yield of the dryer and cost (per kg water evaporated or per kg powder produced)) can then be predicted from the mass and energy balance. Validation tests of over 3 different food concentrates have been performed using different spraydryers (evaporation capacity from 5 kg.h 1 to 6 t.h 1 ), and a good match between measured and predicted parameters (±1 5% error) has been reported (Schuck et al. 29). In this method, the drying rate is evaluated indirectly from the change in the relative humidity by considering that the area below the relative humidity curve is representative of the quantity of evaporated water. However, the relative humidity values measured at the final stage of desorption are very low and close to the detection threshold of the thermo-hygrometric sensor. This could be a source of error for the calculation of mass change values. The aim of the current study was therefore to modify the existing desorption method to confirm whether such an indirect desorption method can predict the drying kinetics accurately during drying by desorption. 2 Materials and methods 2.1 Preparation of the concentrates Concentrates with 4% w/w total solids were prepared from dairy powders dissolved in distilled water. The skim milk powder was produced at Bionov (Rennes, France). The infant formula powder was provided by the Laiterie de Montaigu (Montaigu, France).

35 P. Zhu et al. 2.2 Determination of total solid content Total solid content of the concentrates was measured according to the weight loss profiles after oven drying of 1 g of the sample mixed with sand at 15 C for 7 h. 2.3 Drying-by-desorption method A new desorption cell has been designed to measure, at the same time, mass change and relative humidity change during desorption (Fig. 2). Similar to the original desorption cell (Fig. 1), the new desorption cell is filled with 12±1 ml zeolite particles. A plastic cup containing the concentrate sample (16±1 mg) is placed in the middle of the cell and on the weighing platform of the microbalance (max 22 g, d=.1 mg, model: wxts25du, Mettler Toledo, Viroflay, France). A relative humidity sensor (Rotronic, Bassersdorf, Switwerland) is placed on top of this desorption cell. The drying temperature of the original desorption method was 45 C (Schuck et al. 29; Zhu et al. 211a). However, due to the temperature range limit of the microbalance, the desorption temperature used in this study was fixed at 4 C. The duration of desorption was 3 min. The changes in mass of the concentrate and of the relative humidity are recorded by computer every minute during drying by desorption. The desorption of skim milk was performed twice. The general standard deviations were 1.65 mg and.6% for measured mass and relative humidity, respectively. Using this new desorption cell, mass change can be obtained directly from the microbalance or indirectly from the relative humidity data. Calculation of the mass change from relative humidity was based on the hypothesis that the area between the measured relative humidity curve and the bulk relative humidity (RH b ) curve represents the amount of water evaporated during desorption (Fig. 3). The total water evaporated during desorption can be calculated from the difference between the initial (m )andthefinal(m f )mass weighed before and after desorption drying. Having established the total amount of water evaporated, the instantaneous water evaporated (m t )attimet can thus be calculated from the relative humidity at time t (RH t ) and the cumulated relative humidity (RH cumulated, area between the relative humidity curve and the Relative humidity sensor Microbalance Cup Zeolites Desorption cell Weighing platform Fig. 2 Diagram of the new experimental desorption setup

Modified balance-based desorption method 351 Measured relative humidity (%) 4 3 2 1 RH cumulated 1 2 3 Time (min) Fig. 3 An example of the desorption curve (line: measured relative humidity of 4% w/w total solids skim milk; broken line: bulk relative humidity of zeolite; RH cumulated : cumulated relative humidity) bulk relative humidity), Eq. (1): m t ¼ RH t RH b ðm m f Þ ð1þ RH cumulated Then the mass change (m) during the whole desorption drying process can be calculated from 1 to 3 min, Eq. (2): m ¼ X3 m t t¼1 The drying rate is calculated separately from measured mass change data and calculated mass change during the same desorption process. The mass change data measured by microbalance was taken every 1 min from the recorded raw data for the calculation of the drying rate. Then it is possible to compare the indirect method (calculate from relative humidity) with the direct method (measure by microbalance). ð2þ 3 Results Figure 4 shows the linear correlation between the measured and calculated mass changes obtained by the microbalance and relative humidity sensor, respectively, for distilled water (Fig.4a), 4% w/w total solids skim milk (Fig. 4b) and 4% w/w total solids infant formula (Fig. 4c) concentrates. Good correlation was obtained between the measured and calculated mass for water and the two concentrates under the drying conditions tested. The average differences between the measured and calculated mass values were.58±.42,.3±.55 and 1.3±.75 mg for distilled water, skim milk and infant formula concentrates, respectively. Figure 5 shows the drying rate of desorption calculated from the measured and calculated mass changes of distilled water (Fig. 5a), 4% w/w total solids skim milk (Fig. 5b) and 4% w/w total solids infant formula (Fig. 5c) concentrates. For each of the three trials, the two curves exhibit approximately the same form. It can be seen that for distilled water, there was only one drying period up to 263 min. Desorption of distilled water was at a fast drying rate during this period, with an average value of

352 P. Zhu et al. Calculated mass (mg) Calculated mass (mg) Calculated mass (mg) 2 15 1 5 2 15 1 y =.99x +.96 R² = 1 5 1 15 2 Measured mass (mg) 5 2 15 1 a y =.98x + 1.66 R² = 1 5 1 15 2 Measured mass (mg) 5 y =.99x + 2.51 R² = 1 5 1 15 2 Measured mass (mg) Fig. 4 Correlation between calculated and measured mass during drying by desorption of distilled water (a), 4% w/w total solids skim milk (b) and 4% w/w total solids infant formula (c) b c.61 mg.min 1 and the highest drying rate of.64 mg.min 1. The drying rate of distilled water decreased slightly during desorption and reached approximately.5 mg.min 1 at the end (263 min). In contrast, two different periods can easily be distinguished for desorption of skim milk and infant formula concentrates from the drying rate curve: first a fast drying rate period (approximately from 11 to 11 min for skim milk and from 11 to 12 min for infant formula), and then a slow drying rate period (for duration higher than 11 and 12 min for skim milk and infant formula, respectively). The average drying rates in the fast drying rate period were.56 mg.min 1 whatever the product, and the fastest drying rates were.61 and.58 mg.min 1 for skim milk and infant formula concentrates, respectively. The

Modified balance-based desorption method 353 Drying rate (mg.min -1 ) Drying rate (mg.min -1 ) Drying rate (mg.min -1 ).7.6 a.5.4.3.2.1 1 2 3 Time (min).7.6.5.4.3.2.1 1 2 3 Time (min).7.6.5.4.3.2.1 1 2 3 Time (min) Fig. 5 Drying rate profiles obtained during drying by desorption of distilled water (a), 4% w/w total solids skim milk (b) and 4% w/w total solids infant formula (c) (diamond: by microbalance, cross: by relative humidity sensor) b c average drying rates of the slow drying rate period were.16 and.13 mg.min 1 for skim milk and infant formula concentrates, respectively. 4 Discussion The first 1 min of desorption can be considered as the preheating stage (Zhu et al. 211b). From 11 to 263 min, the desorption of distilled water showed a fast drying rate (Fig. 5a). After this period, no more water was left in the plastic cup. As the thickness of the sample is only about 1 mm, the slight decrease in drying rate may be

354 P. Zhu et al. explained by the reduction in the liquid surface when there is not enough water to cover the bottom of the plastic cup entirely due to the surface tension of water. For the desorption of skim milk and infant formula concentrates (Fig. 5b, c), the existence of the fast and the slow dying rate period has been previously reported and explained by the relative easiness /difficulty to remove water, according to its binding to dairy constituents (Zhu et al. 211b). These two periods have been observed and reported as well for single droplet drying at different drying conditions: higher drying air temperature (7 11 C) to dry a lower initial solids concentration (1 2% w/w) (Lin and Chen 26; Fu et al. 211). This fast drying rate period was once considered as the constant drying rate stage (Zhu et al. 211b). However, the results showed that the drying rate decreased during this period, mainly because of the increases in concentration and viscosity during evaporation of water. The sudden decrease in drying rate in this period may relate to the formation of a crust at the surface of the concentrate along with an increase in viscosity and a decrease of free water to remove (Zhu et al. 211b). The drying temperature used in this study was 4 C for both of the direct and indirect drying-by-desorption methods due to the temperature range limit of the microbalance. For all three trials, the results show a good correlation between the measured and calculated mass changes under the drying conditions tested (Fig. 4). The differences between measured and calculated mass values may be due to the accuracy of the relative humidity sensor, which is limited by temperature variation at the preheating stage and by low relative humidity values at the final stage of desorption drying. In contrast, using a microbalance gave constant accuracy over the entire desorption process, even at the final stage. The difference between measured and calculated mass values was never exceeding 1.3±.75 mg for the concentrates tested in this study. It is reasonable to believe that the original desorption method using a relative humidity sensor can be used to evaluate correctly the mass change of a concentrate during drying by desorption up to 3 min, even at 45 C. It is a simple and reliable method to provide drying kinetics of a concentrate for the SD 2 P software. 5 Conclusion A desorption method has been developed to evaluate the drying kinetics of a concentrate for a spray-drying software (SD 2 P, Schuck et al. 29). This method uses a relative humidity sensor to determine indirectly the mass change of a concentrate during drying by desorption. In order to evaluate the reliability of the original desorption method, a modified drying-by-desorption method was developed in this study. The modified method can measure directly the mass change during the drying by desorption process. Different materials (water, skim milk, infant formulae) were tested using this new method. The results obtained with direct (microbalance) and indirect (relative humidity sensor) measurements were found to be highly consistent under the drying conditions tested. This confirmed that the estimation of mass change by the original desorption method was correct overall under current conditions. The original desorption method using a relative humidity sensor is a simple and reliable method to provide drying kinetics for the SD 2 P software.

Modified balance-based desorption method 355 However, the proposed new desorption method makes it possible to monitor water transfer with constant accuracy over the entire desorption process, thus permitting further reliable study of mass transfer phenomena throughout experiments. Acknowledgements The authors thank the Laiterie de Montaigu (France) for their scientific and financial support for this study. References Fu N, Woo MW, Lin XQ, Zhou Z, Chen XD (211) Reaction engineering approach (REA) to model the drying kinetics of droplets with different initial sizes experiments and analyses. Chem Eng Sci 66:1738 1747 Lin XQ, Chen XD (26) A model for drying of an aqueous lactose droplet using the reaction engineering approach. Dry Technol 24:1329 1334 Patel K, Chen XD, Jeantet R, Schuck P (21) One-dimensional simulation of co-current, dairy spray drying systems pros and cons. Dairy Sci Technol 9:181 21 Schuck P, Roignant M, Brulé G, Davenel A, Famelart MH, Maubois JL (1998) Simulation of water transfer in spray drying. Dry Technol 16(7):1371 1393 Schuck P, Dolivet A, Méjean S, Zhu P, Blanchard E, Jeantet R (29) Drying by desorption: a tool to determine spray drying parameters. J Food Eng 94:199 24 Zhu P, Méjean S, Blanchard E, Jeantet R, Schuck P (211a) Prediction of dry mass glass transition temperature and the spray drying behaviour of a concentrate using a desorption method. J Food Eng 15:46 467 Zhu P, Patel K, Lin S, Méjean S, Blanchard E, Chen XD, Schuck P, Jeantet R (211b) Simulating industrial spray drying operations using a reaction engineering approach and a modified desorption method. Dry Technol 29(4):419 428