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.

The space fractional PDEs (SFPDEs) play an important role in the fractional calculus field. Proposing a high-order, stable and flexible numerical procedure for solving SFPDEs is the main aim of most researchers. This paper devotes to developing a novel spectral algorithm to solve the FitzHugh–Nagumo models with space fractional derivatives.

The fractional derivative is defined based upon the Riesz derivative. First, a second-order finite difference formulation is used to approximate the time derivative. Then, the Jacobi spectral collocation method is employed to discrete the spatial variables. On the other hand, authors assume that the approximate solution is a linear combination of special polynomials which are obtained from the Jacobi polynomials, and also there exists Riesz fractional derivative based on the Jacobi polynomials. Also, a reduced order plan, such as proper orthogonal decomposition (POD) method, has been utilized.

A fast high-order numerical method to decrease the elapsed CPU time has been constructed for solving systems of space fractional PDEs.

The spectral collocation method is combined with the POD idea to solve the system of space-fractional PDEs. The numerical results are acceptable and efficient for the main mathematical model.

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.

This study aims to derive a novel spatial numerical method based on multidimensional local Taylor series representations for solving high-order advection-diffusion (AD) equations.

The parabolic AD equations are reduced to the nonhomogeneous elliptic system of partial differential equations by utilizing the Chebyshev spectral collocation method (ChSCM) in the temporal variable. The implicit-explicit local differential transform method (IELDTM) is constructed over two- and three-dimensional meshes using continuity equations of the neighbor representations with either explicit or implicit forms in related directions. The IELDTM yields an overdetermined or underdetermined system of algebraic equations solved in the least square sense.

The IELDTM has proven to have excellent convergence properties by experimentally illustrating both h-refinement and p-refinement outcomes. A distinctive feature of the IELDTM over the existing numerical techniques is optimizing the local spatial degrees of freedom. It has been proven that the IELDTM provides more accurate results with far fewer degrees of freedom than the finite difference, finite element and spectral methods.

This study shows the derivation, applicability and performance of the IELDTM for solving 2D and 3D advection-diffusion equations. It has been demonstrated that the IELDTM can be a competitive numerical method for addressing high-space dimensional-parabolic partial differential equations (PDEs) arising in various fields of science and engineering. The novel ChSCM-IELDTM hybridization has been proven to have distinct advantages, such as continuous utilization of time integration and optimized formulation of spatial approximations. Furthermore, the novel ChSCM-IELDTM hybridization can be adapted to address various other types of PDEs by modifying the theoretical derivation accordingly.

The purpose of this study is to develop an efficient model to solve the electromagnetic forward problem using a novel semi-analytical approach to compute the electromagnetic fields because of the presence of a scatterer.

The proposed model involves a novel formulation of a complete orthonormal set of radiating/nonradiating polarization currents. Furthermore, an integral equation-based representation is derived, and the appropriate boundary conditions are imposed to get the scattered electromagnetic field. An error term is introduced to evaluate the obtained results.

The proposed model was tested using several examples at different frequencies. The results of this study show that the novel representation exhibits fast convergence behavior and achieves highly accurate results, when compared to the results provided by the transmission line method.

The derived formulations presented in this study are significant in the electromagnetic forward modelling field because of the meaningful physical representation they provide. This is an important aspect that leads to precise calculation of electromagnetic fields for various applications.

This study aims to investigate the stability performance of partial journal bearings of 120° and 180° partial angles with micropolar lubricant.

To investigate the stability characteristics of partial journal bearing, a MATLAB source code is written. To solve the Reynolds’ equation, the finite element method is used. Stability performances of 120° and 180° partial journal bearings are computed for a wide range of non-dimensional micropolar fluid parameters and working eccentricities.

The presented results provide design data for stability parameters in terms of equivalent stiffness, whirl frequency ratio, critical mass and threshold speed of the rotor with respect to eccentricities and material size of the lubricant. The stability of 180° partial journal bearing is found to be higher than 120° partial journal bearing.

In open literature, it is rare to find the stability of a partial journal bearing lubricated with micropolar fluid. Very few researchers have studied the combined effect of eccentricities and micropolar lubricant parameters on the dynamic performance of such bearings. Hence, it is important to study the dynamic stability to explore the complete investigation of the performance of partial journal bearings with micropolar fluid.

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.

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.