L’objectif de ce travail est de proposer de nouvelles stratégies de commande non linéaire pour la stabilisation des Réacteurs Parfaitement Agités Continus (RPAC). Pour cela, nous utilisons d’une part, l’approche thermodynamique entropique. Plus précisément, nous utilisons la notion de disponibilité thermodynamique et les propriétés de la thermodynamique irréversible pour définir une fonction de Lyapunov utilisable pour la stabilisation du système en boucle fermée. Nous proposons aussi une fonction disponibilité réduite afin d’obtenir des lois stabilisantes plus performantes en terme de sollicitation des actionneurs. D’autre part, nous proposons une extension du formalisme (pseudo) hamiltonien à ports dissipatifs aux réacteurs chimiques ouverts. Nous montrons que l’Hamiltonien est lié à l’enthalpie libre de Gibbs dans le cas isotherme et à l’ectropie (opposée de l’entropie) dans le cas non isotherme. Par ce formalisme, la dissipation du système représente la production irréversible d’entropie due à la réaction chimique. Nous appliquons ensuite les techniques de commande passive (modelage de l’énergie) pour la synthèse de lois de commande en choisissant la disponibilité thermodynamique comme fonction hamiltonienne à modeler en boucle fermée. Finalement, nous montrons que les commandes synthétisées par l’approche thermodynamique entropique et la formulation pseudo-hamiltonienne sont, dans certains cas, équivalentes. Certaines propriétés relatives à la stabilisation et l’admissibilité des commandes sont aussi considérées. Les développements théoriques sont mis en oeuvre sur des exemples différents de RPAC : un réacteur académique et l’hydrolyse par catalyse acide de l’oxirane-méthanol en glycérine. / The goal of this thesis is to propose new nonlinear control strategies for the stabilization of perfectly Continuous Stirred Tank Reactors (CSTR). To achieve this goal, we use on the one hand, the entropic thermodynamic approach. More precisely, we use the thermodynamic availability concept and the properties of irreversible thermodynamics to define a Lyapunov function candidate for the stabilization of the closed loop system. We also propose a reduced availability function to design more efficient feedback laws in term of control variable solicitations. On the other hand, we propose an extension of the (pseudo) Hamiltonian formalism associated to dissipative systems to open chemical reactors. We show that the Hamiltonian is linked to the Gibbs free enthalpy in the isothermal case and to ectropy (opposed to entropy) in the non isothermal case. By this formalism, the dissipation of the system represents the irreversible entropy production due to chemical reaction. The Interconnection and Damping Assignment-Passivity Based Control (IDA-PBC) approach is then applied to synthesize feedback laws by choosing the thermodynamic availability as desired closed loop hamiltonian storage function. Finally, we show that feedback laws synthetized by the entropic thermodynamic approach and the pseudo-hamiltonian formulation are equivalent in some cases. Some stabilization properties and the control input admissibility are also considered. Theoretical developments are illustrated on some different CSTR examples : an academic case study and the acid catalyzed hydration of oxirane-methanol to glycerol.

L’objet de cette thèse est la modélisation et la commande par approche thermodynamique des réacteurs catalytiques triphasiques en mode continu et en mode discontinu. Ce type de réacteur consiste en un système fortement non linéaire, multivariable et siège de réactions exothermiques. Nous utilisons les concepts de la thermodynamique irréversible pour la synthèse de lois de commande stabilisante pour ces deux types de réacteurs chimiques. En effet, la stricte concavité de la fonction d’entropie nous a permis de définir une fonction de stockage qui sert de fonction de Lyapunov candidate : la disponibilité thermodynamique. Nous utilisons cette fonction de disponibilité thermodynamique pour la synthèse de lois de commande stabilisante d’un mini-réacteur catalytique triphasique intensifié continu. Une stratégie de contrôle à deux couches (optimisation et contrôle) est utilisée pour contrôler la température et la concentration du produit à la sortie du réacteur en présence de perturbations à l’entrée du réacteur. Les performances du contrôleur mis en place sont comparées en simulation à celles d’un régulateur PI. Dans certains cas, l’utilisation de la fonction de disponibilité thermodynamique s’avère problématique. Une autre étude effectuée sur cette fonction nous permet de déterminer une nouvelle fonction de Lyapunov : la disponibilité thermique. Nous utilisons par la suite la fonction de disponibilité thermique pour la synthèse de lois de commande stabilisante d’un réacteur catalytique triphasique semi-fermé. Un observateur grand gain est utilisé pour estimer la vitesse de réaction à partir des mesures de la température du milieu réactionnel. Cette estimation est injectée ensuite dans le calcul de la loi de commande mise en place. La robustesse du schéma de contrôle est testée en simulation face à des incertitudes de modélisation, des perturbations et des bruits de mesure. / The goal of this thesis is the modeling and thermodynamic based control of three-phase catalytic reactor working in continuous or discontinuous modes. These types of reactors are highly nonlinear, multivariable and exothermal processes. We use the concepts of irreversible thermodynamics for the synthesis of stabilizing control laws for these two types of chemical reactors. Indeed, the strict concavity of the entropy function has allowed us to define a storage function used as a candidate Lyapunov function: the availability function. We use this availability function for the synthesis of control laws for stabilizing a three-phase catalytic continuous intensified mini-reactor. A control strategy with two layers (optimization and control) is used to control the temperature and concentration of the product at the outlet of the reactor in the presence of disturbances. The performances of the controller are compared by simulation to those of a PI controller. In some cases, the use of the availability function may cause some problems. A new Lyapunov candidate function is then derived from the original availability function: the thermal availability. We use this thermal availability for the synthesis of stabilizing control laws for a three-phase catalytic fed-batch reactor. A high gain observer is used to estimate the chemical reaction rate from the measurements of the temperature. This estimate is then used for the control law implementation. The robustness of the control scheme is tested in simulation against modelling uncertainties, disturbances and noise measurements.

Cette thèse s’insère dans le cadre général de l’étude de la dynamique des populations. La population prise en compte étant constituée de cellules souches normales et/ou cancéreuses. Nous proposons et analysons trois modèles mathématiques décrivant la dynamique de cellules souches. Le premier modèle proposé est un modèle d’équations aux dérivées partielles structurées en âge que nous transformons, via la méthode des caractéristiques, en un système d'équations différentielles à retard pour lequel on étudie l'existence et la stabilité des points d'équilibres. On effectue, après, des simulations numériques permettant d'illustrer le comportement des états d'équilibres. Dans le deuxième modèle, on considère que la durée du cycle cellulaire dépend de la population totale de cellules quiescentes. La méthode des caractéristiques nous permet de réduire notre modèle structuré en âge à un système d'équations différentielles avec un retard dépendant de l'état pour lequel on effectue une analyse détaillée de la stabilité. Nous confirmons, ensuite, les résultas analytiquement obtenus par des simulations numériques. Pour le troisième et dernier modèle de cette thèse, on propose un système d'équations différentielles ordinaires décrivant la dynamique de cellules souches saines et cancéreuses et prenant en compte leurs interactions avec les réponses immunitaires. Ce modèle nous a permis de souligner l'ampleur de l'impact que peuvent avoir différentes infections sur la prolifération tumoral que ce soit par le biais de leurs fréquences, leurs durées ou la façon dont ils agissent sur le système immunitaire. / This thesis fits into the general framework of the study of population dynamics. The population particularly considered in this work is comprised of stem cells with both cases of healthy and cancerous cells being investigated. We propose and analyze three mathematical models describing stem cells dynamics. The first model is an age-structured partial differential model that we reduce to a delay differential system using the characteristics method. We investigate the existence and stability of the steady states of the reduced delay differential system. We, then, conduct some numerical simulations to illustrate the behavior of the steady states. In the second model, the duration of the cell cycle is considered to depend upon the total population of quiescent cells. The method of characteristics reduces the age-structured model to a system of differentialequations with a state-dependent delay. We perform a detailed stability analysis of the resulting delay differential system. We confirm the analytical results by numerical simulations. The third and final model, proposed in this thesis, is an ordinary differential equations model describing healthy and cancerous stem cells dynamics and their interactions with immune system responses. Through this model, we show that the frequency, the duration of infections and their action (positive or negative) on immune responses may impact significantly tumor proliferation.

Cette thèse s’intéresse à l’étude de la stabilité des systèmes à commutation qui évoluent sur un domaine de temps non uniforme en introduisant la théorie des échelles de temps. On s’intéresse essentiellement aux systèmes dynamiques linéaires à commutation définis sur une échelle de temps particulière T = P{tσk ,tk+1} = ∪∞k=0[tσk , tk+1]. Le système étudié commute entre un sous-système dynamique continu sur les intervalles ∪∞k=0[tσk , tk+1[ et un sous-système dynamique discret aux instants ∪∞k=0{tk+1} (à temps discret) avec un pas discret qui varie dans le temps. Dans une première partie, des conditions suffisantes sont données pour garantir la stabilité exponentielle de cette classe de systèmes à commutation. Ensuite, des conditions nécessaires et suffisantes de stabilité sont données en déterminant une région de stabilité exponentielle. Dans une deuxième partie, la stabilité de cette classe des systèmes à commutation avec des perturbations nonlinéaires a été traitée en utilisant des majorations de la solution, puis en introduisant l’approche de la fonction de Lyapunov commune. La troisième partie est consacrée au problème du consensus en présence d’interruptions de transmission d’informations où le système multi-agent en boucle fermée peut être représenté comme un système à commutation par une combinaison de modèles de systèmes linéaires à temps continu et de systèmes linéaires à temps discret. / This thesis deals with the stability analysis of switched systems that evolve on non uniform time domain by introducing the time scale theory. We are interested mainly in dynamical linear switched systems defined on particular time scale T = P{tσk ,tk+1} = ∪∞k=0[tσk, tk+1]. The studied system switches between a continuous-time dynamical subsystem on the intervals ∪∞k=0[tσk, tk+1[ and a discrete-time dynamical subsystem on instants ∪∞k=0{tk+1} (a discrete time) with a time-varying discrete step. In a first part, sufficient conditions are given to guarantee the exponential stability of this class of switched systems. Then necessary and sufficient conditions for stability are given by determining a region of exponential stability. In the second part, the stability of this class of switched systems with nonlinear uncertainties, is treated using majoration of the solution, and after that by introducing the approach of a common Lyapunov function. The third part is devoted to the consensus problem under intermittent information transmissions where the closed-loop multi-agent system can be represented as a switched system using a combination of linear continuous-time and linear discrete-time systems.

Dans cette thèse, nous prouvons des résultats concernant le contrôle et la stabilisation d'équations dispersives étudiées sur un intervalle borné. Pour commencer, nous étudions la stabilisation interne du système de Gear-Grimshaw, qui est un système de deux équations de Korteweg-de-Vries (KdV) couplées. Nous obtenons une décroissance exponentielle de l'énergie totale associée au modèle en introduisant une fonction de Lyapunov convenable. Nous prouvons aussi des résultats de contrôlabilité à zéro et exacte pour l'équation de Korteweg-de Vries avec un contrôle distribué à support dans un sous-intervalle du domaine. Pour la contrôlabilité à zéro du système linéarisé, nous utilisons l'approche classique basée sur la dualité qui ramène le problème à l'étude d'une inégalité d'observabilité qui, dans ce travail, est établie à l'aide d'une inégalité de Carleman. Ensuite, utilisant des fonctions plateau, nous prouvons un résultat de contrôlabilité exacte. Dans les deux cas, le résultat concernant le système non linéaire est obtenu à l'aide d'un argument de point fixe. Enfin, dans la lignée du résultat de contrôlabilité au bord obtenu par L. Rosier pour KdV, nous prouvons que le système linéaire de Boussinesq de type KdV-KdV est exactement contrôlable lorsque des contrôles sont appliqués au bord. Notre méthode repose sur l'utilisation de multiplicateurs et l'approche de la dualité mentionnée ci-dessus. Lorsqu'un mécanisme d'amortissement est introduit au bord, nous montrons que le système non linéaire est aussi exactement contrôlable et que l'énergie associée au modèle décroit exponentiellement / This work is devoted to prove a series of results concerning the control and stabilization properties of dispersive models posed on a bounded interval. Initially, we study the internal stabilization of a coupled system of two Korteweg-de Vries equations (KdV), the so-called Gear-Grimshaw system. Defining a convenient Lyapunov function we obtain the exponential decay of the total energy associated to the model. We also prove results of null and exact controllability for the Korteweg-de Vries equation with a control acting internally on a subset of the domain. In the case of the null controllability for the linear model, we use a classical duality approach which reduces the problem to the study of an observability inequality that, in this work, is proved by means of a Carleman inequality. Then, making use of cut-off functions, the exact controllability is also investigated. In both cases, the result for the nonlinear system is obtained by means of fixed-point argument. Finally, in view of the result of the boundary controllability obtained by L. Rosier for the KdV equation, we prove that the linear Boussinesq system of KdV-KdV type is exactly controllable when the controls act in the boundary conditions. Our analysis is performed using multipliers and the duality approach mentioned above. Adding a damping mechanism in the boundary, it is proved that the nonlinear system is also exactly controllable and that the energy associated to the model decays exponentially

Ce mémoire aborde la synthèse de filtres H-infine d'ordre plein et d'ordre réduit pour les systèmes stochastiques à temps continu avec bruits multiplicatifs. Les bruits considérés dans l'équation d'état et dans l'équation de mesures sont des processus de Wiener.<br /><br />Les systèmes stochastiques étudiés dans ce mémoire sont écrits sous la forme d'une équation différentielle stochastique au sens d'Itô dans lesquels la dérive et la diffusion sont linéaires ou bilinéaires. Les systèmes avec plusieurs bruits multiplicatifs et les systèmes dont les mesures sont affectées par des bruits multiplicatifs sont également traités dans ce mémoire. La conception d'une commande H-infine basée sur un observateur d'ordre réduit pour les systèmes stochastiques incertains est étudiée.<br /><br />Le critère de performance considéré est le critère H-infine du signal de perturbation vers le signal d'erreur d'estimation. La stabilité retenue pour ces systèmes stochastiques dans ce travail est la stabilité exponentielle en moyenne quadratique.<br /><br />La méthode utilisée pour trouver les matrices des filtres est basée sur l'utilisation de la théorie de Lyapunov pour les équations différentielles stochastiques, la formule d'Itô et sur la résolution des Inégalités Matricielles Affines couplées à des contraintes bilinéaires qui assurent la stabilité et la performance.

Cette thèse s'intéresse au diagnostic et à l'observation de systèmes linéaires à commutations et à l'application au convertisseur multicellulaire série. L'objectif est de proposer des solutions pour des sous-systèmes non-observables au sens classique et dont des fautes continues ou discrètes peuvent être présentes. Après la présentation d'un état de l'art sur les techniques d'observation et de diagnostic pourles systèmes à commutations, le mémoire est scindé en deux parties. La première partie propose, d'une part, une stratégie d'estimation des états discret et continu d'un système linéaire à commutation soumis à une entrée inconnue. Un observateur hybride basé sur la théorie des modes glissants d'ordre supérieur est développé. D'autre part, deux procédures de diagnostic sont présentées. La première combine un observateur hybride et un diagnostiqueur pour détecter une faute continue. Pour la seconde, un diagnostic actif est défini sur la base de la théorie du test afin de détecter et d'isoler une faute discrète. Dans la seconde partie de ce mémoire, les étapes de la réalisation d'un convertisseur multicellulaire sont détaillées. Ensuite, un chapitre est dédié à la validation des approches théoriques d'observation et de diagnostic sur le convertisseur à trois cellules. Un observateur est synthétisé afin d'estimer les tensions des capacités. Les deux procédures de diagnostic sont appliquées pour la détection d'une variation des valeurs des capacités et le diagnostic de cellules bloquées. Enfin, une commande binaire pour le convertisseur est proposée. L'application de cette stratégie permettra, par la suite, la commande tolérante aux fautes du convertisseur.

Ma thèse de doctorat se concentre principalement sur le comportement en temps long des processus de Markov, les inégalités fonctionnelles et les techniques relatives. Plus spécifiquement, Je vais présenter les taux de convergence sous-exponentielle explicites des processus de Markov dans deux approches : la méthode Meyn-Tweedie et l'hypocoercivité (faible). Le document se divise en trois parties. Dans la première partie, Je vais présenter quelques résultats importants et des connaissances connexes. D'abord, un aperçu de mon domaine de recherche sera donné. La convergence exponentielle (ou sous-exponentielle) des chaînes de Markov et des processus de Markov (à temps continu) est un sujet d'actualité dans la théorie des probabilité. La méthode traditionnelle développée et popularisée par Meyn-Tweedie est largement utilisée pour ce problème. Dans la plupart des résultats, le taux de convergence n'est pas explicite, et certains d'entre eux seront brièvement présentés. De plus, la fonction de Lyapunov est cruciale dans l'approche Meyn-Tweedie, et elle est aussi liée à certaines inégalités fonctionnelles (par exemple, inégalité de Poincaré). Cette relation entre fonction de Lyapounov et inégalités fonctionnelles sera donnée avec les résultats au sens L2. En outre, pour l'exemple de l'équation cinétique de Fokker-Planck, un résultat de convergence exponentielle explicite de la solution sera introduite à la manière de Villani : l'hypocoercivité. Ces contenus sont les fondements de mon travail, et mon but est d'étudier la décroissance sous-exponentielle. La deuxième partie, fait l'objet d'un article écrit en coopération avec d'autres sur les taux de convergence sous-exponentielle explicites des processus de Markov à temps continu. Comme nous le savons, les résultats sur les taux de convergence explicites ont été donnés pour le cas exponentiel. Nous les étendons au cas sous-exponentielle par l'approche Meyn-Tweedie. La clé de la preuve est l'estimation du temps de passage dans un ensemble "petite", obtenue par Douc, Fort et Guillin, mais pour laquelle nous donnons une preuve plus simple. Nous utilisons aussi la construction du couplage et donnons une ergodicité sous exponentielle explicite. Enfin, nous donnons quelques applications numériques. Dans la dernière partie, mon second article traite de l'équation cinétique de Fokker-Planck. Je prolonge l'hypocoercivité à l'hypocoercivité faible qui correspond à inégalité de Poincaré faible. Grâce à cette extension, on peut obtenir le taux de convergence explicite de la solution, dans des cas sous-exponentiels. La convergence est au sens H1 et au sens L2. A la fin de ce document, j'étudie le cas de l'entropie relative comme Villani, et j'obtiens la convergence au sens de l'entropie. Enfin, Je donne deux exemples pour les potentiels qui impliquent l'inégalité de Poincaré faible ou l'inégalité de Sobolev logarithmique faible pour la mesure invariante.

Ce mémoire de thèse présente deux nouvelles approches pour l’analyse et la commande des systèmes non-linéaires complexes, comme les systèmes dynamiques à commutation de la classe des convertisseurs d’énergie électrique. Ces systèmes ont plusieurs modes de fonctionnement et ont un point de fonctionnement désiré qui, en général, n’est le point d’équilibre d’aucun des modes. Dans cette classe de systèmes, la commutation d’un mode de fonctionnement à un autre est commandée selon une loi qui doit être synthétisée. Par conséquent, la synthèse de commande implique l’étude des conditions qui permettent à un cycle limite stable de s’établir au voisinage du point de fonctionnement désiré, puis de la trajectoire de commande qui permet de l’atteindre en respectant les contraintes physiques de comportement (courant maximum supporté par les composants,. . .) ou les contraintes de temps (durée minimum entre deux commutations,. . .). Le cycle limite sera qualifié d’hybride car il est composé de plusieurs dynamiques(deux dans ces travaux).La première méthode développée s’appuie sur les propriétés géométriques des champs de vecteurs et est une extension d’une partie des travaux de thèse de Manon au LAGEP. Une condition nécessaire et suffisante d’existence et de stabilité d’un cycle limite hybride composé d’une séquence de deux modes de fonctionnement dans IR2 est présentée. Ce cycle définit la région finale à atteindre par le système depuis son état initial, par une trajectoire déterminée de manière optimale selon un critère donné (durée totale, énergie dépensée, . . .). La méthode proposée est appliquée aux convertisseurs d’énergie Buck et Buck-Boost alimentant une charge résistive. Une extension à IRn a été proposée et démontrée. Elle est illustrée sur un système non-linéaire dans IR3.La deuxième méthode est développée dans IR2 et basée sur la théorie de Lyapunov, bien connue en automatique pour étudier la stabilité des systèmes non-linéaires et concevoir des commandes stabilisantes.Il s’agit de déterminer par une approche géométrique, une fonction de Lyapunov quadratique commune aux deux modes de fonctionnement du système, qui permette d’obtenir un cycle limite hybride stable le plus proche possible du point de fonctionnement désiré et une commande stabilisante directe des interrupteurs / This PhD-thesis presents two new approaches for the analysis and control of complex nonlinearsystems, such as switching dynamic systems of the class of power converters. These systems have severalmodes of operation and a desired operating point, which, in general, is not the equilibrium point of anymode. In this class of systems, switching from one mode to another is controlled by a switching law tobe designed. Therefore, the synthesis of a control law involves the study of the conditions that allow astable limit cycle to settle near the desired operating point, and of the control trajectory to reach thislimit cycle and stabilize on it, meeting the constraints dues to the physical behavior (maximum currentsupported by the components, . . .) or time constraints (minimum duration between two switchings, . . .).The limit cycle is called hybrid because it is composed of several dynamics (two in this work).The first method is based on the geometric properties of vector fields and is an extension of part ofthe PhD-thesis of Manon at LAGEP. A necessary and sufficient condition of existence and stability of ahybrid limit cycle consisting of a sequence of two operating modes, is presented in IR2. This limit cycledefines the final region to be reached by the system from its initial state, along a trajectory determinedoptimally according to a given criterion (total duration, energy expended, . . .). This method is applied tothe Buck and Buck-Boost power converters with a resistive load. An extension to IRn has been proposedand demonstrated. It is illustrated on a nonlinear system in IR3.The second method is developed in IR2, based on the Lyapunov theory, well known in automatic controlfor studying the stability of nonlinear systems and designing stabilizing control methods. The objective isto design, with a geometric approach, a quadratic Lyapunov function common to both modes of operation,which defines a stable hybrid limit cycle closest to the desired operating point and a direct stabilizingcontrol of the switches.

Le travail présenté dans cette thèse porte sur la modélisation et la commande d'un système thermodynamique non linéaire de dimension infinie, le réacteur tubulaire. Nous abordons le problème de commande sur ce système non linéaire en nous appuyant sur les propriétés thermodynamiques du procédé. Cette approche nécessite l'utilisation d'un modèle ayant comme variables d'état les variables extensives thermodynamiques classiques. Nous utilisons la fonction de disponibilité thermodynamique ainsi qu'une autre fonction déduite de la précédente, la disponibilité réduite, comme fonction de Lyapunov candidate pour résoudre le problème de stabilisation du réacteur autour d'un profil d'équilibre en utilisant comme commande distribuée la température de la double enveloppe. Des simulations illustrent ces résultats ainsi que l'efficacité des commandes en présence de perturbations. Nous nous intéressons aussi à la représentation hamiltonienne à port des systèmes irréversibles de dimension infinie. La structure de Stokes-Dirac pour un modèle réaction diffusion est obtenue en étendant les vecteurs de variables de flux et d'effort. Nous présentons cette démarche pour les équations du système réaction-diffusion en prenant premièrement l'énergie interne comme Hamiltonien puis deuxièmement l'opposé de l'entropie. Nous montrons dans les deux cas qu'en utilisant une extension des couples de variables effort-flux thermodynamiques classiques nous obtenons une structure de Stokes-Dirac. Enfin nous donnons quelques résultats aboutissant à une représentation pseudo hamiltonienne. Enfin nous abordons le problème de commande à la frontière. L'objectif est d'étudier l'existence de solutions associées à un modèle linéarisé de réacteur tubulaire complet commandé à la frontière / The main objective of this thesis consists to investigate the problem of modelling and control of a nonlinear parameter distributed thermodynamic system : the tubular reactor. We address the control problem of this non linear system relying on the thermodynamic properties of the process. This approach requires to use the classical extensive variables as the state variables. We use the thermodynamic availability as well as the reduced thermodynamic availability (this function is formed from some terms of the thermodynamic availabilty) as Lyapunov functions in order to asymptotically stabilize the tubular reactor aroud a steady profile. The distributed temperature of the jacket is the control variable. Some simulations illustrate these results as well as the eficiency of the control in presence of perturbations. Next we study the Port Hamiltonian representation of irreversible infinite dimensional systems. We propose a Stokes-Dirac structure of a reaction-diffusion system by means of the extension of the vectors of the flux and effort variables. We illustrate this approach on the example of the reaction-diffusion system. For this latter we use the internal energy as well as the opposite of the entropy to obtain Stokes-Dirac structures. We propose also a pseudo-Hamiltonian representation for the two Hamiltonians. Finally we tackle the boundary control problem. The objective is to study the existence of solutions associated to a linearized model of the tubular reactor controlled to the boundary