Fluids, Volume 8, Issue 9 (September 2023) – 20 articles

In this work, a three-dimensional CFD model for the gas–solid flow of two different particle sizes in a CFB riser coupled with a kinetic theory (KT) has been developed. The properties of the solid phases are calculated using the proposed multi-particle kinetic theory. The CFD model is implemented in the commercial CFD software CFX4.4. In the current model, one gas phase and two solid phases are used. However, the model is generalised for one carrier phase and N number of solid phases to enable a realistic particle size distribution in the system. The momentum, volume fraction and granular temperature equations are solved for each individual solid phase and implemented into the CFD model through user-defined functions (UDFs). The k-ε turbulence model is used in simulating the circulating fluidised bed model. For verification, simulation results obtained with the new KT model were compared with experimental data, and then the model was used for further analysis. It was found that the proposed multi-particle model can be used to calculate the properties of gas–solid systems with particles of different sizes and/or densities, removing the assumptions of previous models that required all the particles to be of an equal mass, size and density. Full article

Valorizing agri-food industrial waste is essential for a circular economy, yielding high-value products, waste reduction, technological solutions, employment opportunities, and enhanced food security. This work shows the valorization of seeds generated as residues from the agri-food industries of guava pera (Psidium guajava) and Tommy Atkins mango (Mangifera indica L.), through extraction with supercritical carbon dioxide (scCO₂). After the optimization of the initial solid condition of the raw material (i.e., particle size and moisture content), scCO₂ pressure and temperature were established through the response surface methodology (RSM) to obtain an oily extract with the highest content in bioactive compounds of commercial relevance, as well as with a high antioxidant capacity. The total amount of oily extract in guava and mango seeds was 14% and 9%, respectively, while the maximum recovery of supercritical extract was 95% from guava seeds at 38 MPa and 50 °C, and 88% from mango seeds at 37 MPa and 63 °C. Bioactive fractions rich in squalene, γ-tocopherol, α-tocopherol, campesterol, β-sitosterol, and stigmasterol were obtained. The best supercritical extraction conditions, in terms of the bioactive fractions richest in minor compounds, were at 17 MPa and 50 °C for guava seeds and at 23 MPa and 63 °C for mango seeds. At these conditions, the highest antioxidant capacities were also found for the extracts. Thus, these bioactive fractions could be used in a variety of products in the cosmetic, food, pharmaceutical, and medical activities due to the beneficial properties of the identified compounds in health as antioxidants, anti-inflammatories, and cholesterol reducers. Full article

Marangoni bursting describes the spontaneous spread of a droplet of a binary mixture of alcohol/water deposited on a bath of oil, followed by its fast spontaneous fragmentation into a large number of smaller droplets in a self-similar way. Several papers have aimed to describe the physical phenomena underlying this spectacular phenomenon, in which two opposite effects, solutal and thermal Marangoni stresses, play competitive roles. We performed investigations of the Marangoni bursting phenomenon, paying attention to the surface temperature changes during bursting and after it. Fragmentation instabilities were monitored using a thermal camera for various initial alcohol/water compositions and at different stages of the process. We uncovered the role of thermocapillary Marangoni flows within the more viscous oil phase that are responsible for outward and inward shrinking of the periphery circle at the final stage of the phenomenon, enabling a more comprehensive understanding of the thermal Marangoni effect. Simulations of the Marangoni thermocapillary effect in an oil bath by solving coupled Navier–Stokes and heat transport equations using the COMSOL Multiphysics software platform support our experimental observations. Full article

Peristaltic pump technology is widely used wherever relatively low, highly accurately dosed volumetric flow rates are required and where fluid contamination must be excluded. Thus, typical fields of application include food, pharmaceuticals, medical technology, and analytics. In certain cases, when applied in conjunction with polymer-based tubing material, supplied peristaltic flow rates are reported to be significantly lower than the expected set flow rates. Said flow rate reductions are related to (i) the chosen tube material, (ii) tube material fatigue effects, and (iii) the applied pump frequency. This work presents a fast, dynamic, multiphysics, 1D peristaltic pump solver, which is demonstrated to capture all qualitatively relevant effects in terms of peristaltic flow rate reduction within linear peristaltic pumps. The numerical solver encompasses laminar fluid dynamics, geometric restrictions provided by peristaltic pump operation, as well as viscoelastic tube material properties and tube material fatigue effects. A variety of validation experiments were conducted within this work. The experiments point to the high degree of quantitative accuracy of the novel software and qualify it as the basis for elaborating an a priori drive correction. Full article

Highway buses are used in a wide range of commuting services and in the tourist industry. The demand for highway bus transportation has dramatically increased in the recent post-pandemic world, and airborne transmission risks may increase alongside the demand for highway buses, owing to a higher passenger density in bus cabins. We developed a numerical prediction method for the spatial distribution of airborne transmission risks inside bus cabins. For a computational fluid dynamics analyses, targeting two types of bus cabins, sophisticated geometries of bus cabins with realistic heating, ventilation, and air-conditioning were reproduced. The passengers in bus cabins were reproduced using computer-simulated persons. Airflow, heat, and moisture transfer analysis were conducted based on computational fluid dynamics, to predict the microclimate around passengers and the interaction between the cabin climate and passengers. Finally, droplet dispersion analysis using the Eulerian–Lagrangian method and an investigation of the spatial distribution of infection/spread risks, assuming SARS-CoV-2 infection, were performed. Through parametric analyses of passive and individual countermeasures to reduce airborne infection risks, the effectiveness of countermeasures for airborne infection was discussed. Partition installation as a passive countermeasure had an impact on the human microclimate, which decreased infection risks. The individual countermeasure, mask-wearing, almost completely prevented airborne infection. Full article

This study aims to investigate the impact of body heat loss on the thermal and aerodynamic conditions in a whole-body cryotherapy chamber. The underlying hypothesis is that the heat generated by the human body alters the thermal and aerodynamic environment inside the cabin. A numerical study was conducted to test this hypothesis and analyze the thermodynamic exchanges between the human body and the cabin during a 3 min whole-body cryotherapy session. The computational fluid dynamics (CFD) approach was used to study the unsteady heat transfer between the human body and the interior of the cryotherapy cabin. A thermal boundary condition, based on a mathematical model developed from experimental data, was applied to simulate skin cooling kinetics over time. The post-processing of the 3D results, including temperature, velocity fields, and thermal flux maps at the body surface, provided insight into the thermo-convective mechanisms involved in a whole-body cryotherapy session. The study found that body heat loss significantly affects the temperature fields inside the cabin, leading to global modifications of the aeraulic and thermal conditions. These findings suggest that cryotherapy protocols may need to be adjusted or the cabin set temperature optimized to enhance the therapeutic benefits. Full article

Partial electrification of hydraulic circuits to achieve energy savings requires an increase in the angular speed of the positive displacement pumps, with the risk of incomplete filling. In this context, the paper focuses on developing a computational fluid dynamics (CFD) model using SimericsMP+ for two external gear pumps, namely helical and spur type gears. The objective of this study is the analysis of the phenomena occurring on the suction side under conditions of incomplete filling at high speeds. Both CFD models have been validated by conducting experimental tests for measuring the flow rate delivered at various inlet pressures and angular speeds. The experimental results confirm the model’s capability to accurately detect the operating conditions at which the delivered flow rate starts to decrease due to the partial filling of the inter-teeth chambers. Furthermore, this paper investigates the effects of certain geometrical modifications to the spur gear pump. Specifically, the influence of the gear’s width-to-diameter ratio is studied, revealing that a lower ratio leads to slightly better filling. Conversely, increasing the inlet port diameter results in no improvement. Based on this study, the modelling approach appears to be accurate enough to serve as design tool for optimizing pumps to improve their filling capability. Full article

High-speed gas/vapour jets injected into a cross-moving sonic liquid signifies a vital phenomenon which bears useful applications in environmental and energy processes. In the present experimental study, a pulsating jet of supersonic steam was injected into cross-flowing water. Circulation zones of opposite vorticity owing to the interaction between the steam jet and cross-water flow were found. However, a large circulation appeared in front of the nozzle exit. Also, most small circulation regions were observed at higher water-flow rates (>2 m³/s). Among the prime mixing variables (i.e., turbulence kinetic energy (TKE) and Reynolds shear stress (RSS)), the RSS estimations backed a small diffusive phenomenon within a region far from the nozzle exit. Further information extracted from the PIV images indicated the existence of Kelvin–Helmholtz (KH) instabilities. The counter-rotating vortex pairs (CVPs) appeared to be significant in the region close to the nozzle exit, and they exhibited leeward side folds. Moreover, the effects of the operating conditions on the pressure recovery and mixing efficiency as well as the penetration and the separation height were evaluated to determine the optimisation of the phenomenon. By applying extreme difference analysis, the mixing efficiency was found as the most influential parameter. Full article

A novel device of flow rate augmentation is proposed and experimentally examined in a horizontal valveless closed loop pump using a time-dependent stenosis (convergent–divergent channel) in contrast with the commonly used taper tubes of constant opening as flow rectifiers. The stenosis, being a part of the flexible tube of the pump, is formed by a semi-cylindrical surface attached to a compression spring of adjustable pretension compressing the tube against a flat plate. Located at either side of the pump pincher, the shape of the stenosis changes in time, without any external power source, as a function of the fluid pressure and the pretension of the spring. The spring pretension is adjusted by a trial-and-error procedure aiming for net flow rate maximization for each pinching frequency. For the examined pitching frequencies (5 Hz to 11 Hz, for which net flow rate is maximized) and for compression ratios 38% to 75%, the maximum net flow rate was found to be 720% of the non-stenosis case. Important parameters for flow enhancement were found to be the stenosis location along the loop, its opening, the compression ratio at the pincher area and the pinching frequency. Full article

Traditionally, Fourier spectra have been employed to gain a deeper understanding of turbulence flow structures. The investigation of isotropic forced turbulence with passive scalars offers a straightforward means to examine the disparities between velocity and passive scalar spectra. This flow configuration has been extensively studied in the past, encompassing a range of Reynolds and Schmidt numbers. In this present study, direct numerical simulations (DNS) of this flow are conducted at sufficiently high Reynolds numbers, enabling the formation of a wide inertial range. The primary focus of this investigation is to quantitatively assess the variations in scalar spectra with the Schmidt number (Sc). The spectra exhibit a transition from a k^−5/3 scaling for low Sc to a k^−4/3 scaling for high Sc. The emergence of the latter power law becomes evident at Sc = 2, with its width expanding as Sc increases. To gain further insights into the underlying flow structures, a statistical analysis is performed by evaluating quantities aligned with the principal axes of the strain field. The study reveals that enstrophy is primarily influenced by the vorticity aligned with the intermediate principal strain axis, while the scalar gradient variance is predominantly controlled by the compressive strain. To provide a clearer understanding of the differences between enstrophy and scalar gradient variance, joint probability density functions (PDFs) and visualizations of the budget terms for both quantities are presented. These visualizations serve to elucidate the distinctions between the two and offer insights into their respective behaviors. Full article

We present a rare exact solution of the Navier–Stokes equations for the Hagen–Poiseuille flow through a quarter-elliptic tube. Utilizing the separation of variables method, we derive the solution and report expressions for both the volumetric flow rate and the friction factor–Reynolds number product. Full article

This paper proposes a CFD methodology for the simulation of the slipper’s dynamics of a swash-plate axial piston unit under actual operating conditions. The study considers a typical slipper design, including a vented groove at the swash-plate interface. The dynamic fluid–body interaction (DFBI) model is exploited to find the instantaneous position of the slipper, while the morphing approach is adopted to cope with the corresponding mesh distortion. A modular approach is adopted to ensure high-quality mesh on the entire slipper surface and sliding interfaces provide the fluid dynamic connection between neighboring regions. The external forces acting on the slipper are included by means of user-defined lookup tables with the simulation estimating the lift force induced by fluid compression. Moreover, the force produced by the metal-to-metal contact between the slipper and the swash plate is modeled through a specific tool of the software. The pressure signal over an entire revolution of the pump is taken as an input of the simulation and a variable time step is used to manage the high-pressure gradients occurring in the regions of inner and outer dead points of the piston. The weakly compressible characteristic of the fluid is considered by a specific pressure-dependent density approach, and the two-equation eddy-viscosity k-ω SST (shear stress transport) model is used to assess the turbulent behavior of the flow. Furthermore, the transitional model predicts the onset of transition, thus solving different equations depending on whether the flow enters a laminar or turbulent regime. In conclusion, the proposed methodology investigates the motion of the slipper in response to several external forces acting on the component. The numerical results are discussed in terms of variable clearance height, pressure distribution within the gap, and lift forces acting on the slipper under specific pump operations. Full article

The author’s published correlations for subcooled boiling in channels are further studied and developed in this work. The areas explored include choice of equivalent diameters for annuli and partially heated channels, effects of flow direction, micro-gravity, and orientation of heated surface. A new correlation is developed, which is a modification of the author’s earlier correlation. The author’s previous correlations and the new correlation are compared with a very wide range of test data for round tubes, rectangular channels, and annuli. Several other correlations are also compared with the same data. The authors’ correlations provide good agreement with data, the new correlation giving the least deviation. The data included hydraulic diameters from 0.176 to 22.8 mm, reduced pressure from 0.0046 to 0.922, subcooling from 0 to 165 K, mass flux from 59 to 31,500 kgm⁻²s⁻¹, all flow directions, and terrestial to micro gravity. The new correlation has mean absolute deviation (MAD) of 13.3% with 2270 data points from 49 sources. Correlations by others had MAD of 18% to 116%. The results are presented and discussed. Full article

The geometric complexity and high-pressure gradients that characterize the design of the flow field of gear pumps make it very difficult to obtain an accurate CFD simulation of the component. Usually, assumptions are made both in terms of geometrical features and physics being included in the analysis. The contact between the teeth, which is a key factor for the correct functioning of these pumps, represents a critical challenge in 3D CFD simulations, mainly due to the intrinsic limits of the dynamic meshing techniques that can hardly effectively manage a zero or close to zero gap point forming during gear rotation. The geometric complexity and high-pressure gradients that characterize the gear pump flow field make a CFD analysis quite difficult, and the contact between the gear teeth is usually avoided, thus being an extremely important feature. In this paper, a gear pump composed of inlet and outlet pipes was considered, and the contact between the gear’s teeth was modeled in two different ways, one where it is effectively implemented and one where it is avoided using distancing and a proper casing modification. Herein, a new methodology is proposed for the application of the dynamic mesh method in the Simcenter STAR-CCM+ environment using an adaptive remeshing technique. The proposed methodology is compared with the alternative overset meshing method available in the software. The new meshing method is implemented using a user-routing that reproduces the real geometry of the gears while rotating during the pump operation, with teeth contact included. The routine is optimized in order to limit the additional computation and time needed for the remeshing process. The results that can be obtained using the two meshing approaches for the gear pump are compared in terms of computational effort and the accuracy of the results. The two methods showed opposite results in almost all the reported results, with the overset being more precise in the radial pressure evaluation and the dynamic being more reliable in the cavitation/aeration extension cloud. Full article

Pectin is a versatile polysaccharide produced mainly from natural food sources and agro-industrial wastes, adding value to these by-products. For food applications, it is necessary that pectin first interacts with water for technical purposes. As a food additive, pectin acts as a solution thickener and gelling agent for food formulation, even in concentrations of less than 1 (g/100 mL or g/100 g), and it is sufficient to influence food products’ stability, rheology, texture, and sensory properties. Therefore, this review paper attempts to discuss the versability of pectin use, focusing on food application. It starts by showing the chemical structure, the sources’ potential, thickening, and gelling mechanisms and concludes by showing the main applications to the food sector and its rheological properties. Full article

This critical review delves into the impact of Computational Fluid Dynamics (CFD) modeling techniques, specifically 2D, 2.5D, and 3D simulations, on the performance and vortex dynamics of Darrieus turbines. The central aim is to dissect the disparities apparent in numerical outcomes derived from these simulation methodologies when assessing the power coefficient ($C_p$) within a defined velocity ratio ($\lambda$) range. The examination delves into the prevalent turbulence models shaping $C_p$ values, and offers insightful visual aids to expound upon their influence. Furthermore, the review underscores the predominant rationale behind the adoption of 2D CFD modeling, attributed to its computationally efficient nature vis-à-vis the more intricate 2.5D or 3D approaches, particularly when gauging the turbine’s performance within the designated $\lambda$ realm. Moreover, the study meticulously curates a compendium of findings from an expansive collection of over 250 published articles. These findings encapsulate the evolution of pivotal parameters, including $C_p$, moment coefficient ($C_m$), lift coefficient ($C_l$), and drag coefficient ($C_d$), as well as the intricate portrayal of velocity contours, pressure distributions, vorticity patterns, turbulent kinetic energy dynamics, streamlines, and Q-criterion analyses of vorticity. An additional focal point of the review revolves around the discernment of executing 2D parametric investigations to optimize Darrieus turbine efficacy. This practice persists despite the emergence of turbulent flow structures induced by geometric modifications. Notably, the limitations inherent to the 2D methodology are vividly exemplified through compelling CFD contour representations interspersed throughout the review. Vitally, the review underscores that gauging the accuracy and validation of CFD models based solely on the comparison of $C_p$ values against experimental data falls short. Instead, the validation of CFD models rests on time-averaged $C_p$ values, thereby underscoring the need to account for the intricate vortex patterns in the turbine’s wake—a facet that diverges significantly between 2D and 3D simulations. In a bid to showcase the extant disparities in CFD modeling of Darrieus turbine behavior and facilitate the selection of the most judicious CFD modeling approach, the review diligently presents and appraises outcomes from diverse research endeavors published across esteemed scientific journals. Full article

A three-dimensional, transient computational fluid dynamics analysis was conducted on an idealised geometry of a coronary artery fitted with representative geometries of an Absorb bioresorbable vascular scaffold (BVS) or a Xience drug-eluting stent (DES) in order to identify and compare areas of disturbed flow and potential risk sites. A non-Newtonian viscosity model was used with a transient velocity boundary condition programmed with user-defined functions. At-risk areas were quantified in terms of several parameters linked to restenosis: wall shear stress, time-averaged wall shear stress, oscillatory shear index, particle residence time, and shear rate. Results indicated that 71% of the BVS stented surface area had time-averaged wall shear stress values under 0.4 Pa compared to 45% of the DES area. Additionally, high particle residence times were present in 23% and 8% of the BVS and DES areas, respectively, with risk areas identified as being more prominent in close proximity to crowns and link struts. These results suggest an increased risk for thrombosis and neointimal hyperplasia for the BVS compared to the DES, which is in agreement with the outcomes of clinical trials. It is intended that the results of this study may be used as a pre-clinical tool to aid in the design of bioresorbable coronary stents. Full article

The formation of flow-induced, oriented structures in two-phase systems, as in this study, is a phenomenon of considerable interest to the scientific and industrial sectors. The main difficulty in understanding the formation of bands of droplets is the simultaneous interplay of physicochemical, hydrodynamic, and mechanical effects. Additionally, banded structure materials frequently show multiple length scales covering several decades as a result of complex time-dependent stress fields. Here, to facilitate understanding a subset of these structures, we studied water in oil emulsions and focused on the effects of three variables specifically: the confinement factor$(Co=2R/H)$, the viscosity ratio $\left(p\right),$ and the applied shear rate $\left(\dot{\gamma }\right)$. The confinement $\left(Co\right)$ is the ratio between the drop’s diameter ($2R$) and the separation of (the gap between) the circular rotating disks ($H$) containing the emulsion. We carried out (a) observations of the induced structure under different simple shear rates, as well as (b) statistical and morphological analysis of these bands. At low shear rates, the system self-assembles into bands along the direction of the flow and stacked normal to the velocity gradient direction. At higher shear rates is possible to observe bands normal to the vorticity direction. Here, we show that a detailed analysis of the dynamics of the band structures is amenable, as well as measurements of flow field anomalies simultaneously observed. The local emulsion viscosity varies in time, increasing in regions of higher droplet concentration and subsequently inducing velocity components perpendicular to the main flow direction. Thus, the emulsion morphology evolves and changes macroscopically. A relatively plausible explanation is attributed to the competitive effects of coalescence and the rupture of drops, where $p$ values less than one predominate coalescence. Full article

Accurate evaluation of thermo-fluid dynamic characteristics in tanks is critically important for designing liquid hydrogen tanks for small-scale hydrogen liquefiers to minimize heat leakage into the liquid and ullage. Due to the high costs, most future liquid hydrogen storage tank designs will have to rely on predictive computational models for minimizing pressurization and heat leakage. Therefore, in this study, to improve the storage efficiency of a small-scale hydrogen liquefier, a three-dimensional CFD model that can predict the boil-off rate and the thermo-fluid characteristics due to heat penetration has been developed. The prediction performance and accuracy of the CFD model was validated based on comparisons between its results and previous experimental data, and a good agreement was obtained. To evaluate the insulation performance of polyurethane foam with three different insulation thicknesses, the pressure changes and thermo-fluid characteristics in a partially liquid hydrogen tank, subject to fixed ambient temperature and wind velocity, were investigated numerically. It was confirmed that the numerical simulation results well describe not only the temporal variations in the thermal gradient due to coupling between the buoyance and convection, but also the buoyancy-driven turbulent flow characteristics inside liquid hydrogen storage tanks with different insulation thicknesses. In the future, the numerical model developed in this study will be used for optimizing the insulation systems of storage tanks for small-scale hydrogen liquefiers, which is a cost-effective and highly efficient approach. Full article

This comprehensive study focused on the standard taxi sign used in Ireland and its impact on drag production, fuel expenses, and CO${}_2$ emissions. Experimental analysis revealed that the conventional taxi sign significantly increased drag, especially when mounted on streamlined vehicles such as saloon cars, due to flow separation issues on the rear roof and rear windshield. Longitudinal reorientation of the sign offered a 14-fold reduction in drag increase compared to the traditional placement. It was found that positioning the sign in the middle of the roof offered the greatest fuel efficiency. Furthermore, the study estimated that implementing longitudinal repositioning on all Irish taxi signs could save drivers approximately EUR 832 per year and reduce national CO${}_2$ emissions by a substantial 22,464 tonnes annually. Comparative analyses with international taxi signs demonstrated that the Irish sign had significantly larger drag contributions, emphasizing the need for improved aerodynamics. To address the inherent drag issue, the study explored novel appendable devices and proposed alternative taxi sign designs. Among the tested solutions, a magnet-mounted front ramp proved the most effective, reducing total drag by nearly 30%. Additionally, a motorized flip-up taxi sign design demonstrated a remarkable 40% reduction in drag. Finally, a newly proposed taxi sign design, featuring longitudinal positioning and pointed triangular front and rear faces, resulted in a minimal 4.3% increase in vehicle drag compared to the baseline car. Full article