The computational fluid dynamics-population balance model (CFD-PBM) has been presented and used to evaluate the bubble behavior in a large-scale high pressure bubble column with an inner diameter of 300 mm and a height of 6600 mm. In the heterogeneous flow regime, bubbles can be divided into “large bubbles” and “small bubbles” by a critical bubble diameter dc. In this study, large and small bubbles were classified according to different slopes in the experiment only by the method of dynamic gas disengagement, the critical bubble diameter was determined to be 7 mm by the experimental results and the simulation values. In addition, the effects of superficial gas velocity, operating pressure, surface tension and viscosity on gas holdup of large and small bubbles in gas−liquid two-phase flow were investigated using a CFD-PBM coupling model. The results show that the gas holdup of small and large bubbles increases rapidly with the increase of superficial gas velocity. With the increase of p... [more]

In the past two decades, aseptic processing has been implemented in the food industry to sterilize particulate liquid food mixtures. To ensure that particulates in the liquid receive sufficient heating, mathematical modeling is employed to evaluate the temperature and lethality level in the particles. We developed a model for the thermal processing of liquid foods containing cubic particles in a continuous laminar pipe flow system, comprising a tubular heat exchanger. In our simplified approach, heat transfer equations for particulate liquid foods were solved analytically and numerically to evaluate the effect of certain process parameters on the time temperature profiles of particles and the lethality value in the products. A comparison of the particles’ lethality values was made between the experiment and simulation for two different particle residence times in a case study, and the model predictions were in good agreement with experimental data. Based on modeling studies, it was fou... [more]

Summary: Determining the cost of healthcare and social care for patients is a crucial issue for many parties; therefore, both public and private payments play a decisive role in patient care. The article deals with the analysis of the possibilities of the simulation of costs related to Alzheimer’s disease. This disease is highly variable, and the cost items vary considerably. Therefore, it is necessary to use simulation methods. The results of simulation models can then be implemented in sophisticated methods working with activity costs. Findings: Models for health and social care are specific. No significant re-engineering is expected in this area, so the models must be unambiguous and easy to understand for all representatives involved. Modeling of business processes is thus a suitable means for analysis in this area. The process cost calculation is built on two simulation models. The first model is the Business Process Model and Notation (BPMN), the results of which are verified in... [more]

Many pharmaceutical molecules, fine chemicals, and proteins exhibit liquid−liquid phase separation (LLPS, also known as oiling out) during solution crystallization. LLPS is of significant concern in crystallization process development, as oiling out can compromise the effectiveness of a crystallization and can lead to operational problems. A comprehensive methodology that allows a process scientist/engineer to characterize the various phase boundaries relevant to oiling out is currently lacking. In this work, we present a modeling framework useful in predicting the binodal, spinodal, and gelation boundaries starting from the solubility data of a solute that is prone to oiling out. We collate the necessary theoretical concepts from the literature and describe a unified approach to model the phase equilibria of solute−solvent systems from first principles. The modeling effort is validated using experimental data reported in the literature for various solute−solvent systems. The predictiv... [more]

The diffusion coefficient of gases in coal varies with time. This study aims to develop an unsteady dynamic diffusion (UDD) model based on the decay of diffusion coefficient with time and the change of integral. This study conducted a series of gas desorption and diffusion experiments with three different combinations of particle sizes and gas pressures and compared the diffusion coefficients of the three models. The UDD model exhibited good fitting results, and both the UDD and bidisperse models fitted the experimental data better than the unipore model. In addition, the dynamic diffusion coefficient (DDe) decreased rapidly in the initial stage but gradually decreased to a stable level in the later stage. All the effective diffusion coefficients of the three models negatively correlated with the particle size. In the unipore model, the diffusion coefficient of coal samples with three particle sizes increased with gas pressure. In the bidisperse and UDD models, the diffusion coefficien... [more]

Pipeline end termination (PLET) installation is an essential part of offshore pipe-laying operation. Pipe-laying operations are sensitive to pipe-laying barge motion and marine environmental conditions. Monitoring the field environment can provide a reasonable basis for planning pipe-laying. Therefore, the measurement and analysis of sea wave motion is helpful for the control and operational safety of the pipeline and vessels. In this study, an environmental monitoring system was established to measure wave motion during PLET operation. Fourier transforms were used to process images that were acquired by ultra-high-frequency X-band marine radar to extract wave parameters. The resulting wave spectra, as measured each minute, were used to simulate real-time wave data and calculate wave characteristics and regressed wave frequency and direction spectrum throughout the PLET operation. The regressed frequency, spectral density, and direction spectra were compared with the theoretical spectr... [more]

Spray cooling with liquid nitrogen is increasingly utilized as an efficient approach to achieve cryogenic cooling. Effects of injection mass flow rate fluctuations on the evaporation, temperature distribution, and droplet distribution of a spray field were examined by employing a validated Computational Fluid Dynamics (CFD) numerical model. The numerical results indicated that injection fluctuations enhanced the volume-averaging turbulent kinetic energy and promoted the evaporation of the whole spray field. The strengthened mass and heat transfer between the liquid nitrogen droplets and the surrounding vapor created by the fluctuating injection led to a lower temperature of the whole volume. A relatively smaller droplet size and a more inhomogeneous droplet distribution were obtained under the unsteady inlet. The changes of the frequency and the amplitude of the fluctuations had little effects on the overall spray development. The results could enrich the knowledge of the relation betw... [more]

The compact finite difference scheme on non-uniform meshes and the Fourier spectral hybrid method are used to directly simulate the evolution of vortex structures in a laminar boundary layer over a flat plate. To this end, two initial local vibration disturbances, namely, the positive−negative and the negative−positive models, at the wall were adopted. The numerical results show that the maximum amplitudes of vortex structures experience a process of linear growth and nonlinear rapid growth. The vertical disturbance velocity and mean flow shear and the derivative term of the stream-wise disturbance velocity and the span-wise disturbance velocity, are important factors for vortex structure development; the high- and low-speed stripe and the stream-wise vortex are consistent with structures seen in full turbulence. The maximum amplitude of the negative−positive model grows more quickly than that of the negative−positive model, and the detailed vortex structures are different for the two... [more]

The lumpy zone in a blast furnace is composed of piles formed naturally during burden charging. The properties of this zone have significant effects on the blast furnace operation, including heat and mass transfer, chemical reactions and gas flow. The properties of the layers mainly include the angle of repose and porosity distribution. This paper introduces two methods, the Discharging Method and the Lifting Method, to study the influence of the packing method on the angle of repose of the pile. The relationships of the angle of repose and porosity with physical parameters are also investigated. The porosity distribution in the bottom of a pile shows a decreasing trend from the region below the apex to the center. The coordination number of the particles is employed to explain this change. The maximum of the frequency distribution of it was found to show a negative correlation to the static friction coefficient, but becomes insensitive to the parameter as the static friction coefficie... [more]

Nanofiltration of aqueous solutions of succinic acid with the addition of sodium hydroxide or magnesium hydroxycarbonate has been investigated experimentally and modeled with the comprehensively described Donnan−Steric partitioning model. The experimental retentions of acid at the same pH varied between 16% and 78%, while the estimated total volume membrane charge densities were in the range of −35.73 and +875.69 mol/m3. This work presents a novel insight into the modeling of nanofiltration and investigates the relations between the estimated total volume membrane charge densities, ionic strength, and component concentration on the performance of ceramic membrane. In addition, this study takes into consideration other parameters such as pH regulation and viscosities of solutions.

Pore network models are powerful tools to simulate invasion and transport processes in porous media. They are widely applied in the field of geology and the drying of porous media, and have recently also received attention in fuel cell applications. Here we want to describe and discuss how pore network models can be used as a prescriptive tool for future water electrolysis technologies. In detail, we suggest in a first approach a pore network model of drainage for the prediction of the oxygen and water invasion process inside the anodic porous transport layer at high current densities. We neglect wetting liquid films and show that, in this situation, numerous isolated liquid clusters develop when oxygen invades the pore network. In the simulation with narrow pore size distribution, the volumetric ratio of the liquid transporting clusters connected between the catalyst layer and the water supply channel is only around 3% of the total liquid volume contained inside the pore network at th... [more]

This paper proposes a novel computationally efficient, easy-to-implement electromagnetic hysteresis based dynamics model of a kind of intelligent electromagnetic torque controlled coupling (EMTC), which, with drag torque under consideration, first models the electromagnetic hysteresis existing in the primary clutch with the classical Preisach model, and then models the transferred torques in the three friction elements of the center coupling in the slipping and locked modes, respectively. The performance of the model is verified by simulation and experiment jointly, which lays the basis for the development of advanced control algorithm.

This paper presents an experimental evaluation of a closed loop lake water heat pump (LWHPs) system based on the slinky coiled configuration. Initially, a mathematical model is developed in the Engineering Equation Solver (EES) for the heat pump system and the submerged coils in a lake. System performance is determined for the submerged slinky copper coils under the various operating conditions. Afterwards, parametric analysis is performed considering different influencing parameters, such as the lake water temperature, ambient temperature, and mass flow rate of the circulating fluid at constant lake depth of 4 ft. The experimental setup is developed for 3.51 kW cooling capacity after cooling load calculation for a small room. In the current study, slinky copper coils are used to exchange heat with lake water. The experimental setup is installed in Taxila, Pakistan, and the system’s performance is analyzed during selected days. After experimentation based on hourly and daily operation... [more]

The mixing and disinfection performance of a full-scale chlorine contact tank (CCT) is thoroughly investigated by means of numerical simulations for seasonal water supply variations in the water treatment plant (WTP) of Eskisehir in Turkey. Velocity measurements and tracer studies are carried out on a 1:10 scale laboratory model of the CCT to validate the numerical model. A good agreement between numerical and experimental results shows that the numerical model developed can be reliably used for the simulation of turbulent flow and solute transport in the full-scale CCT. Tracer studies indicate that the hydraulic performance of the CCT is classified as “average” according to the baffling factor, while the Morrill, Aral-Demirel (AD), and dispersion indexes indicate low mixing due to the recirculating and short-circuiting effects inside the chambers of the CCT. With respect to the first order modeling of chlorine decay and pathogen inactivation, chlorine concentrations are found to be si... [more]

The performance of a CO2 double-pipe evaporator was studied through experiments and a simulation model that was established by the steady-state distribution parameter method and experimentally verified while using a CO2 transcritical water‒water heat pump system. The effects of different operating parameters on heat transfer performance were studied over a range of evaporation temperatures (−5 to 5 °C), mass velocity (100‒600 kg/m2s), and heat flux (5000‒15,000 W/m2). It was found that the dryout quality increased at a small evaporation temperature, a large mass velocity, and a small heat flux. The simulation yield means relative error (RE) of heat transfer for the evaporation temperature and that of the CO2 pressure drop for the chilled water inlet temperature were 5.21% and 3.78%, respectively. The effect of tube diameter on the performance of CO2 double-pipe evaporator is probed through simulations. At the same time, this paper defines a parameter α , which is the proportio... [more]

The numerical method was used to investigate heat and moisture transport during dehydration of bananas from microwave heating. COMSOL multi-physics software was employed to perform the simulation task. A banana is defined as a porous medium. It has constituents of water, vapor, air as the liquid phase and a solid porous matrix. The numerical results of this study were validated with experimental data. The profiles of moisture, vapor and pressure are discussed in this study. Moreover, the effects of the ripening stages of the banana are examined. A higher heat flux was observed from the beginning period along with the increasing time steps until 50 s. Heat generation decreased during 50 s to 60 s, coinciding with a small rise in temperature, but the temperature gradient remained constant. The temperature distribution of both unripe and ripe banana samples was non-uniform. At the center of the banana, the temperature increased rapidly and reached its highest temperature with the negative... [more]

With the increasing use of Computational Fluid Dynamics to investigate multiphase flow scenarios, modelling surface tension effects has been a topic of active research. A well known associated problem is the generation of spurious velocities (or currents), arising due to inaccuracies in calculations of the surface tension force. These spurious currents cause nonphysical flows which can adversely affect the predictive capability of these simulations. In this paper, we implement the Continuum Surface Force (CSF), Smoothed CSF and Sharp Surface Force (SSF) models in OpenFOAM. The models were validated for various multiphase flow scenarios for Capillary numbers of 10 − 3 −10. All the surface tension models provide reasonable agreement with benchmarking data for rising bubble simulations. Both CSF and SSF models successfully predicted the capillary rise between two parallel plates, but Smoothed CSF could not provide reliable results. The evolution of spurious current were studied f... [more]

In the current research, a high-pressure submerged cavitation jet is investigated numerically. A cavitation model is created considering the effect of shear stress on cavitation formation. As such, this model is developed to predict the cavitation jet, and then the numerical results are validated by high-speed photography experiment. The turbulence viscosity of the renormalization group (RNG) k-ε turbulence model is used to provide a flow field for the cavitation model. Furthermore, this model is modified using a filter-based density correction model (FBDCM). The characteristics of the convergent-divergent cavitation nozzle are investigated in detail using the current CFD simulation method. It is found that shear stress plays an important role in the cavitation formation in the high-pressure submerged jet. In the result predicted by the Zwart-Gerber-Belamri (ZGB) cavitation model, where critical static pressure is used for the threshold of cavitation inception, the cavitation bubble on... [more]

Digital twins are rigorous mathematical models that can be used to represent the operation of real systems. This connection allows for deeper understanding of the actual states of the analyzed system through estimation of variables that are difficult to measure otherwise. In this context, the present manuscript describes the successful implementation of a digital twin to represent a four-stage multi-effect evaporation train from an industrial sugar-cane processing unit. Particularly, the complex phenomenological effects, including the coupling between thermodynamic and fluid dynamic effects, and the low level of instrumentation in the plant constitute major challenges for adequate process operation. For this reason, dynamic mass and energy balances were developed, implemented and validated with actual industrial data, in order to provide process information for decision-making in real time. For example, the digital twin was able to indicate failure of process sensors and to provide est... [more]

Carbon solution loss reaction of coke gasification is one of the most important reasons for coke deterioration and degradation in a blast furnace. It also affects the permeability of gas and fluids, as well as stable working conditions. In this paper, a three dimensional model is established based on the operational parameters of blast furnace B in Bayi Steel. The model is then used to calculate the effects of oxygen enrichment, coke oven gas injection, and steel scrap charging on the carbon solution loss ratio of coke in the blast furnace. Results show that the carbon solution loss ratio of coke gasification for blast furnace B is almost 20% since the results of a model are probably only indicative. The oxygen enrichment and the addition of steel scrap can reduce the carbon solution loss ratio with little effect on the working condition. However, coke oven gas injection increases the carbon solution loss ratio. Therefore, coke oven gas should not be injected into the blast furnace unl... [more]

Metal−organic frameworks (MOFs) are the porous, crystalline structures made of metal−ligands and organic linkers that have applications in gas storage, gas separation, and catalysis. Several experimental and computational tools have been developed over the past decade to investigate the performance of MOFs for such applications. However, the experimental synthesis of MOFs is still empirical and requires trial and error to produce desired structures, which is due to a limited understanding of the mechanism and factors affecting the crystallization of MOFs. Here, we show for the first time a comprehensive kinetic model coupled with population balance model to elucidate the mechanism of MOF synthesis and to estimate size distribution of MOFs growing in a solution of metal−ligand and organic linker. The oligomerization reactions involving metal−ligand and organic linker produce secondary building units (SBUs), which then aggregate slowly to yield MOFs. The formation of secondary building u... [more]

Heat stimulation of coalbed methane (CBM) reservoirs has remarkable promotion to gas desorption that enhances gas recovery. However, coalbed deformation, methane delivery and heat transport interplay each other during the stimulation process. This paper experimentally validated the evolutions of gas sorption and coal permeability under variable temperature. Then, a completely coupled heat-gas-coal model was theoretically developed and applied to a computational simulation of CBM thermal recovery based on a finite element approach of COMSOL with MATLAB. Modeling and simulation results show that: Although different heat-gas-coal interactions have different effects on CBM recovery, thermal stimulation of coalbed can promote methane production effectively. However, CBM thermal recovery needs a forerunner heating time before the apparent enhancement of production. The modeling and simulation results may improve the current cognitions of CBM thermal recovery.

In this work, we employed a computational fluid dynamics (CFD)-based model with a Eulerian multiphase approach to simulate the fluidization hydrodynamics in biomass gasification processes. Air was used as the gasifying/fluidizing agent and entered the gasifier at the bottom which subsequently fluidized the solid particles inside the reactor column. The momentum exchange related to the gas-phase was simulated by considering various viscous models (i.e., laminar and turbulence models of the re-normalisation group (RNG), k-ε and k-ω). The pressure drop gradient obtained by employing each viscous model was plotted for different superficial velocities and compared with the experimental data for validation. The turbulent model of RNG k-Ɛ was found to best represent the actual process. We also studied the effect of air distributor plates with different pore diameters (2, 3 and 5 mm) on the momentum of the fluidizing fluid. The plate with 3-mm pores showed larger turbulent viscosities above th... [more]

This editorial provides a brief overview of the Special Issue “Modeling and Simulation of Energy Systems.” This Special Issue contains 21 research articles describing some of the latest advances in energy systems engineering that use modeling and simulation as a key part of the problem-solving methodology. Although the specific computer tools and software chosen for the job are quite variable, the overall objectives are the same—mathematical models of energy systems are used to describe real phenomena and answer important questions that, due to the hugeness or complexity of the systems of interest, cannot be answered experimentally on the lab bench. The topics explored relate to the conceptual process design of new energy systems and energy networks, the design and operation of controllers for improved energy systems performance or safety, and finding optimal operating strategies for complex systems given highly variable and dynamic environments. Application areas include electric powe... [more]

Monitoring and control of the blast furnace hearth is critical to achieve the required production levels and adequate process operation, as well as to extend the campaign length. Because of the complexity of the draining, the outflows of iron and slag may progress in different ways during tapping in large blast furnaces. To categorize the hearth draining behavior, principal component analysis (PCA) was applied to two extensive sets of process data from an operating blast furnace with three tapholes in order to develop an interpretation of the outflow patterns. Representing the complex outflow patterns in low dimensions made it possible to study and illustrate the time evolution of the drainage, as well as to detect similarities and differences in the performance of the tapholes. The model was used to explain the observations of other variables and factors that are known to be affected by, or affect, the state of the hearth, such as stoppages, liquid levels, and tap duration.