Design and Pinch Analysis of a GFT Process for Production of Biojet Fuel from Biomass and Plastics

: Environmental problems are frequently related to energy use, estimated to grow at 1.6% per year until 2035. The transport sector accounts for 30% of energy demand and aviation is growing around 2.6% per year. Thus, low-emissions policies promote the use of sustainable aviation fuels. This work simulates a gasiﬁcation and Fischer-Tropsch process to obtain biojet fuel from biomass and plastic waste. Syngas obtained through cogasiﬁcation is puriﬁed by amine scrubbing and subjected to a Fischer-Tropsch process to produce hydrocarbons, which are upgraded for optimal fuel properties. Pinch analysis is applied to minimize energy usage, while Rankine cycles and a cooling tower are designed to cover the demand of electricity and cooling water. Results show that mass yields of the process towards biofuels are 13.06%, with an output of 1697.45 kg/h of biojet fuel. Density, kinematic viscosity, pour and ﬂammability points and the lower caloriﬁc value of the biojet fuel comply with the ASTM D7566 standard. Pinch analysis allows to reduce 41.58% and 100% of cooling and heating demands, respectively, using biomass as renewable energy for heating. Moreover, steam generation covers 38.73% of the required electricity. The produced biojet fuel emits 20.14 g CO2eq /MJ and has a minimum selling price of 1.37 EUR/L.


Introduction
Energy is the engine of society and one of the fundamental pillars of human progress. Global energy usage has increased significantly thanks to the large-scale incorporation of fossil fuels [1], which has led to a significant increase in living standards. Today, fossil fuels are still the main source of energy used worldwide, accounting for approximately 90% of energy demand [2]. However, the use of these raw materials as an energy source is the cause of the main environmental problems affecting the planet (climate change) and, due to their geographic location, has caused conflicts based on energy dependence and geopolitical strategies. According to the United Nations, at least 40% of violent conflicts in the last 60 years are related to the exploitation of strategic resources [3].
Renewable energies represent an alternative to conventional energy sources and biomass stands out among them, due to its ubiquity, availability, and low polluting nature [4]. Due to the increase in the price of fossil fuels and materials (in the United States the price of petrol between 2000 and 2021 has risen from 1.52 to 2.26 USD/gal [5]), in recent years, biorefineries present an interesting alternative, being able to efficiently convert renewable raw materials into a wide range of marketable products of biological origin.
Besides, nowadays plastic waste is another major environmental problem [6,7], as it ends up as waste that accumulates in landfills or is dumped into the environment. Plastics

Location and Feedstock
In relation to sustainable aviation fuel, Europe is the continent that is making the most significant commitment [34], given its CO 2 emission reduction targets and its dependence on fossil fuels from other countries. While the geographical location of the plant has little influence on its design, it can play an important role in the economic analysis. The biojet fuel plant is located in Spain, because it is one of the main agricultural producers [35] and, therefore, a promising biojet fuel producer from biomass waste blended with plastics.
Residual biomass: In Spain, straw is the most abundant residual lignocellulosic biomass of agricultural origin [36]. Specifically, in the province of Huesca, wheat and barley straw production amounts to 178,002 and 579,448 tonnes per year (t/year), respectively. Although the total utilisation of wheat and barley residues available in the province is optimistic, other types of residual lignocellulosic biomass are generated and must be considered. Therefore, the mass flow of this raw material is considered achievable.
Plastic waste: Spain produces 1.5 million tonnes of plastic waste per year, of which approximately 20% is exported to other countries [9]. In order to estimate the quantity and composition of the plastic used as raw material, the following assumptions are made:

•
The amount of waste generated per inhabitant is homogeneous throughout Spain; • The average composition of plastic waste is based on the type of polymers demanded in Europe [8]; • A total of 90% of the polyethylene (PE) and polypropylene (PP) that would be exported from the provinces of Huesca, Zaragoza, Navarra, and La Rioja are fed into the chemical process.
For a plant operating factor of 8000 h per year, the raw material quantity and composition by mass (wt%) are shown in Table 1:

Process Simulation and Thermodynamic Models
Aspen Plus and Aspen HYSYS v10 have been used to simulate the process. The amount of raw material and the ratio between lignocellulosic biomass and plastic is based on Table 1. The biomass components (cellulose, hemicellulose, and lignin [37][38][39][40]) are entered as user-defined compounds using their molecular structure [4] for the software to estimate their physical properties. The software already contains polyethylene and polypropylene in its "APV100 POLYMER" database, which is used to simulate the plastic properties.
Furthermore, Aspen Plus has been used to model the physical properties of biomass and plastic applying NRTL [41] as thermodynamic model. Mass balance calculation for the distribution of pyrolysis products has been carried out by using Microsoft Excel. The gasification reactor and downstream processes have been simulated via Aspen HYSYS: using NRTL for gasification (unit 1, see Figure 1), FT process and upgrading (units 3 and 4, see Figure 2) [42]; Acid Gas-Chemical Solvents for amine scrubbing (unit 2, see Figure 1 and unit 5, see Figure 3) [43], and SRK for product separation (unit 6, see Figure 3) [44].

Heat Integration and Utilities
Energy integration between process streams has been studied using pinch analysis, performed via Aspen Energy Analyzer v10 software. Then, Rankine cycles (unit 7) were designed using Peng-Robinson [45], in Aspen HYSYS to cover part of the electrical demand, using the surplus steam generated.
Finally, EBSILON Professional v14.2 software was used to simulate a cooling tower. This software allows us to calculate the mass balance and the electrical power required by the equipment.

Economic Analysis
Total investment was estimated by the percentage method [46], based on equipment and raw materials costs. Minimum selling price (MSP) of biojet fuel and coproducts is assessed considering an internal rate of return (IRR) of 9% [47].

Process Flow Diagram (PFD)
The chemical process starts in unit 1, shown in Figure 1 (refer to Table S9 of Supplementary Information to see the complete list of streams and equipment in PFD), where biomass (S101) and plastic (S102) are introduced into a fluidised bed boiler to be gasified in an oxygen-poor atmosphere (S109) to produce syngas. Steam (S105) is introduced to adjust the amount of hydrogen (H 2 ) generated. On one hand, the syngas leaves the boiler (S112) and is transported to unit 2. On the other hand, ashes and part of the bed-forming solids are removed from the boiler (S110).
Syngas enters unit 2, shown in Figure 1, where it must be purified to reduce its concentration of carbon dioxide (CO 2 ), as it poisons the Fischer-Tropsch reactor catalyst. Water (H 2 O) and other condensable compounds (S203) are first separated from the syngas in a flash chamber. Recovered water is taken to a treatment plant so that it can be reused in the process. Gases separated in this flash chamber (S202) are taken to an absorption column, where an aqueous solution of monoethanolamine (MEA) is used to absorb the CO 2 , leaving at the bottom as rich amine (S208).
Rich amine is taken to a flash chamber (V-203) and a rectification tower (T-202), where the CO 2 is separated and emitted to the atmosphere (S218).
The regenerated amine, known as lean amine, leaves the distillation tower and is taken to a mixer, where its composition is adjusted by the introduction of water and fresh amine, conditioning it back to the absorption column (S206).
The purified gas stream leaves the absorption column at the top (S207) and is fed to a pressure swing adsorption (PSA) bed, where the ratio of H 2 to carbon monoxide (CO) is adjusted. The syngas (S226) is fed to unit 3, and the separated hydrogen (S225) to unit 4.
In unit 3, shown in Figure 2, syngas is brought to Fischer-Tropsch reactor operating conditions. In the reactor, syngas is transformed into water and straight-chain hydrocarbons (S305), which are transported to unit 4.
In unit 4, products obtained from unit 3 are taken to a 3-phase separator, where the gaseous products (S403) are sent to unit 5. The water generated is taken to the water treatment plant (S405), and the hydrocarbons (S404) to the upgrading reactor, to react with the hydrogen separated in the PSA of unit 2 (S225).
Final products are obtained in units 5 and 6. In unit 5 (see Figure 3), fuel gas from the three-phase separator of unit 4 (S403), with a high carbon dioxide content, is purified by procedures similar to those of unit 2. The purified gases (S503) are the final product.
Rich amine (S504) must be regenerated and, for this purpose, it enters a rectification column where CO 2 is separated and emitted to the atmosphere (S507). The lean amine leaves the equipment and is taken to a mixer where water and fresh amine are added in order to replace the losses produced in the unit and adjust its composition, before being recirculated to the absorption column.
Unit 6 consists of a set of columns that allow the separation of the different products obtained. These products are fuel gas (S610), petrol (S611), biojet fuel (S619), diesel (S627), and fuel oil (S629).

Gasification (Unit 1)
Feedstock is fed into the gasification reactor, where the processes of pyrolysis, oxidation, and reforming take place. Pyrolysis is endothermic and starts at around 300 • C for biomass and 450 • C for plastic [22,48]: Cellulose + Hemicellulose + Lignin → Gas + Tar + Char Plastic → Gas + Tar In the case of biomass, the gas stream is mainly composed by carbon monoxide, carbon dioxide, methane (CH 4 ), ethylene (C 2 H 4 ), and hydrogen. Char fraction is composed of carbon (C), ash, extractables and other solids. The tar is made up of hydrocarbons with an average composition of C 10 H 10 and oxygenated compounds with an average composition of C 6 H 10 O 3 [48].
In the case of plastics, the gas stream is considered to be composed of hydrogen, ethylene, and 1-butene (C 4 H 8 ). The produced tar is of two types: primary tar, consisting of aromatic hydrocarbons equivalent to cyclohexane (C 6 H 12 ), and secondary tar, consisting of olefins of average composition C 10 H 20 [48].
Pyrolysis mechanism is very fast and usually considered to occur instantaneously. Thus, a mass distribution of compounds formed after pyrolysis at 600 • C is calculated using the following equation [48]: where a 0 , a 1 , a 2 , and a 3 are parameters from Table 2. Temperature (T) is expressed in K for biomass and in • C for plastic. Average specific heat (c p ) and specific enthalpy (h) of the species involved in the pyrolysis have been obtained from Aspen Plus (see Table 3). The mass of reactants is the mass of feedstock minus extractables and ash of S101 (see mass balances in Supplementary Information Section A): The heat required for pyrolysis reaction is calculated using the equation: where ∆T is the temperature difference between the pyrolysis (600 • C) and the reference temperature (25 • C). The heat required for the reaction is 106 MW, which is the minimum heat to be released during gasification to make the reactor autothermal. After pyrolysis, the gasification takes place, in which the compounds shown in Table 2 react with each other and the gasifying agent to form syngas. In this study, a sufficient amount of air is introduced so the exothermic mechanism of gasification compensates the endothermic mechanisms of pyrolysis and reforming [49,50].
The reactions that take place in gasification and their kinetics are as follows (r is the reaction rate (kmol/(m 3 ·s)); the concentrations of referenced compounds are in kmol/m 3 ; T (K); E A is the activation energy (kJ/kmol) and p O2 , p H2O is the partial pressure of O 2 and H 2 O (Pa)):

•
Oxidation of gaseous compounds [51]: • Tar oxidation [48]: • Char oxidation [48]: Direct r = 1.51 × 10 6 ·e − 108.67 Reverse • Reforming [48]: Reverse Once unit 4 (see Figure 2) has been simulated, the amount of steam to be introduced into the reactor to generate the hydrogen required is known to be 2600 kg/h. On the other hand, 38,200 kg/h of oxygen must be fed to make the reactor autothermal, producing 3.6 MW of heat that must be removed to maintain isothermal conditions.

Amine Scrubbing (Units 2 and 5)
There are two streams with a high concentration of carbon dioxide: the syngas from gasification and the Fischer-Tropsch produced gas. The syngas obtained after gasification (S112, see Figure 1), with 70.4 wt% CO2 that can poison the cobalt (Co) catalyst used in the Fischer-Tropsch (F-T) reactor [25]. On the other hand, the fuel gas obtained from F-T (S403, see Figure 3), with 62.4 wt% CO2 reduces the calorific value of the stream.
MEA is selected as it requires a flow rate several orders of magnitude lower than other amines [26] for keeping CO 2 around 1 n/n% (mole percentage) in the clean gas stream. Before the recycle is introduced, the composition of the lean amine entering the absorption tower (S206/S502, see Increase of pressure improves the absorption and favours the formation of syngas. Thus, taking into account the pressure drop between R-101 (see Figure 1) and T-201, an operating pressure of 65 bar is selected for the absorber. In the case of T-501, the selected pressure is 3 bar, which is the maximum discharge pressure of K-501 compressor for which the gas temperature is less than 100 • C, as this would cause the water in the amine solution to boil.
Regarding the operating temperature, absorption improves at low temperatures, so both absorption towers operate at a temperature of 50 • C. Absorbers (T-201/T-501, see Figures 1 and 3) are designed with 15 and 10 stages, respectively, to obtain less than 1 vol% CO2 in the clean gas.
In unit 2, a flash chamber (V-203, see Figure 1) can be used for CO 2 separation, because there is a large pressure difference between the absorption tower and the regeneration tower (operating at a pressure of 2 bar). This separates 22.5 wt% of the CO 2 (S210).
Regeneration towers (T-202/T-502, see Figures 1 and 3) operate at atmospheric pressure in the condensers and at 2 bar in the reboilers. These towers are designed for high CO 2 and amine recovery while avoiding high reflux and vaporisation ratios in both condenser and reboiler.
Condensers of all rectification towers in the plant operate at full reflux and temperatures higher than 50 • C, so water from cooling tower can be used for refrigeration. Table 4 lists the main operating parameters of the regeneration towers: In the amine syngas scrubbing unit (unit 2, see Figure 1), 98 n/n% of the carbon dioxide is recovered in the distillate, while 100 mol% of amine and 75 mol% of water are recovered in the residue. Sufficient water is generated in the process to replace these losses, so reducing the operating costs in the column is preferable to having a higher water recovery. The separated carbon dioxide (S214) is combined with the carbon dioxide obtained in V-203 (S210) and emitted to the atmosphere (S218). In the case of the amine fuel gas scrubbing unit (unit 5, see Figure 3), the carbon dioxide released is emitted directly.
After the regeneration towers, fresh water (S222/S513) and amine (S221/S514) are added in order to adjust the lean amine flow rate and composition. Then, lean amine (S206/S502) is conditioned and sent back to the absorption tower (T-201/T-501).
Following the absorption tower in unit 2, the clean syngas (S207) is sent to a PSA (V-202A/B, see Figure 1) which is used, on the one hand, to separate part of the hydrogen that is fed to the upgrading reactor (S225) and, on the other hand, to adjust the ratio of carbon monoxide to hydrogen in the syngas (S226) that is sent to the Fischer-Tropsch reactor.

Fischer-Tropsch Process (Unit 3)
The mass distribution of products obtained from this process (defined by Equation (19)) is calculated by the term "chain growth probability" (α) [28], which is a function of the type of catalyst used and the reaction conditions.
where n is the carbon number in the paraffin formed [52]. The chain growth probability must be iterated in order to maximise the production of biojet fuel, as it determines whether the hydrocarbons formed are light or heavy. Once 0.88 has been obtained as the optimal α value, the results obtained are summarised in Figure 4, where the "Products" series has been obtained by applying Equation (19). The width of rectangles in "Fractions" series resembles the carbon number of hydrocarbons considered to be contained in each fraction (C 8 -C 16 for biojet fuel) and their height is the sum of the mass distribution of the hydrocarbons contained (44.76 wt% of biojet fuel). Reactions that take place for the production of straight-chain paraffins are: Therefore, combining the probabilities of obtaining each hydrocarbon with the stoichiometry of each reaction, it is possible to calculate the optimal H 2 /CO mole ratio, which has been determined at 2.08.
The operating temperature of the reactor is calculated taking into account the partial pressure of each species in the reactor, as shown in Equation (20): where n H 2 n CO opt is the optimal H 2 /CO mole ratio in the syngas. Therefore, the optimum operating temperature is 211 • C. Mass composition of products from the Fischer Tropsch reactor (R-301, see Figure 2) is given in Table 5: Table 5. Composition of S305 (see Figure 2). Finally, using the enthalpy of the inlet and outlet streams and eliminating the sensible heat term in Equation (2), the heat given off in the reactor is estimated as 51 MW, which must be removed to maintain isothermal conditions.

Upgrading (Unit 4)
Linear hydrocarbons are taken to the upgrading reactor, where hydrocracking and isomerization take place, increasing the yield of light products and improving their properties for their commercialisation [53].
Inside the reactor, operating at 250 • C and 30 bar, the following reactions (with the correspondent mass yield η) occur [29]: The amount of remaining reactants and formed products are estimated from mass yields and stoichiometric ratio of the different reactions. It is assumed that all isoparaffins formed have a single branch, with the methyl group on the second carbon of the chain [30][31][32][33]54,55].
Finally, using Equation (2) for the upgrading reactor (R-401, see Figure 2), the heat that must be removed to maintain isothermal conditions is estimated as 196 kW.

Product Separation (Unit 6)
The objective of the unit is to maximise the biojet fuel fraction, which comprises hydrocarbons from C 8 to C 16 . However, n-octane (n-C 8 ) and isooctane (i-C 8 ) improve the properties of petrol, which contains hydrocarbons from C 5 to C 11 . Therefore, in the first column (T-601, see Figure 3) the aim is to recover n-C 7 , which forms part of the petrol, and n-C 9 , which forms part of the biojet fuel (both are set to 99.5 n/n%). Under these conditions, the column has a certain flexibility of operation in which 90.7 n/n% of the C 8 is obtained in the distillate.
Light fractions are led to a flash chamber (V-601) to separate petrol from fuel gas, while heavier hydrocarbons are taken to two other rectification towers (T-602 and T603). T-602 is used to separate the biojet fuel from the heavier compounds; so, the recovery of i-C 16 in the distillate and n-C 17 in the residue is close to 100 n/n%. Finally, T-603 is used to separate diesel (C 15 -C 20/25 ) from fuel oil (C 20+ ), setting the recovery of n-C 19 (head) and n-C 21 (bottom) at 90 n/n%. Table 6 lists the main parameters of the rectification columns: Table 6. Main parameters of the rectification columns (T-601/T-602/T-603, see Figure 3). The mass flow between units is shown in Figure 5 and Table 7 (see Figures 1-3   (1) Fuel gas (S503 + S602), petrol (S611), diesel (S627) and fuel oil (S629), see SI E-F for detailed information.

T-601 T-602 T-603
As shown, 103,994 kg/h of biomass and plastic are used as raw material, which is converted into 1207 kg/h of fuel gas, 1590 kg/h of petrol, 5676 kg/h of biojet fuel, 1697 kg/h of diesel, and 3416 kg/h of fuel oil, which means a yield of 13.06% towards products. This yield could be slightly improved by adjusting the supply of oxygen as a gasifying agent, since CO 2 production has been increased to make the reactor autothermal.
All mass balance data sheets of the chemical process can be found in Supplementary Information Sections A-F.

Heat Integration
Aspen Energy Analyzer is used to optimize energy demand. Stream properties were imported from Aspen HYSYS. The heat load released or given off by reactors, which require cooling or heating, is calculated for a difference of 1 • C from the operating temperature (see Supplementary Information Table S1). Prior to the pinch analysis, the heating and cooling demands are, respectively, 55,098 and 296,573 kW.
For the simulation of heat exchanger networks, Aspen Energy Analyzer allows the input of additional parameters such as viscosity (µ), thermal conductivity (k), and density (ρ), which have been obtained directly from Aspen HYSYS. These parameters and the flow velocity through the exchanger (v) are used to calculate the individual heat transfer coefficient of the fluid. This velocity is set to 15 m/s for gases and 3 m/s for liquids or phase changes (see Supplementary Information Table S2).
Hot streams, whose temperature need to decrease and require cooling, form the Hot Composite Curve (HCC). Cold streams, whose temperature need to increase and require heating, form the Cold Composite Curve (CCC). Plotting both curves in a T-H diagram (see Figure 6), there is a maximum approach, represented in green, which is given by the minimum admissible temperature difference (∆T min ) in the heat exchangers [56]. This approach marks the optimal heating and cooling demands, from a thermodynamic point of view [57]. The design of the heat exchanger network and, therefore, the potential energy integration, implies the need to fix this minimum temperature difference (∆T min ), so the heat exchanged is maximised without the size of equipment being excessive [58]. Using the Aspen Energy Analyzer software, the relationship between fixed and operating costs related to heat exchangers and utilities can be estimated, as shown in Figure 7:  Figure 6 shows that the whole CCC lies below the HCC, which means that there exists at least one heat exchanger network that fully satisfies the heating demand of the system. The excess heat from the process could be used for steam generation, which can be sold, so total costs shown in Figure 7 are negative, which indicate that the design of a heat exchanger network will have a positive impact on the economics of the project, which is an improvement to the GFT process. A ∆Tmin of 6 • C is selected as economically optimal (see Figure 7) [57].
Once the ∆Tmin is selected, the Grand Composite Curve (GCC) is drawn in Figure 8. The GCC allows to know pinch temperature (847 • C), and thermodynamically optimal demands (0 MW for heating and 170 MW for cooling): In addition, the GCC is used to select the utilities for the heat exchanger network, as shown in Table 8: (1) Given the high temperatures, it is possible to have steam generation. Steam is superheated 10 • C to avoid condensation during transport due to heat losses, so this heat is assumed to be dissipated.
Thus, knowing the demands, flow rates and properties of all streams and temperature of utilities, the heat exchange network is designed using the Aspen Energy Analyzer software (see Supplementary Information Figure S10).
The steam demanded by the cold streams is generated from the hot streams, so that the final heating and cooling demands of the chemical process are set at 173 MW for cooling (2.12% higher than the optimal) and 0 kW for heating; thus, the network manages to reduce the initial demands by 41.58 and 100%, respectively. By oversizing the necessary cooling water flow by 10%, the required mass flow rate of this utility is 6,830,113 kg/h.

Thermodynamic Cycles (Unit 7)
The excess steam generation, obtained from the heat exchanger network, is used to partially cover the electricity demand of the plant, calculated as 26,146 kW (see Supplementary Information Section L and unit 8 electricity demands in Section 3.4.2. Aspen HYSYS is used to design the different Rankine cycles for each steam pressure. All cycles start in a heat exchanger where steam is generated (E-950, E-G HPS , E-G MPS and E-949, see Figure 9). This steam is taken, if required, to the heat exchanger network (E-939, E-952) and pumped (P-701, P-703) to the operating pressure of steam generators.  Table 9 shows heat used as an input for evaporators: Table 9. Inputs of the evaporators.

Cycle Heat Exchanged (kW) Evaporators
Very Likewise, the inlet and outlet temperatures for the steam at the different pressures are used as initial data: 259-270 • C for HHPS, 249-260 • C for HPS, 174-185 • C for MPS and 124-135 • C for LPS. For each steam type, the pressure at which evaporation occurs (at 260, 250, 175, and 125 • C respectively), is set at the steam outlet after the pressure drop in the evaporator.
In the high and very high-pressure cycles, flow dividers are used as part of the steam is sent for heating. A cooler is used to simulate 10 • C of temperature loss in transport. Moreover, exchangers that demand the steam are placed downstream, where the heat exchanged is introduced as an input, so that Aspen HYSYS can calculate the ratio of steam to be turbined and sent as utility.
The pressure drop in the turbine is limited so steam fraction at the turbine outlet is higher than 90%, or the outlet pressure higher than 0.4 bar.
Condensers (E-701, E-702, E-703, E-704, see Figure 9) use cooling water (CW) from the tower, so the inlet and outlet temperatures are assumed to be 28 and 42 • C, respectively. The data shown in Tables 10 and 11 are obtained after the simulation of the cooling tower (unit 8, see Figure 10):   Cooling water flow required by the cycle, oversized by 10%, reaches 4,448,986 kg/h, which is added to the 6,830,113 kg/h demanded by the chemical process, so that the total flow of water to be cooled by the cooling tower is 11,279,099 kg/h. After simulating the cycles in Aspen HYSYS, the net power required is 10,431 kW, which means that the steam generated does not cover all the electrical demand of the chemical process equipment, so 16,698 kW must be purchased.

Cooling Tower (Unit 8)
This unit (see Figure 10) is simulated using Ebsilon Professional and a mechanical forced draft tower with countercurrent flow [59].
The hot water enters the cooling tower (S801, see Figure 10), is atomised, and then flows down. Air (S802), at ambient conditions, is introduced at the bottom and flows upwards, gets in contact with hot water and leaves the tower as saturated air (S803), cooling the liquid water, which falls into a pool at the bottom. The water losses increase the concentration of salts in the cycle over time, so part of the cooled water (S804) must be purged. The rest is recirculated (S805) and mixed with fresh water (S806) to make up for the losses, restarting the cycle.
Atmospheric pressure and a temperature of 25 • C for all raw materials are used as input. In addition, the average relative humidity of the air is assumed to be 75% [60]. Results of the cooling tower simulation are shown in Table 12: The flow of fresh water taken from the medium represents approximately 3% of the refrigerated water flow. The mass balance data sheets of the thermodynamic cycles and the cooling tower can be found in SI G-H.

CO 2 Equivalent Emissions and Biojet Fuel Properties
On the one hand, emissions from biomass combustion are considered biogenic and not accounted in the biojet fuel life cycle. On the other hand, the gasification reactor is fed with 9 wt% of plastic waste (see Table 1), so 8227 kg/h of the carbon dioxide emitted has a non-renewable origin. These emissions are imputed to 13,586 kg of products; with an impact of 0.61 kg CO2eq /kg Biofuel . It should be noted that emissions could be reduced if the CO 2 captured from amine scrubbing were stored or transformed in different products [61].
In Spain, 0.2 t of carbon dioxide are emitted for each MWh demanded [62]. Thus, 0.25 kg CO2eq /kg Biofuel are caused by the usage of 16.7 MW of electricity. Taking into account a default value of 0.085 kg CO2eq /(t·km) [63] for truck transportation (transport of the waste to the production plant and the fuel to the airports) and a distance of 200 km, transportation produces 0.02 kg CO2eq /kg Biofuel . The sum of CO 2 equivalent emissions is 0.9 kg CO2eq /kg Biofuel ; which is considerably lower than in the case of conventional fuels [64], around 3 kg CO2eq /kg Fuel . (production) and 0.9 kg CO2eq /kg Fuel (consumption).
The low heating value of biojet fuel produced is 43.1 MJ/kg, therefore, around 20 g CO2eq /MJ are emitted; complying with the 65% mandatory reduction stated in the European Directive 2018/2001 [65].
Biojet fuel properties are estimated using the analysis tools offered by Aspen HYSYS, to check whether it meets the requirements specified in ASTM D7566-20c standard [66]. Results are shown in Table 13: All estimated parameters comply with the standard. Other parameters (not estimated) could also comply, due to the nature of the fuel produced. Aromatic and naphthalene content is negligible, since the biojet fuel produced contains straight and branched chain hydrocarbons. Moreover, there is no sulphur (and, therefore, no mercaptan content), which together with water, acids, and carbon dioxide, is one of the main causes of corrosion and high electrical conductivity [67]. Distillation curve temperatures should also be met, since the hydrocarbons that make up the biojet fuel match those that make up conventional paraffin (see Figure 4). Finally, the smoke point requirement should not cause problems, since paraffins and isoparaffins (with higher H/C ratio than other hydrocarbons) reduce the amount of smoke produced [67].
Therefore, the estimated properties of the produced biojet fuel fulfill the standard. Experimental tests should be performed to ensure that it could be blended up to 50% with conventional jet fuel [66].

Economic Analysis
Fixed and variable costs are taken into account to estimate the minimum selling price of the biojet fuel produced and assess if the pinch analysis and the use of Rankine cycles makes the process economically more attractive.

Cost of Equipment
The cost the equipment is calculated using equations of Table 14: Table 14. Cost of equipment [46].
• Plastic and biomass feed hoppers and mixers: They are calculated as vessels. A length to diameter (L/D) ratio of 5 gives the minimum cost. Applying Equation S3 (see SI L) and a residence time of 15 min, the hoppers volume can be calculated, which is oversized by 100%: A total of 3 mm of corrosion allowance (C.A.), carbon steel as construction material, and a design pressure of 3.5 kg/cm 2 are used for the following calculations, as lean amine and raw materials are assumed to be non-corrosive: where: D is the diameter of the vessel (m); V, volume (m 3 ); L, lenght (m); t, thickness (mm); C.A., corrosion allowance (mm); P D , the design pressure (kg/cm 2 ); St, a dimensionless factor depending on the material, pressure and temperature conditions. This value is set to 1055 when using carbon steel and moderate conditions; E, modulus of elasticity of carbon steel, which is 0.85 GPa; W, weight of the equipment (t); X, a complexity factor, considered 2 for simple equipment; • Flash chambers: Thickness, mass, and cost of flash chambers are calculated using the following equations: Cost 2007 = 10, 000 + 29·m 0.85 , where data such as shell thicknesses (t Shell ), corrosion allowances, lengths, diameters, and densities of the materials have been obtained from Aspen HYSYS, using a liquid volume in the vessels of 50%.

• Distillation columns:
Their cost is calculated as the sum of vessel and trays cost: Cost Trays 2007 = 0.92· 290 + 550·D 1.9 ·N Trays .
Sections and heights are sized and trays are selected using Aspen HYSYS, so the fluid dynamics in each plate and all parameters are in a normal operating range. Columns are divided into sections and the cost is calculated as the sum of each individual section.
In addition, the installed cost of PSAs for oxygen (to be fed to F-101) and hydrogen purification (V-202A/B), and the R-301 and R-401 reactors, is calculated using the exponent (exp) method and their respective installation factors: where the base capacities (Cap Base ) are the following [28]: • For oxygen separation PSA, the air flow rate obtained (t/h); • For hydrogen separation PSA, the flow rate of hydrogen separated (kg/h); • For the Fischer-Tropsch reactor, the feed flow rate (kmol/h); • For the upgrading reactor, the feed flow rate (kg/h).
The capacity (Cap) of the equipment and costs are shown in Table 15: Table 15. Cost of PSAs (V-202A/B, see Figure 1) and reactors (R-301/R-401, see Figure 2). For the rest of the equipment, the installation factor (F Installation ) is unique and is calculated from different factors, such as material (f M , carbon steel), equipment placement (f ER ), piping cost (f P ), instrumentation and control (f I ), electricity (f EI ), civil works (f C ), buildings and structures (f S ), and thermal insulation and painting (f L ), whose values are 1, 0.3, 0.8, 0.3, 0.2, 0.3, 0.2, and 0.1 respectively [46].

Equipment
To obtain the cost of installed equipment, the following equation is applied: The equipment costs are calculated and updated using the following expression, where the chemical plant cost index (CEPCI) of the year 2007 (former) is 525.4 and for 2018 (current) is 603.1 [45]: Table 16 lists the cost of all equipment: Costs for cooling water (S806), process water (S001), fuel oil (S629), and electricity have been obtained for the year 2015, while that for biomass (S101) has been obtained for the year 2007, so they must be updated using Equation (47); knowing that the CEPCI for the year 2015 was 556.8 [68]. The material balance for the calculation of the fresh process water flow rate is calculated in Supplementary Information Section I.

Minimum Selling Price
The percentage method (see Table 18) is used to obtain a reference value for the fixed and working capital, based on the total cost of equipment [46]: where [46]: F.C. = ISBL + Ut. + O.S. + S.C. + U.C., W.C. = 0.2·C.I..
Engineering and construction works are assumed to last 3 years, in which the fixed capital is invested (20% first year, 50% second year, 30% last year), from years 0 to 2. Plant lifespan is 20 years, from years 3 to 22 and inflation is considered to be 4% per year. Thus, fixed capital is updated: The working capital is invested at the end of the project in the third year at the updated cost, and that amount is recovered in year 22.
During plant operation years, there is no investment, only the profits or losses from the sale and purchase of products and raw materials. In this scenario, 45 employees will be paid an average of 1800 EUR per worker [70], 1,134,000 EUR in total. Insurance is assumed to be 2% [46] of the fixed capital, 19,281,762 EUR. Both costs are added to the cost of raw materials to be purchased. These costs (C.) and sales profits (S.) are updated annually: During the first 10 years of plant operation, depreciation (D.) for the fixed capital invested in the project is accounted as: This depreciation has a value of 0 for the remaining years of operation, and allows tax deductions (T.), as these are calculated as 35% of the pre-tax profits (P.T.P.): T. = 0.35·P.T.P..
Profit after tax (P.A.T.) is calculated as the difference between profit before tax and tax. Finally, funds generated (F.G.) and cash flow (C.F.) are calculated: Results are shown in Supplementary Information Table S10. From the cash flow obtained for each year, a sensitivity analysis is performed by taking each year's cash flow and updating it for a chosen rate of return (RoR): The sum of cash flows gives the net present value (NPV) of the plant. Bioproducts selling prices need to be higher than conventional petroleum fuels prices for the plant to be profitable. Thus, if an internal rate of return (IRR) of 9% is used, [47], bioproducts MSP are shown in Table 19: These prices are 3.47 times higher than conventional fuels prices (see Table 17), but similar to prices found in other studies for biojet fuel [18,21,28]. EU sets a tax of approximately 25 EUR/t CO2eq [18] to penalise CO 2 emissions, which would increase the price of biojet fuel by around 0.02 EUR/kg, which increases the fossil jet fuel by 0.10 EUR/kg. CO 2 emissions price can significantly increase in the coming years [71], and could help the biojet fuel become economically competitive.
According to Deane et al. [72], this would add near 5 EUR to the ticket's price per passenger of a regular 1000 km flight. Thus, it is worth noting that the use of this biojet fuel will not have a significant impact on the final consumer.

Conclusions
Production of biojet fuel from the cogasification of biomass and plastic (GFT-J process) has been described in detail, considering the design of a heat exchanger network, thermodynamic cycles, cooling tower, and an economic analysis. The following conclusions can be drawn: The addition of plastic to biomass gasification improves the volatility of the mixture and reduces tar and char residues, so the yield of raw materials toward products is 13.06%, where the mass distribution of products is: 8.88% fuel gas; 11.7% petrol; 41.8% biojet fuel; 12.5% diesel, and 25% fuel oil. According to Aspen HYSYS properties estimation, the produced biojet fuel complies with the requirements established by the standard ASTM D7566. Thus, it could be blended up to 50% with traditional jet fuel.
Pinch analysis reveals that steam can be generated using the calorific power of hot streams, so 38.73% of electricity, 41.58% of cooling, and 100% of heating demand can be covered with renewable energy (waste biomass) of the process, and the design of a heat exchanger network reduces the total cost of the GFT process.
Biojet fuel MSP (1.37 EUR/L) is higher than conventional jet fuel price, however, the rise of CO 2 emission taxes could help biojet fuel to be more competitive. The produced biojet fuel emits 20.14 g CO2eq /MJ, a significant reduction compared to 94 g CO2eq /MJ emitted by fossil fuels. Moreover, emissions caused by electricity consumption during biojet fuel production could be further reduced by the change to a more renewable electricity mix.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/en14196035/s1, Figure S1: Simulation of unit 1, Figure S2: Simulation of unit 2, Figure S3: Simulation of unit 3, Figure S4: Simulation of unit 4, Figure S5: Simulation of unit 5, Figure S6: Simulation of unit 6, Figure S7: Simulation of unit 7, Figure S8: Simulation of unit 8, Figure S9: Water treatment plant, Figure S10: Heat exchanger network (unit 9), Table S1: Chemical process stream demands, Table S2: Viscosity, density and thermal conductivity of streams, Table S3: Heat exchangers parameters, Table S4: Power demand of mixers, Table S5: Power demand of process pumps, Table S6: Power demand of the pumps upstream of the reboilers,   Acknowledgments: The authors would like to thank the Fundación Iberdrola for supporting this research under project grant "Smart BioH2".

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following are abbreviations, units, and acronyms used in the article.