Search results

The effects of tube inclination and Grashof number on the fully developed hydrodynamic and thermal fields are investigated numerically for laminar ascending flow of air and water in uniformly heated circular tubes. The effects of the buoyancy induced secondary flow on the hydrodynamic and thermal fields are complex and strongly dependent on the Grashof number, the Prandtl number and the tube inclination. The influence of these parameters on the intensity of the secondary flow, on the distortion of the axial velocity profile and of the temperature field from the corresponding distributions for pure forced flow, as well as on the circumferential variation of the local shear stress and of the local Nusselt number are analysed. The average shear stress is higher than for pure forced flow and it increases with both the tube inclination and with the Grashof number. The average Nusselt number is higher than for pure forced flow and increases with the Grashof number. For a given fluid and Grashof number there exists an optimum tube inclination which maximizes the average Nusselt number. Correlations for the average Nusselt number in terms of Gr and Pr are presented for four different tube inclinations.

This study aims to assess the robustness and accuracy of the face-centred finite volume (FCFV) method for the simulation of compressible laminar flows in different regimes, using numerical benchmarks.

The work presents a detailed comparison with reference solutions published in the literature –when available– and numerical results computed using a commercial cell-centred finite volume software.

The FCFV scheme provides first-order accurate approximations of the viscous stress tensor and the heat flux, insensitively to cell distortion or stretching. The strategy demonstrates its efficiency in inviscid and viscous flows, for a wide range of Mach numbers, also in the incompressible limit. In purely inviscid flows, non-oscillatory approximations are obtained in the presence of shock waves. In the incompressible limit, accurate solutions are computed without pressure correction algorithms. The method shows its superior performance for viscous high Mach number flows, achieving physically admissible solutions without carbuncle effect and predictions of quantities of interest with errors below 5%.

The FCFV method accurately evaluates, for a wide range of compressible laminar flows, quantities of engineering interest, such as drag, lift and heat transfer coefficients, on unstructured meshes featuring distorted and highly stretched cells, with an aspect ratio up to ten thousand. The method is suitable to simulate industrial flows on complex geometries, relaxing the requirements on mesh quality introduced by existing finite volume solvers and alleviating the need for time-consuming manual procedures for mesh generation to be performed by specialised technicians.

This study is concerned with developing laminar flow of an incompressible, Newtonian fluid, having constant viscosity, rotating in circular and rectangular ducts that contain a 180° bend. The Reynolds number ranges from 100 to 400, the rotation number from 0 to 0.4, and the Dean number from 66 to 264. Positive and negative rotation modes are considered. The artificial compressibility method is used for the numerical calculations and new boundary conditions are developed for these flows. It is shown that rotation causes the secondary flow to occur in ducts of any geometry, and that the strength of the secondary flow in the bend due to both rotation and curvature decreases as compared to the no rotation case.

The purpose of this paper is to investigate a novel stabilization scheme to handle convection and pressure oscillation in the process of solving incompressible laminar flows by finite element method (FEM).

The semi-implicit stabilization scheme, characteristic-based polynomial pressure projection (CBP3) consists of the Characteristic-Galerkin method and polynomial pressure projection. Theoretically, the proposed scheme works for any type of element using equal-order approximation for velocity and pressure. In this work, linear 3-node triangular and 4-node tetrahedral elements are the focus, which are the simplest but most difficult elements for pressure stabilizations.

The present paper proposes a new scheme, which can stabilize FEM solution for flows of both low and relatively high Reynolds numbers. And the influence of stabilization parameters of the CBP3 scheme has also been investigated.

The research in this work is limited to the laminar incompressible flow.

The verification and validation of the CBP3 scheme are conducted by several 2 D and 3 D numerical examples. The scheme could be used to deal with more practical fluid problems.

The application of scheme to study complex hemodynamics of patient-specific abdominal aortic aneurysm is also presented, which demonstrates its potential to solve bio-flows.

The paper simulated 2 D and 3 D numerical examples with superior results compared to existing results and experiments. The novel CBP3 scheme is verified to be very effective in handling convection and pressure oscillation.

This study seeks to focus on the annular flow between rectangular and equilateral‐triangular ducts under all possible arrangements. The aim of this work is to obtain accurate prediction of the friction factor of this flow using high‐order finite element method.

Steady and fully developed laminar flow of incompressible Newtonian fluid in an annulus of variable cross‐sectional geometry is investigated numerically. Accurate prediction of the friction factor of this flow was obtained using high‐order finite element method.

The results were in agreement with already published findings in the literature. It was found that a higher annular area ratio will lead to a monotonic increase in fRe value in the case of regular annuli, and will lead to an increase followed by a decrease in fRe value in the case of irregular annuli. Also, it was, found that irregular annuli have lower fRe value than regular annuli, and that the square‐in‐triangle case has the lowest fRe value, whereas the square‐in‐square case has the highest fRe value.

Accurate prediction of the friction factor of the laminar flow in irregular annuli was obtained. Also, the obtained results can be utilized to optimize the annular geometries under consideration. In addition, the obtained results can lead to the design of more efficient heat exchangers.

The purpose of the paper is to propose some modifications to the SIMPLE (semi-implicit method for pressure-linked equations) algorithm. These modifications can ensure the numerical robustness and optimize computational efficiency. They remarkably promote the ability of the SIMPLE algorithm for incompressible DNS (direct numerical simulation) of multiscale problems, such as transitional flows and turbulent flows, by improving the properties of dispersion and dissipation.

The MDCD (minimized dispersion and controllable dissipation) scheme and MMIM (modified momentum interpolation method) are introduced. Six typical test cases are used to validate the modified algorithm, including the linear convective flow, lid-driven cavity flow, laminar boundary layer, Taylor vortex and DHIT (decaying homogenous isotropic turbulence). Particularly, a highly unsteady DNS of separated-flow transition in turbomachinery is precisely predicted by the modified algorithm.

The numerical examples show the distinct superiority of the modified algorithm in both internal flows and external flows. The advantages of the MDCD scheme and MMIM make the SIMPLE algorithm a promising method for DNS.

Some effective modifications to the SIMPLE algorithm are addressed. It is the first attempt to introduce the MDCD approach into the SIMPLE-type algorithms. The new algorithm is especially suitable for the incompressible DNS of convection-dominated flows.

The purpose of this paper is to present a novel numerical method to solve incompressible flows with natural and mixed convections using pseudo‐compressibility formulation.

The present method is derived using the framework of Harten Lax and van Leer with contact (HLLC) method of Toro, Spruce and Spears, that was originally developed for compressible gas dynamics equations. This work generalizes the algorithm described in the previous paper to the case where heat transfer is involved. Here, the solution of the Riemann problem is approximated by a three‐wave system.

A few test cases involving incompressible laminar flows inside 2D square cavity for various Rayleigh and Reynolds numbers are considered for validating the present method. The computed results from the present method are found to be quite promising.

Although pseudo‐compressibility formulation has been found to have superior performance and has the potential to have numerical treatments similar to compressible flow equations, only two numerical methods have been applied so far; namely Jameson method and Roes flux difference splitting method. A new sophisticated numerical method, following the framework of HLLC method, is derived and implemented for solving pseudo‐compressibility‐based incompressible flow equations with heat transfer.

To develop a high‐order compact finite‐difference method for solving flow problems containing shock waves.

A numerical algorithm based on high‐order compact finite‐difference schemes is developed for solving Navier‐Stokes equations in two‐dimensional space. The convective flux terms are discretized by using advection upstream splitting method (AUSM). The developed method is then used to compute some example laminar flow problems. The problems considered have a range of Mach number that corresponds to subsonic incompressible flow to hypersonic compressible flows that contain shock waves and shock/boundary‐layer interaction.

The paper shows that the AUSM flux splitting and high‐order compact finite‐difference methods can be used accurately and robustly in resolving shear layers and capturing shock waves. The highly diffusive nature of conventional flux splitting especially on coarse grids makes them inaccurate for boundary layers even with high‐order discretization.

This paper presents a high‐order numerical method that can accurately and robustly capture shock waves without deteriorating oscillations and resolve boundary layers and shock/boundary layer interaction.

This study aims to feature the application of the artificial compressibility method (ACM) for the numerical prediction of two-dimensional (2D) axisymmetric swirling flows.

The respective academic numerical solver, named IGal2D, is based on the axisymmetric Reynolds-averaged Navier–Stokes (RANS) equations, arranged in a pseudo-Cartesian form, enhanced by the addition of the circumferential momentum equation. Discretization of spatial derivative terms within the governing equations is performed via unstructured 2D grid layouts, with a node-centered finite-volume scheme. For the evaluation of inviscid fluxes, the upwind Roe’s approximate Riemann solver is applied, coupled with a higher-order accurate spatial reconstruction, whereas an element-based approach is used for the calculation of gradients required for the viscous ones. Time integration is succeeded through a second-order accurate four-stage Runge-Kutta method, adopting additionally a local time-stepping technique. Further acceleration, in terms of computational time, is achieved by using an agglomeration multigrid scheme, incorporating the full approximation scheme in a V-cycle process, within an efficient edge-based data structure.

A detailed validation of the proposed numerical methodology is performed by encountering both inviscid and viscous (laminar and turbulent) swirling flows with axial symmetry. IGal2D is compared against the commercial software ANSYS fluent – by using appropriate metrics and characteristic flow quantities – but also against experimental measurements, confirming the proposed methodology’s potential to predict such flows in terms of accuracy.

This study provides a robust methodology for the accurate prediction of swirling flows by combining the axisymmetric RANS equations with ACM. In addition, a detailed description of the convective flux Jacobian is provided, filling a respective gap in research literature.

Numerical solutions are presented for steady two‐dimensional symmetric flow past parabolic bodies in a uniform stream parallel to its axis. For this study, the full Navier‐Stokes equations and energy equation in parabolic coordinates were solved. A second order accurate finite difference scheme on a non‐uniform grid was used. The solution domain does not exclude the leading edge region as it is usually done with boundary layer flows. A wide range of Reynolds number (Re) is studied for different values of Prandtl number (Pr). It is found that the average Nusselt number (Nu) increases as (Pr) increases meanwhile, (Nu) decreases with the increase in (Re).