Robust optimal operation of continuous catalytic reforming process under feedstock uncertainty

The continuous catalytic regenerative (CCR) reforming process is one of the most significant sources of hydrogen production in the petroleum refining process. However, the fluctuations in feedstock composition and flow rate could significantly affect both product distribution and energy consumption. In this study, a robust deviation criterion based multi-objective optimization approach is proposed to perform the optimal operation of CCR reformer under feedstock uncertainty, with simultaneous maximization of product yields and minimization of energy consumption. Minimax approach is adopted to handle these uncertain objectives, and the Latin hypercube sampling method is then used to calculate these robust deviation criteria. Multi-objective surrogate-based optimization methods are next introduced to effectively solve the robust operational problem with high computational cost. The level diagram method is finally utilized to assist in multi-criteria decision-making. Two robust operational optimization problems with different objectives are solved to demonstrate the effectiveness of the proposed method for robust optimal operation of the CCR reforming process under feedstock uncertainty.     

NNI-1型C_5/C_6烷烃异构化催化剂长周期运行分析及建议

中国石化海南炼油化工有限公司0.2 Mt/a C5/C6烷烃异构化装置以连续重整装置的拔头油为原料,使用NNI-1催化剂,采用一次通过流程,不设脱异戊烷塔和稳定塔,经设在连续重整装置内的脱丁烷塔稳定处理后作为汽油调合组分。该装置于2006年9月开工投产,截至2015年3月已连续运行3个周期。长周期运行分析结果表明:前两个周期中NNI-1催化剂具有较高的异构化活性及选择性,C5异构化率为60%左右,C6异构化率为80%左右,C6选择性为15%左右,产品辛烷值基本达到技术指标要求(RON≥78);而在第三周期运行中,催化剂积炭增加等原因导致其异构化活性及选择性降低,异构化产品辛烷值提升能力呈现逐步衰减的趋势,提高反应苛刻度已不能弥补催化剂活性下降造成的产品辛烷值降低。为保证装置长周期运行,建议择机停工对催化剂进行再生,或是直接换用与装置原料性质匹配的异构化催化剂。     

Evaluation of hydrodynamic behavior of the perforated gas distributor of industrial gas phase polymerization reactor using CFD-PBM coupled model

A 2D computational fluid dynamics (Eulerian–Eulerian) multiphase flow model coupled with a population balance model (CFD-PBM) was implemented to investigate the fluidization structure in terms of entrance region in an industrial-scale gas phase fluidized bed reactor. The simulation results were compared with the industrial data, and good agreement was observed. Two cases including perforated distributor and complete sparger were applied to examine the flow structure through the bed. The parametric sensitivity analysis of time step, number of node, drag coefficient, and specularity coefficient was carried out. It was found that the results were more sensitive to the drag model. The results showed that the entrance configuration has significant effect on the flow structure. While the dead zones are created in both corners of the distributors, the perforated distributor generates more startup bubbles, heterogeneous flow field, and better gas–solid interaction above the entrance region due to jet formation.     

Numerical study of entropy generation in Darcy-Forchheimer (D-F) Bödewadt flow of CNTs

Three-dimensional Bödewadt flow (fluid rotates at a large enough distance from the stationary plate) of carbon nanomaterial is examined. Single walled and multi walled CNTs are dissolved in water and gasoline oil baseliquids. Darcy-Forchheimer porous medium is considered. Stationary disk is further stretched linearly in radial direction. Heat transfer effect is examined in presence of radiation and convection. Effect of viscous dissipation is accounted. Entropy generation rate is studied. By using adequate transformation (von Kármán relations), the flow field equations (PDEs) are transmitted into ODEs. Solutions to these ODEs are constructed via implementation of shooting method (bvp4c). In addition to Entropy generation rate, Bejan number, heat transfer rate (Nusselt number), skin friction and temperature of fluid are examined through involved physical parameters. Axial component of velocity intensifies with increment in nanoparticles volume fraction and ratio of stretching rate to angular velocity parameter while it decays with higher porosity parameter. Higher nanoparticles volume fraction and porosity parameter lead to decay in radial as well as tangential component of velocity. However it enhances with higher ratio of stretching rate to angular velocity parameter. Temperature of fluid directly varies with higher ratio of stretching rate to angular velocity parameter, radiation parameter, Eckert number, Biot number and nanoparticles volume fraction. Rate of Entropy generation is reduced with higher estimations of porosity parameter, nanoparticles volume fraction and radiation parameter. Skin friction coefficient decays with higher porosity parameter and ratio of stretching rate to angular velocity parameter. Intensification in porosity parameter, nanoparticles volume fraction and Biot number leads to higher Nusselt number. Prominent impact is shown by multiple-walled CNTs with gasoline oil basefluid than single-walled CNTs with water basefluid.     

Upscaled DEM-CFD model for vibrated fluidized bed based on particle-scale similarities

Simulation based on discrete element method (DEM) coupled with computational fluid dynamics (CFD), coupled DEM-CFD, is a powerful tool for investigating the details of dense particle–fluid interaction problems such as in fluidized beds and pneumatic conveyers. The addition of a mechanical vibration to a system can drastically alter the particle and fluid flows; however, their detailed mechanisms are not well understood. In this study, a DEM-CFD model based on a non-inertial frame of reference is developed to achieve a better understanding of the influence of vibration in a vibrated fluidized bed. Because the high computational cost of DEM-CFD calculations is still a major problem, an upscaled coarse-graining model is also employed. To realize similar behaviors with enlarged model particles, non-dimensional parameters at the particle scale were deduced from the governing equations. The suitability and limitations of the proposed model were examined for a density segregation problem of a binary system. To reduce the computational costs, we show that the ratio between the bed width and model particle size can be reduced to a minimum value of 100; to obtain similar segregation behaviors, the ratio between the bed height and model particle size is considered unchanged.     

Studying local liquid velocity in liquid–solid packed bed using the newly developed X-ray DIR technique

Combination of X-ray Digital Industrial Radiography (DIR) and Particle Tracking Velocimetry (PTV) techniques for local liquid velocity measurement (V_LL) has been newly developed and successfully applied for trickle bed reactor (TBR). The technique was validated against newly developed fiber optical probe technique. This work attempts to highlight the applicability of this newly developed technique on a liquid–solid packed bed reactor. In this work, liquid was represented by water and solids were represented by EPS beads. The EPS beads were chosen because of its low density property. Three superficial liquid velocities (V_SL) were applied to the system. The experiment was replicated four times. The digital industrial radiography (DIR) consists of a complementary metal oxide semiconductor (CMOS) digital detector and X-ray source. Results of this work suggest that the technique has been successfully applied and comparable with previous work that has been done in the literature. It also suggests that there will be a maximum measurable interstitial liquid velocity when it travel inside the packed bed. The measured V_LL can have a maximum range that is between 4 and 4.7 times that of its V_SL. For V_SL=0.42±±2%, the V_LL-Max is in between 1.7 cm/s and 1.9 cm/s, V_SL=0.84±±2%, the V_LL-Max is in between 3.6 cm/s and 4.0 cm/s, and for V_SL=1.11±±2%, the V_LL-Max is in between 4.3 cm/s and 4.8 cm/s.     

Simulation of thermal barrier coating spallation induced by the initiation/growth/coalescence of internal crack and interfacial crack based on a real image model

In the furnace cycle test, the growth of oxide film leads to the propagation and coalescence of multiple cracks near the interface, which should be responsible for the spallation of thermal barrier coatings (TBCs). A TBC model with real interface morphology is created, and the near-interface large pore is retained. The purpose of this work is to clarify the mechanism of TBC spallation caused by successive initiation, propagation, and linkage of cracks near the interface during thermal cycle. The dynamic growth of thermally grown oxide (TGO) is carried out by applying a stress-free strain. The crack nucleation and arbitrary path propagation in YSZ and TGO are simulated by the extended finite element method (XFEM). The debonding along the YSZ/TGO/BC interface is evaluated using a surface-based cohesive behavior. The large-scale pore in YSZ near the interface can initiate a new crack. The ceramic crack can propagate to the YSZ/TGO interface, which will accelerate the interfacial damage and debonding. For the TGO/BC interface, the normal compressive stress and small shear stress at the valley hinder the further crack propagation. The growth of YSZ crack and the formation of through-TGO crack are the main causes of TBC delamination. The accelerated BC oxidation increases the lateral growth strain of TGO, which will promote crack propagation and coalescence. The optimization design proposed in this work can provide another option for developing TBC with high durability.     

Additive manufacturing of bipolar plates for hydrogen production in proton exchange membrane water electrolysis cells

Water electrolysis is a process that can produce hydrogen in a clean way when renewable energy sources are used. This allows managing large renewable surpluses and transferring this energy to other sectors, such as industry or transport. Among the electrolytic technologies to produce hydrogen, proton exchange membrane (PEM) electrolysis is a promising alternative. One of the main components of PEM electrolysis cells are the bipolar plates, which are machined with a series of flow distribution channels, largely responsible for their performance and durability. In this work, AISI 316L stainless steel bipolar plates have been built by additive manufacturing (AM), using laser powder bed fusion (PBF-L) technology. These bipolar plates were subjected to ex-situ corrosion tests and assembled in an electrolysis cell to evaluate the polarization curve. Furthermore, the obtained results were compared with bipolar plates manufactured by conventional machining processes (MEC). The obtained experimental results are very similar for both manufacturing methods. This demonstrates the viability of the PBF-L technology to produce metal bipolar plates for PEM electrolyzers and opens the possibilities to design new and more complex flow distribution channels and to test these designs in initial phases before scaling them to larger surfaces.     

渣油催化临氢热转化钼铁双金属催化剂的研究

Replacement of precious single metal catalysts with cost-effective, highly-dispersed composite catalysts for catalytic hydrothermal conversion of residue holds tremendous promise for the residue upgrading technologies. Organic metals were added to the feed as the oil-soluble precursors, and transformed into the catalytic active phases in this work. Physical properties and structures of the composite catalysts had been investigated by X-ray fluorescence spectroscopy, X-ray photoelectron spectroscopy, X-ray diffraction, scanning electron microscope and transmission electron microscopy. The composite catalysts were found to be highly efficient in the catalytic hydrothermal conversion of both model compound and residue. Increased metal dispersion and synergistic effects of two metals played indispensable roles in such catalytic system. Results showed that under the test conditions in the article, the catalyst had the best catalytic performance when the mass ratio of molybdenum to iron was 1.5.