Analysis of energy dissipation in stirred suspension polymerisation reactors using computational fluid dynamics

This work evaluates the spatial distribution of normalised rates of droplet breakage and droplet coalescence in liquidliquid dispersions maintained in agitated tanks at operation conditions normally used to perform suspension polymerisation reactions. Particularly, simulations are performed with multiphase computational fluid dynamics (CFD) models to represent the flow field in liquidliquid styrene suspension polymerisation reactors for the first time. CFD tools are used first to compute the spatial distribution of the turbulent energy dissipation rates (e) inside the reaction vessel; afterwards, normalised rates of droplet breakage and particle coalescence are computed as functions of e. Surprisingly, multiphase simulations showed that the rates of energy dissipation can be very high near the free vortex surfaces, which has been completely neglected in previous works. The obtained results indicate the existence of extremely large energy dissipation gradients inside the vessel, so that particle breakage occurs primarily in very small regions that surround the impeller and the free vortex surface, while particle coalescence takes place in the liquid bulk. As a consequence, particle breakage should be regarded as an independent source term or a boundary phenomenon. Based on the obtained results, it can be very difficult to justify the use of isotropic assumptions to formulate particle population balances in similar systems, even when multiple compartment models are used to describe the fluid dynamic behaviour of the agitated vessel. (C) 2011 Canadian Society for Chemical Engineering

Performance and composition of bacterial communities in anaerobic fluidized bed reactors for hydrogen production: Effects of organic loading rate and alkalinity

This study evaluated the effects of the organic loading rate (OLR) and pH buffer addition on hydrogen production in two anaerobic fluidized bed reactors (AFBRs) operated simultaneously. The AFBRs were fed with glucose, and expanded clay was used as support material. The reactors were operated at a temperature of 30 degrees C, without the addition of a buffer (AFBR1) and with the addition of a pH buffer (AFBR2, sodium bicarbonate) for OLRs ranging from 19.0 to 140.6 kg COD m(-3) d(-1) (COD: chemical oxygen demand). The maximum hydrogen yields for AFBR1 and AFBR2 were 2.45 and 1.90 mol H-2 mol(-1) glucose (OLR of 84.3 kg COD m(-3) d(-1)), respectively. The highest hydrogen production rates were 0.95 and 0.76 L h(-1) L-1 for AFBR1 and AFBR2 (OLR of 140.6 kg COD m(-3) d(-1)), respectively. The operating conditions in AFBR1 favored the presence of such bacteria as Clostridium, while the bacteria in AFBR2 included Clostridium, Enterobacter, Klebsiella, Veillonellaceae, Chryseobacterium, Sporolactobacillus, and Burkholderiaceae. Copyright (C) 2012, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Gauge and integrable theories in loop spaces

We propose an integral formulation of the equations of motion of a large class of field theories which leads in a quite natural and direct way to the construction of conservation laws. The approach is based on generalized non-abelian Stokes theorems for p-form connections, and its appropriate mathematical language is that of loop spaces. The equations of motion are written as the equality of a hyper-volume ordered integral to a hyper-surface ordered integral on the border of that hyper-volume. The approach applies to integrable field theories in (1 + 1) dimensions, Chern-Simons theories in (2 + 1) dimensions, and non-abelian gauge theories in (2 + 1) and (3 + 1) dimensions. The results presented in this paper are relevant for the understanding of global properties of those theories. As a special byproduct we solve a long standing problem in (3 + 1)-dimensional Yang-Mills theory, namely the construction of conserved charges, valid for any solution, which are invariant under arbitrary gauge transformations. (C) 2012 Elsevier B.V. All rights reserved.

Effect of sulfide concentration on autotrophic denitrification from nitrate and nitrite in vertical fixed-bed reactors

Autotrophic denitrification coupled with sulfide oxidation represents an interesting alternative for the simultaneous removal of nitrate/nitrite and sulfide from wastewaters. The applicability of such bioprocess is especially advantageous for the post treatment of effluents from anaerobic reactors, since they usually produce sulfides, which can be used as endogenous electron donor for autotrophic denitrification. This study evaluated the effect of sulfide concentration on this bioprocess using nitrate and nitrite as electron acceptors in vertical fixed-bed reactors. The results showed that intermediary sulfur compounds were mainly produced when excess of electron donor was applied, which was more evident when nitrate was used. Visual evidences suggested that elemental sulfur was the intermediary compound produced. There was also evidence that the elemental sulfur previously formed was being used when sulfide was applied in stoichiometric concentration relative to nitrate/nitrite. Nitrite was more readily consumed than nitrate. For both electron acceptors and sulfide concentrations tested, autotrophic denitrification was not affected by residual heterotrophic denitrification via endogenic activity, occurring as a minor additional nitrogen removal process. (C) 2012 Elsevier Ltd. All rights reserved.

Closed-loop stable model predictive control of integrating systems with dead time

Model predictive control (MPC) applications in the process industry usually deal with process systems that show time delays (dead times) between the system inputs and outputs. Also, in many industrial applications of MPC, integrating outputs resulting from liquid level control or recycle streams need to be considered as controlled outputs. Conventional MPC packages can be applied to time-delay systems but stability of the closed loop system will depend on the tuning parameters of the controller and cannot be guaranteed even in the nominal case. In this work, a state space model based on the analytical step response model is extended to the case of integrating time systems with time delays. This model is applied to the development of two versions of a nominally stable MPC, which is designed to the practical scenario in which one has targets for some of the inputs and/or outputs that may be unreachable and zone control (or interval tracking) for the remaining outputs. The controller is tested through simulation of a multivariable industrial reactor system. (C) 2012 Elsevier Ltd. All rights reserved.

All-loop finiteness of the two-dimensional noncommutative supersymmetric gauge theory

Within the superfield approach, we discuss the two-dimensional noncommutative super-QED. Its all-order finiteness is explicitly shown. Copyright (C) EPLA, 2012

Running coupling and pomeron loop effects on inclusive and diffractive DIS cross sections

Within the framework of a (1 + 1)-dimensional model which mimics high-energy QCD, we study the behavior of the cross sections for inclusive and diffractive deep inelastic gamma*h scattering cross sections. We analyze the cases of both fixed and running coupling within the mean-field approximation, in which the evolution of the scattering amplitude is described by the Balitsky-Kovchegov equation, and also through the pomeron loop equations, which include in the evolution the gluon number fluctuations. In the diffractive case, similarly to the inclusive one, suppression of the diffusive scaling, as a consequence of the inclusion of the running of the coupling, is observed.

Conformation analysis of a surface loop that controls active site access in the GH11 xylanase A from Bacillus subtilis

Xylanases (EC 3.2.1.8 endo-1,4-glycosyl hydrolase) catalyze the hydrolysis of xylan, an abundant hemicellulose of plant cell walls. Access to the catalytic site of GH11 xylanases is regulated by movement of a short beta-hairpin, the so-called thumb region, which can adopt open or closed conformations. A crystallographic study has shown that the D11F/R122D mutant of the GH11 xylanase A from Bacillus subtilis (BsXA) displays a stable "open" conformation, and here we report a molecular dynamics simulation study comparing this mutant with the native enzyme over a range of temperatures. The mutant open conformation was stable at 300 and 328 K, however it showed a transition to the closed state at 338 K. Analysis of dihedral angles identified thumb region residues Y113 and T123 as key hinge points which determine the open-closed transition at 338 K. Although the D11F/R122D mutations result in a reduction in local inter-intramolecular hydrogen bonding, the global energies of the open and closed conformations in the native enzyme are equivalent, suggesting that the two conformations are equally accessible. These results indicate that the thumb region shows a broader degree of energetically permissible conformations which regulate the access to the active site region. The R122D mutation contributes to the stability of the open conformation, but is not essential for thumb dynamics, i.e., the wild type enzyme can also adapt to the open conformation.

Determination of the intrinsic kinetic parameters of sulfide-oxidizing autotrophic denitrification in differential reactors containing immobilized biomass

Nitrogen removal coupled with sulfide oxidation has potential for the treatment of effluents from anaerobic reactors because they contain sulfide, which can be used as an endogenous electron donor for denitrification. This work evaluated the intrinsic kinetics of sulfide-oxidizing autotrophic denitrification via nitrate and nitrite in systems containing attached cells. Differential reactors were fed with nitrified synthetic domestic sewage and different sulfide concentrations. The intrinsic kinetic parameters of nitrogen removal were determined when the mass transfer resistance was negligible. This bioprocess could be described by a half-order kinetic model for biofilms. The half-order kinetic coefficients ranged from 0.425 to 0.658 mg N-1/2 L-1/2 h(-1) for denitrification via nitrite and from 0.190 to 0.609 mg N-1/2 L-1/2 h(-1) for denitrification via nitrate. In this latter, the lower value was due to the use of electrons donated from intermediary sulfur compounds whose formation and subsequent consumption were detected. (C) 2011 Elsevier Ltd. All rights reserved.

Long-term variability and effects on larval fish distribution of the Gulf of Mexico Loop Current and Rings

Programa de doctorado: Oceanografía

Sviluppo di modelli motore e veicolo per l'analisi di strategie di controllo in applicazioni Software e Hardware In the Loop

This work describes the development of a simulation tool which allows the simulation of the Internal Combustion Engine (ICE), the transmission and the vehicle dynamics. It is a control oriented simulation tool, designed in order to perform both off-line (Software In the Loop) and on-line (Hardware In the Loop) simulation. In the first case the simulation tool can be used in order to optimize Engine Control Unit strategies (as far as regard, for example, the fuel consumption or the performance of the engine), while in the second case it can be used in order to test the control system. In recent years the use of HIL simulations has proved to be very useful in developing and testing of control systems. Hardware In the Loop simulation is a technology where the actual vehicles, engines or other components are replaced by a real time simulation, based on a mathematical model and running in a real time processor. The processor reads ECU (Engine Control Unit) output signals which would normally feed the actuators and, by using mathematical models, provides the signals which would be produced by the actual sensors. The simulation tool, fully designed within Simulink, includes the possibility to simulate the only engine, the transmission and vehicle dynamics and the engine along with the vehicle and transmission dynamics, allowing in this case to evaluate the performance and the operating conditions of the Internal Combustion Engine, once it is installed on a given vehicle. Furthermore the simulation tool includes different level of complexity, since it is possible to use, for example, either a zero-dimensional or a one-dimensional model of the intake system (in this case only for off-line application, because of the higher computational effort). Given these preliminary remarks, an important goal of this work is the development of a simulation environment that can be easily adapted to different engine types (single- or multi-cylinder, four-stroke or two-stroke, diesel or gasoline) and transmission architecture without reprogramming. Also, the same simulation tool can be rapidly configured both for off-line and real-time application. The Matlab-Simulink environment has been adopted to achieve such objectives, since its graphical programming interface allows building flexible and reconfigurable models, and real-time simulation is possible with standard, off-the-shelf software and hardware platforms (such as dSPACE systems).

The massless two-loop two-point function and zeta functions in counterterms of Feynman diagrams

The main part of this thesis describes a method of calculating the massless two-loop two-point function which allows expanding the integral up to an arbitrary order in the dimensional regularization parameter epsilon by rewriting it as a double Mellin-Barnes integral. Closing the contour and collecting the residues then transforms this integral into a form that enables us to utilize S. Weinzierl's computer library nestedsums. We could show that multiple zeta values and rational numbers are sufficient for expanding the massless two-loop two-point function to all orders in epsilon. We then use the Hopf algebra of Feynman diagrams and its antipode, to investigate the appearance of Riemann's zeta function in counterterms of Feynman diagrams in massless Yukawa theory and massless QED. The class of Feynman diagrams we consider consists of graphs built from primitive one-loop diagrams and the non-planar vertex correction, where the vertex corrections only depend on one external momentum. We showed the absence of powers of pi in the counterterms of the non-planar vertex correction and diagrams built by shuffling it with the one-loop vertex correction. We also found the invariance of some coefficients of zeta functions under a change of momentum flow through these vertex corrections.

O (alpha 4 s) loop-by-loop contributions to heavy quark pair production in hadronic collisions

The present state of the theoretical predictions for the hadronic heavy hadron production is not quite satisfactory. The full next-to-leading order (NLO) ${cal O} (alpha_s^3)$ corrections to the hadroproduction of heavy quarks have raised the leading order (LO) ${cal O} (alpha_s^2)$ estimates but the NLO predictions are still slightly below the experimental numbers. Moreover, the theoretical NLO predictions suffer from the usual large uncertainty resulting from the freedom in the choice of renormalization and factorization scales of perturbative QCD.In this light there are hopes that a next-to-next-to-leading order (NNLO) ${cal O} (alpha_s^4)$ calculation will bring theoretical predictions even closer to the experimental data. Also, the dependence on the factorization and renormalization scales of the physical process is expected to be greatly reduced at NNLO. This would reduce the theoretical uncertainty and therefore make the comparison between theory and experiment much more significant. In this thesis I have concentrated on that part of NNLO corrections for hadronic heavy quark production where one-loop integrals contribute in the form of a loop-by-loop product. In the first part of the thesis I use dimensional regularization to calculate the ${cal O}(ep^2)$ expansion of scalar one-loop one-, two-, three- and four-point integrals. The Laurent series of the scalar integrals is needed as an input for the calculation of the one-loop matrix elements for the loop-by-loop contributions. Since each factor of the loop-by-loop product has negative powers of the dimensional regularization parameter $ep$ up to ${cal O}(ep^{-2})$, the Laurent series of the scalar integrals has to be calculated up to ${cal O}(ep^2)$. The negative powers of $ep$ are a consequence of ultraviolet and infrared/collinear (or mass ) divergences. Among the scalar integrals the four-point integrals are the most complicated. The ${cal O}(ep^2)$ expansion of the three- and four-point integrals contains in general classical polylogarithms up to ${rm Li}_4$ and $L$-functions related to multiple polylogarithms of maximal weight and depth four. All results for the scalar integrals are also available in electronic form. In the second part of the thesis I discuss the properties of the classical polylogarithms. I present the algorithms which allow one to reduce the number of the polylogarithms in an expression. I derive identities for the $L$-functions which have been intensively used in order to reduce the length of the final results for the scalar integrals. I also discuss the properties of multiple polylogarithms. I derive identities to express the $L$-functions in terms of multiple polylogarithms. In the third part I investigate the numerical efficiency of the results for the scalar integrals. The dependence of the evaluation time on the relative error is discussed. In the forth part of the thesis I present the larger part of the ${cal O}(ep^2)$ results on one-loop matrix elements in heavy flavor hadroproduction containing the full spin information. The ${cal O}(ep^2)$ terms arise as a combination of the ${cal O}(ep^2)$ results for the scalar integrals, the spin algebra and the Passarino-Veltman decomposition. The one-loop matrix elements will be needed as input in the determination of the loop-by-loop part of NNLO for the hadronic heavy flavor production.