Experimental investigations and model-based optimization of CZZ/H-FER 20 bed compositions for the direct synthesis of DME from CO 2 -rich syngas

Experimental kinetics studies and model-based optimization were performed for the direct synthesis of dimethyl ether (DME) in a wide range of process conditions, including various CuO/ZnO/ZrO 2 :H-FER 20 bed compositions. Thus, in order to improve the DME productivity using fluctuating CO/CO 2 gas feeds such as those that could be used in power-to-fuel technologies, a series of experiments was carried out at 30 bar by varying the CO 2 /CO x inlet ratio (0.4 – 0.9), temperature (483 – 513 K), gas hourly space velocity (2.78, 3.57 s − 1 ), and the weight ratio of the catalyst components CuO/ZnO/ZrO 2 and H-FER 20 (from 1 to 70.4). A lumped kinetic model for the methanol synthesis from the literature was extended for the direct DME synthesis. With only eight fitted parameters, the model adequately simulates an extensive array of data points measured with a six-fold parallel reactor setup. The resulting model was applied to optimize the CuO/ZnO/ZrO 2 :H-FER 20 ratio and the obtained optimization results were validated experimentally, confirming the simulated performance enhancement. For a CO 2 -rich feed (CO 2 /CO x = 0.9), the optimum volumetric CuO/ZnO/ZrO 2 :H-FER 20 ratio under kinetic operating conditions is 91.5/8.5 vol%. Extrapolation based on the new model for industrially relevant operating cases show an optimal H-FER 20 volume fraction below 5 vol%, resulting in a CO x conversion of 47% at a DME selectivity of 88%.


Introduction
A future with a CO 2 -neutral carbon economy necessitates the development of mature power-to-fuel technologies. 1,2ynthetic fuels and chemicals produced from sustainable and economically viable hydrogen production can compensate for fluctuating and over-potential power generation from renewable energy sources. 3In addition, adapted synthesis routes can be used to recycle CO 2 yielding CO 2 -neutral or even CO 2 -consuming fuels and chemicals. 4One of the most promising options for a flexible and carbon neutral production of chemical energy carriers is the use of CO 2 -rich syngas and its further conversion using efficient and longterm stable catalysts.Besides other synthetic hydrocarbons, dimethyl ether (DME) is a particularly interesting product due to its promising physical and chemical properties. 5DME is the simplest ether, has no C-C bonds and contains oxygen.Therefore, the exhaust gas after the combustion of DME contains less carbon monoxide, unburned hydrocarbons and soot than, for example, during the combustion of butane or diesel. 6[9][10][11] The industrially applied two-step synthesis requires two individual reactors, one for methanol (MeOH) production, where CO and CO 2 are hydrogenated to MeOH ((R1) and (R2)), linked by the WGS reaction (R3).
CO hydrogenation to MeOH CO + 2H 2 ⇌ CH 3 OH ΔH°298 K = −90.4kJ mol −1 (R1) whereas a minor loading of γ-Al 2 O 3 showed a significant decrease in DME selectivity.A catalyst bed variation using CO 2 -rich syngas was studied by Peinado et al. 19 with CZA/γ-Al 2 O 3 (CO/CO 2 /H 2 : 1/1.9/7.7 − CO 2 /CO x = 0.66) showing increased CO x conversions (up to 24 mol%) for a CZA/γ-Al 2 O 3 weight ratio of 9.However, DME was only converted from CO x with a selectivity of 1%.Increasing the amount of γ-Al 2 O 3 in the bed improved the DME selectivity, but at the expense of CO x conversion.Bae et al. 20 also investigated a bed variation of the CZA/γ-Al 2 O 3 catalyst system for a CO/CO 2 syngas mixture of CO/CO 2 /H 2 : 41/21/38 − CO 2 /CO x = 0.34.Here, CO conversion (29.3 mol%) is increasing up to a CZA : γ-Al 2 O 3 weight ratio of 5, however, accompanied by a decrease in DME selectivity (7.9 mol%).In a recent study, Delgado Otalvaro et al. 21found similar results with a variable CO 2 /CO x feed, where an increasing CZA/γ-Al 2 O 3 ratio promotes CO x conversion and a DME yield at a low CO 2 content in the feed; however, a high CO 2 feed content significantly reduces DME selectivity.
A catalyst bed variation of H-ZSM5 in combination with a fixed CZA catalyst mass 24 in a pure CO 2 syngas feed (CO 2 /CO x = 1.0) revealed a nearly stable CO 2 conversion (26 mol%) while DME selectivity remained constant up to a CZA : HZSM-5 mass ratio of 5.
Arena et al. 26 reported the advantages of CZZ compared to CZA systems in the scope of CO 2 hydrogenation.Extensive studies of a CuO/ZnO/ZrO 2 : H-FER 20 (CZZ/FER) catalyst system revealed that the catalyst system showed promising activity for flexible conversion of both CO-rich and CO 2 -rich syngas to DME, as well as superior activity compared to a commercial CZA system, 27,28 confirming the aforementioned    27,29 suggesting that even at higher CZZ ratios, better CO x conversion and thus higher DME yields are possible.Based on the current state of knowledge, a detailed understanding of the dependency of the process variable conversion and selectivity at large CZZ/ FER ratios is of great interest.Besides, the dependence on the syngas composition is important as both CO and CO 2 can be converted 27 and the transition from WGS (R3) to the rWGS reaction is observable to reveal the respective correlations between the catalyst bed composition and CO/CO 2 hydrogenation ((R1) and (R2)).
To the best of our knowledge, such a broad range of operating conditions has not been studied or modeled yet.Our hypothesis is that important impulses for further process development will be derived from the knowledge of syngas composition dependencies.Herein, we investigate the catalyst bed composition of the system CZZ/FER at different process parameters, i.e. the CO 2 content in the synthesis gas, temperature and space time.To describe and simulate a wide operating range, a kinetic model was developed to calculate an optimized CZZ/FER ratio within the studied operating conditions.The model is also applied for predicting the process performance at industrially relevant operating cases.
For CO 2 hydrogenation to DME, Ren et al. 24 showed that the mixing method of the catalyst components (CZA/HZSM-5) has a minor impact on conversion and the DME yield but significantly influences catalyst stability.Similarly, in a COrich environment as reported by García-Trenco et al., 23 slightly increased conversion is observed with a physically mixed bed compared to a hybrid mixed bed configuration.It has also been reported 30,31 that due to closer contact within a hybrid catalyst, Cu and even Zn ion migration from the MeOH-forming catalyst into the dehydration component may occur, reducing both acidity and the available Cu and Cu/Zn surface areas.Also, increased Cu migration into the zeolite pores was found by Fierro et al. 32 at elevated water vapour concentrations.Accordingly, a higher CO 2 feed content should increase Cu migration due to enhanced water production via the rWGS reaction (R3).In order to reduce these transport effects, the catalytic components used in this work were separately pressed, sieved and subsequently physically mixed, as described afterwards.

Catalyst preparation
The CZZ catalyst was prepared by a continuous coprecipitation method from metal nitrate and sodium bicarbonate solutions resulting in pH 7 using a micro jet mixer. 28The suspension was then aged at 313 K for 2 h.The solids were filtered, dried at 383 K for 16 h and calcined at 623 K with 3 K min −1 for 4 h.Characterization techniques used were described in a previous work. 28A commercial ferrierite-type zeolite H-FER 20 (FER) from Zeolyst International with a Si/Al ratio of 20 was used as a dehydration catalyst.Before use, FER was calcined at 823 K for 4 h in air.

Catalyst characterization
Selected chemical and physical data of the CZZ catalyst are shown in Table 4.The XRD analysis of the CZZ catalyst precursor and the calcined pre-catalyst can be found in chapter S9. † Physico-chemical properties of FER are taken from a study by Kim et al. 33 and shown in Table 5.

Catalytic activity study
Each catalyst component was finely powdered, pressed and sieved into fractions of 250-500 μm, followed by physical mixing of the components in the required mass ratio.The studied catalytic bed configurations are summarized in Table 2. Within this sieve fraction, the bulk density amounts to 882.5 kg m −3 for CZZ and to 415.0 kg m −3 for FER.
The catalyst beds were 100 ± 1 mm in length.The physical mixtures were diluted with silicon carbide (SiC, Hausen Mineraliengroßhandel GmbH) of the same grain size as the catalyst components for isothermal operation.To ensure an adequate grain distribution of the catalytic components CZZ, FER and SiC, the beds were filled as five-fold stacks.The catalyst bed volume was constant at 4.20 ml.At the top and bottom of the beds, additional SiC layers were placed.
Direct DME synthesis was performed in a self-constructed parallel reactor system "MURSS" (multi-reactor-screeningsystem, see Fig. 1 and S2 †) with stainless steel fixed bed reactors (inner diameter: 17.4 mm; length: 600 mm).Each reactor can be independently heated with heating cartridges and brass jaws over a length of 400 mm.A catalyst bed temperature profile is measured using axial thermocouples.Gas supply is controlled using mass flow controllers (Bronkhorst Hi-Tec).Feed gases, carbon monoxide (CO, 99.97 vol%), nitrogen (N 2 , 99.9999 vol%), hydrogen (H 2 , 99.9999 vol%) and a mixture carbon of CO 2 /N 2 (50 : 50 ± 1.0 vol%), were supplied by Air Liquid Germany GmbH.The total inlet gas flow is distributed via a capillary system, with fine flow adjustment by flow meters before each reactor.A volume fraction of approx.5 vol% of the reactor outlet gas flows and the bypass flow is directed to a multi-position valve (Valco), which is connected to a FTIR CX4000 (Gasmet Technologies Oy) equipped with a micro gas chromatograph (Inficon Micro  Reduction of the CZZ catalyst volume fraction was performed at 2 bar with 5 vol% H 2 diluted in N 2 , while temperature was increased from 373 K to 473 K with a ramp of 10 K h −1 , followed by further heating to a final reduction temperature of 493 K with 50/50 vol% H 2 /N 2 at a rate of 10 K h −1 .After another 60 min, the reactor was purged with N 2 , cooled to 483 K and, subsequently the pressure increased to 30 bar.The feed gas composition was H 2 /CO x /N 2 = 45/20/35 vol%, and the considered CO and CO 2 feed concentrations respectively CO 2 /CO x inlet ratios are shown in Table 3.Before kinetic measurements, the catalyst beds were exposed to a run-in period of 75 h time on stream (ToS) by varying the temperature (493, 503, and 513 K), CO 2 /CO x inlet ratio (0.4, 0.7, and 0.9) and gas hourly space velocity (GHSV) (2.78 and 3.57 s −1 ) at 30 bar.These measurements were used for the validation of the optimized catalyst bed ratio.The kinetic measurements were performed at each feed gas composition shown in Table 3, with temperatures between 483 and 513 K and two GHSV values of 2.78 and 3.57 s −1 with regard to the catalyst volume.

Experimental error analysis
Before starting the kinetic investigations, the experimental error of the novel system (Fig. 1) was estimated to verify that the data quality in parallel operation is in a similar range to when using a single PFR system.The error caused by an unequal flow distribution via the capillary system is in the range of ±1.2-1.5%.The complete error analysis and the error estimation can be found in S3. †

Performance indicators
In all experiments, the carbon balance (eqn (S1) †) presented a maximum deviation of ±3%.The CO x conversion is calculated as follows (eqn (1)): where n ˙i is the respective molar flowrate.The role (i.e.reactant or product) of CO and CO 2 on DME and MeOH formation is represented by the CO and CO 2 conversion (eqn (2) and ( 3)): The oxygenate/hydrocarbon-based selectivity (eqn ( 4)) is used to obtain a selectivity independent of the possibly produced CO 2 and CO: Here, ν x corresponds to the number of carbon atoms in each C x O y H z product.Formation of DME and MeOH with respect to each catalyst bed is calculated as volumetric DME and MeOH productivities (eqn ( 5) and ( 6)):

Kinetic model development
Kinetic description of the direct DME synthesis is based on the six-parameter model for the MeOH synthesis from Lacerda de Oliveira Campos et al., 34 which has been extended to include the MeOH dehydration step.
In the MeOH synthesis model, only CO 2 hydrogenation (eqn (7)) and the rWGS reaction (eqn ( 8)) are considered, as theoretical studies suggest that direct CO hydrogenation is not significant at a moderate or high CO 2 content. 35,36The reaction rates of CO 2 hydrogenation (r CO 2 mol s −1 ) and rWGS (r rWGS mol s −1 ) are shown as follows:

Reaction Chemistry & Engineering Paper
where m CZZ is the mass of the CZZ catalyst (kg), k 0,CO 2 hyd and k 0,rWGS are the reaction rate constants, E A,CO 2 hyd and E A,rWGS are the activation energies, ϕ Zn is the zinc coverage on the catalyst surface, θ b is the coverage of free Cu/Zn sites, θ c is the coverage of free sites available only for H 2 and H 2 O, f j is the fugacity of gas component j (bar), and K P,CO 2 hyd. 0(bar −2 ) and K P,rWGS 0 are the global equilibrium constants.The zinc coverage varies depending on the gas composition. 37,38We investigated theoretical approaches for the zinc coverage estimation (Kuld et al. 38 and Ovesen et al. 39 ), but the additional zinc coverage estimation did not improve the simulation results.Still, in our operating region of interest (CO 2 /CO x between 0.40 and 0.90), we assumed that the zinc coverage change is small enough for a simplified constant value to be considered.Therefore, we followed the approach of Lacerda de Oliveira Campos et al. 34 of considering a constant zinc coverage value of 0.50.
The fugacity coefficients are calculated with the Peng-Robinson equation, using binary interaction parameters and other parameters from the literature. 40,41The free sites θ b and θ c are calculated as follows: Here, K 1 (bar −1.5 ) and K 2 (bar −0.5 ) are the adsorption constants.Arvidsson et al. 42 performed density functional theory (DFT) calculations for MeOH dehydration on different zeolites and found that the associative mechanism is dominant at temperatures lower than 573 K.In our work, we assume that the reaction mechanism of MeOH dehydration over FER is also the associative path.This mechanism consists of two elementary reactions: the MeOH adsorption (eqn (11)) and the associative reaction (eqn ( 12)), where Z-H is a free acid site and Z-H-CH 3 OH represents the adsorbed MeOH.Considering the associative reaction (eqn ( 12)) as the rate determining step (RDS), 42 the reaction rate of MeOH dehydration (r Dehyd , mol s −1 ) is calculated in eqn (13).
The mathematical derivation of this reaction rate is provided in the ESI.† Here, m FER is the mass of FER (kg), k Dehyd is the MeOH dehydration rate constant, θ d is the coverage of free zeolite sites, and K P,Dehyd 0 is the global equilibrium constant.It should be noted that the rate of MeOH dehydration over FER can be affected by mass transfer limitations, due to the microporous nature of the zeolites.However, in this model, we assume that there are no significant mass transfer limitations, due to the small particle size range of FER (250-500 μm).The coverage of the free zeolite sites (θ d ) is calculated as follows: where K 3 (bar −1 ) is an adsorption constant.The kinetic model has initially nine parameters to be estimated: the reaction rate constants (k 0,CO 2 hyd , k 0,rWGS , and k 0,Dehyd ), the activation energies (E A,CO 2 hyd , E A,rWGS , and E A,Dehyd ), and the adsorption constants (K 1 , K 2 , and K 3 ).Using our experimental database, it was found that E A,Dehyd is statistically not significant, probably because the MeOH dehydration is at quasiequilibrium for operating points with high amounts of FER.
Therefore, E A,Dehyd is excluded from the model (E A,Dehyd = 0) and eight remaining parameters (k 0,CO 2 hyd , k 0,rWGS , k 0,Dehyd , E A,CO 2 hyd , E A,rWGS , K 1 , K 2 , and K 3 ) are to be estimated.The estimation of the kinetic parameters is done by solving an optimization problem, where the objective function is the minimization of the normalized squared errors of the prediction (y ˆi out,j ) of the carbonaceous compounds (CO, CO 2 , MeOH and DME), the so-called chisquared (χ 2 ) regression method.The normalization with experimental y i out,i squared values is performed to prevent overweighting of high conversion points and underweighting of low conversion points.(15)   Here, N p is the number of data points in the training set, y i out,j is the experimental output mole fraction of gaseous component j at the operating point i, and y ˆi out,j is the simulated output mole fraction of gaseous component j of the operating point i.

Reaction Chemistry & Engineering Paper
In order to figure out the best parameter set and to calculate the confidence intervals, a 5-fold cross validation (CV) method was used. 43,44The 240 experimental points were randomly divided into five groups of 48 points each, and the optimization problem was solved five times, each one with four of the five groups as a training set, resulting in five sets of parameters.The group of parameters with which the model has the lowest χ 2 value for the total 240 points is chosen as the best one.The confidence interval (CI) of each parameter is derived by calculating the standard deviation of each parameter considering the five parameter sets and multiplying by the t-student factor (0.05 two-tail significance and 232 degrees of freedom).The optimization problems are solved with the Matlab function fminsearch (varying the starting values).
The 5-fold cross validation (CV ( 5) ) value is a mean value of the χ 2 five different parameter sets (eqn ( 16)).The closer it is to the χ 2 value of the best fit, the better the model should simulate experiments outside the training set.
The mean squared error (MSE) and the mean error (ME) are also statistical indicators of the model quality.They are calculated for each component j according to eqn ( 17) and (18).

Model validation
The group of parameters with the best fit and the respective CI are shown in Table 6, with all eight parameters being statistically significant.The statistical evaluation of the chosen parameter set is shown in In Fig. 2, parity plots are shown correlating the experiments and simulated results using the best-fit parameter set shown in Table 6.
The simulated concentrations are mostly located inside the ±10% lines of the experimental values, including all H 2 points, all CO points, all CO 2 points, 74% of the MeOH points, and 87% of the DME points.In the parity plot of H 2 O, 86% of the points are within ±20% lines.The experiments seem to be affected by H 2 O concentrations above 1.5 vol% as a result of condensation effects, most likely due to the pressure-affected analysis stage as well as dead volumes or cold spots.Therefore, condensation of mainly water, but also of MeOH at higher product concentrations may occur.In our study, however, only 12.5% of all points are within this concentration range.
In Fig. 3, experimental and simulated DME and MeOH productivities are shown as a function of CZZ volume fraction (ζ CZZ ) for the following operating conditions: CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) at 30 bar, 2.78 s −1 and 483-513 K in steps of 10 K in between (further CO 2 /CO x inlet ratios are shown in the ESI †).
The dashed lines correspond to the simulated values and the dots correspond to the experimentally measured productivities for MeOH (framed dots) and DME (unframed dots), respectively.The model simulates the trends adequately, with slight underestimations of the DME productivity at 513 K. Here, non-isothermal bed temperature can be excluded from axial temperature measurements with maximum ΔT values of 1.1 K (Fig. S13 †).
In Fig. 4, CO and CO 2 conversions are shown, each as a function of ζ CZZ for the following operating conditions: CO 2 / CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) at 30 bar, 2.78

DME productivity vs. CZZ/FER variation
From the results shown in Fig. 3, it can be observed that an increase in temperature leads to a largely uniform increase in the productivity of DME and MeOH, indicating a strong kinetically controlled operating regime, additionally evidenced by the wide gap from the thermodynamic equilibrium (Fig. S10 †).With a ζ CZZ of 0 vol%, no MeOH is formed and the subsequent MeOH dehydration cannot take place and, as a result, no DME is produced.With increasing ζ CZZ , both MeOH and consequently DME are formed.Further increasing ζ CZZ up to ca. 90 vol% enhances MeOH and DME productivities, whereby different temperatures or CO 2 /CO x inlet ratios slightly shift the optimum towards more or less CZZ content (see Fig. 7).This indicates that, for ζ CZZ up to ca. 90 vol%, FER is in excess, which means that the MeOH dehydration is running close to equilibrium.In this operating region, the rate of the reaction system is controlled only by the MeOH synthesis (over CZZ).A further increase of ζ CZZ leads to a marked decrease in DME productivity and, accordingly, to a sharp increase in non-dehydrated MeOH.This suggests that FER is not in excess anymore, which means that neither the MeOH synthesis nor the MeOH dehydration are in equilibrium, and the rates of the reaction system are determined by both CZZ and FER masses.
High MeOH dehydration activity, as shown here with a FER catalyst component, was similarly reported by Peláez et al. 17 with a variation of the CZA/γ-Al 2 O 3 catalyst system.In their study, a highest DME yield was achieved with a CO 2 -free syngas at a CZA/γ-Al 2 O 3 ratio of 92.5 : 7.5 wt%, being close to the experimental CZZ/FER optimum of 95 : 5 wt% reported in this work.In this context, it should be noted that changing the CO x source from CO to CO 2 increases water production (R3) typically leading to a rapid decrease in dehydration performance when γ-Al 2 O 3 is used. 19,21,27,29In comparison, when the dehydration of MeOH takes place over a solid acid component with fewer Lewis acid sites (known to adsorb water) such as FER, 15,27 HZSM5 (ref.45-47) or SAPO, 48 high DME productivity can be maintained even at an increased CO 2 content in syngas or in the presence of increased water content.In our study, increasing the CO 2 feed content leads to decreased DME productivity (Fig. 3a-c), being more evident at higher temperatures.Since the thermodynamic influence can be largely ruled out in this operating range (Fig. S10 †), it is a kinetically controlled phenomenon.A temperaturecontrolled decrease in DME and MeOH productivities by increasing the CO 2 content is explained by Sahibzada et al. 49 arguing that the lower concentrations of products under differential conditions (low conversion) are less inhibiting MeOH production from CO 2 .

CO and CO 2 conversion vs. CZZ/FER variation
For a CO 2 /CO x inlet ratio of 0.6 (Fig. 4a), an increase in ζ CZZ leads to an increase in CO 2 conversion up to a ζ CZZ of approx.30 vol%.At ζ CZZ values below 32 vol%, the CO 2 conversion is higher than the CO conversion.Below a ζ CZZ of approx.23 vol%, CO is produced in small amounts.A further increase of ζ CZZ has a beneficial effect on CO conversion up to a ζ CZZ value of about 95 vol%.Analogously, the CO 2 conversion Fig. 3 Comparison of the experimental DME (dots) and MeOH (framed dots) productivities and the simulated values (dashed lines) for a CZZ volume fraction variation at 30 bar, 2.78 s −1 , and 483-513 K and the three CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) with 20 vol% total CO x in the feed.
Fig. 4 Comparison of the experimental CO (dots) and CO 2 (framed dots) conversion and the simulated values (dashed lines) for a CZZ volume fraction variation at 30 bar, 2.78 s −1 , and 483-513 K and the three CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) with 20 vol% total CO x in the feed.
decreases slightly.This effect is mainly caused by the fact that the increased MeOH formation leads to increased MeOH dehydration (R4) and thus water formation.The produced water allows higher CO conversion via the accelerated WGS reaction (R3) to CO 2 and H 2 , partially regenerating the consumed CO 2 .At ζ CZZ above approx.95 vol%, less DME and consequently water are formed, resulting in a reduced WGS reaction rate and thus decreasing CO conversion, leading to a slightly increased CO 2 conversion.
An increase in the CO 2 /CO x inlet ratio to 0.75 shows a similar pattern for CO and CO 2 conversion (Fig. 4b) with shifts toward more CO formation at lower ζ CZZ values and less CO conversion at higher ζ CZZ values, these effects being even more pronounced at higher reaction temperatures.At ζ CZZ values above ca.35%, the CO 2 conversion shows almost constant values at the respective measured temperatures, showing that the WGS reaction rate is lowered by less CO in the syngas feed (R3).
Further increasing the CO 2 /CO x inlet ratio to 0.9 (Fig. 4c) leads to a pattern where CO is unexceptionally formed over the entire CZZ/FER variation range via the rWGS reaction (R3).The CO 2 conversion is favoured by an increased ζ CZZ value, but to a smaller extent than the CO conversions at lower CO 2 /CO x inlet ratios (Fig. 4a and b).A further potential increase in CO 2 conversion with an optimized CZZ/FER ratio is probably limited by increased water formation from the rWGS reaction occupying the active sites (Cu/Zn) at the CZZ catalyst. 35,36n summary, the increase of CO 2 conversion with optimized ζ CZZ is limited due to water formation by CO 2 hydrogenation and MeOH dehydration, accelerating the WGS reaction and preventing the enhancement of CO 2 conversion, whereby this effect being more pronounced at increased temperature.An increase of CO 2 conversion at higher ζ CZZ values should be possible by in situ water removal, e.g. with a membrane reactor system 50,51 or by doping the CZZ catalyst aiming at stronger binding of formate to the Cu sites, which could inhibit the rWGS reaction and thus accelerate CO 2 hydrogenation on the Cu/Zn sites. 36OH and DME selectivity vs. CZZ/FER variation In Fig. 5  CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) at 30 bar, 2.78 s −1 , and 483-513 K, (further CO 2 /CO x inlet ratios are shown in Fig. 8).A similar pattern can be seen for each CO 2 /CO x inlet ratio investigated: high and constant DME selectivities for ζ CZZ up to ca. 90 vol% and a sharp decrease of DME when ζ CZZ is increased above 90 vol%.
As MeOH dehydration is slightly exothermic (R4), an increase in temperature should directly decrease DME  selectivity.Nevertheless, temperature is affecting DME selectivity in a more pronounced way by influencing the rWGS reaction (R3) and CO 2 hydrogenation (R2) and thus the water concentration/selectivity (see Fig. S9 †).The patterns of the respective water selectivities (eqn (S15) †) clearly display that the changes of DME selectivities follow the inverse trend of the water selectivities.

Optimized CZZ/FER ratio
The validated kinetic model was used to find the optimum FER amount with respect to the highest DME productivity.For a CO 2 -rich feed (CO 2 /CO x = 0.9), 30 bar, and 2.78 s −1 , an optimal ζ CZZ value of 91.5 vol% was found, with a temperature variation (493-513 K) having only marginal influence under these operating conditions (see Fig. 7c).This simulated optimum ζ CZZ value was experimentally validated, following the same run-in period as the catalyst beds in Table 2. Fig. 6 shows the DME productivities of the run-in period (unframed dots) and the experimental validation of the model-based optimized catalyst bed (framed dots).
Compared to the catalyst bed with the highest experimental DME productivity (ζ CZZ = 90 vol%), the simulated optimal value (ζ CZZ = 91.5 vol%) achieved a slight relative increase in DME productivity at all measured temperatures: 493 K (2.9%), 503 K (3.2%), and 513 K (3.8%).Since the present model adequately simulated the experiments, it is applied to investigate the optimum FER volume fraction (ζ FER ) within a broader GHSV range, including industrially relevant operating conditions.In Fig. 7, the optimized FER amounts are shown at 30 bar and CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) and under a variation of temperature and GHSV.The FER optimum value is mostly affected by the GHSV, with less FER being necessary for lower GHSV (achieving ζ FER values lower than 5 vol%).Increasing the temperature slightly increases the FER optimum amount, while increasing the CO 2 content in the feed marginally decreases the FER optimum amount.Since lower GHSV leads to increased CO x conversion and thus enhanced DME production, higher water accumulation is expected.The lower FER amount required for an optimized DME productivity at lower GHSV is therefore a hint that the activity of FER is less affected by water accumulation than the activity of CZZ.This hypothesis is supported by the fact that the optimum FER amount is slightly lower for higher CO 2 /CO x inlet ratios (cf.Fig. 7a-c), which also increase water formation.
In Fig. 8, CO x conversion is plotted against the GHSV under the same conditions considered in Fig. 7.The GHSV variation shows a clearly kinetically dominated range between 1 and 4 s −1 , while an increase in the thermodynamic control is expected as GHSV decreases, with CO x conversion achieving its equilibrium at 513 K.
By using process parameters of 513 K and a CO 2 /CO x inlet ratio of 0.6, a GHSV of 0.4 s −1 leads to 95.8% of the thermodynamically possible CO x conversion (X CO x ,equil = 31.9%)at a DME selectivity of 90.1%, while CO 2 /CO x = 0.75 affords 95.4% of the possible conversion (X CO x ,equil = 24.4%) at a DME selectivity of 87.4%, and CO 2 /CO x = 0.9 affords 94.8% of the possible conversion (X CO x ,equil = 17.6%) at a DME selectivity of 83.4%.A GHSV of 0.4 s −1 corresponds to a seven times lower volume flow, chosen in the kinetic measurements presented here.Further simulations were performed at 60 bar and without N 2 dilution (H 2 /CO/CO 2 in the feed = 70/12/18 vol%), i.e. closer to industrial conditions.Although this requires a greater extrapolation of our model, the preliminary and interesting finding is that in this case probably even lower amounts of FER (less than 3 vol%) already ensure efficient DME production.From simulation at 513 K, such an optimized CZZ/FER bed ratio enables a CO x conversion of more than 47% at a GHSV below 0.4 s −1 (Fig. S11 and S12 †), which corresponds to 90.3% of the thermodynamically possible CO x conversion (X CO x ,equil = 52.2%)at a DME selectivity of 88.9%.Further experimental validation in this operating range would be valuable to back up these promising simulated results.
Other CZZ catalysts with different properties (cf.Table 4) as well as ferrierites (cf.Table 5) with other Si/Al ratios would require a refit of the parameter set (Table 6) and possibly a new evaluation of the result interpretation.

Conclusions
Kinetic experiments for direct DME synthesis were carried out under various operating conditions, with particular emphasis on the CZZ/FER bed composition.A new kinetic model was developed and validated with a vast array of experimental data.It was applied in model-based optimization to determine the optimal FER volume fraction with respect to DME productivity under the variation of temperature, GHSV and CO 2 /CO x inlet ratios.An optimum CZZ/FER catalyst bed ratio (91.5/8.5 vol%) for a CO 2 -rich feed was calculated with the model and then validated experimentally verifying that the results obtained from the optimization are accurate.Extrapolations of the model to process conditions closer to industrially relevant conditions showed that the optimal FER volume fraction actually decreases at lower GHSV.From these findings, the necessity for additional validation beyond the scope of the present study is derived, which is the subject of current investigations.Our experiments, together with the results from modelling also underline that water formation from MeOH dehydration accelerates the WGS reaction regenerating CO 2 .This at first prevents a significant increase of CO 2 conversion with an optimized catalyst bed, which is why alternative reactor concepts where water can be separated in situ, e.g.membrane reactors, appear promising.In the outlook for a further possibility to increase CO 2 conversion, a catalyst modification (e.g., with promoters) could be purposeful, which strengthens binding of the intermediate formate to the Cu sites, which otherwise inhibit the rWGS reaction and thus accelerate CO 2 hydrogenation
This article is licensed under a Creative Commons Attribution 3.0 Unported Licence.

GC
Fusion).The GC is equipped with two thermal conductivity (TCD) detectors connected to RT-Molsieve 5A, 0.25 mm (10 m) and RT-Q-Bond, 0.25 mm (12 m) columns.Each operating point was held for 200 min and the catalyst beds were measured consecutively.After finishing the variation loops of CO 2 /CO x values for each temperature, the reactor was purged with N 2 for one hour (the general sequence of the process parameters is shown in Fig.S1 †).

Fig. 1 A
Fig. 1 A schematic flow chart of the parallel reactor system.

s − 1 ,
and 483-513 K (further CO 2 /CO x inlet ratios are shown in Fig.S7 †).The CO and CO 2 conversions are adequately predicted by the kinetic model.At a CO 2 /CO x feed of 0.90, there are some deviations from the experimentally measured CO conversion, probably due to the fact that the initial CO concentration is low, which leads to an amplification of the small deviations in the final CO concentration when calculating the corresponding CO conversion.

Fig. 2
Fig. 2 Parity plots of CO, CO 2 , H 2 , MeOH, DME and H 2 O comparing the 240 simulated results with the respective experimental results.

Fig. 5
Fig. 5 Comparison of the experimental DME (dots) and MeOH (framed dots) selectivity and the simulated values (dashed lines) for a CZZ volume fraction variation at 30 bar, 2.78 s −1 , and 483-513 K and the three CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c).

Fig. 6
Fig.6Comparison of the experimental DME productivity in the run-in period and the model-optimized CZZ/FER volume ratio of 91.5/8.5 at CO 2 /CO x = 0.9, 30 bar, 2.78 s −1 and three different temperatures.

Fig. 7
Fig. 7 Contour plots of the optimal FER volume fraction in view of maximum DME productivity under the variation of temperature and GHSV at 30 bar and CO 2 /CO x inlet ratios: 0.6 (a), 0.75 (b) and 0.9 (c).

Fig. 8
Fig.8Simulated CO x conversion under the variation of GHSV with the optimized CZZ/FER catalyst beds at 30 bar and 483-513 K for the CO 2 /CO x inlet ratios of 0.6 (a), 0.75 (b) and 0.9 (c) and the respective thermodynamic equilibrium at 513 K.

Table 1
A literature summary of optimal catalyst bed compositions in the direct DME synthesis activity of Zr-promoted Cu/ZnO systems in CO 2 -rich environments.A comparison of the dehydration properties of γ-Al 2 O 3 and FER at different CO 2 /CO x syngas ratios shows that DME selectivity is less affected for both in CO-rich and CO 2 -rich syngas when FER is used, a Based on DME yields.bBased on CO x conversion.favorable

Table 3
CO 2 /CO x inlet ratios and respective feed gas compositions used in the direct DME synthesis

Table 4
Selected properties of the CZZ catalyst component

Table 5
33T surface areas and total acidity properties of the acid dehydration catalyst FER at low-temperature (LT) and high-temperature (HT) taken from the study by Kim et al.33CatalystS BET /m 2 g −1 NH 3 -TPD peak positon/K Acid amount/mmol NH 3 per g cat

Table 7
, with a total χ 2 of 3.406 for the best fit and a CV (5) of 3.439.The CV(5)value is close to the best χ 2 , which suggests that the model can adequately simulate experiments outside the training set.The MSE and ME of all optimized species are significantly low, with MeOH and DME showing mean errors of only 7% and 5%, respectively.The MSE and ME values of the training and the validation sets are significantly close, which also points to the good performance of the model outside the training region.

Table 6
Fitted parameters of k 0,i , E a,i and K i within the training set

Table 7
Statistical evaluation of the best parameter set