Supercritical Technology-Based Date Sugar Powder Production: Process Modeling and Simulation

: Date palm fruits ( Phoenix dactylifera ) contain high levels of fructose and glucose sugars. These natural sugar forms are healthy, nutritional and easily assimilate into human metabolism. The successful production of soluble date sugar powder from nutritious date fruits would result in a new food product that could replace the commercial reﬁned sugar. In this work, a novel process technology based on the supercritical extraction of sugar components from date pulp was modeled and simulated using Aspen Plus software. The process model consisted of three main steps that were individually simulated for their optimal working conditions as follows: (a) freeze-drying of the date pulp at − 42 ◦ C and 0.0001 bar; (b) supercritical extraction of the sugar components using a 6.77 wt.% water mixed CO 2 solvent system at a pressure of 308 bar, temperature of 65 ◦ C, and CO 2 ﬂow rate of 31,000 kg/h; and (c) spray-drying of the extract using 40 wt.% Gum Arabic as the carrier agent and air as drying medium at 150 ◦ C. The overall production yield of the process showed an extraction efﬁciency of 99.1% for the recovery of total reducing sugars from the date fruit. The solubility of the as-produced date sugar powder was improved by the process selectivity, elimination of insoluble ﬁber contents, and the addition of Gum Arabic. The solubility of the ﬁnal date sugar product was estimated as 0.89 g/g water. Results showed a maximum product ﬂow (glucose + fructose) of 73.70 kg/h was achieved at an operating pressure of 308 bar. and simulation studies were performed using Aspen Plus software with the sensitivity analysis to identify the optimum values of critical parameters for each stage. Freeze-drying studies revealed a minimum duty requirement of 13.23 kW at a temperature of − 42 ◦ C and pressure of 0.0001 bar to reduce the moisture content of the fruit pulp from 19.62% to 0.26%. The Sc-CO 2 assisted extraction resulted in a total sugar recovery of 99.1% with ﬂow rates of 34.14 kg/h and 39.13 kg/h for glucose and fructose sugars, respectively. Optimization studies showed that a pressure of 308 bar, temperature of 65 ◦ C, Sc-CO 2 ﬂow rate of 31,000 kg/h, and Sc-CO 2 : water ratio of 0.07 resulted in the maximum extraction efﬁciency. Further, spray-drying simulation of the date extract into a powdered product using Gum Arabic as the carrier agent showed a mean particle size of 50.4 µ m for the date sugar powder with a solubility of 0.89 g/g water (at 25 and 1 bar) and density of 0.0023 kg/m 3 . Sensitivity analysis revealed that the extraction pressure and solvent ﬂow rates were the signiﬁcant parameters for the extraction step while operating temperature had no salient effects on the recovery of the sugar. For the spray-drying process, drying air ﬂow rate and its inlet temperature had positive effects, while feed ﬂow rate negatively impacted moisture removal. Results showed that a mass ﬂow rate of 8000 kg/h and an inlet temperature of 150 ◦ C for the drying air was optimal to dry a feed stream mixture of date extract and Gum Arabic ﬂowing at 435 kg/h containing 40 wt.% of Gum Arabic. Thus, this work conﬁrmed the technical feasibility of IDSPP for the production of soluble date sugar powder and provided the necessary data with optimum operating conditions required for the process. Experimental validations and detailed characterizations of the product are further required to realize the outcomes of this study.


Introduction
Recently, the increasing awareness of the health effects of refined sugar has shifted consumers' focus towards natural and organic sugars. This has provided a tremendous scope for the global organic sugar market, which is expected to grow by 14% between 2019 and 2027 for a net worth of USD 4500 million by 2027 [1]. To meet the growing demands for organic sugar, a wide variety of plant sources and their products are being analyzed for the development of natural fruit-sugar products that benefit human health. Among them, date palm (Phoenix dactylifera) is a promising plant widely cultivated in the arid regions of the Persian Gulf and North Africa [2]. Date palm fruits contain high amounts of fructose and glucose sugars with very little sucrose content [3]. These natural sugar forms are simple, healthy, and assimilate easily into human metabolism. The fruit contains numerous health benefits, including nutritional, functional, and therapeutic values [4]. Date sugar is commercially produced in syrup form through a hot water extraction method. However, there are various drawbacks associated with this method, such as low extraction performance, high energy input, and degradation of thermolabile compounds [5].
Recovery of sugar components from the date fruit biomatrix is mainly hindered by the tough skin and gummy consistency of the fruit pulp. Moreover, the rigid cell membrane layers of the biomatrix retard the movement of solvent in and out of the plant cells. As a result, most of the valuable sugar components are lost with spent biomass in the conventional hot water extraction process, decreasing the extraction efficiency [6]. Thus, an innovative process technology is required to overcome these limitations and efficiently recover the sugar components from date fruits. Non-conventional techniques such as microwave-, ultrasonic-, supercritical-, enzymatic-, and magnetic-/electrical-assisted extraction are widely explored for the effective extraction of various bioactive compounds from different plants and their parts [7]. These advanced approaches are advantageous because they are less costly to operate, require a shorter treatment time, consume low amounts of energy, and do not produce any toxic byproducts.
Compared to other non-conventional extraction techniques, supercritical CO 2 (Sc-CO 2 ) extraction has gained increased popularity due to its sustainability, versatility, and absence of solvent contamination [8]. Due to the properties that can be readily modified with changes in pressure and temperature, supercritical fluids have a wide range of applications. Indeed, the properties of both liquid and gaseous matter are combined in the supercritical phase [9]. One of the notable benefits of supercritical solvent-based extraction is that the solvent can be fully separated from the desired solutes, thus ensuring product purity. Various studies have investigated the extraction of high-value components from fruits and vegetable tissues using supercritical fluids, specifically Sc-CO 2 [9][10][11][12]. However, usage of Sc-CO 2 for plant carbohydrates extraction is partly hampered by its low polarity, which limits the solubilization of the targeted compounds into Sc-CO 2 . To overcome this problem, high-density polar solvents are used as a co-solvent to enhance the solubility of polar sugars in Sc-CO 2 . Studies on carbohydrate solubility in Sc-CO 2 have shown that using polar co-solvents such as water significantly improves the solubility and separation of carbohydrates from mixtures at an optimal concentration of the co-solvent [13].
Further, a pre-drying step of the date pulp (prior to supercritical extraction) is essential for the maximum recovery of sugars and other nutritive analytes from the fruit matrix. Among the different pre-drying methods that could be applied to produce dehydrated food, freeze-drying is recognized as one of the best drying techniques to preserve food quality [14][15][16]. Freeze-drying, also known as lyophilization, is the process of removing water from the food by an initial freezing phase followed by a sublimation phase to remove the ice as vapor. The removal of water from food by sublimation protects the food product against the loss of essential components and damage caused by chemical reactions associated with the withdrawal or vaporization of liquid water [17].
A post-drying stage of the extract (obtained from the Sc-CO 2 extraction) to a solid powder form ensures product stability and longevity of the sugar. The surface area of the final date sugar product plays a vital role in determining its dissolution rate. Spray-drying is one of the most widely used solvent evaporation techniques in the food industry to obtain products with a large surface area [18,19]. However, the rubbery nature of the low molecular weight sugars in the extract (owing to their lower transition temperature) causes product instability and particle agglomeration during the drying process. Hence, a carrier agent is highly required in the post-drying step to overcome these limitations. Gum Arabic (GA) is a widely studied carrier agent due to its low viscosity and high solubility [18]. It is a natural prebiotic and is considered natural, edible, and a safe source of dietary fiber [20]. GA is a complex polysaccharide that benefits gut health, digestion, and cardiovascular health.
A careful evaluation of the scientific literature shows that freeze-drying and spraydrying are the best drying techniques for the pre-and post-drying stages of the main extraction [14,21]. Hence, the same methods were adopted in this study to develop a novel process technology based on Sc-CO 2 extraction for the enhanced recovery and production of soluble date sugar powder from the date fruit. Based on these literature findings, an Integrated Date Sugar Production Process (IDSPP) involving three main stages is conceptualized as follows: (i) Freeze-drying of date pulp, (ii) Sc-CO 2 -assisted extraction of date syrup with water as a co-solvent [8], and (iii) spray-drying of the extract for the production of soluble powder date sugar.
This study focuses on the modeling and simulation of the IDSPP process using Aspen Steady-State Modeler software to establish the feasibility of the process and determine the optimum working conditions for the process to maximize the sugar yield of date fruits. The simulation results would provide the necessary process data for the new technology and the recommended values for various critical process parameters for the experimental development and realization of the IDSPP process. Continuous-model simulation was performed for each stage of the process as a discrete block to determine the optimal values of the critical operating parameters for the respective stage of the IDSPP plant.

IDSPP Flowsheet
The process flowsheet for the as-proposed IDSPP technology is shown in Figure 1. Freeze-drying, supercritical extraction, and spray-drying are the major unit operations in the process. Operating pressure, freezing temperature, and drying temperature are the main parameters examined for the freeze-drying (FD) stage to minimize the moisture content in the fruit. Regarding the supercritical extraction (SCE) stage, variables such as operating pressure, temperature, and solvent:co-solvent ratio were considered to optimize the extraction efficiency. For the final spray-drying (SD) stage, the mass flow rate and temperature of air (drying medium), and feed flow rate were studied to understand their effects on the size and solubility of the final date sugar product.

Materials
The lack of experimental data for the IDSPP at its design stage creates greater reliability on numerical simulation to provide useful information for the process. Commercial process simulation software, Aspen Plus V10 (2019), was used to solve the flowsheets encountered in this work. The software was used to understand the compositions of various phases and streams involved in each stage of the IDSPP. The simulations were run on a standalone personal computer (PC) using a Microsoft Windows 10 (2020) operating system.

Data
The essential properties and data for the Barhe variant of the date palm fruits were obtained from a study conducted by Rambabu et al., [2]. The properties (Table S1, Supplementary Materials) and structure of the GA ( Figure S1, Supplementary Materials) were included in the Aspen database through the user interface [22]. The necessary user-defined functions (UDFs) pertaining to the solubility behavior of the analytes into the Sc-CO 2 /water system (for the extraction stage) were developed using Phyton (V 3.9.9) scripting and patched to Aspen software.

Method
The simulation was performed through a model-based representation of chemical, physical, biological, and technical processes and unit operations in the software. The simulation prerequisites involved a detailed understanding of the chemical and physical properties of pure components and mixtures of various streams, as well as mathematical models that allowed stoichiometry-based process calculations. The simulator described processes in flow diagrams where the unit operations are positioned and connected by streams. The solver computed the necessary mass and energy balances for each node to find a stable operating point. The thermodynamic model suitable to present the behavior of each participating phase was selected after careful analysis of a wide range of scientific literature.
For the modeling, a basis of 100 kg/h of date fruit pulp, including the moisture content, was selected. All necessary components were found and defined in the Aspen database except for the date fruit, that was specified as a non-conventional solid. The date fruit was specified using the HCOALGEN enthalpy model and the DCHARIGT density model. Electrolyte Wizard was selected for the ice formation and drying of water, and the ice formation reaction was chosen and added to the chemistry profile to be used as an input to the units involving water phase change [23]. The Particle Size Distribution (PSD) was defined to describe the appearance of the non-conventional substream. A log-normal distribution function was used to describe the PSD with a geometric mean deviation of 27.5 mm, and a mass median diameter (D 50 ) value of 27.5 mm. The IDEAL property method was selected that accommodates both Raoult's law and Henry's law and is recommended for vaccum systems where gas is assumed to behave ideally.
The FD process was modeled using a Gibbs reactor which uses Gibbs free energy minimization with phase splitting to calculate the chemical equilibrium between any number of solid components and the fluid phases. The amount of sublimated water was calculated using a UDF that was developed using the results ofİzli [14]. Accordingly, the Two Term model (as shown in Equation (1)) that best described the FD behavior was used for the freeze-dryer simulation.
where MR is the moisture ratio (M t /M 0 ), M t represents the moisture at a specific time, t (g water/g dry solid); M 0 represents the initial moisture content (g water/g dry solid); k 0 and k 1 are the drying constants with values 0.002291 min −1 and 0.4768 min −1 , respectively; a FD and b FD are FD coefficients with values 0.835 and 0.4768, respectively. These values were experimentally estimated using various thin-layer drying models for the freeze-drying of dates [14]. For SCE stage modeling, the frozen dates were demonstrated as a mixture of specific components to simulate the extraction process and evaluate the exact yield of these components from the overall fruit matrix. Among the various represented components, glucose and fructose were highlighted as the targeted entities for the extraction by Sc-CO 2 . The non-polar portion of the date fruit matrix was represented by fatty acids such as linoleic, palmitic, and oleic acids. The non-extractable material of the biomatrix was represented as cellulose due to its negligible solubility in Sc-CO 2 [12]. The phase equilibrium behavior of the glucose-fructose, Sc-CO 2 , and water mixture was thermodynamically developed. The process simulation was based on the concept of stages in the extraction. Each stage was assumed as a component separator module at fixed temperature and pressure. The Soave-Redlich-Kwong (SRK) equation of state was used to calculate the thermodynamic properties. The SRK model is suitable for high-temperature and high-pressure extractions, such as supercritical conditioning [24]. Separation unit modules were used to manipulate the separation process of CO 2 -water-glucose-fructose based on the solubility values at a given temperature and pressure. The separation was based on an empirical solubility equation (Equation (2)) proposed by Stahl and Schilz that was fed to manipulator block modules [25].
where y gf is the solubility of the glucose/fructose in the Sc-CO 2 (mol solute/mol CO 2 ); P is the total pressure of the system (bar); and A, B, C, are the empirical constants with values −4.3855, −6.2459 × 10 −3 , and 2.0307 × 10 −5 for glucose, and 4.5455, −1.2617 × 10 −3 , and 1.1561 × 10 −5 for fructose, respectively [25]. The solubilities of glucose and fructose in water were calculated according to the work of Saldaña et al., [26]. Additionally, Chrastil's equation [27] was used to model the solubility of the palmitic, oleic, and linoleic acids into Sc-CO 2 . Equation (3) presents this model, which is based on the formation of acid-complex upon association of the solute and solvent molecules.
where S is the solubility of the solute in the supercritical solvent (g/L); ρ s is the density of the pure solvent (g/L); and k, a, and b are empirical constants. Further, the solubility of water in CO 2 was obtained through regression of experimental literature data [28] and is represented by Equation (4).
where S w is the solubility of water in CO 2 (kg H 2 O/kg CO 2 ), and T is the temperature ( • C).
Modeling and simulation of the SD stage were performed through the built-in spraydryer unit block model available in Aspen Plus. The necessary parameters for the spraydryer model were specified for the drying process based on the IDEAL property method [29]. The properties for water, air, glucose, and fructose were defined with the conventional values available in the Aspen database. The properties of GA were explicitly included in the software. A rotary atomizer was chosen and modeled with a wheel diameter of 10 cm, four droplet intervals, four 30 mm blades, and a wheel speed of 18,000 rpm [19]. The characteristic droplet size distribution and nozzle exit velocity were calculated based on the geometry of the nozzle. The median droplet size (d 50 , in µm) was calculated according to the model (Equation (5)) described by Masters [30].
where n is the number of blades; m liq is the mass flow rate of the feed (kg/h); N is the rotational speed (rpm); D R is the diameter of the atomizer (m); and h is the height of the vanes (m). The constant K and exponents p, q, r, s were determined from the correlations for the vane liquid loading [30]. Subsequently, the mass flow ratio/vane (M s ) to the height of the vanes was determined as 3587.1 kg/h m. The validity limits were expressed in terms of the speed of the circumference of the disc (U) using Equation (6).
Thus, the speed of the circumference of the disc was evaluated as 94.25 m/s. As the values of the M s and U values were in the same category, the design was considerably adequate according to the validity table [30]. The values of the constants and exponents were obtained from the same category. Based on these calculations and validations, the median droplet size (d 50 ) was estimated as 50.3 µm. Additionally, the critical moisture content for the drying process was evaluated as 2.67 kg water/kg dry solid (Equation (1)). The equilibrium moisture content was taken as 0.03 kg water/kg dry solid. Table 1 summarizes the values for various parameters used for the simulation of each stage. The sensitivity analysis tool of the software was used to study the effect of different critical parameters on the final extracted amount of glucose and fructose.

Results and Discussions
Steady-state modeling and simulation of the IDSPP were conducted successfully using Aspen-Plus Software. Each sub-process was modeled and simulated separately. In addition, sensitivity analysis was performed for each block to identify the optimum operating conditions of the critical parameters for the respective stage.

Freeze-Drying
In the FD method, the biological qualities of the components, flavors, and colors of the food are usually preserved well [31]. For the analysis of IDSPP, vacuum-based FD was considered for the effective removal of moisture at a minimal drying time. The FD model accounted for the phase changes and water vapor diffusion within the date pulp. Modeling for the non-condensable species was accounted through the species transport model [32]. Technically, the feed sample was categorized into four layers, namely frozen front, vapor-ice interface, dry zone, and gas phase. A source term was used to adjust the distribution of gas-phase species due to water vapor generation [33]. Figure 2 depicts the process flow diagram of the FD process. Based on the simulation results, it was observed that the moisture content of the date pulp was reduced from 19.62% to 0.26%, calculated on a dry weight basis. The mass density of the freeze-dried dates was estimated as 1,639.05 kg/m 3 . This relatively low-density value of the freeze-dried date pulp signified higher porosity and thus a potential for improved cell rupture in the SCE stage [32]. Further, sensitivity analysis was performed by varying the temperature and pressure for the FD to examine the effect on the moisture content and duty consumption. The results for the optimization studies of FD are displayed in Figure S2

Supercritical Extraction
An Aspen flowsheet for the SCE process is illustrated in Figure 3. For the proper design of the SCE stage, it was essential to employ appropriate mathematical representations. The SCE stage involved extraction from various solid natural matrices that were thermodynamically represented through solubility control, external mass transfer control, and/or internal mass transfer conventional models. Notably, three manipulators were used in the modeling of the SCE stage. The first one, "EXTRCALC" was used to estimate the quantity of the analytes to be extracted by the Sc-CO 2 /water solvent system based on the UDF patched to the software. The UDF contained the necessary data and equations for the solubility behavior of glucose, fructose, and other essential components of date fruit into the solvent mixture. The second one, "CO2OUTCA" calculated the amount of CO 2 lost due to the opening of the extraction cell to fix the overall mass balance of the system. The third manipulator, "SEPARATO" was used to estimate the extracted mass lost with the exhaust stream of CO 2 gas.

Operating Pressure
The extraction pressure of the SCE process is an important parameter that impacts the overall extraction efficiency [34]. Figure 4 shows the pressure effects on the extraction of reducing sugars from the freeze-dried date pulp. Simulation studies revealed that the extraction pressure had varying effects on the recovery of glucose and fructose sugars. For fructose sugar, the recovery efficiency decreased with incremental values in the operating pressure (Figure 4a). On the contrary, the extraction yield increased for glucose sugar with an increase in the operating pressure (Figure 4b). This opposing behavior can be ascribed to the different signs of the empirical coefficients in Equation (2) for glucose and fructose solubility in Sc-CO 2 . Consequently, the trade-off optimization tool of Aspen software was utilized to estimate the optimum pressure required for the maximum extraction of the total reducing sugars. Results showed a maximum product flow (glucose + fructose) of 73.70 kg/h was achieved at an operating pressure of 308 bar.

Extraction Temperature
Similar to pressure, the operating temperature of the SCE process is an imperative parameter that influences the efficiency of the separation process. Simulation studies showed that an increase in the extraction temperature had no significant effects on the recovery of reducing sugars from the fruit pulp. This insignificant temperature effect was anticipated as the solubility equation (Equation (2)) for the monosaccharide sugar molecules into Sc-CO 2 did not involve any temperature effects. However, the extraction of non-polar analytes such as linoleic acid, palmitic acid, and oleic acid by the solvent showed a significant increase with temperature for T > 70 • C ( Figure S3, Supplementary Materials). The increase in the temperature enhanced the kinetic energy of the Sc-CO 2 fluid, which in turn improved the solvent diffusion and bioactive compound recovery [35]. However, it is worthy of specifying that higher extraction temperatures would cause the thermal degradation of the recovered sugar molecules, reducing the overall yield [36]. Based on the simulation results, an extraction temperature of 65 • C was recommended for the SCE stage.

Solvent Flow Rate and Composition
Results on the sensitivity analysis for the solvent:co-solvent ratio showed that a cosolvent concentration of 6.77 wt.% produced the maximum recovery of reducing sugars from the date pulp biomatrix. The simulation run resulted in an overall extraction yield of 99.1% for the reducing sugars with flow rates of 34.14 kg/h and 39.13 kg/h for glucose and fructose, respectively. Moreover, the sensitivity analysis showed that the total extracted quantity of glucose and fructose increased with the flow rate of Sc-CO 2, as shown in Figure 5. Almost all of the reducing sugars were extracted into the solvent system for an Sc-CO 2 flow rate of 31,000 kg/h. Further increases in the Sc-CO 2 flow rate had no significant effects on the sugar extraction. This quantity of Sc-CO 2 defines the solvent limit required for the maximum removal of targeted sugar components from the date palm fruit. Based on the optimization studies for the solvent:co-solvent ratio and the Sc-CO 2 flow rate, the optimum flow rate of water was estimated as 2100 kg/h for the maximum recovery of glucose and fructose components. Figure 6 shows the effect of the water flow rate on the extraction of reducing sugars from the fruit biomatrix. Similar to Sc-CO 2 , the extraction yield increased with incremental levels of water flow rate. The improvement in the extraction of monosaccharides was due to the enhanced permeability and polarity of water molecules at the extraction conditions.
In general, the mass transfer mechanism for the recovery of plant bioactive compounds using supercritical technology involves four major steps [37]. Firstly, the analytes are desorbed from the biomatrix into the solvent film. Second, the analytes diffuse into the solvent, while in the third, the analytes dissolve into the solvent. The final step involves the bulk mass transfer of the solvent from the biomatrix to the external. Based on this mechanism, the increase in the extract yield with the flow rate of the solvent system signified that the external film transfer diffusion (second step) quantified by a diffusion coefficient, and the dissolution potential of sugars into the Sc-CO 2 /water system (third step) quantified by a partitioning coefficient, played the limiting roles for the recovery of the targeted analytes [38].

Spray-Drying
In general, beneficial features such as high drying rate, shorter drying time, and broad working temperature make SD an ideal choice for drying of food products to powder form [32]. The complicated three-dimensional nature of spray-dryers often makes it difficult to describe their exact behavior using empirical models. In this work, Aspen's builtin simulator was used to characterize the SD process of the SCE extract (glucose and fructose solution) mixed with GA as the carrying agent. The model employed an Eulerian-Lagrangian approach to compute the particle motion paths, and heat and mass transfer phenomena occurring between the drying medium and the feed particles [32]. The choice of GA was to overcome the stickiness effect and increase the glass transition temperature of the mixture. Based on the previously reported experimental studies of SD for date fruit powder, the mass ratio of GA to SCE extract was fixed as 0.4 [18]. Figure 7 represents the process flow diagram for the SD stage of the IDSPP technology. Sensitivity analysis of the SD stage involved the study of the effect of different parameters on the final moisture content of the sugar powder. Results showed that the final moisture content of the powdered product decreased with an increase in the mass flow rate of the drying air ( Figure S4, Supplementary Materials). Complete removal of excess moisture was observed at an airflow rate of 8000 kg/h resulting in powdered sugar with the equilibrium moisture content. High flow rates of the drying air reduced the droplet size due to their improved atomizing effect. This resulted in a finer particle size of the product with a larger surface area that facilitated the rapid and higher drying rate of the product with complete removal of moisture [30]. It is worth specifying that the blank simulation without GA addition resulted in low drying efficiency. This can be attributed to the possible stickiness of the feed extract, as reported in several other related experimental works [39]. Figure 8a shows the effect of the inlet air temperature on the moisture removal from the feed solution for the powdered product formation. It was observed that an increase in the air temperature resulted in improved moisture elimination. This was clearly due to the improved heat and mass transfer processes between the air and feed particles owing to the greater driving force for moisture evaporation [40]. Moreover, higher operating temperatures would eliminate insufficient drying of the particles, thereby reducing the probability of semi-wet/sticky solids formation that would typically stick to the wall and thus result in higher yield [39]. Simulation results showed the attainment of equilibrium moisture content at a minimum air inlet temperature of 150 • C for the formation of dry powdered products of date sugar. Further, the results for the dependency of product moisture content on the feed flow rate are shown in Figure 8b. It was evident that lower feed flow rates resulted in a relatively drier product than higher flow rates. In the SD stage, a higher feed rate may result in slower heat and mass transfer processes. An increase in feed flow rate reduces the contact time between the drying medium and the atomized particles leading to low heat transfer rates resulting in low moisture removal. In general, feed flow rate not only affects the moisture content and yield of the product but also significantly influences the product solubility and water activity [41]. Based on the simulation results, a feed flow rate of 435 kg/h was recommended for the SD stage of the IDSPP.

Product Solubility
Product solubility is a critical benchmark for food products, especially food-grade sugars, to be used as an additive in food and beverages. In general, several parameters of the product, such as composition, water activity, morphology, equilibrium moisture content, etc., affect the product solubility [42]. The solubility of the product (glucose and fructose mixture) in an aqueous medium was simulated over a temperature range from 0 to 100 • C and at a pressure of 1 bar based on the UNIQUAC (universal quasichemical) thermodynamic model. The solid-liquid phase equilibrium of the powdered date sugar in pure water was simulated using activity coefficients [43]. The solubility of the sugar product at a temperature of 25 • C and 1 bar was estimated as 0.89 g/g water. Comparing this value with the solubility of the table sugar (2.1 g/g water), it was observed that sucrose (in table sugar) had a higher solubility than the extracted date sugar. However, it is expected that experimental validation would result in higher solubility of the date sugar powder due to the good emulsification properties and high solubility of GA in aqueous solutions [22]. Figure 9 depicts the PSD for the final powder product obtained for the optimal operating conditions of the SD step. The mean particle size and density were estimated as 504 µm and 0.0023 kg/m 3 , respectively, demonstrating the presence of significant levels of bioactive compounds in the product [42]. In addition, powders obtained using lower drying temperatures tend to possess a high percentage of smaller particles in the micro-range.  Experimental validation of IDSPP technology, with extensive analysis of extracts (including other vital analytes such as fibers, minerals, amino acids, etc.), determination of the extraction equilibrium and kinetics [44][45][46][47], and establishment of the necessary economics for production along with regulatory waste management approaches [48][49][50] would potentially realize effective production of date sugar that could ideally replace refined sugar.

Conclusions
In summary, this work conceptualized a novel IDSPP technology comprising of three main stages, namely freeze-drying, supercritical extraction, and spray-drying to process 100 kg/h of date fruit pulp to produce a soluble date sugar product. Stage-wise modeling and simulation studies were performed using Aspen Plus software with the sensitivity analysis to identify the optimum values of critical parameters for each stage. Freeze-drying studies revealed a minimum duty requirement of 13.23 kW at a temperature of −42 • C and pressure of 0.0001 bar to reduce the moisture content of the fruit pulp from 19.62% to 0.26%. The Sc-CO 2 assisted extraction resulted in a total sugar recovery of 99.1% with flow rates of 34.14 kg/h and 39.13 kg/h for glucose and fructose sugars, respectively. Optimization studies showed that a pressure of 308 bar, temperature of 65 • C, Sc-CO 2 flow rate of 31,000 kg/h, and Sc-CO 2 : water ratio of 0.07 resulted in the maximum extraction efficiency. Further, spray-drying simulation of the date extract into a powdered product using Gum Arabic as the carrier agent showed a mean particle size of 50.4 µm for the date sugar powder with a solubility of 0.89 g/g water (at 25 and 1 bar) and density of 0.0023 kg/m 3 . Sensitivity analysis revealed that the extraction pressure and solvent flow rates were the significant parameters for the extraction step while operating temperature had no salient effects on the recovery of the sugar. For the spray-drying process, drying air flow rate and its inlet temperature had positive effects, while feed flow rate negatively impacted moisture removal. Results showed that a mass flow rate of 8000 kg/h and an inlet temperature of 150 • C for the drying air was optimal to dry a feed stream mixture of date extract and Gum Arabic flowing at 435 kg/h containing 40 wt.% of Gum Arabic. Thus, this work confirmed the technical feasibility of IDSPP for the production of soluble date sugar powder and provided the necessary data with optimum operating conditions required for the process. Experimental validations and detailed characterizations of the product are further required to realize the outcomes of this study.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/pr10020257/s1, Figure S1: Chemical structure of Gum Arabic, Figure S2: Sensitivity analysis for FD: Temperature and pressure effects on the final moisture content and the net heat duty, Figure S3: Effect of temperature on the recovery of non-polar analytes of date fruit pulp, Figure S4: Effect of drying air flow rate on the moisture content of the date sugar powder product, Table S1: Properties of Gum Arabic used for the ASPEN simulation.