Concentration Wave for a Class of Reaction Chromatography System with Pulse Injections

By using fluid dynamics theory with the effects of adsorption and reaction, the chromatography model with a reaction A → B was established as a system of two hyperbolic partial differential equations (PDE’s). In some practical situations, the reaction chromatography model was simplified a semi-coupled system of two linear hyperbolic PDE’s. In which, the reactant concentration wave model was the initial-boundary value problem of a self-closed hyperbolic PDE, while the resultant concentration wave model was the initial-boundary value problem of hyperbolic PDE coupling reactant concentration. The general explicit expressions for the concentration wave of the reactants and resultants were derived by Laplace transform. The δ-pulse and wide pulse injections were taken as the examples to discuss detailedly, and then the stability analysis between the resultant solutions of the two modes of pulse injection was further discussed. It was significant for further analysis of chromatography, optimizing chromatographic separation, determining the physical and chemical characters.


Introduction
With the appearance of diverse production chromatography (such as the reaction chromatography), the chromato-graphy technology has been widely applied in chemistry, chemical engineering, biological engineering and pharmaceutical engineering, etc., while the demand of chromatography theory is increasing higher.The relationships among the chromatographic input-output and the system conditions play the very important role in chromatography model [1]- [6].
In fact, the mathematical model of chromatography system is a initial-boundary value problem of hyperbolic partial differential equations system [7]- [11], which is hard and challenging mathematics problem to chromatography scientists.In the other hand, the practical application and demand for chromatography is also difficult to understand deeply by mathematicians.The relative works of partial differential equations in the practical chromatography are still not enough.
If the chromatographic process contains reactions, it is labeled as reaction chromatography.An important example is the catalyst for the column packing, accompanied the catalytic [2]- [6] in the adsorption process, and the isomerization reaction is the common situation.
In this paper, a chromatography model with a reaction A B → was established, which is a initial-boundary value problem for the semi-coupled system of two linear hyperbolic partial differential equations.Then the general explicit expressions of concentration waves for reactant and resultant were derived using Laplace transform.It was significant for further analysis between input and output of chromatography, optimizing chromatographic separation, determining the physical and chemical characters.Finally, the δ-pulse and wide pulse injections were taken as the examples to discuss detailedly, and then the stability analysis between the resultant solutions of the two modes of pulse injection was further discussed.The results provided proper theory models for further chromatographic data analysis.

Reaction Chromatography Model
Set the concentrations of the reactant A and the resultant B in the mobile phase and in the stationary phase as , , , c c f f respectively.Reaction rate was r k .And the linear velocity of the mobile phase was u.The volume shares in chromatographic column in the mobile phase and in the stationary phase as , ε µ , respectively.
Denoted that F µ ε = , then the mass conservation equations between reactant and resultant in the catalytic chromatographic process was shown as below: where, 1 r k f − was the reactant reduction rate, and The concentration wave Equation (1) were a system of two nonlinear hyperbolic partial differential equations, which was a hard mathematical problem.But in some practical situations, the problem can be simplified [2].Assume 1 c was small, or the adsorption coefficient 1 b was small, that was, 1 1 1 b c  .While considering the assumed reaction rate r k is relatively minor, then 2 c was also small, that was, 2 1 c  , 2 2 1 b c  .In fact, in the quantitative analysis using high performance liquid chromatography (HPLC), the concentrations of most analytes, such as the reactant A and the resultant B here, were all very small [2] [3].Thus the adsorption isotherm above can be approximated as a linear and regarded as follows: and denoted concretely: they were positive constant, thus Equation ( 1) can be simplified to the following semi-coupled system of two linear hyperbolic partial differential equations.In which, the reactant concentration wave model was the initialboundary value problem of a self-closed hyperbolic partial differential equations, while the resultant concentration wave model was the initial boundary value problem of hyperbolic partial differential equations coupling reactant concentration.
Chromatographic process started from the boundary, and there were many types of the boundary conditions, such as the injection methods of δ-pulse, wide pulse, head-on, etc.; whose corresponding boundary condition were not zero.The initial state of chromatography columns were typically empty, that the initial conditions corresponding to 0. However, in practical problems, there were some important chromatograph whose corresponding initial conditions is not zero, such as simulated moving bed chromatography.Therefore, it is necessary to study the general initial-boundary value problem with both the initial and boundary values were not 0.That was, 1 2 , c c satisfied the following the general initial-boundary value problems.

Explicit Solution of Concentration Wave
Firstly, solved the initial-boundary value problem (6) for concentration wave of of reactant 1 c .According to Laplace transform of t, noted that: Then solved the ordinary differential Equation (8) about  ( ) 1 , c x p , we got: e e e , 0, , ) To sum ( 9) and (10) up, ( ) ( ) , e , 0 .
Then solved the initial-boundary value problem (7) for the concentration wave of resultant 2 c .Similarly, according to Laplace transform of t, noted that: The above problem (7) satisfied the following ordinary differential equation: Solved the ordinary differential Equation (12) about  ( ) 2 , c x p , we got: 0, , ) , d , 0 0, 0, To sum (13), ( 14) and (15) up, Using the expression (11) of 1 c and the relation Equation ( 16) of 1 c and 2 c , the explicit solution ex- pressions of 2 c were derived by dividing into the following three cases.In the case of 1 1 e , 0 , 1 e , .
In the case of 1 2 λ λ < , we had ( ) e e d , .
Particularly, when the initial-boundary problem ( 6) and ( 7) satisfied the following conditions ( ) ( ) ( ) the explicit solution of reactant and resultant concentration wave 1 2 , c c were obtained as follows.Following (11), we had ( ) ( ) x According to the expressions (17), ( 18) and (19), we had the explicit solution expressions of 2 c as follows.When 1

Solutions and Stability for δ-Pulse and Wide Pulse Injections
In this section, we derived the solutions of reactant and resultant concentration waves in wide pulse and δ-pulse injections detailedly.And the stability analysis between the resultant solutions of the two modes of pulse injection was further discussed.

δ-Pulse Injection
Chromatographic process started from the boundary, and there were many types of the boundary conditions, such as the methods of δ-pulse, wide pulse, head-on, etc; whose corresponding boundary condition was not zero.Where, δ-pulse and wide pulse were the most common way of chromatography injection method.Firstly, initial state of chromatography column in the δ-pulse method, which injection function was a kind of δ-function, was typically empty.So in the case of δ-Pulse, 1 c satisfied the following initial-boundary problem.
( ) ( ) ( ) where k is a constant represented the injection size, which is equal to 10 p c t in wide pulse method in Section 4.2.According to the behavior of the δ-function, we had ( ) The solution of concentration wave for reactant was obtained by Laplace transform as similar with Section 3.
The concentration wave corresponding to δ-pulse injection of reactant and resultant can be expressed as follows.
If there was no reaction terms, that was, 0 α = , we got ( ) ( ) As for the solution of concentration wave for resultant, the initial and boundary values were both 0. From the expression (21), ( 22) and (23), we had the explicit solution expressions of 2 c .When

Wide Pulse Injection
Wide pulse was the another most common way of chromatography injection method, its initial state of chromatography column was typically empty, so the initial condition was the follows, ( ) ( ) 0, 0.
The corresponding injection function was given as follows, ( ) ( ) where, p t was the injection time, 10 c was the injection rate, both of them are constant.In this paper, Wide pulse was taken as an another example, the solution of concentration wave for reactant and resultant were derived detailedly.
Similarly, we had the explicit solution expressions of 1 c and 2 c as follows, ( )

Stability Analysis between Wide Pulse and δ-Pulse Injections
Note that, the boundary condition in wide pulse injection tended to the condition in δ-pulse injection.We also showed that the mentioned limit relationship was still valid for the solutions in the two modes of pulse injection.
The main result of this work is the following theorem: Theorem 1.If = and when ( ) ∃ (a sufficiently small constant), so that which was converging to the solution (27) in δ-pulse injection.
In the case of In the case of For any fixed 0 x , p t ∃ (a sufficiently small constant), so that 0 3) When 1 2 λ λ < , the solution in wide pulse method (34) was equivalent to the following.
In the case of     39) and (43), we can conclude that this Theorem is true.□

Conclusion
The chromatography model with a reaction A B → was established and can be simplified a semi-coupled system of two linear hyperbolic PDE's in some practical situations.In which, the reactant concentration wave model was the initial-boundary value problem of a self-closed hyperbolic PDE, while the resultant concentration wave model was the initial-boundary value problem of hyperbolic PDE coupling reactant concentration.The general explicit expressions for the concentration wave of the reactants and resultants were derived by Laplace transform.The δ-pulse and wide pulse injections were taken as the examples to discuss detailedly, and it was proved that the continuous dependence of solutions was in accordance with the dependence under corresponding boundary conditions.It was significant for further analysis of chromatography in nonlinear case, optimizing chromatographic separation, determining the physical and chemical characters.

=
solution of concentration wave for resultant in wide pulse injection converges to the resultant solution in δ-pulse injection., from (32), For any fixed 0 x .