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 study aims to numerically investigate heat transport in a trapezoidal cavity using hybrid nanoparticles (Ag-$Al_2O_3$). Unlike previous studies, this one covers magnetohydrodynamics, joule heating with viscous dissipation, heat absorption and generation. The left and right sides of the chasm are frigid. The upper wall heats, whereas the bottom wall remains adiabatic.

After reducing the system of dimensional equations to dimensionless equations, the authors use the Galerkin finite element method to solve them numerically. Geometric parameters affect heating efficiency; thus, the authors use flow metrics such as the Reynold number Re, magnetic parameter M, volume fraction coefficient, heat absorption and Eckert number Ec. The authors use the finite volume method to solve the governing equations after converting them to dimensionless form. The authors also try the artificial neural network method to predict the innovative cavity’s heat response in future scenarios. Transition state charts, regression analysis, MSE and error histograms accelerate, smooth and accurately converge solutions.

As the magnetic parameter and Eckert number increase, the enclosure emits more heat. As Reynold and volume fraction coefficients rise, the Nusselt number falls. It rose as magnetic, Eckert and heat absorption characteristics increased. The average Nusselt number rises with Reynolds and volume fraction coefficients. The magnetic, Eckert and heat absorption characteristics have inverse values.

This study numerically investigates heat transport in a trapezoidal cavity using hybrid nanoparticles (Ag-$Al_2O_3$). Unlike previous studies, this one covers MHD, joule heating with viscous dissipation, heat absorption and generation. The left and right sides of the chasm are frigid. The upper wall heats, whereas the bottom wall remains adiabatic.

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).

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.

To improve the thermal performance of base fluid, nanoparticles of three types are dispersed in the base fluid. A novel theory of non-Fourier heat transfer is used for design and development of models. The thermal performance of sample fluids is compared to determine which types of combination of nanoparticles are the best for an optimized enhancement in thermal performance of fluids. This article aims to: (i) investigate the impact of nanoparticles on thermal performance; and (ii) implement the Galerkin finite element method (GFEM) to thermal problems.

The mathematical models are developed using novel non-Fourier heat flux theory, conservation laws of computational fluid dynamics (CFD) and no-slip thermal boundary conditions. The models are approximated using thermal boundary layer approximations, and transformed models are solved numerically using GFEM. A grid-sensitivity test is performed. The accuracy, correction and stability of solutions is ensured. The numerical method adopted for the calculations is validated with published data. Quantities of engineering interest, i.e. wall shear stress, wall mass flow rate and wall heat flux, are calculated and examined versus emerging rheological parameters and thermal relaxation time.

The thermal relaxation time measures the ability of a fluid to restore its original thermal state, called thermal equilibrium and therefore, simulations have shown that the thermal relaxation time associated with a mono nanofluid has the most substantial effect on the temperature of fluid, whereas a ternary nanofluid has the smallest thermal relaxation time. A ternary nanofluid has a wider thermal boundary thickness in comparison with base and di- and mono nanofluids. The wall heat flux (in the case of the ternary nanofluids) has the most significant value compared with the wall shear stresses for the mono and hybrid nanofluids. The wall heat and mass fluxes have the highest values for the case of non-Fourier heat and mass diffusion compared to the case of Fourier heat and mass transfer.

An extensive literature review reveals that no study has considered thermal and concentration memory effects on transport mechanisms in fluids of cross-rheological liquid using novel theory of heat and mass [presented by Cattaneo (Cattaneo, 1958) and Christov (Christov, 2009)] so far. Moreover, the finite element method for coupled and nonlinear CFD problems has not been implemented so far. To the best of the authors’ knowledge for the first time, the dynamics of wall heat flow rate and mass flow rate under simultaneous effects of thermal and solute relaxation times, Ohmic dissipation and first-order chemical reactions are studied.

This study aims to undertake a comprehensive examination of heat transfer by convection in porous systems with top and bottom walls insulated and differently heated vertical walls under a magnetic field. For a specific nanofluid, the study aims to bring out the effects of different segmental heating arrangements.

An existing in-house code based on the finite volume method has provided the numerical solution of the coupled nondimensional transport equations. Following a validation study, different explorations include the variations of Darcy–Rayleigh number (Ra_m = 10–10⁴), Darcy number (Da = 10^–5–10^–1) segmented arrangements of heaters of identical total length, porosity index (ε = 0.1–1) and aspect ratio of the cavity (AR = 0.25–2) under Hartmann number (Ha = 10–70) and volume fraction of φ = 0.1% for the nanoparticles. In the analysis, there are major roles of the streamlines, isotherms and heatlines on the vertical mid-plane of the cavity and the profiles of the flow velocity and temperature on the central line of the section.

The finding of a monotonic rise in the heat transfer rate with an increase in Ra_m from 10 to 10⁴ has prompted a further comparison of the rate at Ra_m equal to 10⁴ with the total length of the heaters kept constant in all the cases. With respect to uniform heating of one entire wall, the study reveals a significant advantage of 246% rate enhancement from two equal heater segments placed centrally on opposite walls. This rate has emerged higher by 82% and 249%, respectively, with both the segments placed at the top and one at the bottom and one at the top. An increase in the number of centrally arranged heaters on each wall from one to five has yielded 286% rate enhancement. Changes in the ratio of the cavity height-to-length from 1.0 to 0.2 and 2 cause the rate to decrease by 50% and increase by 21%, respectively.

Further research with additional parameters, geometries and configurations will consolidate the understanding. Experimental validation can complement the numerical simulations presented in this study.

This research contributes to the field by integrating segmented heating, magnetic fields and hybrid nanofluid in a porous flow domain, addressing existing research gaps. The findings provide valuable insights for enhancing thermal performance, and controlling heat transfer locally, and have implications for medical treatments, thermal management systems and related fields. The research opens up new possibilities for precise thermal management and offers directions for future investigations.

This study aims to delve into the coupled mixed convective heat transport process within a grooved channel cavity using CuO-water nanofluid and an inclined magnetic field. The cavity undergoes isothermal heating from the bottom, with variations in the positions of heated walls across the grooved channel. The aim is to assess the impact of heater positions on thermal performance and identify the most effective configuration.

Numerical solutions to the evolved transport equations are obtained using a finite volume method-based indigenous solver. The dimensionless parameters of Reynolds number (1 ≤ Re ≤ 500), Richardson number (0.1 ≤ Ri ≤ 100), Hartmann number (0 ≤ Ha ≤ 70) and magnetic field inclination angle (0° ≤ γ ≤ 180°) are considered. The solved variables generate both local and global variables after discretization using the semi-implicit method for pressure linked equations algorithm on nonuniform grids.

The study reveals that optimal heat transfer occurs when the heater is positioned at the right corner of the grooved cavity. Heat transfer augmentation ranges from 0.5% to 168.53% for Re = 50 to 300 compared to the bottom-heated case. The magnetic field’s orientation significantly influences the average heat transfer, initially rising and then declining with increasing inclination angle. Overall, this analysis underscores the effectiveness of heater positions in achieving superior thermal performance in a grooved channel cavity.

This concept can be extended to explore enhanced thermal performance under various thermal boundary conditions, considering wall curvature effects, different geometry orientations and the presence of porous structures, either numerically or experimentally.

The findings are applicable across diverse fields, including biomedical systems, heat exchanging devices, electronic cooling systems, food processing, drying processes, crystallization, mixing processes and beyond.

This work provides a novel exploration of CuO-water nanofluid flow in mixed convection within a grooved channel cavity under the influence of an inclined magnetic field. The influence of different heater positions on thermomagnetic convection in such a cavity has not been extensively investigated before, contributing to the originality and value of this research.

This paper aims to investigate the tribological benefits of a biomimetic teardrop surface texture inspired by snakeskin compared to conventional surface textures with the help of geometrical and flow parameters using computational fluid dynamics techniques.

The lubricant is assumed to be Newtonian, and the flow is laminar with constant viscosity and isothermal property. The governing equations, continuity and Navier–Stokes equation, are discretised by the finite volume method, and cavitation modelling is included. The discretisation for the momentum equations is carried out using the second-order difference method for the SIMPLEC algorithm of pressure–velocity coupling.

The results indicate that biomimetic teardrop surface texturing performs better than conventional shapes surface textures in improving tribological performance. Furthermore, the parallel texture orientation along with the flow generates a high-pressure distribution relative to other orientations. Surface texture area density also highly influences the load-carrying capacity, which is optimum at 29%. Zigzag pattern arrangement performs better compared to linear pattern arrangement of texturing.

The paper proposes that this unique biomimetic teardrop shape can give better tribological performance than conventional shapes.

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

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.