Nothing Special   »   [go: up one dir, main page]

WO2006075085A2 - Optimisation d’une reaction chimique dans un reacteur ouvert du type a plaques - Google Patents

Optimisation d’une reaction chimique dans un reacteur ouvert du type a plaques Download PDF

Info

Publication number
WO2006075085A2
WO2006075085A2 PCT/FR2006/000049 FR2006000049W WO2006075085A2 WO 2006075085 A2 WO2006075085 A2 WO 2006075085A2 FR 2006000049 W FR2006000049 W FR 2006000049W WO 2006075085 A2 WO2006075085 A2 WO 2006075085A2
Authority
WO
WIPO (PCT)
Prior art keywords
reaction
reactor
fluid
medium
reaction chamber
Prior art date
Application number
PCT/FR2006/000049
Other languages
English (en)
Other versions
WO2006075085A3 (fr
Inventor
Sébastien ELGUE
Fabrice Chopard
Michel Roger Cabassud
Patrick Alain Cognet
Laurent Emile Georges Prat
Christophe Jacques Jean Gourdon
Original Assignee
Alfa Laval Vicarb
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Alfa Laval Vicarb filed Critical Alfa Laval Vicarb
Priority to MX2007008567A priority Critical patent/MX2007008567A/es
Priority to CA2594822A priority patent/CA2594822C/fr
Priority to CN2006800022907A priority patent/CN101288078B/zh
Priority to JP2007550813A priority patent/JP5107054B2/ja
Priority to AU2006205731A priority patent/AU2006205731B2/en
Priority to BRPI0606194-0A priority patent/BRPI0606194A2/pt
Priority to EP06709060A priority patent/EP1836627A2/fr
Publication of WO2006075085A2 publication Critical patent/WO2006075085A2/fr
Priority to US11/773,667 priority patent/US8065038B2/en
Priority to NO20074160A priority patent/NO20074160L/no
Publication of WO2006075085A3 publication Critical patent/WO2006075085A3/fr

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J19/0006—Controlling or regulating processes
    • B01J19/0033—Optimalisation processes, i.e. processes with adaptive control systems
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J19/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J19/24—Stationary reactors without moving elements inside
    • B01J19/248—Reactors comprising multiple separated flow channels
    • B01J19/249—Plate-type reactors
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2219/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J2219/24—Stationary reactors without moving elements inside
    • B01J2219/2401—Reactors comprising multiple separate flow channels
    • B01J2219/245—Plate-type reactors
    • B01J2219/2451—Geometry of the reactor
    • B01J2219/2453—Plates arranged in parallel
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2219/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J2219/24—Stationary reactors without moving elements inside
    • B01J2219/2401—Reactors comprising multiple separate flow channels
    • B01J2219/245—Plate-type reactors
    • B01J2219/2461—Heat exchange aspects
    • B01J2219/2462—Heat exchange aspects the reactants being in indirect heat exchange with a non reacting heat exchange medium
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
    • B01J—CHEMICAL OR PHYSICAL PROCESSES, e.g. CATALYSIS OR COLLOID CHEMISTRY; THEIR RELEVANT APPARATUS
    • B01J2219/00—Chemical, physical or physico-chemical processes in general; Their relevant apparatus
    • B01J2219/24—Stationary reactors without moving elements inside
    • B01J2219/2401—Reactors comprising multiple separate flow channels
    • B01J2219/245—Plate-type reactors
    • B01J2219/2474—Mixing means, e.g. fins or baffles attached to the plates
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F28—HEAT EXCHANGE IN GENERAL
    • F28D—HEAT-EXCHANGE APPARATUS, NOT PROVIDED FOR IN ANOTHER SUBCLASS, IN WHICH THE HEAT-EXCHANGE MEDIA DO NOT COME INTO DIRECT CONTACT
    • F28D9/00—Heat-exchange apparatus having stationary plate-like or laminated conduit assemblies for both heat-exchange media, the media being in contact with different sides of a conduit wall
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/10—Analysis or design of chemical reactions, syntheses or processes

Definitions

  • the present invention relates to a method for optimizing a chemical reaction in an open reactor of the plate type consisting of a stack of plates delimiting between them at least one block comprising a reaction chamber and two lateral chambers of circulation of a utility fluid for heat exchange.
  • a major disadvantage of these reactors is the difficulty of evacuating the heat produced by a reaction, which limits the possibilities of control and optimization of the reaction.
  • the present invention is intended in particular to provide a simple and effective solution to the control and optimization of chemical reactions through the optimization of an open-type reactor that does not have the aforementioned drawbacks of batch reactors or reactors. type
  • an open reactor consisting of a stack of plates delimiting between them at least one block comprising a reaction chamber formed between two lateral heat exchange chambers with heat exchange fluid exchange utility, feed means of the reaction chamber by a continuous flow of one or more reagents and means supplying the two lateral chambers with a continuous flow of utility fluid;
  • this dynamic model comprising balance equations matter and energy and constraints
  • a continuous or open reactor comprising a superposition of plates defining at least one reaction chamber between two fluid circulation chambers, solves the problems of heat dissipation found in discontinuous or closed reactors of the type " batch ", by appropriate selection of the temperatures and flows of the utility fluid or fluids, which circulate on both sides of the reaction chamber, and also by the selection of the direction of passage of these fluids (co-current, at countercurrent or cross flow) relative to the direction of passage of the reaction fluid in the reaction chamber.
  • the dynamic reactor model includes material and energy balance equations and stress equations, which are characteristic of the reactor and the chemical reaction.
  • the modeling of the reaction chamber comprises a decomposition of this chamber into successive elementary cells each containing a perfectly stirred fluid medium.
  • the fluid circulation chambers are even modeled by decomposition into elementary cells whose number is equal to the number of elementary cells of the reaction chamber.
  • the modeling of the heat exchanges also takes into account the direction of passage of the fluid in each elementary cell of the reaction chamber relative to the direction of passage (co-current, countercurrent or cross-flow) or utility fluids in the corresponding elementary cells of the lateral chambers located on either side of this elementary cell of the reaction chamber. It also includes heat balance equations relating to the reaction fluid, the plates delimiting the lateral heat exchange chambers, and utilizing fluids circulating in the lateral chambers, and material balance equations of the reaction medium and of the utility fluids in the reaction chambers. elementary cells of the different chambers. Each reaction chamber is connected to inlet zones and outlet zones which are associated with an elementary cell of the reaction chamber and which are for example connecting elements between the reagent feed means and a pressure point. supplying a reaction chamber, and connecting elements between an outlet point of a reaction chamber and means for collecting the reaction product. In the case where two reaction chambers are interconnected in series, the output element of the first chamber coincides with the input element of the second chamber.
  • the dynamic model of the reactor advantageously comprises a model of the inlet and outlet zones of the reactor, the inlet and outlet zones of each reaction chamber and the transition zones between the different blocks of the reactor.
  • the dynamic model also offers the possibility of a feed point per elementary cell of the reaction chamber.
  • the integration software tool applied to the system of equations of the dynamic model of the reactor, makes it possible to solve these equations and to simulate the behavior of the reactor called "state of the reactor" for the chemical reaction studied.
  • the determination and the optimization of the dimensional parameters of the reactor are obtained, for given objectives and constraints, by adjustment of variable characteristics of the reactor comprising the number of feed points in the reaction chamber, the directions of passage of the operating fluids. relative to the reaction medium, the nature and distribution of the utility fluids, the total volume of the reaction chamber and / or the volumes of the side chambers with utility fluid circulation, and the determination and optimization of the functional parameters of the reactor are obtained, for given objectives and constraints, by adjustment of variable reaction data or operating parameters comprising the temperatures, pressures, compositions and / or flow rates of the reaction medium and the utility fluid (s).
  • the method consists, for example, in adjusting the feed rates of one or more reagents in different elementary cells in order to obtain a better yield of the reaction, and in defining limit values of the operating parameters which satisfy safety constraints and / or environmental.
  • the dynamic model of the reactor comprises a dynamic model of the reaction chamber formed from equations of material balance, enthalpy balance, pressure balance and volume constraints. . These equations make it possible to define the state and evolutions of the fluid medium in each elementary cell.
  • the characteristic quantities temperature, pressure, flow, composition, etc.
  • the physical properties density, viscosity, etc.
  • the physicochemical phenomena reaction, transfer of matter, etc.
  • reaction medium in two-phase medium (liquid-liquid or liquid-gas), where the reaction medium is a mixture of a continuous phase and a dispersed phase, and in the case where the reaction medium is a mixture of two miscible liquids, the state and the evolution of the reaction medium are defined in each elementary cell from equations for estimating the density, the specific heat, the thermal conductivity and the viscosity of the medium. reaction, and material balance equations, enthalpy balance, pressure balance and volume constraints.
  • the characteristic quantities temperature, pressure, flow, composition, etc.
  • the physical properties density, viscosity, etc.
  • the physicochemical phenomena chemical reactions (reaction, transfer of material, etc.) are homogeneous in each cell, the temperatures and the pressures of the phases are identical, the reaction medium is assimilated to a pseudo-homogeneous medium and the physical properties of the two phases are determined by laws of association to connect the physical properties of the medium to the properties of each of the phases and their proportions, the distribution of the phases being assumed homogeneous in each elementary cell.
  • the Sauter diameter makes it possible to determine the size of the droplets of the dispersed phase.
  • the evolution of the size of the droplets of the dispersed phase can further be followed in the different elementary cells of the reaction chamber.
  • Other models can complement and specify the dynamic model of the reactor, such as:
  • the process may also include an estimation of the physical properties of the reactor components and the reaction and utility fluids.
  • the method according to the invention has the advantage of predicting the behavior of a plate reactor for given operating conditions, and also makes it possible to define the optimal sizing or configuration and the optimal operation or operating conditions of the reactor for a reaction. or a given set of reactions. This process also makes it possible to describe the evolutions of this behavior over time, the representation of this evolution offering many possibilities in terms of optimal control of the reactor, safety studies and management of unproductive transitional periods.
  • FIG. 1 is a diagrammatic exploded perspective view of an open type plate reactor used in the process according to the invention
  • FIG. 2 is a schematic view of the reaction chamber of the reactor of FIG. 1, broken down into elementary cells;
  • FIG. 3 is a schematic view of an elementary cell of the reaction chamber of FIG. 2;
  • FIG. 4 is a schematic view of an elementary reaction cell associated with two elementary cells for utilitarian fluid circulation
  • FIGS. 5a, 5b and 5c are diagrams illustrating the main steps of the optimization method according to the invention.
  • FIGS. 6 and 7 are graphs showing the influence of the concentration and the flow rate of the reagents, and the flow rate of utilitarian fluid on the temperature and the yield of the oxidation reaction of sodium thiosulfate;
  • FIGS. 8 to 10 are graphs showing the influence of the number and the arrangement of the feed points of the reactants on the temperature of the hydrolysis reaction of acetic anhydride
  • FIGS. 11 and 12 are graphs showing the influence of the feed rate of the reactants at several feed points on the temperature and the yield of the neopentyl glycol production reaction.
  • Figure 1 shows schematically an open or continuous reactor comprising a stack of plates 20, 21, 22, 23 which delimit between them parallel chambers.
  • transition plates delimit between them a reaction chamber 24 comprising a point 25a for supplying a reagent and one or more points 25b for supplying another reagent or reagents, and a point 26 for leaving the reaction product or products.
  • the plates 20 and 23 are adiabatic plates and delimit, with the plates 21, 22, a left lateral chamber 27 and a right lateral chamber 28 for the circulation of heat transfer fluid, each lateral chamber 27, 28 comprising an inlet 29 and a utility fluid outlet 30.
  • the utility fluid flowing in the left side chamber 27 may be the same or different from that flowing in the right side chamber 28, the utility fluid commonly used being water or ethylene glycol.
  • the stack of plates 20, 21, 22, 23 forms a reaction block that the plates 20 and 23 can thermally isolate from the external environment or from one or more other reaction blocks connected in series with this one. this.
  • the chambers 24, 27, 28 are generally equipped with means determining the pathways for the passage of fluids within these chambers.
  • the reaction chamber 24 comprises, for example, baffle inserts which impose on the reaction medium a zigzag circulation inside the reaction chamber 24, as shown in FIGS. 1 and 2.
  • the chambers 24, 27, 28 are fed continuously by appropriate means not shown. Under optimum operating conditions, the overall heat transfer coefficient of a plate reactor can vary between 2000 and 5000 Wm -2 K -1 .
  • the control of the temperature in the reaction chamber 24 is carried out for example by a appropriate selection of the temperature and / or flow rate of the utility fluid or fluids circulating on either side of the reaction chamber 24, the composition, the temperature and / or the flow rate of the reaction medium in the reaction chamber , and / or by the meaning utilizing fluid passage in each side chamber relative to that of the reaction medium in the reaction chamber.
  • a commercial fluid flows cocurrently if it flows in the two lateral chambers 27, 28 in the same direction as the reaction medium in the reaction chamber 24, against the current if it flows in the opposite direction, and cross flow if it flows in a side chamber in the same direction as the reaction medium and in the opposite direction in the other side chamber, as shown in Figure 1.
  • the reaction medium and the utility fluids intersect in because of the zigzag horizontal circulation of the reaction medium and the vertical circulation of utility fluid, but any other circulation configuration can be adopted such as a parallel configuration.
  • An essential characteristic of the process according to the invention consists in establishing a dynamic model of this reactor for a given reaction.
  • the modeling includes a decomposition of the reaction chamber into perfectly stirred elementary cells by means of experimental residence time distribution studies allowing flow analysis with the detection of dead volumes and preferential passages, and characterization of the mixture in the reaction chamber 24.
  • this number depends in particular on the feed rate of the reaction chamber 24 and the viscosity of the reaction medium.
  • the number of elementary cells in the case of a plate reactor comprising three reaction blocks is for example 70 for a reaction medium consisting solely of water supplied at a flow rate of 25 Lh -1 , 100 for water fed at 50 Lh -1 , 130 for water fed at 80 Lh -1 and 25 for a reaction medium composed of ethylene glycol fed at a flow rate of 50 Lh -1 .
  • the elementary cells of the reaction chamber are 12 in number. These 12 cells are numbered from 1 to 12 according to the horizontal zig-zag circulation of the reaction medium in the reaction chamber 24. from the supply point (s) 25a, 25b to the exit point 26 of this chamber 24. This circulation is imposed by the abovementioned inserts present in the reaction chamber 24.
  • the 12 cells form 4 rows of 3 cells.
  • a reagent is introduced into the reaction chamber 24 from the feed point 25a in the cell 1.
  • the reaction medium passes horizontally from the cell 1 to the cell 2, then to the 3, it leaves the cell 3 and enters the cell 4 located vertically beneath it, then passes horizontally in the cells 5 and 6.
  • the reaction medium leaves the cell 6 and enters the cell 7 located vertically beneath it, then passes horizontally in the cells 8 and 9.
  • the reaction medium exits the cell 9 and enters the cell 10 located vertically beneath it, then passes horizontally in the cells 11 and 12 to exit the reaction chamber 24 by the point 26.
  • Other feed points 25b of one or more other reagents may further be provided on the first cell or other elementary cells than the first.
  • reaction medium can flow in the reaction chamber 24 in parallel with the utility fluid flowing in the side chambers 27, 28.
  • the side chambers 27 and 28 are decomposed into elementary cells for utilitarian fluid circulation, so that each elementary reaction cell is located between two elementary cells of utility fluid circulation.
  • the left lateral chambers 27 and 28 respectively comprise series A, B, C, D, E, F, G, H, I, J, K, L and A ', B', C, D '. , E ', F', G ', H', I ', J', K ', L' of elementary cells for utilitarian fluid circulation.
  • the alphabetical order of the letters of the cells takes into account the increasing order of the numbers of the elementary cells of reaction, and not of the circulation of utilitarian fluid which is up and down in the left chamber 27 and from bottom to top in the right chamber 28.
  • the decomposition of the reaction chamber 24 into elementary cells makes it possible to determine, for the reaction medium circulating in each elementary reaction cell, the thermal influence of the commercial fluid circulating in the two corresponding lateral elementary cells, through the two transition plates 21, 22 delimiting the reaction chamber 24.
  • Heat is for example exchanged between the elementary cell E of the utility fluid circulation and the elementary reaction cell 5 through the plate 21 and between the elementary cell E utilitarian fluid circulation and the elementary reaction cell 5 through the plate 22.
  • the thermal characteristics of the cell 5 are determined by those of the cell located just before it in the direction of circulation of the reaction medium, that is to say the cell 4. From the in the same way, the thermal characteristics of the cells E and E 'are determined respectively by those of the cells B and H', located before E and E 'in the direction of circulation of the utilitarian fluid.
  • the number of equations of the mathematical model is also a function of the number of utility zones in the utility fluid circulation chambers.
  • Utility area means a utility fluid circulation chamber portion that includes its own utility fluid supply and outlet points.
  • a utility fluid circulation chamber that includes three utility zones thus has three feed points and three utility fluid exit points.
  • n [i n thisUlblock + 3 ) X n block + 1 I * ( 16 + “corn) + 3 X " * "» / "** X n ⁇ bblock '
  • n [(" this muock + 3 ) xn b iock + l] x (21 + 2 xn com ) + 3 xn zonelblock xn bhck ;
  • n con s is the number of constituents of the chemical reaction
  • nbiock the number of reaction blocks
  • n zon e / biock the number of utility zones per reaction block
  • n ce ⁇ / biock the number of elementary cells by reaction block.
  • the characteristic equations are for example: for each elementary cell of the reaction chamber: a global material balance, partial material balance, an enthalpy balance, a pressure balance, a volume model, a molar enthalpy model, a volume constraint, a molar volume model, a heat balance of the thermal environment, thermal balances of transition plates, thermal balances of fluids utilities, thermal balances of adiabatic plates and material balances of utilitarian fluids;
  • the dynamic model of the reactor according to the invention is established on the basis of a dynamic model of the reaction chamber and of a dynamic model of heat exchanges, each dynamic model comprising equations related to an elementary reaction cell 5, schematically represented in Figure 3 by a rectangle.
  • the arrows 32 and 33 respectively represent the direction of circulation of the reaction medium in the elementary cell 5 of the reaction chamber 24, and the utility fluid in the lateral lateral cells E and E 'of the lateral chambers 27, 28, the direction of passage. fluid utilitarian in the two chambers 27, 28 being in this example co-current with respect to the circulation of the reaction medium in the reaction chamber 24, contrary to Figure 1.
  • the indices p and u respectively define the reaction medium and the utility fluid.
  • F p represents the molar flow rate of the reaction medium in an elementary cell in mol.s -1 , T p and T u respectively the temperature of the reaction medium and of the utility fluid in their elementary cell at 0 K, H p the molar enthalpy of the reaction medium in an elementary cell in J.
  • the dynamic model of the reaction chamber comprises for example the following equations:
  • the exponent cell denotes the physical cell, the exponent ⁇ the order of reaction, the index j the number of the reaction, the index i the constituent considered, t the time in s, T the temperature in 0 K, P the pressure, and ⁇ P the pressure drop in Pa which is defined by the following equation:
  • the dynamic model of the utility fluid comprises for example the following equations:
  • ncr denotes the number of cells per row of cells in the utility fluid circulation chamber (in Fig. 1, ncr is equal to three for chambers 27 and 28), and denotes the last cell in the utility fluid circulation chamber ( in Fig. 1, in the chamber 27, denotes L, nc-1 denotes K and nc-2 denotes J, and in the chamber 28, denotes A ', nc-1 denotes B' and nc-2 denotes C).
  • the dynamic model also comprises a model of behavior of the reagent (s) in the connecting element 34 (FIG. 2) between the reagent supplying means and the feed point 25a on the cell 1, and a model of behavior of the reagent reaction medium in the connecting element 35 between the outlet point 26 on the cell 12 and the reaction product collection means or a feed or inlet point on a cell of a second reaction chamber (not shown) (second reaction block).
  • the connecting elements 34, 35 in the dynamic modeling of the plate reactor is important because the cumulative volume of these elements is not negligible compared to the volume of the reaction chamber 24, and may have an influence on the yield of the chemical reaction.
  • the heat exchange is limited with the utilitarian fluids.
  • the model of the heat exchanges of the reaction medium in the elementary reaction cells comprises characteristic equations related to an elementary reaction cell 5 schematically represented in FIG. 4 by a rectangle which is delimited by the two plates 21 and 22, and two E cells.
  • E 'lateral utilitarian circulation circumscribed respectively by the plates 20, 21 and 22, 23 in Figure 1.
  • the arrows 32 and 33 respectively represent the direction of circulation of the reaction medium in the elementary cell 5 of the reaction chamber 24 and of utilitarian fluid in the elementary lateral cells E 1 E 'of the lateral chambers 27, 28, the direction of passage of utilitarian fluid in the two chambers 27, 28 being co-current with respect to the circulation of the reaction medium in the chamber of reaction 24, unlike FIG.
  • the utilitarian fluid temperature in a lateral elementary cell depends on that of the previous lateral elementary cell in the direction of utilitarian fluid flow and also the presence of utility areas in the reaction block. It is thus possible to represent all the thermal configurations that can be envisaged for the different directions of passage of utilitarian fluid with respect to the reaction medium.
  • the utilitarian fluids in the left chamber 27 and in the right chamber 28, respectively, are defined by the indices u L and u R , the other designations and representations being the same as previously described.
  • the heat exchange model of the reaction medium includes for example the following equations:
  • F p and F u respectively represent the volume flow rate of the reaction medium in mol. s "1 and utilitarian fluid in m 3 .s “ 1 , e the thickness of the plates 21, 22 delimiting the reaction chamber 24 in m, ⁇ the thermal conductivity in J. s "1 " m “1 . K “1 , h the heat transfer coefficient in J. s ' 1 .m “ 2 . K “1 , T the temperature in K, Cp the specific heat of the medium in J. kg -1 . K “1 , p the density of the medium in kg m '3 , ⁇ q the heat generated by the reactions in J. m ' 3 ' s 1 , ⁇ Hr molar reaction heat in J. mol " 1 , A heat exchange area in m 2 , and pu L and PUR respectively indicate the plate 21 between the reaction chamber 24 and the left chamber 27 and the plate 22 between the reaction chamber 24 and the right chamber 28.
  • the aforementioned model of heat exchange is completed to include the thermal inertia of the plate reactor.
  • the thermal inertia model makes it possible to take into account the structure and the composition of the various elements of the plate reactor. It comprises, for example, the following equations, which differ from the preceding equations in that they take into account the thermal environment of the different fluids: energy balance equation of the environment of the reaction medium in an elementary reaction cell (FIG. s "1 ):
  • thermodynamic equation of the utilitarian fluid environment in the lateral cell of the left chamber 27 in J. s- 1 , the same equation is used for the lateral cell of the right chamber 28:
  • the heat balance equation of the environment of the reaction medium in the elementary reaction cell which depends on the temperature of the thermal environment of the reaction medium, made it possible to take into account the thermal inertia due to the structure of the reaction medium. bedroom reaction 24 and inserts (eg polyetheretherketone) present in this chamber. For a given elementary cell, the heat transfer area and the mass of this thermal environment are calculated from the geometric characteristics of the inserts and the structure of the reaction chamber 24.
  • thermal equilibrium equations of the environment of the utilitarian fluids in the elementary cells of the left lateral chambers 27 and right 28, which depend on the temperature of the thermal environment of the utilitarian fluids, made it possible to take into account the thermal inertia due to adiabatic plates 20 and 23 (eg stainless steel) in the dynamic model.
  • the characteristics of the thermal environment are calculated from the geometric characteristics of these adiabatic plates.
  • the dynamic model of the reactor comprises all the above-mentioned characteristic and balance equations (energy, thermal, mass, etc.) established for each phase, the continuous phase being liquid and the dispersed phase being liquid. or gas.
  • the dynamic models of the reaction medium in the reaction chamber 24, utilitarian fluids in the lateral chambers 27, 28, of the behavior of the reaction medium in the input and output elements 35 of the reaction chamber or chambers 24, charge losses, thermal exchanges of the reaction medium in the elementary cells 1 to 12 and the thermal inertia of the reactor are adapted for the establishment of the dynamic model of a plate-type open reactor for a reaction in two-phase medium given.
  • n represents the molar number of the mol phase
  • M the molar mass of the phase in kg. mol "1 ,
  • m the total mass of the diphasic medium in kg, i the index of the continuous or dispersed phase, ⁇ the thermal conductivity of the medium in J. m- 1 . s ' 1 . K '1 and ⁇ the viscosity of the medium in Pa s.
  • the accuracy of the dynamic model of the reactor for a reaction in a two-phase medium is based on the stability of the reaction medium and therefore the equilibrium between the two phases.
  • the equilibrium rules used in the mathematical model depend on the size of the droplets of the dispersed phase in the continuous phase.
  • the size of the droplets and their distribution are determined by means of a software tool for evaluating the population of the droplets.
  • the software tool for evaluating the population makes it possible to predict the evolution over time of an initial droplet size distribution and therefore of the Sauter diameter in the different elementary reaction cells of the chamber 24 of the reactor. with plates. This evolution depends on the rates of appearance and disappearance of the droplets, which are directly related to the phenomena of coalescence and rupture. The phenomena of coalescence and rupture are introduced into the program from correlations using the physical and transport properties of the two-phase medium.
  • Material transfer modeling is based on the assumption that the dispersed phase is in the form of spherical droplets with a mean diameter of Sauter diameter.
  • the transfer of matter between the two phases is modeled from Whitman's double film theory.
  • the mass transfer characteristics are calculated for each elementary reaction cell.
  • all mass transfer characteristics mass transfer coefficient, interfacial area, solubility
  • Sauter diameter the physical properties of a two-phase liquid-liquid medium and the Sauter diameter. Consequently, for a given liquid-liquid diphasic medium, an accurate representation of the behavior of the plate reactor requires only a proper estimate of the Sauter diameter for each elementary cell of the reaction chamber.
  • FIG. 5a is a diagram illustrating the essential steps of the optimization method according to the invention, in which the dynamic model 36 of the plate reactor comprises a plurality of characteristic equations and reactor balance equations such as previously described.
  • An integration software tool 37 is used to simulate the behavior of the plate reactor from the dynamic model of the reactor for a given reaction.
  • This tool makes it possible to solve the mixed systems of algebraic and differential equations of the dynamic model, by coupling with a database which includes the hydrodynamic description of the reactor 38, the physical properties 39 of the reaction medium, utility fluid (s) and reactor component materials (inserts, plates, etc.), the kinetic and heat reaction models of the chemical reaction 40, the heat transfers 41 between the different fluids, and the transfers of liquid-liquid or liquid-gas material for the reactions in two-phase medium.
  • parameters or inputs that occur in the different models are needed by the tool for the integration of the dynamic model. These parameters include the operating conditions of the reaction such as the number of reactions, the number of reagents, the nature of the reagents, their flow rate, their temperature, their pressure, their duration of feed, the characteristics of the reactions, the stoichiometric coefficients.
  • a description of a basic configuration 44 of the reactor comprising, for example, the number of elementary cells per block the number of elementary cells per row, the number of reaction blocks, the number of utility zones per reaction block, the dimensions of the plate reactor, the temperature of the external environment, the heat exchange areas, the number of plates, their thickness, their density, their heat capacity, their thermal conductivity, the dimensions of inserts, their density, their heat capacity, and the nature of utility fluids, and other integration parameters such as tolerances on absolute errors or the final time of integration.
  • the skilled person knows this kind of tool and database.
  • the database is for example Bipphy ® .
  • the integration software tool is for example DISCo, capable of solving dynamic models with more than 5000 characteristic equations and balance sheet with great speed and accuracy.
  • DlSCo Do Integrate by Software Components
  • DlSCo Do Integrate by Software Components
  • This tool has many advantages such as a global treatment of the system, without discrimination between equations or variables, the possibility of integrating equations and therefore additional variables into the basic system, thus offering great flexibility to the mathematical model, procedures automatic detection of events and calculation of consistent initial conditions, particularly useful in the search and management of events (startup procedure, dynamics) and the possibility of treatment of the system by hollow matrix, which leads to a reduction important computing time.
  • the integration software tool 37 makes it possible to simulate the behavior of the chemical reactor called "state of the reactor" 45 for a given chemical reaction. This simulation makes it possible, for example, to predict the yield of the chemical reaction, the evolution of the temperature during the reaction, etc.
  • the operator first seeks to define an optimal sizing or configuration of the reactor 46. For this, he determines or evaluates the objective (s) 47 to be achieved, such as the improvement of the productivity (yield reaction), the respect of the environment (quantity of weak rejects), or a reduced cost, the constraints 48 to respect which may concern safety for example with a maximum temperature of the reaction medium not to be exceeded, productivity or the profitability with maximum or minimum amounts of predetermined reagents, or the environment with a limitation of toxic or harmful or hardly recyclable or biodegradable discharges, and the variable characteristics of reactor structure.
  • the objective (s) 47 to be achieved such as the improvement of the productivity (yield reaction), the respect of the environment (quantity of weak rejects), or a reduced cost
  • the constraints 48 to respect which may concern safety for example with a maximum temperature of the reaction medium not to be exceeded, productivity or the profitability with maximum or minimum amounts of predetermined reagents, or the environment with a limitation of toxic or harmful or hardly recyclable or biodegradable discharges, and
  • variable characteristics 49 comprise, for example, the number of reaction blocks, the number of utility zones per block, the number and the position of the feed points, the dimensions of the reaction chambers, the co-current flow direction. , against current or cross flows of utilizing fluids relative to the circulation of the reaction medium in the reaction chamber, and the number, nature and geometry of the inserts in the reaction chamber.
  • the objective to be achieved may be a mathematical function comprising several quantitative criteria.
  • the operator can also define acceptable variations of these objectives and constraints, that is, allowable margins.
  • the method according to the invention makes it possible to optimize the optimal sizing or configuration 46 of the reactor from the state of reactor 45 previously defined, the operating conditions 43 mentioned above and of these objectives 47, constraints 48 and variables 49. optimization is obtained by a plurality of successive simulations 50 which make it possible to achieve as many objectives or objective functions as possible while respecting the constraints and even the authorized margins of the objectives and constraints, and by acting on all the variable structure characteristics plate reactor.
  • optimization of the design of the reactor is carried out by means of a second software tool 51 such as, for example, SQP ® (Successive Quadratic Programming) which makes it possible to define the parameters of the design or optimal configuration 46 of the reactor for a given reaction.
  • SQP ® Successessive Quadratic Programming
  • variable reaction data 53 which comprise, for example, the flow rate, the temperature, the pressure and the composition of the various reagents, the feed order of these reagents in the reaction chamber 24, and the flow and temperature of utility fluids.
  • the method according to the invention makes it possible to optimize the operation of the reactor by determining the optimal operating conditions of the reactor for a given reaction, based on the design or optimal configuration of the reactor previously defined, the aforementioned operating conditions and objectives, constraints and variables mentioned above.
  • the optimization is obtained by a plurality of successive simulations 50 which make it possible to achieve as many of the above objectives or objective functions as possible while respecting the constraints and even the authorized margins of the objectives and constraints, and by acting on all the data. reaction variables.
  • optimization of the operation of the reactor is also carried out by means of the aforementioned software tool 51.
  • the optimal operating conditions make it possible, for example, to carry out comparative studies in terms of reactor control, and safety and sensitivity of the chemical reaction. Examples of optimization of the behavior of a plate reactor for a given dimensioning and for a given reaction are illustrated in FIGS. 6 to 12.
  • Figures 6 and 7 are graphs for monitoring the oxidation reaction of sodium thiosulfate.
  • the axis 149 of the abscissa represents the volume of the reaction chamber in liters.
  • the axis 150 of the ordinates on the left represents the temperature in degrees Celsius and the axis 151 of the ordinates on the right represents the efficiency of the reaction on a scale of 1.
  • the curve 152 is the evolution of the yield of the oxidation reaction
  • curves 153, 154 and 155 are respectively the evolution of the temperature of the reaction medium, utility fluids and plates defining the reaction chamber.
  • the oxidation of sodium thiosulfate is carried out with hydrogen peroxide according to the reaction below, the reaction products being sodium trithionate, sodium sulfate and water.
  • the reaction is highly exothermic and carried out in a homogeneous liquid medium in a plate reactor comprising three reaction blocks.
  • the aforementioned experimental studies of residence time distribution make it possible to assimilate the plate reactor to a succession of 91 elementary cells, the reactor comprising three parallel reaction blocks formed of 27 elementary cells each, and separated from each other by a transition element formed of 10 elementary cells.
  • the number of equations in the dynamic model of the reactor for, this reaction is 1920.
  • the reaction medium is not very concentrated in reagents
  • the sodium thiosulfate is injected at 40 Lh -1 with a concentration of 0.62 mol.l -1
  • the hydrogen peroxide at 10 Lh -1 with a concentration of of 1.28 mol L -1
  • the utility fluids are water at 14 0 C fed at a flow rate of 1.8 m 3 . h -1 in the first reaction block and 0.9 m 3 h -1 in the other two reaction blocks.
  • the operator sets a single objective, improving the efficiency of the oxidation reaction, a variable reaction data which is the flow rates of the utility fluids, and a temperature constraint of the reaction medium which must remain below 100 ° C.
  • the reaction medium is more concentrated in reagent
  • the sodium thiosulfate is injected at 40 Lh -1 with a concentration of 0.75 mol L -1
  • the hydrogen peroxide at 10 L -1 -1 with a concentration of of 1. 59 mol L -1
  • the utility fluids are water at 14 ° C supplied at a rate of 2.3 m 3 ⁇ h -1 in the first reaction block and 1.2 m 3 .h "1 of the other two blocks reactions.
  • the process according to the invention makes it possible to optimize the plate reactor for the oxidation reaction of sodium thiosulfate by acting on variable reaction data which are the flow rates of the utility fluids.
  • the oxidation reaction of sodium thiosulfate is then carried out in a plate reactor built to operate with the optimized flow rates of the utility fluids, and measurements of physical parameters such as temperature and pressure are carried out during this implementation. to validate the operation of the reactor.
  • the diamond-shaped points 156 are experimental measurements of the reaction efficiency; the square-shaped points 157 and the triangle 158 are respectively experimental measurements of the temperature of the reaction medium and utility fluids.
  • the simulation also makes it possible to precisely locate the temperature peak 158 in FIG. 7 (at more than 90 ° C.), which does not appear in the experimental measurements because of the small number of sensors and their arrangement, the maximum temperature measured being less than 60 ° C.
  • FIGS. 8 to 10 are graphs for monitoring the hydrolysis reaction of acetic anhydride.
  • the axis 159 of the abscissa represents the number of the elementary reaction cell.
  • the plate reactor comprises 30 elementary cells separated into three reaction blocks, the passage of a reaction chamber from one block to the next being schematically represented by the vertical dashed lines 160.
  • the ordinate represents the temperature in degrees Celsius.
  • the curves 162, 163 and 164 are respectively the evolution of the temperature of the reaction medium, Utility fluids and transition plates delineating the reaction chamber.
  • the operator sets a goal of maximum productivity, a temperature constraint of the reaction medium which must remain below 80 ° C., and variable characteristics of the reactor structure which are the number and the arrangement of the feed points of the reagents for the reaction medium. optimization of the design of the reactor.
  • the acetic anhydride was entirely injected into the cell 1 of the reaction chamber at a temperature of 30 ° C. and a flow rate of 40 L. h -1 .
  • the water was also injected into the cell. 1 at a temperature of 30 ° C and a flow rate of 10 Lh "1 .
  • Utility fluids have a temperature of 15 0 C and a flow rate of 10800 Lh "1 because the reaction is very exothermic It is seen that the chosen values of the optimization variables do not allow to respect the temperature constraint.
  • the acetic anhydride was entirely injected into the cell 1 of the reaction chamber, at a temperature of 30 ° C. and a flow rate of 40 Lh -1,
  • the water was injected into the cell 1 (first cell of the first block), at a temperature of 30 ° C. and a flow rate of 5 Lh "1 , and in the cell 11 (first cell of the second reaction block), at a temperature of 30 ° C. and a flow rate of 5 Lh "1.
  • Utility fluids have a temperature of 15 ° C and a flow rate of 10800 Lh " 1 .
  • the water was therefore introduced into the reaction chamber twice and at two different feed points. This makes it possible to reduce a little the maximum temperature of the reaction medium which goes from 138 ° C. in FIG.
  • the acetic anhydride was entirely injected into the cell 1 of the reaction chamber, at a temperature of 30 ° C. and a flow rate of 40 Lh -1, The water was injected into the cell 1, a temperature of 30 ° C. and a flow rate of 3.33 Lh -1 , in the cell 11, at a temperature of 30 ° C. and a flow rate of 3.33 Lh -1 , and in the cell 21 (first cell of the third block), at a temperature of 30 ° C. and a flow rate of 3.33 Lh -1 .
  • Utility fluids have a temperature of 15 ° C and a flow rate of 10800 Lh "1 .
  • the water was here introduced into the reaction chamber in three times and at three different feed points. This makes it possible to reduce the maximum temperature of the reaction medium below the constraint of 80 ° C.
  • the optimization process has made it possible to identify an appropriate dimensioning meeting the objective of productivity and respecting the temperature constraint.
  • Figures 11 and 12 are graphs for monitoring the production reaction of neopentyl glycol.
  • the axis 165 of the abscissa represents the volume of the reaction chamber in liters.
  • the reaction chamber comprises 40 elementary cells separated into four blocks of elementary cells each having a volume of
  • the axis 166 of the ordinates on the left represents the temperature in degrees Celsius and the axis 167 of the ordinates on the right represents the yield of the reaction in percentage.
  • Curve 168 is the evolution of the yield of the reaction
  • curves 169, 170 and 171 are respectively the evolution of the temperature of the reaction medium, utilitarian fluids and plates delimiting the reaction chamber.
  • neopentyl glycol 2 (formaldehyde) + (2-methyl propanal) + (sodium hydroxide) -> (neopentyl glycol) + (sodium formate)
  • the dashed vertical lines 172 and 173 respectively represent the feed points of sodium hydroxide and 2-methyl propanal.
  • Sodium hydroxide is therefore introduced at a single time when the volume of the reaction chamber is close to 0.05L, which corresponds to the second cell of the reaction chamber (0.035L per elementary cell), and the 2-methylpropanal is introduced in three times, in the third cell, the fifteenth cell and the twenty-eighth cell.
  • the utility fluids have a temperature of 50 0 C and flow rate of 1 5 m 3 .h -1.
  • the reagents are injected at a temperature of 18 ° C, and the 2-methylpropanal has a concentration of 2 mol. kg -1 for a total flow rate of 9.14 Lh -1 (three injections at identical flow rates of approximately 3.05 Lh -1 ), the conversion rate achieved is 89.4% and the maximum temperature of the reaction medium at The course of the reaction is 66.6 ° C.
  • the operator sets a goal of increasing the yield of the production reaction of neopentyl glycol, variable data for optimization of the operation of the reaction which are the feed rates of 2-methylpropanal for each feed point, and a constraint to be respected which is a limit temperature of the reaction medium not to be exceeded, which is 65 ° C.
  • the optimization process has thus made it possible to achieve a goal of yielding a reaction and of meeting a temperature constraint by acting on variable reaction data which are the feed rates of a reagent at several points of the reaction chamber.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Organic Chemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Physical Or Chemical Processes And Apparatus (AREA)

Abstract

Procédé d'optimisation d'une réaction chimique dans un réacteur ouvert du type à plaques, qui consiste à : - établir un modèle dynamique du réacteur pour cette réaction, - appliquer un outil logiciel d'intégration à ce modèle dynamique pour résoudre ses équations, - déterminer et optimiser un ensemble de paramètres du réacteur à partir d'évaluations des objectifs à atteindre, des contraintes à respecter et des données variables du réacteur et de la réaction, - construire un réacteur conforme à l'ensemble optimisé de paramètres, - et faire des mesures de paramètres physiques lors de la mise en œuvre de la réaction dans Ie réacteur construit, pour valider son fonctionnement.

Description

Optimisation d'une réaction chimique dans un réacteur ouvert du type à plaques
La présente invention concerne un procédé d'optimisation d'une réaction chimique dans un réacteur ouvert du type à plaques constitué d'un empilement de plaques délimitant entre elles au moins un bloc comprenant une chambre de réaction et deux chambres latérales de circulation d'un fluide utilitaire d'échange de chaleur.
On cherche actuellement à optimiser le déroulement et à améliorer le rendement des réactions chimiques en général, qui sont réalisées le plus souvent dans des réacteurs du type « batch », c'est à dire des réacteurs à fonctionnement discontinu, comprenant des cuves dans lesquelles on dépose et on mélange des quantités déterminées de réactifs que l'on soumet ensuite à des conditions optimales d'environnement (pression, température, pH, etc.) pour favoriser les réactions.
Un inconvénient majeur de ces réacteurs est la difficulté d'évacuer la chaleur produite par une réaction, ce qui limite les possibilités de contrôle et d'optimisation de la réaction.
La présente invention a notamment pour but d'apporter une solution simple et efficace au contrôle et à l'optimisation des réactions chimiques grâce à l'optimisation d'un réacteur du type ouvert qui ne présente pas les inconvénients précités des réacteurs discontinus ou réacteurs du type
« batch ».
Elle propose à cet effet un procédé d'optimisation d'une réaction chimique dans un réacteur, caractérisé en ce qu'il consiste à :
- utiliser un réacteur ouvert constitué d'un empilement de plaques délimitant entre elles au moins un bloc comprenant une chambre de réaction formée entre deux chambres latérales d'échange thermique à circulation de fluide utilitaire d'échange de chaleur, des moyens d'alimentation de la chambre de réaction par un débit continu d'un ou de plusieurs réactifs et des moyens d'alimentation des deux chambres latérales par un débit continu de fluide utilitaire ;
- établir un modèle dynamique du réacteur pour une réaction donnée, à partir d'une modélisation de la chambre de réaction et des échanges thermiques entre cette chambre de réaction et les chambres latérales à circulation de fluide utilitaire, ce modèle dynamique comprenant des équations de bilans de matière et d'énergie et de contraintes ;
- appliquer un outil logiciel d'intégration à ce modèle dynamique pour résoudre les équations précitées ; - déterminer et optimiser un ensemble de paramètres dimensionnels et/ou fonctionnels du réacteur à partir d'évaluations des objectifs à atteindre, des contraintes à respecter, et des données variables du réacteur et/ou de la réaction ;
- construire un réacteur à plaques du type précité conforme audit ensemble optimisé de paramètres ; et
- faire des mesures de paramètres physiques tels que la température et la pression lors de la mise en œuvre de ladite réaction dans le réacteur construit, pour valider son fonctionnement.
L'utilisation d'un réacteur continu ou ouvert comprenant une superposition de plaques délimitant au moins une chambre de réaction entre deux chambres de circulation de fluide utilitaire, permet de résoudre les problèmes de dissipation de chaleur constatés dans les réacteurs discontinus ou fermés du type « batch », par une sélection appropriée des températures et débits du ou des fluides utilitaires, qui circulent de part et d'autre de la chambre de réaction, et également par la sélection des sens de passage de ces fluides (à co-courant, à contre-courant ou à courants croisés) par rapport au sens de passage du fluide réactionnel dans la chambre de réaction.
Cela permet aussi d'optimiser le dimensionnement ou la configuration optimale du réacteur ainsi que le fonctionnement ou les conditions opératoires optimales du réacteur, en fournissant des paramètres de contrôle de la réaction chimique, et de décrire avec précision la distribution de la température dans le réacteur.
La détermination précise de la loi de variation de la température dans ce réacteur, qui constitue la différence essentielle entre l'invention et les procédés connus d'optimisation des réactions dans des réacteurs discontinus, est rendue possible par la modélisation dynamique du réacteur pour la réaction étudiée, et par l'optimisation d'un ensemble de paramètres dimensionnels et/ou fonctionnels du réacteur en fonction des objectifs à atteindre et des contraintes à respecter. Un réacteur à plaques du type précité est ensuite construit conformément à l'ensemble de paramètres dimensionnels et/ou fonctionnels préalablement optimisé, la réaction chimique est mise en œuvre dans ce réacteur et des mesures expérimentales de paramètres physiques sont faites pour vérifier l'optimisation de la réaction chimique. Le procédé selon l'invention peut également consister, après la construction du réacteur, à :
- établir un nouveau modèle dynamique de réacteur à partir du réacteur construit ;
- appliquer l'outil logiciel d'intégration à ce nouveau modèle dynamique pour résoudre les équations précitées ;
- optimiser à nouveau ledit ensemble de paramètres dimensionnels et /ou fonctionnels du réacteur ; et
- modifier si nécessaire ledit réacteur construit, conformément au nouvel ensemble optimisé de paramètres. Le modèle dynamique du réacteur comprend des équations de bilan de matière et d'énergie et des équations de contraintes, qui sont caractéristiques du réacteur et de la réaction chimique.
Selon une autre caractéristique de l'invention, la modélisation de la chambre de réaction comprend une décomposition de cette chambre en cellules élémentaires successives contenant chacune un milieu fluide parfaitement agité. Les chambres de circulation de fluide utilitaire sont de même modélisées par décomposition en cellules élémentaires dont le nombre est égal au nombre de cellules élémentaires de la chambre de réaction.
La modélisation des échanges thermiques tient également compte du sens de passage du fluide dans chaque cellule élémentaire de la chambre de réaction par rapport au sens de passage (à co-courant, à contre-courant ou à courants croisés) du ou des fluides utilitaires dans les cellules élémentaires correspondantes des chambres latérales situées de part et d'autre de cette cellule élémentaire de la chambre de réaction. Elle comprend également des équations de bilan thermique relatives au fluide réactionnel, aux plaques délimitant les chambres latérales d'échange thermique, et aux fluides utilitaires circulant dans les chambres latérales, et des équations de bilan de matière du milieu réactionnel et des fluides utilitaires dans les cellules élémentaires des différentes chambres. Chaque chambre de réaction est reliée à des zones d'entrée et des zones de sortie qui sont associées à une cellule élémentaire de la chambre de réaction et qui sont par exemple des éléments de liaison entre les moyens d'alimentation en réactifs et un point d'alimentation d'une chambre de réaction, et des éléments de liaison entre un point de sortie d'une chambre de réaction et des moyens de collecte du produit de la réaction. Dans le cas où deux chambres de réaction sont reliées entre elles en série, l'élément de sortie de la première chambre est confondu avec l'élément d'entrée de la seconde chambre.
Le modèle dynamique du réacteur comprend avantageusement un modèle des zones d'entrée et de sortie du réacteur, des zones d'entrée et de sortie de chaque chambre de réaction et des zones de transition entre les différents blocs du réacteur.
Le modèle dynamique offre en outre la possibilité d'un point d'alimentation par cellule élémentaire de la chambre de réaction. L'outil logiciel d'intégration, appliqué au système d'équations du modèle dynamique du réacteur, permet de résoudre ces équations et de simuler le comportement du réacteur appelé « état du réacteur » pour la réaction chimique étudiée.
La détermination et l'optimisation des paramètres dimensionnels du réacteur sont obtenues, pour des objectifs et des contraintes donnés, par ajustement de caractéristiques variables du réacteur comprenant le nombre de points d'alimentation dans la chambre de réaction, les sens de passage des fluides utilitaires par rapport au milieu réactionnel, la nature et la distribution des fluides utilitaires, Ie volume total de la chambre de réaction et/ou les volumes des chambres latérales à circulation de fluide utilitaire, et la détermination et l'optimisation des paramètres fonctionnels du réacteur sont obtenues, pour des objectifs et des contraintes donnés, par ajustement de données variables de réaction ou de paramètres opératoires comprenant les températures, les pressions, les compositions et/ou les débits du milieu réactionnel et du ou des fluides utilitaires. Le procédé consiste par exemple à ajuster les débits d'alimentation d'un ou de plusieurs réactifs dans différentes cellules élémentaires afin d'obtenir un meilleur rendement de la réaction, et à définir des valeurs limites des paramètres opératoires satisfaisant à des contraintes de sécurité et/ou environnementales. Dans le cas d'une réaction chimique en milieu homogène, le modèle dynamique du réacteur comprend un modèle dynamique de la chambre de réaction formé à partir d'équations de bilan de matière, de bilan enthalpique, de bilan de pression et de contraintes de volume. Ces équations permettent de définir l'état et les évolutions du milieu fluide dans chaque cellule élémentaire.
Pour cela, on considère que dans les cellules élémentaires à contenu parfaitement agité de la chambre de réaction, les grandeurs caractéristiques (température, pression, débit, composition, etc.), les propriétés physiques (masse volumique, viscosité, etc.) et les phénomènes physico-chimiques (réaction, transfert de matière, etc.) sont homogènes dans chaque cellule et les propriétés physiques des fluides utilitaires sont homogènes et constantes pour une température donnée.
Dans le cas d'une réaction chimique en milieu diphasique (liquide- liquide ou liquide-gaz), où le milieu réactionnel est un mélange d'une phase continue et d'une phase dispersée, et dans le cas où le milieu réactionnel est un mélange de deux liquides miscibles, l'état et l'évolution du milieu réactionnel sont définis dans chaque cellule élémentaire à partir d'équations d'estimation de la masse volumique, de la chaleur spécifique, de la conductivité thermique et de la viscosité du milieu réactionnel, et d'équations de bilan de matière, de bilan enthalpique, de bilan de pression et de contraintes de volume.
On considère que dans les cellules élémentaires à contenu parfaitement agité de Ia chambre de réaction, les grandeurs caractéristiques (température, pression, débit, composition, etc.), les propriétés physiques (masse volumique, viscosité, etc.) et les phénomènes physico-chimiques (réaction, transfert de matière, etc.) sont homogènes dans chaque cellule, les températures et les pressions des phases sont identiques, on assimile le milieu réactionnel à un milieu pseudo-homogène et les propriétés physiques des deux phases sont déterminées par des lois d'association permettant de relier les propriétés physiques du milieu aux propriétés de chacune des phases ainsi qu'à leurs proportions, la répartition des phases étant supposée homogène dans chaque cellule élémentaire.
On peut en outre appliquer au mélange de phases des règles d'équilibre entre les deux phases en fonction de la dimension des gouttelettes de la phase dispersée, des propriétés thermodynamiques des phases en présence, et/ou des conditions d'écoulement.
Le diamètre de Sauter permet de déterminer la dimension des gouttelettes de la phase dispersée. L'évolution de la taille des gouttelettes de la phase dispersée peut en outre être suivie dans les différentes cellules élémentaires de la chambre de réaction. D'autres modélisations peuvent venir compléter et préciser le modèle dynamique du réacteur, telles que :
- une modélisation du comportement hydrodynamique du réacteur faite à partir d'études expérimentales de distribution de temps de séjour des fluides dans les différentes parties du réacteur, et permettant de définir le nombre de cellules élémentaires à contenu parfaitement agité ;
- une modélisation de la réaction à partir d'équations de vitesse de la réaction, de taux de production de constituants de la réaction et de génération de chaleur, dans chaque cellule élémentaire de la ou des chambres de réaction ;
- une modélisation du transfert de matière, dans une réaction en milieu diphasique, entre les deux phases du milieu réactionnel, à partir de propriétés physiques du milieu diphasique et de la taille des gouttelettes de la phase dispersée ; - une modélisation du transfert thermique à partir de définitions des coefficients de transfert thermique entre les plaques délimitant les chambres latérales et les fluides circulant dans la chambre de réaction et dans les chambres latérales, et d'estimations des coefficients de film du fluide réactionnel et du fluide utilitaire ; et - une modélisation des pertes de charge du fluide réactionnel dans le réacteur, faite à partir de mesures expérimentales.
Le procédé peut comprendre également une estimation des propriétés physiques des composants du réacteur et des fluides réactionnels et utilitaires. Le procédé selon l'invention présente l'avantage de prédire le comportement d'un réacteur à plaques pour des conditions opératoires données, et permet également de définir le dimensionnement ou la configuration optimale et le fonctionnement ou les conditions opératoires optimales du réacteur pour une réaction ou un ensemble de réactions donné. Ce procédé permet aussi de décrire les évolutions de ce comportement au cours du temps, la représentation de cette évolution offrant de nombreuses possibilités en termes de contrôle optimal du réacteur, d'études de sécurité et de gestion des périodes transitoires non productives.
D'autres avantages et caractéristiques de l'invention apparaîtront à la lecture de la description suivante faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels :
- la figure 1 est une vue schématique en perspective éclatée d'un réacteur ouvert du type à plaques utilisé dans le procédé selon l'invention ;
- la figure 2 est une vue schématique de la chambre de réaction du réacteur de la figure 1 , décomposée en cellules élémentaires ;
- la figure 3 est une vue schématique d'une cellule élémentaire de la chambre de réaction de la figure 2 ;
- la figure 4 est une vue schématique d'une cellule élémentaire de réaction associée à deux cellules élémentaires de circulation de fluide utilitaire ; - les figures 5a, 5b et 5c sont des diagrammes illustrant les principales étapes du procédé d'optimisation selon l'invention ;
- les figures 6 et 7 sont des graphes représentant l'influence de la concentration et du débit des réactifs, et du débit de fluide utilitaire sur la température et le rendement de la réaction d'oxydation du thiosulfate de sodium ;
- les figures 8 à 10 sont des graphes représentant l'influence du nombre et de la disposition des points d'alimentation des réactifs sur la température de la réaction d'hydrolyse de l'anhydride acétique ;
- les figures 11 et 12 sont des graphes représentant l'influence du débit d'alimentation des réactifs en plusieurs points d'alimentation sur la température et le rendement de la réaction de production du néopentyl glycol.
On se réfère d'abord à la figure 1 qui représente schématiquement un réacteur ouvert ou continu comprenant un empilement de plaques 20, 21 , 22, 23 qui délimitent entre elles des chambres parallèles. Les plaques
21 et 22, appelées plaques de transition, délimitent entre elles une chambre de réaction 24 comprenant un point 25a d'alimentation d'un réactif et un ou plusieurs points 25b d'alimentation d'un autre ou de plusieurs autres réactifs, et un point 26 de sortie du ou des produits de la réaction. Les plaques 20 et 23 sont des plaques adiabatiques et délimitent avec les plaques 21 , 22 une chambre latérale gauche 27 et une chambre latérale droite 28 de circulation de fluide utilitaire d'échange de chaleur, chaque chambre latérale 27, 28 comprenant une entrée 29 et une sortie 30 de fluide utilitaire.
Le fluide utilitaire qui circule dans la chambre latérale gauche 27 peut être identique ou différent de celui qui circule dans la chambre latérale droite 28, le fluide utilitaire couramment utilisé étant l'eau ou l'éthylène glycol.
L'empilement des plaques 20, 21 , 22, 23 forme un bloc de réaction que les plaques 20 et 23 permettent d'isoler thermiquement de l'environnement extérieur ou d'un ou de plusieurs autres blocs de réaction montés en série avec celui-ci.
Les chambres 24, 27, 28 sont généralement équipées de moyens déterminant les trajets de passage des fluides à l'intérieur de ces chambres. La chambre de réaction 24 comprend par exemple des inserts en chicane qui imposent au milieu réactionnel une circulation en zigzag à l'intérieur de la chambre de réaction 24, comme cela est représenté aux figures 1 et 2.
Les chambres 24, 27, 28 sont alimentées en continu par des moyens appropriés non représentés. Dans des conditions optimales de fonctionnement, le coefficient de transfert thermique global d'un réacteur à plaques peut varier entre 2000 et 5000 W.m'2.K"1. Le contrôle de la température dans la chambre de réaction 24 est réalisé par exemple par une sélection appropriée de la température et/ou du débit du ou des fluides utilitaires qui circulent de part et d'autre de la chambre de réaction 24, de la composition, de la température et/ou du débit du milieu réactionnel dans la chambre de réaction, et/ou par le sens de passage de fluide utilitaire dans chaque chambre latérale par rapport à celui du milieu réactionnel dans la chambre de réaction. Un fluide utilitaire circule à co-courant s'il circule dans les deux chambres latérales 27, 28 dans le même sens que le milieu réactionnel dans la chambre de réaction 24, à contre-courant s'il circule dans le sens inverse, et à courants croisés s'il circule dans une chambre latérale dans le même sens que le milieu réactionnel et dans le sens inverse dans l'autre chambre latérale, comme représenté en figure 1. Dans cet exemple, le milieu réactionnel et les fluides utilitaires se croisent en raison de la circulation horizontale en zigzag du milieu réactionnel et de la circulation verticale de fluide utilitaire, mais toute autre configuration de circulation peut être adoptée telle qu'une configuration en parallèle.
Une caractéristique essentielle du procédé selon l'invention consiste à établir un modèle dynamique de ce réacteur pour une réaction donnée. La modélisation comprend une décomposition de la chambre de réaction en cellules élémentaires parfaitement agitées, au moyen d'études expérimentales de distribution de temps de séjour permettant une analyse de l'écoulement avec la mise en évidence de volumes morts et de passages préférentiels, et une caractérisation du mélange dans la chambre de réaction 24.
L'analyse de la distribution de temps de séjour expérimentale permet de la caractériser en termes de temps de séjour moyen (t ) et de variance (σ2). A partir de ces caractéristiques, il est possible d'estimer le nombre de cellules élémentaires représentatif du degré de mélange et de l'hydrodynamique du réacteur réel. Ce nombre est calculé par exemple au moyen de la formule suivante :
Figure imgf000012_0001
En pratique, ce nombre dépend notamment du débit d'alimentation de la chambre de réaction 24 et de la viscosité du milieu réactionnel. Le nombre de cellules élémentaires dans le cas d'un réacteur à plaques comprenant trois blocs de réaction est par exemple de 70 pour un milieu réactionnel composé uniquement d'eau alimentée à un débit de 25 Lh"1, de 100 pour de l'eau alimentée à 50 Lh"1, de 130 pour de l'eau alimentée à 80 Lh"1 et de 25 pour un milieu réactionnel composé d'éthylène glycol alimenté à un débit de 50 Lh"1.
Dans l'exemple des figures 1 et 2, les cellules élémentaires de la chambre de réaction sont au nombre de 12. Ces 12 cellules sont numérotées de 1 à 12 suivant la circulation horizontale en zig-zag du milieu réactionnel dans la chambre de réaction 24, depuis le ou les points d'alimentation 25a, 25b jusqu'au point de sortie 26 de cette chambre 24. Cette circulation est imposée par les inserts précités présents dans la chambre de réaction 24.
Les 12 cellules forment 4 rangées de 3 cellules. Un réactif est introduit dans la chambre de réaction 24 depuis le point d'alimentation 25a dans la cellule 1. Le milieu réactionnel passe horizontalement de la cellule 1 à la cellule 2, puis à la 3, il sort de la cellule 3 et pénètre dans la cellule 4 située verticalement sous elle, puis passe horizontalement dans les cellules 5 et 6. Le milieu réactionnel sort de la cellule 6 et pénètre dans la cellule 7 située verticalement sous elle, puis passe horizontalement dans les cellules 8 et 9. Enfin, le milieu réactionnel sort de la cellule 9 et pénètre dans la cellule 10 située verticalement sous elle, puis passe horizontalement dans les cellules 11 et 12 pour sortir de la chambre de réaction 24 par le point 26. D'autres points d'alimentation 25b d'un ou de plusieurs autres réactifs peuvent en outre être prévus sur Ia première cellule ou d'autres cellules élémentaires que la première.
En variante, le milieu réactionnel peut circuler dans la chambre de réaction 24 en parallèle par rapport au fluide utilitaire circulant dans les chambres latérales 27, 28.
De même, les chambres latérales 27 et 28 sont décomposées en cellules élémentaires de circulation de fluide utilitaire, afin que chaque cellule élémentaire de réaction soit située entre deux cellules élémentaires de circulation de fluide utilitaire.
En figure 1 , les chambres latérales gauche 27 et droite 28 comprennent respectivement des séries A, B, C, D, E, F, G, H, I, J, K, L et A', B', C, D', E', F', G', H', I', J', K', L' de cellules élémentaires de circulation de fluide utilitaire. L'ordre alphabétique des lettres des cellules tient compte de l'ordre croissant des chiffres des cellules élémentaires de réaction, et non pas de la circulation de fluide utilitaire qui est de haut en bas dans la chambre de gauche 27 et de bas en haut dans la chambre de droite 28. La décomposition de la chambre de réaction 24 en cellules élémentaires permet de déterminer pour le milieu réactionnel circulant dans chaque cellule élémentaire de réaction, l'influence thermique du fluide utilitaire circulant dans les deux cellules élémentaires latérales correspondantes, à travers les deux plaques 21 , 22 de transition délimitant la chambre de réaction 24. De la chaleur est par exemple échangée entre la cellule élémentaire E de circulation de fluide utilitaire et la cellule élémentaire de réaction 5 à travers la plaque 21 et entre la cellule élémentaire E' de circulation de fluide utilitaire et la cellule élémentaire de réaction 5 à travers la plaque 22. Dans l'exemple de réalisation représenté en figure 1 , les caractéristiques thermiques de la cellule 5 sont déterminées par celles de la cellule située juste avant elle dans le sens de circulation du milieu réactionnel, c'est-à-dire la cellule 4. De la même façon, les caractéristiques thermiques des cellules E et E' sont déterminées respectivement par celles des cellules B et H', situées avant E et E' dans le sens de circulation du fluide utilitaire.
Les cellules élémentaires de réaction sont considérées comme des réacteurs agités continus, ce qui permet de généraliser les équations caractéristiques et de bilan du réacteur et de la réaction chimique du modèle mathématique à un nombre d'équations qui est fonction du nombre de cellules élémentaires de réaction, et qui a pour expression : n = [nombre de cellules] * (nombre d'équations par cellule)
Le nombre d'équations du modèle mathématique est également fonction du nombre de zones d'utilité dans les chambres de circulation de fluide utilitaire. On entend par zone d'utilité, une partie de chambre de circulation de fluide utilitaire qui comprend ses propres points d'alimentation et de sortie de fluide utilitaire. Une chambre de circulation de fluide utilitaire qui comprend trois zones d'utilité présente donc trois points d'alimentation et trois points de sortie de fluide utilitaire. Le nombre d'équations du modèle mathématique a alors pour expression : n = [nombre de cellules] * (nombre d'équations par cellule) +
(nombre d'équations par zone) * [nombre de zones d'utilité] En fonction du modèle considéré, milieu homogène ou milieu diphasique, ce nombre d'équations n'est pas le même, le nombre total d'équations du modèle étant défini à partir des équations suivantes : - Modèle homogène :
n = [inceUlblock + 3) X nblock + 1I* (16 + «corn ) + 3 X «*»»/«** X n ιbblock '
- Modèle diphasique
n = [("cemuock + 3) x n biock + l]x (21 + 2 x ncom ) + 3 x nzonelblock x nbhck ;
dans lesquels ncons est le nombre de constituants de la réaction chimique, nbiock le nombre de blocs de réaction, nzone/biock le nombre de zones d'utilité par bloc de réaction et nceιι/biock le nombre de cellules élémentaires par bloc de réaction.
Pour une réaction chimique en milieu homogène ou monophasique, les équations caractéristiques sont par exemple : - pour chaque cellule élémentaire de la chambre de réaction : un bilan de matière global, des bilans de matière partiels, un bilan enthalpique, un bilan de pression, un modèle de volume, un modèle d'enthalpie molaire, une contrainte de volume, un modèle de volume molaire, un bilan thermique de l'environnement thermique, des bilans thermiques des plaques de transition, des bilans thermiques des fluides utilitaires, des bilans thermiques des plaques adiabatiques et des bilans de matière des fluides utilitaires ;
- pour chaque zone d'utilité de circulation de fluide utilitaire : des contraintes de température des fluides utilitaires et une contrainte de débit de fluides utilitaires.
Pour une réaction chimique en milieu diphasique ou hétérogène, les équations caractéristiques sont par exemple :
- pour chaque cellule élémentaire de la chambre de réaction : des bilans de matière globaux de la phase continue et de la phase dispersée, des bilans de matière partiels de la phase continue et de la phase dispersée, un bilan enthalpique, un bilan de pression, des modèles de volume de la phase continue et de la phase dispersée, des modèles d'enthalpie molaire de la phase continue et de la phase dispersée, des modèles de volume molaire de la phase continue et de la phase dispersée, des contraintes de volume de la phase dispersée et de la phase continue, un bilan thermique de l'environnement thermique, des bilans thermiques des plaques de transition, des bilans thermiques des fluides utilitaires, des bilans thermiques des plaques adiabatiques et des bilans de matières des fluides utilitaires ; - pour chaque zone d'utilité de circulation de fluide utilitaire : des contraintes de température des fluides utilitaires et une contrainte de débit des fluides utilitaires.
Le modèle dynamique du réacteur selon l'invention est établi à partir d'un modèle dynamique de la chambre de réaction et d'un modèle dynamique des échanges thermiques, chaque modèle dynamique comprenant des équations rapportées à une cellule élémentaire de réaction 5, schématiquement représentée en figure 3 par un rectangle.
Les flèches 32 et 33 représentent respectivement le sens de circulation du milieu réactionnel dans la cellule élémentaire 5 de la chambre de réaction 24, et du fluide utilitaire dans les cellules latérales élémentaires E et E' des chambres latérales 27, 28, le sens de passage du fluide utilitaire dans les deux chambres 27, 28 étant dans cet exemple à co- courant par rapport à la circulation du milieu réactionnel dans la chambre de réaction 24, contrairement à la figure 1.
On définit par les indices p et u respectivement le milieu réactionnel et le fluide utilitaire. Fp représente le débit molaire du milieu réactionnel dans une cellule élémentaire en mol.s"1, Tp et Tu respectivement la température du milieu réactionnel et de fluide utilitaire dans leur cellule élémentaire en 0K, Hp l'enthalpie molaire du milieu réactionnel dans une cellule élémentaire en J. mol'1, Vp le volume d'une cellule élémentaire de réaction en m3, up le nombre de moles dans le milieu réactionnel dans une cellule élémentaire et x la fraction molaire d'un réactif ou d'un constituant (i) dans une cellule élémentaire, les exposants k-1 et k indiquant le numéro de la cellule élémentaire de réaction (1 à 12 pour l'exemple de réalisation de la chambre de réaction 24 des figures 1 et 2) ou la lettre de Ia cellule élémentaire de circulation de fluide utilitaire (A à L et A' à L' pour les chambres gauche 27 et droite 28 de la figure 1).
Si on considère que dans les cellules élémentaires à contenu parfaitement agité de la chambre de réaction 24, la réaction chimique est réalisée dans un milieu homogène monophasique liquide dans lequel on peut appliquer la loi d'Arrhenius, le modèle dynamique de la chambre de réaction comprend par exemple les équations suivantes :
- équation de bilan de masse ou de matière global du milieu réactionnel dans une cellule élémentaire (en mol.s"1) :
^• = Fp'* + Fp" - F; + άnp k x V;
Figure imgf000017_0001
- équation de bilan de masse ou de matière partiel des constituants (i) du milieu réactionnel dans une cellule élémentaire (en mol.s'1) :
" P>, J = pf*χf* -FX, + Δ<; X F; dt + 1<P XP,< - équation de bilan de pression du milieu réactionnel dans une cellule élémentaire (en Pa) : pk _ pk-\ _ \ pk . rp ~ P ^- p '
- équation de modèle de volume du milieu réactionnel dans une cellule élémentaire (en m3) :
F; - «7F; = 0 avec mV; = up* x Vp ml Jc ,
- équation d'enthalpie molaire du milieu réactionnel dans une cellule élémentaire (en J. mol"1) :
Hp k - mHp k = 0 ; - contrainte de volume ou équation de débit molaire du milieu réactionnel dans une cellule élémentaire (en m3 ou mol. s"1) : Fp = 0 pendant la phase de remplissage de toutes les cellules ycell élémentaires de réaction p ~uis u Pk n y ml r Jrc = 0 «• v Pk = Vce" une fois les
P cellules remplies ; dans lesquelles Ff est le débit d'alimentation en mol. s'1, Δn représente la vitesse de production en mol.m"3.s~1, a ie coefficient stoechiométrique, r la vitesse de réaction en mol. m'3. s"1, k° le facteur pré-exponentiel, Ea l'énergie d'activation en J. mol"1, R Ia constante des gaz parfaits, C la concentration d'un réactif en mol. m"3, xf la fraction molaire dans le débit d'alimentation, V le volume en m3, mV et mH le modèle de volume et le modèle d'enthalpie en m3 et en J. mol"1, respectivement, Vml le volume molaire en m3.mol"1, l'exposant cell désigne la cellule physique, l'exposant β l'ordre de réaction, l'indice j le numéro de la réaction, l'indice i le constituant considéré, t le temps en s, T la température en 0K, P la pression, et ΔP la perte de charge en Pa qui est définie par l'équation suivante :
Δ? = — p v dh 2 μ dans laquelle le facteur de frottement ou de Fanning f est donné par / = 5.0464 Re"05328 pour un écoulement laminaire (nombre de Reynolds Re < 2200) du milieu réactionnel dans la cellule élémentaire de réaction ou / = 2.17347 Re"042316 pour un écoulement turbulent (Re > 2200) ; L étant la longueur du parcours du milieu réactionnel dans la chambre de réaction 24 en m, p la densité du milieu réactionnel en kg. m"3, v la vitesse du milieu réactionnel en m.s"1 et dh le diamètre hydraulique équivalent d'une cellule élémentaire de réaction en m.
Si on considère que le réacteur à plaques est constitué d'une succession de réacteurs agités continus et que les propriétés physiques du fluide utilitaire sont constantes et homogènes pour une température donnée, le modèle dynamique de fluide utilitaire comprend par exemple les équations suivantes :
- équation de continuité de débit de fluide utilitaire (en m3.s"1) : F;" = F:'" ;
- équation de continuité de température de fluide utilitaire (en K) : ne rpout _ ι=nc-ncι +1 ncr dans lesquelles les exposants in et out désignent respectivement l'entrée et la sortie de la chambre de circulation de fluide utilitaire, ncr désigne le nombre de cellule par rangée de cellules dans la chambre de circulation de fluide utilitaire (en figure 1 , ncr est égal à trois pour les chambres 27 et 28), et ne désigne la dernière cellule dans la chambre de circulation de fluide utilitaire (en figure 1 , dans la chambre 27, ne désigne L, nc-1 désigne K et nc-2 désigne J, et dans la chambre 28, ne désigne A', nc-1 désigne B' et nc-2 désigne C).
Le modèle dynamique comprend également un modèle de comportement du ou des réactifs dans l'élément de liaison 34 (figure 2) entre les moyens d'alimentation des réactifs et le point d'alimentation 25a sur la cellule 1 , et un modèle de comportement du milieu réactionnel dans l'élément de liaison 35 entre le point de sortie 26 sur la cellule 12 et les moyens de collecte des produits de la réaction ou un point d'alimentation ou d'entrée sur une cellule d'une deuxième chambre de réaction non représentée (deuxième bloc de réaction). La prise en compte des éléments de liaison 34, 35 dans la modélisation dynamique du réacteur à plaques est importante car le volume cumulé de ces éléments n'est pas négligeable par rapport au volume de la chambre de réaction 24, et peut avoir une influence sur le rendement de la réaction chimique. Par ailleurs, dans ces éléments de liaison, les échanges thermiques sont limités avec les fluides utilitaires.
Plusieurs éléments du réacteur à plaques participent aux échanges thermiques, tels que les plaques de transition 21 , 22, les chambres de circulation de fluide utilitaire 27, 28, les plaques adiabatiques 20, 23 et l'environnement thermique du milieu réactionnel (présence d'inserts dans la chambre de réaction 24).
Le modèle des échanges thermiques du milieu réactionnel dans les cellules élémentaires de réaction comprend des équations caractéristiques rapportées à une cellule élémentaire de réaction 5 schématiquement représentée en figure 4 par un rectangle qui est délimitée par les deux plaques 21 et 22, et deux cellules E, E' latérales de circulation de fluide utilitaire délimitées respectivement par les plaques 20, 21 et 22, 23 en figure 1. Les flèches 32 et 33 représentent respectivement le sens de circulation du milieu réactionnel dans la cellule élémentaire 5 de la chambre de réaction 24 et de fluide utilitaire dans les cellules latérales élémentaires E1 E' des chambres latérales 27, 28, le sens de passage de fluide utilitaire dans les deux chambres 27, 28 étant à co-courant par rapport à la circulation du milieu réactionnel dans la chambre de réaction 24, contrairement à la figure 1.
La température de fluide utilitaire dans une cellule élémentaire latérale dépend de celle de la cellule élémentaire latérale précédente dans le sens d'écoulement de fluide utilitaire et également de la présence de zones d'utilité dans le bloc de réaction. On peut ainsi représenter toutes les configurations thermiques envisageables pour les différents sens de passage de fluide utilitaire par rapport au milieu réactionnel.
On définit par les indices uL et uR, les fluides utilitaires respectivement dans la chambre gauche 27 et dans la chambre droite 28, les autres désignations et représentations étant les mêmes que précédemment décrites.
Si on considère que les propriétés physiques des plaques de transition et adiabatiques et du fluide utilitaire sont constantes pour une cellule élémentaire à une température donnée et que la température de chaque plaque est calculée au milieu de la paroi, le modèle des échanges thermiques du milieu réactionnel comprend par exemple les équations suivantes :
- équation de bilan énergétique du milieu réactionnel dans une cellule élémentaire de réaction (en J. s"1) : d(uk x Hk)
= F^H? - FkHp k + Aqp k x Vk + hp*Λ Ap k (^ -Tk) + h^ Ap k (Tj - Tk ) dt 1 Vk n hp k flR -- - , A Àk n -- A Ace" v x p
Figure imgf000021_0001
Figure imgf000021_0002
et Aqp k = ∑ (AHrj χ rj k) ;
- équation de bilan thermique de fluide utilitaire dans la cellule latérale de la chambre gauche 27 (en J. s"1, la même équation est utilisée pour la cellule latérale de la chambre droite 28) :
ok Vk Cnk — = Fk ok Cnk (Tk~λ - Tk } + hk Ak (Tk - Tk ) '
- équation de bilan thermique de la plaque 21 disposée entre une cellule élémentaire de réaction et la cellule latérale de la chambre gauche 27 (en J. s"1, la même équation est utilisée pour la plaque 22 disposée entre une cellule élémentaire de réaction et la cellule latérale de la chambre droite 28) :
P nk puj V pkuL C^Pn pkuL β P". L - - hnpkuL,p A-Αpk - T l pkuL ~)) + + n hpkuL,uL Λ Ace" - T l pkuL ï)
Figure imgf000022_0001
Figure imgf000022_0002
'
Figure imgf000022_0003
dans lesquelles Fp et Fu représentent respectivement le débit volumique du milieu réactionnel en mol. s"1 et de fluide utilitaire en m3.s"1 , e l'épaisseur des plaques 21 , 22 délimitant la chambre de réaction 24 en m, λ la conductivité thermique en J. s"1. m"1. K"1, h le coefficient de transfert thermique en J. s'1. m"2. K"1, T la température en K, Cp la chaleur spécifique du milieu en J. kg"1. K"1, p la masse volumique du milieu en kg. m'3, Δq la chaleur générée par les réactions en J. m'3. s'1, ΔHr la chaleur molaire de réaction en J. mol"1, A l'aire d'échange thermique en m2, et les indices puL et PUR désignent respectivement la plaque 21 entre la chambre de réaction 24 et la chambre gauche 27 et la plaque 22 entre la chambre de réaction 24 et la chambre droite 28.
Le modèle précité des échanges thermiques est complété pour inclure l'inertie thermique du réacteur à plaques. Le modèle d'inertie thermique permet de tenir compte de la structure et de la composition des différents éléments du réacteur à plaques. Il comprend par exemple les équations suivantes, qui diffèrent des équations précédentes en ce qu'elles prennent en compte l'environnement thermique des différents fluides : - équation de bilan énergétique de l'environnement du milieu réactionnel dans une cellule élémentaire de réaction (en J. s"1) :
P P = _ F/>kHf'k +Fk'1Hk'1 -FkHk + Aak y. Vk +hk Ak(Tk -Tk) • dt ~ x " "" " p
+ hk Àk(Tk - Tk) + hk Âk(Tk - Tk~) - équation de bilan thermique de l'environnement de fluide utilitaire dans la cellule latérale de la chambre gauche 27 (en J. s"1, la même équation est utilisée pour la cellule latérale de la chambre droite 28) :
nk Vk r-nk (Tk~λ — Tk Λ
Pu, VuL LPuL KJ- U1 1U1 )
Figure imgf000023_0001
. j Ic Λ cell /ηπk ηrt k \ , j k jcell /ηnk mk \
+ KuL,uL A (TP"L ~ TuL ) + KL,KL A (4,α - TUL ) - équation de bilan thermique de l'environnement du milieu réactionnel dans une cellule élémentaire de réaction (en J. s"1) :
A1Y k
P nk te / V tek , C^Pn tke „ ^ "' - ~ n hpk,lep Λ Âpk ( \XT pk - 1TIek,, ) ï > •
- équation de bilan thermique de l'environnement de la plaque adiabatique 20 de la chambre gauche 27 (en J. s"1 , la même équation est utilisée pour la plaque adiabatique 23 de la chambre droite 28) : k yk le«L _ 1 k 4k (Tk —Tk Λ •
P te.L V*,L
Figure imgf000023_0002
—jp- ~ KL ,,e,L Ap U ,,L * ,e,L ) , dans lesquelles les coefficients de transfert thermique ont pour expression (en J. m'2. K'1. s"1):
, k 1 , k _ ____J____ l,k _ et
L . l i P-11P 1 hU - 1
• + βP«L • + hp λp,puL K λp,'ep K K dans lesquelles les indices tep et teu désignent respectivement l'environnement thermique du milieu réactionnel et du fluide utilitaire, et Hf l'enthalpie molaire du débit d'alimentation en J. mol"1.
Les échanges thermiques entre le milieu réactionnel et les fluides utilitaires ont été pris en compte séparément pour chaque plaque de transition 21 , 22, afin d'inclure l'inertie thermique due à ces plaques (par exemple en acier inoxydable) dans le modèle dynamique.
L'équation de bilan thermique de l'environnement du milieu réactionnel dans la cellule élémentaire de réaction, qui dépend de la température de l'environnement thermique du milieu réactionnel, a permis de prendre en compte l'inertie thermique due à la structure de la chambre de réaction 24 et aux inserts (par exemple en polyétheréthercétone) présents dans cette chambre. Pour une cellule élémentaire donnée, la surface de transfert thermique et la masse de cet environnement thermique sont calculées à partir des caractéristiques géométriques des inserts et de la structure de la chambre de réaction 24.
Les équations de bilan thermique de l'environnement des fluides utilitaires dans les cellules élémentaires des chambres latérales gauche 27 et droite 28, qui dépendent de la température de l'environnement thermique des fluides utilitaires, ont permis de prendre en compte l'inertie thermique due aux plaques adiabatiques 20 et 23 (par exemple en acier inoxydable) dans le modèle dynamique. Les caractéristiques de l'environnement thermique (surface de transfert thermique, masse) sont calculées à partir des caractéristiques géométriques de ces plaques adiabatiques.
Pour une réaction chimique réalisée en milieu diphasique, le modèle dynamique du réacteur comprend toutes les équations caractéristiques et de bilan (énergétique, thermique, de masse, etc.) précitées établies pour chaque phase, la phase continue étant liquide et la phase dispersée étant liquide ou gaz. Les modèles dynamiques du milieu réactionnel dans la chambre de réaction 24, des fluides utilitaires dans les chambres latérales 27, 28, de comportement du milieu réactionnel dans les éléments d'entrée 34 et de sortie 35 de la ou des chambres de réaction 24, des pertes de charge, des échanges thermique du milieu réactionnel dans les cellules élémentaires 1 à 12 et de l'inertie thermique du réacteur, sont adaptés pour l'établissement du modèle dynamique d'un réacteur ouvert du type à plaques pour une réaction en milieu diphasique donnée.
Pour cela, les températures et les pressions dans chacune des phases sont considérées comme identiques et on assimile le milieu réactionnel à un milieu pseudo-homogène dont les propriétés physiques du milieu diphasique correspondant sont déterminées par des lois d'association qui relient les propriétés physiques du milieu diphasique aux propriétés de chacune des phases ainsi qu'à leurs proportions respectives, la répartition des phases étant supposée homogène dans chaque cellule élémentaire. Des exemples de lois d'association sont présentés ci- dessous :
- estimation de la masse volumique du milieu diphasique ;
Figure imgf000025_0001
- estimation de la chaleur spécifique du milieu diphasique
Cp = ^n1M1Cp1 ;
- estimation de la conductivité thermique du milieu diphasique
Figure imgf000025_0002
- estimation de la viscosité du milieu diphasique
Figure imgf000025_0003
dans lesquelles n, représente le nombre de mole de la phase en mol, M, la masse molaire de la phase en kg. mol"1, m la masse totale du milieu diphasique en kg, i l'indice de la phase continue ou dispersée, λ la conductivité thermique du milieu en J. m"1. s'1. K'1 et μ la viscosité du milieu en Pa. s.
La précision du modèle dynamique du réacteur pour une réaction en milieu diphasique est basée sur la stabilité du milieu réactionnel et donc de l'équilibre entre les deux phases. Les règles d'équilibre utilisées dans le modèle mathématique sont fonction de la dimension des gouttelettes de la phase dispersée dans la phase continue. La dimension des gouttelettes et leur distribution sont déterminées au moyen d'un outil logiciel d'évaluation de la population des gouttelettes.
L'outil logiciel d'évaluation de la population permet de prédire l'évolution au cours du temps d'une distribution de tailles de gouttelettes initiale et donc par conséquent du diamètre de Sauter dans les différentes cellules élémentaires de réaction de la chambre 24 du réacteur à plaques. Cette évolution dépend des vitesses d'apparition et de disparition des gouttelettes, qui sont directement liés aux phénomènes de coalescence et de rupture. Les phénomènes de coalescence et de rupture sont introduits dans le programme à partir de corrélations faisant appel aux propriétés physiques et de transport du milieu diphasique.
La modélisation du transfert de matière repose sur l'hypothèse que la phase dispersée se présente sous la forme de gouttelettes sphériques dont le diamètre moyen est le diamètre de Sauter. Le transfert de matière entre les deux phases est modélisé à partir de la théorie du double film de Whitman.
Dans le modèle dynamique du procédé d'optimisation selon l'invention, les caractéristiques de transfert de masse sont calculées pour chaque cellule élémentaire de réaction. Grâce à des corrélations de transfert de masse, toutes les caractéristiques de transfert de masse (coefficient de transfert de masse, surface interfaciale, solubilité) sont reliées aux propriétés physiques d'un milieu diphasique liquide-liquide et au diamètre de Sauter. En conséquence, pour un milieu diphasique liquide- liquide donnée, une représentation précise du comportement du réacteur à plaques nécessite seulement une estimation correcte du diamètre de Sauter pour chaque cellule élémentaire de la chambre de réaction.
On se réfère désormais à la figure 5a qui est un diagramme illustrant les étapes essentielles du procédé d'optimisation selon l'invention, dans laquelle le modèle dynamique 36 du réacteur à plaques comprend une pluralité d'équations caractéristiques et de bilan du réacteur telles que décrites précédemment.
Un outil logiciel d'intégration 37 est utilisé pour simuler le comportement du réacteur à plaques à partir du modèle dynamique du réacteur pour une réaction donnée. Cet outil permet de résoudre les systèmes mixtes d'équations algébriques et différentielles du modèle dynamique, par couplage avec une base de données qui comprend la description hydrodynamique du réacteur 38, les propriétés physiques 39 du milieu réactionnel, du ou des fluides utilitaires et des matériaux des composants du réacteur (inserts, plaques,...), les modèles cinétiques et de chaleur de réaction de la réaction chimique 40, les transferts thermiques 41 entre les différents fluides, et les transferts de matière 42 liquide-liquide ou liquide-gaz pour les réactions en milieu diphasique.
De nombreux paramètres ou entrées qui interviennent dans les différents modèles sont nécessaires à l'outil pour l'intégration du modèle dynamique. Ces paramètres comprennent les conditions opératoires 43 de la réaction tels que le nombre de réactions, le nombre de réactifs, la nature des réactifs, leur débit, leur température, leur pression, leur durée d'alimentation, les caractéristiques des réactions, les coefficients stoechiométriques, les ordres réactionnels, le facteur préexponentiel des constantes de vitesse, les énergies d'activation, et les débits et les températures des fluides utilitaires, un descriptif d'une configuration de base 44 du réacteur comprenant par exemple le nombre de cellules élémentaires par bloc de réaction, le nombre de cellules élémentaires par rangées, le nombre de blocs de réaction, le nombre de zones d'utilité par bloc de réaction, les dimensions du réacteur à plaques, la température du milieu extérieur, les aires d'échange thermique, le nombre de plaques, leur épaisseur, leur masse volumique, leur capacité calorifique, leur conductivité thermique, les dimensions des inserts, leur masse volumique, leur capacité calorifique, et la nature des fluides utilitaires, et d'autres paramètres d'intégration tels que les tolérances sur les erreurs absolues ou le temps final de l'intégration. L'homme du métier connaît ce genre d'outil et de base de données.
La base de donnée est par exemple Bipphy®. L'outil logiciel d'intégration est par exemple DISCo, capable de résoudre des modèles dynamiques avec plus de 5000 équations caractéristiques et de bilan avec une grande rapidité et précision. DlSCo (Do Integrate by Software Components) est un outil d'intégration de systèmes d'équations algébro-différentielles basé sur la méthode de Gear de différentiation rétrograde selon un schéma de prédicteur-correcteur. Cet outil présente de nombreux avantages tels qu'un traitement global du système, sans discrimination entre équations ou variables, la possibilité d'intégrer au système de base des équations et donc des variables supplémentaires, offrant ainsi une grande flexibilité au modèle mathématique, des procédures automatiques de détection des événements et de calcul de conditions initiales cohérentes, particulièrement utiles dans la recherche et la gestion d'événements (procédure de démarrage, dynamiques) et la possibilité d'un traitement du système par matrice creuse, ce qui conduit à une réduction importante du temps de calcul.
L'outil logiciel d'intégration 37 permet de simuler le comportement du réacteur chimique appelé « état du réacteur » 45 pour une réaction chimique donnée. Cette simulation permet par exemple de prévoir le rendement de la réaction chimique, l'évolution de la température au cours de la réaction, etc.
En figure 5b, l'opérateur cherche d'abord à définir un dimensionnement ou une configuration optimale du réacteur 46. Pour cela, il détermine ou évalue le ou les objectif(s) 47 à atteindre tels que l'amélioration de la productivité (rendement de la réaction), le respect de l'environnement (quantité de rejets faibles), ou un coût réduit, les contraintes 48 à respecter qui peuvent concerner la sécurité par exemple avec une température maximale du milieu réactionnel à ne pas dépasser, la productivité ou la rentabilité avec des quantités maximales ou minimales de réactifs prédéterminées, ou l'environnement avec une limitation des rejets toxiques ou néfastes ou difficilement recyclables ou biodégradables, et les caractéristiques variables 49 de structure du réacteur .
Ces caractéristiques variables 49 comprennent par exemple le nombre de blocs de réaction, le nombre de zones d'utilité par bloc, le nombre et la position des points d'alimentation, les dimensions des chambres de réaction, le sens de passage à co-courant, contre courant ou courants croisés des fluides utilitaires par rapport à la circulation du milieu réactionnel dans la chambre de réaction, et le nombre, la nature et la géométrie des inserts dans la chambre de réaction.
En variante, l'objectif à atteindre peut être une fonction mathématique comprenant plusieurs critères quantitatifs. A titre d'exemple, la fonction objectif d'une ligne de production peut être définie par rapport à un critère de productivité et un critère environnemental, par l'équation suivante : f = qp . kp - qx . /cx dans laquelle qp est la quantité de produit recherché de la réaction, kp le coût unitaire de ce produit, qx la quantité de produit non valorisable de la réaction, et kx le coût de traitement ou d'élimination de ce produit.
L'opérateur peut également définir des variations acceptables de ces objectifs et de ces contraintes, c'est-à-dire des marges autorisées. Le procédé selon l'invention permet d'optimiser le dimensionnement ou la configuration optimale 46 du réacteur à partir de l'état du réacteur 45 précédemment défini, des conditions opératoires 43 précitées et de ces objectifs 47, contraintes 48 et variables 49. L'optimisation est obtenue par une pluralité de simulations 50 successives qui permettent d'atteindre le plus possible d'objectifs ou de fonctions objectifs en respectant les contraintes, voire les marges autorisées des objectifs et des contraintes, et en agissant sur toutes les caractéristiques variables de structure du réacteur à plaques.
L'optimisation du dimensionnement du réacteur est réalisée au moyen d'un second outil logiciel 51 tel par exemple que SQP® (Successive Quadratic Programming) qui permet de définir les paramètres du dimensionnement ou de la configuration optimale 46 du réacteur pour une réaction donnée.
En figure 5c, l'opérateur cherche ensuite à déterminer les conditions opératoires optimales 52 pour optimiser le fonctionnement du réacteur pour la réaction donnée. Le ou les objectif(s) 47 à atteindre et les contraintes 48 à respecter restent les mêmes que lors de l'optimisation du dimensionnement du réacteur. Seules les variables changent, celles-ci étant des données variables de réaction 53 qui comprennent par exemple le débit, la température, la pression et la composition des différents réactifs, l'ordre d'alimentation de ces réactifs dans la chambre de réaction 24, et le débit et la température des fluides utilitaires.
Le procédé selon l'invention permet d'optimiser le fonctionnement du réacteur en déterminant les conditions opératoires optimales du réacteur pour une réaction donnée, à partir du dimensionnement ou de la configuration optimale du réacteur précédemment défini, des conditions opératoires précitées et des objectifs, contraintes et variables précités. L'optimisation est obtenue par une pluralité de simulations 50 successives qui permettent d'atteindre le plus possible d'objectifs ou de fonctions objectifs précitées en respectant les contraintes, voire les marges autorisées des objectifs et des contraintes, et en agissant sur toutes les données variables de réaction.
L'optimisation du fonctionnement du réacteur est également réalisée au moyen de l'outil logiciel 51 précité. Les conditions opératoires optimales permettent par exemple de réaliser des études comparatives en termes de contrôle du réacteur, et de sécurité et de sensibilité de la réaction chimique. Des exemples d'optimisation du comportement d'un réacteur à plaques pour un dimensionnement donné et pour une réaction donnée sont illustrés en figures 6 à 12.
Les figures 6 et 7 sont des graphes du suivi de la réaction d'oxydation du thiosulfate de sodium. L'axe 149 des abscisses représente le volume de la chambre de réaction en litre. L'axe 150 des ordonnés à gauche représente la température en degré Celsius et l'axe 151 des ordonnés à droite représente le rendement de la réaction sur une échelle de l La courbe 152 est l'évolution du rendement de la réaction d'oxydation, les courbe 153, 154 et 155 sont respectivement l'évolution de la température du milieu réactionnel, des fluides utilitaires et des plaques délimitant la chambre de réaction.
L'oxydation du thiosulfate de sodium est réalisée par du peroxyde d'hydrogène selon la réaction ci-dessous, les produits de la réaction étant du trithionate de sodium, du sulfate de sodium et de l'eau.
2 Na2S2O3 + 4 H2O2 → Na2S3O6 + Na2SO4 + 4 H2O
La réaction est très exothermique et réalisée en milieu homogène liquide dans un réacteur à plaques comprenant trois blocs de réaction. Les études expérimentales précitées de distribution de temps de séjour permettent d'assimiler le réacteur à plaques à une succession de 91 cellules élémentaires, le réacteur comprenant trois blocs de réaction en parallèle formés de 27 cellules élémentaires chacun, et séparés les uns des autres par un élément de transition formé de 10 cellules élémentaires. Le nombre d'équations dans le modèle dynamique du réacteur pour, cette réaction est de 1920.
En figure 6, le milieu réactionnel est peu concentré en réactifs, le thiosulfate de sodium est injecté à 40 L.h"1 avec une concentration de 0,62 mol. L"1, le peroxyde d'hydrogène à 10 L.h"1 avec une concentration de 1 ,28 mol. L"1, et les fluides utilitaires sont de l'eau à 140C alimentée à un débit de 1 ,8 m3. h"1 dans le premier bloc de réaction et de 0,9 m3. h"1 dans les deux autres blocs de réactions.
L'opérateur fixe un unique objectif, l'amélioration du rendement de la réaction d'oxydation, une donnée variable de réaction qui est les débits des fluides utilitaires, et une contrainte de température du milieu réactionnel qui doit rester inférieure à 1000C.
En figure 7, le milieu réactionnel est plus concentré en réactif, le thiosulfate de sodium est injecté à 40 L.h"1 avec une concentration de 0,75 mol. L"1, le peroxyde d'hydrogène à 10 L.h"1 avec une concentration de 1 ,59 mol. L"1, et les fluides utilitaires sont de l'eau à 14°C alimentée à un débit de 2,3 m3.h"1 dans le premier bloc de réaction et de 1 ,2 m3.h"1 dans les deux autres blocs de réactions. Le procédé selon l'invention permet d'optimiser le réacteur à plaques pour la réaction d'oxydation du thiosulfate de sodium en agissant sur des données variables de réaction qui sont les débits des fluides utilitaires.
La réaction d'oxydation du thiosulfate de sodium est ensuite mise en œuvre dans un réacteur à plaques construit pour fonctionner avec les débits optimisés des fluides utilitaires, et des mesures de paramètres physiques tels que la température et la pression sont réalisées lors de cette mise en œuvre pour valider le fonctionnement du réacteur.
En figures 6 et 7, les points en forme de losange 156 sont des mesures expérimentales du rendement de la réaction, les points en forme de carré 157 et de triangle 158 sont respectivement des mesures expérimentales de la température du milieu réactionnel et des fluides utilitaires.
Les résultats expérimentaux sont proches des courbes correspondantes, ce qui permet de valider la simulation de cette réaction chimique dans le réacteur à plaques et donc la précision du modèle dynamique.
La simulation permet d'autre part de localiser précisément le pic de température 158 en figure 7 (à plus de 900C), qui n'apparaît pas dans les mesures expérimentales à cause du petit nombre de capteurs et leur disposition, la température maximale mesurée étant inférieure à 600C.
Les figures 8 à 10 sont des graphes du suivi de la réaction d'hydrolyse de l'anhydride acétique. L'axe 159 des abscisses représente le numéro de la cellule élémentaire de réaction. Dans cet exemple de réalisation, le réacteur à plaques comprend 30 cellules élémentaires séparées en trois blocs de réaction, le passage d'une chambre de réaction d'un bloc à la suivante étant schématiquement représenté par les lignes verticales en pointillés 160. L'axe 161 des ordonnées représente la température en degré Celsius. Les courbes 162, 163 et 164 sont respectivement l'évolution de la température du milieu réactionnel, des fluides utilitaires et des plaques de transition délimitant la chambre de réaction.
L'hydrolyse de l'anhydride acétique produit de l'acide acétique selon la réaction ci-dessous. CH3-CO-O-CO-CH3 + H2O -> 2 CH3-CO-OH
L'opérateur fixe un objectif de productivité maximale, une contrainte de température du milieu réactionnel qui doit rester inférieure à 800C, et des caractéristiques variables de structure du réacteur qui sont le nombre et la disposition des points d'alimentation des réactifs pour l'optimisation du dimensionnement du réacteur.
En figure 8, l'anhydride acétique a été entièrement injecté dans la cellule 1 de la chambre de réaction, à une température de 300C et un débit de 40 L. h"1. L'eau a aussi été injectée dans la cellule 1 à une température de 30°C et un débit de 10 L.h"1. Les fluides utilitaires ont une température de 150C et un débit de 10800 L.h"1 car la réaction est très exothermique. On voit que les valeurs choisies des variables d'optimisation ne permettent pas de respecter la contrainte de température.
En figure 9, l'anhydride acétique a été entièrement injecté dans la cellule 1 de la chambre de réaction, à une température de 300C et un débit de 40 L.h"1. L'eau a été injectée dans la cellule 1 (première cellule du premier bloc), à une température de 300C et un débit de 5 L.h"1, et dans la cellule 11 (première cellule du deuxième bloc de réaction), à une température de 300C et un débit de 5 L.h"1. Les fluides utilitaires ont une température de 15°C et un débit de 10800 L.h"1. L'eau a donc été introduite dans la chambre de réaction en deux fois et en deux points d'alimentation différents. Cela permet de diminuer un peu la température maximale du milieu réactionnel qui passe de 138°C en figure 8 à 1040C en figure 9. Le nombre et la disposition des points d'alimentation des réactifs du réacteur ne sont toutefois pas corrects puisque la contrainte de température n'est toujours pas respectée. En figure 10, l'anhydride acétique a été entièrement injecté dans la cellule 1 de la chambre de réaction, à une température de 300C et un débit de 40 Lh"1. L'eau a été injectée dans la cellule 1 , à une température de 3O0C et un débit de 3,33 Lh"1, dans la cellule 11 , à une température de 30°C et un débit de 3, 33 Lh"1, et dans la cellule 21 (première cellule du troisième bloc), à une température de 300C et un débit de 3,33 Lh"1. Les fluides utilitaires ont une température de 15°C et un débit de 10800 L.h"1.
L'eau a ici été introduite dans la chambre de réaction en trois fois et en trois points d'alimentation différents. Cela permet de diminuer la température maximale du milieu réactionnel en dessous de la contrainte de 800C. Le procédé d'optimisation a permis d'identifier un dimensionnement approprié répondant à l'objectif de productivité et respectant la contrainte de température.
Les figures 11 et 12 sont des graphes du suivi de la réaction de production du néopentyl glycol. L'axe 165 des abscisses représente le volume de la chambre de réaction en litre. Dans cet exemple de réalisation, la chambre de réaction comprend 40 cellules élémentaires séparées en quatre blocs de 10 cellules élémentaires ayant chacune un volume de
0,035L environ. L'axe 166 des ordonnés à gauche représente la température en degré Celsius et l'axe 167 des ordonnés à droite représente le rendement de la réaction en pourcentage.
La courbe 168 est l'évolution du rendement de la réaction, les courbe 169, 170 et 171 sont respectivement l'évolution de la température du milieu réactionnel, des fluides utilitaires et des plaques délimitant la chambre de réaction.
L'équation de production du néopentyl glycol est la suivante : 2 (Formaldéhyde) + (2-méthyl propanai) + (hydr oxyde de sodium) — > (néopentyl glycol) + (formate de sodium)
Aux figures 11 et 12, les lignes verticales en pointillés 172 et 173 représentent respectivement les points d'alimentation de l'hydroxyde de sodium et du 2-méthyl propanal. L'hydroxyde de sodium est donc introduit en une seule fois lorsque le volume de la chambre de réaction est proche de 0.05L, ce qui correspond à la deuxième cellule de la chambre de réaction (0.035L par cellule élémentaire), et le 2-méthyl propanai est introduit en trois fois, dans la troisième cellule, la quinzième cellule et la vingt-huitième cellule.
En figure 11 , la réaction étant exothermique, les fluides utilitaires ont une température de 500C et un débit de 1 ,5 m3.h'1. Les réactifs sont injectés à une température de 18°C, et le 2-méthyl propanai a une concentration de 2 mol. kg"1 pour un débit total de 9,14 Lh"1 (trois injections à des débits identiques de 3,05 Lh"1 environ). Le taux de conversion atteint est de 89,4% et la température maximale du milieu réactionnel au cours de la réaction est 66,6°C.
L'opérateur fixe un objectif d'augmentation du rendement de la réaction de production du néopentyl glycol, des données variables d'optimisation du fonctionnement de la réaction qui sont les débits d'alimentation du 2-méthyl propanai pour chaque point d'alimentation, et une contrainte à respecter qui est une température limite du milieu réactionnel à ne pas dépasser, qui est de 65°C.
Le résultat est représenté en figure 12, où les fluides utilitaires ont une température de 500C et un débit de 1 ,5 m3. h"1, les réactifs sont alimentés à une température de 18°C, et le 2-méthyl propanai a une concentration de 2 mol. kg"1 pour un débit total de 9,14 L.h"1. Les trois injections de 2-méthyl propanai ont des débits différents : de 5,91 L.h"1 pour la première alimentation en cellule 3, de 2,73 L.h"1 pour la deuxième alimentation en cellule 11 et de 0,50 L.h'1 pour la troisième alimentation en cellule 21. Le taux de conversion atteint est de 95,6% et la température maximale du milieu réactionnel au cours de la réaction atteint 650C et respecte la contrainte fixée.
Le procédé d'optimisation a donc permis d'atteindre un objectif de rendement d'une réaction et de respecter une contrainte de température en agissant sur des données variables de réaction qui sont les débits d'alimentation d'un réactif en plusieurs points de la chambre de réaction.

Claims

REVENDICATIONS
1. Procédé d'optimisation d'une réaction chimique dans un réacteur, caractérisé en ce qu'il consiste à : - utiliser un réacteur ouvert constitué d'un empilement de plaques (20, 21 , 22, 23) délimitant entre elles au moins un bloc comprenant une chambre de réaction (24) formée entre deux chambres latérales (27, 28) d'échange thermique à circulation de fluide utilitaire d'échange de chaleur, des moyens d'alimentation de la chambre de réaction (24) par un débit continu d'un ou de plusieurs réactifs et des moyens d'alimentation des deux chambres latérales (27, 28) par un débit continu de fluide utilitaire ;
- établir un modèle dynamique (36) du réacteur pour une réaction donnée, à partir d'une modélisation de la chambre de réaction (24) et des échanges thermiques entre cette chambre de réaction (24) et les chambres latérales (27, 28) à circulation de fluide utilitaire, ce modèle dynamique (36) comprenant des équations de bilans de matière et d'énergie et de contraintes ;
- appliquer un outil logiciel d'intégration (37) à ce modèle dynamique (36) pour résoudre les équations précitées ; - déterminer et optimiser un ensemble de paramètres dimensionnels et/ou fonctionnels du réacteur à partir d'évaluations des objectifs (47) à atteindre, des contraintes (48) à respecter, et des données variables du réacteur (49) et/ou de la réaction (53) ;
- construire un réacteur à plaques du type précité conforme audit ensemble optimisé de paramètres ; et
- faire des mesures de paramètres physiques tels que la température et Ia pression lors de la mise en œuvre de ladite réaction dans le réacteur construit, pour valider son fonctionnement.
2. Procédé selon la revendication 1 , caractérisé en ce qu'il consiste, après la construction du réacteur, à : - établir un nouveau modèle dynamique de réacteur à partir du réacteur construit en procédant comme précité ;
- appliquer l'outil logiciel d'intégration (37) à ce nouveau modèle dynamique pour résoudre les équations précitées ; - optimiser à nouveau ledit ensemble de paramètres dimensionnels et /ou fonctionnels du réacteur ; et
- modifier si nécessaire ledit réacteur construit, conformément au nouvel ensemble optimisé de paramètres.
3. Procédé selon la revendication 1 ou 2, caractérisé en ce que la modélisation de la chambre de réaction (24) comprend une décomposition de cette chambre en une série de cellules élémentaires successives contenant chacune un milieu fluide parfaitement agité.
4. Procédé selon la revendication 3, caractérisé en ce que le milieu fluide dans chaque cellule élémentaire est homogène et son état et ses évolutions sont définis par des équations de bilan de matière, de bilan enthalpique, de bilan de pression et de contraintes de volume.
5. Procédé selon la revendication 3, caractérisé en ce que le milieu fluide contenu dans chaque cellule élémentaire de la chambre de réaction (24) est un milieu diphasique comprenant une phase continue et une phase dispersée.
6. Procédé selon Ia revendication 5, caractérisé en ce que l'état et l'évolution du milieu diphasique sont définis dans chaque cellule élémentaire à partir d'équations d'estimation de la masse volumique, de la chaleur spécifique, de la conductivité thermique et de la viscosité du milieu diphasique, et d'équations de bilan de matière, de bilan enthalpique, de bilan de pression et de contraintes de volume.
7. Procédé selon l'une des revendications 3 à 6, caractérisé en ce que la modélisation des échanges thermiques comprend une décomposition de chacune des chambres latérales (27, 28) en une série de cellules élémentaires contenant chacune un milieu fluide parfaitement agité, le nombre de cellules élémentaires de chaque chambre latérale (27, 28) étant égal au nombre de cellules élémentaires de la chambre de réaction (24), et une définition du sens de circulation de fluide utilitaire dans chaque cellule élémentaire par rapport au sens de passage du fluide réactionnel dans la cellule élémentaire correspondante de la chambre de réaction (24)
8. Procédé selon la revendication 7, caractérisé en ce que la modélisation des échanges thermiques comprend également des équations de bilan thermique relatives au fluide réactionnel, aux plaques (20, 21 , 22, 23) délimitant les chambres latérales (27, 28) d'échanges thermiques, et aux fluides utilitaires circulant dans les chambres latérales (27, 28), et des équations de bilan de matière des fluides utilitaires dans les différentes cellules élémentaires précitées des chambres latérales (27, 28).
9. Procédé selon l'une des revendications précédentes, caractérisé en ce que le modèle dynamique du réacteur comprend également des modèles des zones d'entrée et de sortie du réacteur, des zones d'entrée et de sortie de chaque chambre de réaction (24) et des zones de transition entre les différents blocs du réacteur.
10. Procédé selon l'une des revendications 3 à 9, caractérisé en ce qu'il consiste à modéliser le comportement hydrodynamique du réacteur à partir d'études expérimentales de distribution de temps de séjour des fluides dans les différentes parties du réacteur, pour définir le nombre de cellules élémentaires à contenu parfaitement agité.
11. Procédé selon l'une des revendications 3 à 10, caractérisé en ce qu'il comprend également une modélisation de la réaction à partir d'équations de vitesse de Ia réaction, de taux de production de constituants de la réaction et de génération de chaleur, dans chaque cellule élémentaire de la ou des chambres de réaction (24).
12. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend également une modélisation du transfert de matière, dans une réaction en milieu diphasique, entre les deux phases du milieu réactionnel, à partir de propriétés physiques du milieu diphasique et de la taille des gouttelettes de la phase dispersée.
13. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend également une modélisation du transfert thermique à partir de définitions des coefficients de transfert thermique entre les plaques délimitant les chambres latérales (27, 28) et les fluides circulant dans la chambre de réaction (24) et dans les chambres latérales (27, 28), et d'estimations des coefficients de film du fluide réactionnel et de fluide utilitaire.
14. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend également une modélisation des pertes de charge du fluide réactionnel dans le réacteur, faite à partir de mesures expérimentales.
15. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il comprend également une estimation des propriétés physiques des composants du réacteur et des fluides réactionnels et utilitaires.
16. Procédé selon l'une des revendications 3 à 15, caractérisé en ce qu'il consiste à optimiser le dimensionnement du réacteur, pour des objectifs (47) et des contraintes (48) données, par ajustement de paramètres tels que le nombre de séries de cellules élémentaires, le nombre de points d'alimentation de réactifs dans la chambre de réaction (24), le sens de passage du fluide utilitaire par rapport au milieu réactionnel, la nature et Ia distribution du fluide utilitaire, le volume de la chambre de réaction (24) et/ou les volumes des chambres latérales (27, 28) de circulation de fluide utilitaire.
17. Procédé selon l'une des revendications précédentes, caractérisé en ce qu'il consiste à optimiser le fonctionnement du réacteur, pour des objectifs (47) et des contraintes (48) données, par ajustement de paramètres opératoires comprenant les températures, les pressions, les compositions et/ou les débits des fluides réactionnels et des fluides utilitaires.
18. Procédé selon la revendication 17, caractérisé en ce que, le réacteur ou la chambre de réaction (24) comprenant plusieurs zones d'entrée de fluides, il consiste à ajuster les débits, les températures, les pressions et/ou les compositions d'injection des fluides dans ces différentes zones.
19. Procédé selon l'une des revendications 5 à 18, caractérisé en ce que le milieu réactionnel diphasique est considéré comme un milieu pseudo-homogène et les propriétés physiques des deux phases sont déterminées par des lois d'association reliant les propriétés physiques du milieu diphasique aux propriétés de chacune des phases ainsi qu'à leurs proportions.
20. Procédé selon la revendication 19, caractérisé en ce qu'il consiste à appliquer au mélange de phases des règles d'équilibre entre les deux phases en fonction de la dimension des gouttelettes de la phase dispersée, des propriétés thermodynamiques des phases en présence, et/ou des conditions d'écoulement.
21. Procédé selon la revendication 19 ou 20, caractérisé en ce qu'il consiste à suivre l'évolution de la taille des gouttelettes de la phase dispersée dans les cellules élémentaires de la ou des chambres de réaction (24).
22. Procédé selon l'une des revendications 16 à 21 , caractérisé en ce qu'il consiste également à définir des valeurs limites des paramètres opératoires satisfaisant à des contraintes de sécurité et/ou environnementales.
PCT/FR2006/000049 2005-01-14 2006-01-10 Optimisation d’une reaction chimique dans un reacteur ouvert du type a plaques WO2006075085A2 (fr)

Priority Applications (9)

Application Number Priority Date Filing Date Title
MX2007008567A MX2007008567A (es) 2005-01-14 2006-01-10 Optimizacion de una reaccion quimica en un reactor abierto del tipo de placas.
CA2594822A CA2594822C (fr) 2005-01-14 2006-01-10 Optimisation d'une reaction chimique dans un reacteur ouvert du type a plaques
CN2006800022907A CN101288078B (zh) 2005-01-14 2006-01-10 最优化板型敞口反应器中的化学反应的方法
JP2007550813A JP5107054B2 (ja) 2005-01-14 2006-01-10 開放プレート型反応器における化学反応の最適化
AU2006205731A AU2006205731B2 (en) 2005-01-14 2006-01-10 Optimisation of a chemical reaction in an open plate-type reactor
BRPI0606194-0A BRPI0606194A2 (pt) 2005-01-14 2006-01-10 processo de otimização de uma reação quìmica em um reator
EP06709060A EP1836627A2 (fr) 2005-01-14 2006-01-10 Optimisation d'une reaction chimique dans un reacteur ouvert du type à plaques
US11/773,667 US8065038B2 (en) 2005-01-14 2007-07-05 Optimizing a chemical reaction in a plate-type open reactor
NO20074160A NO20074160L (no) 2005-01-14 2007-08-13 Optimalisering av en kjemisk reaksjon i en åpen reaktor av platetypen

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0500424 2005-01-14
FR0500424A FR2880967A1 (fr) 2005-01-14 2005-01-14 Optimisation d'une reaction chimique dans un reacteur ouvert du type a plaques

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US11/773,667 Continuation US8065038B2 (en) 2005-01-14 2007-07-05 Optimizing a chemical reaction in a plate-type open reactor

Publications (2)

Publication Number Publication Date
WO2006075085A2 true WO2006075085A2 (fr) 2006-07-20
WO2006075085A3 WO2006075085A3 (fr) 2008-06-19

Family

ID=35045034

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/FR2006/000049 WO2006075085A2 (fr) 2005-01-14 2006-01-10 Optimisation d’une reaction chimique dans un reacteur ouvert du type a plaques

Country Status (12)

Country Link
US (1) US8065038B2 (fr)
EP (1) EP1836627A2 (fr)
JP (1) JP5107054B2 (fr)
CN (1) CN101288078B (fr)
AU (1) AU2006205731B2 (fr)
BR (1) BRPI0606194A2 (fr)
CA (1) CA2594822C (fr)
FR (1) FR2880967A1 (fr)
MX (1) MX2007008567A (fr)
NO (1) NO20074160L (fr)
RU (1) RU2417119C2 (fr)
WO (1) WO2006075085A2 (fr)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008051143A1 (fr) * 2006-10-27 2008-05-02 Perstorp Specialty Chemicals Ab Procédé de production de diol, triol, polyol dans un ou plusieurs réacteur(s) à plaque
JP2011517454A (ja) * 2008-04-03 2011-06-09 サントル ナスィオナル ド ラ ルシェルシュ スィアンティフィク（セ．エン．エル．エス．） イソシアネートの連続的オリゴマー化の方法

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070059713A1 (en) 2005-09-09 2007-03-15 Lee Jun E SSB-DNA polymerase fusion proteins
US8150560B2 (en) * 2006-06-23 2012-04-03 Saudi Arabian Oil Company Methods for heat exchanger network energy efficiency assessment and lifetime retrofit
US8150559B2 (en) * 2006-06-23 2012-04-03 Saudi Arabian Oil Company Systems and program product for heat exchanger network energy efficiency assessment and lifetime retrofit
US7729809B2 (en) 2006-06-23 2010-06-01 Saudi Arabian Oil Company System, method, and program product for targeting and identification of optimal process variables in constrained energy recovery systems
US8116920B2 (en) 2009-10-08 2012-02-14 Saudi Arabian Oil Company System, method, and program product for synthesizing non-thermodynamically constrained heat exchanger networks
US8150668B2 (en) * 2009-02-11 2012-04-03 Livermore Software Technology Corporation Thermal fluid-structure interaction simulation in finite element analysis
FR2971783B1 (fr) 2011-02-17 2013-02-15 Rhodia Operations Procede d'hydroxylation de phenols et d'ethers de phenols
EP2874029A1 (fr) * 2013-11-15 2015-05-20 Bayer Technology Services GmbH Procédé de fonctionnement d'une installation conçue pour réaliser au moins une réaction chimique
US10238992B2 (en) * 2016-05-03 2019-03-26 Saudi Arabian Oil Company Processes for analysis and optimization of multiphase separators, particularly in regard to simulated gravity separation of immiscible liquid dispersions
CN107563051B (zh) * 2017-08-30 2019-04-02 南京大学 微界面强化反应器气泡尺度构效调控模型建模方法
FR3071750B1 (fr) * 2017-09-29 2022-01-28 Atmostat Reacteur-echangeur a vieillissement maitrise
KR102276219B1 (ko) * 2019-02-15 2021-07-12 씨제이제일제당 (주) 생물반응기의 운전 조건을 결정하는 장치 및 방법
CN113283199B (zh) * 2021-06-28 2022-04-22 中国人民解放军国防科技大学 含相变的空气预冷器设计方法、装置、计算机系统及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002085511A1 (fr) * 2001-04-25 2002-10-31 Alfa Laval Vicarb Dispositif perfectionne d'echange et/ou de reaction entre fluides

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5223864A (en) * 1989-09-06 1993-06-29 J. D. Moller Optische Werke Gmbh Phoropter
JP2548880B2 (ja) * 1992-04-23 1996-10-30 インターナショナル・ビジネス・マシーンズ・コーポレイション 空間分割型確率シミュレーション・システムおよび物質エネルギー・システムをシミュレートする方法
JP3573602B2 (ja) * 1997-07-11 2004-10-06 株式会社荏原製作所 ポンプ試験装置
US6524482B2 (en) * 2001-04-19 2003-02-25 Ibc Advanced Technologies, Inc. Use of ion binding ligands attached to solid supports and membranes for ion removal from a biological system
JP2002356682A (ja) * 2001-05-30 2002-12-13 Mitsubishi Materials Corp 石炭ガス化をシミュレーションするためのモデリング方法および石炭ガス化のシミュレーションシステム
JP2004051654A (ja) * 2002-07-16 2004-02-19 Mitsubishi Chemicals Corp 処理槽の設計方法及びその処理槽並びに処理槽を用いた重合樹脂粒子の製造方法及びその重合樹脂粒子
DE10307233A1 (de) * 2003-02-14 2004-08-26 Basf Ag Verfahren und Vorrichtung zur Simulation von stoffübergangskontrollierten Trennprozessen

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002085511A1 (fr) * 2001-04-25 2002-10-31 Alfa Laval Vicarb Dispositif perfectionne d'echange et/ou de reaction entre fluides

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ANDERSSON R ET AL: "Development of a multi-scale simulation method for design of novel multiphase reactors" CHEMICAL ENGINEERING SCIENCE, OXFORD, GB, vol. 59, no. 22-23, novembre 2004 (2004-11), pages 4911-4917, XP004668399 ISSN: 0009-2509 *
CALVET J-P ET AL: "MODELLING, OPTIMISATION AND CONTROL OF A CONTINUOUS MULTISTAGE REACTOR FOR AN INDUSTRIAL REACTIVE SYSTEM" JOURNAL A, SOFT VISION, BRUSSELS, BE, vol. 35, no. 1, 1 avril 1994 (1994-04-01), pages 5-15, XP000460882 ISSN: 0771-1107 *
ZHANG J ET AL: "Design and optimisation of batch and semi-batch reactors" CHEMICAL ENGINEERING SCIENCE, OXFORD, GB, vol. 59, no. 2, janvier 2004 (2004-01), pages 459-478, XP004484748 ISSN: 0009-2509 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2008051143A1 (fr) * 2006-10-27 2008-05-02 Perstorp Specialty Chemicals Ab Procédé de production de diol, triol, polyol dans un ou plusieurs réacteur(s) à plaque
JP2011517454A (ja) * 2008-04-03 2011-06-09 サントル ナスィオナル ド ラ ルシェルシュ スィアンティフィク（セ．エン．エル．エス．） イソシアネートの連続的オリゴマー化の方法

Also Published As

Publication number Publication date
AU2006205731A1 (en) 2006-07-20
JP2008529757A (ja) 2008-08-07
CA2594822C (fr) 2014-07-08
AU2006205731B2 (en) 2009-11-12
US8065038B2 (en) 2011-11-22
FR2880967A1 (fr) 2006-07-21
CN101288078B (zh) 2012-06-13
RU2007124076A (ru) 2009-02-20
BRPI0606194A2 (pt) 2009-06-09
MX2007008567A (es) 2008-03-10
EP1836627A2 (fr) 2007-09-26
RU2417119C2 (ru) 2011-04-27
JP5107054B2 (ja) 2012-12-26
CN101288078A (zh) 2008-10-15
CA2594822A1 (fr) 2006-07-20
WO2006075085A3 (fr) 2008-06-19
US20080140376A1 (en) 2008-06-12
NO20074160L (no) 2007-08-13

Similar Documents

Publication Publication Date Title
WO2006075085A2 (fr) Optimisation d’une reaction chimique dans un reacteur ouvert du type a plaques
Dixon Local transport and reaction rates in a fixed bed reactor tube: Endothermic steam methane reforming
Filtvedt et al. Chemical vapor deposition of silicon from silane: Review of growth mechanisms and modeling/scaleup of fluidized bed reactors
Mehrabi et al. Modeling of heat transfer and fluid flow characteristics of helicoidal double-pipe heat exchangers using adaptive neuro-fuzzy inference system (ANFIS)
Finaud-Guyot et al. 1D–2D coupling for river flow modeling
CN104634708B (zh) 基于计算流体力学流化床内颗粒密度和粒径分布预测方法
Ceballos et al. Slow invasion of a nonwetting fluid from multiple inlet sources in a thin porous layer
Hosseini et al. Smart models to predict the minimum spouting velocity of conical spouted beds with non-porous draft tube
Seitz Modeling fixed-bed reactors for thermochemical heat storage with the reaction system CaO/Ca (OH) 2
Mustafa et al. Numerical analysis and moth flame optimization of passive t-micromixer with twist and bend mixing channel
Wang et al. Numerical simulation of biomass steam gasification in an internally interconnected fluidized bed using a two-grid MP-PIC model
Gamba et al. Residence time distribution determination of a continuous stirred tank reactor using computational fluid dynamics and its application on the mathematical modeling of styrene polymerization
Xiao et al. 2D numerical modelling of meandering channel formation
Wang et al. Multi-objective optimization of rectangular cooling channel design using Design of Experiments (DOE)
Sui et al. A universal ANN-based approach predicting PCHEs’ off-design performance across various operating conditions of sCO2 RCBCs
Zoljalali et al. Developing MLP-ICA and MLP algorithms for investigating flow distribution and pressure drop changes in manifold microchannels
Ponomarev et al. Modelling non-homogeneous flow and residence time distribution in a single-screw extruder by means of Markov chains
Xu et al. Enhancement of mixing by different baffle arrays in cavity flows
Artikov et al. Computer modeling of the heat exchange process in the calefacient chamber of the evaporator
Biswas et al. Two-dimensional Simulation of Droplet Splitting in Multi-Furcating Microchannel
Vismara et al. Computational fluid dynamics calculations of high efficiency heat exchangers operating in laminar motion consisting of extruded plastic profiles with inserts
Krause Development of a design for the reduction chamber of a swept, open, moving particle reactor for solar thermochemical fuel production
Petersson Experiments and simulations of gas flow through a NiCd battery
Rautenbach Modelling of flow through porous packing elements of a CO2 absorption tower
Aminova et al. Investigation of simultaneous chemical-transformation and heat-exchange process occurring in the process of synthesis of divinyl rubber on the basis of a cobaltous catalyst in a cascade of two continuous-operation reactors

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200680002290.7

Country of ref document: CN

121 Ep: the epo has been informed by wipo that ep was designated in this application
REEP Request for entry into the european phase

Ref document number: 2006709060

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2006709060

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 5194/DELNP/2007

Country of ref document: IN

WWE Wipo information: entry into national phase

Ref document number: MX/a/2007/008567

Country of ref document: MX

Ref document number: 2594822

Country of ref document: CA

Ref document number: 2007550813

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2006205731

Country of ref document: AU

WWE Wipo information: entry into national phase

Ref document number: 2007124076

Country of ref document: RU

ENP Entry into the national phase

Ref document number: 2006205731

Country of ref document: AU

Date of ref document: 20060110

Kind code of ref document: A

WWP Wipo information: published in national office

Ref document number: 2006205731

Country of ref document: AU

WWP Wipo information: published in national office

Ref document number: 2006709060

Country of ref document: EP

ENP Entry into the national phase

Ref document number: PI0606194

Country of ref document: BR

Kind code of ref document: A2