Some examples of modeling in chemical engineering: thermal treatment of polymers, phase equilibrium calculations and process design

Presentation submitted to PSE Seminar, Chemical Engineering Department, Center for Advanced Process Design-making (CAPD), Carnegie Mellon University, Pittsburgh (USA), October 2012.

Clay-supported graphene materials: application to hydrogen storage

The present work refers to clay–graphene nanomaterials prepared by a green way using caramel from sucrose and two types of natural clays (montmorillonite and sepiolite) as precursors, with the aim of evaluating their potential use in hydrogen storage. The impregnation of the clay substrates by caramel in aqueous media, followed by a thermal treatment in the absence of oxygen of these clay–caramel intermediates gives rise to graphene-like materials, which remain strongly bound to the silicate support. The nature of the resulting materials was characterized by different techniques such as XRD, Raman spectroscopy and TEM, as well as by adsorption isotherms of N2, CO2 and H2O. These carbon–clay nanocomposites can act as adsorbents for hydrogen storage, achieving, at 298 K and 20 MPa, over 0.1 wt% of hydrogen adsorption excess related to the total mass of the system, and a maximum value close to 0.4 wt% of hydrogen specifically related to the carbon mass. The very high isosteric heat for hydrogen sorption determined from adsorption isotherms at different temperatures (14.5 kJ mol−1) fits well with the theoretical values available for hydrogen storage on materials that show a strong stabilization of the H2 molecule upon adsorption.

MINLP Model for the Synthesis of Heat Exchanger Networks with Handling Pressure of Process Streams

This paper introduces a new mathematical model for the simultaneous synthesis of heat exchanger networks (HENs), wherein the handling pressure of process streams is used to enhance the heat integration. The proposed approach combines generalized disjunctive programming (GDP) and mixed-integer nonlinear programming (MINLP) formulation, in order to minimize the total annualized cost composed by operational and capital expenses. A multi-stage superstructure is developed for the HEN synthesis, assuming constant heat capacity flow rates and isothermal mixing, and allowing for streams splits. In this model, the pressure and temperature of streams must be treated as optimization variables, increasing further the complexity and difficulty to solve the problem. In addition, the model allows for coupling of compressors and turbines to save energy. A case study is performed to verify the accuracy of the proposed model. In this example, the optimal integration between the heat and work decreases the need for thermal utilities in the HEN design. As a result, the total annualized cost is also reduced due to the decrease in the operational expenses related to the heating and cooling of the streams.

Chemical Reactors. Problem 2: Constant pressure adiabatic stirred batch reactor with variable heat capacities (2011-2012 Course)

Wording of problem 2 (week 3, 17/10/11).

Chemical Reactors. Problem 2: Constant pressure adiabatic stirred batch reactor with variable heat capacities. Solution (2011-2012 Course)

Chemical Reaction Engineering. Course 2011-12. Solution of problem 2: constant pressure adiabatic stirred batch reactor with variable heat capacities.

Review and extension of the McCabe-Thiele method covering multiple feeds, products and heat transfer stages

The McCabe-Thiele method is a classical approximate graphical method for the conceptual design of binary distillation columns which is still widely used, mainly for didactical purposes, though it is also valuable for quick preliminary calculations. Nevertheless, no complete description of the method has been found and situations such as different thermal feed conditions, multiple feeds, possibilities to extract by-products or to add or remove heat, are not always considered. In the present work we provide a systematic analysis of such situations by developing the generalized equations for: a) the operating lines (OL) of each sector, and b) the changeover line that provides the connection between two consecutive trays of the corresponding sectors separated by a lateral stream of feed, product, or a heat removal or addition.

Simultaneous synthesis of work exchange networks with heat integration

The optimal integration of work and its interaction with heat can represent large energy savings in industrial plants. This paper introduces a new optimization model for the simultaneous synthesis of work exchange networks (WENs), with heat integration for the optimal pressure recovery of process gaseous streams. The proposed approach for the WEN synthesis is analogous to the well-known problem of synthesis of heat exchanger networks (HENs). Thus, there is work exchange between high-pressure (HP) and low-pressure (LP) streams, achieved by pressure manipulation equipment running on common axes. The model allows the use of several units of single-shaft-turbine-compressor (SSTC), as well as stand-alone compressors, turbines and valves. Helper motors and generators are used to respond to any demand and excess of energy. Moreover, between the WEN stages the streams are sent to the HEN to promote thermal recovery, aiming to enhance the work integration. A multi-stage superstructure is proposed to represent the process. The WEN superstructure is optimized in a mixed-integer nonlinear programming (MINLP) formulation and solved with the GAMS software, with the goal of minimizing the total annualized cost. Three examples are conducted to verify the accuracy of the proposed method. In all case studies, the heat integration between WEN stages is essential to improve the pressure recovery, and to reduce the total costs involved in the process.

Analysis of the Connecting Zone between Consecutive Sections in Distillation Columns Covering Multiple Feeds, Products and Heat Transfer Stages

In the present work, we provide a systematic analysis about all tine streams involved in the zone connecting two consecutive sections for the design of distillation columns with different thermal feed conditions, product extractions and heat additions or withdrawals. This analysis allows a better understanding of what happens on a feed or side draw (of mass or energy) stage, what compositions are or are not in equilibrium, and the impact on internal liquid and vapor flows.

Quiescent thermal emission from neutron stars in low-mass X-ray binaries

Context. We monitored the quiescent thermal emission from neutron stars in low-mass X-ray binaries after active periods of intense activity in X-rays (outbursts). Aims. The theoretical modeling of the thermal relaxation of the neutron star crust may be used to establish constraints on the crust composition and transport properties, depending on the astrophysical scenarios assumed. Methods. We numerically simulated the thermal evolution of the neutron star crust and compared them with inferred surface temperatures for five sources: MXB 1659−29, KS 1731−260, XTE J1701−462, EXO 0748−676  and IGR J17480−2446. Results. We find that the evolution of MXB 1659−29, KS 1731−260 and EXO 0748−676 can be well described within a deep crustal cooling scenario. Conversely, we find that the other two sources can only be explained with models beyond crustal cooling. For the peculiar emission of XTE J1701−462 we propose alternative scenarios such as residual accretion during quiescence, additional heat sources in the outer crust, and/or thermal isolation of the inner crust due to a buried magnetic field. We also explain the very recent reported temperature of IGR J17480−2446 with an additional heat deposition in the outer crust from shallow sources.

Thermal Runaway in Tubular Reactor Systems with n-th Order Kinetics

This paper addresses the problem of predicting the critical parameters that characterize thermal runaway in a tubular reactor with wall cooling, introducing a new view of the n-th order kinetics reactions. The paper describes the trajectories of the system in the temperature-(concentration)n plane, and deduces the conditions for the thermal risk.