Next Article in Journal
Testing and Analysis of Torsional Vibration of Ship Transmission Shafting Based on Five-Point Smoothing Algorithm
Next Article in Special Issue
Antioxidant and Sensorial Properties: Meat Analogues versus Conventional Meat Products
Previous Article in Journal
Supercritical Fluid Extraction from Zataria multiflora Boiss and Impregnation of Bioactive Compounds in PLA for the Development of Materials with Antibacterial Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Study of a Microbial Spray-Drying Process Based on Real-Time Sampling

1
College of Mechanical & Electrical Engineering, Shaanxi University of Science & Technology, Xi’an 710016, China
2
College of Land and Environment, Shenyang Agricultural University, Shenyang 110866, China
*
Author to whom correspondence should be addressed.
Processes 2022, 10(9), 1789; https://doi.org/10.3390/pr10091789
Submission received: 6 August 2022 / Revised: 1 September 2022 / Accepted: 2 September 2022 / Published: 5 September 2022
(This article belongs to the Special Issue Bio-Active Compounds in Food Production)

Abstract

:
The process of microbial spray-drying has inherent defects such as short time, complexity, and non-visualization of particle trajectory. However, there has been a lack of effective methods for real-time sampling, rehydration, and non-destructive storage of dried particles, as well as mathematical modeling of the drying process of yeast particles based on sampling and measurement data. Therefore, firstly, a real-time sampling system was developed which completed real-time sampling, rehydration, and non-destructive storage of spray-dried particles, and realized the real reproduction of the changes of yeast particles in the process of spray drying. The laws that the motion trajectory of microbial particles during spray drying are divided into the first cycle region and the reflux cycle region were concluded, and the partition mechanism was explored. Then, based on the sampling data and the law of heat and mass transfer, a mathematical model of porous media was established to predict the variation of moisture content and rehydration survival rate of dried microbial particles with drying time. Finally, the mathematical model was tested by a microbial spray drying experiment, and the maximum errors between the predicted value and the test value of moisture content and rehydration survival rate were X max 1 = 0.027 ( d . b ) and Q v max 1 = 1.06 % , respectively, both were less than 5% which proved the correctness of the mathematical model of porous media and laid a foundation for the study of the damage mechanism of microbial spray drying.

Graphical Abstract

1. Introduction

Microbial desiccation is the core research topic of livelihood fields (such as food, medicine, and etc.). As the most widely used microbial drying process, spray drying has the advantages of high efficiency, low cost, and continuous production. It not only occupies a dominant position in the market, but also is a hot spot in current research [1]. Nevertheless, due to the inherent defects of the microbial spray drying process, such as short time, complexity, and invisible trajectory of spray particles [2], it is very difficult to directly sample and study the drying process of microbial particles. Therefore, there has been a lack of effective means for real-time sampling, rehydration, and non-destructive storage of dried particles, and further establishing mathematical models to accurately predict the moisture content and rehydration survival rate of microbial particles in the spray drying process [3]. At present, the common research method is “single droplet simulation experiment + mathematical modeling” [4]. On the one hand, the single droplet simulation experiment could not reproduce environmental components such as short residence time, complex movement trajectory, and high-speed collision of spray-dried particles, and the experimental data could not truly reflect the changes of microorganisms in the spray-drying process. On the other hand, most mathematical models ignore changes in the trajectory of a dry particle and simplify the trajectory to top-down homogeneous motion [5]. Therefore, the mathematical models established by this combination method could not predict the spray-drying process of microorganisms well.
Dr. Gong Pi-min [6], (Harbin Institute of Technology), proposed a vertical arrangement of a single-cup straight cylinder sampling system. Although this system could complete the real-time sampling of spray-dried particles, it suffered from the disadvantages of poor sampling stability and nondestructive storage, which could not reproduce the real microbial changes during spray-drying. Previous research on modeling has mainly focused on regression analysis to build empirical models, such as Akpinar [7], equal to the thin layer model proposed in 2006, there was little in-depth exploration of heat and mass transfer and microbial cell damage of the drying process, so the prediction accuracy was not high [8]. In 2013, L. Spreutels [9], a Canadian scholar, established an image-only mathematical model [10], based on the dynamic drying curve method based on the single droplet simulation experiment, although the spray particles were assumed to fall at a uniform speed in the spray drying process. However, the model was based on the single drop test, so the rehydration survival of yeast particles during drying could not be accurately predicted as a function of drying time [11].
In response to the above problems, firstly a novel microbial spray-drying real-time sampling system (submitted for a Chinese invention patent) was developed which could implement real-time sampling, rehydration, and nondestructive storage on spray-dried particles. Then, on the basis of real-time sampling data, the heat and mass transfer model of porous media proposed by the research team in the previous study was used to complete the mathematical modeling of the spray-drying process. Finally, through the comparison with and the analysis of the experimental data, it was proved that the model could accurately predict the moisture content and rehydration survival rate of dried particles with time in the process of microbial spray drying.

2. A Brief Introduction of the New Sampling System

2.1. A Brief Introduction of the Sampling Arrangement of Double Wet and Dry Samples in Step-Ladder Type

In the novel sampling scheme, a dry-pellet sampling cup and a rehydration sample sampling cup, 50 mm apart, were arranged on the same ladder step and were responsible for collecting the dry pellets and real-time rehydration of the dry pellets, respectively, with nondestructive storage. The drying tower was then bisected equally by height and perimeter into 6 elevations and 7-segment arcs, respectively, with each ladder step again arranged at the intermediate value of height and radians, completing the 7-stage sampling arrangement. In order to prevent the spray particles from directly entering the sampling cup and affecting the sampling stability, arc-shaped baffles were installed at 30 mm above the mouth of each sampling cup. Finally, the design of the stepped-type dry and wet double-sample sampling arrangement was completed. This sampling scheme, with a step ladder-type design, avoids interference between the individual sampling cups in a vertical-type arrangement scheme, which, combined with the installation of arcuated baffles, greatly improves the stability and accuracy of sampling.

2.2. A Brief Review of Two Cup Surface Sampling Cups

In terms of sampling principle, the sampling cup attracts dry particles loaded with microorganisms and carboxyl iron powder (CIP) into the sampling cup through the magnetic field of a strong magnet at the bottom of the outer cup, so as to achieve the purpose of stable sampling [12].
At constant temperature and humidity, the sampling cup was filled with −40 °C antifreeze for mobility in the freezing layer, and then put into the −60 °C refrigerator for quick freezing to form a solid state, so as to maintain the low-temperature environment in the inner cup. At the same time, the thermal and magnetic exchange between inside and outside of the sampling cup was isolated by pasting the thermal insulation magnetic patch evenly on the outer surface of the outer cup, and the low-temperature environment in the inner cup was further maintained, while the strong magnetic field was prevented from penetrating the outer cup and affecting the sampling process.
In order to improve the sampling stability and the capacity of the water to hold a constant temperature at the same time, the sampling cup was improved to design a straight cylinder inner cup to be a convex curved inner cup by rotating the first 20 mm of the Bezier cubic curve [13], with control points P 0 = ( 0 , 0 ) , P 1 = ( 12 , 9 ) , P 2 = ( 24 , 16 ) and P 3 = ( 48 , 8 ) . The curve group is shown in Figure 1 below, and the equation of curve group is shown in Equations (1)–(3) below [14]:
P ( μ ) = i = 0 n B i n ( μ ) P i
B i n ( μ ) = n ! i ! ( n 1 ) ! μ i ( 1 μ ) n i
P ( μ ) = μ 0 ( 1 μ ) 3 P 0 + 3 μ ( 1 μ ) 2 P 1 + 3 μ 2 ( 1 μ ) P 2 + μ 3 ( 1 μ ) 0 P 3 = ( 1 μ ) 3 P 0 + 3 μ ( 1 μ ) 2 P 1 + 3 μ 2 ( 1 μ ) P 2 + μ 3 P 3
In the Equations (1)–(3): n—degree of polynomial; i—non-negative integer; P i —control point representing curve equation; and B i n ( μ ) —The basis function of the nth order curve equation.
The curved inner cup structure changed the air flow field, wind speed, and pressure field in the cup and reduced the inflow of external hot air (as shown in Figure 2). Finally, the purpose of improving the constant temperature water holding capacity of the sampling cup without increasing the depth of the sampling cup to ensure the sampling capacity and stability was achieved.
In order to achieve non-destructive installation, the sampling cup was fixed by connecting a magnet and spray-drying tower to achieve the purpose of non-destructive installation. In addition, the four exposed surfaces of the magnet were affixed with thermal insulation magnetic patches to prevent them from affecting the spray-drying sampling process.
In order to measure temperature and humidity in real time, the sampling cup recorded the sampling temperature and humidity at the sampling point in real time through a pasted wireless temperature and humidity sensor placed in the outer cup mouth. The curve of temperature and humidity change was drawn.
The overall structure diagram of the sampling cup is shown in Figure 3 below. The dimensions of the sampling cup are: the upper surface of the inner cup is 20 mm high, the minimum inner diameter is 18 mm; the lower part of the inner cup is 10 mm high and the inner diameter is 24 mm; the height of the outer cup is 65 mm, the outer diameter is 42 mm; and the wall thickness is 2 mm.
A large number of experiments have proven that the combination of the double cup curved sampling cup and step-ladder type sampling scheme could achieve stable real-time sampling, rehydration, and non-destructive storage of spray-drying particles, and the sampling data could accurately reflect the changes of microorganisms in the process of spray drying.

3. Experimental Design and Data Analysis

In this test, Baijiu yeast was sampled by a novel sampling system using three parameters: spray-dried material flow, hot air wind speed, and temperature as test variables, and the temperature and humidity at each sampling site were recorded. Then, the dry-based water content, moisture status and particle size of the dried pellets in each dried pellet sampling cup were determined, respectively, and the rehydration survival rate of yeasts in each rehydrated sample sampling cup were determined. Finally, the determination parameters were analyzed to obtain the movement law of microbial particles during spray drying.

3.1. Materials and Methods

The schematic diagram of the spray-drying test device is shown in Figure 4. The left figure is the overall structure diagram, and the right figure is the structure diagram of sampling and temperature and humidity detection in the drying tower. The parameters of the spray-drying tower used in the experiment are tower height h = 900 mm, cylinder diameter d = 400 mm, conical taper at the bottom a = 60°, and tower bottom diameter d = 150 mm.
As shown in Figure 4 above, the temperature and humidity measurements of the sampling point were completed by the temperature and humidity sensors arranged at the mouth of the sampling cup in cooperation with the data acquisition system. The dry base moisture content, water state, particle size, and rehydration survival rate of samples in the sampling cup were determined by water meter, differential scanning calorimeter (DSC), Winner3001 dry powder laser particle size meter, and hemocytometry, respectively.

3.2. Test Procedure

The purpose of this test is to preliminarily analyze the changes of microorganisms represented by Baijiu Yeast in the drying process by measuring the values of various evaluation indexes of samples in different sampling cups under various test variables, so as to lay a foundation for the establishment of a mathematical model for predicting the drying process [15].

3.2.1. Preparation of Bacterial Solution to Be Dried

According to the research team’s prior studies: ① The optimized N medium formulation were yeast extract 15 g/L, peptone 15 g/L, beef extract 15 g/L, fructose 40 g/L, KH2PO4 3 g/L, and inositol 0.1 g/L. ② Yeast desiccant protectant formulations were sucrose 100 g/L, trehalose 100 g/L, lactose 120 g/L, and sorbitol 80 g/L. ③ The optimal vehicles when dried were reconstituted skimmed milk powder (RSM) at a concentration of 30% in culture broth [16]. ④ Carboxyiron powder (CIP) was added as follows: CIP: RSM = 1:8 [17].
The steps to formulate the broth to be dried are: ① Baijiu Yeast was inoculated in 12 groups of 200 mL of novel n medium at an inoculation ratio of 0.5%, then placed into a constant temperature shaking incubator at 28 °C and incubated at a rotation rate of 180 r/min for 9 h. ② Aliquots were centrifuged for 10 min in a centrifuge at 4000 r/min to remove supernatant, remove bacteria, and wash twice by centrifugation in normal saline. ③ The collected bacteria were re-suspended in 20 mL PBS buffer. ④ To the cell suspension, 180 mL of culture solution containing 30% RSM was added, and then CIP was added at a CIP: RSM ratio of 1:8 to mix well. ⑤ According to the above formula of protective agents, protective agents were added, and then placed in a 4 °C freezer for refrigeration.

3.2.2. Spray-Drying Test Design

The experimental procedures for spray-drying were: ① Before each test, the components of the spray-drying system were cleaned in full accordance with the cleaning requirements of the spray-drying system. Then the spray-drying system was preheated by predetermined test inlet hot air temperature and flow rate for 10 min to obtain stable test conditions. ② Set the feed pump parameters by the predetermined mass flow and then start the feed pump, starting the test. ③ The temperature and humidity values of each sampling point displayed by the collection system were recorded. After the test was completed, the samples in the sampling cup were immediately taken out and placed in the refrigerator at 4 °C after numbering. ④ The dry-basis water content, moisture status, particle size, and rehydration survival of each sample were determined separately.
First, the yeast particles in the sample cup were put into the crucible of the moisture content measuring machine and dried at 120 °C for 7 h. Then, the mass m3 and m4 of Baijiu yeast granules before and after drying were read out, and then the CPI mass was deducted proportionally and recorded. Finally, calculation of the dry-basis moisture content of the sample particles was performed according to the following Equation (4):
Moisture   content   ( % ) = m 3 m 4 m 4 × 100 %
The method for measuring the particle size of dry particles is as follows: first, the sampling hopper was inserted into the sampling port of the dry powder sampler of winner 3001 dry powder laser particle sizer. Then, the dry particles mixed in the feed tank were filled up through the feed hopper, and after switching to 1.0–300 μm, the operation button was clicked. Finally, the particle size distribution and average value of dried yeast particles were automatically measured by the laser particle size analyzer, and a report was automatically generated. The average particle diameter of the dried particle sample in the first sampling cup was the input parameter—initial diameter of yeast particles (D0).
The measurement method of moisture state is as follows: first, the scanning temperature of DSC was set to 0–160 °C, and the temperature rise rate is 5 °C/min. Then, 0.1 g of dry particle sample was put into an aluminum crucible for sealing and then was placed into a DSC scanner for scanning. Finally, with the blank aluminum crucible as the blank control, the bound water content and free water content of the dried particles were automatically determined. The content value of the combined water in the water state, i.e., the water content at the transfer point, was taken as the input parameter.
The moisture content, moisture state, and particle size of the dried particles could be determined directly by a moisture meter, DSC, and laser granulometry, respectively, according to standard methods. The rehydration survival rate (Qv) could be calculated according to the following Equation (3) after counting the number of colorless bacteria (Nus) and the total number of bacteria (NT) in the four top corner counting areas of the blood cell counting board by methylene blue staining + blood cell counting method [18]. The equation is as follows:
Q v = N u s N T t h r e e f i e l d s × 100 %
The 16 groups of experimental parameters formed by the combination of 3 test variables and their representative patterns when mapping is shown in Table 1 below (The intersection diagram of the three groups of data in Table 1 is a representative icon of the results obtained after the experiment of the experimental parameters composed of the three groups of data in the subsequent mapping). Three repetitions shall be made for each group of experiments. a total of 48 groups of spray-drying experiments shall be completed, and 1392 test data shall be measured and recorded.

3.3. Test Data Analysis

The spray-drying experiment of Baijiu Yeast according to the above experimental design was carried out, and the temperature at the sampling point was recorded. Then, the moisture content, particle size, and the survival rate of yeast in rehydrated samples were determined respectively. Finally, the analysis with Origin is shown in Figure 5, Figure 6, Figure 7 and Figure 8.
As shown in Figure 5 above, with increasing sampling height, the measured sampling point temperatures at the cup mouth of the sampling cup all showed a decreasing trend of speed reduction. The rate of temperature drops at the sampling site decreased obviously after sampling height of 15 cm and stabilized near the outlet temperature at sampling site temperature when the sampling height was 45 cm.
As shown in Figure 7, when the sampling height was less than 45 cm, the dry-based water content of yeasts all decreased rapidly with increasing sampling height. After sampling height was greater than 45 cm, under each drying parameter, with the increase of sampling height, the dry-basis water content of yeast, however, showed different variation rules as follows. ① In curve 2 of group-c and all curves of group-d, the dry-basis moisture content of yeast decreased with the increase of sampling height. ② In all the curves except curve 2. In group-b and group-c, the moisture content of dry basis at the sampling height of 75 cm is lower than that at 90 cm. ③ In all curves of group-a, the moisture content of dry basis is lower than 75 cm at the sampling height of 60 cm and then lower than 90 cm.
Based on the analysis of Figure 5 and Figure 7, it could be concluded that in the constant rate drying stage, yeast particles were rapidly dehydrated in a high-temperature environment, which makes the temperature in this stage decrease rapidly, and the sampling height was 45 cm, which is close to the outlet temperature of spray-drying. In the subsequent stage, yeast particles were dried at a low rate in a medium- and low- temperature environment close to the outlet temperature.
It was assumed that the time for the droplet to move from one sampling point to another is equal, which is 2 s. If quadratic interpolation was performed at the time point without data, the change law of moisture content and rehydration survival rate of yeast dry particles with time could be obtained, as shown in Figure 9 and Figure 10.
As shown in Figure 9 above, the moisture content of dry particles decreases with the increase of drying time at a constant rate first and then at a reduced rate.
Combined with the analysis of Figure 7 and Figure 9, yeast particles was not entirely top-down moving. Firstly, in group-c curve 2 and all curves in group-d, yeast particles move from top to bottom without reflux. Secondly, in all the curves except curve 2 in group-b and group-c, yeast particles refluxed from the sampling height of 90 cm to 75 cm, and then decreased along the tower wall. Finally, in all curves of group-a, the yeast particles refluxed from the sampling height of 90 cm to 60 cm, and then decreased along the tower wall. The motion track is shown in Figure 11.
By comparing Figure 8 and Figure 10, it could be seen that in the high-speed and constant-speed drying stage, the surface moisture of bacteria particles was sufficient to achieve gas-liquid two-phase equilibrium, and the contact temperature of yeast is the wet bulb temperature, lower than the hot air temperature. At this stage, the main cause of cell inactivation was dehydration damage caused by rapid dehydration. The survival rate of yeast is high and was less affected by drying process parameters. In the stage of reduced rate drying, the moisture on the particle surface could not meet the requirements of gas–liquid two-phase equilibrium, and the moisture inside the particle then migrates to the surface through the pore structure. With the increase of migration resistance, the particle drying rate also decreases. At this stage, when the bacteria were exposed to hot air, the main causes of cell inactivation are heat damage and dehydration damage. The survival rate of yeast was greatly affected by the outlet temperature.

4. Mathematical Modeling of the Spray-Drying Process

Based on the test data measured by the new sampling system after sampling, the spray-drying process of Baijiu yeast was modelled to establish a mathematical model that could accurately predict the changes of particle size, water content and rehydration survival rate of Baijiu Yeast. Finally, it was generalized to apply the prediction of all-microbial drying, which is the core purpose of this study.

4.1. Common Data Values and Assumptions

To establish the mathematical prediction model based on the experimental parameters, reasonable assumptions must be made in combination with the experimental data, and then the physical model of the spray-drying process and yeast particles should be simplified. Therefore, based on the experimental results, five hypotheses were made as follows.
(1) According to the analysis of water content and rehydration survival rate of yeast particles, the water content in yeast particles could be divided into free water (Class E), intracellular bound water (Class D), and residual water (Class r), as shown in Figure 12 below. The figure on the right is the schematic diagram of yeast spray-drying particles, and the figure on the left is the abstract model diagram of drying particles.
Since cell membrane is the most typical semi-permeable membrane, which has no effect on the entry and exit of small molecules such as water molecules, in this study, the water in the outer protective agent of cell membrane and the flowing water that cells do not combine with cell structures were uniformly classified as Class E free water [19]. According to the pore network theory, it was assumed that the pore is uniformly distributed in a microscopic pore network (pore throat) with a diameter of 0.1 μ m , and its loss was mainly obstructed by capillary forces in porous media and saturated vapor pressure on the particle surface. The presence of Class E water maintained the gas–liquid two-phase balance on the surface of dried particles, and kept the bacteria in the “wet bulb temperature” at the early stage of drying. The yeast was almost not damaged by heat (survival rate >95%).
Class D water is the bound water in the bacterial body (excluding Class r bound water remaining after drying), which binds to each cellular substructure through hydrogen bounding or intermolecular forces, etc., and acts to maintain the normal function of each cellular substructure. Its dispersion and loss were mainly impeded by the adsorption of cellular substructures, capillary forces in the porous medium and saturated vapor pressure at the particle surface. The loss of Class D water leads to damage to cellular substructures, so the loss of Class D water is the main cause of dry inactivation of yeasts. It was assumed that the two types of water are lost in the order of E and D, and the remaining water (Class r) after drying is all Class D water. Because   t m d and t c are much smaller than t d r , it was considered that the mass transfer inside and outside the particles is in a quasi-steady state.
(2) Since the radius of the grid channel of the protective agent outside the yeast body is 10 μ m , the capillary effect has negligible effect on the saturation pressure of the particle gas phase.
(3) Since t e is much smaller than t d e , gas-phase heat exchange around the particle could be assumed to be quasi steady state. Since t d is much smaller than t d r , it could be assumed that the temperature inside the particles is uniformly distributed.
(4) Since t d r and t h have the same order of magnitude, it could be assumed that the heating of the particles continues throughout the drying cycle and that the particle surface temperature is equal to the air temperature.
Based on the above assumptions, the following equations could be used to describe the drying of Baijiu yeast particles in spray-drying, and to predict the changes of moisture content and rehydration survival rate of yeast particles with drying time. Furthermore, the new sampling scheme could be used to sample and measure the spray-drying process of other microorganisms, to change the parameters of the mathematical model and extend it to the drying of all microorganisms.
The meanings and values of known or unknown variables used in the modeling process have shown in Table 2 below [20].

4.2. Mathematical Modeling

4.2.1. Initial Conditions

The initial conditions were expressed as follows.
D ( t = 0 ) = D 0
T g ( t = 0 ) = T g 0
u = u 0
X ( t = 0 ) = X 0
X e ( t = 0 ) = X e 0
X d ( t = 0 ) = X d 0
In the Equations (6)–(11), D —diameter of yeast particles, mm, which could be measured by a laser particle sizer; T —yeast particle surface (hot air) temperature, °C; u —wind speed; X —water content, %; Each letter in the lower corner corresponds to each type of water, and 0 represents the initial value; X d 0 —initial dry-basis content of bound water, %, determined by DSC scanning.
After conversion, Equation (6) corresponds to the spray flow (Particle size), Equation (7) corresponds to the hot air temperature, Equation (8) corresponds to the wind speed, and Equations (9)–(11) correspond to the moisture content at the transfer point. The initial conditions of the model were found by entering the equations one at a time.

4.2.2. Particle Size Modeling

According to the experimental results and pore theory, it could be considered that the diameter of yeast particles is a linear function of the water content of the particles, and it is isotropic. The real-time diameter D of the dried particles could be expressed as a function of the real-time moisture content X 0 with the coefficients D 0 , X 0 , and the shrinkage coefficient κ as follows:
D = D 0 ( 1 κ X 0 X X 0 ) = κ D 0 X X 0 + D 0 κ D 0
The shrinkage coefficient κ of particles in Equation (12) could be estimated as κ = 0.46 based on a large number of experimental results in Figure 6.

4.2.3. Mass Transfer Modeling

The mass balance of water could be expressed as:
X = X e + X d + X r
d X e d t = J e
d X d d t = J d
d X r d t = 0
In Equations (13)–(16), J X —Class X water loss rate, where X is representative of three categories of water, kg / s ; t —time, s .

Class E Water

Class E water is in a pore channel network of pore throat diameter 0.1 μ m made up of cytoplasm and cell membrane. Moreover, the dissipation rate J e of Class E water could be characterized by vapor diffusion in the particle dry zone (internal mass transfer coefficient k i ) and combined vapor diffusion convection (external mass transfer coefficient k e ) in the boundary layer around the particle as follows.
J e = a M w k e k i k e + k i ( C s a t C )
In Equation (17), a —specific surface area, m 2 / kg ; M w —relative molecular mass of water, kg / mol ; k e —external mass transfer coefficient of the particles; k i —internal mass transfer coefficient of the particle; C s a t —the saturated concentration of water in air, mol / m 3 ; C —the concentration of water vapor in the air far from the particle, mol / m 3 ; As available by the ideal gas assumption, C s a t and C could be expressed as follows [21]:
C s a t = P s a t R g T g
C = Y M w P t o t R g T g ( M a + Y M w )
In Equations (18) and (19), T g —yeast particle surface (hot air) temperature, °C, which could be measured by temperature sensor in the experiment; T —internal temperature of yeast particles, °C; P t o t —total pressure in drying tower, Pa; P s a t —water saturation pressure, Pa. The relationship between P s a t and the surface (hot air) temperature of yeast particles could be expressed by Clausius Clapeyron law as follows:
P s a t = P s a t r e f exp ( L k M w R g T r e f T g T g T r e f )
The external mass transfer coefficient ( k e ) in Equation (17) could be expressed as follows by the relation to the Sherwood’s number ( s h ) of a fluid flowing around a particle [19]:
s h = k e D δ i = 2 + 0.6 Re 0.5 S c 1 / 3
In Equation (21), Reynolds number Re and Schmidt’s number S c were defined as follows [22]:
Re = u D ν
S c = ν δ i
The expressions for the external mass transfer coefficient ( k e ) available by carrying Equations (22) and (23) into Equation (21) are as follows:
k e = 2 δ i + 0.6 u 1 / 2 δ i 2 / 3 D 1 / 2 v 1 / 6
The intra-particle mass transfer coefficient ( k i ) in Equation (17), under the quasi steady state assumption, could be expressed in terms of the mass balance of water vapor inside the drying zone of the particles as follows:
d d t ( r d C d r ) = 0
In Equation (25), the distance of the r—Class e water from the center of the particle, mm; C —concentration of water vapor in the voids of particles, mol / m 3 ; Its boundary conditions are:
C r = D 2 = C s u r f
C r = D h 2 = C s a t
In Equations (26) and (27): C s u r f —concentration of water vapor in the contact air at the outer surface of the particle, mol / m 3 ; C s a t —the saturated concentration of water in the air, mol / m 3 .
In Equation (25), d C / d r could be expressed by the mass flux j away from the yeast particle water vapor as follows:
d C d r r = 1 2 D h = j τ ε δ i
Take Equations (26)–(28) into Equation (25):
ε δ i τ d C d r r = 0.5 D h = ε δ i τ 2 D 1 ln ( D ) ln ( D h ) ( C s a t T C s u r f ) = k i ( C s a t T C s u r f )
The expression of internal mass transfer coefficient ( k i ) could be obtained by sorting out Equation (29) as follows:
k i = ε δ i τ 2 D 1 ln D D h
The moisture content and the diameter of the spheres could be considered to be directly proportional, assuming that all water is uniformly distributed and lost layer-by-layer within the spherical region where it was located (Figure 12), so the proportional relationship could be written as when starting drying (t = 0) and during drying (t = t):
4 π 3 D 2 3 = α X e 0
4 π 3 D h 2 3 = α X e
Dividing Equation (32) by Equation (31) was obtained:
D D h = ( X e 0 X e ) 1 / 3

Class D Water

J d could be expressed by the first-order dynamic model [23],
J d = k d ( X d X d 0 X e X e 0 )
In Equation (34), k d —the kinetic constant related to the particle surface temperature T g , s 1 .
In summary, the dispersion rates and the water content models for the two types of water available after finishing E and D are as follows:
J e = 2 ε δ i a M w ( 2 δ i + 0.6 u 1 / 2 δ i 2 / 3 D 1 / 2 v 1 / 6 ) ( ( 2 δ i + 0.6 u 1 / 2 δ i 2 / 3 D 1 / 2 v 1 / 6 ) τ D ln ( ( X e 0 X e ) 1 / 3 ) 2 ε δ i P s a t r e f exp ( L k M w R g T r e f T g T g T r e f ) M a + Y M w ( P s a t r e f exp ( L k M w R g T r e f T g T g T r e f ) P t o t ) R g T g ( M a + Y M w ) = d X e d t
J d = k d ( X d X d 0 X e X e 0 ) = k d X d = d X d d t

4.2.4. Heat Transfer Modeling

The macroscopic thermal equilibrium equation for yeast particles was as follows [24]:
a h ( T g T ) = ( J e + J d ) L k + ( 1 + X ) c p , w l d T d t
The relaxation time t h * for the particle heating was defined as follows [25]:
t h * = ( a h ( 1 + X ) c p , w l ) 1
The expression of the heat transfer rate was obtained by bringing Equation (38) into Equation (37), as shown below:
d T d t = 1 t h * ( T g T ( J e + J d ) L k a h )
In Equations (37)–(39): L k —vaporization latent heat of water, j / k g ; h —heat exchange coefficient, W / k m 2 ;
In Equations (38) and (39), the heat exchange coefficient h between the particle and the surrounding air, could be expressed by the association formula of the Nussle number (Nu) of the fluid flow around the particle as follows [26]:
N u = h D λ a = 0.3 + 0.62 Re 0.5 P r 1 / 3 ( 1 + ( 0.4 Pr ) 2 / 3 ) 1 / 4
In Equation (40): Re —Reynolds number; Pr —Prandtl number. The definition formula is:
Pr = μ c p , a λ a
The heat transfer rate calculation model could be obtained by successively bringing Equations (38)–(41) into Equation (37):
d T d t = [ 0.3 a λ a ( 1 + X ) c p , w l + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 ( 1 + X ) c p , w l ] [ T g T ( J e + J d ) L k D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 0.3 a λ a D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 ]
Since the water in yeast particles was divided into two types, E and D, which are lost in turn, Equation (42) should be divided into two parts as follows:
d T d t = [ 0.3 a λ a ( 1 + X e ) c p , w l + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 ( 1 + X e ) c p , w l ] [ T g T J e L k D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 0.3 a λ a D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 ]
d T d t = [ 0.3 a λ a ( 1 + X d ) c p , w l + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 ( 1 + X d ) c p , w l ] [ T g T J d L k D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 0.3 a λ a D ( 1 + ( 0.4 λ a μ c p , a ) 2 / 3 ) 1 / 4 + 0.62 a Re 0.5 μ 1 / 3 c p , a 1 / 3 λ a 2 / 3 ]

4.2.5. Rehydration Survival Rate Modeling

Thermal inactivation of microbial cells is usually expressed using the first-order kinetic equation as follows [27]:
d Q v d t = k r Q v
In Equation (45), k r —deactivation rate constant. Since dehydration damage and thermal damage are the main factors causing yeast inactivation in spray drying, available from the above heat transfer modeling, the change of temperature is also closely related to the moisture content within the yeast particles. So in this study, it was assumed that thermal inactivation is the most dominant inactivation way of yeast, and k r could be correlated with the internal temperature T of the particles using the Arrhenius equation as follows [28]:
k r = k 0 e E d R T
In Equation (46), k 0 —pre exponential factor, E d —activation energy of inactivated cells, R—gas constant.
Carrying Equation (46) into Equation (45) was available:
d Q v d t = Q v k 0 e E d R T

4.3. Model Description

The initial parameters of spray drying were input according to Equations (4)–(9), and then the values of the value parameters in Table 2 were brought into Equations (10), (33), (34), (41), (42) and (45). Then, according to the mathematical model operation logic shown in Figure 13, the function curve of moisture content and rehydration survival rate of dried particles in spray drying process with time could be obtained.

4.4. Determination of Model Parameters

As shown in Table 2, there are also three modeling parameters ε , X d 0 , and k d that will change with the change of the spray-drying process and microbial species. Therefore, this study plans to take Baijiu Yeast as an example to determine the changing model parameters, to guide the determination of other microbial spray-drying modeling parameters.
ε is the porosity before the drying of yeast pellets, which was determined by a large number of trials in a modelling study by debase on the hot-air drying of yeast pellets ε = 0.58 [18].
As shown in Figure 12, it was assumed that the three types of water E, D, and r are uniformly distributed in the circle from the outside to the inside in that order, and are lost in turn. Based on the ease of water dispersal and the degree of damage to yeast in each category, X e 0 = 1.8–2.2 (d.b), X d 0 = 0.3–0.4 (d.b), and X r = 0.05–0.11 (d.b) could be tentatively determined, which is the yeast particles dry-based water content at the water conversion point of Class E and Class D and was about 0.3–0.4 (d.b) [18].
The moisture content (moisture state) of the bound water on a dry basis, also known as the moisture content at the transition point, could be accurately measured by differential scanning calorimeter (DSC) in the experiment and used as the input parameter of the model. The value of bound water dry-basis moisture content X d 0 of yeast dry granules obtained under various spray-drying parameters is shown in Table 3 below, and X d 0 = 0.35 (d.b) was taken after comparative analysis.
Because b < 0.4 (D.B), k d could be determined by linear regression ln ( X X r ) , and there is X = X d + X r , that is, the change of X d with time could be expressed as follows by Equation (34) [29]:
X d = X d 0 exp ( k d t )
Take the logarithm of Equation (48) to obtain:
k d = ln ( X d ) ln ( X d 0 ) t = ln ( X X r ) ln ( X d 0 ) t

5. Results and Analysis

Firstly, based on four kinds of inlet hot air temperatures, a group of typical process parameters were selected for spray-drying of Baijiu Yeast. Then, the test data were brought into the above established mathematical model to obtain the corresponding prediction function, then Origin was used to make a diagram. Finally, the predicted function curve of the mathematical model of the Particle size, moisture content, and rehydration survival rate of yeast particles with time has been compared with the experimental results, as shown in Figure 14, Figure 15 and Figure 16 respectively.
As shown in Figure 14, the coincidence rate of the function curve of dry particle size predicted by the mathematical model with its dry-basis moisture content and the experimental data is as high as 96%, which proves that the mathematical model could better predict the historical change of dry particles during spray drying.
According to the analysis in Figure 15 and Figure 16, the prediction curves and experimental values of the heat and mass transfer mathematical model of porous media based on real-time sampling data for the moisture content and rehydration survival rate of yeast particles in the spray drying process were obtained when the inlet temperature was between 95 °C and 125 °C. When the inlet hot air temperature is 115 °C—the hot air flow is 60 m 3 / h —the material flow is 600 mL / h , the maximum errors X max = 0.03 ( d . b ) and Q v max = 6.43 % were taken when the sampling time was 12 s. In the dry DC circulation area, the moisture content and rehydration survival rates predicted by the model almost coincide with the experimental values. After drying, the predicted water content and rehydration survival rate were lower than the experimental values. However, they were close to overlapping at the end.
The predicted and experimental values of the water bearing rate and rehydration survival rate of yeast particles during the drying process show that the predicted values in the DC circulation area almost coincide with the end point, while the predicted values in the recirculation circulation area are less than the experimental values. The main reason is that in the reflux circulation area, the sampling cup collects both reflux dry particles and falling dry particles, resulting in the moisture content and rehydration survival rate of dry particles measured in the test being higher than the predicted value. Apart from the recirculation area, the maximum error values of the model for predicting the moisture content and rehydration survival rate of yeast dry particles during spray-drying are X max 1 = 0.027 ( d . b ) and Q v max 1 = 1.06 % , respectively, which are less than 5%. It was proven that the mathematical model of heat and mass transfer in porous media-based on real-time sampling data could accurately predict the moisture content and rehydration survival rate of microbial particles in the process of spray-drying, and provide guidance for revealing the changes of microbial particles in the process of spray-drying. At the same time, it was proven again that the inference of particle trajectory in the process of spray drying was correct.

6. Conclusions

In this study, a sampling system for real-time sampling, rehydration, and non-destructive storage of microbial dry particles in the spray-drying process (which has been submitted for a Chinese invention patent) was introduced. With the cooperation of later detection, the system could reveal the changes of microorganisms and their particles in the high-speed, complex, and invisible spray-drying process, which lays a foundation for later mathematical modeling to predict and explore the microscopic damage mechanism of microorganisms in the spray drying process. Based on the measured data, the mathematical model of Baijiu yeast granule drying was established according to the law of heat and mass transfer and the theory of porous media. Comparing the predicted value of the mathematical model with the experimental results, the maximum error values of moisture content and rehydration survival rate are only X max 1 = 0.027 ( d . b ) and Q v max 1 = 1.06 % , respectively, both are less than 5%, which proves that the model could accurately predict the changes of moisture content and rehydration survival rate of microbial dried particles with drying time in the process of spray-drying.
Based on the analysis of the experimental data and the predicted values of the model, it was concluded that the spray-drying process could be divided into DC circulation area and reflux circulation area. In the DC circulation area, the yeast particles have enough kinetic energy to overcome the resistance of the laterally flowing air in the vortex area, move along the Original spray path, and dehydrate rapidly. In the recirculation area, the kinetic energy of yeast particles was weakened, which is not enough to overcome the resistance of the lateral air flow in the vortex area, but was pulled up along the wall direction and further reduced for dehydration. The moisture content Xz = 0.3–0.4 (d.b) of yeast particles transferred from the DC circulation area to the reflux area. After entering the reflux circulation area, the yeast begins to inactivate rapidly.
Due to the lack of yeast activity detection equipment, it is impossible to quantitatively detect the activity of yeast in dry granules. Therefore, it is impossible to conduct mathematical modeling and prediction on the change of yeast activity with time in the drying process and further verify the correctness of the mathematical model. In follow-up studies, the author will seek better experimental conditions and conduct in-depth research.

Author Contributions

Conception and writing, F.-k.X.; Review and editing, Y.-j.Y.; Language polish, Y.-y.X.; Supervision and data curation, J.-y.L.; Review, Z.Z.; Language polish, L.-b.T. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China [grant number 51876109]; the National Key R&D Program of China [grant number 2017YFD0400902-1]; and the Key Project of International Science and Technology Cooperation Program for Shaanxi Province [grant number 2020KWZ-015]; Scientific research plan project of youth innovation team of Shaanxi Provincial Department of Education [grant number 22JP012].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

This is a research article, researched by the author himself, without conflict of interest with any individual or collective, hereby declares.

References

  1. Tay, J.B.J.; Chua, X.; Ang, C.; Subramanian, G.S.; Tan, S.Y.; Lin, E.M.J.; Wu, W.-Y.; Goh, K.K.T.; Lim, K. Effects of Spray-Drying Inlet Temperature on the Production of High-Quality Native Rice Starch. Processes 2021, 9, 1557. [Google Scholar] [CrossRef]
  2. Toro-Sierra, J.; Schumann, J.; Kulozik, U. Impact of spray-drying conditions on the particle size of microparticulated whey protein fractions. Dairy Sci. Technol. 2013, 93, 487–503. [Google Scholar] [CrossRef]
  3. Koopmann, I.K.; Möller, S.; Elle, C.; Hindersin, S.; Kramer, A.; Labes, A. Optimization of Astaxanthin Recovery in the Downstream Process of Haematococcus pluvialis. Foods 2022, 11, 1352. [Google Scholar] [CrossRef]
  4. Al Zaitone, B.; Al-Zahrani, A.; Ahmed, O.; Saeed, U.; Taimoor, A.A. Spray Drying of PEG6000 Suspension: Reaction Engineering Approach (REA) Modeling of Single Droplet Drying Kinetics. Processes 2022, 10, 1365. [Google Scholar] [CrossRef]
  5. Teixeira, P.C.; Castro, M.H.; Malcata, F.X.; Kirby, R.M. Survival of Lactobacillus delbrueckii ssp. bulgaricus Following Spray-Drying. J. Dairy Sci. 1995, 78, 1025–1031. [Google Scholar] [CrossRef]
  6. Gong, P.; Sun, J.; Lin, K.; Di, W.; Zhang, L.; Han, X. Changes process in the cellular structures and constituents of Lactobacillus bulgaricus sp1.1 during spray drying. LWT—Food Sci. Technol. 2018, 102, 30–36. [Google Scholar] [CrossRef]
  7. Akpinar, E.K. Determination of suitable thin layer drying curve model for some vegetables and fruits. J. Food Eng. 2006, 73, 75–84. [Google Scholar] [CrossRef]
  8. Akpinar, S. Treatment of Restless Legs Syndrome with Levodopa Plus Benserazide. Arch. Neurol. 1982, 39, 739. [Google Scholar] [CrossRef]
  9. Spreutels, L.; Debaste, F.; Legros, R.; Haut, B. Experimental characterization and modeling of Baker’s yeast pellet drying. Food Res. Int. 2013, 52, 275–287. [Google Scholar] [CrossRef]
  10. Huang, S.; Vignolles, M.; Chen, X.D.; Loir, Y.; Jan, G.; Schuck, P.; Jeantet, R. Spray drying of probiotics and other food-grade bacteria: A review. Trends Food Sci. Technol. 2017, 63, 1–17. [Google Scholar] [CrossRef]
  11. Peighambardoust, S.H.; Tafti, A.G.; Hesaria, J. Application of spray drying for preservation of lactic acid starter cultures: A review. Trends Food Sci. Technol. 2011, 22, 215–224. [Google Scholar] [CrossRef]
  12. Park, B.J.; Kang, H.S.; Choi, H.J. Magnetic carbonyl iron nanoparticle based magnetorheological suspension and its characteristics. Mater. Lett. 2009, 63, 1350–1352. [Google Scholar] [CrossRef]
  13. Forrest, A.R. Interactive Interpolation and Approximation by Bezier Polynomials. Comput. Aided Des. 1990, 22, 527–537. [Google Scholar] [CrossRef]
  14. Gordon, W.J.; Riesenfeld, R.F. Bernstein-Bezier Methods for the Computer-Aided Design of Free-Form Curves and Surfaces. J. Assoc. Comput. Mach. 1974, 21, 293–310. [Google Scholar] [CrossRef]
  15. Bushnaq, H.; Krishnamoorthy, R.; Abu-Zahra, M.; Hasan, S.W.; Taher, H.; Alomar, S.Y.; Ahmad, N.; Banat, F. Supercritical Technology-Based Date Sugar Powder Production: Process Modeling and Simulation. Processes 2022, 10, 257. [Google Scholar] [CrossRef]
  16. Gong, P.; Di, W.; Yi, H.; Sun, J.; Zhang, L.; Han, X. Improved viability of spray-dried Lactobacillus bulgaricus sp1.1 embedded in acidic-basic proteins treated with transglutaminase. Food Chem. 2019, 281, 204–212. [Google Scholar] [CrossRef]
  17. Gong, P.; Zhang, L.; Zhang, J.; Han, X.; Yi, H.; Di, W.; Chen, S. Sampling method for studying the activity of lactic acid bacteria during spray drying. Dry. Technol. 2018, 36, 1236–1242. [Google Scholar] [CrossRef]
  18. Spreutels, L.; Haut, B.; Chaouki, J.; Bertrand, F.; Legros, R. Conical spouted bed drying of Baker’s yeast: Experimentation and multi-modeling. Food Res. Int. 2014, 62, 137–150. [Google Scholar] [CrossRef]
  19. Türker, M.; Kanarya, A.; Yüzgeç, U.; Kapucu, H.; Şenalp, Z. Drying of baker’s yeast in batch fluidized bed. Chem. Eng. Processing Process Intensif. 2006, 45, 1019–1028. [Google Scholar] [CrossRef]
  20. Pauli, M.; Kayser, T.; Adamiuk, G.; Wiesbeck, W. Modeling of Mutual Coupling between Electromagnetic and Thermal Fields in Microwave Heating. In Proceedings of the 2007 IEEE/MTT-S International Microwave Symposium, Honolulu, HI, USA, 3–8 June 2007. [Google Scholar]
  21. Martínez-Sánchez, A.; López-Cañavate, M.; Guirao-Martínez, J.; Roca, M.; Aguayo, E. Aloe vera Flowers, a Byproduct with Great Potential and Wide Application, Depending on Maturity Stage. Foods 2020, 9, 1542. [Google Scholar] [CrossRef]
  22. Koumbogle, K.; Gitzhofer, F.; Abatzoglou, N. Moisture Transport Coefficients Determination on a Model Pharmaceutical Tablet. Processes 2022, 10, 254. [Google Scholar] [CrossRef]
  23. Ursu, M.S.; Aprodu, I.; Milea, Ș.A.; Enachi, E.; Râpeanu, G.; Bahrim, G.E.; Stănciuc, N. Thermal Degradation Kinetics of Anthocyanins Extracted from Purple Maize Flour Extract and the Effect of Heating on Selected Biological Functionality. Foods 2020, 9, 1593. [Google Scholar] [CrossRef]
  24. Grasa, G.; Abanades, J. Narrow fluidised beds arranged to exchange heat between a combustion chamber and a CO2 sorbent regenerator. Chem. Eng. Sci. 2007, 62, 619–626. [Google Scholar] [CrossRef]
  25. Olivares-Quiroz, L.; García-Colín, L.S. Evidence of α fluctuations in myoglobin’s denaturation in the high temperature region: Average relaxation time from an Adam–Gibbs perspective. Biophys. Chem. 2009, 144, 123–129. [Google Scholar] [CrossRef]
  26. Parapanov, R.; Nussle, S.; Crausaz, M.; Senn, A.; Hausser, J.; Vogel, P. Testis size, sperm characteristics and testosterone concentrations in four species of shrews (Mammalia, Soricidae). Anim. Reprod. Sci. 2009, 114, 269–278. [Google Scholar] [CrossRef]
  27. Li, X.; Lin, S.X.Q.; Chen, X.D.; Chen, L.; Pearce, D. Inactivation Kinetics of Probiotic Bacteria during the Drying of Single Milk Droplets. Dry. Technol. 2006, 24, 695–701. [Google Scholar] [CrossRef]
  28. Huang, H.; Brooks, M.S.-L.; Huang, H.-J.; Chen, X.D. Inactivation Kinetics of Yeast Cells during Infrared Drying. Dry. Technol. 2009, 27, 1060–1068. [Google Scholar] [CrossRef]
  29. Schwaab, M.; Pinto, J.C. Optimum reference temperature for reparameterization of the Arrhenius equation. Part 1: Problems involving one kinetic constant. Chem. Eng. Sci. 2007, 62, 2750–2764. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of Bezier cubic curve.
Figure 1. Schematic diagram of Bezier cubic curve.
Processes 10 01789 g001
Figure 2. Comparison of pressure field and wind speed field between straight cylinder sampling cup and curved sampling cup.
Figure 2. Comparison of pressure field and wind speed field between straight cylinder sampling cup and curved sampling cup.
Processes 10 01789 g002
Figure 3. Diagram of a new double surface sampling cup structure. (1) Heat and magnetic insulation tape; (2) Connecting magnet; (3) Outer cup; (4) Frozen layer; (5) Wireless temperature and humidity sensors; (6) Upper part of inner cup; (7) Lower part of inner cup; (8) Powerful magnets; (9) Saline.
Figure 3. Diagram of a new double surface sampling cup structure. (1) Heat and magnetic insulation tape; (2) Connecting magnet; (3) Outer cup; (4) Frozen layer; (5) Wireless temperature and humidity sensors; (6) Upper part of inner cup; (7) Lower part of inner cup; (8) Powerful magnets; (9) Saline.
Processes 10 01789 g003
Figure 4. Schematic of the parallel flow drying column structure. (1) Feed pump, (2) & (6) Filter, (3) & (9) Fan feeders, (4) Steam heaters, (5) Electric heaters, (7) Nebulizers, (8) Drying column, (10) Air pump, (11) Gas sweep device, (12) Cyclone separator.
Figure 4. Schematic of the parallel flow drying column structure. (1) Feed pump, (2) & (6) Filter, (3) & (9) Fan feeders, (4) Steam heaters, (5) Electric heaters, (7) Nebulizers, (8) Drying column, (10) Air pump, (11) Gas sweep device, (12) Cyclone separator.
Processes 10 01789 g004
Figure 5. Sampling point temperature of each sampling cup under different test parameters.
Figure 5. Sampling point temperature of each sampling cup under different test parameters.
Processes 10 01789 g005
Figure 6. Variation of yeast particle size with moisture content.
Figure 6. Variation of yeast particle size with moisture content.
Processes 10 01789 g006
Figure 7. Moisture of yeast in each sampling cup under different test parameters.
Figure 7. Moisture of yeast in each sampling cup under different test parameters.
Processes 10 01789 g007
Figure 8. Rehydration survival of yeast in each sampling cup under different test parameters.
Figure 8. Rehydration survival of yeast in each sampling cup under different test parameters.
Processes 10 01789 g008
Figure 9. Moisture content of yeast under different test parameters with time.
Figure 9. Moisture content of yeast under different test parameters with time.
Processes 10 01789 g009
Figure 10. Rehydration survival rate of yeast under different test parameters with time.
Figure 10. Rehydration survival rate of yeast under different test parameters with time.
Processes 10 01789 g010
Figure 11. Schematic diagram of yeast particle trajectory under three groups of parameters.
Figure 11. Schematic diagram of yeast particle trajectory under three groups of parameters.
Processes 10 01789 g011
Figure 12. Schematic diagram of Baijiu yeast granule drying model. C . Water vapor concentration in the air far from the pellet, C s u r f . Water vapor concentration in the air in contact with the external pellet surface, C s a t . Saturation concentration of water in pellet, D . Real time diameter of yeast particles, D h . Diameter of bound water, D r . Diameter of the residual water, D 0 . Yeast particle initial diameter.
Figure 12. Schematic diagram of Baijiu yeast granule drying model. C . Water vapor concentration in the air far from the pellet, C s u r f . Water vapor concentration in the air in contact with the external pellet surface, C s a t . Saturation concentration of water in pellet, D . Real time diameter of yeast particles, D h . Diameter of bound water, D r . Diameter of the residual water, D 0 . Yeast particle initial diameter.
Processes 10 01789 g012
Figure 13. Liquor yeast diagram of mathematical model operation logic.
Figure 13. Liquor yeast diagram of mathematical model operation logic.
Processes 10 01789 g013
Figure 14. Comparison between model curve and experimental value of yeast particle size changing with moisture content.
Figure 14. Comparison between model curve and experimental value of yeast particle size changing with moisture content.
Processes 10 01789 g014
Figure 15. Model curve of yeast particle dry based moisture content as a function of time contrasting with experimental values.
Figure 15. Model curve of yeast particle dry based moisture content as a function of time contrasting with experimental values.
Processes 10 01789 g015
Figure 16. Comparison between model curve and experimental value of yeast granule rehydration survival rate with time.
Figure 16. Comparison between model curve and experimental value of yeast granule rehydration survival rate with time.
Processes 10 01789 g016
Table 1. Test variable group setting and representative pattern.
Table 1. Test variable group setting and representative pattern.
Inlet Hot Air Temperature (°C)95105115125
Hot air inflow (m3/h)60+
70×
Material inflow (mL/h)600700600700600700600700
Table 2. Meanings and values of modeling parameters.
Table 2. Meanings and values of modeling parameters.
VariableMeaningValue
ρ w l Density of liquid water (25 °C)1000 kg m−3
ρ a Density of water vapor (25 °C)1.2 kg m−3
M w Relative molecular weight of water1.8 × 10−3 kg mol−1
M a Molecular weight of air2.8 × 10−3 kg mol−1
R g Universal gas constant8.314 J (mol K)−1
P s a t r e f Reference water saturation pressure4243 Pa
T r e f Reference temperature25 °C
L k Latent heat of vaporization of water2436 Mj kg−1
δ i Diffusion coefficient of water vapor in air2.5 × 10−5 m2 s−1
ν Kinematic viscosity of air1.7 × 10−5 m2 s−1
μ Aerodynamic viscosity1.8 × 10−5 Pa s
c p , wl Specific heat capacity of liquid water4185 J (kg K)−1
c p , a Specific heat capacity of air1000 J (kg K)−1
λ a Air thermal conductivity2.62 × 10−2 W (m K)−1
D 0 Initial diameter of yeast spray particles 100 400   μ m (Determination by particle size)
DDiameter of yeast particles 100 300   μ m
P t o t Total pressure in drying tower101,325 Pa
Y Moisture content of air on a dry basis0.01 (d.b.)
u Hot wind speed6–8 m s−1
X 0 Initial moisture content of yeast granules2.2–2.5 (d.b.)
T g 0 Initial hot air temperature40–125 °C
T 0 Initial yeast spray particle temperature28 °C
t d r drying time14–22 s
t c Contact time between gas flow and particles t c 2 R   / u 0.001 s
  t m d Diffusion time of water vapor in particles   t m d R 2 / δ 1 0.01 s
t h d Thermal diffusion time inside yeast particles t h d R 2 / k w l 2 s
t h Heating time of yeast particles in gas flow t h ( 1 + X ) C p , w l / a h 14 22 s
τ Tortuosity of dry porous shell4.45
T g surface temperature of yeast particlesUnknown
ε Porosity of yeast granulesUnknown
k d Kinetic constantUnknown
Table 3. Dry-basis content of bound water (d.b).
Table 3. Dry-basis content of bound water (d.b).
Inlet Hot Air Temperature (°C)95105115125
Hot air inflow (m3/h)600.35 ± 0.0120.36 ± 0.010.34 ± 0.0170.35 ± 0.0090.35 ± 0.0130.35 ± 0.0100.36 ± 0.0040.34 ± 0.010
700.35 ± 0.0040.36 ± 0.0110.35 ± 0.0120.34 ± 0.0140.36 ± 0.0080.34 ± 0.0170.35 ± 0.0110.36 ± 0.007
Material inflow (mL/h)600700600700600700600700
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xiong, F.-k.; Yuan, Y.-j.; Xu, Y.-y.; Li, J.-y.; Zhao, Z.; Tan, L.-b. Modeling Study of a Microbial Spray-Drying Process Based on Real-Time Sampling. Processes 2022, 10, 1789. https://doi.org/10.3390/pr10091789

AMA Style

Xiong F-k, Yuan Y-j, Xu Y-y, Li J-y, Zhao Z, Tan L-b. Modeling Study of a Microbial Spray-Drying Process Based on Real-Time Sampling. Processes. 2022; 10(9):1789. https://doi.org/10.3390/pr10091789

Chicago/Turabian Style

Xiong, Feng-kui, Yue-jin Yuan, Ying-ying Xu, Jing-yu Li, Zhe Zhao, and Li-bin Tan. 2022. "Modeling Study of a Microbial Spray-Drying Process Based on Real-Time Sampling" Processes 10, no. 9: 1789. https://doi.org/10.3390/pr10091789

APA Style

Xiong, F. -k., Yuan, Y. -j., Xu, Y. -y., Li, J. -y., Zhao, Z., & Tan, L. -b. (2022). Modeling Study of a Microbial Spray-Drying Process Based on Real-Time Sampling. Processes, 10(9), 1789. https://doi.org/10.3390/pr10091789

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop