Analytical and numerical analysis of time fractional dual-phase-lag heat conduction during short-pulse laser heating

Numerical Algorithms - In this study, we analytically and numerically investigate the non-Fourier heat conduction behavior within a finite medium based on the time fractional dual-phase-lag model....     

Thermoelastic damping effect analysis in micro flexural resonator of atomic force microscopy

In the design of high-Q micro/nano-resonators, dissipation mechanisms may have damaging effects on the quality factor (Q). One of the major dissipation mechanisms is thermoelastic damping (TED) that needs an accurate consideration for prediction. Aim of this paper is to evaluate the effect of TED on the vibrations of thin beam resonators. In particular, we will focus on cantilever beam resonator used in atomic force microscopy (AFM). AFM resonator is actually a cantilever with a spring attached to its free end. The end spring is considered to capture the effect of surface stiffness between tip and sample surface. The coupled governing equations of motion of thin beam with consideration of TED effects are derived. In general, there are four elastic equations that are coupled with thermal conduction equation. Based on accurate assumptions, these equations are simplified and the various boundary conditions have been used in order to validate the computational procedure. In order to accurately determine TED effects, the coupled thermal conduction equation is solved for the temperature field by considering three-dimensional (3-D) heat conduction along the length, width and thickness of the beam. Weighted residual Galerkin technique is used to obtain frequency shift and the quality factor of the thin beam resonator. The obtained results for quality factor, frequency shift and sensitivity change due to thermo-elastic coupling are presented graphically. Furthermore, the effects of beam aspect ratio, stress-free temperature on the quality factor and the influence of the surface stiffness on the frequencies and modal sensitivity of the AFM cantilever with and without considering thermo-elastic damping effects are discussed.     

Thermoelastic damping in flexural vibration of bilayered microbeams with circular cross-section

Predicting of thermoelastic damping (TED) is crucial in the design of micro-resonators with composite structures. Several analytical models were developed to evaluate TED in bilayered and three-layered microbeams in the past. However, the previous models focus on the microbeams with rectangular cross-section. This paper aims to study the TED in a bilayered microbeam with circular cross-section. The temperature field is approximated by using sine series and Bessel series in the circular cross-section. An analytical full model for TED in flexural vibration of bilayered microbeam is presented in the form of an infinite series. A simple model is also developed by retaining only the first term. The simulation results of the present model have a good agreement with those of the finite element method (FEM). The results indicate that well-resolved two peaks of TED are typically observed in the bilayered microbeams in which the value of k₂/C₂ is much less or higher than that of k₁/C₁. The present model is not a rapidly converging infinite series for most of the bilayered microbeams. The present model is more suitable for the long slender beams.     

Efficient and accurate spectral method for the time-fractional dual-phase-lag heat transfer model and its parameter estimation

In the current paper, a heat transfer model is suggested based on a time-fractional dual-phase-lag (DPL) model. We discuss the model in two parts, the direct problem and the inverse problem. Firstly, for solving it, the finite difference/Legendre spectral method is constructed. In the temporal direction, we employ the weighted and shifted Grünwald approximation, which can achieve second order convergence. In the spatial direction, the Legendre spectral method is used, it can obtain spectral accuracy. The stability and convergence are theoretically analyzed. For the inverse problem, the Bayesian method is used to construct an algorithm to estimate the four parameters for the model, namely, the time-fractional order α, the time-fractional order β, the delay time τ_T, and the relaxation time τ_q. Next, numerical experiments are provided to demonstrate the effectiveness of our scheme, with the values of τ_q and τ_T for processed meat employed. We also make a comparison with another method. The data obtained for the direct problem are used in the parameter estimation. The paper provides an accurate and efficient numerical method for the time-fractional DPL model.     

Dual-phase-lag analysis of CNT–MoS2–ZrO2–SiO2–Si nano-transistor and arteriole in multi-layered skin

The numerical simulation of non-Fourier dual-phase-lag heat conduction in three-dimensional bio and nano media has been investigated with two developed solvers using the OpenFOAM C++ finite-volume libraries. Despite the rapid growth in the use of non-Fourier heat transfer model, most researchers conduct simulations over simple geometries due to numerical restrictions. The developed solvers have been verified for four benchmark cases. After that, two complex geometries have been selected to show the accuracy of the solvers for complicated geometries. The first case is the temperature distribution in the human hand skin in the presence of a circular arteriole, while the second case is the thermal analysis of a newly constructed CNT–MoS₂ transistor with 1 nm gate length. Both cases are multilayered structures with curved boundaries. The computed temperature rises in the blood vessel and in the transistor are 63°C after 15 s and 16 °C after 10 ps, respectively. It is found that considering the temperature-dependent electrical properties of ZrO₂ and CNT is essential to present an accurate mathematical model.     

考虑记忆效应及尺寸效应窄长薄板的磁-热弹性耦合动态响应

引入记忆依赖微分的双相滞后热弹性理论能较完善地描述非Fourier导热现象,然而迄今尚未发现该理论综合考虑微尺度效应和磁、热、弹等多场耦合效应对材料力学行为的影响.通过考虑记忆依赖效应和非局部效应修正了双相滞后广义热弹性理论,基于改进后的理论研究了受周期性变化热源作用时窄长薄板的磁-热弹性耦合问题.首先建立问题的控制方程;然后结合边界条件与初值条件,利用Laplace变换和反变换技术对该问题进行求解;最后分别考察了磁场、相位滞后、时间延迟因子、核函数、非局部效应、时间对各无量纲量的影响,为微尺度材料的动态响应提供了有力参考依据.     

Finite integral transform-based analytical solutions of dual phase lag bio-heat transfer equation

A finite integral transform (FIT)-based analytical solution to the dual phase lag (DPL) bio-heat transfer equation has been developed. One of the potential applications of this analytical approach is in the field of photo-thermal therapy, wherein the interest lies in determining the thermal response of laser-irradiated biological samples. In order to demonstrate the applicability of the generalized analytical solutions, three problems have been formulated: (1) time independent boundary conditions (constant surface temperature heating), (2) time dependent boundary conditions (medium subjected to sinusoidal surface heating), and (3) biological tissue phantoms subjected to short-pulse laser irradiation. In the context of the case study involving biological tissue phantoms, the FIT-based analytical solutions of Fourier, as well as non-Fourier, heat conduction equations have been coupled with a numerical solution of the transient form of the radiative transfer equation (RTE) to determine the resultant temperature distribution. Performance of the FIT-based approach has been assessed by comparing the results of the present study with those reported in the literature. A comparison of DPL-based analytical solutions with those obtained using the conventional Fourier and hyperbolic heat conduction models has been presented. The relative influence of relaxation times associated with the temperature gradients (τ_T) and heat flux (τ_q) on the resultant thermal profiles has also been discussed. To the best of the knowledge of the authors, the present study is the first successful attempt at developing complete FIT-based analytical solution(s) of non-Fourier heat conduction equation(s), which have subsequently been coupled with numerical solutions of the transient form of the RTE. The work finds its importance in a range of areas such as material processing, photo-thermal therapy, etc.     

Resonance and chaos of micro and nano electro mechanical resonators with time delay feedback

The resonance and chaos of micro (nano) electro mechanical resonators with time delay feedback is concerned in the paper. Based on the experimental results, a lumped single degree-of-freedom (1DOF) model is studied and the effects of time delay displacement and velocity feedback on the system are investigated. In order to have a deep insight into the system, the amplitude frequency response curve of the system is firstly obtained using the multiple scales method. The Melnikov function method is then extended to the two time delay systems, and the analytically required condition for chaos was obtained. Finally, the fourth-order Runge–Kutta method, point-mapping method and spectrum diagram are used to simulate the evolution of the dynamic behavior of the time delay control system. Also, the stability of this time delay control system is studied thoroughly. The results show that time delay feedback is a good method for the control system and that reasonable selection of control system parameters can effectively suppress the vibration level for micro/nano-electro-mechanical resonator systems.     

Vibration of nanobeams of different boundary conditions with multiple cracks based on nonlocal elasticity theory

The aim of this paper is to study the free vibration of nanobeams with multiple cracks. The analysis procedure is based on nonlocal elasticity theory. This theory states that stress at a point is a function of strains at all points in the continuum. The nonlocal elasticity theory becomes significant for small length scale in micro and nanostructures. The effects of nonlocality, crack location and crack parameter are investigated on the natural frequencies of the cracked nanobeam. In this study, analytical solutions are given for cracked Euler–Bernoulli nanobeams of different boundary conditions.     

Size-dependent dynamics of a FG Nanobeam near nonlinear resonances induced by heat

This study explores heat-induced nonlinear vibration of a functionally graded (FG) capacitive nanobeam within the framework of nonlocal strain gradient theory (NLSGT). The elastic FG beam, which is firstly deflected by a DC voltage, is driven to vibrate about its deflected position by a periodic heat load. The nano-structure, which consists of a clamped-clamped nanobeam, is modeled assuming Euler–Bernoulli beam assumption which accounts for the nonlinear von-Karman strain and the electrostatic and intermolecular forcing. To simulate the static and dynamic responses, a model reduction procedure is carried out by employing the Galerkin method. The method of Averaging as a regular semi-analytic perturbation method is applied to obtain governing equations of the steady-state responses. With the purpose of establishing the validity of the solution, a Shooting technique in conjunction with the Floquet theory is used to capture the periodic motions and then examine their stability. The nonlinear resonance frequency of the FG nanobeam near its fundamental natural frequency (primary resonance) and near principal parametric resonance is investigated while the emphasis is placed on studying the effect of various parameters including DC voltage, amplitude of the periodic heat source, material index, damping ratio, and small scale parameters. The main objective of this study is to model a miniature structure which can be used as either a sensitive remote temperature sensor or a high-efficiency thermal energy harvester.     

Surface effects on nonlinear free vibration of functionally graded nanobeams using nonlocal elasticity

In this paper, to consider all surface effects including surface elasticity, surface stress, and surface density, on the nonlinear free vibration analysis of simply-supported functionally graded Euler–Bernoulli nanobeams using nonlocal elasticity theory, the balance conditions between FG nanobeam bulk and its surfaces are considered to be satisfied assuming a cubic variation for the component of the normal stress through the FG nanobeam thickness. The nonlinear governing equation includes the von Kármán geometric nonlinearity and the material properties change continuously through the thickness of the FG nanobeam according to a power-law distribution of the volume fraction of the constituents. The multiple scale method is employed as an analytical solution for the nonlinear governing equation to obtain the nonlinear natural frequencies of FG nanobeams. The effect of the gradient index, the nanobeam length, thickness to length ratio, mode number, amplitude of deflection to radius of gyration ratio and nonlocal parameter on the frequency ratios of FG nanobeams is investigated.     

Non-linear forced vibration analysis of nanobeams subjected to moving concentrated load resting on a viscoelastic foundation considering thermal and surface effects

Analytical solution for the steady-state response of an Euler–Bernoulli nanobeam subjected to moving concentrated load and resting on a viscoelastic foundation with surface effects consideration in a thermal environment is investigated in this article. At first, based on the Eringen's nonlocal theory, the governing equations of motion are derived using the Hamilton's principle. Then, in order to solve the equation, Galerkin method is applied to discretize the governing nonlinear partial differential equation to a nonlinear ordinary differential equation; solution is obtained employing the perturbation technique (multiple scales method). Results indicate that by increasing of various parameters such as foundation damping, linear stiffness, residual surface stress and the temperature change, the jump phenomenon is postponed and with increasing the amplitude of the moving force and the nonlocal parameter, the jump phenomenon occurs earlier and its frequency and the peak value of amplitude of vibration increases. In addition, it is seen that the non-linear stiffness and the critical velocity of the moving load are important factors in studying nanobeams subjected to moving concentrated load. Presence of the non-linear stiffness of Winkler foundation resulting nanobeam tends to instability and so, the jump phenomenon occurs. But, presence of the linear stiffness will lead to stability of the nanobeam. In the next sections of the paper, frequency responses of the nanobeam made of temperature-dependent material properties under multi-frequency excitations are investigated.     

Geometrically nonlinear analysis of Timoshenko piezoelectric nanobeams with flexoelectricity effect based on Eringen's differential model

This study analyzes the nonlinear free vibration and post-buckling of nanobeams with flexoelectric effect based on Eringen's differential model. The nanobeam is modeled based on Timoshenko beam's theory. The von-Kármán strain–displacement relation together with the electrical Gibbs free energy and Hamilton's principle are employed to derive equations of motion. The nonlinear free vibration frequencies are obtained for pinned–pinned (P–P) and clamped–clamped (C–C) boundary conditions. Multiple scales method is employed to obtain the closed-form solution for the nonlinear governing equations. By employing this methodology, the natural frequencies of nanobeams are obtained and their post-buckling behavior is examined. The influence of nonlocal parameter, amplitude ratio, and input voltage on the top surface and flexoelectricity constant on nonlinear free vibration and post-buckling characteristics of nanobeam is investigated. In this paper, it is concluded that the flexoelectricity has a significant effect on free vibration of the beams in nano-scale and its effect has to be considered in designing nano-electro-mechanical systems (NEMS) such as nano- generators and nano-sensors.     

非Fourier温度场分布的奇摄动解

应用非Fourier热传导定律构建了单层材料中温度场模型，即一类在无界域上带小参数的奇摄动双曲方程，通过奇摄动展开方法，得到了该问题的渐近解.首先应用奇摄动方法得到了该问题的外解和边界层矫正项，通过对内解和外解的最大模估计和关于时间导数的最大模估计以及线性抛物方程理论，得到了内外解的存在唯一性，从而得到了解的形式渐近展开式.通过余项估计，给出了渐近解的L2估计，得到了渐近解的一致有效性，从而得到了无界域上温度场的分布.通过奇摄动分析，给出了非Fourier 温度场与Fourier 温度场的关系，描述了非Fourier温度场的具体形态.     

Modelling non-Fourier heat conduction with periodic thermal oscillation using the finite integral transform

The non-Fourier heat conduction in a finite medium subjected to a periodic heat flux is modelled using the finite integral transform technique and an analytic solution is obtained. An analogy between thermal oscillation and oscillation of mechanical and electrical systems is drawn. A transition criterion from the non-Fourier heat conduction formulation to the Fourier formulation is obtained and a simple analytical expression of the phase and amplitude of thermal oscillation is derived.     

Long-time dynamics of an extensible plate equation with thermal memory

This work is concerned with long-time dynamics of solutions of extensible plate equations with thermal memory. The problem corresponds to a model of thermoelasticity based on a theory of non-Fourier heat flux. By considering the case where rotational inertia is present we show that the thermal dissipation is sufficient to stabilize the system and guarantees the existence of a finite-dimensional global attractor. In addition, the existence of exponential attractors is also considered.     

A superstatistical model for anomalous heat conduction and diffusion

We propose a superstatistical model for anomalous heat conduction and diffusion, which is formulated by the thermal conductivity distribution, overall temperature and heat flux distributions. Our model obeys Fourier's law and the continuity equation at the individual level. The evolution of the thermal conductivity distribution is described by an advection-diffusion equation. We show that the superstatistical model predict anomalous behaviors including the time-dependent effective thermal conductivity and slow long-time asymptotics. The time-dependence of the effective thermal conductivity is determined by the mean square displacement (MSD), which coincides with existing investigations. The superstatistical structure can also be extended into other non-Fourier models including the Cattaneo and fractional-order heat conduction models.     

Finite element modeling of the nonlinear oscillations in axisymmetric acoustic resonators

Oscillation of a gas in closed resonators has gained considerable interest in the past years. In this paper, the nonlinear equations governing the behavior of the gas oscillations inside the resonator are formulated in a weak form and then modeled using the finite element method. The pressure ratios, predicted by the proposed model, are in close agreement with the exact solutions available for simple geometries such as cylindrical, exponential and linearly varying area resonators. The presented comparisons validate the accuracy of the finite element model and emphasize its potential for predicting the performance or resonators of more complex geometries which are necessary for generating high pressures from the standing waves. Also, gas flow through the boundaries of the resonator is implemented in the proposed model. The presented finite element model presents an invaluable tool for designing a new class of acoustic compressors which can be used, for example, in refrigeration and vibration control applications.     

Numerical Algorithm with Fourth-Order Spatial Accuracy for Solving the Time-Fractional Dual-Phase-Lagging Nanoscale Heat Conduction Equation

Nanoscale heat transfer cannot be described by the classical Fourier law due to the very small dimension, and therefore, analyzing heat transfer in nanoscale is of crucial importance for the design and operation of nano-devices and the optimization of thermal processing of nano-materials. Recently, time-fractional dual-phase-lagging (DPL) equations with temperature jump boundary conditions have showed promising for analyzing the heat conduction in nanoscale. This article proposes a numerical algorithm with high spatial accuracy for solving the time-fractional dual-phase-lagging nano-heat conduction equation with temperature jump boundary conditions. To this end, we first develop a fourth-order accurate and unconditionally stable compact finite difference scheme for solving this time-fractional DPL model. We then present a fast numerical solver based on the divide-and-conquer strategy for the obtained finite difference scheme in order to reduce the huge computational work and storage. Finally, the algorithm is tested by two examples to verify the accuracy of the scheme and computational speed. And we apply the numerical algorithm for predicting the temperature rise in a nano-scale silicon thin film. Numerical results confirm that the present difference scheme provides ${\rm min}\{2−α, 2−β\}$ order accuracy in time and fourth-order accuracy in space, which coincides with the theoretical analysis. Results indicate that the mentioned time-fractional DPL model could be a tool for investigating the thermal analysis in a simple nanoscale semiconductor silicon device by choosing the suitable fractional order of Caputo derivative and the parameters in the model.     

Buckling and stability analysis of a piezoelectric viscoelastic nanobeam subjected to van der Waals forces

A study on the buckling and dynamic stability of a piezoelectric viscoelastic nanobeam subjected to van der Waals forces is performed in this research. The static and dynamic governing equations of the nanobeam are established with Galerkin method and under Euler–Bernoulli hypothesis. The buckling, post-buckling and nonlinear dynamic stability character of the nanobeam is presented. The quasi-elastic method, Leibnitz’s rule, Runge–Kutta method and the incremental harmonic balanced method are employed for obtaining the buckling voltage, post-buckling characteristics and the boundaries of the principal instability region of the dynamic system. Effects of the electrostatic load, van der Waals force, creep quantity, inner damping, geometric nonlinearity and other factors on the post-buckling and the principal region of instability are investigated.