International Journal of Numerical Methods for Heat & Fluid Flow最新文献

Purpose

The purpose of this paper is to study the two-dimensional micropolar fluid flow with conjugate heat transfer and mass transpiration. The considered nanofluid has graphene nanoparticles.

Design/methodology/approach

Governing nonlinear partial differential equations are converted to nonlinear ordinary differential equations by similarity transformation. Then, to analyze the flow, the authors derive the dual solutions to the flow problem. Biot number and radiation effect are included in the energy equation. The momentum equation was solved by using boundary conditions, and the temperature equation solved by using hypergeometric series solutions. Nusselt numbers and skin friction coefficients are calculated as functions of the Reynolds number. Further, the problem is governed by other parameters, namely, the magnetic parameter, radiation parameter, Prandtl number and mass transpiration. Graphene nanofluids have shown promising thermal conductivity enhancements due to the high thermal conductivity of graphene and have a wide range of applications affecting the thermal boundary layer and serve as coolants and thermal management systems in electronics or as heat transfer fluids in various industrial processes.

Findings

Results show that increasing the magnetic field decreases the momentum and increases thermal radiation. The heat source/sink parameter increases the thermal boundary layer. Increasing the volume fraction decreases the velocity profile and increases the temperature. Increasing the Eringen parameter increases the momentum of the fluid flow. Applications are found in the extrusion of polymer sheets, films and sheets, the manufacturing of plastic wires, the fabrication of fibers and the growth of crystals, among others. Heat sources/sinks are commonly used in electronic devices to transfer the heat generated by high-power semiconductor devices such as power transistors and optoelectronics such as lasers and light-emitting diodes to a fluid medium, thermal radiation on the fluid flow used in spectroscopy to study the properties of materials and also used in thermal imaging to capture and display the infrared radiation emitted by objects.

Originality/value

Micropolar fluid flow across stretching/shrinking surfaces is examined. Biot number and radiation effects are included in the energy equation. An increase in the volume fraction decreases the momentum boundary layer thickness. Nusselt numbers and skin friction coefficients are presented versus Reynolds numbers. A dual solution is obtained for a shrinking surface.

Purpose

The purpose of this study is two-fold. First, it aims to differentiate the response of a stretching jet encountering a quadratic air resistance from the classical jet shape formed in a frictionless medium. Second, it investigates how the resulting jet forms with and without air resistance, seeking evidence that supports the similarity flows frequently studied for stretching/moving thin bodies under the boundary layer approximation.

Design/methodology/approach

This study extends the established electrohydrodynamic stretching jet theory, used to model electrospinning or jet printing in the absence of air resistance, to encompass the impact of the retarding force on the jet stretching in both the cone and final regimes before it impinges on a substrate.

Findings

A close examination of the nonlinear governing equations reveals that the jet rapidly thins near the nozzle because of the combined action of viscous and electrical forces. In this region, the exponentially decaying jet receives further support from the air resistance, resulting in a closer alignment with the observed experimental jet. This exponential decay, accelerated by the inversely quadratic speed of the liquid particles, serves as clear evidence for the existence of a similarity flow over an exponentially stretching sheet. Furthermore, in the final regime, the jet stretching exhibits an algebraic decay in the absence of air friction, while with air resistance, it decays exponentially to reach a limiting speed. In the former case, a square root dependence of the stretching jet speed leads to the emergence of a similarity flow over a thin stretching jet, while in the latter case, a Sakiadis’ similarity flow appears over a continuously moving flat surface.

Practical implications

The analysis goes beyond jet hydrodynamics, delving into the interplay of electrostatic forces (including Coulomb’s law) and quadratic air drag, drawing upon experimental data on glycerol liquid presented in earlier publications.

Originality/value

Finally, the asymptotic behavior of the stretching jet under the combined influence of electrostatic pull and its electric currents because of bulk conduction and surface convection is validated through a comprehensive numerical simulation of the nonlinear system.

Purpose

Aerodynamic shape optimisation is a complex problem usually governed by transonic non-linear aerodynamics, a high dimensional design space and high computational cost. Consequently, the use of a numerical simulation approach can become prohibitive for some applications. This paper aims to propose a computationally efficient multi-fidelity method for the optimisation of two-dimensional axisymmetric aero-engine nacelles.

Design/methodology/approach

The nacelle optimisation approach combines a gradient-free algorithm with a multi-fidelity surrogate model. Machine learning based on artificial neural networks (ANN) is used as the modelling technique because of its ability to handle non-linear behaviour. The multi-fidelity method combines Reynolds-averaged Navier Stokes and Euler CFD calculations as high- and low-fidelity, respectively.

Findings

Ratios of low- and high-fidelity training samples to degrees of freedom of n_LF/n_DOFs = 50 and n_HF/n_DOFs = 12.5 provided a surrogate model with a root mean squared error less than 5% and a similar convergence to the optimal design space when compared with the equivalent CFD-in-the-loop optimisation. Similar nacelle geometries and aerodynamic flow topologies were obtained for down-selected designs with a reduction of 92% in the computational cost. This highlights the potential benefits of this multi-fidelity approach for aerodynamic optimisation within a preliminary design stage.

Originality/value

The application of a multi-fidelity technique based on ANN to the aerodynamic shape optimisation problem of isolated nacelles is the key novelty of this work. The multi-fidelity aspect of the method advances current practices based on single-fidelity surrogate models and offers further reductions in computational cost to meet industrial design timescales. Additionally, guidelines in terms of low- and high-fidelity sample sizes relative to the number of design variables have been established.

Purpose

The purpose of this study is to develop a new numerical algorithm to simulate the phase-field model.

Design/methodology/approach

First, the derivative of the temporal direction is discretized by a second-order linearized finite difference scheme where it conserves the energy stability of the mathematical model. Then, the isogeometric collocation (IGC) method is used to approximate the derivative of spacial direction. The IGC procedure can be applied on irregular physical domains. The IGC method is constructed based upon the nonuniform rational B-splines (NURBS). Each curve and surface can be approximated by the NURBS. Also, a map will be defined to project the physical domain to a simple computational domain. In this procedure, the partial derivatives will be transformed to the new domain by the Jacobian and Hessian matrices. According to the mentioned procedure, the first- and second-order differential matrices are built. Furthermore, the pseudo-spectral algorithm is used to derive the first- and second-order nodal differential matrices. In the end, the Greville Abscissae points are used to the collocation method.

Findings

In the numerical experiments, the efficiency and accuracy of the proposed method are assessed through two examples, demonstrating its performance on both rectangular and nonrectangular domains.

Originality/value

This research work introduces the IGC method as a simulation technique for the phase-field crystal model.

Purpose

Despite the reputation of the metal-based porous media for their ability to augment heat transfer as widely witnessed in the literature and practically operating heat exchanging applications, the coexisting penalty of the increased pressure drop demanding increased pumping power poses a major concern that invites the need for an alternate solution to handle this unsought outcome. Therefore, this study aims at providing a better solution to the existing cost and benefit scenarios to benefit a plethora of engineering applications including energy transfer, energy storage and energy conversion.

Design/methodology/approach

This work highlights on the property of stacked woven wire mesh porous media such as their stacking types, porous conditions and thickness scenarios that can potentially result in distinct trade-off scenarios. A vertical channel is numerical modelled by using REV scaled modelling technique using Darcy-Forchheimer and local thermal non-equilibrium models to illustrate the possibilities of this variety of trade off scenarios between the desirable heat transfer and the unsought flow resistance.

Findings

This work illustrates the advantages of wire mesh-based porous medium and its distinct potential in controlling the existing trade-offs between the cost and benefit aspects. It is found that by varying the features of wire mesh porous media, the interplay between the conflictingly existing characteristics can be much easily handled specific to distinct requirements associated with variety of engineering applications.

Originality/value

The study emphasizes on a new solution or methodology to handle the penalty of pressure drop associated with metal-based porous media. Through this study, a novel approach to control the ultimately costing pumping power at the benefit of increased heat transfer is provided considering various requirements that could be associated with any thermal management systems. Various possibilities and potentials of wire mesh porous media are illustrated highlighting on their benefit of ease with which the mentioned goals can be achieved.

Purpose

Sustainable urban rail transit requires noise barriers. However, these barriers’ durability varies due to the differing aerodynamic impacts they experience. The purpose of this paper is to investigate the aerodynamic discrepancies of trains when they meet within two types of rectangular noise barriers: fully enclosed (FERNB) and semi-enclosed with vertical plates (SERNBVB). The research also considers the sensitivity of the scale ratio in these scenarios.

Design/methodology/approach

A 1:16 scaled moving model test analyzed spatiotemporal patterns and discrepancies in aerodynamic pressures during train meetings. Three-dimensional computational fluid dynamics models, with scale ratios of 1:1, 1:8 and 1:16, used the improved delayed detached eddy simulation turbulence model and slip grid technique. Comparing scale ratios on aerodynamic pressure discrepancies between the two types of noise barriers and revealing the flow field mechanism were done. The goal is to establish the relationship between aerodynamic pressure at scale and in full scale.

Findings

The aerodynamic pressure on SERNBVB is influenced by the train’s head and tail waves, whereas for FERNB, it is affected by pressure wave and head-tail waves. Notably, SERNBVB's aerodynamic pressure is more sensitive to changes in scale ratio. As the scale ratio decreases, the aerodynamic pressure on the noise barrier gradually increases.

Originality/value

A train-meeting moving model test is conducted within the noise barrier. Comparison of aerodynamic discrepancies during train meets between two types of rectangular noise barriers and the relationship between the scale and the full scale are established considering the modeling scale ratio.

Purpose

This study aims to examine the plane wave reflection problem in micropolar orthotropic magneto-thermoelastic half space, considering the influence of impedance as a boundary in a nonlocal elasticity.

Design/methodology/approach

This study presents the novel formulation of governing partial differential equations for micropolar orthotropic medium with impact of nonlocal thermo-elasticity under magnetic field.

Findings

This study provides the numerical results validation for a particular numerical data and expression for the amplitude ratios of reflected waves and identifies the existence of four different waves, namely, quasi longitudinal displacement qCLD-wave, quasi thermal wave qCT-wave, quasi transverse displacement qCTD-wave and quasi-transverse micro-rotational qCTM-wave. The study derives the velocity equation giving the speed and phase velocity of these waves. The study also shows that the small-scale size effect gives significant impact on phase velocity.

Research limitations/implications

The graphical analysis examines the variation of speeds and coefficients of attenuation of these waves due to frequency, magnetic field and nonlocal parameters. Also, significant conclusions on the variation of reflection coefficient against nonlocal parameter, frequency, impedance parameter and angle of incidence are provided graphically.

Practical implications

The creation of more effective micropolar orthotropic anisotropic materials which are very useful in the daily life and their applications in earth science are greatly impacted by the findings of this study.

Originality/value

The authors of the submitted document initiated and produced it collectively, with equal contributions from all members.

Purpose

This study aims to investigate numerically the impact of the three-dimensional convective nanoliquid flow on a rotating frame embedded in the non-Darcy porous medium in the presence of activation energy. The cross-diffusion effects, i.e. Soret and Dufour effects, and heat generation are included in the study. The convective heating condition is applied on the bounding surface.

Design/methodology/approach

The control model consisted of a system of partial differential equations (PDE) with boundary constraints. Using suitable similarity transformation, the PDE transformed into an ordinary differential equation and solved numerically by the Runge–Kutta–Fehlberg method. The obtained results of velocity, temperature and solute concentration characteristics plotted to show the impact of the pertinent parameters. The heat and mass transfer rate and skin friction are also calculated.

Findings

It is found that both Biot numbers enhance the heat and mass distribution inside the boundary layer region. The temperature increases by increasing the Dufour number, while concentration decreases by increasing the Dufour number. The heat transfer is increased up to 8.1% in the presence of activation energy parameter (E). But, mass transfer rate declines up to 16.6% in the presence of E.

Practical implications

The applications of combined Dufour and Soret effects are in separation of isotopes in mixture of gases, oil reservoirs and binary alloys solidification. The nanofluid with porous medium can be used in chemical engineering, heat exchangers and nuclear reactor.

Social implications

This study is mainly useful for thermal sciences and chemical engineering.

Originality/value

The uniqueness in this research is the study of the impact of activation energy and cross-diffusion on rotating nanoliquid flow with heat generation and convective heating condition. The obtained results are unique and valuable, and it can be used in various fields of science and technology.

Purpose

This study aims to investigate magnetohydrodynamic (MHD) conjugate pure mixed convection considering interior heat production and resistive heating inside a square closed/open cavity featuring a rotating cylinder for aiding (clockwise) and opposing (counterclockwise) flow configurations. Moreover, the impacts of altering cylinder size and conductivity on the system’s overall performance to determine optimum conditions are examined in this investigation.

Design/methodology/approach

The closed chamber is differentially heated by keeping high and low temperatures at the vertical boundaries. In contrast, the open cavity has a heated left wall and an open right boundary. The Galerkin finite element method is used to solve the Navier–Stokes and the thermal energy equations, which construct the present study’s mathematical framework. Numerical simulations are conducted for the specified ranges of several controlling parameters: Reynolds (31.62 ≤ Re ≤ 1000), Grashof (10³ ≤ Gr ≤ 10⁶) and Hartmann numbers (0 ≤ Ha ≤ 31.62), and volumetric heat generation coefficient (Δ = 0, 3).

Findings

When Gr, Re and Ha simultaneously increase, the average Nusselt number along the warmed boundary rises accordingly. Conversely, interior heat production lowers heat transmission within the computational domain, which is also monitored regarding mean fluid temperature, overall entropy production and thermal performance criterion. Finally, the open cavity confirms better thermal performance than the closed cavity.

Originality/value

Comprehending the impacts of the magnetic field, Joule heating, internal heat generation and enclosed or open boundary on pure MHD combined free-forced convective flow offers valuable understandings of temperature fluctuations, velocity propagations, heat transport and irretrievable energy loss in numerous engineering applications.

Purpose

Cilia serves numerous biological functions in the human body. Malfunctioning of nonmotile or motile cilia will have different kinds of consequences for human health. More specifically, the directed and rhythmic beat of motile cilia facilitates the unidirectional flow of fluids that are crucial in both homeostasis and the development of ciliated tissues. In cilia-dependent hydrodynamic flows, tapering geometries look a lot like the structure of biological pathways and vessels, like airways and lymphatic vessels. In this paper, the Carreau fluid model through the cilia-assisted tapered channel (asymmetric) under the influence of induced magnetic field and convective heat transfer is investigated.

Design/methodology/approach

Lubrication theory is a key player in the mathematical formulation of momentum, magnetic field and energy equations. The formulated nonlinear and coupled differential equations are solved with the aid of the homotopy perturbation method (HPM). The graphical results are illustrated with the help of the computational software “Mathematica.”

Findings

The impact of diverse emerging physical parameters on velocity, induced magnetic field, pressure rise, current density and temperature profiles is presented graphically. It is observed that the cilia length parameter supported the velocity and current density profiles, while the Hartman number and Weissenberg number were opposed. A promising effect of emerging parameters on streamlines is also perceived.

Originality/value

The study provides novel aspects of cilia-driven induced magnetohydrodynamics flow of Carreau fluid under the influence of induced magnetic field and convective heat transfer through the asymmetric tapered channel.