Integration of different models in the design of chemical processes: Application to the design of a power plant

With advances in the synthesis and design of chemical processes there is an increasing need for more complex mathematical models with which to screen the alternatives that constitute accurate and reliable process models. Despite the wide availability of sophisticated tools for simulation, optimization and synthesis of chemical processes, the user is frequently interested in using the ‘best available model’. However, in practice, these models are usually little more than a black box with a rigid input–output structure. In this paper we propose to tackle all these models using generalized disjunctive programming to capture the numerical characteristics of each model (in equation form, modular, noisy, etc.) and to deal with each of them according to their individual characteristics. The result is a hybrid modular–equation based approach that allows synthesizing complex processes using different models in a robust and reliable way. The capabilities of the proposed approach are discussed with a case study: the design of a utility system power plant that has been decomposed into its constitutive elements, each treated differently numerically. And finally, numerical results and conclusions are presented.

Removal of paracetamol on biomass-derived activated carbon: Modeling the fixed bed breakthrough curves using batch adsorption experiments

The remediation of paracetamol (PA), an emerging contaminant frequently found in wastewater treatment plants, has been studied in the low concentration range (0.3–10 mg L−1) using as adsorbent a biomass-derived activated carbon. PA uptake of up to 100 mg g−1 over the activated carbon has been obtained, with the adsorption isotherms being fairly explained by the Langmuir model. The application of Reichemberg and the Vermeulen equations to the batch kinetics experiments allowed estimating homogeneous and heterogeneous diffusion coefficients, reflecting the dependence of diffusion with the surface coverage of PA. A series of rapid small-scale column tests were carried out to determine the breakthrough curves under different operational conditions (temperature, PA concentration, flow rate, bed length). The suitability of the proposed adsorbent for the remediation of PA in fixed-bed adsorption was proven by the high PA adsorption capacity along with the fast adsorption and the reduced height of the mass transfer zone of the columns. We have demonstrated that, thanks to the use of the heterogeneous diffusion coefficient, the proposed mathematical approach for the numerical solution to the mass balance of the column provides a reliable description of the breakthrough profiles and the design parameters, being much more accurate than models based in the classical linear driving force.

The evolution of mass tourism destinations: new approaches beyond deterministic models in Benidorm (Spain)

Although deterministic models of the evolution of mass tourism coastal resorts predict an almost inevitable decline over time, theoretical frameworks of the evolution and restructuring policies of mature destinations should be revised to reflect the complex and dynamic way in which these destinations evolve and interact with the tourism market and global socio-economic environment. The present study examines Benidorm because its urban and tourism model and large-scale tourism supply and demand make it one of the most unique destinations on the Mediterranean coast. The investigation reveals the need to adopt theories and models that are not purely deterministic. The dialectic interplay between external factors and the internal factors inherent in this destination simultaneously reveals a complex and diverse stage of maturity and the ability of destinations to create their own future.

Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

Effect of the porous structure in carbon materials for CO2 capture at atmospheric and high-pressure

Activated carbons prepared from petroleum pitch and using KOH as activating agent exhibit an excellent behavior in CO2 capture both at atmospheric (∼168 mg CO2/g at 298 K) and high pressure (∼1500 mg CO2/g at 298 K and 4.5 MPa). However, an exhaustive evaluation of the adsorption process shows that the optimum carbon structure, in terms of adsorption capacity, depends on the final application. Whereas narrow micropores (pores below 0.6 nm) govern the sorption behavior at 0.1 MPa, large micropores/small mesopores (pores below 2.0–3.0 nm) govern the sorption behavior at high pressure (4.5 MPa). Consequently, an optimum sorbent exhibiting a high working capacity for high pressure applications, e.g., pressure-swing adsorption units, will require a poorly-developed narrow microporous structure together with a highly-developed wide microporous and small mesoporous network. The appropriate design of the preparation conditions gives rise to carbon materials with an extremely high delivery capacity ∼1388 mg CO2/g between 4.5 MPa and 0.1 MPa. Consequently, this study provides guidelines for the design of carbon materials with an improved ability to remove carbon dioxide from the environment at atmospheric and high pressure.

KOH activation of carbon materials obtained from the pyrolysis of ethylene tar at different temperatures

A complete study of the importance of the pyrolysis temperature (up to 1500 °C) of a petroleum residue (ethylene tar) in the activation with KOH of the resultant pyrolysis products (covering from the own ethylene tar to pitches and well developed cokes) has been carried out. The trend in the porosity found for activated carbons is as follows: the pore volume increases with the pyrolysis temperature reaching a maximum value (1.39 cm3/g) at about 460 °C, just at the transition temperature between a fluid pitch and a solid coke. It is the pitch with highest mesophase content that develops the maximum porosity when activated with KOH. The amount of H2, CO and CO2 produced during the reaction of the mesophase pitch and coke with KOH has been quantified, and a trend as described for the pore volume was found with the pyrolysis temperature. Therefore, there is a relationship between the reactivity of the precursor with KOH and the porosity developed by the activated carbon. Since the reactions that produce H2 initiate at temperatures as low as 300 °C, it seems that KOH is modifying the conditions under which the pyrolysis occurs, and this fact is critical in the development of porosity.

CeO2-promoted Ni/activated carbon catalysts for the water–gas shift (WGS) reaction

The low temperature water–gas shift (WGS) reaction has been studied over carbon-supported nickel catalysts promoted by ceria. To this end, cerium oxide has been dispersed (at different loadings: 10, 20, 30 and 40 wt.%) on the activated carbon surface with the aim of obtaining small ceria particles and a highly available surface area. Furthermore, carbon- and ceria-supported nickel catalysts have also been studied as references. A combination of N2 adsorption analysis, powder X-ray diffraction, temperature-programmed reduction with H2, X-ray photoelectron spectroscopy and TEM analysis were used to characterize the Ni–CeO2 interactions and the CeO2 dispersion over the activated carbon support. Catalysts were tested in the low temperature WGS reaction with two different feed gas mixtures: the idealized one (with only CO and H2O) and a slightly harder one (with CO, CO2, H2, and H2O). The obtained results show that there is a clear effect of the ceria loading on the catalytic activity. In both cases, catalysts with 20 and 10 wt.% CeO2 were the most active materials at low temperature. On the other hand, Ni/C shows a lower activity, this assessing the determinant role of ceria in this reaction. Methane, a product of side reactions, was observed in very low amounts, when CO2 and H2 were included in the WGS feed. Nevertheless, our data indicate that the methanation process is mainly due to CO2, and no CO consumption via methanation takes place at the relevant WGS temperatures. Finally, a stability test was carried out, obtaining CO conversions greater than 40% after 150 h of reaction.

Functionalization of carbon nanotubes using aminobenzene acids and electrochemical methods. Electroactivity for the oxygen reduction reaction

Functionalized carbon nanotubes (CNTs) using three aminobenzene acids with different functional groups (carboxylic, sulphonic, phosphonic) in para position have been synthesized through potentiodynamic treatment in acid media under oxidative conditions. A noticeable increase in the capacitance for the functionalized carbon nanotubes mainly due to redox processes points out the formation of an electroactive polymer thin film on the CNTs surface along with covalently bonded functionalities. The CNTs functionalized using aminobenzoic acid rendered the highest capacitance values and surface nitrogen content, while the presence of sulfur and/or phosphorus groups in the aminobenzene structure yielded a lower functionalization degree. The oxygen reduction reaction (ORR) activity of the functionalized samples was similar to that of the parent CNTs, independently of the functional group present in the aminobenzene acid. Interestingly, a heat treatment in N2 atmosphere with a very low O2 concentration (3125 ppm) at 800 °C of the CNTs functionalized with aminobenzoic acid produced a material with high amounts of surface oxygen and nitrogen groups (12 and 4% at., respectively), that seem to modulate the electron-donor properties of the resulting material. The onset potential and limiting current for ORR was enhanced for this material. These are promising results that validates the use of electrochemistry for the synthesis of novel N-doped electrocatalysts for ORR in combination with adequate heat treatments.

A multi-method study of the transformation of the carbonaceous skeleton of a polymer-based nanoporous carbon along the activation pathway

The change in the carbonaceous skeleton of nanoporous carbons during their activation has received limited attention, unlike its counterpart process in the presence of an inert atmosphere. Here we adopt a multi-method approach to elucidate this change in a poly(furfuryl alcohol)-derived carbon activated using cyclic application of oxygen saturation at 250 °C before its removal (with carbon) at 800 °C in argon. The methods used include helium pycnometry, synchrotron-based X-ray diffraction (XRD) and associated radial distribution function (RDF) analysis, transmission electron microscopy (TEM) and, uniquely, electron energy-loss spectroscopy spectrum-imaging (EELS-SI), electron nanodiffraction and fluctuation electron microscopy (FEM). Helium pycnometry indicates the solid skeleton of the carbon densifies during activation from 78% to 93% of graphite. RDF analysis, EELS-SI, and FEM all suggest this densification comes through an in-plane growth of sp2 carbon out to the medium range without commensurate increase in order normal to the plane. This process could be termed ‘graphenization’. The exact way in which this process occurs is not clear, but TEM images of the carbon before and after activation suggest it may come through removal of the more reactive carbon, breaking constraining cross-links and creating space that allows the remaining carbon material to migrate in an annealing-like process.

Platinum supported on highly-dispersed ceria on activated carbon for the total oxidation of VOCs

Catalysts consisting in platinum supported on cerium oxide highly dispersed on activated carbon, with a Pt loading of 1 wt.% and ceria loadings of 5, 10 and 20 wt.% have been prepared by impregnation method and characterized by several techniques (N2 adsorption at 77 K, ICP, XRD, H2-TPR and XPS). Their catalytic behavior has been evaluated in the total oxidation of ethanol and toluene after reduction at 473 K. The obtained results show that the prepared catalysts have better performances than platinum supported on bulk CeO2. The best catalytic performance was obtained for 10 wt.% ceria loading, likely due to an optimum synergistic interaction between highly dispersed cerium oxide and platinum particles. Pt-10Ce/C achieves total conversion of ethanol and toluene to CO2 at 433 K and 453 K, respectively, and shows no deactivation during a test for 100 h. Under humid conditions (relative humidity, RH, of 40 and 80%), the activity was only slightly influenced due to the hydrophobic character of the activated carbon support, which prevents the adsorption of water.

Hydrothermal carbonization of lignocellulosic biomass: Effect of process conditions on hydrochar properties

Although hydrothermal carbonization of biomass components is known to be mainly governed by reaction temperature, consistent reports on the effect and statistical significance of process conditions on hydrochar properties are still lacking. The objective of this research was to determine the importance and significance of reaction temperature, retention time and solid load on the properties of hydrochar produced from an industrial lignocellulosic sludge residue. According to the results, reaction temperature and retention time had a statistically significant effect on hydrochar ash content, solid yield, carbon content, O/C-ratio, energy densification and energy yield as reactor solid load was statistically insignificant for all acquired models within the design range. Although statistically significant, the effect of retention time was 3–7 times lower than that of reaction temperature. Predicted dry ash-free solid yields of attained hydrochar decreased to approximately 40% due to the dissolution of biomass components at higher reaction temperatures, as respective oxygen contents were comparable to subbituminous coal. Significant increases in the carbon contents of hydrochar led to predicted energy densification ratios of 1–1.5 with respective energy yields of 60–100%. Estimated theoretical energy requirements of carbonization were dependent on the literature method used and mainly controlled by reaction temperature and reactor solid load. The attained results enable future prediction of hydrochar properties from this feedstock and help to understand the effect of process conditions on hydrothermal treatment of lignocellulosic biomass.