The effects of processing intensity, time and particle surface energy on mixing of binary cohesive blends (size ratio 1:2, fine concentration at 10 %) in high intensity vibration system were investigated via DEM simulations. Results show that both increasing processing intensity from 50 to 100 Gs and reducing surface energy from 50 to 0.5 J/m2 lead to a faster mixing rate. Mixing Bond number (〖Bo〗_m) was introduced to capture the effective mixing rate, Rm; higher 〖Bo〗_m corresponding to lower mixing rate. The coefficient of variation, Cv, formed the basis for the mixing quality and Rm, while the mixing action is quantified by the product of Rm and mixing time (Pr,t). Simulation results show that Cv values drop initially, and then rise with Pr,t. Hence, low Pr,t indicates inadequate mixing intensity, while high Pr,t most likely indicates mixture segregation, and therefore too high or too low Pr,t values should be avoided.

Bubble breakup plays an important role in gas-liquid flows, but detailed studies are still scarce. In this work, the breakup behavior of a single bubble in a stirred tank was experimentally studied with a high-speed camera, focusing on the effect of gas density, liquid properties, agitation speed and mother bubble size. The bubble breakup time, breakup probability, breakup rate and daughter bubble size distribution were determined. The internal flow phenomenon inside a deformed bubble was studied in detail, which accounted for the effect of gas density or operating pressure. The results showed that with increasing gas density, agitation speed, mother bubble size and decreasing surface tension, the bubble breakup rate and probability of equal-size distribution significantly increased. With increasing liquid viscosity, the bubble breakup rate decreased especially in the high viscosity range. An M-shaped daughter bubble size distribution was observed, which was consistent with our previous bubble breakup model.

Reactor corrosion and salt deposition problems severely restrict the industrialization of supercritical water oxidation. Transpiring wall reactor can effectively weaken these two problems through a protective water film formed on its internal surface. In this work, the effects of key structural parameters on water film properties of transpiring wall reactor were explored by numerical simulation, and established models were validated by comparing simulation and experimental values. The results show that transpiration water layer, transpiring wall porosity and inner diameter hardly affected organic matter degradation. Increasing transpiration water layer and transpiring wall porosity reduced reactor center temperatures in the middle and lower zones of the reactor. Increasing transpiration water layer, transpiring wall porosity and inner diameter decreased water film temperatures but increased water film coverage rates. Increasing reactor length affected slightly on the volume of the upper supercritical oxidation zone but enlarged the subcritical zone.

We present the development and application of a two-phase stirred reactor model for heavy oil upgrading in the presence of supercritical water (SCW), with coupled phase-specific thermolysis reaction kinetics and multicomponent hydrocarbon water phase equilibrium. We demonstrate the inference of oil and water phase kinetics parameters for a compact lumped reaction kinetics model through the application of this model to two different sets of batch reactor experiments reported in the literature. We infer that, though SCW can suppress the formation of newer polynuclear aromatics (PNA) from distillate range species, it is broadly ineffective in deterring the combination of pre-existing PNA fragments in the oil feed. Quantification of the conversion to distillate liquids before the onset of coke formation helps arrive at a clearer conclusion on whether the use of SCW in the batch reactor leads to better product outcomes for different oil feeds and operating conditions.

We investigate the drying process of monodisperse colloidal film over a wide range of Péclet number (Pe) by using the Brownian dynamics simulation. We analyze the detailed process in three aspects; accumulation front, normal stress, and microstructure. The evolution of particle distribution is quantified by tracking the accumulation front. The accumulated particles contribute to the continuous increase of the normal stress at the interface. At the substrate, the normal stress first stays constant and then increases as the accumulation front touches the substrate. We quantitatively analyze the stress development by a scaled normal stress difference between the two boundaries. At all tested Pe, the stress difference increases to the maximum, followed by a decrease during drying. Interestingly, a mismatch is observed between the stress difference maximum and the initial stress increase at the substrate. The microstructural analysis reveals that this mismatch is related to the microstructural development at the substrate.

Centrifugal pumps are used in several industrial processes. It is common the operation of this equipment with gas-liquid mixtures, which is the case of the electrical submersible pumping artificial lift method used in the oil industry. The increase of free gas fraction inside the pump may lead to unstable operation and problems such as surging and gas locking phenomena to occur. In this study a drift-flux model is proposed for the gas-liquid flow subjected to centrifugal fields using the impeller as an example. The model is closed with experimental data of bubble diameter, displacements and velocities acquired via high-speed camera at several different rotational speeds and gas mass flow rates using water as the continuous medium. From the modeling and the forces balance in the bubbles, a quantitative criterion for the start point of surging and gas locking conditions was proposed.

Hydrophobic deep eutectic solvents (DESs) have emerged as excellent extractants. Their performance depends on the heterogeneous hydrogen bond environment formed by multiple hydrogen bond donors and acceptors. An understanding of this heterogeneous hydrogen bond environment can be used to develop principles for designing high-performance DESs for extraction and other separation applications. We investigate the structure and dynamics of hydrogen bonds in eight hydrophobic DESs formed by decanoic acid, menthol, thymol, and Lidocaine using molecular dynamics simulations. The results show the diversity of hydrogen bonds in the eight DESs and their impact on diffusivity and molecular association. Each DES possesses four-six types of hydrogen bonds and one or two of them overwhelm the others in quantity and lifetime. The dominating hydrogen bonds determine whether the DESs are governed by intra- or inter-component associations. The component diffusivity presents an inverse relationship with the hydrogen bond strength.

First principles studies combined with the microkinetic analysis were performed to study the reliability and reaction mechanisms of single-atom doped graphene (SADGr) materials in catalyzing NOx reduction with CO. By screening the 3d transition metals (Sc-Zn) and group-IV elements (Si and Ge), it was found that the Ti and Co doped graphene sheets (TiGr and CoGr) respectively own excellent catalytic activities in the NO/NO2-to-N2O and the N2O-to-N2 processes at low temperatures. Therefore, the TiGr/CoGr composite can be a promising catalyst in NOx reduction with CO. It was further revealed the combination of adsorption energy and electronegativity was a good descriptor to predict the activation energies. The obtained results can provide useful information for rational design of carbon-based single-atom catalysts for NOx reduction by CO at low temperatures.

Product formulations for industrial processes are typically developed at lab-scale. However, the mixing conditions are not easily mimicked in the lab. A rotational device is proposed in this work as a fast lab-scale formulation development, which enables mimicking the mixing conditions in the industrial process. The geometrical configurations of the rotational device are from rheometry devices (plate-plate and cone-plate). The main advantages of this method are the small amounts of raw materials, and shorter testing times. This methodology is applied to an industrial case study, the Reaction Injection Moulding (RIM) process. The mixing length scales evolution in the rotational rheometer were matched to those in RIM machines. The main novelty of this work is the introduction of a protocol that bridges the processing conditions at lab using small amounts of raw materials to high throughput continuous flow reactors.

In this study, effect of swirling addition on the liquid mixing behavior of multi-orifice-impinging transverse jet mixer has been investigated by planar laser induced fluorescence as well as large eddy simulation (LES). In the case of swirling addition into the jet flow, there exists an optimized swirling jet angle or optimized jet-to-cross velocity ratio for the fixed mixer configuration. A larger swirling jet angle will make the flow dominated by the swirling, resulting in a slower mixing process. Interaction of swirling crossflow with no-swirling injected streams, or with swirling injected streams in the opposite direction is beneficial for the mixing. LES predictions show that many small vortices are produced homogenously due to intensified impingement in the case of opposite swirling directions, leading to a relative fast mixing process in several milliseconds. Whereas the mixing is restrained when the swirling directions of two flows are the same.

The study proposed an isotopes-tagging method for investigation of reactions under the atmosphere of product gas. To illustrate this method, the calcination kinetics of calcium carbonate Ca13CO3 in CO2 atmospheres were investigated by monitoring 13CO2 produced using a micro fluidized bed reaction analyzer (MFBRA). The results demonstrated that the presence of CO2 in the reaction atmosphere increases the apparent activation energy. The increase in the apparent activation energy is, however, significantly overestimated by the TGA because of the excessive suppression by stagnated product gas inside the sample crucible. Comparatively, the results from the MFBRA are due primarily to the thermal equilibrium limitation, because the gas diffusion in the MFBRA is essentially eliminated. It is thus concluded that the MFBRA is quite capable of acquiring the real kinetics of reactions in such inhibitory atmospheres. atmosphere.

As an attractive alternative technology for the separation of long chain olefin and paraffin, a novel silver-based deep eutectic solvent (Ag-DES) was prepared and utilized for 1-octene/n-octane separation. Comprehensive reactive extraction separation experiments were performed to highlight the Ag-DES concentration and operating temperature discriminations using compounds with different ratio of 1-octene/n-octane. The novel Ag-DES showed optimal separation performance regarding 1-octene/n-octane and possessed the highest levels separation selectivity in the range 3.75-16.74 with excellent circulation stability in our best knowledge. Furthermore, FT-Raman measurements and quantum chemistry calculation were conducted to elucidate the interaction mechanism of Ag-DES in the separation of 1-octene and n-octane, which revealed that both chemical complexation and strong physical attraction existed in the complex of Ag-DES with 1-octene rather than n-octane. This study lends important insight for the development of Ag-DES reactive extraction separation process for the energy-efficient long chain α-olefin purification from F-T synthesis products.

The existing methods of flexibility index are mainly based on mixed-integer linear or nonlinear programming methods, making it difficult to readily deal with complex mathematical models. In this article, a novel solution strategy is proposed for finding a reliable upper bound of the flexibility index where the process model is implemented in a black box that can be directly executed by a commercial simulator, and also avoiding the need for calculating derivatives. Then, the flexibility index problem is formulated as a sequence of univariate derivative-free optimization (DFO) models. An external DFO solver based on trust-region methods can be called to solve this model. Finally, after calculating the critical point of the model parameters, the vertex enumeration method and two gradient approximation methods are proposed to evaluate the impact of process parameters and to evaluate the flexibility index. A reaction model is studied to show the efficiency of the proposed algorithm.

The low effectiveness factor of catalyst pellet caused by high internal diffusion limitation is a common issue in fixed-bed reactor. Nevertheless, hierarchical structured catalyst provides a promising solution for the contradiction between reaction activity and diffusion efficiency in large catalyst pellets. Herein, we studied the effect of pore structure parameters of the meso-macroporous catalyst on Fischer-Tropsch synthesis performances through experiment and pellet scale reaction-diffusion simulation. The pellet simulation firstly elucidated the reason for the significant improvement on activity and product selectivity for the meso-macroporous catalyst observed in our experiment. Further optimization via pellet simulation indicated the critical influences of wax filling degree and that the perfect matching between reaction and mass transfer rates by increasing macropore size and adjusting porosity within pellet enables the C5+ space-time yield to the maximum. This work could provide a theoretical guideline for the engineering design of the hierarchical structured catalyst pellet.

The complexity of lipid feedstocks and the lack of data on physical properties hinder the simulation of oleochemical processing units. In this work, an iterative lumping approach is proposed to define an adequate number of key components such that diversification between lipid feedstocks becomes possible, while keeping the determination of physical properties as required for process modelling manageable. As a case study, the iterative lumping approach is used for simulation and optimization of a fatty acid distillation plant. For predicting vapour-liquid equilibria of fatty acids, the best results were acquired using the property method UNIQ-HOC. Using the iterative lumping approach, 11 key components were selected to represent the feedstock. The process model properly predicts the product composition, yield, purity and heat duty. The most important process parameters are found to be side-reflux-ratio, reboiler-outlet-temperature and heat-duty of the pitch-distiller. For optimization, an increase of the side-reflux-ratio and reboiler-outlet-temperature, is recommended.

In this work, methods based upon nonequilibrium thermodynamics are elucidated to predict stationary states of chemical reactions in nonequilibrium plasma, and limits for energy conversion efficiency. Two example reactions are used: CO2 splitting and NH3 synthesis, with emphasis on CO2 splitting. Expectations from the theoretical framework are compared to experimental results for both reactions, and reasonable agreement is obtained. The conclusion is that the probability of observing either reactants or products increases with the amount of energy dissipated by that side of the reaction as heat through collisions with hot electrons. The side of the reaction that dissipates more energy as heat has a higher probability of occurrence. Furthermore, endergonic chemical reactions in nonequilibrium plasma, such as CO2 splitting at low temperature, require an intrinsic energy dissipation to satisfy the 2nd law of thermodynamics – a sufficient and necessary waste. This intrinsic dissipation limits the maximum theoretical energy conversion efficiency

Aging effects of off-gas streams including dry air and humid air on reduced silver exchanged mordenite (Ag0Z) were studied. Aged Ag0Z was prepared by exposing Ag0Z to dry air and humid air at different aging temperatures, time, and water vapor concentrations. Iodine loading capacity on the aged Ag0Z was obtained through a continuous-flow adsorption system. Significant iodine loading capacity losses were observed after the Ag0Z was exposed to dry air and humid air. Physical and chemical analyses were conducted to observe the physical and chemical changes of Ag0Z after being aged. From iodine adsorption data and sample analyses, it was found that iodine loading capacity on the aged Ag0Z in dry air and humid air decreases with increasing aging temperatures, time and water vapor concentrations. The pseudo reaction model describes experimental data well and the oxidation of Ag0 is the rate determining step in the aging process.

Biotechnological application of multiple enzymes in different phases for target compounds synthesis poses a significant challenge for industrial process development. At the same time, a growing demand for natural flavors and fragrances opens up possibilities for novel biotechnological processes to replace current chemical synthesis routes, with additional advantages such as avoiding harsh reaction conditions and toxic chemicals, and less by-products in the system. Within complex biotechnological processes, the key for unfolding their industrial application potential in bioprocess engineering lies in their mathematical modeling. In this contribution, a multi-enzyme cascade reaction in a two-phase system implemented in a miniplant-scale reactor setup is mathematically modeled for the example of the flavoring agent cinnamyl cinnamate. Using our validated model and a mathematical optimization tool based on a genetic algorithm, optimization runs are performed to demonstrate the potential of computer-aided process development for complex biotechnological processes.