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.

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.

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.

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.

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.

This paper aims to develop a model that reflects the current transformer (CT) core materials nonlinearity. The model enables simulation and analysis of the CT excitation current that includes the inductive magnetizing current and the resistive excitation current.

A nonlinear CT model is established with the magnetizing current as the solution variable. This model presents the form of a nonlinear differential equation and can be solved discretely using the Runge–Kutta method.

By simulating variations in the excitation current for different primary currents, loads and core materials, the results demonstrate that enhancing the permeability of the B–H curve leads to a more significant improvement in the CT ratio error at low primary currents.

The proposed model has three obvious advantages over the previous models with the secondary current as the solution variable: (1) The differential equation is simpler and easier to solve. Previous models contain the time differential terms of the secondary current and excitation flux or the integral term of the flux, making the iterative solution complicated. The proposed model only contains the time differential of the magnetizing current. (2) The accuracy of the excitation current obtained by the proposed model is higher. Previous models calculate the excitation current by subtracting the secondary current from the converted primary current. Because these two currents are much greater than the excitation current, the error of calculating the small excitation current by subtracting two large numbers is greatly enlarged. (3) The proposed model can calculate the distorted waveform of the excitation current and error for any form of time-domain primary current, while previous models can only obtain the effective value.

The high-pressure accumulator has been widely used in the hydraulic system. Failure pressure prediction is crucial for the safe design and integrity assessment of the accumulators. The purpose of this study is to accurately predict the burst pressure and location for the accumulator shells due to internal pressure.

This study concentrates the non-linear finite element simulation procedure, which allows determination of the burst pressure and crack location using extensive plastic straining criterion. Meanwhile, the full-scale hydraulic burst test and the analytical solution are conducted for comparative analysis.

A good agreement between predicted and measured the burst pressure that was obtained, and the predicted failure point coincided very well with the fracture location of the actual shell very well. Meanwhile, the burst pressure of the shells increases with wall thickness, independent of the length. It can be said that the non-linear finite element method can be employed to predict the failure behavior of a cylindrical shell with sufficient accuracy.

This paper can provide a designer with additional insight into how the pressurized hollow cylinder might fail, and the failure pressure has been predicted accurately with a minimum error below 1%, comparing the numerical results with experimental data.

As to different angles of attack and nonlinear problems caused by high temperatures in coexisting hypersonic aircraft, people mainly rely on fluid software for research but lack analysis of flow mechanisms. Owing to computational difficulties, few people use numerical algorithms to combine them for discussion. Hence, this study aims to make a deep inquiry into the laminar flow and heat transfer of compressible Newtonian fluid in hypersonic aircraft with small attack angles.

In this paper, on the basis of mass, momentum and energy conservation laws, the governing equations of the hypersonic boundary layer are established. Viscosity, specific heat capacity and thermal conductivity are considered nonlinear functions concerning temperature. In virtue of the MacCormack finite difference method, the stationary numerical solutions are solved directly, and the validity of the algorithm is verified.

The results demonstrate that at Mach number 5, compared to the 0° attack angle, the maximum temperature near-wall at the 3° attack angle increases by about 25%. An enjoyable phenomenon is discovered, where the position corresponding to the maximum wall shear force shifts back as the attack angle and Mach number increase. The relationship between the near-wall maximum temperature versus attack angle and Mach number is fitted through numerical calculation results.

Empirical formulas can be used to estimate heat transfer characteristics at small attack angles, which will guide the design of aircraft thermal protection systems.