Vapour reactive distillation process for hydrogen production by hi decomposition from hi-i2-h2o solutions

In this contribution, a sequential and hierarchical approach for the feasibility analysis and the preliminary design of reactive distillation columns is extended to systems involving vapour phase chemical reaction and is successfully applied to the HI vapour phase decomposition to produce H2. The complex phase and physico chemical behaviour of the quaternary HI-H2-I2-H2O system is represented by the Neumann’s thermodynamic model and instantaneous vapour phase chemical equilibrium is assumed. Then, from minimal information concerning the physicochemical properties of the system, three successive steps lead to the design of the unit and the specification of its operating conditions: the feasibility analysis, the synthesis and the design step. First, the analysis of reactive condensation curve map method (rCCM), assuming infinite internal liquid and vapour flow rate and infinite reflux ratio, is used to assess the feasibility of the process. It determines the column structure and estimates the attainable compositions. These results are used as inputs data for the synthesis step. Based on the boundary value design method (BVD), considering finite internal liquid and vapour flow rate and finite reflux ratio while neglecting all thermal effects and assuming a constant heat of vaporisation, the synthesis step provides more precise information about the process configuration (minimum reflux ratio, number of theoretical stages, localisation and number of reactive plates, position of the feed plate). Finally, the BVD method results are used to initialise rigorous simulations, based on an equilibrium stage model with energy balance, to estimate the reflux ratio taking into account thermal effect on the process. The resulting design configuration consists in a single feed and entirely reactive distillation column. The column operates under a pressure of 22 bars. The feed of the reactive distillation column, coming from the Bunsen reaction section [xHI=0.10; xI2=0.39 xH2O=0.51], is at its boiling temperature. The residue consists in pure iodine. Water and produced hydrogen are recovered at the distillate. The column operates at a reflux ratio of 5 and is composed of 11 theoretical plates including the reboiler and the partial condenser with the feed at the stage 10 (counted downwards). The obtained HI dissociation yield is 99.6%

Simulation of inviscid compressible multi-phase flow with condensation

Condensation of vapours in rapid expansions of compressible gases is investigated. In the case of high temperature gradients the condensation will start at conditions well away from thermodynamic equilibrium of the fluid. In those cases homogeneous condensation is dominant over heterogeneous condensation. The present work is concerned with development of a simulation tool for computation of high speed compressible flows with homogeneous condensation. The resulting flow solver should preferably be accurate and robust to be used for simulation of industrial flows in general geometries

A numerical model for the fractional condensation of pyrolysis vapours

Experimentation on the fast pyrolysis process has been primarily focused on the pyrolysis reactor itself, with less emphasis given to the liquid collection system (LCS). More importantly, the physics behind the vapour condensation process in LCSs has not been thoroughly researched mainly due to the complexity of the phenomena involved. The present work focusses on providing detailed information of the condensation process within the LCS, which consists of a water cooled indirect contact condenser. In an effort to understand the mass transfer phenomena within the LCS, a numerical simulation was performed using the Eulerian approach. A multiphase multi-component model, with the condensable vapours and non-condensable gases as the gaseous phase and the condensed bio-oil as the liquid phase, has been created. Species transport modelling has been used to capture the detailed physical phenomena of 11 major compounds present in the pyrolysis vapours. The development of the condensation model relies on the saturation pressures of the individual compounds based on the corresponding states correlations and assuming that the pyrolysis vapours form an ideal mixture. After the numerical analysis, results showed that different species condense at different times and at different rates. In this simulation, acidic components like acetic acid and formic acids were not condensed as it was also evident in experimental works, were the pH value of the condensed oil is higher than subsequent stages. In the future, the current computational model can provide significant aid in the design and optimization of different types of LCSs

Theory of Sorption Hysteresis in Nanoporous Solids: II. Molecular condensation

Motivated by the puzzle of sorption hysteresis in Portland cement concrete or cement paste, we develop in Part II of this study a general theory of vapor sorption and desorption from nanoporous solids, which attributes hysteresis to hindered molecular condensation with attractive lateral interactions. The classical mean-field theory of van der Waals is applied to predict the dependence of hysteresis on temperature and pore size, using the regular solution model and gradient energy of Cahn and Hilliard. A simple "hierarchical wetting" model for thin nanopores is developed to describe the case of strong wetting by the first monolayer, followed by condensation of nanodroplets and nanobubbles in the bulk. The model predicts a larger hysteresis critical temperature and enhanced hysteresis for molecular condensation across nanopores at high vapor pressure than within monolayers at low vapor pressure. For heterogeneous pores, the theory predicts sorption/desorption sequences similar to those seen in molecular dynamics simulations, where the interfacial energy (or gradient penalty) at nanopore junctions acts as a free energy barrier for snap-through instabilities. The model helps to quantitatively understand recent experimental data for concrete or cement paste wetting and drying cycles and suggests new experiments at different temperatures and humidity sweep rates.Comment: 26 pages, 10 fig

On the Modeling of Droplet Evaporation on Superhydrophobic Surfaces

When a drop of water is placed on a rough surface, there are two possible extreme regimes of wetting: the one called Cassie-Baxter (CB) with air pockets trapped underneath the droplet and the one characterized by the homogeneous wetting of the surface, called the Wenzel (W) state. A way to investigate the transition between these two states is by means of evaporation experiments, in which the droplet starts in a CB state and, as its volume decreases, penetrates the surface's grooves, reaching a W state. Here we present a theoretical model based on the global interfacial energies for CB and W states that allows us to predict the thermodynamic wetting state of the droplet for a given volume and surface texture. We first analyze the influence of the surface geometric parameters on the droplet's final wetting state with constant volume, and show that it depends strongly on the surface texture. We then vary the volume of the droplet keeping fixed the geometric surface parameters to mimic evaporation and show that the drop experiences a transition from the CB to the W state when its volume reduces, as observed in experiments. To investigate the dependency of the wetting state on the initial state of the droplet, we implement a cellular Potts model in three dimensions. Simulations show a very good agreement with theory when the initial state is W, but it disagrees when the droplet is initialized in a CB state, in accordance with previous observations which show that the CB state is metastable in many cases. Both simulations and theoretical model can be modified to study other types of surface.Comment: 23 pages, 7 figure

Contrail microphysics in the near wake of a realistic wing through RANS simulations

This paper focuses on Steady Reynolds Average Navier Stokes simulations (RANS) of ice particles growth in the near field of a wing-injector configuration. The multiphysics multiphase flow solver CEDRE, enriched with a microphysical model, has been developed in order to simulate the impact of a more real aircraft geometry in contrail formation studies. As a first evaluation case, a simplified aircraft description, i.e. a NACA0012 2D wing with two injectors, has been used. Ice formation has been simulated by assuming water condensation and instantaneous freezing on activated soot particles, initially emitted by aircraft engines. Our investigation focuses on the near field, extending from the nozzle exit until eight wing spans. Although the main goal is to address the question of ice formation, the aerodynamic flow field has been investigated and numerical results compared with existing experimental data. The first results indicate that the exhaust jet is correctly wrapped around the vortex and that the pattern of dilution qualitatively matches observations in the near field. Sensitivity studies to humidity and to the initial soot particle radius have also been performed