Fluids

The research aims to design and construct a new mixed vertical boiler (fire tube – water tube) with three gas passes. The strength of this technological innovation is in the best use of the thermic transmission receiving fluid (hot water, steam, thermal oil), this due to its multipurpose function of three steps using alternative fuels (Diesel, Liquid Petroleum Gas LPG, natural gas), by improving the thermal efficiency of the boiler its temperature is reduced with gases at low temperatures, which in turn also reduce environmental pollution. The methodology focuses on calculating the transfer area with the calculation method that will allow dimensioning the boiler, considering the calculation of losses and the fluid speed, with two defined procedures, the first for fire tube and water tube boilers. And another alternative. The results obtained allowed optimizing the thermal efficiency level, achieving very significant thermal efficiency results: With LPG 92.4% for hot water and 92.42% to generate steam in the same way with natural gas 90.25% for hot water and 90.24% to generate steam as well with Diesel 2; 89.21% for hot water and 89.31% to generate steam. Full article

The measurement of the apparent contact angle on structured surfaces is much more difficult to obtain than on smooth surfaces because the pinning of liquid to the roughness has a tremendous influence on the three phase contact line. The results presented here clearly show an apparent contact angle variation along the three phase contact line. Accordingly, not only one value for the apparent contact angle can be provided, but a contact angle distribution or an interval has to be given to characterize the wetting behavior. For measuring the apparent contact angle distribution on regularly structured surfaces, namely micrometric pillars and grooves, an experimental approach is presented and the results are provided. A short introduction into the manufacturing process of such structured surfaces, which is a combination of Direct LASER Writing (DLW) lithography, electroforming and hot embossing shows the high quality standard of the used surfaces. Full article

In this study, a numerical modelling of thermal radiation and turbulent thermogravitational convection in a large-scale chamber containing a thermally-generating element is conducted. The lower border of the cabinet is maintained under adiabatic conditions, while on the other walls the convective boundary conditions (Robin boundary condition) are used. The managing equations with corresponding restrictions are transformed using the stream function–vorticity formulation and then solved by employing a finite difference method. The influence of both the height and wall emissivity of the heated source on fluid motion and the heat transmission in a large-scale chamber is investigated. Our results of the calculations on non-uniform grids with algebraic transformation are in excellent agreement with other available experimental and numerical outcomes for turbulent thermal convection in enclosures. The computations indicate that the average total Nusselt number is enhanced up to 2 times with an increase in the heater height. The results show that the surface emissivity of the heat source has a great influence on the total thermal transference coefficient. Furthermore, a growth of the heater surface emissivity has no significant effect on the flow structure. Full article

The dynamic behaviours of a filament in a viscoelastic uniform flow were investigated by an immersed boundary-lattice Boltzmann method. The effects of the Reynolds numbers ($Re$, ranging from 10 to 200) and the Weissenberg number ($Wi$, ranging from 0 to 1.2) on the filament flapping motion and the drag and lift coefficients on the filament were studied. It was found that a higher inertial effect (larger $Re$) promotes the flapping motion of the filament. In addition, the major effect of the viscoelasticity of the Giesekus fluid is to decrease the critical Reynolds number for the flapping motion of the filament and to promote the flapping motion. The drag coefficient on the filament in a Giesekus uniform flow decreases with the increase of $Wi$ at low $Re$ ($Re<100$), and experiences oscillations with similar amplitudes at all $Wi$ at a sufficiently high $Re$ ($Re>100$). In contrast, the viscoelasticity of the FENE-CR fluid increases the critical Reynolds number at lower $Wi$ ($Wi<0.8$), and shows little influence on the critical Reynolds number at higher $Wi$ ($Wi\ge 0.8$). In addition, the viscoelasticity of the FENE-CR fluid hinders the flapping motion of the filament, and increases the drag coefficient on the filament at low $Re$ ($Re<100$). Full article

The numerical investigation of the interpenetrating flow dynamics of a gas injected into a homogeneous porous media saturated with liquid is presented. The analysis is undertaken as a function of the inlet velocity, liquid–gas viscosity ratio (D) and physical properties of the porous medium, such as porous geometry and surface wettability. The study aims to improve understanding of the interaction between the physical parameters involved in complex multiphase flow in porous media (e.g., CO₂ sequestration in aquifers). The numerical simulation of a gaseous phase being introduced through a 2D porous medium constructed using seven staggered columns of either circular- or square-shaped micro-obstacles mimicking the solid walls of the pores is performed using the multiphase Lattice Boltzmann Method (LBM). The gas–liquid fingering phenomenon is triggered by a small geometrical asymmetry deliberately introduced in the first column of obstacles. Our study shows that the amount of gas penetration into the porous medium depends on surface wettability and on a set of parameters such as capillary number (C_a), liquid–gas viscosity ratio (D), pore geometry and surface wettability. The results demonstrate that increasing the capillary number and the surface wettability leads to an increase in the effective gas penetration rate, disregarding porous medium configuration, while increasing the viscosity ratio decreases the penetration rate, again disregarding porous medium configuration. Full article

Flexible propulsors are ubiquitous in aquatic and flying organisms and are of great interest for bioinspired engineering. However, many animal models, especially those found in the deep sea, remain inaccessible to direct observation in the laboratory. We address this challenge by conducting an integrative study of the giant larvacean, an invertebrate swimmer and “fluid pump” of the mesopelagic zone. We demonstrate a workflow involving deep sea robots, advanced imaging tools, and numerical modeling to assess the kinematics and resulting fluid transport of the larvacean’s beating tail. A computational model of the tail was developed to simulate the local fluid environment and the tail kinematics using embedded passive (elastic) and active (muscular) material properties. The model examines how varying the extent of muscular activation affects the resulting kinematics and fluid transport rates. We find that muscle activation in two-thirds of the tail’s length, which corresponds to the observed kinematics in giant larvaceans, generates a greater average downstream flow speed than other designs with the same power input. Our results suggest that the active and passive material properties of the larvacean tail are tuned to produce efficient fluid transport for swimming and feeding, as well as provide new insight into the role of flexibility in biological propulsors. Full article

The problem of Rayleigh–Bénard’s natural convection subjected to a temporally periodic cooling condition is solved numerically by the Lattice Boltzmann method with multiple relaxation time (LBM-MRT). The study finds its interest in the field of thermal comfort where current knowledge has gaps in the fundamental phenomena requiring their exploration. The Boussinesq approximation is considered in the resolution of the physical problem studied for a Rayleigh number taken in the range 10³ ≤ Ra ≤ 10⁶ with a Prandtl number equal to 0.71 (air as working fluid). The physical phenomenon is also controlled by the amplitude of periodic cooling where, for small values of the latter, the results obtained follow a periodic evolution around an average corresponding to the formulation at a constant cold temperature. When the heating amplitude increases, the physical phenomenon is disturbed, the stream functions become mainly multicellular and an aperiodic evolution is obtained for the heat transfer illustrated by the average Nusselt number. Full article

When using point measurement for environmental or sediment laden flows, there is well-recognised risk for not having aligned measurements that causes misinterpretation of the measured velocity data. In reality, these kinds of mismeasurement mainly happen due to the misinterpretation of bed orientation caused by the complexity of its determination in natural flows, especially in bedload laden or rough bed flows. This study proposes a novel bed realignment method to improve the measured data benchmarking by three-dimensional (3D) bed profile orientation and implemented it into different sets of experimental data. More specifically, the effects of realignment on velocity profile and streamwise turbulence structure measurements were investigated. The proposed technique was tested against experimental data collected over a water-worked and an experimentally arranged well-packed beds. Different from the well-packed rough bed, the water-worked bed has been generated after long sediment transport and settling and hence can be used to verify the proposed bed-alignment technique thoroughly. During the flow analysis, the corrected velocity, turbulence intensity and Reynolds stress profiles were compared to the theoretical logarithmic law, exponential law and linear gravity (universal Reynolds stress distribution) profiles, respectively. It has been observed that the proposed method has improved the agreement of the measured velocity and turbulence structure data with their actual theoretical profiles, particularly in the near-bed region (where the ratio of the flow measurement vertical distance to the total water depth, z/h, is limited to ≤0.4). Full article

The unsteady Ekman problem involves finding the response of the near-surface currents to wind stress forcing under linear dynamics. Its solution can be conveniently framed in the frequency domain in terms of a quantity that is known as the transfer function, the Fourier transform of the impulse response function. In this paper, a theoretical investigation of a fairly general transfer function form is undertaken with the goal of paving the way for future observational studies. Building on earlier work, we consider in detail the transfer function arising from a linearly-varying profile of the vertical eddy viscosity, subject to a no-slip lower boundary condition at a finite depth. The horizontal momentum equations, rendered linear by the assumption of horizontally uniform motion, are shown to transform to a modified Bessel’s equation for the transfer function. Two self-similarities, or rescalings that each effectively eliminate one independent variable, are identified, enabling the dependence of the transfer function on its parameters to be more readily assessed. A systematic investigation of asymptotic behaviors of the transfer function is then undertaken, yielding expressions appropriate for eighteen different regimes, and unifying the results from numerous earlier studies. A solution to a numerical overflow problem that arises in the computation of the transfer function is also found. All numerical code associated with this paper is distributed freely for use by the community. Full article

The influence of viscosity and surface tension on oxygen transfer was investigated using planar laser-induced fluorescence with inhibition (PLIF-I). The surface tension and the viscosity were modified using Triton X-100 and polyacrylamide, respectively. Changes in the hydrodynamic parameters of millimetric bubbles were identified, and transfer parameters were calculated. The results revealed a decrease in the mass transferred in the presence of a contaminant. For modified viscosity, the decrease in mass transferred was allowed for by current correlations, but the presence of surfactant led to a sharp decrease in the liquid side mass transfer coefficient, which became even lower when polymer was added. An explanation for the gap between classical correlations and experimental values of k_L is discussed, and a hypothesis of the existence of an accumulation of contaminant in the diffusion layer is proposed. This led to the possibility of a decrease in the diffusion coefficient and oxygen saturation concentration in the liquid film, explaining the discrepancy between models and experience. Adapted values of D_O2 and [O₂] * in this layer were estimated. This original study unravels the complexity of mass transfer from an air bubble in a complex medium. Full article

A short review of wavelet-based adaptive methods for modeling and simulation of incompressible turbulent flows is presented. Wavelet-based computational modeling approaches of different fidelities are recast into an integrated hierarchical adaptive eddy-capturing turbulence modeling framework. The wavelet threshold filtering procedure and the associated wavelet-filtered Navier–Stokes equations are briefly discussed, along with the adaptive wavelet collocation method that is used for numerical computations. Depending on the level of wavelet thresholding, the simulation is possibly supplemented with a localized closure model. The latest advancements in spatiotemporally varying wavelet thresholding procedures along with the adaptive-anisotropic wavelet-collocation method make the development of a fully adaptive approach feasible with potential applications for complex turbulent flows. Full article

The cross-flow over a surface-mounted elastic plate and its vibratory response are studied as a fundamental two-dimensional configuration to gain physical insight into the interaction of viscous flow with flexible structures. The governing equations are numerically solved on a deforming mesh using an arbitrary Lagrangian-Eulerian finite-element method. The turbulent flow is resolved using the unsteady Reynolds-averaged Navier–Stokes equations at a Reynolds number of $2.5\times {10}^4$ based on the plate height. The material properties of the plate are selected so that the structural frequency is close to the frequency of vortex shedding from the free edge of a rigid plate, which is studied initially as the reference case. The results show that the plate tip oscillates back and forth in response to unsteady fluid loading at twice the frequency of vortex shedding, which is attributable to the sequential formation of a primary vortex from the free edge and a secondary vortex near the base of the plate. The effects of the plate elasticity and density on the structural response are considered, and results are compiled in terms of the reduced velocity $U^{*}$ and the density ratio ${\rho }^{*}$. The standard deviation of tip displacement increases with reduced velocity in the range $7.1⩽U^{*}⩽18.4$, irrespective of whether the elasticity or the density of the plate is varied. However, the average deflection of the plate in the streamwise direction displays different scaling with $U^{*}$ and ${\rho }^{*}$, but scales almost linearly with the Cauchy number ∼$U^{*2}/{\rho }^{*}$. Interestingly, the synchronization between plate motion and vortex shedding ceases at $U^{*}=18.4$, and the excitation mechanism in the latter case resembles flutter instability, rather than vortex-induced vibration found at lower $U^{*}$. Full article

There is a lack of detailed information about wind flow and distribution of wind pressure around atypically shaped high-rise buildings. The national standard EN 1991-1-4 Eurocode 1 used to determine the effects of wind on the territory of Slovakia (and indeed other countries of the European Union) does not have a procedure for determining the effects of wind on objects of triangular shape. This presents a problem for designers and engineers, as there exist no generally binding/valid rules to follow when performing the wind effect analysis. This paper shows the procedure of identification and results of the external wind pressure coefficient for the triangularly shaped high-rise building. Two methods of calculation have been chosen for this purpose. First, experimental measurements were performed on a scaled model of the building cross-section in the wind tunnel. Subsequently, software simulations were performed on the same scaled model in the CFD (computational fluid dynamics) program ANSYS CFX. Results of wind pressure were obtained for two directions of wind flow measured in 16 sampling points distributed irregularly around the circumference of the model. Results were mutually compared and verified. At the end, the wind flow effects on a real-size triangular high-rise building in the built-up area performed by software simulation are shown. Full article

Many thermohydraulic issues about the safety of light water reactors are related to complicated two-phase flow phenomena. In these phenomena, computational fluid dynamics (CFD) analysis using the volume of fluid (VOF) method causes numerical diffusion generated by the first-order upwind scheme used in the convection term of the volume fraction equation. Thus, in this study, we focused on an interface compression (IC) method for such a VOF approach; this technique prevents numerical diffusion issues and maintains boundedness and conservation with negative diffusion. First, on a sufficiently high mesh resolution and without the IC method, the validation process was considered by comparing the amplitude growth of the interfacial wave between a two-dimensional gas sheet and a quiescent liquid using the linear theory. The disturbance growth rates were consistent with the linear theory, and the validation process was considered appropriate. Then, this validation process confirmed the effects of the IC method on numerical diffusion, and we derived the optimum value of the IC coefficient, which is the parameter that controls the numerical diffusion. Full article

Thermal hydraulics, in certain components of nuclear reactor systems, involve complex flow scenarios, such as flows assisted by free jets and stratified flows leading to turbulent mixing and thermal fluctuations. These complex flow patterns and thermal fluctuations can be extremely critical from a reactor safety standpoint. The component-level lumped approximations (0D) or one-dimensional approximations (1D) models for such components and subsystems in safety analysis codes cannot capture the physics accurately, and may introduce a large degree of modeling uncertainty. On the other hand, high-fidelity computational fluid dynamics codes, which provide numerical solutions to the Navier–Stokes equations, are accurate but computationally intensive, and thus cannot be used for system-wide analysis. An alternate way to improve reactor safety analysis is by building reduced-order emulators from computational fluid dynamics (CFD) codes to improve system scale models. One of the key challenges in developing a reduced-order emulator is to preserve turbulent mixing and thermal fluctuations across different-length scales or time-scales. This paper presents the development of a reduced-order, non-linear, “Markovian” statistical surrogate for turbulent mixing and scalar transport. The method and its implementation are demonstrated on a canonical problem of differentially heated channel flow, and high-resolution direct numerical simulations (DNS) data are used for emulator or surrogate development. This statistical surrogate model relies on Kramers–Moyal expansion and emulates the turbulent velocity signal with a high degree of accuracy. Full article

The ability to accurately predict the dynamics of fast moving and deforming interfaces is of interest to a number of applications including ink printing, drug delivery and fuel injection. In the current work we present a new compressible framework within OpenFOAM which incorporates mitigation strategies for the well known issue of spurious currents. The framework incorporates the compressible algebraic Volume-of-Fluid (VoF) method with additional interfacial treatment techniques including volume fraction smoothing and sharpening (for the calculation of the interface geometries and surface tension force, respectively) as well as filtering of the capillary forces. The framework is tested against different benchmarks: A 2D stationary droplet, a high velocity impact droplet case (500 m/s impact velocity) against a dry substrate and, with the same impact conditions, against a liquid film. For the 2D static droplet case, our results are consistent with what is observed in the literature when these strategies are implemented within incompressible frameworks. For the high impact droplet cases we find that accounting for both compressibility and correct representation of the interface is very important in numerical simulations, since pressure waves develop and propagate within the droplet interacting with the interface. While the implemented strategies do not alter the dynamics of the impact and the droplet shape, they have a considerable effect on the lamella formation. Our numerical method, although currently implemented for droplet cases, can also be used for any fast moving interface with or without considering the impact on a surface. Full article

If a gas volume is distributed into many microbubbles of a sub-millimetre size, the total gas/liquid surface becomes very large. This increases overall heat and/or mass transport across the sum of surfaces. The paper discusses several applications in which the use of microbubbles increases efficiency of various processes, especially in wastewater treatment and in growing microorganisms such as algae, yeast, bacteria, or primitive fungi. The problem of microbubble generation by percolation in aerator is their coalescence into larger bubbles, whatever small are the pores in the aerator in which the microbubbles are generated. The solution of this size discrepancy question was found in agitating the gas flow by a fluidic oscillator prior to its injection through the aerator. The oscillator is a no-moving-part device, simple, inexpensive, resistant to external effects like acceleration or heat, and with long maintenance-free working life. Full article

Coronavirus disease 2019 (Covid-19) has been identified as being transmitted among humans with droplets from breath, cough, and sneezes. Understanding the droplets’ behavior can be critical information to avoid disease transmission, especially while designing a device deals with human air respiratory. Although various studies have provided enormous computational fluid simulations, most cases are too specific and quite challenging to combine with other similar studies directly. Therefore, this paper proposes a systematic approach to predict the droplet behavior for coughing cases using machine learning. The approach consists of three models, which are droplet generator, mask model, and free droplet model modeled using feedforward neural network (FFNN). The evaluation has shown that the three FFNNs models’ accuracies are relatively high, with R-values of more than 0.990. The model has successfully predicted the evaporation effect on the diameter reduction and the completely evaporated state, which can be considered unlearned cases for machine learning models. The predicted horizontal distance pattern also agrees with the data in the literature. In summary, the proposed approach has demonstrated the capability to predict the diameter pattern according to the experimental or previous work data at various mask face types. Full article

In this paper, we describe and evaluate an algebraic model that has been adopted in a diagnostic wind flow model. Our numerical model is based on the Röckle’s wind modelling approach and we intend to reproduce the steady-state flow patterns of recirculation vortices that are generated in the near-wake region behind a vehicle. The evaluation of the practicality of our proposed model is performed by comparing the wind tunnel experiments of the flow around a vehicle conducted by the Loughborough University. We also compare the numerical results of our model with the CFD model. The proposed model reproduces the flow patterns behind a vehicle and it has significant advantages, such as low numerical costs. We expect that further improvements in the algebraic model when considering the vehicle’s shape will improve its practicality for the numerical analysis of flow fields around a vehicle. Full article