Search results

The boundary integral method (BIM) is very attractive to practicing engineers as it reduces the dimensionality of the problem by one, thereby making the procedure computationally inexpensive compared to its peers. The principal feature of this technique is the limitation of all its computations to only the boundaries of the domain. Although the procedure is well developed for the Laplace equation, the Poisson equation offers some computational challenges. Nevertheless, the literature provides a couple of solution methods. This paper revisits an alternate approach that has not gained much traction within the community. The purpose of this paper is to address the main bottleneck of that approach in an effort to popularize it and critically evaluate the errors introduced into the solution by that method.

The primary intent in the paper is to work on the particular solution of the Poisson equation by representing the source term through a Fourier series. The evaluation of the Fourier coefficients requires a rectangular domain even though the original domain can be of any arbitrary shape. The boundary conditions for the homogeneous solution gets modified by the projection of the particular solution on the original boundaries. The paper also develops a new Gauss quadrature procedure to compute the integrals appearing in the Fourier coefficients in case they cannot be analytically evaluated.

The current endeavor has developed two different representations of the source terms. A comprehensive set of benchmark exercises has successfully demonstrated the effectiveness of both the methods, especially the second one. A subsequent detailed analysis has identified the errors emanating from an inadequate number of boundary nodes and Fourier modes, a high difference in sizes between the particular solution and the original domains and the used Gauss quadrature integration procedures. Adequate mitigation procedures were successful in suppressing each of the above errors and in improving the solution accuracy to any desired level. A comparative study with the finite difference method revealed that the BIM was as accurate as the FDM but was computationally more efficient for problems of real-life scale. A later exercise minutely analyzed the heat transfer physics for a fin after validating the simulation results with the analytical solution that was separately derived. The final set of simulations demonstrated the applicability of the method to complicated geometries.

First, the newly developed Gauss quadrature integration procedure can efficiently compute the integrals during evaluation of the Fourier coefficients; the current literature lacks such a tool, thereby deterring researchers to adopt this category of methods. Second, to the best of the author’s knowledge, such a comprehensive error analysis of the solution method within the BIM framework for the Poisson equation does not currently exist in the literature. This particular exercise should go a long way in increasing the confidence of the research community to venture into this category of methods for the solution of the Poisson equation.

This paper presents a Monotonic Unbounded Schemes Transformer (MUST) approach to bound/monotonize (remove undershoots and overshoots) unbounded spatial differencing schemes automatically, and naturally. Automatically means the approach (1) captures the critical cell Peclet number when an unbounded scheme starts to produce physically unrealistic solution automatically, and (2) removes the undershoots and overshoots as part of the formulation without requiring human interventions. Naturally implies, all the terms in the discretization equation of the unbounded spatial differencing scheme are retained.

The authors do not formulate new higher-order scheme. MUST transforms an unbounded higher-order scheme into a bounded higher-order scheme.

The solutions obtained with MUST are identical to those without MUST when the cell Peclet number is smaller than the critical cell Peclet number. For cell Peclet numbers larger than the critical cell Peclet numbers, MUST sets the nodal values to the limiter value which can be derived for the problem at-hand. The authors propose a way to derive the limiter value. The authors tested MUST on the central differencing scheme, the second-order upwind differencing scheme and the QUICK differencing scheme. In all cases tested, MUST is able to (1) capture the critical cell Peclet numbers; the exact locations when overshoots and undershoots occur, and (2) limit the nodal value to the value of the limiter values. These are achieved by retaining all the discretization terms of the respective differencing schemes naturally and accomplished automatically as part of the discretization process. The authors demonstrated MUST using one-dimensional problems. Results for a two-dimensional convection–diffusion problem are shown in Appendix to show generality of MUST.

The authors present an original approach to convert any unbounded scheme to bounded scheme while retaining all the terms in the original discretization equation.

The present numerical investigation examines the magnetohydrodynamic (MHD) double diffusion natural convection of power-law non-Newtonian nano-encapsulated phase change materials (NEPCMs) in a trapezoidal cavity.

The governing Navier-Stokes, energy and concentration equations based on the Cartesian curvilinear coordinates are solved using the collocated grid arrangement’s finite volume method. The in-house FORTRAN code is validated with the different benchmark problems. The NEPCM nanoparticles consist of a core-shell structure with Phase Change Material (PCM) at the core. The enclosure, shaped as a trapezoidal hollow, features a warmed (T_h) left wall and a cold (T_c) right wall. Various parameters are considered, including the power law index (0.6 ≤ n ≤ 1.4), Hartmann number (0 ≤ Ha ≤ 30), Rayleigh number (10⁴ ≤ Ra ≤ 10⁵) and fixed variables such as buoyancy ratio (Br = 0.8), Prandtl number (Pr = 6.2), Lewis number (Le = 5), fusion temperature (Θ_f = 0.5) and volume fraction (ϕ = 0.04).

The findings indicate a decrease in local Nusselt (Nu) and Sherwood (Sh) numbers with increasing Hartmann numbers (Ha). Additionally, for a shear-thinning fluid (n = 0.6) results in the maximum local Nu and Sh values. As the Rayleigh number (Ra) increases from 10⁴ to 10⁵, the structured vortex in the streamline pattern is disturbed. Furthermore, for different Ra values, an increase in n from 0.6 to 1.4 leads to a 67.43% to 76.88% decrease in average Nu and a 70% to 77% decrease in average Sh.

This research is for two-dimensioal laminar flow only.

PCMs represent a class of practical substances that behave as a function of temperature and have the innate ability to absorb, release and store heated energy in the form of hidden fusion enthalpy, or heat. They are valuable in these systems as they can store significant energy at a relatively constant temperature through their latent heat phase change.

As per the literature review and the authors’ understanding, an examination has never been conducted on MHD double diffusion natural convection of power-law non-Newtonian NEPCMs within a trapezoidal enclosure. The current work is innovative since it combines NEPCMs with the effect of magnetic field Double diffusion Natural Convection of power-law non-Newtonian NEPCMs in a Trapezoidal enclosure. This outcome can be used to improve thermal management in energy storage systems, increasing safety and effectiveness.

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.

The performance of a bladeless Tesla turbine is closely tied to momentum diffusion, kinetic energy transfer and wall shear stress generation on its rotating disks. The surface roughness adds complexity of flow analysis in such a domain. This paper aims to assess the effect of roughness on flow structures and the application of roughness models in flow cross sections with submillimeter height, including both stationary and rotating walls.

This research starts with the examination of flow over a rough flat plate, and then proceeds to study flow within minichannels, evaluating the effect of roughness on flow characteristics. An in-house test stand validates the numerical solutions of minichannel. Finally, flow through the minichannel with corotating walls was analyzed. The k-ω SST turbulent model and Aupoix's roughness method are used for numerical simulations.

The findings emphasize the necessity of considering the constricted dimensions of the flow cross section, thereby improving the alignment of derived results with theoretical estimations. Moreover, this study explores the effects of roughness on flow characteristics within the minichannel with stationary and rotating walls, offering valuable insights into this intricate phenomenon, and depicts the appropriate performance of chosen roughness model in studied cases.

The originality of this investigation is the assessment and validation of flow characteristics inside minichannel with stationary and corotating walls when the roughness is implemented. This phenomenon, along with the effect of roughness on the transportation of kinetic energy to the rough surface of a minichannel in an in-house test setup, is assessed.

The purpose of this study is to investigate the sealing performance of different deep groove mechanical seals by considering the changing law of dynamic pressure effect and temperature gradient caused by high speed and high pressure.

A thermohydrodynamic lubrication model (THD) of the mechanical seal was constructed and solved using the commercial software FLUENT. The pressure and temperature distributions of the fluid under different groove types, as well as the sealing performance under different pressures, rotational speeds and sealing gaps, are obtained.

The annular groove (AG) can effectively reduce the temperature, and the T-type spiral groove (STG) can effectively inhibit the leakage. The increase of pressure and rotational speed leads to the enhancement of dynamic pressure effect and the increase of leakage, while the sealing gap increases and the leakage increases while taking away more heat. The choice of groove type is very important to the impact of sealing performance.

In consideration of the beneficial effect of deep grooves on cooling performance, the viscous temperature equation and the impact of the thermodynamic lubrication model are evaluated in conjunction with the sealing performance of four distinct groove types. This approach provides a theoretical basis for the optimal design of mechanical seals.

The peer review history for this article is available at: https://publons.com/publon/10.1108/ILT-05-2024-0184/

This study aims to investigate the erosion wear rate of a stainless steel automobile exhaust manifold, both computationally and physically.

The experiment was performed on a motorcycle exhaust manifold as well as on a 3D model, created using SolidWorks 2022 CAD software. The analysis was later achieved using ANSYS 19.2 simulation software using Fluent – code.

The analysis of solid particle erosion in the exhaust manifold revealed that erosion wear is concentrated predominantly at the extrados of the manifold, with the most significant wear occurring at the lowermost bend. The erosion wear rate increases with larger particulate sizes and varies among bends, with negligible wear observed in straight pipes. The SEM analysis further confirmed surface degradation, with rugged textures, pits and grooves indicating abrasive wear. Spine-like structures and fractured soot particles suggest erosive and abrasive forces caused by high-speed contact of exhaust gas compounds. Energy dispersive X-ray spectroscopy revealed significant carbon abundance, indicating carbonaceous compounds from fuel combustion, along with notable amounts of oxygen and iron, typical of oxidized metallic constituents. The discrete phase modeling (DPM) analysis highlighted peak particulate matter deposition at the first bend exit, with maximum concentrations observed at specific angles. This deposition is influenced by centrifugal force, leading to increased PM concentration at outer bend walls. Velocity magnitude contours showed asymmetrical flow profiles, with high turbulence levels and secondary flow induced by centrifugal effects in bend areas. Dynamic pressure contours revealed varying pressures at intrados and extrados, with maximum pressure observed at the intrados of the manifold’s bends. These findings provide valuable insights into erosion wear, particulate dispersion and flow dynamics within the exhaust manifold.

The study investigated an automobile exhaust manifold model using ANSYS Fluent code and DPM to analyze erosion wear rate phenomena and its various constituents. This analysis was conducted in comparison with a physically eroded sample. The study offers insights into the mechanism underlying the exhaust manifold of an automobile.

This study aims to focus on understanding how different jet angles and Reynolds numbers influence the phase change materials’ (PCMs) melting process and their capacity to store energy. This approach is intended to offer novel insights into enhancing thermal energy storage systems, particularly for applications where heat transfer efficiency and energy storage are critical.

The research involved an experimental and numerical analysis of PCM with a melting temperature range of 22 °C–26°C under various conditions. Three different jet angles (45°, 90° and 135°) and two container angles (45° and 90°) were tested. Additionally, two different Reynolds numbers (2,235 and 4,470) were used to explore the effects of jet outlet velocities on PCM melting behaviour. The study used a circular container and analysed the melting process using the hot air inclined jet impingement (HAIJI) method.

The obtained results showed that the average temperature for the last time step at Ф = 90° and Re = 4,470 is 6.26% higher for Ф = 135° and 14.23% higher for Ф = 90° compared with the 45° jet angle. It is also observed that the jet angle, especially for Ф = 90°, is a much more important factor in energy storage than the Reynolds number. In other words, the jet angle can be used as a passive control parameter for energy storage.

This study offers a novel perspective on the effective storage of waste heat transferred with air, such as exhaust gases. It provides valuable insights into the role of jet inclination angles and Reynolds numbers in optimizing the melting and energy storage performance of PCMs, which can be crucial for enhancing the efficiency of thermal energy storage systems.

The lubricating fluid stored in the porous matrix will spontaneously exude to supplement the lubricating film in the damaged area, thus ensuring the long-term self-lubricating function of the porous surface. To reveal the repair mechanism of oil film, it is necessary to understand the flow characteristics of oil in micropores. The purpose of this study guides the design of micropore structure to realize the rapid exudation of oil to the porous surface and the rapid repair of the lubricating film.

In this paper, cylindrical orifice, convergent orifice and divergent orifice were studied. The numerical model of lubricating oil exudation in micropores was established. The distribution characteristics of oil pressure, velocity and three-phase contact line in the process of oil exudation were investigated. The effects of different orifice shapes and orifice structure parameters on the pinning and spreading characteristics of oil droplet were analyzed. Then the internal mechanisms of oil droplet formation and spread on the orifice surface were summarized.

The results show that during the process of oil exudation, the three-phase contact line of the oil drop is pinned once at the edge of the cylindrical and convergent orifice. Compared with the three orifice structures, the inlet pressure of the oil drop is low, and the oil velocity at the pinning point is stable in the divergent orifice. Resulting in favorable oil exudation. It is easier for oil droplet to depin by appropriately reducing the wall wetting angle, increasing the aperture or controlling the wall inclination angle. Ensure the self-healing and long-lasting lubrication film of porous oil-bearing surfaces.

The effect of pore structure on the flow behavior of lubricating fluid has always been concerned. But the mechanism by which different orifice shape affect the pinning behavior of oil droplets is not yet clear, which is crucial for understanding the self-healing mechanism of oil films on porous surfaces. It is meaningful to analyze the mechanism of oil exudation and spreading on the porous surface of oil in the special orifice, to optimize the design of the orifice structure.

Orifice shape has influence on internal flow field parameters. There is no report on the influence of orifice shape on the film formation process of oil seepage and diffusion from pores. The effects of different orifice shapes and orifice structure parameters on the characteristics of oil droplet pinning and diffusion were studied.

The peer review history for this article is available at: https://publons.com/publon/10.1108/ILT-04-2024-0118/

This paper aims to investigate the effect of different wing height layouts on the aerodynamic performance and flow structure of high-speed train, in a train-wing coupling method with multiple tandem wings installed on the train roof.

The improved delayed detached eddy simulation method based on shear stress transport k- ω turbulence model has been used to conduct computational fluid dynamics simulation on the train with three different wing height layouts, at a Reynolds number of 2.8 × 10⁶. The accuracy of the numerical method has been validated by wind tunnel experiments.

The wing height layout has a significant effect on the lift, while its influence on the drag is weak. There are three distinctive vortex structures in the flow field: wingtip vortex, train body vortex and pillar vortex, which are influenced by the variation in wing height layout. The incremental wing layout reduces the mixing and merging between vortexes in the flow field, weakening the vorticity and turbulence intensity. This enhances the pressure difference between the upper and lower surfaces of both the train and wings, thereby increasing the overall lift. Simultaneously, it reduces the slipstream velocity at platform and trackside heights.

This paper contributes to understanding the aerodynamic characteristics and flow structure of a high-speed train coupled with wings. It provides a reference for the design aiming to achieve equivalent weight reduction through aerodynamic lift synergy in trains.