Mathematical modeling of drying process of unripe banana slices

Unripe banana flour (UBF) production employs bananas not submitted to maturation process, is an interesting alternative to minimize the fruit loss reduction related to inappropriate handling or fast ripening. The UBF is considered as a functional ingredient improving glycemic and plasma insulin levels in blood, have also shown efficacy on the control of satiety, insulin resistance. The aim of this work was to study the drying process of unripe banana slabs (Musa cavendishii, Nanicão) developing a transient drying model through mathematical modeling with simultaneous moisture and heat transfer. The raw material characterization was performed and afterwards the drying process was conducted at 40 ºC, 50 ºC e 60 ºC, the product temperature was recorded using thermocouples, the air velocity inside the chamber was 4 m·s-1. With the experimental data was possible to validate the diffusion model based on the Fick\'s second law and Fourier. For this purpose, the sorption isotherms were measured and fitted to the GAB model estimating the equilibrium moisture content (Xe), 1.76 [g H2O/100g d.b.] at 60 ºC and 10 % of relative humidity (RH), the thermophysical properties (k, Cp, ?) were also measured to be used in the model. Five cases were contemplated: i) Constant thermophysical properties; ii) Variable properties; iii) Mass (hm), heat transfer (h) coefficient and effective diffusivity (De) estimation 134 W·m-2·K-1, 4.91x10-5 m-2·s-1 and 3.278?10-10 m·s-2 at 60 ºC, respectively; iv) Variable De, it presented a third order polynomial behavior as function of moisture content; v) The shrinkage had an effect on the mathematical model, especially in the 3 first hours of process, the thickness experienced a contraction of about (30.34 ± 1.29) % out of the initial thickness, finding two decreasing drying rate periods (DDR I and DDR II), 3.28x10-10 m·s-2 and 1.77x10-10 m·s-2, respectively. COMSOL Multiphysics simulations were possible to perform through the heat and mass transfer coefficient estimated by the mathematical modeling.

Liquid-liquid extraction in a pilot scale rotating disc contactor

A study of the hydrodynamics and mass transfer characteristics of a liquid-liquid extraction process in a 450 mm diameter, 4.30 m high Rotating Disc Contactor (R.D.C.) has been undertaken. The literature relating to this type of extractor and the relevant phenomena, such as droplet break-up and coalescence, drop mass transfer and axial mixing has been revjewed. Experiments were performed using the system C1airsol-350-acetone-water and the effects of drop size, drop size-distribution and dispersed phase hold-up on the performance of the R.D.C. established. The results obtained for the two-phase system C1airso1-water have been compared with published correlations: since most of these correlations are based on data obtained from laboratory scale R.D.C.'s, a wide divergence was found. The hydrodynamics data from this study have therefore been correlated to predict the drop size and the dispersed phase hold-up and agreement has been obtained with the experimental data to within +8% for the drop size and +9% for the dispersed phase hold-up. The correlations obtained were modified to include terms involving column dimensions and the data have been correlated with the results obtained from this study together with published data; agreement was generally within +17% for drop size and within +14% for the dispersed phase hold-up. The experimental drop size distributions obtained were in excellent agreement with the upper limit log-normal distributions which should therefore be used in preference to other distribution functions. In the calculation of the overall experimental mass transfer coefficient the mean driving force was determined from the concentration profile along the column using Simpson's Rule and a novel method was developed to calculate the overall theoretical mass transfer coefficient Kca1, involving the drop size distribution diagram to determine the volume percentage of stagnant, circulating and oscillating drops in the sample population. Individual mass transfer coefficients were determined for the corresponding droplet state using different single drop mass transfer models. Kca1 was then calculated as the fractional sum of these individual coefficients and their proportions in the drop sample population. Very good agreement was found between the experimental and theoretical overall mass transfer coefficients. Drop sizes under mass transfer conditions were strongly dependant upon the direction of mass transfer. Drop Sizes in the absence of mass transfer were generally larger than those with solute transfer from the continuous to the dispersed phase, but smaller than those with solute transfer in the opposite direction at corresponding phase flowrates and rotor speed. Under similar operating conditions hold-up was also affected by mass transfer; it was higher when solute transfered from the continuous to the dispersed phase and lower when direction was reversed compared with non-mass transfer operation.

Drying of portland cement raw material slurries:predicting the performance of a counter-current spray drier

An extensive review of literature has been carried out concerning the drying of single drops, sprays of droplets and the prediction of spray drier performances. The experimental investigation has been divided into two broad parts mainly: (1) Single Drop Experiments, and (2) Spray Drying and Residence Time Distribution Experiments. The thermal conductivity of slurry cakes from five different sources have been experimentally determined using a modified Lee's Disc Apparatus and the data collected was correlated by the polynominal... Good agreement was observed between the experimental thermal conductivity values and the predicted ones. The fit gave a variance ... for the various samples experimented on. A mathematical model for estimating crust mass transfer coefficient at high drying temperatures was derived.

Studies of the drying of particulate slurries

The literature relating to the drying characteristics of pure liquid drops and particulate slurry drops has been reviewed. The experimental investigation was, therefore, divided into three parts: Pure water drops, Aqueous sodium sulphate decahydrate drops, and, Slurry drops from nine detergent formulations. The value of the constant,'Ψ, reported by Ranz and Marshall, was found to be temperature dependent. In the temperature range o 26.5≤T≤118.5°C,Ψ , for pure water drops, varied between 0.38 and 0.47. A revised correlation of the mass transfer coefficients is therefore proposed.  A mathematical model for estimating the variation of crust thickness, for aqueous sodium sulphate drops, with time is proposed: β = R _ {R3 - ( 1.5G/πCo ) ( ΔHD - ΔHU) Δ} 1/3 Experimental crust thickness evaluated from stereoscan micrographs showed good agreement with theoretical prediction. It has been shown that drying characteristics of detergent drops can be evaluated from the porosity:thickness ratio, {ε/\β}. Formulations having large {ε/β I-ratios dry better than those with smaller values. The agreement between the experimental and theoretical mass transfer coefficients shows, in addition to the above correlation, that the overall mass transfer coefficient can be predicted from the expression1/K = 1/K + β/DMε 1.5 The crust is the controlling resistance to transfer in particulate slurry drops. For aqueous sodium sulphate drops, the crust provides 64.2% of the total resistance while for detergents with thicker, but less porous crusts, the value is 97.5%.

Liquid-phase adsorption of n-paraffins on type 5A molecular sieves

As a basis for the commercial separation of normal paraffins a detailed study has been made of factors affecting the adsorption of binary liquid mixtures of high molecular weight normal paraffins (C12, C16, and C20) from isooctane on type 5A molecular sieves. The literature relating to molecular sieve properties and applications, and to liquid-phase adsorption of high molecular weight normal paraffin compounds by zeolites, was reviewed. Equilibrium isotherms were determined experimentally for the normal paraffins under investigation at temperatures of 303oK, 323oK and 343oK and showed a non-linear, favourable- type of isotherm. A higher equilibrium amount was adsorbed with lower molecular weight normal paraffins. An increase in adsorption temperature resulted in a decrease in the adsorption value. Kinetics of adsorption were investigated for the three normal paraffins at different temperatures. The effective diffusivity and the rate of adsorption of each normal paraffin increased with an increase in temperature in the range 303 to 343oK. The value of activation energy was between 2 and 4 kcal/mole. The dynamic properties of the three systems were investigated over a range of operating conditions (i.e. temperature, flow rate, feed concentration, and molecular sieve size in the range 0.032 x 10-3 to 2 x 10-3m) with a packed column. The heights of adsorption zones calculated by two independent equations (one based on a constant width, constant velocity and adsorption zone and the second on a solute material balance within the adsorption zone) agreed within 3% which confirmed the validity of using the mass transfer zone concept to provide a simple design procedure for the systems under study. The dynamic capacity of type 5A sieves for n-eicosane was lower than for n-hexadecane and n-dodecane corresponding to a lower equilibrium loading capacity and lower overall mass transfer coefficient. The values of individual external, internal, theoretical and experimental overall mass transfer coefficient were determined. The internal resistance was in all cases rate-controlling. A mathematical model for the prediction of dynamic breakthrough curves was developed analytically and solved from the equilibrium isotherm and the mass transfer rate equation. The experimental breakthrough curves were tested against both the proposed model and a graphical method developed by Treybal. The model produced the best fit with mean relative percent deviations of 26, 22, and 13% for the n-dodecane, n-hexadecane, and n-eicosane systems respectively.

A study of liquid-liquid extraction in a sieve plate column

The literature relating to sieve plate liquid extraction columns and relevant hydrodynamic phenomena have been surveyed. Mass transfer characteristics during drop formation, rise and coalescence, and related models were also reviewed. Important design parameters i.e. flooding, dispersed phase hold-up, drop size distribution, mean drop size, coalescence/flocculation zone height beneath a plate and jetting phenomena were investigated under non-mass transfer and mass transfer conditions in a 0.45m diameter, 2.3m high sieve plate column. This column had provision for four different plate designs, and variable plate spacing and downcomer heights, and the system used was Clairsol `350' (dispersed) - acetone - deionised water (continuous) with either direction of mass transfer. Drop size distributions were best described by the functions proposed by Gal-or, and then Mugele-Evans. Using data from this study and the literature, correlations were developed for dispersed phase hold-up, mean drop size in the preferred jetting regime and in the non-jetting regime, and coalescence zone height. A method to calculate the theoretical overall mass transfer coefficient allowing for the range of drop sizes encountered in the column gave the best fit to experimental data. This applied the drop size distribution diagram to estimate the volume percentage of stagnant, circulating and oscillating drops in the drop population. The overall coefficient Kcal was then calculated as the fractional sum of the predicted individual single drop coefficients and their proportion in the drop population. In a comparison between the experimental and calculated overall mass transfer coefficients for cases in which all the drops were in the oscillating regime (i.e. 6.35mm hole size plate), and for transfer from the dispersed(d) to continuous(c) phase, the film coefficient kd predicted from the Rose-Kintner correlation together with kc from that of Garner-Tayeban gave the best representation. Droplets from the 3.175mm hole size plate, were of a size to be mainly circulating and oscillating; a combination of kd from the Kronig-Brink (circulating) and Rose-Kintner (oscillating) correlations with the respective kc gave the best agreement. The optimum operating conditions for the SPC were identified and a procedure proposed for design from basic single drop data.

Concentration polarization model of spiral-wound membrane modules with application to batch-mode RO desalination of brackish water

Batch-mode reverse osmosis (batch-RO) operation is considered a promising desalination method due to its low energy requirement compared to other RO system arrangements. To improve and predict batch-RO performance, studies on concentration polarization (CP) are carried out. The Kimura-Sourirajan mass-transfer model is applied and validated by experimentation with two different spiral-wound RO elements. Explicit analytical Sherwood correlations are derived based on experimental results. For batch-RO operation, a new genetic algorithm method is developed to estimate the Sherwood correlation parameters, taking into account the effects of variation in operating parameters. Analytical procedures are presented, then the mass transfer coefficient models are developed for different operation processes, i.e., batch-RO and continuous RO. The CP related energy loss in batch-RO operation is quantified based on the resulting relationship between feed flow rates and mass transfer coefficients. It is found that CP increases energy consumption in batch-RO by about 25% compared to the ideal case in which CP is absent. For continuous RO process, the derived Sherwood correlation predicted CP accurately. In addition, we determined the optimum feed flow rate of our batch-RO system.

Estudo da absorção das espécies pesadas do gás natural em octanol: efeitos da temperatura e vazão

In the oil industry, natural gas is a vital component of the world energy supply and an important source of hydrocarbons. It is one of the cleanest, safest and most relevant of all energy sources, and helps to meet the world's growing demand for clean energy in the future. With the growing share of natural gas in the Brazil energy matrix, the main purpose of its use has been the supply of electricity by thermal power generation. In the current production process, as in a Natural Gas Processing Unit (NGPU), natural gas undergoes various separation units aimed at producing liquefied natural gas and fuel gas. The latter should be specified to meet the thermal machines specifications. In the case of remote wells, the process of absorption of heavy components aims the match of fuel gas application and thereby is an alternative to increase the energy matrix. Currently, due to the high demand for this raw gas, research and development techniques aimed at adjusting natural gas are studied. Conventional methods employed today, such as physical absorption, show good results. The objective of this dissertation is to evaluate the removal of heavy components of natural gas by absorption. In this research it was used as the absorbent octyl alcohol (1-octanol). The influence of temperature (5 and 40 °C) and flowrate (25 and 50 ml/min) on the absorption process was studied. Absorption capacity expressed by the amount absorbed and kinetic parameters, expressed by the mass transfer coefficient, were evaluated. As expected from the literature, it was observed that the absorption of heavy hydrocarbon fraction is favored by lowering the temperature. Moreover, both temperature and flowrate favors mass transfer (kinetic effect). The absorption kinetics for removal of heavy components was monitored by chromatographic analysis and the experimental results demonstrated a high percentage of recovery of heavy components. Furthermore, it was observed that the use of octyl alcohol as absorbent was feasible for the requested separation process.

Estudo da absorção das espécies pesadas do gás natural em octanol: efeitos da temperatura e vazão

In the oil industry, natural gas is a vital component of the world energy supply and an important source of hydrocarbons. It is one of the cleanest, safest and most relevant of all energy sources, and helps to meet the world's growing demand for clean energy in the future. With the growing share of natural gas in the Brazil energy matrix, the main purpose of its use has been the supply of electricity by thermal power generation. In the current production process, as in a Natural Gas Processing Unit (NGPU), natural gas undergoes various separation units aimed at producing liquefied natural gas and fuel gas. The latter should be specified to meet the thermal machines specifications. In the case of remote wells, the process of absorption of heavy components aims the match of fuel gas application and thereby is an alternative to increase the energy matrix. Currently, due to the high demand for this raw gas, research and development techniques aimed at adjusting natural gas are studied. Conventional methods employed today, such as physical absorption, show good results. The objective of this dissertation is to evaluate the removal of heavy components of natural gas by absorption. In this research it was used as the absorbent octyl alcohol (1-octanol). The influence of temperature (5 and 40 °C) and flowrate (25 and 50 ml/min) on the absorption process was studied. Absorption capacity expressed by the amount absorbed and kinetic parameters, expressed by the mass transfer coefficient, were evaluated. As expected from the literature, it was observed that the absorption of heavy hydrocarbon fraction is favored by lowering the temperature. Moreover, both temperature and flowrate favors mass transfer (kinetic effect). The absorption kinetics for removal of heavy components was monitored by chromatographic analysis and the experimental results demonstrated a high percentage of recovery of heavy components. Furthermore, it was observed that the use of octyl alcohol as absorbent was feasible for the requested separation process.

Evaluation of thermo-physical properties and drying kinetics of carrots in a convective hot air drying system

This paper presents an experimental study on the evolution of carrot properties along convective drying by hot air at different temperatures (50ºC, 60ºC and 70ºC). The thermo-physical properties calculated were: specific heat, thermal conductivity, diffusivity, enthalpy, heat and mass transfer coefficients. Furthermore, the data of drying kinetics were treated and adjusted according to the three empirical models: Page, Henderson & Pabis and Logarithmic. The sorption isotherms were also determined and fitted using the GAB model. The results showed that, generally, the thermo-physical properties presented a decline during the drying process, and the decrease was faster for the temperature of 70ºC. It was possible to verify that the Page model presented the best prediction ability for the representation of kinetics of the drying process. The GAB model used to fit the sorption isotherms showed a good prediction capacity and, at a given water activity, despite some variations, the amount of water sorbed increased with the decrease of drying temperature.

Flow boiling heat transfer of R134a and R245fa in a 2.3 mm tube

This paper presents new experimental flow boiling heat transfer results in micro-scale tubes. The experimental data were obtained in a horizontal 2.3 mm I.D stainless steel tube with heating length of 464 mm, R134a and R245fa as working fluids, mass velocities ranging from 50 to 700 kg m(-2) s(-1), heat flux from 5 to 55 kW m(-2), exit saturation temperatures of 22, 31 and 41 degrees C, and vapor qualities ranging from 0.05 to 0.99. Flow pattern characterization was also performed from images obtained by high-speed filming. Heat transfer coefficient results from 1 to 14 kW m(-2) K(-1) were measured. It was found that the heat transfer coefficient is a strong function of heat flux, mass velocity and vapor quality. The experimental data were compared against ten flow boiling predictive methods from the literature. Liu and Winterton [3], Zhang et al. [5] and Saitoh et al. [6] worked best for both fluids, capturing most of the experimental heat transfer trends. (C) 2010 Elsevier Ltd. All rights reserved.

Experimental and numerical investigation of dynamic heat transfer parameters in packed bed

Dynamic experiments in a nonadiabatic packed bed were carried out to evaluate the response to disturbances in wall temperature and inlet airflow rate and temperature. A two-dimensional, pseudo-homogeneous, axially dispersed plug-flow model was numerically solved and used to interpret the results. The model parameters were fitted in distinct stages: effective radial thermal conductivity (K (r)) and wall heat transfer coefficient (h (w)) were estimated from steady-state data and the characteristic packed bed time constant (tau) from transient data. A new correlation for the K (r) in packed beds of cylindrical particles was proposed. It was experimentally proved that temperature measurements using radially inserted thermocouples and a ring-shaped sensor were not distorted by heat conduction across the thermocouple or by the thermal inertia effect of the temperature sensors.

Numerical and ExperimentalModelling of Gas Flow and Heat Transfer in the Air Gap of An Electric Machine

This work deals with the cooling of high-speed electric machines, such as motors and generators, through an air gap. It consists of numerical and experimental modelling of gas flow and heat transfer in an annular channel. Velocity and temperature profiles are modelled in the air gap of a high-speed testmachine. Local and mean heat transfer coefficients and total friction coefficients are attained for a smooth rotor-stator combination at a large velocity range. The aim is to solve the heat transfer numerically and experimentally. The FINFLO software, developed at Helsinki University of Technology, has been used in the flow solution, and the commercial IGG and Field view programs for the grid generation and post processing. The annular channel is discretized as a sector mesh. Calculation is performed with constant mass flow rate on six rotational speeds. The effect of turbulence is calculated using three turbulence models. The friction coefficient and velocity factor are attained via total friction power. The first part of experimental section consists of finding the proper sensors and calibrating them in a straight pipe. After preliminary tests, a RdF-sensor is glued on the walls of stator and rotor surfaces. Telemetry is needed to be able to measure the heat transfer coefficients at the rotor. The mean heat transfer coefficients are measured in a test machine on four cooling air mass flow rates at a wide Couette Reynolds number range. The calculated values concerning the friction and heat transfer coefficients are compared with measured and semi-empirical data. Heat is transferred from the hotter stator and rotor surfaces to the coolerair flow in the air gap, not from the rotor to the stator via the air gap, althought the stator temperature is lower than the rotor temperature. The calculatedfriction coefficients fits well with the semi-empirical equations and precedingmeasurements. On constant mass flow rate the rotor heat transfer coefficient attains a saturation point at a higher rotational speed, while the heat transfer coefficient of the stator grows uniformly. The magnitudes of the heat transfer coefficients are almost constant with different turbulence models. The calibrationof sensors in a straight pipe is only an advisory step in the selection process. Telemetry is tested in the pipe conditions and compared to the same measurements with a plain sensor. The magnitudes of the measured data and the data from the semi-empirical equation are higher for the heat transfer coefficients than thenumerical data considered on the velocity range. Friction and heat transfer coefficients are presented in a large velocity range in the report. The goals are reached acceptably using numerical and experimental research. The next challenge is to achieve results for grooved stator-rotor combinations. The work contains also results for an air gap with a grooved stator with 36 slots. The velocity field by the numerical method does not match in every respect the estimated flow mode. The absence of secondary Taylor vortices is evident when using time averagednumerical simulation.

Modeling and simulation of forced-air cooling of strawberries using variable convective coefficient

In the forced-air cooling process of fruits occurs, besides the convective heat transfer, the mass transfer by evaporation. The energy need in the evaporation is taken from fruit that has its temperature lowered. In this study it has been proposed the use of empirical correlations for calculating the convective heat transfer coefficient as a function of surface temperature of the strawberry during the cooling process. The aim of this variation of the convective coefficient is to compensate the effect of evaporation in the heat transfer process. Linear and exponential correlations are tested, both with two adjustable parameters. The simulations are performed using experimental conditions reported in the literature for the cooling of strawberries. The results confirm the suitability of the proposed methodology.

Solid mass falling into the reservoirs shallow water equations validity from the engineering point of view

Water waves generated by a solid mass is a complex phenomenon discussed in this paper by numerical and experimental approaches. A model based on shallow water equations with shocks (Saint Venant) has developed. It can reproduce the amplitude and the energy of the wave quite well, but because it consistently generates a hydraulic jump, it is able to reproduce the profile, in the case of high relative thickness of slide, but in the case of small relative thickness it is unable to reproduce the amplitude of the wave. As the momentum conservation is not verified during the phase of wave creation, a second technique based on discharge transfer coefficient α, is introduced at the zone of impact. Numerical tests have been performed and validated this technique from the experimental results of the wave's height obtained in a flume.