Fluids

This paper presents a comparison of two turbulence models implemented in two different frameworks (Eulerian and Lagrangian) in order to simulate the motion in calm water of a displacement hull. The hydrodynamic resistance is calculated using two open-source Computational Fluid Dynamics (CFD) software packages: OpenFOAM and DualSPHysics. These two packages are employed with two different numerical treatments to introduce turbulence closure effects. The methodology includes rigorous validation using a Wigley hull with experimental data taken from the literature. Then, the validated frameworks are applied to model a ship hull with a 30 m length overall (LOA), and their results discussed, outlining the advantages and disadvantages of the two turbulence treatments. In conclusion, the resistance calculated with OpenFOAM offers the best compactness of results and a shorter simulation time, whereas DualSPHysics can better capture the free-surface deformations, preserving similar accuracy. Full article

Meteorological factors, specifically wind direction and magnitude, influence the dispersion of atmospheric pollutants due to road traffic by affecting their spatial and temporal distribution. In this study, we are interested in the effect of the evolution of horizontal wind components, i.e., in the plane perpendicular to the altitude axis. A two-dimensional numerical model for solving the coupled traffic flow/pollution problem, whose pollutants are generated by vehicles, is developed. The numerical solution of this model is computed via an algorithm combining the characteristics method for temporal discretization with the finite-element method for spatial discretization. The numerical model is validated through a sensitivity study on the diffusion coefficient of road traffic and its impact on traffic density. The distribution of pollutant concentration, computed based on a source generated by traffic density, is presented for a single direction and different magnitudes of the wind velocity (stationary, Gaussian, linearly increasing and decreasing, sudden change over time), taking into account the stretching and tilting of plumes and patterns. The temporal evolution of pollutant concentration at various relevant locations in the domain is studied for two wind velocities (stationary and sudden change). Three regimes were observed for transport pollution depending on time and velocity: nonlinear growth, saturation, and decrease. Full article

This study investigates the morphology of a free-falling liquid jet by using a computational approach with an experimental validation. Numerical simulations are developed by means of the Finite Element Method (FEM) for solving the viscous fluid flow and the level set method in order to track the interface between the fluid and air. Experiments are conducted in order to capture the shape of a free-falling jet of viscous fluid via circular orifice, where the shape is measured optically. The numerical results are found to be in agreement with the experimental data, demonstrating the validity of the proposed approach. Furthermore, we analyze the role of the surface tension by implementing linear as well as nonlinear surface energy models. All computational codes are developed with the aid of open-source packages from FEniCS and made publicly available. The combination of experimental and numerical techniques provides a comprehensive understanding of the morphology of free-falling jets and may be extended to multiphysics problems rather in a straightforward manner. Full article

The hydrodynamic processes of landfill leachate migration in the base of a solid waste landfill can have a critical impact on the natural environment. In the case of improper operation of municipal solid waste placement facilities, highly contaminated leachate may penetrate into groundwater and subsequently into surface water. This work addresses fundamental issues of multicomponent fluid propagation in a multilayer porous medium, taking into account temperature inhomogeneities caused by waste decomposition with heat release. The regimes of diffusion and convection of leachate water penetrating into soil layers in the base of municipal solid waste facilities are numerically studied. Archival data from a set of field and laboratory measurements in the area of the operating landfill are used to model the features of pollutant propagation and determine migration parameters. The process of pollutant propagation and migration is described by quantitative values of dry residue content in leachate. Factors that have a significant impact on the migration of leachate are considered. The main ones are convective transfer, diffusion, and properties of the geological composition of the landfill base, which are taken into account in the mathematical formulation of the problem. The calculations show that leachate self-heating can substantially intensify leachate filtration and has to be taken into account in the assessment of leachate migration rate. Full article

A growing number of scholars are drawn to using numerical approaches powered by computer simulations as a potential solution to industrial problems. Replicating the compaction process in powder metallurgy with accuracy is one such issue. The Drucker-Prager-Cap model requires parameter calibration as the most used method for simulating powder compaction. This paper addresses this issue and presents a new technique for doing so. Utilizing Abaqus software 2020, the compaction process was simulated for the benchmark powder, which is the alloy Ag57.6-Cu22.4-Sn10-In10. The difference between simulation results and experimental data was reduced by applying the Particle Swarm Optimization technique in Python. The suggested approach may accurately forecast the Drucker-Prager-Cap model parameters, as demonstrated by comparing the optimized parameters utilizing the research’s method with their experimental values. The findings revealed how well the suggested approach in this study calibrated the DPC model, yielding three parameters—Young’s modulus, material cohesion, and hydrostatic pressure yield stress—with respective RMSEs of 1.95, 0.12, and 324.64 concerning their experimental values. Full article

Measurements of transmission as well as thermophysical properties have been carried out for different concentrations of SiO₂ nanoparticles (0, 1, 2, 5, 10, and 20 wt.%) in pure water (W) and ethylene glycol–water mixture (EG/W) in a weight ratio of 25/75, from 298 to 323 K at 100 kPa. In particular, the density, specific heat capacity, and thermal diffusivity are measured by a density meter, differential scanning calorimetry, and the laser flash method. In the case of 20 wt.% SiO₂, transmission in the visible range is reduced by 9.3%. Simultaneously, the density rises linearly to 12.3% (in W) and 11.3% (in EG/W). The specific heat capacity decreases to 15.9% (in W) and 17.3% (in EG/W), while the thermal diffusivity rises to 16.4% (in W) and 20.4% (in EG/W). While the density measurements are in very good agreement with the literature, the measured values of the specific heat capacity deviate more than 5%, especially for concentrations below 5 wt.% SiO₂. Moreover, it is shown that the thermal conductivity increases less for fluids in small gaps compared to other authors, which could be due to the suppression of the Brownian motion. Based on the measurement results, temperature- and concentration-dependent correlations for the investigated thermophysical properties are developed using two adjustable parameters. In general, these correlations show deviations of less than 4% from the experimental results, which will help to fill the gaps in the variation of experimental results due to size, SiO₂ nanofluid production, and different measurement devices, and thus optimize solar thermal applications with SiO₂ nanofluid. Full article

In this article, we address resistance forces and torques acting onto a body with arbitrary shape moving in an unsteady Stokes flow. We start analyzing the functional form of the expressions for forces and torques, which depend on the frequency parameter and on the position of the body in the domain of the fluid, and determining the asymptotic limits for high and low frequencies. In this regard, we show that, for high frequencies (hence short times), forces and torques are obtained by the associated hydrodynamic problems considering ideal potential flows, independently of the geometry of the problem. Afterwards, with the aim of obtaining expressions for forces and torques valid in the entire range of frequencies, we extend to the unsteady case the reflection method, largely employed in the theory of the steady Stokes flows. In this way, general expressions are provided in terms of the Faxén operators of the body and the Green function associated with the geometry of the confinement, that are valid, to the leading order, at any frequency, independently of the geometry of the problem. Finally, as the application of the general expressions, explicit relations for the resistance forces acting onto a spherical body with no-slip boundary conditions near a plane wall with full-slip boundary conditions are obtained, valid over the entire frequency range, provided that the distance between the plane and the sphere is larger than one sphere radius. Full article

This paper aims to numerically assess the cavitation damage of hydrodynamic machines and hydraulic components and its development in time, based on cavitation erosion tests with samples of used materials. The theoretical part of this paper is devoted to the numerical simulation of unsteady multiphase flow by means of computational fluid dynamics (CFD) and to the prediction of the erosive effects of the collapses of cavitation bubbles in the vicinity of solid surfaces. Compressible unsteady Reynolds-averaged Navier–Stokes equations (URANS) are solved together with the Zwart cavitation model. To describe the destructive collapses of vapor bubbles, the modeling of cavitation bubble dynamics along selected streamlines or trajectories is applied. The hybrid Euler–Lagrange approach with one-way coupling and the full Rayleigh–Plesset equation (R–P) are therefore utilized. This paper also describes the experimental apparatus with a rotating disc used to reach genuine hydrodynamic cavitation and conditions similar to those of hydrodynamic machines. The simulations are compared with the obtained experimental data, with good agreement. The proposed methodology enables the application of the results of erosion tests to the real geometry of hydraulic machines and to reliably predict the locations and magnitude of cavitation erosion, so as to select appropriate materials or material treatments for endangered parts. Full article

In this paper, we used Fourier Neural Operator (FNO) networks to solve reaction–diffusion equations. The FNO is a novel framework designed to solve partial differential equations by learning mappings between infinite-dimensional functional spaces. We applied the FNO to the Surface Quasi-Geostrophic (SQG) equation, and we tested the model with two significantly different initial conditions: Vortex Initial Conditions and Sinusoidal Initial Conditions. Furthermore, we explored the generalization ability of the model by evaluating its performance when trained on Vortex Initial Conditions and applied to Sinusoidal Initial Conditions. Additionally, we investigated the modes (frequency parameters) used during training, analyzing their impact on the experimental results, and we determined the most suitable modes for this study. Next, we conducted experiments on the number of convolutional layers. The results showed that the performance of the models did not differ significantly when using two, three, or four layers, with the performance of two or three layers even slightly surpassing that of four layers. However, as the number of layers increased to five, the performance improved significantly. Beyond 10 layers, overfitting became evident. Based on these observations, we selected the optimal number of layers to ensure the best model performance. Given the autoregressive nature of the FNO, we also applied it to solve the Gray–Scott (GS) model, analyzing the impact of different input time steps on the performance of the model during recursive solving. The results indicated that the FNO requires sufficient information to capture the long-term evolution of the equations. However, compared to traditional methods, the FNO offers a significant advantage by requiring almost no additional computation time when predicting with new initial conditions. Full article

This study presents a numerical evaluation of a Horizontal Rectilinear Earth–air Heat Exchanger (EAHE), considering the climatic and soil conditions of Viamão, Brazil, a subtropical region. The Constructal Design method, combined with the Exhaustive Search, was utilized to define the system constraints, degree of freedom, and performance indicators. The degree of freedom was characterized by the aspect ratio between the vertical and horizontal lengths of the elliptical cross-section duct (H/L). The performance indicators for the EAHE configurations were assessed based on thermal potential (TP) and pressure drop (PD). The Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) was applied for multi-objective evaluation, and a methodology for EAHE is proposed. The problem was solved using FLUENT software (version 2024 R2), which employs the Finite Volume Method to solve the conservation equations for mass, momentum, and energy. The (H/L)_T,o = 6.0 configuration showed a 16.4% increase in thermal performance for heating and 15.9% for cooling compared to the conventional circular duct. Conversely, the (H/L)_F,o = 1.0 configuration reduced pressure loss by 65.33%. The integration of Constructal Design with TOPSIS facilitated the identification of optimized geometries that achieve a balance between performance indicators and those that specifically prioritize thermal or fluid dynamic aspects, being this approach an original scientific contribution of the present work. Full article

An adaptive Free-Form Deformation parameterization method based on a spring analogy is presented for aerodynamic shape optimization problems. The proposed method effectively incorporates the gradients of the objective and constraint functions, achieving automatic control point adjustment based on variances in design variable components. To evaluate the performance of the adaptive FFD parameterization method, two 2D airfoil optimization design problems are examined. The optimization of the RAE2822 airfoil with 12, 18 and 24 design variables demonstrates superior results for the adaptive method compared to uniform parameterization. The adaptive method requires fewer iterations and achieves lower objective function values. Additionally, the optimization design from NACA0012 to RAE2822 airfoil with 18 design variables shows that the adaptive parameterization method achieves a lower drag coefficient while satisfying the optimization objective. This validates the method’s capability to finely adjust airfoil shapes and capture more optimal design points by exerting stronger control over local shapes. The proposed adaptive FFD parameterization method proves highly effective for optimizing aerodynamic shapes, offering stability and efficiency in the early stages of optimization, even with a limited number of design variables. Full article

We present a new formulation of the central moment lattice Boltzmann (LB) method based on a minimal continuous Fokker-Planck (FP) kinetic model, originally proposed for stochastic diffusive-drift processes (e.g., Brownian dynamics), by adapting it as a collision model for the continuous Boltzmann equation (CBE) for fluid dynamics. The FP collision model has several desirable properties, including its ability to preserve the quadratic nonlinearity of the CBE, unlike that based on the common Bhatnagar-Gross-Krook model. Rather than using an equivalent Langevin equation as a proxy, we construct our approach by directly matching the changes in different discrete central moments independently supported by the lattice under collision to those given by the CBE under the FP-guided collision model. This can be interpreted as a new path for the collision process in terms of the relaxation of the various central moments to “equilibria”, which we term as the Markovian central moment attractors that depend on the products of the adjacent lower order moments and a diffusion coefficient tensor, thereby involving of a chain of attractors; effectively, the latter are nonlinear functions of not only the hydrodynamic variables, but also the non-conserved moments; the relaxation rates are based on scaling the drift coefficient by the order of the moment involved. The construction of the method in terms of the relevant central moments rather than via the drift and diffusion of the distribution functions directly in the velocity space facilitates its numerical implementation and analysis. We show its consistency to the Navier-Stokes equations via a Chapman-Enskog analysis and elucidate the choice of the diffusion coefficient based on the second order moments in accurately representing flows at relatively low viscosities or high Reynolds numbers. We will demonstrate the accuracy and robustness of our new central moment FP-LB formulation, termed as the FPC-LBM, using the D3Q27 lattice for simulations of a variety of flows, including wall-bounded turbulent flows. We show that the FPC-LBM is more stable than other existing LB schemes based on central moments, while avoiding numerical hyperviscosity effects in flow simulations at relatively very low physical fluid viscosities through a refinement to a model founded on kinetic theory. Full article

This paper reviews recent biomechanical studies on heart valve flutter. The function of the heart valves is essential for maintaining effective blood circulation. Heart valve flutter is a kind of small vibration phenomenon like a flag fluttering in the wind, which is related to many factors such as a thrombus, valve calcification, regurgitation, and hemolysis and material fatigue. This vibration phenomenon is particularly prevalent in valve replacement patients. The biomechanical implications of flutter are profound and can lead to micro-trauma of valve tissue, accelerating its degeneration process and increasing the risk of thrombosis. We conducted a systematic review along with a critical appraisal of published studies on heart valve flutter. In this review, we summarize and analyze the existing literature; discuss the detection methods of frequency and amplitude of heart valve flutter, and its potential effects on valve function, such as thrombosis and valve degeneration; and discuss some possible ways to avoid flutter. These findings are important for optimizing valve design, diagnosing diseases, and developing treatment strategies. Full article

Motivated by existing techniques for implementing roughness on cylinders to control flow disturbances, we performed delayed detached eddy simulations (DDES) at Re = 6×106 that generated unsteady turbulent flow around a rectangular cylinder with a controlled wrinkled surface and a 1:4 aspect ratio. A systematic study of the roughness effect was carried out by implementing different configurations of equally spaced grooves and bumps on the top-surface of the cylinder. Our results suggest that groove geometries reduce energy dissipation at higher rates than the smooth reference case, whereas bumped cylinders produce relative pressures characterized by a sawtooth pattern along the middle-upper part of the cylinder. Moreover, cylinders with triangular bumps increase mean drag and lift forces by up to 8% and 0.08 units, respectively, while circular bumps increase vorticity and pressure disturbances on the wrinkled surface. All of these effects impact energy dissipation, vorticity, pressure coefficients, and flow velocity along the wrinkled surface. Both the surface-manufactured cylinders and the proposed visualization techniques could be replicated in a variety of engineering developments involving flow characterization in the presence of roughness. Full article

Many designs for natural-convection heat sinks and semi-empirical correlations have been proposed in the recent years, but they are only valid in a limited range of Elenbaas numbers El and were mostly tested for laminar flows. To alleviate those limits, parametric studies with 2D and quasi-3D models were carried out, in ranges of Grashof numbers up to 1.55×1011 and Elenbaas numbers up to 3.42×107. Ansys Fluent’s laminar, transition-SST, SST k-ω and k-ϵ models were applied. In addition, when used in this valid range, i.e., mean Elenbaas numbers, with the simplified quasi-3D model, the transition-SST model could predict better results, overestimating the heat flux by 10 to 15% compared to semi-empirical correlations. The 2D model was not deemed satisfying, regarding turbulence models. Consequently, a quasi-3D model was developed: it appeared to be an efficient trade-off between computational time and prediction accuracy, in particular for turbulence models. New grouping factors were also found, to ensure proper dimensioning of natural-convection heat sinks. They corresponded to non-dimensional parameters that dictated the physical behaviour of the heat sink with respect to the semi-empirical correlations. Typically, the ratio of the spacing to the optimal spacing predicted by Bar-Cohen’s correlation turned out to be an appropriate grouping factor with a threshold of 1, above which the fins could safely be considered as isolated, thus greatly simplifying all further calculations. Full article

This review provides a comprehensive examination of CFD modeling procedures for SCPP, with an emphasis on the detailed methodologies and a case study of the Manzanares prototype in Spain. The introduction delineates the historical context and physical modeling principles of solar chimneys, while highlighting their potential in industrial applications. The governing equations are meticulously discussed, covering assumptions in both 2D and 3D CFD modeling, the continuity and momentum equations, and the selection and accuracy of turbulence models, particularly the k-ε equations. The review also delves into heat transfer modeling, encompassing the energy equation and radiation modeling. Analytical evaluations of turbine pressure drop ratios and performance metrics for power generation efficiency are critically analyzed. The establishment of boundary conditions in solar chimney applications, including sky temperature assessments and distinctions between 2D and 3D boundary conditions, is extensively explored. Mesh generation techniques for both 2D and 3D CFD models are presented, supported by case studies. Parametric studies and experimental investigations are scrutinized to elucidate their impact on the performance of solar chimneys. The temperature–entropy diagram for an idealized Brayton cycle is introduced as a conceptual framework for efficiency analysis. Validation of the CFD codes, both 2D and 3D, against experimental data is performed to ensure model accuracy. The review further examines energy balance approaches in modeling solar chimneys, presenting state-of-the-art CFD results and discussing their implications in both 2D and 3D contexts. The synthesis of these findings culminates in a comprehensive conclusion, offering insights into the future directions and potential advancements in the CFD modeling of solar chimneys. This work aims to serve as a definitive reference for researchers and practitioners in the field, providing a robust foundation for the development and optimization of SCPP technology. Full article

This is a comprehensive overview on our research work to link interdisciplinary modeling and simulation techniques to improve the predictability and reliability simulations (PARs) of compressible turbulence with shock waves for general audiences who are not familiar with our nonlinear approach. This focused nonlinear approach is to integrate our “nonlinear dynamical approach” with our “newly developed high order entropy-conserving, momentum-conserving and kinetic energy-preserving methods” in the quantification of numerical uncertainty in highly nonlinear flow simulations. The central issue is that the solution space of discrete genuinely nonlinear systems is much larger than that of the corresponding genuinely nonlinear continuous systems, thus obtaining numerical solutions that might not be solutions of the continuous systems. Traditional uncertainty quantification (UQ) approaches in numerical simulations commonly employ linearized analysis that might not provide the true behavior of genuinely nonlinear physical fluid flows. Due to the rapid development of high-performance computing, the last two decades have been an era when computation is ahead of analysis and when very large-scale practical computations are increasingly used in poorly understood multiscale data-limited complex nonlinear physical problems and non-traditional fields. This is compounded by the fact that the numerical schemes used in production computational fluid dynamics (CFD) computer codes often do not take into consideration the genuinely nonlinear behavior of numerical methods for more realistic modeling and simulations. Often, the numerical methods used might have been developed for weakly nonlinear flow or different flow types other than the flow being investigated. In addition, some of these methods are not discretely physics-preserving (structure-preserving); this includes but is not limited to entropy-conserving, momentum-conserving and kinetic energy-preserving methods. Employing theories of nonlinear dynamics to guide the construction of more appropriate, stable and accurate numerical methods could help, e.g., (a) delineate solutions of the discretized counterparts but not solutions of the governing equations; (b) prevent numerical chaos or numerical “turbulence” leading to FALSE predication of transition to turbulence; (c) provide more reliable numerical simulations of nonlinear fluid dynamical systems, especially by direct numerical simulations (DNS), large eddy simulations (LES) and implicit large eddy simulations (ILES) simulations; and (d) prevent incorrect computed shock speeds for problems containing stiff nonlinear source terms, if present. For computation intensive turbulent flows, the desirable methods should also be efficient and exhibit scalable parallelism for current high-performance computing. Selected numerical examples to illustrate the genuinely nonlinear behavior of numerical methods and our integrated approach to improve PARs are included. Full article

Cavitation is a phase-change phenomenon from the liquid to the gas phase due to an increased flow velocity. As it causes severe erosion and noise, it is harmful to hydraulic machinery such as pumps, valves, and screw propellers. However, it can be utilized for water treatment, in chemical reactors, and as a mechanical surface treatment, as radicals and impacts at the point of cavitation bubble collapse can be utilized. Mechanical surface treatment using cavitation impacts is called “cavitation peening”. Cavitation peening causes less pollution because it uses water to treat the mechanical surface. In addition, cavitation peening improves on traditional methods in terms of fatigue strength and the working life of parts in the automobile, aerospace, and medical fields. As cavitation bubbles are utilized in cavitation peening, the study of cavitation bubbles has significant value in improving this new technique. To achieve this, many numerical analyses combined with field experiments have been carried out to measure the stress caused by bubble collapse and rebound, especially when collapse occurs near a solid boundary. Understanding the mechanics of bubble collapse can help to avoid unnecessary surface damage, enabling more accurate surface preparation, and improving the stability of cavitation peening. The present study introduces three cavitation bubble types: single, cloud, and vortex cavitation bubbles. In addition, the critical parameters, governing equations, and high-speed camera images of these three cavitation bubble types are introduced to support a broader understanding of the collapse mechanism and characteristics of cavitation bubbles. Then, the results of the numerical and experimental analyses of non-spherical cavitation bubbles are summarized. Full article

This article presents a comparison between the performance obtained by using a spatial discretization of the Euler equations based on a high-order discontinuous Galerkin (dG) method and different sets of variables. The sets of variables investigated are as follows: (1) conservative variables; (2) primitive variables based on pressure and temperature; (3) primitive variables based on the logarithms of pressure and temperature. The solution is advanced in time by using a linearly implicit high-order Rosenbrock-type scheme. The results obtained using the different sets are assessed across several canonical unsteady test cases, focusing on the accuracy, conservation properties and robustness of each discretization. In order to cover a wide range of physical flow conditions, the test-cases considered here are (1) the isentropic vortex convection, (2) the Kelvin–Helmholtz instability and (3) the Richtmyer–Meshkov instability. Full article