Search results

The purpose of this paper is to investigate the two-phase flow problems involving gas–liquid mixture.

The governed equations consist of a range of conservation laws modeling a classification of two-phase flow phenomena subjected to a velocity nonequilibrium for the gas–liquid mixture. Effects of the relative velocity are accounted for in the present model by a kinetic constitutive relation coupled to a collection of specific equations governing mass and volume fractions for the gas phase. Unlike many two-phase models, the considered system is fully hyperbolic and fully conservative. The suggested relaxation approach switches a nonlinear hyperbolic system into a semilinear model that includes a source relaxation term and characteristic linear properties. Notably, this model can be solved numerically without the use of Riemann solvers or linear iterations. For accurate time integration, a high-resolution spatial reconstruction and a Runge–Kutta scheme with decreasing total variation are used to discretize the relaxation system.

The method is used in addressing various nonequilibrium two-phase flow problems, accompanied by a comparative study of different reconstructions. The numerical results demonstrate the suggested relaxation method’s high-resolution capabilities, affirming its proficiency in delivering accurate simulations for flow regimes characterized by strong shocks.

While relaxation methods exhibit notable performance and competitive features, as far as we are aware, there has been no endeavor to address nonequilibrium two-phase flow problems using these methods.

This paper aims to derive a reduced-order model for the heat transfer across the interface between a millimetric thermocapillary liquid bridge from silicone oil and the surrounding ambient gas.

Numerical solutions for the two-fluid model are computed covering a wide parametric space, making a total of 2,800 numerical flow simulations. Based on the computed data, a reduced single-fluid model for the liquid phase is devised, in which the heat transfer between the liquid and the gas is modeled by Newton’s heat transfer law, albeit with a space-dependent Biot function Bi(z), instead of a constant Biot number Bi.

An explicit robust fit of Bi(z) is obtained covering the whole range of parameters considered. The single-fluid model together with the Biot function derived yields very accurate results at much lesser computational cost than the corresponding two-phase fully-coupled simulation required for the two-fluid model.

Using this novel Biot function approach instead of a constant Biot number, the critical Reynolds number can be predicted much more accurately within single-phase linear stability solvers.

The Biot function for thermocapillary liquid bridges is derived from the full multiphase problem by a robust multi-stage fit procedure. The derived Biot function reproduces very well the theoretical boundary layer scalings.

The purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.

In the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.

The larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.

The concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

This paper aims to establish a lattice Boltzmann (LB) method for solid-liquid phase transition (SLPT) from the pore scale to the representative elementary volume (REV) scale. By applying this method, detailed information about heat transfer and phase change processes within the pores can be obtained, while also enabling the calculation of larger-scale SLPT problems, such as shell-and-tube phase change heat storage systems.

Three-dimensional (3D) pore-scale enthalpy-based LB model is developed. The computational input parameters at the REV scale are derived from calculations at the pore scale, ensuring consistency between the two scales. The approaches to reconstruct the 3D porous structure and determine the REV of metal foam were discussed. The implementation of conjugate heat transfer between the solid matrix and the solid−liquid phase change material (SLPCM) for the proposed model is developed. A simple REV-scale LB model under the local thermal nonequilibrium condition is presented. The method of bridging the gap between the pore-scale and REV-scale enthalpy-based LB models by the REV is given.

This coupled method facilitates detailed simulations of flow, heat transfer and phase change within pores. The approach holds promise for multiscale calculations in latent heat storage devices with porous structures. The SLPT of the heat sinks for electronic device thermal control was simulated as a case, demonstrating the efficiency of the present models in designing and optimizing SLPT devices.

A coupled pore-scale and REV-scale LB method as a numerical tool for investigating phase change in porous materials was developed. This innovative approach allows for the capture of details within pores while addressing computations over a large domain. The LB method for simulating SLPT from the pore scale to the REV scale was given. The proposed method addresses the conjugate heat transfer between the SLPCM and the solid matrix in the enthalpy-based LB model.

Water jet propulsion is widely used in various military and civilian fields due to its advantages of simple structure and high propulsion efficiency. The process of mooring involves utilizing specially designed equipment to secure a ship at a designated berth. During the process of water jet propulsion, the single propeller operates within a complex and turbulent three-dimensional flow. Hence, studying the coupling between the water jet propeller and the hull is critical to comprehending the characteristics of the device and the distribution of the flow field in detail.

Firstly, we conducted computational fluid dynamics (CFD)-based self-propulsion calculations to evaluate the interaction between the hull and the propeller. We subsequently analyzed the propeller's performance and the forces acting on the hull to understand how the presence or absence of the hull influenced the water jet propeller. Finally, we performed calculations and analysis of the cavitation characteristics of the coupling between the hull and the water jet propeller, considering different rotational speeds and water depths at the bottom of the pool.

The study demonstrated that the presence of the hull boundary layer under the hull-propeller coupling condition led to reduced uniformity of propeller inlet flow and lower efficiency of the propulsion pump. However, it also increased the bias toward low-flow conditions. Additionally, increasing the impeller speed led to a gradual increase in the cavitation volume within the water jet propeller, resulting in a gradual decrease in the propeller's performance.

This research provides the technical support required for effective design and operation of water jet propulsion systems. This paper involves studying and analyzing the performance and flow field of the coupling between the hull and propeller under mooring conditions with a specified hull model.

This paper aims to study the effect of concentric hot circular cylinder inside egg-cavity porous-copper nanofluid on natural convection phenomena.

The finite element method–based Galerkin approach is applied to solve numerically the set of governing equations with appropriate boundary conditions.

The effects of different range parameters, such as Darcy number (10–3 = Da = 10–1), Rayleigh number (103 = Ra = 106), nanoparticle volume fraction (0 = ϑ = 0.06) and eccentricity (−0.3 = e = 0.1) on the fluid flow represent by stream function and heat transfer represent by temperature distribution, local and average Nusselt numbers.

A comparison between oval shape and concentric circular concentric cylinder was investigated.

In the current numerical study, heat transfer by natural convection was identified inside the new design of egg-shaped cavity as a result of the presence of a circular inside it supported by a porous medium filled with a nanofluid. After reviewing previous studies and considering the importance of heat transfer by free convection inside tubes for many applications, to the best of the authors’ knowledge, the current work is the first study that deals with a study and comparison between the common shape (concentric circular tubes) and the new shape (egg-shaped cavity).

The purpose of this paper is to introduce a double-stator switched reluctance machine (DS-SRM) for electric vehicles (EVs) and to propose multi-mode operations for this machine.

Analysis of flux linkage distributions and torque characteristics using finite element method (FEM). Building a dynamic simulation model based on electromagnetic characteristics, mathematical equations and mechanical motion equations of the DS-SRM drive system. The paper proposes multi-mode operations (inner-stator excitation mode, outer-stator excitation mode and double-stator excitation mode) based on motor working regions. It also conducts simulation and experimental results to verify the effectiveness of the proposed multi-mode operations strategies and control schemes.

There is almost no electromagnetic coupling between the inner and outer stators due to the specially designed rotor structure and optimized windings polarity configuration. Analysis of flux linkage distributions and torque characteristics verified the independence of inner and outer stators. Proposal of multi-mode operations and corresponding control rules achieved the smooth switching between different modes.

The paper introduced the DS-SRM for EVs and proposed multi-mode operations, along with control rules, to optimize its performance. The specially designed rotor structure, optimized winding polarity configuration, and the proposed multi-mode operations contribute to the originality of the research.

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.

This study aims to perform three-dimensional numerical computations for blood flow through a double stenosed carotid artery. Pulsatile flow with Womersley number (Wo) of 4.65 and Reynolds number (Re) of 425, based on the diameter of normal artery and average velocity of inlet pulse, was considered.

Finite volume method based ANSYS Fluent 20.1 was used for solving the governing equations of three-dimensional, laminar, incompressible and non-Newtonian blood flow. A high-quality grid with sufficient refinement was generated using ICEM CFD 20.1. The time-averaged flow field was captured to investigate the effect of severity and eccentricity on the lumen flow characteristics.

The results show that an increase in interspacing between blockages brings shear layer instability within the region between two blockages. The velocity profile and wall shear stress distribution are found to be majorly influenced by eccentricity. On the other hand, their peak magnitude is found to be primarily influenced by severity. Results have also demonstrated that the presence of eccentricity in stenosis would assist in flow development.

Variation in severity and interspacing was considered with a provision of eccentricity equal to 10% of diameter. Eccentricity refers to the offset between the centreline of stenosis and the centreline of normal artery. For the two blockages, severity values of 40% and 60% based on diameter reduction were permuted, giving rise to four combinations. For each combination, three values of interspacing in the multiples of normal artery diameter (D), viz. 4D, 6D and 8D were considered.

This study aims to minimize the pressure drop across wavy microchannels using secondary branches without compromising its capacity to transfer the heat. The impact of secondary flows on the pressure drop and heat transfer capabilities at different Reynolds numbers are investigated numerically for different wavy microchannels. Finally, different channels are evaluated using performance evaluation criteria to determine their effectiveness.

To investigate the flow and heat transfer capabilities in wavy microchannels having secondary branches, a 3D conjugate heat transfer model based on finite volume method is used. In conventional wavy microchannel, secondary branches are introduced at crest and trough locations. For the numerical simulation, a single symmetrical channel is used to minimize computational time and resources and the flow within the channels remains single-phase and laminar.

The findings indicate that the suggested secondary channels notably improve heat transfer and decrease pressure drop within the channels. At lower flow rates, the secondary channels demonstrate superior performance in terms of heat transfer. However, the performance declines as the flow rate increased. With the same amplitude and wavelength, the introduction of secondary channels reduces the pressure drop compared with conventional wavy channels. Due to the presence of secondary channels, the flow splits from the main channel, and part of the core flow gets diverted into the secondary channel as the flow takes the path of minimum resistance. Due to this flow split, the core velocity is reduced. An increase in flow area helps in reducing pressure drop.

Many complex and intricate microchannels are proposed by the researchers to augment heat dissipation. There are challenges in the fabrication of microchannels, such as surface finish and achieving the required dimensions. However, due to the recent developments in metal additive manufacturing and microfabrication techniques, the complex shapes proposed in this paper are feasible to fabricate.

Wavy channels are widely used in heat transfer and micro-fluidics applications. The proposed wavy microchannels with secondary channels are different when compared to conventional wavy channels and can be used practically to solve thermal challenges. They help achieve a lower pressure drop in wavy microchannels without compromising heat transfer performance.