Search results

This study aims to analyse the non-linear losses of a porous media (stack) composed by parallel plates and inserted in a resonator tube in oscillatory flows by proposing numerical correlations between pressure gradient and velocity.

The numerical correlations origin from computational fluid dynamics simulations, conducted at the microscopic scale, in which three fluid channels representing the porous media are taken into account. More specifically, for a specific frequency and stack porosity, the oscillating pressure input is varied, and the velocity and the pressure-drop are post-processed in the frequency domain (Fast Fourier Transform analysis).

It emerges that the viscous component of pressure drop follows a quadratic trend with respect to velocity inside the stack, while the inertial component is linear also at high-velocity regimes. Furthermore, the non-linear coefficient b of the correlation ax + bx² (related to the Forchheimer coefficient) is discovered to be dependent on frequency. The largest value of the b is found at low frequencies as the fluid particle displacement is comparable to the stack length. Furthermore, the lower the porosity the higher the Forchheimer term because the velocity gradients at the stack geometrical discontinuities are more pronounced.

The main novelty of this work is that, for the first time, non-linear losses of a parallel plate stack are investigated from a macroscopic point of view and summarised into a non-linear correlation, similar to the steady-state and well-known Darcy–Forchheimer law. The main difference is that it considers the frequency dependence of both Darcy and Forchheimer terms. The results can be used to enhance the analysis and design of thermoacoustic devices, which use the kind of stacks studied in the present work.

This study aims to computational numerical simulations to clarify and explore the influences of periodic cellular lattice (PCL) morphological parameters – such as lattice structure topology (simple cubic, body-centered cubic, z-reinforced body-centered cubic [BCCZ], face-centered cubic and z-reinforced face-centered cubic [FCCZ] lattice structures) and porosity value ( ) – on the thermal-hydraulic characteristics of the novel trussed fin-and-elliptical tube heat exchanger (FETHX), which has led to a deeper understanding of the superior heat transfer enhancement ability of the PCL structure.

A three-dimensional computational fluid dynamics (CFD) model is proposed in this paper to provide better understanding of the fluid flow and heat transfer behavior of the PCL structures in the trussed FETHXs associated with different structure topologies and high-porosities. The flow governing equations of the trussed FETHX are solved by the CFD software ANSYS CFX® and use the Menter SST turbulence model to accurately predict flow characteristics in the fluid flow region.

The thermal-hydraulic performance benchmarks analysis – such as field synergy performance and performance evaluation criteria – conducted during this research successfully identified demonstrates that if the high porosity of all PCL structures decrease to 92%, the best thermal-hydraulic performance is provided. Overall, according to the obtained outcomes, the trussed FETHX with the advantages of using BCCZ lattice structure at 92% porosity presents good thermal-hydraulic performance enhancement among all the investigated PCL structures.

To the best of the authors’ knowledge, this paper is one of the first in the literature that provides thorough thermal-hydraulic characteristics of a novel trussed FETHX with high-porosity PCL structures.

The purpose of this paper is to present a procedure to design an experimental setup meant to validate an innovative approach for simulating, via computational fluid dynamics, a high-pressure gas release from a rupture (e.g. on an offshore oil and gas platform). The design is based on a series of scaling exercises, some of which are anything but trivial.

The experimental setup is composed of a wind tunnel, the instrumented scaled (1:10) mock-up of an offshore platform and a gas release system. A correct scaling approach is necessary to define the reference speed in the wind tunnel and the conditions of the gas release to maintain similarity with respect to the real-size phenomena. The scaling of the wind velocity and the scaling of the gas release were inspired by the approach proposed by Hall et al. (1997): a dimensionless group was chosen to link release parameters, wind velocity and geometric scaling factor.

The theoretical scaling approaches for each different part of the setup were applied to the design of the experiment and some criticalities were identified, such as the existence of a set of case studies with some release parameters laying outside the applicability range of the developed scaling methodology, which will be further discussed.

The resulting procedure is one of a kind because it involves a multi-scaling approach because of the different aspects of the design. Literature supports for the different scaling theories but, to the best of the authors’ knowledge, fails to provide an integrated approach that considers the combined effects of scaling.

Efficient thermal management of central processing unit (CPU) cooling systems is vital in the context of advancing information technology and the demand for enhanced data processing speeds. This study aims to explore the thermal performance of a CPU cooling setup using a cylindrical porous metal foam heat sink.

Nanofluid flow through the metal foam is simulated using the Darcy–Brinkman–Forschheimer equation, accounting for magnetic field effects. The temperature distribution is modeled through the local thermal equilibrium equation, considering viscous dissipation. The problem’s governing partial differential equations are solved using the similarity method. The CPU’s hot surface serves as a solid wall, with nanofluid entering the heat sink as an impinging jet. Verification of the numerical results involves comparison with existing research, demonstrating strong agreement across numerical, analytical and experimental findings. Ansys Fluent® software is used to assess temperature, velocity and streamlines, yielding satisfactory results from an engineering standpoint.

Investigating critical parameters such as Darcy number (10⁻⁴ ≤ Da_D ≤ 10⁻²), aspect ratio (0.5 ≤ H/D ≤ 1.5), Reynolds number (5 ≤ Re_D,bf ≤ 3500), Eckert number (0 ≤ EC_bf ≤ 0.1) , porosity (0.85 ≤ ε ≤ 0.95), Hartmann number (0 ≤ Ha_D,bf ≤ 300) and the volume fraction of nanofluid (0 ≤ φ ≤ 0.1) reveals their impact on fluid flow and heat sink performance. Notably, Nusselt number will reduce 45%, rise 19.2%, decrease 14.1%, and decrease 0.15% for Reynolds numbers of 600, with rising porosity from 0.85 to 0.95, Darcy numbers from 10⁻⁴ to 10⁻², Eckert numbers from 0 to 0.1, and Hartman numbers from 0 to 300.

Despite notable progress in studying thermal management in CPU cooling systems using porous media and nanofluids, there are still significant gaps in the existing literature. First, few studies have considered the Darcy–Brinkman–Forchheimer equation, which accounts for non-Darcy effects and the flow and geometric interactions between coolant and porous medium. The influence of viscous dissipation on heat transfer in this specific geometry has also been largely overlooked. Additionally, while nanofluids and impinging jets have demonstrated potential in enhancing thermal performance, their utilization within porous media remains underexplored. Furthermore, the unique thermal and structural characteristics of porous media, along with the incorporation of a magnetic field, have not been fully investigated in this particular configuration. Consequently, this study aims to address these literature gaps and introduce novel advancements in analytical modeling, non-Darcy flow, viscous dissipation, nanofluid utilization, impinging jets, porous media characteristics and the impact of a magnetic field. These contributions hold promising prospects for improving CPU cooling system thermal management and have broader implications across various applications in the field.

The purpose of this paper is to present a fast and bare bones implementation of a numerical method for quickly simulating turbulent thermal flows on GPUs. The work also validates earlier research showing that the lattice Boltzmann method (LBM) method is suitable for complex thermal flows.

A dual lattice hydrodynamic (D3Q27) thermal (D3Q7) multiple-relaxation time LBM model capable of thermal DNS calculations is implemented in CUDA.

The model has the same computational performance compared to earlier publications of similar LBM solvers. The solver is validated against three benchmark cases for turbulent thermal flow with available data and is shown to be in excellent agreement.

The combination of a D3Q27 and D3Q7 stencil for a multiple relaxation time -LBM has, to the authors’ knowledge, not been used for simulations of thermal flows. The code is made available in a public repository under a free license.

Endwall film cooling protects vane endwall by coolant coverage, especially at the leading edge (LE) region and vane-pressure side (PS) junction region. Strong flow impingement and complex vortexaa structures on the vane endwall cause difficulties for coolant flows to cover properly. This work aims at a full-scale arrangement of film cooling holes on the endwall which improves coolant efficiency in the LE region and vane-PS junction region.

The endwall film holes are grouped in four-holes constructal patterns. Three ways of arranging the groups are studied: based on the pressure field, the streamlines or the heat transfer field. The computational analysis is done with the k-ω SST model after validating the turbulence model properly.

By clustering the film cooling holes in four-holes patterns, the ejection of the coolant flow is stronger. The four-holes constructal patterns also improve the local coolant coverage in the “tough” regions, such as the junction region of the PS and the endwall. The arrangement based on streamlines distribution can effectively improve the coolant coverage and the arrangement based on the heat transfer distribution (HTD) has benefits by reducing high-temperature regions on the endwall.

A full-scale endwall film cooling design is presented considering interactions of different film cooling holes. A comprehensive model validation and mesh independence study are provided. The cooling holes pattern on the endwall is designed as four-holes constructal patterns combined with several arrangement choices, i.e. by pressure, by heat transfer and by streamline distributions.

Reynolds-averaged Navier–Stokes (RANS) models often perform poorly in shock/turbulence interaction regions, resulting in excessive wall heat load and incorrect representation of the separation length in shockwave/turbulent boundary layer interactions. The authors suggest that this can be traced back to inadequate numerical treatment of the inviscid fluxes. The purpose of this study is an extension to the well-known Harten, Lax, van Leer, Einfeldt (HLLE) Riemann solver to overcome this issue.

It explicitly takes into account the broadening of waves due to the averaging procedure, which adds numerical dissipation and reduces excessive turbulence production across shocks. The scheme is derived based on the HLLE equations, and it is tested against three numerical experiments.

Sod’s shock tube case shows that the scheme succeeds in reducing turbulence amplification across shocks. A shock-free turbulent flat plate boundary layer indicates that smooth flow at moderate turbulence intensity is largely unaffected by the scheme. A shock/turbulent boundary layer interaction case with higher turbulence intensity shows that the added numerical dissipation can, however, impair the wall heat flux distribution.

The proposed scheme is motivated by implicit large eddy simulations that use numerical dissipation as subgrid-scale model. Introducing physical aspects of turbulence into the numerical treatment for RANS simulations is a novel approach.

This study aims to develop an experimentally validated three-dimensional numerical model for predicting different flow patterns produced with a gas dynamic virtual nozzle (GDVN).

The physical model is posed in the mixture formulation and copes with the unsteady, incompressible, isothermal, Newtonian, low turbulent two-phase flow. The computational fluid dynamics numerical solution is based on the half-space finite volume discretisation. The geo-reconstruct volume-of-fluid scheme tracks the interphase boundary between the gas and the liquid. To ensure numerical stability in the transition regime and adequately account for turbulent behaviour, the k-ω shear stress transport turbulence model is used. The model is validated by comparison with the experimental measurements on a vertical, downward-positioned GDVN configuration. Three different combinations of air and water volumetric flow rates have been solved numerically in the range of Reynolds numbers for airflow 1,009–2,596 and water 61–133, respectively, at Weber numbers 1.2–6.2.

The half-space symmetry allows the numerical reconstruction of the dripping, jetting and indication of the whipping mode. The kinetic energy transfer from the gas to the liquid is analysed, and locations with locally increased gas kinetic energy are observed. The calculated jet shapes reasonably well match the experimentally obtained high-speed camera videos.

The model is used for the virtual studies of new GDVN nozzle designs and optimisation of their operation.

To the best of the authors’ knowledge, the developed model numerically reconstructs all three GDVN flow regimes for the first time.

Coronary artery disease (CAD) is reported as one of the most common sources of death all over the world. The presence of stenosis (plaque) in the coronary arteries results in the restriction of blood supply, leading to myocardial infarction. The current study investigates the influence of multi stenosis on hemodynamic properties in a patient-specific left coronary artery.

A three-dimensional model of the patient-specific left coronary artery was reconstructed based on computed tomography (CT) scan images using MIMICS-20 software. The diseased model of the left coronary artery was investigated, having the narrowing of 90% and 70% of area stenosis (AS) at the left anterior descending (LAD) and left circumflex (LCX), respectively.

The results indicate that the upstream region of stenosis experiences very high pressure for 90% AS during the systolic period of the cardiac cycle. The pressure drops maximum as the flow travels into the stenotic zone, and the high flow velocities were observed across the 90% AS. The higher wall shear stresses occur at the stenosis region, and it increases with the increase in the flow rate. It is found that the maximum wall shear stress across 90% AS is at the highest risk for rupture. A recirculation region immediately after the stenosis results in the further development of stenosis.

The current study provides evidence that there is a strong effect of multi-stenosis on the blood flow in the left coronary artery.

This paper aims to comprehensively clarify the research status of thermal transport of supercritical aviation kerosene, with particular interests in the effect of cracking on heat transfer.

A brief review of current research on supercritical aviation kerosene is presented in views of the surrogate model of hydrocarbon fuels, chemical cracking mechanism of hydrocarbon fuels, thermo-physical properties of hydrocarbon fuels, turbulence models, flow characteristics and thermal performances, which indicates that more efforts need to be directed into these topics. Therefore, supercritical thermal transport of n-decane is then computationally investigated in the condition of thermal pyrolysis, while the ASPEN HYSYS gives the properties of n-decane and pyrolysis products. In addition, the one-step chemical cracking mechanism and SST k-ω turbulence model are applied with relatively high precision.

The existing surrogate models of aviation kerosene are limited to a specific scope of application and their thermo-physical properties deviate from the experimental data. The turbulence models used to implement numerical simulation should be studied to further improve the prediction accuracy. The thermal-induced acceleration is driven by the drastic density change, which is caused by the production of small molecules. The wall temperature of the combustion chamber can be effectively reduced by this behavior, i.e. the phenomenon of heat transfer deterioration can be attenuated or suppressed by thermal pyrolysis.

The issues in numerical studies of supercritical aviation kerosene are clearly revealed, and the conjugation mechanism between thermal pyrolysis and convective heat transfer is initially presented.