The initial step of works on the designing of a mathematical model of a commercial process to manufacture a mixture of butyric acid hydroxyesters, for use in ecological solvent compositions for the paint and varnish industry, is presented in the paper. The works referred to in this paper relate the mathematical model of the condensation reactor. A unit process kinetic model was designed based on an analysis of the general differential form of the equation of the reactor mass balance and an analysis of the mechanism of the Tishchenko aldol condensation reaction, as used in the synthesis of butyric acid hydroxyesters, confirmed by the results of GC/MS analysis. The forms of the kinetic equations obtained were verified experimentally by determining the values of the respective reaction rate constants.
Mathematical modelling is an extremely valuable tool that enables the dramatic effect of cost reduction to be achieved in the area of experimental works on the designing and development of chemical technology, taking change of scale into account. The ability to simulate the course of commercial processes in silico enables the real-time operation of a chemical plant to be simulated in the aspects of its entire mass balance, energy balance, and financial balance. The present paper relates to the designing of the kinetic model of the unit process of condensation.
The Tishchenko reaction, which is based on a series of methods to obtain esters from the respective aldehydes were developed, raises much interest due to its numerous useful advantages. After discovering its unique features, enabling the effect of stereoselectivity of its course to be obtained, the Tishchenko reaction was rediscovered and appreciated in the 1990’s, specifically in connection with its features and applications in the synthesis of pure enantiomers, that had been demonstrated not long before.
A number of modified Tishchenko [
In connection with the high reactivity of the aldehydes as their specific feature, it is necessary to select the suitable reaction conditions in planning the course of the synthesis in order to reduce the impact of side processes [
The reaction which is the subject of this characterization is also known as the Claisen-Tishchenko reaction in connection with the series of experiments with benzaldehyde and sodium alcoholates, carried out by Claisen in 1887. As mentioned earlier, the Tishchenko aldol condensation runs in the alkaline catalyst environment and concerns aldehydes having a hydrogen atom α, determining the course of the aldol condensation. The reaction leads to the formation of monoesters of 1,3-diols.
Owing to its commercial importance, the above mentioned modification of the Tishchenko reaction has been investigated rather extensively [
The system of two isomeric esters of isobutyric acid is formed in the Tishchenko aldol condensation reaction, during which trimerisation of three molecules of isobutyric aldehyde occurs in accordance with the following general stoichiometric equation:
The process of trimerisation starts with a reaction of two molecules of 2-methylbu- tanal, running according to the aldol condensation mechanism, leading to the formation of the aldol: 3-hydroxy-2,2,4-trimethylpentanal [
Based on the analysis of the possible reaction paths during the process of isobutanal condensation in the presence of an alkaline catalyst as well as experimental findings signalled in patent literature [
1. Aldol condensation of 2 molecules of isobutanal.
1)
2. Condensation and disproportioning of isobutanal and aldol obtained in Reaction 1).
2)
3. Transesterification with diester formation.
3)
4. Partial hydrolysis of diester.
4)
5. Hydrolysis of monoester.
5)
6. The Canizzaro reaction of the aldol obtained in step 3).
6)
Acceptance of the proposed reaction mechanism for further studies is empirically justified by the findings of GC/MS analyses of samples of the reaction mixtures, collected from various steps of the process [
Isobutanal (2-methylpropanal)-pure, from GA ZA Kędzierzyn.
3-Hydroxy-2,2,4-trimethylpentyl isobutyrate-pure, from Sigma-Aldrich.
Sodium hydroxide, 40% solution-from POCH Gliwice S.A.
Acetic acid 80% pure, from POCH Gliwice S.A.
A thermostated (60˚C) reaction vessel with a weighed amount of post-reaction mixture was filled with a weighed amount of isobutyric aldehyde, and 2.7 g of 40% NaOH solution was added to the obtained mixture. The resulting reaction mixture was stirred vigorously until the time t, was then placed in an ice bath whereafter 2.5 g of 80% acetic acid solution was added. The mixture was cooled down to a temperature in the range 15˚C - 20˚C and analyzed by GC to determine the concentrations of the respective reacting substances.
The experiment was carried out for the times t: 2, 4, 8, 15, 30, 60 min, for reaction systems obtained by mixing 3-hydroxy-2,2,4-trimethylpentyl isobutyrate with isobutanal in quantitative ratios shown in
Condensation of isobutyric aldehyde is performed in continuous-flow tank reactor equipped with efficient propeller or turbine type stirrer, providing fulfillment of kinetic criterion:
where tmix denotes mixing time and t1/2 half-life time of reaction. In such conditions the isobutyric aldehyde condensation reactor may be approximated by an ideal continuous- flow tank reactor model (
where ρi is density of i―component, w―velocity field, ji―dispersion flow field of i, Mi―molecular weight of i, νri―stoichimetric factor of i in reaction r, r―reaction rate, R―total number of reaction involved in process.
Integration of (3) over volume V of reactor results Equation (4):
Considering Gauss-Ostrogradsky theorem, volume integral of divergence appearing in (4) may be converted to the surface integral over the closed oriented boundary F:
where
Since the mass flows coming from the surface surrounding inside of reactor flow only through the surfaces f0 and f1, the surface integral (11) may be expressed by an algebraic sum of simple integrals as follows:
Assuming the field of i-compound velocity on surfaces f0 and f1 equals the barycentric velocity and then in consequence negligibility of dispersion flow of i, what could be expressed by criteria:
then
In view of (10) and (14), assuming the process runs in condition of non-compressible fluid, the mass balance of continuous-flow tank reactor may be expressed by model equation, applicable for engineering purposes:
where
Considering the presence of the kinetic member in model Equation (9), in order to determine its form, an analysis of the isobutanal condensation reaction kinetics was performed as it follows below.
Acceptance of the proposed reaction mechanism for further studies, step 1)-6), is empirically justified by the findings of GC/MS analyses of samples of the reaction mixtures, collected from various steps of the process.
collected after 1 min of reaction.
Significant feature of the shape of this chromatographic curve occurs to be an appearance of peaks characteristic for the presence of mentioned above system hemiacetale/dioxanole and transitional aldole reaction product of 2 moles of isobutanal: 2-hy- droxy-2,2,4-trimethylpentanal. As it has been already stated, both undergo consequent Tishchenko reaction, what finds reflection in the shape of the chromatogram of the sample of post reaction mixture (
The shape of post reaction mixture sample proves the decay of transitional by products of Tishchenko reaction appearing at initial stages of reaction, hemiacetale/dioxa- nole system and 2-hydroxy-2,2,4-trimethylpentanal, which underwent esterification through the reaction run, according to the Tishchenko reaction mechanism scheme signaled earlier.
where:
C―reagents concentration column matrix, mol・dm−3.
νT―transpose of chemical matrix of the reaction system.
r―condensation process component reactions rates equations column matrix, mol・dm3・min−1.
Referring to the results of preliminarily carried out investigations of the reaction mixture compositions, showing negligibly low values of concentration of some components of these mixtures, it may be concluded the composition of the reaction mixture in conditions of the industrial reactor run will be determined only by the process described by step 1), 2) and 5). Therefore, the general kinetic Equation (10) may be formulated as follows:
where:
a―isobutanal concentration, mol・dm−3.
b―3-hydroxy-2,2,4-trimethylpentyl isobutyrate (the sum of both isomers) concentration, mol・dm−3.
c―concentration of 2,2,4-trimethyl-1,3-penthanodiole, mol・dm−3.
r1―sum of reaction (4) and (5) rate, mol・dm−3・min−1.
r2―reaction (8) rate, mol・dm−3・min−1.
where:
k1―isobutanal trimerisation reaction rate constant.
k2―hydrolysis reaction, step 4) rate constant.
Solution of the first member of ordinary differential equation system (11) using the separation of variables method results the particular integral of irrational form:
Solution of ordinary differential equation (ODE):
was obtained by Laplace transformation:
where:
and s is complexed variable defined as: s = x + iy.
Result of solution of step 1) obtained after inversion of transformation 2) may be expressed by following equation:
where:
The solution procedure of ODE (15) referred to function b(t) met some inconveniences connected with the properties of Laplace transformation and limited possibility of extraction of complex forms of transform inversions. In order to convert imaginary function
And
what allowed to introduce substitution:
where erfi[f(t)] denotes imaginary error function, expandable to convergent Maclaurin’s series [
Introduction of the series (27) made possible to output results of calculations using Equation (18) in the real range. Because of the difficulties with performance of the Laplace transformation (17) on expression (18) the analytical solution of (11) vs. concentration of diole was not been obtained.
The Equation (11) was also solved numerically using Runge-Kutta method [
The course of the curve of changes of the concentration of the hydroxyester indicates the presence, in the environment of the point of time t = 45 min, of a stationary point which is characterized by a maximum concentration of the hydroxyester. After the time is exceeded, the course of the model curve shows a decrease in the concentrations of the hydroxyester, which is both due to the reaction getting slower at lower concentrations of isobutanal and due to the course of the side reaction of hydrolysis of which the rate does not depend on the concentration of isobutanal. The last reaction is accompanied by an increase in the concentration of the diol. In the aspect of process kinetics and the purity of final product, leading to a reduction in capital expenditure and process costs connected with the utilization of the diol, the optimal time range of 10 - 20 min of the residence of the reaction mixture in the reaction zone is of commercial value.
In order to verify the adopted kinetic model, as described by the system of differential Equation (11), the measurement points obtained during kinetic tests were approximated by means of curves obtained by solving Equation (11) numerically, by minimizing the sums of squared deviations between the values predicted by the model and those obtained experimentally.
Figures 5-9 show the graphs of dependence of the concentrations of isobutanal, hydroxyester, and diol vs. time, on which the values of the measurement points obtained during the experiments were plotted along with the curves which are approximations of experimental results, obtained by means of general integrals of Equation (11).
The analysis of the mechanism of isobutanal condensation that was carried out indicates the possibility of assumption that the process may run in a single step of trimerization of three molecules of the aldehyde to 2,2,4-trimethyl-1,3-pentanediol 2-methyl- propionate, and that the side reaction of product hydrolysis produces 2,2,4-trimethyl- 1,3-pentanediol. Based on experimental findings, the kinetic equation of the aldol homocondensation of isobutanal was formulated as the dependence of reaction rate vs. time. The set of kinetic equations was solved by independent mathematical and numerical methods. Measurement points, obtained during kinetic studies, were approximated with experimentally obtained values in order to verify the adopted kinetic model. A high adequacy of the adopted kinetic model of the isobutanal condensation process was obtained as the result. Reaction rate constants were found for competitive reactions, determining the process yields, i.e., those of trimerization and hydrolysis, which were k1 = 0.00122 (dm3)2・mol−2・min−1; k2 = 0.0035 min−1, respectively.
This research work was financed as a R&D project No. 13856, with the support of National Research and Development Centre, Poland.
Tic, W.J. (2016) Kinetic Study of the Aldol Condensation of Isobutyraldehyde Catalyzed by Sodium Hy- droxide. Advances in Chemical Engineering and Science, 6, 456-469. http://dx.doi.org/10.4236/aces.2016.64041