Study of Industrial Naphtha Catalytic Reforming Reactions via Modelling and Simulation

Aminu Zakari Yusuf 1,2, B. O. Aderemi 3, Raj Patel 1 and Iqbal M. Mujtaba 1,* 1 Chemical Engineering Department, Faculty of Engineering & Informatics, University of Bradford, Bradford BD7 1DP, UK; a.y.zakari@bradford.ac.uk (A.Z.Y.); R.Patel@bradford.ac.uk (R.P.) 2 Nigerian National Petroleum Cooperation (NNPC), Abuja 900211, Nigeria 3 Chemical Engineering Department, Ahmadu Bello University Zaria, Zaria 810107, Nigeria; benjaminaderemi@gmail.com * Correspondence: I.M.Mujtaba@bradford.ac.uk


Introduction
Naphtha is a low octane cut from atmospheric distillation of crude oil from the topping plant usually in the range between 30 • C-200 • C and constitutes typically 15-30% by weight of the crude oil [1].This cut contains some impurities such as sulphur, oxygen and nitrogen compounds.This low run naphtha, sometimes referred to as straight run, requires reforming to a more valuable higher run naphtha called reformate which is blended into a pool of gasoline for use as a motor spirit.SR naphtha can be Light naphtha a fraction with boiling point ranging from 30 • C to 90 • C, containing the C 5 and C 6 hydrocarbons or Heavy naphtha which is the fraction with boiling range from 90 • C to 200 • C. The medium naphtha is a fraction of this heavy cut that boils below 150 • C and includes mostly C 7 -C 9 hydrocarbons [2].Prior to reforming, the SR is treated in a hydrotreating unit.Reforming is a form of reaction where the components paraffins, naphthenes and aromatics are rearranged to produce higher gasolines.These components typically contain 40-70 wt% paraffins, 20-50 wt% naphthenes, 5-20 wt% aromatics [3].The various reforming reactions are dehydrogenation reactions of the paraffins to naphthenes and dehydrogenation of the naphthenes to aromatics, dehyrocylization of paraffins to naphthenes, isomerization of normal paraffins to iso-paraffins, cracking of paraffins and naphthenes and hydrodealkylation of aromatics.Modelling and simulation of catalytic reforming has been considered by different researchers to determine the plant performance or product distribution behaviour.Hu carried out [4] steady state modelling of semi regenerative catalytic reformer with 17 lumps and monitored the plant performance on stream.Modelling and parameter estimation was carried out by Stijepovic with benchmarking with industrial data and a good prediction of concentration and hydrogen was obtained [5].Optimization of catalyst distribution was carried by [6] to determine the optimal catalyst distribution in the reactors.The dynamics of a process tells a lot about the behaviour of the system and the effect of perturbation due to changes in operational conditions such as temperature.It gives an indication on the transient behaviour and the time for the system to reach steady state.This information helps whenever there is need for plant scale up or process modification.Few works have been carried out on the dynamics of the reaction [7][8][9][10].The modelling and simulation works carried out by [11,12] (steady state) and [10] (dynamic) use the same reactor configurations, same feed throughput and feed stock properties.It requires rigorous modelling method to model and simulate dynamically series of reactors concurrently together having three or four heaters and reactors in series.It is a difficult task and in [10], one reactor was simulated as a first contribution In this work, first a modelling and simulation of catalytic reforming unit by [11,12] was performed on gPROMs software and was compared with that performed on MATLAB ODE45 by the earlier authors.The validated model is then used to study the steady state and dynamic behaviour of an industrial commercial plant of Kaduna Refining and Petrochemical Company, using four reactors was performed sequentially using the gPROMs software (Version 5.0.0Process System Enterprise Limited, London, UK).gPROMs software was found capable and robust in performing complex dynamic modelling and simulations.

Process Description and Process Variables
Naphtha catalytic reforming reactions include Dehydrogenation, Dehydrocyclization, Isomerization and Hydrocracking reactions which are performed on both the metallic and acidic parts of the catalyst.A simple network reaction is given in Figure 1.
Processes 2019, 7, x FOR PEER REVIEW 2 of 29 been considered by different researchers to determine the plant performance or product distribution behaviour.Hu carried out [4] steady state modelling of semi regenerative catalytic reformer with 17 lumps and monitored the plant performance on stream.Modelling and parameter estimation was carried out by Stijepovic with benchmarking with industrial data and a good prediction of concentration and hydrogen was obtained [5].Optimization of catalyst distribution was carried by [6] to determine the optimal catalyst distribution in the reactors.The dynamics of a process tells a lot about the behaviour of the system and the effect of perturbation due to changes in operational conditions such as temperature.It gives an indication on the transient behaviour and the time for the system to reach steady state.This information helps whenever there is need for plant scale up or process modification.Few works have been carried out on the dynamics of the reaction [7][8][9][10].The modelling and simulation works carried out by [11,12] (steady state) and [10] (dynamic) use the same reactor configurations, same feed throughput and feed stock properties.It requires rigorous modelling method to model and simulate dynamically series of reactors concurrently together having three or four heaters and reactors in series.It is a difficult task and in [10], one reactor was simulated as a first contribution In this work, first a modelling and simulation of catalytic reforming unit by [11] and [12] was performed on gPROMs software and was compared with that performed on MATLAB ODE45 by the earlier authors.The validated model is then used to study the steady state and dynamic behaviour of an industrial commercial plant of Kaduna Refining and Petrochemical Company, using four reactors was performed sequentially using the gPROMs software (Version 5.0.0Process System Enterprise Limited, London, UK).gPROMs software was found capable and robust in performing complex dynamic modelling and simulations.

Process Description and Process Variables
Naphtha catalytic reforming reactions include Dehydrogenation, Dehydrocyclization, Isomerization and Hydrocracking reactions which are performed on both the metallic and acidic parts of the catalyst.A simple network reaction is given in Figure 1.M represents the metal catalyst site, A represents the Acid catalyst site.The main reactions involved in reforming processes are catalysed essentially either by the acid support or the metal functions.

Metallic function Dehydrocyclization
Metallic + acidic functions Isomerization Metallic + acidic functions Hydrogenolysis Metallic function

Hydrocracking Metallic + acidic function
A Naphtha catalytic reforming unit (CRU) typically is made up of a preheat system, that is, a long vertical heat exchanger with large surface area to allow for feed preheating prior to further heating in the charge heater usually three or four heaters and reactors in series and a low pressure M represents the metal catalyst site, A represents the Acid catalyst site.The main reactions involved in reforming processes are catalysed essentially either by the acid support or the metal functions.A Naphtha catalytic reforming unit (CRU) typically is made up of a preheat system, that is, a long vertical heat exchanger with large surface area to allow for feed preheating prior to further heating in the charge heater usually three or four heaters and reactors in series and a low pressure separator.The hydrogen is flashed and separated in the separator where part of it is recycled to the meet the feed again continuously before it enters the first heater, while the heavier liquid (reformate) is sent to a stabilization column.Reformate is obtained as a bottom product from the stabilization column while the liquefied petroleum gases (LPG) are obtained from the top of the stabilization after condensation in a trim cooler and condenser.It is important to mention that the flashed hydrogen from the low pressure separator goes through a drying section, a dryer consisting some alumina balls in order to reduce the amount of water vapour in the gas to a maximum of 25 ppm before recycling back to the feed.Also the bottom of the low pressure separator, liquid reformate prior to stabilization is re-contacted with hydrogen from a dry gas compressor at high pressure (high pressure separator) in order to strip more of hydrogen if still contained in the reformate.This high pressure separator otherwise called re-contactor is the source of hydrogen to other units in refinery where it is distributed at high pressure.These reforming reactions globally are endothermic in nature, multiple charge heaters are used to raise the temperature of the reaction to the desired value (400-500 • C) [13].During the reaction in the catalytic bed in the reactor, one of the most significant reaction which is the dehydrogenation of naphthenes to aromatics, is very fast and more endothermic.This results in a high decrease of temperature in the reactor.The temperature of the exit stream from the first reactor is raised and sent to the next reactor.The rate of reaction of the feed decreases whenever it moves through the reactors in series.The reaction becomes less endothermic and the differential temperature along the reactor falls and the quantity of heat required to reheat the stream also decrease.The treated naphtha fed to the reactor is combined with a gas stream containing about 60 to 90 mole % hydrogen whose pressure depends on compressor capacity.Thus, total reactor charge is the naphtha feed plus the hydrogen recycle gas stream.The reactor feed must be raised to the proper temperature for the reforming reactions to occur when the charge contacts the catalyst.As shown in Figure 2 total reactor charge is heated, at first by exchange with effluent from the last reactor and is finally brought up to the first reactor inlet temperature in the first charge heater which is also called the first inter heater.The reactor effluent-to-feed exchanger recovers the heat from the reactor effluent and provides it to the reactor feed.Thus, it is one key to energy conservation in a catalytic reformer.The reactor effluent which may be as high as 750 to 790 K, must be cooled to 315 to 320 K for flash separation of hydrogen from reformate.This heat exchange is carried out in several banks of heat exchangers, arranged for parallel flow.Some of the newer units use one large single-pass vertical exchanger for this purpose.Catalytic reforming is a gas phase process [13].After passing through the reactor effluent-to-feed exchanger and the charge heater, the total reactor charge is 100% vapour at the appropriate reaction temperature and is ready to contact the reforming catalyst.The most commonly used catalyst type is platinum on alumina support.Recent developments in catalysis have produced some bi-metallic and tri-metallic catalysts, such as Pt-Re or Pt-Ir on alumina support [14].The flow scheme in Figure 2 shows the four reactors in series.
The hydrogen circulation to the naphtha charge is called the hydrogen recycle.The hydrogen recycle is maintained in order of control the rates of deactivation of the catalyst.Excess hydrogen from the separator is sent to feed gas or to hydrogen consuming units in the refinery such as hydro-treaters and petrochemical plants.The separator liquid, comprised mostly of the desired reformate product but also containing hydrogen, methane, ethane, propane and butanes is pumped to the reformate stabilizer.The stabilized reformate is sent to storage for gasoline blending while the C 3 and C 4 are sent to LPG storage tanks.The reformer that is modelled in this work is a semi regenerative type of reformer consisting of four reactors in series.When the catalyst gets deactivated after time and has reached end of run (EOR) where the quality of the product is not obtained, the plant is shut down and the deactivated catalyst is regenerated and this is called a catalyst cycle.A cycle may be terminated for a number of reasons but the most common is poor catalyst performance due to deactivation of the catalyst by coke formation on its surface as a result of severity of operation mainly high temperature.This usually leads to a loss of reformate yield.The regeneration procedure restores the catalyst to near fresh-catalyst condition.The temperature, pressure, space velocity, H 2 /Hc ratio and feed quality are the variables affecting rate of the reaction and the sensitivity of the product distribution to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table 1 [15].
distribution to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in

Process Classification
Naphtha catalytic reforming process is typically classified based on the method and procedure of operation and regeneration of catalyst.They are semi catalyst regeneration, cyclic catalyst regeneration and continuous catalyst regeneration.Different refineries have different processes but most widely used is the semi catalyst regeneration, then continuous catalyst regeneration and the less commonly used is cyclic catalyst regeneration.Nowadays, more catalytic reformers are developed and designed for continuous catalyst regeneration and the previous semi catalyst regeneration units are being revamped to operate in continuous regeneration mode.A semi catalyst regeneration process typically consists of three or four reactors in series in a fixed bed catalyst system and operates within a time line where the catalyst gets deactivated, a cycle length and it is regenerated.Usually a period of 6 months to one year.The advantage of this process type is the simplicity of the mode of operation but one major drawback is that the plant has to be shut down for catalyst regeneration after end of run (EOR).Another drawback is the application of high pressure and hydrogen hydrocarbon (HHR) to avoid quick catalyst deactivation.The cyclic catalyst regeneration unit is similar to the semi catalyst regeneration mode but it has an extra swing reactor which is used as a spare reactor whenever the catalyst of one of the reactors is deactivated.It serves as a reactor for the unit while it is regenerated and this prevents the total shut down of the plant.The drawback of this process is the complexity and challenges of reactor swing from one to the other although it has the advantage of use of lower pressure than the semi catalyst regeneration mode.The continuous catalyst regeneration mode is a moving bed mode where the catalyst is continuously regenerated.It does not require shut down for catalyst regeneration and this implies less down time for the unit and more production time.The major challenge of this process is the complexity of mechanism of engaging and disengaging of the product from the last reactor to the regenerator.The hydrogen from the catalyst to be regenerated has to be displaced with nitrogen before it is sent to the regenerator to avoid explosion since the regenerator operates at a high temperature of about 600-700 • C. The catalyst after cooling must be swept of nitrogen before is sent to the reactor.An important aspect in this type of process is disharmony at times between the catalyst suppliers and the technology owners where some catalysts used may not be suitable for the technology of the process.The major operational draw backs of the previous modes have solved in this mode of operation.It uses low pressure and low hydrogen hydrocarbon ratio and gives higher hydrogen purity and higher yield of aromatics.In this work, a semi catalyst regeneration mode of CRU, in Kaduna refining and petrochemical company is modelled.

Kinetic Models
In this work, various models have been developed for both steady state and dynamic process of the catalytic reforming reaction with different lumps to investigate the behaviour of the reaction and products distribution.The reaction steps and equations of the model are given in this section.Table 2 shows the reaction steps and rate constants for the various reactions used in the model.The kinetic rate equations for the components are given in Equations ( 1)-( 25) where SV is space velocity and 1  sv is residence time.
where a where K i is the kinetic constants for the reactions, T 0 and P 0 are the reference temperatures and pressures.E Ai (activation energy) and w (exponential values of pressure) are taken from [10] and are given in Table 3. k i , E Ai , w, T 0 and P 0 are also obtained from [10].

Mathematical Model
The kinetic rate equations from Equations ( 1)-( 25) which determine the rate of the reaction of each component is used in the mathematical modelling Equations of ( 27)-(32) to determine the concentration, temperature and pressure distributions of the components describing the system.The rates of the equations are defined as r j in Equations ( 27) and (28) and r j γ i summations of the rate equations and their stoichiometric coefficients.These model equations, (27)-(32), in cooperated with Equations (1)-( 25) are solved simultaneously on gPROMs.Equations ( 27) and (28) are equations representing the dynamic behaviour of the system and change with time dt.At steady state, the system does not change with time and dt becomes zero.So the equations can be used for both steady state and dynamic behaviour depending on the application.

RON Model
The estimation of the Research octane number, RON of hydrocarbons for the feed and products can be done using different methods of prediction.RON can be calculated for each pure component by using a polynomial equation that correlates to the normal boiling point [16].Equation (33) represents the model equation for the RON estimation.RON = Aa + BbT + CcT 2 + DdT 3 + EeT 4  (33) where T = TBP × 0.01, TBP represents the normal boiling point ( • C) and Aa, Bb, Cc, Dd and Ee, are coefficients.The RON of a hydrocarbon mixture is calculated by assuming that the mixture consists of paraffins, naphthenic hydrocarbons and aromatics.The equation is expressed as the sum of the RON for each pure component multiplied by the volume fractions of the components.Equation ( 27) is the material balance equation where ρb is the bulk density and ε the bed voidage.Equation ( 28) is the heat balance of the reaction to determine the temperature behaviour where −∆H Rj is the heat of the reaction of jth component and Cpcat the heat capacity of the catalyst.Equation ( 29) is the pressure profile of the reaction where Dp is particle diameter, G the mass flux, ε the bed voidage, ρ the gas density and A the reactor cross sectional area.H 0 ri , the standard heat of formation of the components.The constants A, B, C, D used for the heat capacity, Cp, in Equations ( 31) and (32) are taken from [16].These equations are solved simultaneously using gPROMs to determine the behaviour of the system.gPROMs is a robust mathematical software that can solve the equations for all the four reactors in series dynamically given the behaviour of the paraffins, naphthenes and aromatics as well as the dynamics of the temperature in the reactors.All the components behaviours were determined.

Model Validation
The model validation was carried out by modelling and simulating the commercial data of [11,12] on gPROMs to ascertain the capability and ruggedness of the mathematical software in complex modelling.The modelling of [11] was performed using MATLAB software and the result is compared with that obtained with gPROMs.The configuration of the commercial reformer and Feed Properties are given in Table 4; Table 5 respectively.The commercial reformer of [11] was simulated with the properties of Tables 4 and 5 respectively with four semi regenerative reactors in series.The reformer's through put is 30 MBPD at inlet temperature of 495 • C, pressure of 10.5 kg/cm 2 using MATLAB mathematical tool on ODE45.
The simulation was performed on gPROMs software to validate its capability in modelling this problem.The result obtained from gPROMs is compared with that obtained with MATLAB from [11] and are shown in Table 6.The result showed a good and comparable result between the actual and that simulated with gPROMs.This is an indication of the capability of the mathematical tool in solving complex model equations and hence can be used as a validation for modelling other catalytic reformers using same kinetic models.

Results and Discussion
The model was simulated based on the kinetic data from [11,12] and validated with Kaduna Refining & Petrochemical Company, KRPC Commercial Catalytic Reformer data from Nigeria.Tables 7  and 8 show the Configuration and Feed Properties of KRPC Reformer.The results show a good agreement as shown in Figure 3 and Table 9.

Steady State Simulation Results
The time derivatives of the model equation are set to zero to simulate the process at steady state.The modelling results from the model equations of ( 27)-(32) gives the concentration and temperature distribution only across the reactor height.In this work, four reactors are simulated with the reactor configurations given in Table 4.The total reactor height is 25.232 m which is total sum of the heights of the four reactors as in Table 4 and is denoted as height of reactor in Section 7.1 therefore the height in Figures 4-9 refer to reactors total height of 25.232 m.In Section 7.2, the different reactors heights for the four reactors were used differently to show variation with both the height and time as shown in Figures 10-14.
increases, from Table 2 we can see that the rates increase from 0.0, 0.58, 1.33,1.81,2.54 from P6 to P11 however the tendency of paraffins to hydrocrack increases concurrently.Kinetically, the rate of dehydrocyclization increases with low pressure and high temperature.To sum up, the dehydrocyclization of P6 paraffins to benzene is more difficult than that of C7 paraffin to toluene, which itself is more difficult than that of C8 paraffin to xylenes.Accordingly, the most suitable fraction to feed a reforming process is the C7-C10 fraction.Figure 5 shows the behaviour of the naphthenes along the reactor height.N6, N7, N8, N9, N10 are naphthenes with hydrocarbon number 6, 7, 8, 9, 10 respectively.The fastest reaction is aromatization that is, dehydrogenation of naphthenes, resulting in the large temperature drop due to its highly endothermic nature.The drastic reduction is due to the aromatization to aromatics.Thermodynamically the reaction is highly endothermic and is favoured by high temperature and low pressure.In addition, the higher the number of carbon atoms the higher the aromatics production at equilibrium from N8, N9, N10.This can be seen from Table 2 in equations ( 34)-( 39) where the reaction rates increase from 4.02, 9.03, 21.5, 24.5 and 24.5 respectively for N6, N7, N8, N9 and N10.From a kinetic point of view, the rate of reaction increases with temperature.The naphthenes also undergo hydrocracking to lighter hydrocarbons leading to their decrease as shown in reactions 49-57 in Table 2. Figure 6 shows the behaviour of aromatics along the reactor height.A6, A7, A8, A9, A10 are aromatics with hydrocarbon number 5, 6, 7, 8, 9, 10 respectively.The dehydrogenation of naphthenes as seen from Figure 5 increases the aromatics shown in Figure 6.The sharp increase is as a result of drastic aromatization of the naphthenes favoured by the metallic sites of the catalyst.The rate of benzene formation is lower due to lower carbon number of N6 while A7, A8, A9 increase rapidly due to higher carbon number of N7, N8 and N9.This can be seen from reactions 58-71 in Table 2. Figure 7 shows the temperature profile of the reactions along the reactor height.The sharp drop in first and second reactor is due to the more endothermic dehydrogenation reaction of naphthenes to aromatics.Typically, dehydrogenation and isomerization reactions take place in the first reactor, dehydrogenation, isomerization, dehydrogenation and cracking in the second followed by dehydrogenation and cracking in the third and fourth reactors.Figure 6 shows the behaviour of aromatics along the reactor height.A6, A7, A8, A9, A10 are aromatics with hydrocarbon number 5, 6, 7, 8, 9, 10 respectively.The dehydrogenation of naphthenes as seen from Figure 5 increases the aromatics shown in Figure 6.The sharp increase is as a result of drastic aromatization of the naphthenes favoured by the metallic sites of the catalyst.The rate of benzene formation is lower due to lower carbon number of N6 while A7, A8, A9 increase rapidly due to higher carbon number of N7, N8 and N9.This can be seen from reactions 58-71 in Table 2. Figure 7 shows the temperature profile of the reactions along the reactor height.The sharp drop in first and second reactor is due to the more endothermic dehydrogenation reaction of naphthenes to aromatics.Typically, dehydrogenation and isomerization reactions take place in the first reactor, dehydrogenation, isomerization, dehydrogenation and cracking in the second followed by dehydrogenation and cracking in the third and fourth reactors.Figure 4 shows the concentration profile of paraffins along the reactor height.P 5 , P 6 , P 7 , P 8 , P 9 , P 10 are paraffins with hydrocarbon number 5, 6, 7, 8, 9, 10 respectively.Normal Paraffin components basically undergo three major reactions during naphtha reforming.They undergo hydrocracking to lighter paraffins that is, methane, ethane, propane and butane as in shown in reactions 8 to 33 in Table 2, isomerization to iso-paraffins which are not modelled in this work, it is a slow reaction with slow reaction rate, dehydrocyclization to naphthenes, a very slow reaction which leads to a decrease in the paraffins.These reactions are reactions 1-7 from Table 2. Dehydrocyclization reaction becomes easier as the molecular weight of the paraffins increases as in P 8 , P 9 and P 10 while P 7 shows little increase.P 5 and P 6 increase due to hydrocracking.The main effects of hydrocracking are decrease of paraffins (C 5 +) in the reformate, decrease in hydrogen production, increase of LPG production and hydrogenolysis.The isomerization reactions are fast, slightly exothermic and do not affect the number of carbon atoms.The thermodynamic equilibrium of iso-paraffins to paraffins depends mainly on the temperature and pressure has no effect.The paraffins isomerization results in a slight increase of the octane number.These reactions are promoted by the acidic function of the catalyst support.The paraffin dehydrocyclization step becomes easier as the molecular weight of the paraffin increases, from Table 2 we can see that the rates increase from 0.0, 0.58, 1.33,1.81,2.54 from P 6 to P 11 however the tendency of paraffins to hydrocrack increases concurrently.Kinetically, the rate of dehydrocyclization increases with low pressure and high temperature.To sum up, the dehydrocyclization of P 6 paraffins to benzene is more difficult than that of C 7 paraffin to toluene, which itself is more difficult than that of C 8 paraffin to xylenes.Accordingly, the most suitable fraction to feed a reforming process is the C 7 -C 10 fraction.Figure 8 shows the increase in hydrogen yield along the height of the reactor.Figure 9 shows the increase in RON along the reactor height.This is due to the increase in the aromatics along the reactor height.The RON increases from reactor one to four due to increase in aromatization reaction  Figure 8 shows the increase in hydrogen yield along the height of the reactor.Figure 9 shows the increase in RON along the reactor height.This is due to the increase in the aromatics along the reactor height.The RON increases from reactor one to four due to increase in aromatization reaction which gives a higher research octane number, a parameter for antiknock in the gasoline engine.Straight run naphtha has a low RON and cannot be used as gasoline and hence the reforming reaction of the components to give a gasoline with higher research octane number.and isomerization reactions continue, thereby increasing slightly the temperature with time.The temperature here reaches steady state at 0.12 h which is longer than the first value.In reactor three there is lesser temperature drop due to exothermic reactions of isomerization and hydrocracking, reaching a steady state at 0.20 h.Reactor four has the highest exothermic reactions and one can see a sharp rise in temperature.It reaches a steady state at 0.25 h.Figures [15][16][17][18][19][20] give the dynamic behaviours of the paraffins.Unlike in Figure 4 where it shows the steady state behaviour of the paraffins changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4.  Figures [15][16][17][18][19][20] give the dynamic behaviours of the paraffins.Unlike in Figure 4 where it shows the steady state behaviour of the paraffins changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4. Figure 5 the behaviour of the naphthenes along the reactor height.N 6 , N 7 , N 8 , N 9 , N 10 are naphthenes with hydrocarbon number 6, 7, 8, 9, 10 respectively.The fastest reaction is aromatization that is, dehydrogenation of naphthenes, resulting in the large temperature drop due to its highly endothermic nature.The drastic reduction is due to the aromatization to aromatics.Thermodynamically the reaction is highly endothermic and is favoured by high temperature and low pressure.In addition, the higher the number of carbon atoms the higher the aromatics production at equilibrium from N 8 , N 9 , N 10 .This can be seen from Table 2 in equations (34)-(39) where the reaction rates increase from 4.02, 9.03, 21.5, 24.5 and 24.5 respectively for N 6 , N 7 , N 8 , N 9 and N 10 .From a kinetic point of view, the rate of reaction increases with temperature.The naphthenes also undergo hydrocracking to lighter hydrocarbons leading to their decrease as shown in reactions 49-57 in Table 2.
Figure 6 shows the behaviour of aromatics along the reactor height.A 6 , A 7 , A 8 , A 9 , A 10 are aromatics with hydrocarbon number 5, 6, 7, 8, 9, 10 respectively.The dehydrogenation of naphthenes as seen from Figure 5 increases the aromatics shown in Figure 6.The sharp increase is as a result of drastic aromatization of the naphthenes favoured by the metallic sites of the catalyst.The rate of benzene formation is lower due to lower carbon number of N 6 while A 7 , A 8 , A 9 increase rapidly due to higher carbon number of N 7 , N 8 and N 9 .This can be seen from reactions 58-71 in Table 2.
Figure 7 shows the temperature profile of the reactions along the reactor height.The sharp drop in first and second reactor is due to the more endothermic dehydrogenation reaction of naphthenes to aromatics.Typically, dehydrogenation and isomerization reactions take place in the first reactor, dehydrogenation, isomerization, dehydrogenation and cracking in the second followed by dehydrogenation and cracking in the third and fourth reactors.
Figure 8 shows the increase in hydrogen yield along the height of the reactor.Figure 9 shows the increase in RON along the reactor height.This is due to the increase in the aromatics along the reactor height.The RON increases from reactor one to four due to increase in aromatization reaction which gives a higher research octane number, a parameter for antiknock in the gasoline engine.Straight run naphtha has a low RON and cannot be used as gasoline and hence the reforming reaction of the components to give a gasoline with higher research octane number.

Dynamic Simulation of the Four Reactors
The modelling results from the model Equations of ( 27)-(32) gives the concentration and temperature distributions only across the reactor height and the time simultaneously where their behaviours are predicted with changing both time height.In this work, four reactors are simulated with the reactor configurations given in Table 5.
Figures 10-13 shows the 3D profiles of the dynamics of the temperatures in the four reactors respectively.It is a combination of the change along both the reactor height and time.For the variation with time, the interesting aspect could be seeing clearly in Figure 14 that the exit temperatures falls and rise with time before reaching a final temperature value unlike the steady state where it shows only a drop in temperature without seeing such interesting aspect as in Figure 7. From Figure 14, temperature decreases along the height simultaneously with time.The dynamic simulation was performed over 10 h but the time to reach steady state was 0.25 h maximum, hence the representation in Figure 14 stops at 0.4 h for clarity purpose.The drop in temperature in first reactor is mainly due to the high endothermic reaction of naphthenes conversion to aromatics.Dehydrogenation is the main reaction in this reactor, the fastest reaction among all the reactions and reaches a steady state at about 0.05 h.This can also be attributed to the lowest residence time of 0.179 h for the reactions in the first reactor.There are also some endothermic reactions in the second reactor where dehydrogenation and isomerization reactions continue, thereby increasing slightly the temperature with time.The temperature here reaches steady state at 0.12 h which is longer than the first value.In reactor three there is lesser temperature drop due to exothermic reactions of isomerization and hydrocracking, reaching a steady state at 0.20 h.Reactor four has the highest exothermic reactions and one can see a sharp rise in temperature.It reaches a steady state at 0.25 h.
Figures 15-20 give the dynamic behaviours of the paraffins.Unlike in Figure 4 where it shows the steady state behaviour of the paraffins changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4.      Figures 21-25 give the dynamic behaviours of the naphthenes.Unlike in Figure 5 where it shows the steady state behaviour of the naphthenes changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state  Figures 21-25 give the dynamic behaviours of the naphthenes.Unlike in Figure 5 where it shows the steady state behaviour of the naphthenes changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state Figures 21-25 give the dynamic behaviours of the naphthenes.Unlike in Figure 5 where it shows the steady state behaviour of the naphthenes changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4.          Figures 26-30 give the dynamic behaviours of the aromatics.Unlike in Figure 6 where it shows the steady state behaviour of the aromatics changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4.   Figures 26-30 give the dynamic behaviours of the aromatics.Unlike in Figure 6 where it shows the steady state behaviour of the aromatics changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4. Figures 26-30 give the dynamic behaviours of the aromatics.Unlike in Figure 6 where it shows the steady state behaviour of the aromatics changing with the reactor height, in these Figures, the time to reach steady state in all the reactors where studied.The dynamics were run for 10 h but for clarity purpose it stops at 1 h.It takes 0.18 h in reactor 1 to attain steady state, 0.20 h to attain steady state in reactor 2, 0.22 h to attain steady state for reactor 3 while it takes 0.28 h to attain steady state in reactor 4.

Conclusion
In this paper, the catalytic naphtha reforming process of Kaduna refining and petrochemical company, Kaduna were studied via steady state and dynamic modelling using gPROMs software.The basic model of the process is taken from the literature and was developed in gPROMs platform and was validated using literature data.Note, while [10] simulated one of four reactors (in series) this

Figure 4 .
Figure 4. Concentration profile of paraffins with total height of reactors.

Figure 4 .Figure 5 .
Figure 4. Concentration profile of paraffins with total height of reactors.Processes 2019, 7, x FOR PEER REVIEW 14 of 29

Figure 5 . 29 Figure 6 .
Figure 5. Concentration profile of naphthenes with total height of reactors.

Figure 7 .
Figure 7. Temperature profile in the four reactors.

Figure 6 .
Figure 6.Concentration profile of aromatics with total height of reactors.

Figure 6 .
Figure 6.Concentration profile of aromatics with total height of reactors.

Figure 7 .
Figure 7. Temperature profile in the four reactors.

Figure 7 .
Figure 7. Temperature profile in the four reactors.

Figure 9 .
Figure 9. RON along the reactor height.

Figure 9 .
Figure 9. RON along the reactor height.

Figure 9 .
Figure 9. RON along the reactor height.

Figure 10 .
Figure 10.A 3D representation of Reactor 1 Temperature variation with reactor height and time.

FigureFigure 11 .
Figure 3D representation of Reactor 1 Temperature variation with reactor height and time.Processes 2019, 7, x FOR PEER REVIEW 18 of 29

Figure 12 .Figure 11 .
Figure 12.A 3D representation of Reactor 3 Temperature variation with reactor height and time.

Figure 12 .
Figure 12.A 3D representation of Reactor 3 Temperature variation with reactor height and time.A 3D representation of Reactor 3 Temperature variation with reactor height and time.

Figure 13 .
Figure 13.A 3D representation of Reactor 4 Temperature variation with reactor height and time.

Figure 14 .
Figure 14.Dynamic behaviour of Temperature in the four reactors.

Figure 13 .Figure 13 .
Figure 13.A 3D representation of Reactor 4 Temperature variation with reactor height and time.

Figure 14 .
Figure 14.Dynamic behaviour of Temperature in the four reactors.

Figure 14 .
Figure 14.Dynamic behaviour of Temperature in the four reactors.

Figure 15 .
Figure 15.Dynamic behaviour of P5 in the four reactors.

Figure 15 .
Figure 15.Dynamic behaviour of P5 in the four reactors.

Figure 16 .
Figure 16.Dynamic behaviour of P6 in the four reactors.

Figure 20 .
Figure 20.Dynamic behaviour of P10 in the four reactors.

Figure 19 .Figure 19 .
Figure 19.Dynamic behaviour of P 9 in the four reactors.

Figure 20 .
Figure 20.Dynamic behaviour of P10 in the four reactors.

Figure 20 .
Figure 20.Dynamic behaviour of P 10 in the four reactors.

Figure 22 .
Figure 22.Dynamic behaviour of N7 in the four reactors.

Figure 21 .Figure 21 .
Figure 21.Dynamic behaviour of N 6 in the four reactors.

Figure 22 .
Figure 22.Dynamic behaviour of N7 in the four reactors.

Figure 22 .Figure 21 .
Figure 22.Dynamic behaviour of N 7 in the four reactors.

Figure 22 .
Figure 22.Dynamic behaviour of N7 in the four reactors.

Figure 23 .
Figure 23.Dynamic behaviour of N 8 in the four reactors.

Figure 23 .
Figure 23.Dynamic behaviour of N8 in the four reactors.

Figure 24 .
Figure 24.Dynamic behaviour of N9 in the four reactors.

Figure 25 .
Figure 25.Dynamic behaviour of N10 in the four reactors.

Figure 24 . 29 Figure 23 .
Figure 24.Dynamic behaviour of N 9 in the four reactors.

Figure 24 .
Figure 24.Dynamic behaviour of N9 in the four reactors.

Figure 25 .
Figure 25.Dynamic behaviour of N10 in the four reactors.

Figure 25 .
Figure 25.Dynamic behaviour of N 10 in the four reactors.

Figure 26 .
Figure 26.Dynamic behaviour of A6 in the four reactors.

Figure 29 .
Figure 29.Dynamic behaviour of A9 in the four reactors.

Figure 28 .Figure 28 .
Figure 28.Dynamic behaviour of A 8 in the four reactors.

Figure 29 .
Figure 29.Dynamic behaviour of A9 in the four reactors.

Figure 29 .
Figure 29.Dynamic behaviour of A 9 in the four reactors.

Figure 30 .
Figure 30.Dynamic behaviour of A10 in the four reactors.

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (

Table 1 .
Effect of variables on product distribution, coke and Research o

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (

Table 1 .
Effect of variables on product distribution, coke and Research o

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON). to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15]. distribution

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).
[15]esses 2019, 7, x FOR PEER REVIEW distribution to these variables are studied in the result section of this work.A summary o of these variables are shown in Table1[15].

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RO [15]esses 2019, 7, x FOR PEER REVIEW distribution to these variables are studied in the result section of this wo of these variables are shown in Table1[15].

Table 1 .
Effect of variables on product distribution, coke and Research o [15]/Hc ratio Processes 2019, 7, x FOR PEER REVIEW 4 of 29distribution to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table 1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON). to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15]. distribution

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

,
[15]FOR PEER REVIEW distribution to these variables are studied in the result section of this work.A summary o of these variables are shown in Table1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RO to these variables are studied in the result section of this wo of these variables are shown in Table1[15]. distribution

Table 1 .
[15]ct of variables on product distribution, coke and Research o to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15]. distribution

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON). to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15]. distribution

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RO

Table 1 .
[15]ct of variables on product distribution, coke and Research o variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON).variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON). to these variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15]. distribution

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RON).

Table 1 .
[15]ct of variables on product distribution, coke and Research o variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON).variables are studied in the result section of this work.A summary of the effect of these variables are shown in Table1[15].

Table 1 .
[15]ct of variables on product distribution, coke and Research octane number (RON).Processes 2019, 7, x FOR PEER REVIEW distribution to these variables are studied in the result section of this work.A summary o of these variables are shown in Table1[15].

Table 1 .
Effect of variables on product distribution, coke and Research octane number (RO [15]esses 2019, 7, x FOR PEER REVIEW distribution to these variables are studied in the result section of this wo of these variables are shown in Table1[15].

Table 1 .
Effect of variables on product distribution, coke and Research o

Table 2 .
[10]e of reaction steps and rate constants of the Model (adapted from[10]).

Table 8 .
Feed Properties KRPC Commercial Catalytic Reformer.Comparison between Model Simulation and Industrial Plant data of KRPC.

Table 9 .
Relative errors between Industrial and Simulated results of KRPC plant.

Table 9 .
Relative errors between Industrial and Simulated results of KRPC plant.
Dynamic behaviour of P 6 in the four reactors.Dynamic behaviour of P7 in the four reactors.Dynamic behaviour of P 8 in the four reactors.
Figure 16.Figure 17.Dynamic behaviour of P7 in the four reactors.8 Figure 17.Dynamic behaviour of P 7 in the four reactors.