Next Article in Journal
Microwave-Assisted Synthesis of Tri-Substituted 1,3,5-Triazines from Metformin Using Benzotriazole Chemistry
Next Article in Special Issue
Co-Hydroprocessing of Fossil Middle Distillate and Bio-Derived Durene-Rich Heavy Ends under Hydrotreating Conditions
Previous Article in Journal
Advances for Biorefineries: Glycerol Hydrogenolysis to 1,3-Propylene Glycol
Previous Article in Special Issue
Editorial for Special Issue “Hydrogen Production and Storage”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Modeling for the Gas-Phase Hydrogenation of the LOHC γ-Butyrolactone–1,4-Butanediol on a Copper-Zinc Catalyst

1
CEA, Liten, DTCH, Laboratoire Réacteurs et Procédés (LRP), Univ Grenoble Alpes, F-38000 Grenoble, France
2
Univ Lyon, Université Claude Bernard Lyon 1, CNRS, LAGEPP UMR 5007, 43 Boulevard du 11 Novembre 1918, F-69100 Villeurbanne, France
*
Authors to whom correspondence should be addressed.
Reactions 2022, 3(4), 499-515; https://doi.org/10.3390/reactions3040033
Submission received: 20 July 2022 / Revised: 25 August 2022 / Accepted: 8 September 2022 / Published: 26 September 2022
(This article belongs to the Special Issue Hydrogen Production and Storage, 2nd Edition)

Abstract

:
Liquid organic hydrogen carriers (LOHCs) are an interesting alternative for hydrogen storage as the method is based on the reversibility of hydrogenation and dehydrogenation reactions to produce liquid and safe components at room temperature. As hydrogen storage involves a large amount of hydrogen and pure compounds, the design of a three-phase reactor requires the study of gas and liquid-phase kinetics. The gas-phase hydrogenation kinetics of LOHC γ-butyrolactone/1,4-butanediol on a copper-zinc catalyst are investigated here. The experiments were performed with data, taken from the literature, in the temperature and pressure ranges 200–240 °C and 25–35 bar, respectively, for a H2/γ-butyrolactone molar ratio at the reactor inlet of about 90. The best kinetic law takes into account the thermodynamic chemical equilibrium, is based on the associative hydrogen adsorption and is able to simulate temperature and pressure effects. For this model, the confidence intervals are at most 28% for the pre-exponential factors and 4% for the activation energies. Finally, this model will be included in a larger reactor model in order to evaluate the selectivity of the reactions, which may differ depending on whether the reaction takes place in the liquid or gas phase.

1. Introduction

The threat of climate change as well as geopolitical tensions are gradually reducing or even prohibiting the use of fossil fuels. Hydrogen, produced from renewable energy sources, is widely studied as a replacement [1]. Among the hydrogen storage technologies, the most mature are compression or liquefaction and one of the most studied is physisorption. However, they have drawbacks regarding safety, due to high compression pressures and the risk of leakage and evaporation of flammable hydrogen [2,3,4]. With regard to chemical storage of hydrogen in carrier molecules, this technology can involve a wide variety of chemical reactions and storage solutions. The most common ways involve power-to-X systems, by CO2 and/or N2 hydrogenation to fuels capable of being stored under reasonable conditions over extended periods [1]. However, these systems often lack in reversibility because hydrogen is hardly retrieved [1,5]. Another solution consists in binding hydrogen to metals under a solid state [6]. Anyway, the most promising solutions in terms of hydrogen gravimetric density are currently not reversible and the reversible ones are limited in terms of hydrogen content. Another way is the use of a liquid organic hydrogen carrier (LOHC), which involves reversible reactions, is easy to transport and stable over long periods [1,4]. Indeed, a LOHC consists in a couple of hydrogen-lean|hydrogen-rich molecules, which are respectively catalytically hydrogenated and dehydrogenated into the other molecule in order to store and release hydrogen. These organic liquids are fluids sharing similar properties with hydrocarbon fuels or common chemicals. Their stability at ambient temperature and atmospheric pressure in the liquid phase makes them compatible with the existing fuel infrastructures and thus LOHCs can be easily transported in tankers, trucks and/or pipelines [1,4,7].
Multiple LOHC systems have already been identified. The main ones highlighted in the literature are aromatics|cycloalkanes, such as toluene|methylcyclohexane (TOL|MCH), dibenzyltoluene|perhydro-dibenzyltoluene (H0-DBT|H18-DBT), benzyltoluene|perhydro-benzyltoluene (H0-BT|H12-BT) and couples based on N-heterocycles, such as N-ethylcarbazole|perhydro-N-ethylcarbazole (H0-NEC|H12-NEC). The first cycloalkane-based LOHCs were developed because of their good hydrogen capacities (5.8–6.2%), their wide availabilities and the amount of knowledge on the processes of hydrogenation and dehydrogenation. However, reaction enthalpies, which dictate how much energy is lost during the storage process, are high for these LOHCs at around 65–69 kJ·molH2−1, although NEC LOHC has a slightly lower reaction enthalpy (52 kJ·molH2−1). Safety and toxicity profiles have also been improved from the early light aromatics LOHCs to NEC, BT and DBT LOHCs, due to their higher flash points and lesser carcinogenic properties. The liquid state criterion, which is a key point for the process designs, is hardly met for some LOHCs at some conditions. For example, H12-NEC is solid at ambient conditions. Moreover, DBT LOHCs have reaction intermediates with high viscosities. NEC has also reversibility issues with high degradation when subject to harsh dehydrogenation temperatures. Finally, the use of expensive noble metal catalysts and relatively severe reaction conditions hinder the economical or environmental viability of the aromatic-based LOHCs [1,3,7,8,9,10,11,12]. CEA has recently identified a bioavailable LOHC: γ-butyrolactone|1,4-butanediol (GBL|BDO) [13]. The latter, described for the first time by Onoda et al. [14], has low reaction enthalpies of 42 kJ·molH2−1 in liquid phase and 31 kJ·molH2−1 in gas phase, thus making it the most energy efficient. Its lower cycle operating temperature range (160–220 °C), compared to that of cycloalkane-based LOHCs _(300–350 °C), makes its process less energy consuming. BDO is partly produced from biomass, by hydrogenating bio-based succinic acid into GBL or by producing BDO directly from a fermentation of sugars over microorganisms [14,15,16,17], although most of it is currently produced from acetylene and formaldehyde [15,16,17]. It has a H2 capacity of 4.5 wt%, which is lower than that of the other LOHCs. The melting point of pure BDO at 20 °C could limit its use in colder weather or regions [16], but as for the stoichiometric ratio in feed, it remains liquid at reaction operating temperature and pressure, so reaction product separation is easier. Although GBL and BDO have good safety and toxicity profiles, GBL is known as a drug precursor because of its psychotropic properties and chemical resemblance to γ-hydroxybutyric acid (GHB) [18]. They are commonly used as solvents and monomers and their uses are allowed in industrial processes in most countries. See Table S1 in Supplementary Materials for LOHC comparisons [1,16,19,20,21,22].
For all LOHCs, except for TOL|MCH, hydrogenation and dehydrogenation are performed in three-phase reactors. Although LOHCs are usually compared with respect to their mass H2 capacities, the volume H2 capacity is therefore also useful for the design of suitable reactors as it indicates the volume of hydrogen involved during the reactions. For all LOHCs, the ratio H2/LOHC is significant (613–856 LH2·L−1 LOHC) and could act as a barrier to liquid–solid transfer by forming dry zones on the surface of the catalyst pellets. Thus, the knowledge of the reaction kinetics in both the liquid and gas phases is necessary to accurately design a three-phase reactor. As the main quality of LOHCs is the perfect reversibility of their reactions, the formation of co-products in both phases and reactions must also be assessed. This study is therefore a part of a larger study on GBL|BDO. Here, the paper is focused on the kinetics of the hydrogenation in the gas phase. The kinetics of the hydrogenation in the liquid phase will be the subject of another article.
The reaction pathway describing the equilibrium between GBL and BDO has been debated within the literature. Most have studied the reaction on copper-based catalysts in the liquid phase with temperature and pressure ranging from 200 to 250 °C and from 30 to 90 bar [14,15,23,24], respectively. Two parallel reactions, resulting in a dehydration towards tetrahydrofuran (THF) and n-Butanol (BuOH), occur at high temperature and are catalyzed by acid sites present in the catalyst support [23,24,25,26,27]; however, their mechanisms are unclear. Reaction intermediates, such as 2-Hydroxytetrahydrofuran, have been proposed as they could be involved in the formation of secondary products [19,24,27]. Ichikawa et al. proposed a model taking into account its formation without measuring it [27]. Chaudhari et al. proposed a mechanism where both GBL and BDO can lead up to the formation of THF and BuOH [23]. However, no article seems to indicate that the latter are formed from intermediates, despites hydrogen imbalances between the reactants and the secondary products (except for the dehydration of BDO into THF). For instance, the experimental data in Schlander’s Ph.D. thesis [25], performed in the gas phase, seem to indicate that GBL concentrations decrease faster than those of BDO as the secondary products build up in the mixture. As for Chaudhari et al., they contend THF is likely to be formed from the hydrogenation of GBL, rather than the dehydration of BDO [23]; we have also used Schlander’s reaction pathway, shown in Figure 1, for the kinetic modeling.
Kinetic models for the hydrogenation of GBL usually involve a simple power law model such as that proposed by Schlander [15,23,25]. However, it (Model 1 in §2.2.1) proved unable to correctly predict the experimentally described pressure effect. Hermann et al. [24] proposed a Langmuir–Hinshelwood type model because an intermediate is thought to be formed during the hydrogenation of GBL, but this approach was first ruled out because of the number of parameters needed to be estimated, giving the model too many degrees of liberty and thus making it hard to extrapolate. This is why we decided to revise the kinetic study of the GBL gas-phase hydrogenation. The aim was to build a kinetic model, able to predict the effect of temperature and pressure, to be included in a larger reactor model including all phases involved. As the data from Schlander’s works [25,26] were obtained under operating conditions that can mimic the composition of the gas phase in a trickle bed, the published data are used here to determine the appropriate kinetic law.
The GBL hydrogenation was done in gas phase on a CuZnO catalyst in an isothermal tubular reactor. The fixed bed was loaded on a reactor length of 80 mm with crushed catalyst and glass of 400 µm-diameter particles. The catalyst loading was about 8 g. The molar ratio H2/GBL at the reactor entrance was about 90. Temperature and pressure ranges from 200 to 240 °C and from 25 to 35 bar, respectively, were screened. The experimental data presented in Figure 2 are plotted as a function of the residence time tV, which is the ratio between the reactor volume VR and the total volume flow rate Q. The bimetallic catalyst was prepared by co-precipitation and contained 15 mol% Cu and 85 mol% ZnO [28]. See Section S2 in Supplementary Materials, and references [25,26,28] for further details. Regarding the potential transfer limitation, Schlander showed that the apparent kinetics could not be explained either by internal or external mass transfer (Section S3 in Supplementary Materials).
In this paper, based on these previous bibliographic results, a reactor model was first described and then several kinetic laws were tested and kinetic parameters estimated in order to predict the effect of temperature and pressure on the gas-phase hydrogenation of GBL into BDO. The kinetic laws were established in interaction with the thermodynamics to take into account the reversibility of the reactions.

2. Materials and Methods

2.1. Reactor Model

To establish mass balance equations suitable for simulation of the experiments, some criteria were checked. GBL was highly diluted by hydrogen in the feed. As the H2/GBL molar ratio at the reactor inlet was about 90, the GBL conversion induced at most a variation of about 1.8% of the volumetric flow rate along the catalytic bed and a negligible change in the gas phase composition for all experiments considered in the kinetics. Consequently, the variation of the gas velocity along the reactor could be neglected as well as the possible temperature gradients induced by the reaction. So, the catalytic bed was considered as isothermal. Moreover, as the reactor diameter (dt) to particle diameter (dp) ratio was greater than 10 (Equation (1)), then wall effects causing a radial velocity profile resulting in deviations from plug flow could be neglected. In addition, as the reactor length (LB) to particle diameter ratio was greater than 50 (Equation (2)), the axial dispersion resulting in consequences on the residence time distribution could also be neglected.
d t d p = 3.1 · 10 2 > 10
L B d p = 1.7 · 10 3 > 50
Consequently, an isothermal plug flow model was applied. Steady state was supposed to be reached. The material balance followed the evolution of a compound j within an elementary volume of catalyst V c a t a :
d F j d V c a t a = v i j   r P l ¯
where i is the reaction index and j is the compound index. The apparent kinetic rate is:
r p l ¯ = η i r s i
As the diffusion has been estimated as a non-limiting step, the efficiency factor η i is considered equal to 1 (See Section S3 in Supplementary Materials). Capillary condensation has also been assessed in Section S3 in Supplementary Materials.
The experimental data used the scale tV, which is the residence time. Thus
d x j d t V = r j
with xj is the molar fraction of condensable compound j equal to the ratio between the compound partial pressure Pj and the partial pressure of GBL at the reactor entrance PGBL,0 and rj the production rates of compound j.

2.2. Kinetic Modeling

Model 1 and Model 2 were built upon power laws with Model 2 including equilibrium constants. Model 3 and Model 4 included in addition an inhibition parameter on the hydrogen pressure. All the models were based on the reaction scheme described in Figure 1. All the models presented in this article (whether kinetic or thermodynamic) were constructed on the basis of the ideal gas hypothesis (see Section S4 in Supplementary Materials).

2.2.1. Model 1

Model 1 is a simple power law model first investigated in modeling the hydrogenation of GBL into BDO. It is the model proposed by J.H. Schlander in his thesis [25]. To avoid the interpretation bias due to the numerical methods, this model has been assessed here with the same method as used for the other models. That explains the difference between the values of the kinetic parameters estimated in this paper and given in the thesis.
Based on the reactions given in Table 1, the production rates of all active compounds are as follows:
r G B L = k B D O x B D O k G B L x G B L k T H F x G B L k B u O H x G B L
r B D O = k G B L x G B L k B D O x B D O
r T H F = k T H F x G B L
r B u O H = k B u O H x G B L
where kj are the kinetic parameters.

2.2.2. Model 2

In order to study the effect of pressure into the model, the terms of the production rates were modified to introduce the hydrogen pressure and the chemical equilibrium. kTHF and kBuOH were grouped under a single parameter kSP to better fit the single kinetic following of both side reactions. The reaction mechanism for this power law model is given in Table 2. Moreover, because of the large excess of hydrogen, the total pressure within the reactor was assumed to be the partial pressure of hydrogen. The activity of hydrogen for side reactions was neglected because of its low effect on kinetics.
r G B L = k 1 K x B D O x G B L P H 2 2 k S P x G B L
r B D O = k x G B L P H 2 2 1 K x B D O
r S P = k S P x G B L

2.2.3. Model 3

Model 3 is based on power law kinetics combined with an inhibition parameter added to the hydrogen partial pressure, which would describe slower kinetics due to a harder accessibility of the reactants to the active copper sites. Here, the adsorption mechanism of hydrogen is associative, meaning that each hydrogen molecule takes one site on the catalyst. The activity of hydrogen for side reactions was neglected because of its low effect on kinetics. The following reaction steps were considered for the construction of Model 3 (Table 3):
In this model, the production rates also depend on the surface coverage of hydrogen on the catalyst θ H 2 as well as the surface coverage of free sites θ s , as follows:
r G B L = k 1 K x B D O θ s 2 x G B L θ H 2 2 k S P x G B L
r B D O = k x G B L θ H 2 2 1 K x B D O θ s 2
The production rate rSP for the side products (SP) is the same as Equation (12) of Model 2.
A steady state is assumed for the adsorption of H2 onto the catalyst and leads to expressions of the surface coverage of hydrogen on the catalyst θ H 2 (Equation (15)) as well as the surface coverage of free sites   θ s (Equation (16)).
θ H 2 = K H 2 P H 2 1 + K H 2 P H 2
θ s = 1 θ H 2
Moreover, because of the large excess of hydrogen, the total pressure within the reactor is assumed to be the partial pressure of hydrogen. The production rates used for Model 3 are:
r G B L = k 1 K x B D O 1 K H 2 P H 2 1 + K H 2 P H 2 2 x G B L K H 2 P H 2 1 + K H 2 P H 2 2 k S P x G B L  
r B D O = k x G B L K H 2 P H 2 1 + K H 2 P H 2 2 1 K x B D O 1 K H 2 P H 2 1 + K H 2 P H 2 2

2.2.4. Model 4

Model 4 is based on Model 3, assuming a dissociative H2 sorption instead of an associative H2 sorption. For the dissociative H2 sorption, each hydrogen atom takes up one site on the catalyst. Model 4 was built similarly to Model 3. The activity of hydrogen for side reactions was neglected because of its low effect on kinetics.
The proposed reaction mechanism is given in Table 4. The induced production rates for the dissociative H2 sorption are written in Equations (19) and (20).
r G B L = k 1 K x B D O θ s 4 x G B L θ H 4 k S P x G B L  
r B D O = k x G B L θ H 4 1 K x B D O θ s 4
As for Model 3, the production rate rSP for the side products (SP) is the same as Equation (12) of Model 2 and a steady state is assumed for the adsorption of H2 onto the catalyst, and leads to the expressions of the surface coverage of hydrogen on the catalyst θ H (Equation (21)) as well as the surface coverage of free sites θ s (Equation (22)).
θ H = K H 2 P H 2 1 + K H 2 P H 2
θ s = 1 θ H
Moreover, because of the large excess of hydrogen, the total pressure within the reactor is assumed to be the partial pressure of hydrogen. The terms of production rates used for the simulation are:
r G B L = k 1 K x B D O 1 K H 2 P H 2 1 + K H 2 P H 2 4 x G B L K H 2 P H 2 1 + K H 2 P H 2 4 k S P x G B L  
r B D O = k x G B L K H 2 P H 2 1 + K H 2 P H 2 4 1 K x B D O 1 K H 2 P H 2 1 + K H 2 P H 2 4

2.3. Thermodynamic Equilibrium Modeling

The overall reaction equilibrium constant K is calculated through the standard Gibbs free energy change of reaction, ΔrG°:
K T = exp Δ r G ° T R T  
The standard Gibbs free energy change of reaction, ΔrG°, is calculated by Hess’s law, using formation enthalpies ΔfH°j and entropies ΔfH°j for GBL and BDO. Thermodynamics properties were extracted from the NIST database and completed by the ProPhyPlus database and can be seen in Table 5.
Δ r G ° T = Δ f G ° B D O T Δ f G ° G B L T  
Δ f G ° j T = Δ f H ° j T T Δ f S ° j T  
Formation enthalpies and entropies are calculated at a given temperature T according to Kirchhoff’s laws:
Δ f H ° j T = Δ f H ° j 298 K + Δ C p j · T T r e f  
Δ f S ° j T = Δ f S ° j 298 K + Δ C p j · L n T T r e f  
Finally, the heat capacities of the component Cpj are calculated according to Hess’s law, between the component and its hypothetical decomposition into elements:
Δ C p j = C p j v a p l = C ; H 2 ; O 2 υ l C p l

2.4. Numerical Methods

All numerical operations were done using the software MATLAB. Differential equation systems were solved using a modified Rosenbrock formula of order 2 (function ode23s in MATLAB).
Parameters were estimated using a trust–region–reflective least squares algorithm (function lsqnonlin in MATLAB). The minimized objective function J(a) was the sum of gaps between the experimental data and the modeled data squared, represented through Equation (31).
J a = i = 1 n x i e x p x i m o d 2
For Model 1 and 2, parameters were estimated for each set of reaction conditions independently. For Model 3 and 4, parameters were estimated for the experiments at 200 and 240 °C independently and for the two experiments at 220 °C simultaneously in order to obtain kinetic constants independent of pressure conditions. For the last parameter estimation for Model 3, wherein both pre-exponential and activation energy parameters were simultaneously estimated, all the experimental data were used.
Because the uncertainties on experimental data were not available, the precision of estimated parameters only takes into account the uncertainties due to the modeling. Uncertainties on estimated parameters were calculated as follows:
k ^ t N M p C o v k ^ < k < k ^ + t N M p C o v k ^
The covariance matrix C o v k ^ was calculated through the Jacobian matrix given by the function lsqnonlin of MATLAB and tN-Mp is the Student number.

3. Results and Discussion

3.1. Kinetic Model Comparison

For Model 1, the model proposed by J.H. Schlander [25], the estimated kinetic constants and their uncertainties for a 97.5% confidence interval are given in Table 6. The confidence intervals for kTHF and kBuOH show that these parameter values obviously are not significant.
This is due to the very large number of parameters to be estimated compared to the amount of experimental data available. For example, kTHF and kBuOH should not be evaluated separately because all the secondary products produced by the two parallel dehydrations are drawn in a single dataset from the experiments. Moreover, kGBL and kBDO can easily be made independent of one another by a proper modeling of the reaction equilibrium, which is usually reached in the experimental section. The redundancy of kinetic constants makes the model not sensitive towards the estimated parameters, since a default in one kinetic constant can be balanced by another. This gives very large uncertainties on the estimated parameters, especially for the side reactions, thus making this model impossible to extrapolate. Moreover, the model does not take into account the effect of pressure. These unreliable results motivated the development of Model 2 (Table 7).
Although we took into account the hydrogen pressure in a new model (Model 2), this model displayed a significant gap for the k kinetic parameter at 220 °C between 25 and 35 bar. Reaction rates are slower at higher pressure, which is not taken into account in that approach. Moreover, the confidence intervals for the k constant did not overlap for the experiments at 220 °C, meaning that the effect of pressure was not correctly modeled. So, it was impossible to fit an Arrhenius law for this model. Additionally, there was a non-significant k kinetic parameter at 240 °C. Indeed, uncertainties were low for every constant except for the k constant at 240 °C. This is explainable by the fact that the equilibrium was immediately reached right after the first point, and thus made the regression of this k constant happen between just the two first points. Therefore, simple power law models (Model 1 and 2) are inefficient when accounting for the drop of reaction rates as the pressure rises.
Model 3 (Table 8) was able to fit both experiments at 220 °C with a single set of constants in accordance to Arrhenius’s law where kinetic parameters were only dependent on temperatures. Parameter estimation showed that the adsorption constant K H 2   was almost constant in the studied range of temperatures. This result is not surprising, considering that activation energies for adsorption constants are usually low, thus making the dependency of adsorption constants towards temperature weak.
Regarding Model 4, the estimated kinetic constants and their uncertainties for a 97.5% confidence interval are given in Table 9. The estimated constants for kSP and K H 2   kinetic parameters were quite similar to those of Model 3. For the k parameter, some differences in terms of value appeared for all the temperatures, while as for Model 3, the k parameter remained non-significant at 240 °C.
Comparing Model 3 and 4 shows that, while kinetic constants were slightly different, the accuracy of the models was exactly the same, whether the adsorption was dissociative or associative. The model needed to slow down the reaction rates between the two pressure points of 25 and 35 bar regardless of the mathematical expression of this phenomenon. To evaluate which phenomenon, hydrogen dissociative or associative adsorption, was the limiting one taking place, further experiments with a wider range of pressures need to be performed.

3.2. Temperature Extrapolation

In order to design a LOHC process, a kinetic model able to predict the kinetics of the GBL hydrogenation into BDO at a given temperature and pressure, between 200 and 240 °C and between 25 and 35 bar, is needed. This paragraph is devoted to the extrapolation of Model 3, which was one of the only models successful in describing the slower kinetics at higher pressures. Extrapolation of Model 4 is seen to be unnecessary due to its similarity with Model 3. Extrapolation of pressures is built in the mathematical expression of the different models. For the extrapolation of temperatures, kinetic constants are modeled by Arrhenius’s law:
k i = A i exp E a i R T
Pre-exponential factors Ai and activation energies Eai were re-estimated by replacing the kinetic parameters by their respective Arrhenius’s laws and using all the experimental data at once. Because of the added parameters, the problem is highly dependent on initial conditions. To find a good set of initial parameters, exponential fitting of k and kSP constants into an Arrhenius’s law were made. However, the pre-exponential factor for the k parameter given by the exponential fitting proved to be several orders of magnitude off a good initialization point. In order to find a better set of initial parameters, we noted that the uncertainty calculation on the k constant at 240 °C displayed a lot of play in the constant k, because the equilibrium was obtained right after the first point (Figure 3). Therefore, the k constant of experiment (b) was manually adjusted in order to fit perfectly into an Arrhenius’s law. This adjustment was made through an exponential fit made with the estimated k constants for 200 and 220 °C, in furtherance of calculating the k constant of experiment (b) upon the exponential fit. Initial constants found with the retrofitted k constant showed decent results in the model fitting. Exponential fitting of kSP constants into an Arrhenius’s law proved to be good enough and no adjustment was made here.
Model testing showed that the adsorption constant K H 2 was almost constant in the studied range of temperatures. Therefore, the K H 2   constant was manually adjusted to 1.03 (which fitted all confidence intervals) and treated as a constant across the studied range of temperatures. Table 10 displays the initial parameters used.
Estimation of parameters using the aforementioned method led to a good fit of the extrapolated model to the experimental data (Figure 4).
The estimated pre-exponential factors and activation energies are given in Table 11 along with parity diagrams in Figure 5. For all the parameters, the confidence intervals were narrow and showed good reliability of the parameter values.

4. Conclusions

In this paper, a kinetic model has been developed for the gas-phase hydrogenation of γ-butyrolactone in the context of a hydrogen storage process. The data were gathered from experiments conducted by J.H. Schlander on the valorization of succinic acid. Indeed, several works addressing the hydrogenation of both maleic anhydride and succinic acid towards 1,4-butanediol proved to be useful to study the γ-butyrolactone–1,4-butanediol LOHC couple. The newly identified LOHC pair is definitely showing high potential of becoming a competitive LOHC system. With its much improved energy efficiency over cycloalkane-based LOHCs, and better availability than N-heterocycle LOHCs while requiring cheaper catalysts, this pair of molecules could meet the demand for a LOHC of renewable origin, readily available and with low operating costs. Currently, the γ-butyrolactone–1,4-butanediol LOHC lacks reliable kinetic models and data in the literature, thus preventing the development of efficient three-phase reactors and processes. Therefore, this paper proposes for the first time a temperature and pressure extrapolated kinetic model for the gas phase hydrogenation of this LOHC. The model will be included in a more general model for the design of a three-phase hydrogenation reactor currently under study.
Carrying out the hydrogenation of pure GBL in the liquid phase poses difficult problems such as potential dry spots on the catalyst surface. While good selectivity and conversion can be achieved with reasonable reaction conditions in the liquid phase, the presence of dry zones could lead to a drop in selectivity through the production of by-products in the gas phase, making it important to properly assess the impact of gas and liquid phase kinetics for LOHC processes.

Supplementary Materials

The following supporting information can be downloaded at: https://mdpi.longhoe.net/article/10.3390/reactions3040033/s1. The document consists of 4 sections. Section S1: Comparison of different LOHC systems. This section contains a table with the properties of the LOHCs. Section S2: Experimental data extracted from the PhD thesis of J.H. Schlander. This section gives additional information on the data used for parameter estimation, as well as a table with the catalyst properties, a scheme of the active part of the reactor and a figure of the experimental data. Section S3: Intraparticular diffusion and capillary condensation. This section gives results and calculation methods on transfer and capillary condensation. Section S4: Ideal gas hypothesis. This section gives results and calculation methods on the thermodynamic properties of the gas phase, as well as a figure with the calculated fugacity coefficients.

Author Contributions

Conceptualization, A.C. and I.P.; methodology, A.C. and I.P.; software, V.G. and I.C.; validation, all; writing—original draft preparation, V.G.; writing—review and editing, A.C. and I.P.; funding acquisition, A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CEA.

Data Availability Statement

All the produced data are in the text.

Acknowledgments

The authors would like to thank the Carnot Institute for the financial support, Françoise COUENNE and Christian JALLUT (LAGEPP) for their advice on parameters estimation, as well as Sebastien THOMAS, Anne-Cécile ROGER (ICPEES) and Chaimae QASMI (CEA) for the fruitful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Abbreviation
BDO1,4-butanediol
GBLγ-butyrolactone
H0-BTbenzyltoluene
H0-DBTdibenzyltoluene
H0-NECN-ethylcarbazole
H12-BTperhydro-benzyltoluene
H12-NECperhydro-N-ethylcarbazole
H18-DBTperhydro-dibenzyltoluene
LOHCliquid organic hydrogen carrier
MCHmethylcyclohexane
TOLtoluene
Latine letters
A i Pre-exponential factor (s−1)
a Set of estimated parameters
C o v Covariance matrix
C p j Specific heat capacity of elements (j: H2, O2, Cgraphite) (J·K−1·mol−1)
C p j v a p Specific heat capacity of component j in vapour phase (j: GBL, BDO) (J·K−1·mol−1)
Δ C p j Specific heat capacity of component j (j: GBL, BDO) (J·K−1·mol−1)
d p Particule size of the catalyst (m)
d t Reactor diameter (m)
E a Activation energy (J·mol−1)
Δ f H ° j Standard enthalpy of formation for component j (j: GBL, BDO) (J·mol−1)
F j Molar flux for compound j (j: GBL, BDO, THF, BuOH) (mol·s−1)
Δ f G ° j Standard Gibbs free energy of formation for component j (j: GBL, BDO) (J·mol−1)
Δ r G ° Standard Gibbs free energy change of reaction (J·mol−1)
J a Objective function for the least-square solver (-)
K Equilibrium constant for the hydrogenation of GBL into BDO (Bar−2)
K H 2 Adsorption constant of H2 on the catalyst (Bar−1)
k Kinetic parameter for the hydrogenation of GBL into BDO (s−1)
k S P Kinetic parameter for side reactions (s−1)
k ^ Estimated parameter
L B Reactor length (m)
m c a t a Catalyst loading (kg)
M p Number of parameters
N Number of experimental data points
P Pressure (Bar)
P H 2 Partial pressure of hydrogen (Bar)
Q Gas volumetric flow rate (m3·s−1)
R Molar gas constant (J·K−1·mol−1)
r j Production rate for compound j (j: GBL, BDO, THF, BuOH, side products) (s−1)
r P l ¯ Observed reaction rate for reaction i (mol·s−1·m−3catalyst)
r s i Intrinsic reaction rate for reaction i (mol·s−1·m−3catalyst)
T Temperature (K)
t m V Modified residence time (kgcatalyst·s·mR−3)
t N M p Student factor for a N M p degree of freedom.
t V Residence time (s)
V c a t a Volume occupied by the catalyst (m3)
V R Volume of the reactive section (m3)
x j Molar fraction for condensable component j (j: GBL, BDO, Side products) (-)
x k e x p Molar fraction extracted from experimental data for a residence time k (-)
x k m o d Molar fraction calculated by a model for a residence time k (-)
Greek letters
η i Efficiency factor for reaction i (-)
θ H Surface coverage of hydrogen on the catalyst (-)
θ H 2 Surface coverage of dihydrogen on the catalyst (-)
θ s Surface coverage of free sites on the catalyst (-)
υ i j Stoichiometric coefficient for reaction i and for compound j (j: GBL, BDO, THF, BuOH, side products) (-)

References

  1. Aakko-Saksa, P.T.; Cook, C.; Kiviaho, J.; Repo, T. Liquid organic hydrogen carriers for transportation and storing of renewable energy—Review and discussion. J. Power Sources 2018, 396, 803–823. [Google Scholar] [CrossRef]
  2. Makridis, S.S. Hydrogen storage and compression. In Methane and Hydrogen for Energy Storage, 1st ed.; Rupp, C., David, S.-K.T., Eds.; IET Digital Library: Stevenage, UK, 2016; pp. 1–28. [Google Scholar] [CrossRef]
  3. Müller, K. Technologies for the Storage of Hydrogen Part 1: Hydrogen Storage in the Narrower Sense. ChemBioEng Rev. 2019, 6, 72–80. [Google Scholar] [CrossRef]
  4. Dalebrook, A.F.; Gan, W.; Grasemann, M.; Moret, S.; Laurenczy, G. Hydrogen storage: Beyond conventional methods. Chem. Commun. 2013, 49, 8735–8751. [Google Scholar] [CrossRef] [PubMed]
  5. Huang, Z.; Barnett, K.J.; Chada, J.P.; Brentzel, Z.J.; Xu, Z.; Dumesic, J.A.; Huber, G.W. Hydrogenation of γ-Butyrolactone to 1,4-Butanediol over CuCo/TiO2 Bimetallic Catalysts. ACS Catal. 2017, 7, 8429–8440. [Google Scholar] [CrossRef]
  6. Huang, Y.; Cheng, Y.; Zhang, J. A Review of High Density Solid Hydrogen Storage Materials by Pyrolysis for Promising Mobile Applications. Ind. Eng. Chem. Res. 2021, 60, 2737–2771. [Google Scholar] [CrossRef]
  7. Teichmann, D.; Arlt, W.; Wasserscheid, P.; Freymann, R. A future energy supply based on Liquid Organic Hydrogen Carriers (LOHC). Energy Environ. Sci. 2011, 4, 2767–2773. [Google Scholar] [CrossRef]
  8. He, T.; Pei, Q.; Chen, P. Liquid organic hydrogen carriers. J. Energy Chem. 2015, 24, 587–594. [Google Scholar] [CrossRef]
  9. Meille, V.; Pitault, I. Liquid Organic Hydrogen Carriers or Organic Liquid Hydrides: 40 Years of History. Reactions 2021, 2, 94–101. [Google Scholar] [CrossRef]
  10. Makaryan, I.A.; Sedov, I.V.; Maksimov, A.L. Hydrogen Storage Using Liquid Organic Carriers. Russ. J. Appl. Chem. 2020, 93, 1815–1830. [Google Scholar] [CrossRef]
  11. Preuster, P.; Papp, C.; Wasserscheid, P. Liquid Organic Hydrogen Carriers (LOHCs): Toward a Hydrogen-free Hydrogen Economy. Accounts Chem. Res. 2017, 50, 74–85. [Google Scholar] [CrossRef]
  12. Shi, L.; Qi, S.; Qu, J.; Che, T.; Yi, C.; Yang, B. Integration of hydrogenation and dehydrogenation based on dibenzyltoluene as liquid organic hydrogen energy carrier. Int. J. Hydrogen Energy 2019, 44, 5345–5354. [Google Scholar] [CrossRef]
  13. Chappaz, A.; Bengaouer, A. Use of Hydrogenated Organic Liquids, in Particular in Devices for Converting Energy. European Patent EP3725738, April 2020. [Google Scholar]
  14. Onoda, M.; Nagano, Y.; Fujita, K.-I. Iridium-catalyzed dehydrogenative lactonization of 1,4-butanediol and reversal hydrogenation: New hydrogen storage system using cheap organic resources. Int. J. Hydrogen Energy 2019, 44, 28514–28520. [Google Scholar] [CrossRef]
  15. Küksal, A.; Klemm, E.; Emig, G. Reaction kinetics of the liquid-phase hydrogenation of succinic anhydride on CuZnO-catalysts with varying copper-to-zinc ratios in a three-phase slurry reactor. Appl. Catal. A: Gen. 2002, 228, 237–251. [Google Scholar] [CrossRef]
  16. Gräfje, H.; Körnig, W.; Weitz, H.-M.; Reiß, W.; Steffan, G.; Diehl, H.; Bosche, H.; Schneider, K.; Kieczka, H.; Pinkos, R. Butanediols, Butenediol, and Butynediol. In Ullmann’s Encyclopedia of Industrial Chemistr; Wiley: Hoboken, NJ, USA, 2022; pp. 1–12. [Google Scholar] [CrossRef]
  17. Global and China 1,4-Butanediol (BDO) Industry Report, 2017–2021. 2017. Available online: https://www.reportlinker.com/p0318800/Global-and-China-1-4-Butanediol-BDO-Industry-Report.html (accessed on 23 April 2021).
  18. Satta, R.; Dimitrijevic, N.; Manev, H. Drosophila metabolize 1,4-butanediol into γ-hydroxybutyric acid in vivo. Eur. J. Pharmacol. 2003, 473, 149–152. [Google Scholar] [CrossRef]
  19. Schwarz, W.; Schossig, J.; Rossbacher, R.; Höke, H. Butyrolactone. In Ullmann’s Encyclopedia of Industrial Chemistr; Wiley: Hoboken, NJ, USA, 2000; pp. 1–7. [Google Scholar] [CrossRef]
  20. Aslam, R.; Khan, M.H.; Ishaq, M.; Müller, K. Thermophysical Studies of Dibenzyltoluene and Its Partially and Fully Hydrogenated Derivatives. J. Chem. Eng. Data 2018, 63, 4580–4587. [Google Scholar] [CrossRef]
  21. Davies, M.K.; Archer, A.L.; Maczek, A.O.S.; Manzanares-Papayanopoulos, A.E.; A McLure, I. Shear Viscosities of Methycyclohexane, Perfluoromethylcyclohexane, and Their Mixtures in the Vicinity of the Upper Critical Mixing Temperature. 1. Critical Isopleth and Coexistence Curve. J. Chem. Eng. Data 2006, 51, 1502–1508. [Google Scholar] [CrossRef]
  22. Ozokwelu, E.D.; Staff, U.B. (Eds.) Toluene. In Kirk-Othmer Encyclopedia of Chemical Technology; Wiley: Hoboken, NJ, USA, 2006; pp. 1–24. [Google Scholar] [CrossRef]
  23. Chaudhari, R.; Jaganathan, R.; Vaidya, S.; Chaudhari, S.; Naik, R.; Rode, C. Hydrogenation of diethyl maleate in a fixed-bed catalytic reactor: Kinetics, reactor modelling and pilot plant studies. Chem. Eng. Sci. 1999, 54, 3643–3651. [Google Scholar] [CrossRef]
  24. Herrmann, U.; Emig, G. Liquid Phase Hydrogenation of Maleic Anhydride and Intermediates on Copper-Based and Noble Metal Catalysts. Ind. Eng. Chem. Res. 1997, 36, 2885–2896. [Google Scholar] [CrossRef]
  25. Schlander, J.H. Gasphasenhydrierung von Maleinsäuredimethylester zu 1.4-Butandiol, γ-Butyrolacton und Tetrahydrofuran an Kupfer-Katalysatoren. Institut für Chemische Verfahrenstechnik. Ph.D. Thesis, Universität Karlsruhe, Karlsruhe, Germany, 2000. [Google Scholar]
  26. And, J.H.S.; Turek, T. Gas-Phase Hydrogenolysis of Dimethyl Maleate to 1,4-Butanediol and γ-Butyrolactone over Copper/Zinc Oxide Catalysts. Ind. Eng. Chem. Res. 1999, 38, 1264–1270. [Google Scholar] [CrossRef]
  27. Ichikawa, N.; Sato, S.; Takahashi, R.; Sodesawa, T.; Inui, K. Dehydrogenative cyclization of 1,4-butanediol over copper-based catalyst. J. Mol. Catal. A. Chem. 2004, 212, 197–203. [Google Scholar] [CrossRef]
  28. Herman, R.; Klier, K.; Simmons, G.; Finn, B.; Bulko, J.; Kobylinski, T. Catalytic synthesis of methanol from CO/H2I. Phase composition, electronic properties, and activities of the Cu/ZnO/M2O3 catalysts. J. Catal. 1979, 56, 407–429. [Google Scholar] [CrossRef]
  29. National Institute of Standards and Technlogy, U.S. Department of Commerce. Butyrolactone. n.d. Available online: https://webbook.nist.gov/cgi/cbook.cgi?ID=C96480&Units=SI&Mask=1#Refs (accessed on 9 June 2021).
  30. National Institute of Standards and Technlogy, U.S. Department of Commerce. 1,4-Butanediol. n.d. Available online: https://webbook.nist.gov/cgi/cbook.cgi?ID=C110634&Units=SI&Mask=1#Thermo-Gas (accessed on 9 June 2021).
Figure 1. Reaction scheme used for the construction of the kinetic model.
Figure 1. Reaction scheme used for the construction of the kinetic model.
Reactions 03 00033 g001
Figure 2. GBL (times cross), BDO (circle) and side products (SP) (star) molar fraction vs. time on stream (tV): (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C. The data come from the Ph.D. thesis of J.H. Schlander [25], but were rearranged as a function of tV instead of the modified dwell time tmV.
Figure 2. GBL (times cross), BDO (circle) and side products (SP) (star) molar fraction vs. time on stream (tV): (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C. The data come from the Ph.D. thesis of J.H. Schlander [25], but were rearranged as a function of tV instead of the modified dwell time tmV.
Reactions 03 00033 g002
Figure 3. Exponential fitting of k constants of Model 3 between 200 and 240 °C before adjustment.
Figure 3. Exponential fitting of k constants of Model 3 between 200 and 240 °C before adjustment.
Reactions 03 00033 g003
Figure 4. Simulation of the extrapolation of Model 3 with the adjusted parameters. GBL (experimental: times cross; model: blue line), BDO (experimental: circle; model: red line) and side products (experimental: star; model: green line) molar fraction vs. time on stream (tV): (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C.
Figure 4. Simulation of the extrapolation of Model 3 with the adjusted parameters. GBL (experimental: times cross; model: blue line), BDO (experimental: circle; model: red line) and side products (experimental: star; model: green line) molar fraction vs. time on stream (tV): (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C.
Reactions 03 00033 g004
Figure 5. Parity diagram of the extrapolation of Model 3 with the adjusted parameters. GBL (blue times cross), BDO (red circle) and side products (green star). Experimental molar fraction vs. modeled molar fraction: (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C. Dotted lines represent ± 15% deviations.
Figure 5. Parity diagram of the extrapolation of Model 3 with the adjusted parameters. GBL (blue times cross), BDO (red circle) and side products (green star). Experimental molar fraction vs. modeled molar fraction: (a) 35 bar, 200 °C; (b) 25 bar, 240 °C; (c) 35 bar, 220 °C; (d) 25 bar, 220 °C. Dotted lines represent ± 15% deviations.
Reactions 03 00033 g005
Table 1. Reaction mechanism for the construction of the first power law model: Model 1.
Table 1. Reaction mechanism for the construction of the first power law model: Model 1.
ReactionsKinetic ParametersHypothesis
G B L + 2   H 2 B D O k G B L ,   k B D O Main reaction
G B L + 2 H 2 T H F + H 2 O k T H F Side reaction
G B L + 3 H 2 B u O H + H 2 O k B u O H Side reaction
Table 2. Reaction mechanism used for the construction of the second power law model: Model 2.
Table 2. Reaction mechanism used for the construction of the second power law model: Model 2.
ReactionsKinetic or Thermodynamic ParametersHypothesis
G B L + 2   H 2 B D O k ,   K Rate-determining step
G B L + n   H 2 S P k S P Side reaction
Table 3. Reaction mechanism used for the construction of the kinetic model with inhibition of hydrogen partial pressure: Model 3.
Table 3. Reaction mechanism used for the construction of the kinetic model with inhibition of hydrogen partial pressure: Model 3.
ReactionsKinetic or Thermodynamic ParametersHypothesis
H 2 + s H 2 s K H 2 Steady-state
G B L + 2   H 2 s B D O + 2   s k ,   K Rate-determining step
G B L + n   H 2 S P k S P Side reaction
Table 4. Reaction mechanism used to build the kinetic model based on a dissociative adsorption of H2: Model 4.
Table 4. Reaction mechanism used to build the kinetic model based on a dissociative adsorption of H2: Model 4.
ReactionsKinetic or Thermodynamic ParametersHypothesis
H 2 + 2   s 2   H s K H 2 Steady-state
G B L + 4   H s B D O + 4   s k ,   K Rate-determining step
G B L + n   H 2 S P k S P Side reaction
Table 5. Thermodynamic data used for the modeling of the equilibrium constant. Data were extracted from the databases of NIST and ProPhyPlus.
Table 5. Thermodynamic data used for the modeling of the equilibrium constant. Data were extracted from the databases of NIST and ProPhyPlus.
Component Δ f H ° 298 K   k J · m o l 1   Δ f S ° 298 K   J · K 1 · m o l 1   Δ C p   J · K 1 · m o l 1   C p v a p   J · K 1 · m o l 1   C p s o l J · K 1 · m o l 1  
GBL−365−314 *−62.328.1-
BDO−427−500 *−82.323123.697 *-
H2---28.8-
O2---29.1-
C graphite----8.23
* These data come from the ProPhyPlus database [29,30].
Table 6. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 1.
Table 6. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 1.
Experiment k G B L   ( s 1 ) k B D O   ( s 1 ) k T H F   ( s 1 ) k B u O H   ( s 1 )
240 °C|25 bar(1.50 ± 1.28) × 10−1(7.47 ± 6.50) × 10−15.09 × 10−3 ± 1.01 × 10−25.04 × 10−4 ± 1.01 × 10−2
220 °C|25 bar(3.73 ± 0.321) × 10−2(1.00 ± 0.098) × 10−19.91 × 10−4 ± 7.62 × 10−39.74 × 10−4 ± 7.56 × 10−3
220 °C|35 bar(2.31 ± 0.157) × 10−2(3.05 ± 0.265) × 10−2(1.05 ± 9.68) × 10−39.07 × 10−4 ± 9.62 × 10−3
200 °C|35 bar(2.49 ± 2.72) × 10−2(2.94 ± 3.66) × 10−25.13 × 10−4 ± 1.505.47 × 10−4 ± 1.5
Table 7. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 2.
Table 7. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 2.
Experiment k   s 1 k S P   s 1
240 °C|25 bar(0.325 ± 2.10) × 10−3(5.65 ± 0.244) × 10−3
220° C|25 bar(6.16 ± 0.581) × 10−5(2.07 ± 0.120) × 10−3
220 °C|35 bar(1.88 ± 0.152) × 10−5(2.10 ± 0.136) × 10−3
200 °C|35 bar(1.06 ± 0.0843) × 10−5(1.08 ± 0.226) × 10−3
Table 8. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 3.
Table 8. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 3.
Temperature (°C) k   s 1 k S P   s 1 K H 2   B a r 1
240(2.31 ± 5.92) × 10−1(5.61 ± 0.215) × 10−31.035 ± 0.0381
220(2.88 ± 0.185) × 10−2(2.04 ± 0.193) × 10−31.018 ± 0.0337
200(1.33 ± 0.135) × 10−2(1.04 ± 0.229) × 10−31.036 ± 0.0840
Table 9. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 4.
Table 9. Estimated kinetic parameters and their uncertainty for a 97.5% confidence interval for Model 4.
Temperature (°C) k   s 1 k S P   s 1 K H 2   B a r 1
2400.387 ± 1.30(5.61 ± 0.220) × 10−31.032 ± 0.0372
220(5.34 ± 1.37) × 10−2(2.00 ± 0.223) × 10−31.032 ± 0.0388
200(2.32 ± 0.263) × 10−2(1.04 ± 0.229) × 10−31.037 ± 0.0837
Table 10. Initial parameters used for the parameters estimation of the extrapolated Model 3.
Table 10. Initial parameters used for the parameters estimation of the extrapolated Model 3.
A H y d r o   s 1 E a H y d r o   J · m o l 1 A S P   s 1 E a S P   J · m o l 1 K H 2   B a r 1
2.50 × 10674 × 9412.23 × 10684 × 7741.03
Table 11. Estimated parameters for Model 3 with 95% confidence intervals.
Table 11. Estimated parameters for Model 3 with 95% confidence intervals.
A H y d r o   s 1 E a H y d r o J · m o l 1   A S P   s 1 E a S P   J · m o l 1 K H 2   B a r 1
(1.97 ± 0.38) × 106(73.9 ± 0.9) × 103(1.67 ± 0.46) × 106(83.5 ± 3.3) × 1031.05 ± 0.06
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gautier, V.; Champon, I.; Chappaz, A.; Pitault, I. Kinetic Modeling for the Gas-Phase Hydrogenation of the LOHC γ-Butyrolactone–1,4-Butanediol on a Copper-Zinc Catalyst. Reactions 2022, 3, 499-515. https://doi.org/10.3390/reactions3040033

AMA Style

Gautier V, Champon I, Chappaz A, Pitault I. Kinetic Modeling for the Gas-Phase Hydrogenation of the LOHC γ-Butyrolactone–1,4-Butanediol on a Copper-Zinc Catalyst. Reactions. 2022; 3(4):499-515. https://doi.org/10.3390/reactions3040033

Chicago/Turabian Style

Gautier, Vincent, Isabelle Champon, Alban Chappaz, and Isabelle Pitault. 2022. "Kinetic Modeling for the Gas-Phase Hydrogenation of the LOHC γ-Butyrolactone–1,4-Butanediol on a Copper-Zinc Catalyst" Reactions 3, no. 4: 499-515. https://doi.org/10.3390/reactions3040033

Article Metrics

Back to TopTop