The new generation of internal combustion engines is facing various research challenges which often include modern fuels and different operating modes. A robust modeling framework is essential for predicting the dynamic behavior of such complex phenomena. In this article, the implementation, verification, and validation of a Eulerian multi-fluid model for spray applications within the OpenFOAM toolbox are presented. Due to its open-source nature and broad-spectrum of available libraries and solvers, OpenFOAM is an ideal platform for academic research. The proposed work utilizes advanced interfacial momentum transfer models to capture the behavior of deforming droplets at a high phase fraction. Furthermore, the WAVE breakup model is employed for the transfer of mass from larger to smaller droplet classes. The work gives detailed instructions regarding the numerical implementation, with a dedicated section dealing with the implementation of the breakup model within the Eulerian multi-flu... [more]

Positive displacement machines have been identified as appropriate expanders for small-scale power generation systems such as Organic Rankine Cycles (ORCs). Screw expanders can operate with good efficiency in working fluids under both dry and two-phase conditions. Detailed understanding of the fluid expansion process is required to optimise the machine design and operation for specific applications, and accurate design tools are therefore essential. Using experimental data for air expansion, both CFD and chamber models have been applied to investigate the influence of port flow and leakage on the expansion process. Both models are shown to predict pressure variation and power output with good accuracy. The validated chamber model is then used to identify the optimal volume ratio and rotational speed for experimental conditions.

Windcatcher louvers are designed to capture air flowing outside a building in order to increase its natural ventilation. There are no studies that have designed the shape of the louver to increase the natural ventilation efficiency of the building. This study aimed to conduct a computational fluid dynamics simulation and mock-up test of a Clark Y airfoil-type windcatcher louver designed to increase the natural ventilation in a building. The following test results were obtained. The optimal angle of attack of the airfoil was calculated via a numerical analysis, which demonstrated that the wind speed was at its highest when the angle of attack was 8°; further, flow separation occurred at angles exceeding 8°, at which point the wind speed began to decrease. The results of the mock-up test demonstrated that the time required to reduce the concentration of fine particles in the indoor air was 120 s shorter when the windcatcher was installed than when it was not, which indicating that the ti... [more]

Hydrothermal carbonization (HTC) represents one of the emerging and most promising technologies for upgrading biomass. Among the residual biomass waste, olive pomace and olive mill wastewater may be seen as valuable energy sources, especially for the Mediterranean countries, given the key role of the olive oil industry in those regions. This paper deals with the thermo-fluid dynamic performance of the HTC process of olive pomace. Computational Fluid Dynamics (CFD) modeling is employed in this study to numerically simulate such a process in batch reactor with the aim of understanding the complex fluid dynamics, heat transfer and reaction kinetics phenomena occurring under hydrothermal conditions. A parametric analysis is performed to evaluate the temperature fields inside the reactor and the output mass yields as a function of the power input required by the process. Velocity flow fields and the spatial distribution of the mixture during the process are also investigated to understand t... [more]

Annular seals of turbomachinery usually suffer from various degrees of eccentricities and disturbances due to the rotor−stator misalignment and radial loads, while the discussion of annular seal under both large static eccentricities and dynamic disturbances is relatively limited. In this paper, the applicability of linear assumption and reliability of nonlinear dynamic model for eccentric annular seals under large eccentricities and disturbances is discussed based on the investigation of seals with various rotor motions through computational fluid dynamics (CFD). After the validation of transient CFD methods by comparison with experimental and bulk theory results, the dynamic behaviors of annular seal are analyzed by adopting both direct transient simulations and the nonlinear Muszynska model. The results show that the nonlinear dynamic model based on rotor circular whirls around seal center can predict the fluid excitations of different types of rotor motions well under small static... [more]

The design of gear pumps and motors is focused on more efficient units which are possible to achieve using advanced numerical simulation techniques. The flow that appears inside the gear pump is very complex, despite the simple design of the pump itself. The identification of fluid flow phenomena in areas inside the pump, considering the entire range of operating parameters, is a major challenge. This paper presents the results of simulation studies of leakages in axial and radial gaps in an external gear pump carried out for different gap shapes and sizes, as well as various operating parameters. To investigate the processes that affect pump efficiency and visualize the fluid flow phenomena during the pump’s operation, a CFD model was built. It allows for a detailed analysis of the impact of the gears’ eccentricity on leakages and pressure build-up on the circumference. Performed simulations made it possible to indicate the relationship between leakages resulting from the axial and ra... [more]

The present study focuses on investigating the aerodynamic interaction between a three-element wing and wheel in ground effect, following the Formula One regulation change set for 2022, among which is the simplification of the front wing. This was accomplished by conducting a three-dimensional computational fluid dynamics analysis, using a Detached-Eddy Simulation approach, on a simplified one-quarter model of a Formula One racing car. The main goal was to examine how changing the front wing pressure distribution, by changing the incidence of the second flap, affected the wheel wake. The flow investigation indicated that the wheel wake is influenced by the flap configuration, which is mainly due to the fact that different flap configurations produce different upwash flow fields, leading to a variation of the separation point on top of the tire. As the separation point moves rearwards, the downwash generated in the central region (for a vertical plane) of the wheel wake increases increm... [more]

The main challenge of improving spark ignition (SI) engines to achieve ever increasing thermal efficiencies and near-zero pollutant emissions today concerns developing turbulent combustion under homogeneous ultra-lean premixed mixtures (HULP). This continuous shift of the lean operation limit entails questions on the applicability limits of the combustion models used to date for SI engine design and optimization. In this work, an assessment of flamelet-based models, widely used in RANS SI engines simulations of premixed turbulent combustion, is carried out using an open-source 3D-CFD platform to clarify the applicability limits on HULP mixtures. Two different consolidated approaches are selected: the Coherent Flame Model (CFM) and the Flame Area Model (FAM). Both methodologies are embedded by the authors into the same numerical structure and compared against measurements over a simplified and controlled flame configuration, which is representative of engine-like conditions. The experim... [more]

Zeolite 13X molecular sieve with high sorption capacity and significant sorption rate has been considered a promising candidate for seasonal heat storage. In this study, a code is developed to simulate the adsorption process between zeolite and water in all ranges of partial pressures, temperatures, and sorbate loadings. The results from the proposed code were compared with experiments and good agreement was observed. After validation, the developed model was used to study the effective parameters involved in the adsorption process of binder-free Zeolite 13X. A parametric study considering various temperatures and water content in the inflow air was conducted and the influence of different factors on the outlet temperature and adsorption enthalpy has been studied. This parametric study gives a good insight into the measures which can be taken for achieving the desired released energy or having the outlet temperature in the preferred range. The simulations have been conducted in a varie... [more]

This study aimed to analyze pharyngeal airflow characteristics and their relationship with the skeletal movement of the maxilla and mandible after bimaxillary orthognathic surgery in patients with skeletal class III (mandibular prognathism) malocclusion. Cone-beam computed tomography (CBCT) was conducted before surgery (T0), immediately after surgery (T1), and at least six months after surgery (T2). Digital imaging and communications in medicine files were transferred to InVivo (Anatomage) software to measure the skeletal changes after surgery. The changes in the maxillary and mandibular position, tongue position, and hyoid bone position were analyzed. Patient-specific models were reconstructed using 3D-Doctor software. The models after converting to the stereolithography (STL) file for Ansys integrated computer engineering and manufacturing code for computational fluid dynamics (ICEM CFD), commercial software were used for calculating the geometry, pressure drop and adjusted pressure... [more]

Heterogeneous chemical looping combustion (CLC) reactions and conjugate transports in large oxygen carrier particles were numerically investigated with computational fluid dynamics (CFD) approaches, in which a simplified noncatalytic reaction model was implemented for reducing intraparticle modelling computation. Volumic gas-solid reactions were treated as surface reactions based on the equivalent internal surface in the particle model. In large porous particles such as fixed bed CLC reactors, the heterogeneous reactions are often limited by intraparticle diffusion. Comprehensive analyses were conducted on transports across the particle surface and their influences on reactions inside the single particles. A threshold Reynolds number of external convections was found for the enhancement of intraparticle reactions. The heterogeneous reactions, intraparticle diffusions and interstitial transports in a fixed bed CLC reactor randomly packed with 597 spheres were thoroughly analysed with th... [more]

Understanding the influence of gas flow maldistribution in honeycombs can be beneficial for the process design in various technical applications. Although recent studies have investigated the effect of maldistribution by comparing the results of numerical simulations with experimental measurements, an exhaustive 3D full-field comparison is still lacking. Such full-field comparisons are required to identify and eliminate possible limitations of numerical and experimental tools. For that purpose, spatially resolved flow patterns were simulated by computational fluid dynamics (CFD) and measured experimentally by non-invasive NMR velocimetry (MRV). While the latter might suffer from a misinterpretation of artefacts, the reliability of CFD is linked to correctly chosen boundary conditions. Here, a full-field numerical and experimental analysis of the gas flow within catalytic honeycombs is presented. The velocity field of thermally polarized methane gas was measured in a regular 3D-printed... [more]

Hydrogen-fueled fuel cells are considered one of the key strategies to tackle the achievement of fully-sustainable mobility. The transportation sector is paying significant attention to the development and industrialization of proton exchange membrane fuel cells (PEMFC) to be introduced alongside batteries, reaching the goal of complete de-carbonization. In this paper a multi-phase, multi-component, and non-isothermal 3D-CFD model is presented to simulate the fluid, heat, and charge transport processes developing inside a hydrogen/air PEMFC with a serpentine-type gas distributor. Model results are compared against experimental data in terms of polarization and power density curves, including an improved formulation of exchange current density at the cathode catalyst layer, improving the simulation results’ accuracy in the activation-dominated region. Then, 3D-CFD fields of reactants’ delivery to the active electrochemical surface, reaction rates, temperature distributions, and liquid w... [more]

Silicone fluids belong to the group of pseudoplastic non-Newtonian fluids with complex rheological characteristics. They are considered in basic and applied researches and in a wide range of industrial applications due to their favorable physical and thermal properties. One of their specific field of applications in the automotive industry is the working fluid of viscous torsional vibration dampers. For numerical studies in the design and development phase of this damping product, it is essential to have thorough rheological knowledge and mathematical description about the silicone oil viscosity. In the present work, adopted rheological measurement results conducted on polydimethylsiloxane manufactured by Wacker Chemie with initial viscosity of 1000 Pas (AK 1 000 000 STAB silicone oil) are processed. As a result of the parameter identification by nonlinear regression, the temperature-dependent parameter curves of the Carreau−Yasuda non-Newtonian viscosity model are generated. By implem... [more]

Aiming at providing cost-effective approach for upgrading the existing municipal wastewater treatment plants in the cold region of China to meet more stringent discharge standards of nitrogen removal, a full-scale sewage treatment plant with the CASS process was selected through focusing on biological process, key equipment and hydrodynamics in bioreactors by the activated sludge model 1 (ASM1) and computational fluid dynamics (CFD) model. Influent COD fractions and the key characteristic parameters (YH and bH) of the activated sludge were determined through the respirometry at temperatures of 10 °C and 20 °C, respectively. The layout of submerged agitator installation in the bioreactor of the CASS process was optimized through CFD simulation. The calculation of the average relative deviation (less than 20%) between simulated data and the operation data, demonstrated that the ASM1 model could be reasonably used in the wastewater treatment plant simulation. The upgrade solution based on... [more]

A multistage centrifugal pump was developed for high head and high flow rate applications. A double-suction impeller and a twin-volute were installed at the first stage followed by an impeller, diffuser and return vanes for the next four stages. An initial design feasibility study was conducted using three-dimensional computational fluid dynamics tools to study the performance and the hydraulic losses associated with the design. Substantial losses in head and efficiency were observed at the interface between the first stage volute and the second stage impeller. An inducer-type guide vane (ITGV) was installed at this location to mitigate the losses by reducing the circumferential velocity of the fluid exiting the volute. The ITGV regulated the pre-swirl of the fluid entering the second stage impeller. The pump with and without ITGV is compared at the design flow rate. The pump with ITGV increased the stage head by 63.28% and stage efficiency by 47.17% at the second stage. As a result, t... [more]

A symmetrical NACA 0018 airfoil is often used in such applications as small-to-medium scale vertical-axis wind turbines and aerial vehicles. A review of the literature indicates a large gap in experimental studies of this airfoil at low and moderate Reynolds numbers in the previous century. This gap has limited the potential development of classical turbulence models, which in this range of Reynolds numbers predict the lift coefficients with insufficiently accurate results in comparison to contemporary experimental studies. Therefore, this paper validates the aerodynamic performance of the NACA 0018 airfoil and the characteristics of the laminar separation bubble formed on its suction side using the standard uncalibrated four-equation Transition SST turbulence model and the unsteady Reynolds-averaged Navier-Stokes (URANS) equations. A numerical study was conducted for the chord Reynolds number of 160,000, angles of attack between 0 and 11 degrees, as well as for the free-stream turbule... [more]

Three products hyrdrocyclone screen (TPHS) has been proposed for particle separation based on size. In TPHS, a cylindrical screen was embedded in a conventional hydrocyclone (CH) to combine the centrifugal classification and screening to particle separation based on size. The industrial application of TPHS indicates its better device performance than CH. Although, the earlier studies reveal some common understanding for TPHS, the information of the absent air column remains unknown. Hence, the combination of physical experiment and numerical simulation was considered involving a 75 mm TPHS for this knowledge gap. First, both the computational fluid dynamics (CFD) simulation with Reynolds stress mode and the physical experiment with a high-definition camera illustrate the development process details of a flow field in TPHS. That is, the water was imported along the tangential inlet into TPHS; then under the effects of the feed chamber wall and gravity, the liquid phase spiraled downward... [more]

A small Counter-Rotating Wind Turbine (CRWT) has been proposed and its performance has been investigated numerically. Results of a parametric study have been presented in this paper. As parameters, the axial distance between rotors and a tip speed ratio of each rotor have been selected. Performance parameters have been compared with reference to a Single Rotor Wind Turbine (SRWT). Simulations were carried out with Computational Fluids Dynamics (CFD) solver and a Large Eddy Scale approach to model turbulences. An Actuator Line Model has been chosen to represent rotors in the computational domain. Summing up the results of simulation tests, it can be stated that when constructing a CRWT turbine, rotors should be placed at a distance of at least 0.5 D (where D is rotor outer diameter) or more. One can then expect a noticeable power increase compared to a single rotor turbine. Placing the second rotor closer than 0.5 D guarantees a significant increase in power, but in such configurations,... [more]

A fire has been analyzed in a real underground coal mine, using a sublevel method, during an entire year. The study was focused on the collapsed area, reproducing a real mixture formed by coal, waste, and air gap. The analysis was done by means of an experimental analysis, a computational fuid dynamic model (CFD), and simulations using a mine ventilation software. Three scenarios were determined and studied regarding their influence on the evolution of the fire: (a) development of the fire without taking any action, (b) sealing off the affected areas, and (c) sealing and reducing the ventilation in the affected area and surrounding drifts. The study revealed the behavior of the fire in a real mine and the effectiveness of the main fire-fighting measures over time, verifying that none of the measures taken could eliminate the fire-induced in the collapsed area.

A free-piston engine is a novel internal combustion engine which has the advantages of a variable compression ratio and multi-fuel adaptability. This paper focuses on numerical simulation for combustion process and spray angle optimization of an opposed-piston free-piston engine. The working principle and spray-guided central combustor structure of the engine are discussed. A three-dimensional computational fluid dynamic model with moving mesh is presented based on the tested piston motion of the prototype. Calculation conditions, spray models, and combustion models were set-up according to the same prototype. The effects of spray angle on fuel evaporation rate, mixture distribution, heat release rate, in-cylinder pressure, in-cylinder temperature, and emissions were simulated and analyzed in detail. The research results indicate that the performance of the engine was very sensitive to the spray angle. The combustion efficiency and the indicated thermal efficiencies of 97.5% and 39.7%... [more]

Sustainable development can no longer neglect the growth of those technologies that look at the recovery of any energy waste in industrial processes. For example, in almost every industrial plant it happens that pressure energy is wasted in throttling devices for pressure and flow control needs. Clearly, the recovery of this wasted energy can be considered as an opportunity to reach not only a higher plant energy efficiency, but also the reduction of the plant Operating Expenditures (OpEx). In recent years, it is getting common to replace throttling valves with turbine-based systems (tuboexpander) thus getting both the pressure control and the energy recovery, for instance, producing electricity. However, the wide range of possible operating conditions, technical requirements and design constrains determine highly customized constructions of these turboexpanders. Furthermore, manufacturers are interested in tools enabling them to rapidly get the design of their products. For these reas... [more]

In the case of a severe accident, natural convection plays an important role in the atmosphere mixing of nuclear reactor containments. In this case, the natural convection might not in the steady-state condition. Hence, instead of steady-state simulation, the transient simulation should be performed to understand natural convection in the accident scenario within a nuclear reactor containment. The present study, therefore, was aimed at the transient 3-D numerical simulations of natural convection of air around a cylindrical containment with unsteady thermal boundary conditions (BCs) at the vessel wall. For this purpose, the experiment series was done in the CIGMA facility at Japan Atomic Energy Agency (JAEA). The upper vessel or both the upper vessel and the middle jacket was cooled by subcooled water, while the lower vessel was thermally insulated. A 3-D model was simulated with OpenFOAM®, applying the unsteady Reynolds-averaged Navier−Stokes equations (URANS) model. Different turbule... [more]

Wind turbines without pitch control are more preferable from economical point of view but aerodynamic stall affects them more and after a critical wind speed local boundary layer separation occurs. Consequently, their power production is relatively low. In this study, air ducts added on the blade and using the airflow from them the kinetic energy of the low-momentum fluid behind the surface was increased and delay of separation of the boundary layer from the surface was examined The Response Surface Optimization method was utilized in order to get the best possible design under the constraints and targets arranged for the parameters termed the diameter, slope, number and angle of attack of the air ducts. By using computational fluid dynamics analysis, optimum parameter values were obtained and air-ducted and air-duct free blade designs were compared. An improvement in power coefficient between 3.4−4.4% depending on wind speed was achieved with the new design. Due to increase in viscous... [more]

Semiconductor light sources are currently the fastest growing and most energy efficient group of light sources used in lighting technology. Their lighting parameters, such as luminous flux, correlated color temperature and color rendering index depend on the value of the forward current, as well as the temperature of the junction. LED source manufacturers usually specify, in data sheets, the effect of junction temperature and forward current on the luminous flux for individual light sources. The difficulty, however, is the correct determination of temperature and then lighting parameters, by simulation methods for multi-source lighting systems. Determining the junction temperature which affects lighting parameters is particulary important in the case of LED panels and luminaires, where thermally coupled LED sources shaping the output lighting parameters are in close proximity to each other. Additionally, other factors influencing the temperature distribution of sources, such as the des... [more]