Search results

The main objective of this paper is to investigate the response of human skin to an intense temperature drop at the surface. In addition, this paper aims to evaluate the efficiency of finite difference and finite volume methods in solving the highly nonlinear form of Pennes’ bioheat equation.

One-dimensional linear and nonlinear forms of Pennes’ bioheat equation with uniform grids were used to study the behavior of human skin. The specific heat capacity, thermal conductivity and blood perfusion rate were assumed to be linear functions of temperature. The nonlinear form of the bioheat equation was solved using the Newton linearization method for the finite difference method and the Picard linearization method for the finite volume method. The algorithms were validated by comparing the results from both methods.

The study demonstrated the capacity of both finite difference and finite volume methods to solve the one-dimensional and highly nonlinear form of the bioheat equation. The investigation of human skin’s thermal behavior indicated that thermal conductivity and blood perfusion rate are the most effective properties in mitigating a surface temperature drop, while specific heat capacity has a lesser impact and can be considered constant.

This paper modeled the transient heat distribution within human skin in a one-dimensional manner, using temperate-dependent physical properties. The nonlinear equation was solved with two numerical methods to ensure the validity of the results, despite the complexity of the formulation. The findings of this study can help in understanding the behavior of human skin under extreme temperature conditions, which can be beneficial in various fields, including medical and engineering.

A novel method for modelling permanent magnets is investigated based on numerical approximations with rational functions. This study aims to introduce the AAA algorithm and other recently developed, cutting-edge mathematical tools, which provide outstandingly fast and accurate numerical computation of potentials and vector fields.

First, the AAA algorithm is briefly introduced along with its main variants and other advanced mathematical tools involved in the modelling. Then, the analysis of a circular Halbach array with a one-pole pair is carried out by means of the AAA-least squares method, focusing on vector potential and flux density in the bore and validating results by means of classic finite element software. Finally, the investigation is completed by a finite difference analysis.

AAA methods for field analysis prove to be strikingly fast and accurate. Results are in excellent agreement with those provided by the finite element model, and the very good agreement with those from finite differences suggests future improvements. They are also easy programming; the MATLAB code is less than 200 lines. This indicates they can provide an effective tool for rapid analysis.

AAA methods in magnetostatics are novel, but their extension to analogous physical problems seems straightforward. Being a meshless method, it is unlikely that local non-linearities can be considered. An aspect of particular interest, left for future research, is the capability of handling inhomogeneous domains, i.e. solving general interface problems.

The authors use cutting-edge mathematical tools for the modelling of complex physical objects in magnetostatics.

In recent times, there has been a growing interest in buoyancy-induced heat transfer within confined enclosures due to its frequent occurrence in heat transfer processes across diverse engineering disciplines, including electronic cooling, solar technologies, nuclear reactor systems, heat exchangers and energy storage systems. Moreover, the reduction of entropy generation holds significant importance in engineering applications, as it contributes to enhancing thermal system performance. This study, a numerical investigation, aims to analyze entropy generation and natural convection flow in an inclined square enclosure filled with Ag–MgO/water and Ag–TiO2/water hybrid nanofluids under the influence of a magnetic field. The enclosure features heated slits along its bottom and left walls. Following the Boussinesq approximation, the convective flow arises from a horizontal temperature difference between the partially heated walls and the cold right wall.

The governing equations for laminar unsteady natural convection flow in a Newtonian, incompressible mixture is solved using a Marker-and-Cell-based finite difference method within a customized MATLAB code. The hybrid nanofluid’s effective thermal conductivity and viscosity are determined using spherical nanoparticle correlations.

The numerical investigations cover various parameters, including nanoparticle volume concentration, Hartmann number, Rayleigh number, heat source/sink effects and inclination angle. As the Hartmann and Rayleigh numbers increase, there is a significant enhancement in entropy generation. The average Nusselt number experiences a substantial increase at extremely high values of the Rayleigh number and inclination.

This numerical investigation explores advanced applications involving various combinations of influential parameters, different nanoparticles, enclosure inclinations and improved designs. The goal is to control fluid flow and enhance heat transfer rates to meet the demands of the Fourth Industrial Revolution.

In a 90° tilted enclosure, the addition of 5% hybrid nanoparticles to the base fluid resulted in a 17.139% increase in the heat transfer rate for Ag–MgO nanoparticles and a 16.4185% increase for Ag–TiO2 nanoparticles compared to the base fluid. It is observed that a 5% nanoparticle volume fraction results in an increased heat transfer rate, influenced by variations in both the Darcy and Rayleigh numbers. The study demonstrates that the Ag–MgO hybrid nanofluid exhibits superior heat transfer and fluid transport performance compared to the Ag–TiO2 hybrid nanofluid. The simulations pertain to the use of hybrid magnetic nanofluids in fuel cells, solar cavity receivers and the processing of electromagnetic nanomaterials in enclosed environments.

Intumescent coatings are nowadays a dominant passive system used to protect structural materials in case of fire. Due to their reactive swelling behaviour, intumescent coatings are particularly complex materials to be modelled and predicted, which can be extremely useful especially for performance-based fire safety designs. In addition, many parameters influence their performance, and this challenges the definition and quantification of their material properties. Several approaches and models of various complexities are proposed in the literature, and they are reviewed and analysed in a critical literature review.

Analytical, finite-difference and finite-element methods for modelling intumescent coatings are compared, followed by the definition and quantification of the main physical, thermal, and optical properties of intumescent coatings: swelled thickness, thermal conductivity and resistance, density, specific heat capacity, and emissivity/absorptivity.

The study highlights the scarce consideration of key influencing factors on the material properties, and the tendency to simplify the problem into effective thermo-physical properties, such as effective thermal conductivity. As a conclusion, the literature review underlines the lack of homogenisation of modelling approaches and material properties, as well as the need for a universal modelling method that can generally simulate the performance of intumescent coatings, combine the large amount of published experimental data, and reliably produce fire-safe performance-based designs.

Due to their limited applicability, high complexity and little comparability, the presented literature review does not focus on analysing and comparing different multi-component models, constituted of many model-specific input parameters. On the contrary, the presented literature review compares various approaches, models and thermo-physical properties which primarily focusses on solving the heat transfer problem through swelling intumescent systems.

The presented literature review analyses and discusses the various modelling approaches to describe and predict the behaviour of swelling intumescent coatings as fire protection for structural materials. Due to the vast variety of available commercial products and potential testing conditions, these data are rarely compared and combined to achieve an overall understanding on the response of intumescent coatings as fire protection measure. The study highlights the lack of information and homogenisation of various modelling approaches, and it underlines the research needs about several aspects related to the intumescent coating behaviour modelling, also providing some useful suggestions for future studies.

The purpose of this paper is to analyze numerically the blowup in finite time of the solutions to a one-dimensional, bidirectional, nonlinear wave model equation for the propagation of small-amplitude waves in shallow water, as a function of the relaxation time, linear and nonlinear drift, power of the nonlinear advection flux, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of three types of initial conditions.

An implicit, first-order accurate in time, finite difference method valid for semipositive relaxation times has been used to solve the equation in a truncated domain for three different initial conditions, a first-order time derivative initially equal to zero and several constant wave speeds.

The numerical experiments show a very rapid transient from the initial conditions to the formation of a leading propagating wave, whose duration depends strongly on the shape, amplitude and width of the initial data as well as on the coefficients of the bidirectional equation. The blowup times for the triangular conditions have been found to be larger than those for the Gaussian ones, and the latter are larger than those for rectangular conditions, thus indicating that the blowup time decreases as the smoothness of the initial conditions decreases. The blowup time has also been found to decrease as the relaxation time, degree of nonlinearity, linear drift coefficient and amplitude of the initial conditions are increased, and as the width of the initial condition is decreased, but it increases as the viscosity coefficient is increased. No blowup has been observed for relaxation times smaller than one-hundredth, viscosity coefficients larger than ten-thousandths, quadratic and cubic nonlinearities, and initial Gaussian, triangular and rectangular conditions of unity amplitude.

The blowup of a one-dimensional, bidirectional equation that is a model for the propagation of waves in shallow water, longitudinal displacement in homogeneous viscoelastic bars, nerve conduction, nonlinear acoustics and heat transfer in very small devices and/or at very high transfer rates has been determined numerically as a function of the linear and nonlinear drift coefficients, power of the nonlinear drift, viscosity coefficient, viscous attenuation, and amplitude, smoothness and width of the initial conditions for nonzero relaxation times.

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 study is to design a closed hydrostatic guideway has the ability to resist large-side load, pitch moments and yaw moments, has good stiffness and damping characteristics, and provides certain beneficial guidance for the design of large-span closed hydrostatic guideway on the basis of providing a large vertical load bearing capacity.

The Reynolds’ equation and flow continuity equation are solved simultaneously by the finite difference method, and the perturbation method and the finite disturbance method is used for calculating the dynamic characteristics. The static and dynamic characteristics, including recess pressure, flow of lubricating oil, carrying capacity, pitch moment, yaw moment, dynamic stiffness and damping, are comprehensively analyzed.

The designed closed hydrostatic guideway has the ability to resist large lateral load, pitch moment and yaw moment and has good stiffness and damping characteristics, on the basis of being able to provide large vertical carrying capacity, which can meet the application requirements of heavy two-plate injection molding machine (TPIMM).

This paper researches static and dynamic characteristics of a large-span six-slider closed hydrostatic guideway used in heavy TPIMM, emphatically considering pitch moment and yaw moment. Some useful guidance is given for the design of large-span closed hydrostatic guideway.

This paper aims to study numerically the steady natural convective heat transfer of a hybrid nanosuspension (Ag-MgO/H2O) within a partially heated/cooled trapezoidal region with linear temperature profiles at inclined walls under an effect of uniform Lorentz force. This investigation is useful for researchers studying in the area of cavity flows to know features of the flow structures and nature of hybrid nanofluid characteristics. In addition, a detailed entropy generation analysis has been performed to highlight possible regimes with minimal entropy generation rates.

The governing equations formulated using the Oberbeck–Boussinesq approach and single-phase nanoliquid model are transformed to a non-dimensional form by using non-dimensional variables. The obtained equations with appropriate boundary conditions are resolved by the finite difference technique. The developed code has been validated comprehensively. Analysis has been performed for a wide range of governing parameters, including Rayleigh number (Ra = 105), Prandtl number (Pr = 6.82), Hartmann number (Ha = 0–100), magnetic field inclination angle (φ = 0–?/2) and nanoparticles volume fraction (φ_hnf = 0 and 2%).

It has been shown that inclined magnetic field can be used to manage the energy transport performance. An inclusion of nanoparticles without Lorentz force influence allows forming more stable convective regime with descending heat plume in the central zone, while such a regime was performed for clear fluid only for moderate and high Hartmann numbers. Moreover, the average overall entropy generation can be decreased with a growth of the Hartmann number, while an addition of hybrid nanoparticles allows reducing this parameter for Ha = 30 and 50. The average Nusselt number can be increased with a growth of the nanoparticles concentration for low values of the magnetic field intensity.

Governing equations written using the conservation laws and dimensionless non-primitive variables have been resolved by the finite difference approach. The created numerical code has been verified by applying the grid independence test and computational outcomes of other researchers. The comprehensive analysis for various key parameters has been performed.

This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

The purpose of this paper is to assess the feasibility of analytical models, specifically the radial basis function method, Akbari–Ganji method and Gaussian method, in conjunction with the finite element method. The aim is to examine the impact of processing parameters on temperature history.

Through analytical investigation and finite element simulation, this research examines the influence of processing parameters on temperature history. Simufact software with a thermomechanical approach was used for finite element simulation, while radial basis function, Akbari–Ganji and Gaussian methods were used for analytical modeling to solve the heat transfer differential equation.

The accuracy of both finite element and analytical methods was validated with about 90%. The findings revealed direct relationships between thermal conductivity (from 100 to 200), laser power (from 400 to 800 W), heat source depth (from 0.35 to 0.75) and power absorption coefficient (from 0.4 to 0.8). Increasing the values of these parameters led to higher temperature history. On the other hand, density (from 7,600 to 8,200), emission coefficient (from 0.5 to 0.7) and convective heat transfer (from 35 to 90) exhibited an inverse relationship with temperature history.

The application of analytical modeling, particularly the utilization of the Akbari–Ganji, radial basis functions and Gaussian methods, showcases an innovative approach to studying directed energy deposition. This analytical investigation offers an alternative to relying solely on experimental procedures, potentially saving time and resources in the optimization of DED processes.