Military, space, and high-speed civilian applications will continue contributing to the renewed interest in compressible, high-speed turbulent boundary layers. To further complicate matters, these flows present complex computational challenges ranging from the pre-processing to the execution and subsequent post-processing of large-scale numerical simulations. Exploring more complex geometries at higher Reynolds numbers will demand scalable post-processing. Modern times have brought application developers and scientists the advent of increasingly more diversified and heterogeneous computing hardware, which significantly complicates the development of performance-portable applications. To address these challenges, we propose Aquila, a distributed, out-of-core, performance-portable post-processing library for large-scale simulations. It is designed to alleviate the burden of domain experts writing applications targeted at heterogeneous, high-performance computers with strong scaling performance. We provide two implementations, in C++ and Python; and demonstrate their strong scaling performance and ability to reach 60% of peak memory bandwidth and 98% of the peak filesystem bandwidth while operating out of core. We also present our approach to optimizing two-point correlations by exploiting symmetry in the Fourier space. A key distinction in the proposed design is the inclusion of an out-of-core data pre-fetcher to give the illusion of in-memory availability of files yielding up to 46% improvement in program runtime. Furthermore, we demonstrate a parallel efficiency greater than 70% for highly threaded workloads. Full article

In irrigation and drainage channels, vertical drops are generally used to transfer water from a higher elevation to a lower level. Downstream of these structures, measures are taken to prevent the destruction of the channel bed by the flow and reduce its destructive kinetic energy. In this study, the effect of use steps and grid dissipators on hydraulic characteristics regarding flow pattern, relative downstream depth, relative pool depth, and energy dissipation of a vertical drop was investigated by numerical simulation following the symmetry law. Two relative step heights and two grid dissipator cell sizes were used. The hydraulic model describes fully coupled three-dimensional flow with axial symmetry. For the simulation, critical depths ranging from 0.24 to 0.5 were considered. Values of low relative depth obtained from the numerical results are in satisfactory agreement with the laboratory data. The simultaneous use of step and grid dissipators increases the relative energy dissipation compared to a simple vertical drop and a vertical drop equipped with steps. By using the grid dissipators and the steps downstream of the vertical drop, the relative pool depth increases. Changing the pore size of the grid dissipators does not affect the relative depth of the pool. The simultaneous use of steps and grid dissipators reduces the downstream Froude number of the vertical drop from 3.83–5.20 to 1.46–2.00. Full article

The Augmented Lagrangian Smoothed Particle Hydrodynamics (ALSPH) method is a novel incompressible Smoothed Particle Hydrodynamics (SPH) approach that solves Navier–Stokes equations by an iterative augmented Lagrangian scheme through enforcing the divergence-free coupling of velocity and pressure fields. This study aims to systematically investigate the time step size and the number of inner iteration parameters to boost the performance of the ALSPH method. Additionally, the effects of computing spatial derivatives with two alternative schemes on the accuracy of numerical results are also scrutinized. Namely, the first scheme computes spatial derivatives on the updated particle positions at each iteration, whereas the second one employs the updated pressure and velocity fields on the initial particle positions to compute the gradients and divergences throughout the iterations. These two schemes are implemented to the solution of a flow over a circular cylinder at Reynolds numbers of 200 in two dimensions. Initially, simulations are performed in order to determine the optimum time step sizes by utilizing a maximum number of five iterations per time step. Subsequently, the optimum number of inner iterations is investigated by employing the predetermined optimum time step size under the same flow conditions. Finally, the schemes are tested on the same flow problem with different Reynolds numbers using the best performing combination of the aforementioned parameters. It is observed that the ALSPH method can enable one to increase the time step size without deteriorating the numerical accuracy as a consequence of imposing larger ALSPH penalty terms in larger time step sizes, which, overall, leads to improved computational efficiency. When considering the hydrodynamic flow characteristics, it can be stated that two spatial derivative schemes perform very similarly. However, the results indicate that the derivative operation with the updated particle positions produces slightly lower velocity divergence magnitudes at larger time step sizes. Full article

A 3D computational fluid dynamics method is used in the current study to investigate the hybrid nanofluid (HNF) flow and heat transfer in an annulus with hot and cold rods. The chief goal of the current study is to examine the influences of dissimilar Reynolds numbers, emissivity coefficients, and dissimilar volume fractions of nanoparticles on hydraulic and thermal characteristics of the studied annulus. In this way, the geometry is modeled using a symmetry scheme. The heat transfer fluid is a water, ethylene–glycol, or water/ethylene–glycol mixture-based Cu-Al₂O₃ HNF, which is a Newtonian NF. According to the findings for the model at Re = 3000 and ${\varphi }_1$ = 0.05, all studied cases with different base fluids have similar behavior. ϕ₁ and ϕ₂ are the volume concentration of Al₂O₃ and Cu nanoparticles, respectively. For all studied cases, the total average Nusselt number (Nu_ave) reduces firstly by an increment of the volume concentrations of Cu nanoparticles until ${\varphi }_2$ = 0.01 or 0.02 and then, the total Nu_ave rises by an increment of the volume concentrations of Cu nanoparticles. Additionally, for the case with water as the base fluid, the total Nu_ave at ${\varphi }_2$ = 0.05 is higher than the values at ${\varphi }_2$ = 0.00. On the other hand, for the other cases, the total Nu_ave at ${\varphi }_2$ = 0.05 is lower than the values at ${\varphi }_2$ = 0.00. For all studied cases, the case with water as the base fluid has the maximum Nu_ave. Plus, for the model at Re = 4000 and ${\varphi }_1$ = 0.05, all studied cases with different base fluids have similar behavior. For all studied cases, the total Nu_ave reduces firstly by an increment of the volume concentrations of Cu nanoparticles until ${\varphi }_2$ = 0.01 and then, the total Nu_ave rises by an increment of the volume concentrations of Cu nanoparticles. The Nu_ave augments are found by an increment of Reynolds numbers. Higher emissivity values should lead to higher radiation heat transfer, but the portion of radiative heat transfer in the studied annulus is low and therefore, has no observable increment in HNF flow and heat transfer. Full article

This paper reports a numerical study on heat transfer improvement in a double backward-facing expanding channel using different convectional fluids. A finite volume method with the k-ε standard model is used to investigate the effects of step, Reynolds number and type of liquid on heat transfer enhancement. Three types of conventional fluids (water, ammonia liquid and ethylene glycol) with Reynolds numbers varying from 98.5 to 512 and three cases for different step heights at a constant heat flux (q = 2000 W/m²) are examined. The top wall of the passage and the bottom wall of the upstream section are adiabatic, while the walls of both the first and second steps downstream are heated. The results show that the local Nusselt number rises with the augmentation of the Reynolds number, and the critical effects are seen in the entrance area of the first and second steps. The maximum average Nusselt number, which represents the thermal performance, can be seen clearly in case 1 for EG in comparison to water and ammonia. Due to the expanding of the passage, separation flow is generated, which causes a rapid increment in the local skin friction coefficient, especially at the first and second steps of the downstream section for water, ammonia liquid and EG. The maximum skin friction coefficient is detected in case 1 for water with Re = 512. Trends of velocities for positions (X/H1 = 2.01, X/H2 = 2.51) at the first and second steps for all the studied cases with different types of convectional fluids are indicated in this paper. The presented findings also include the contour of velocity, which shows the recirculation zones at the first and second steps to demonstrate the improvement in the thermal performance. Full article

A nonlocal fractional model of Brinkman type fluid (BTF) containing a hybrid nanostructure was examined. The magnetohydrodynamic (MHD) flow of the hybrid nanofluid was studied using the fractional calculus approach. Hybridized silver (Ag) and Titanium dioxide (TiO₂) nanoparticles were dissolved in base fluid water (H₂O) to form a hybrid nanofluid. The MHD free convection flow of the nanofluid (Ag-TiO₂-H₂O) was considered in a microchannel (flow with a bounded domain). The BTF model was generalized using a nonlocal Caputo-Fabrizio fractional operator (CFFO) without a singular kernel of order $\alpha$ with effective thermophysical properties. The governing equations of the model were subjected to physical initial and boundary conditions. The exact solutions for the nonlocal fractional model without a singular kernel were developed via the fractional Laplace transform technique. The fractional solutions were reduced to local solutions by limiting $\alpha \to 1$ . To understand the rheological behavior of the fluid, the obtained solutions were numerically computed and plotted on various graphs. Finally, the influence of pertinent parameters was physically studied. It was found that the solutions were general, reliable, realistic and fixable. For the fractional parameter, the velocity and temperature profiles showed a decreasing trend for a constant time. By setting the values of the fractional parameter, excellent agreement between the theoretical and experimental results could be attained. Full article

The sloshing behavior of systems is influenced by different factors related to the liquid level and tank specifications. Different approaches are applicable for the assessment of sloshing behavior in a tank. In this paper, a new numerical model based on the differential quadrature method and boundary element approaches is adopted to investigate the sloshing behavior of a tank with an elastic thin-walled beam. The model is developed based on small slope considerations of the free surface. The main assumption of fluid modeling is homogeneity, isotropy, inviscid, and only limited compressibility of the liquid. Indeed, the formulation is represented based on the reduced-order method and then is employed for simulating the coupling between structure and fluid in symmetric test cases. The results are verified with the ANSYS and literature for symmetric rigid structural walls and then the code is employed to study the behavior of fluid-structure interaction in a symmetric tank with new and efficient immersed structure. Full article

Micromixers are significant segments inside miniaturized scale biomedical frameworks. Numerical investigation of the effects of galloping cylinder characteristics inside a microchannel Newtonian, incompressible fluid in nonstationary condition is performed. Governing equations of the system include the continuity equation, and Navier–Stokes equations are solved within a moving mesh domain. The symmetry of laminar entering the channel is broken by the self-sustained motion of the cylinder. A parameter study on the amplitude and frequency of passive moving cylinder on the mixing of tiny particles in the fluid is performed. The results show a significant increase to the index of mixing uses of the galloping body in biomedical frameworks in the course of micro-electromechanical systems (MEMS) devices. Full article

In this article, the magnetohydrodynamic (MHD) flow of Casson nanofluid with thermal radiation over an unsteady shrinking surface is investigated. The equation of momentum is derived from the Navier–Stokes model for non-Newtonian fluid where components of the viscous terms are symmetric. The effect of Stefan blowing with partial slip conditions of velocity, concentration, and temperature on the velocity, concentration, and temperature distributions is also taken into account. The modeled equations of partial differential equations (PDEs) are transformed into the equivalent boundary value problems (BVPs) of ordinary differential equations (ODEs) by employing similarity transformations. These similarity transformations can be obtained by using symmetry analysis. The resultant BVPs are reduced into initial value problems (IVPs) by using the shooting method and then solved by using the fourth-order Runge–Kutta (RK) technique. The numerical results reveal that dual solutions exist in some ranges of different physical parameters such as unsteadiness and suction/injection parameters. The thickness of the velocity boundary layer is enhanced in the second solution by increasing the magnetic and velocity slip factor effect in the boundary layer. Increment in the Prandtl number and Brownian motion parameter is caused by a reduction of the thickness of the thermal boundary layer and temperature. Moreover, stability analysis performed by employing the three-stage Lobatto IIIA formula in the BVP4C solver with the help of MATLAB software reveals that only the first solution is stable and physically realizable. Full article

In this paper, the unsteady magnetohydrodynamic (MHD) flow of hybrid nanofluid (HNF) composed of $Cu-A{l}_2{O}_3$ /water in the presence of a thermal radiation effect over the stretching/shrinking sheet is investigated. Using similarity transformation, the governing partial differential equations (PDEs) are transformed into a system of ordinary differential equations (ODEs), which are then solved by using a shooting method. In order to validate the obtained numerical results, the comparison of the results with the published literature is made numerically as well as graphically and is found in good agreements. In addition, the effects of many emerging physical governing parameters on the profiles of velocity, temperature, skin friction coefficient, and heat transfer rate are demonstrated graphically and are elucidated theoretically. Based on the numerical results, dual solutions exist in a specific range of magnetic, suction, and unsteadiness parameters. It was also found that the values of $f^{″}\left(0\right)$ rise in the first solution and reduce in the second solution when the solid volume fraction ${\varphi }_{Cu}$ is increased. Finally, the temporal stability analysis of the solutions is conducted, and it is concluded that only the first solution is stable. Full article