Study of thermal conductivity and thermal rectification in exponential mass graded lattices

Concept of exponential mass variation of oscillators along the chain length of N oscillators is proposed in the present Letter. The temperature profile and thermal conductivity of one-dimensional (1D) exponential mass graded harmonic and anharmonic lattices are studied on the basis of Fermi-Pasta-Ulam (FPU) β model. Present findings conclude that the exponential mass graded chain provide higher conductivity than that of linear mass graded chain. The exponential mass graded anharmonic chain generates the thermal rectification of 70-75% which is better than linear mass graded materials, so far. Thus instead of using linear mass graded material, the use of exponential mass graded material will be a better and genuine choice for controlling the heat flow at nano-scale.     

Magnetically induced thermal rectification

We consider far from equilibrium heat transport in chaotic billiard chains with noninteracting charged particles in the presence of nonuniform transverse magnetic field. If half of the chain is placed in a strong magnetic field, or if the strength of the magnetic field has a large gradient along the chain, heat current is shown to be asymmetric with respect to exchange of the temperatures of the heat baths. Thermal rectification factor can be arbitrarily large for sufficiently small temperature of one of the baths.     

Nontrivial properties of heat flow: Analytical study of some anharmonic lattice microscopic models

We address the analytical investigation of nontrivial properties of heat flow, starting from microscopic models of matter. We present an integral representation for the expression of the heat flow, by taking as our crystal model a self-consistent chain of anharmonic oscillators, precisely, a chain of oscillators with harmonic interparticle interactions, anharmonic on-site potentials, thermal reservoirs at the boundaries, and still with some residual inner stochastic baths. We propose an approximative scheme to make the integral formalism analytically treatable: we test the approximations in harmonic models and analyze some anharmonic systems. For the case of graded anharmonic models with weak interparticle interactions, we derive an expression for the thermal conductivity, and show the existence of thermal rectification and also nonnegative differential thermal resistance. The details of our formalism allow us to note the ingredients behind these phenomena, and the generality of our results (i.e., the results will be valid for other interactions and regimes) shows that rectification in graded materials is a ubiquitous property.     

FGMs稳态热传导分析的重心Lagrange插值配点法

提出一种求解二维功能梯度材料（FGMs）稳态热传导问题的重心Lagrange插值配点法.基于Chebyshev节点构造二维重心Lagrange插值函数及其偏导数,然后基于配点法将其直接代入FGMs热传导问题的控制方程和边界条件,得到系统离散方程.重心Lagrange插值配点法是一种真正的无网格方法,很好地融合了重心Lagrange插值和配点格式的优势,具有高效、稳定、高精度和易于数值实现的优点.采用重心Lagrange插值配点法分别对指数型、二次型和三角型FGMs热传导问题进行数值模拟.结果表明：该方法具有较高的计算效率和计算精度,对材料梯度参数的变化不敏感.可以进一步拓展到FGMs瞬态问题和FGMs的热力耦合分析.     

考虑界面接触热阻的一维复合结构的热整流机理

建立了考虑变截面、变热导率及界面接触热阻效应的组合热整流结构的温度场及热整流系数的理论模型和有限元解.数值算例证明了本文模型及算法的可靠性,进而通过参数影响研究确定了若干几何及材料参数对结构热整流系数的影响规律,揭示界面接触热阻对热整流效果的影响机理.研究结果表明长度比、截面半径变化率、热导率、边界条件温差和界面接触热阻等因素必须通过优化设计才能得到最大的热整流系数,同时界面接触热阻的引入也为调控热整流系数提供了一条新的途径.     

Three-dimensional analysis for transient coupled thermoelastic response of a functionally graded rectangular plate

The generalized coupled thermoelasticity based on the Lord-Shulman theory is considered to study the transient thermoelastic response of functionally graded rectangular plates. The state equations of functionally graded rectangular plate subjected to time-dependent thermal loads were established by using of state space approach, in which three displacement components, three stress components, the temperature and the heat flux were chosen as state variables. By giving simply supported boundary conditions and assuming that the material properties of the plate have an exponential law distribution along the thickness-coordinate, the equations were solved by the numerical Laplace transformations and shooting methods for transient thermal responses of a three dimensional functionally graded rectangular plate due to a thermal shock on its top surface. Effects of the volume fraction distributions of material constituents on the thermal responses, including the temperature change, the displacement and the stresses distributions were investigated.     

Negative differential thermal resistance in a double-chain system

Using nonequilibrium molecular dynamics simulations, we study the phenomenon of negative differential thermal resistance (NDTR) in a double-chain system. We investigate the dependence of NDTR on the external potential, inter- and intra-chain interaction and the system size. It is reported that the NDTR can occur in a small double-chain system with weak external potential and weak inter- and intra-chain interaction. We also present the influence of the external potential, inter- and intra-chain interaction and the system size on the heat current of the system through the phonon spectral analysis.     

A general modelling method for functionally graded materials with an arbitrarily oriented crack

In order to analytically solve crack problems regarding functionally graded materials (FGMs), some ideal assumptions are often made. They are: (1) the properties of FGMs are usually assumed to be described by very particular functions; (2) the crack is assumed to be vertical to (or parallel to) the gradient of FGMs. However, these assumptions may not be practical for actual FGMs. Since the controlling differential equations with general mechanical properties are very difficult to solve and the arbitrarily oriented crack causes great trouble in the analytical procedure, a general piecewise-exponential model (GPE model) is proposed to investigate the fracture behaviour of FGMs with general mechanical properties and an arbitrarily oriented crack. “General mechanical properties” means that the mechanical properties in the GPE model are not required to be very particularly pre-defined functions but arbitrary functions determined by fitting the experimental results of FGMs. The studied FGMs are divided into some sub-layers with each layer’s properties varying exponentially so that the general mechanical properties can be approximated by a series of exponential functions and hence the stresses and displacements of each layer which may contain a mixed-mode crack can be solved analytically. By use of integral transform methods, principle of superposition, residual theorem and theory of singular integral equations, the mixed-mode crack problem can be turned into solving a group of singular integral equations from which mixed-mode stress intensity factors (SIFs) can be obtained. Finally, the influences of the nonhomogeneous and geometric parameters on the mixed-mode SIFs are analysed.     

A meshless model for transient heat conduction analyses of 3D axisymmetric functionally graded solids

A meshless numerical model is developed for analyzing transient heat conductions in three-dimensional （3D） axisymmetric continuously nonhomogeneous functionally graded materials （FGMs）. Axial symmetry of geometry and boundary conditions reduces the original 3D initial-boundary value problem into a two-dimensional （2D） problem. Local weak forms are derived for small polygonal sub-domains which surround nodal points distributed over the cross section. In order to simplify the treatment of the essential boundary conditions, spatial variations of the temperature and heat flux at discrete time instants are interpolated by the natural neighbor interpolation. Moreover, the using of three-node triangular finite element method （FEM） shape functions as test functions reduces the orders of integrands involved in domain integrals. The semi-discrete heat conduction equation is solved numerically with the traditional two-point difference technique in the time domain. Two numerical examples are investigated and excellent results are obtained, demonstrating the potential application of the proposed approach.     

热质的运动与传递-微尺度导热中的热质动能效应

基于热质(热量的当量动质量)的概念,通过建立和分析热质的运动方程得到了反映热质动能变化的稳态导热微分方程,表明Fourier导热定律只有在热质的动能变化相对热质势能变化很小而可以忽略时才成立;在高热流密度和低温的情况下热质的动能变化不可忽略,这种动能效应表现为热流密度和温度梯度不再成线性关系.动能效应也导致Fourier导热定律不能通过热流和温度梯度准确地获得物体的导热系数,本文基于热质运动方程给出了导热系数动能效应的修正式.最后针对高热流密度和低温一维稳态导热进行了分子动力学模拟验证.     

Buckling and free vibration analysis of functionally graded cylindrical shells subjected to a temperature-specified boundary condition

Linear thermal buckling and free vibration analysis are presented for functionally graded cylindrical shells with clamped-clamped boundary condition based on temperature-dependent material properties. The material properties of functionally graded materials (FGM) shell are assumed to vary smoothly and continuously across the thickness. With high-temperature specified on the inner surface of the FGM shell and outer surface at ambient temperature, 1D heat conduction equation along the thickness of the shell is applied to determine the temperature distribution; thereby, the material properties based on temperature distribution are made available for thermal buckling and free vibration analysis. First-order shear deformation theory along with Fourier series expansion of the displacement variables in the circumferential direction are used to model the FGM shell. Numerical studies involved the understanding of the influence of the power-law index, r/h and l/r ratios on the critical buckling temperature. Free vibration studies of FGM shells under elevated temperature show that the fall in natural frequency is very drastic for the mode corresponding to the lowest natural frequency when compared to the lowest buckling temperature mode.     

Exponential DNA replication by laminar convection

It is shown that laminar thermal convection can drive a chain reaction of DNA replication. The convection is triggered by a constant horizontal temperature gradient, moving molecules along stationary paths between hot and cold regions. This implements the temperature cycling for the classical polymerase chain reaction (PCR). The amplification is shown to be exponential and reaches 100,000-fold gains within 25 min. Besides direct applications, the mechanism might have implications for the molecular evolution of life.     

Periodic composites: quasi-uniform heat conduction,Janus thermal illusion,and illusion thermal diodes

Manipulating thermal conductivities at will plays a crucial role in controlling heat flow. By developing an effective medium theory including periodicity, here we experimentally show that nonuniform media can exhibit quasi-uniform heat conduction. This provides capabilities in proposing Janus thermal illusion and illusion thermal rectification. For the former, we study, via experiment and theory, a big periodic composite containing a small periodic composite with circular or elliptic particles. As a result, we reveal the Janus thermal illusion that describes the whole periodic system with both invisibility illusion along one direction and visibility illusion along the perpendicular direction, which is fundamentally different from the existing thermal illusions for misleading thermal detection. Further, the Janus illusion helps to design two different periodic systems that both work as thermal diodes but with nearly the same temperature distribution, heat fluxes and rectification ratios, thus being called illusion thermal diodes. Such thermal diodes differ from those extensively studied in the literature, and are useful for the areas that require both thermal rectification and thermal camouflage. This work not only opens a door for designing novel periodic composites in thermal camouflage and heat rectification, but also holds for achieving similar composites in other disciplines like electrostatics, magnetostatics, and particle dynamics.     

冲击相边界传播过程中梯度材料的形成

 采用一种基于简单混合物模型的本构，研究了不可逆相变材料中一维的冲击相边界的传播规律，发现在突加载荷并连续卸载的应力边界条件下，样品内部有梯度材料形成。分析了边界条件对材料梯度变化的影响，并且用基于特征线的数值方法计算了不同的边界条件下材料的梯度分布，从而提出利用冲击相边界的传播制备梯度材料的可能性。最后对这种方法的利弊做了探讨。     

The multiple scattering of non-homogeneous shear waves from two cavities in functionally graded materials

The propagation and multiple scattering of non-homogeneous shear waves resulting from two cavities embedded in exponential functional graded materials (FGMs) were investigated, and the dynamic stress around the two cavities derived. The non-homogeneous scattering fields of shear waves around the cavities are analytically expressed by using the wave function expansion method. The interaction of non-homogeneous scattering fields between the two cavities is described accurately. The dynamic stresses around the cavities under different geometrical and physical parameters are graphically illustrated and analyzed. Analysis shows that the non-homogeneous properties of FGMs exhibit a significant effect on the dynamic stress around the cavities. The effect of the non-homogeneous properties of FGMs on the dynamic stress is also dependent on the gradation direction, the distance between the two cavities, the relative position of the two cavities and the incident frequency of waves. A comparison with other existing studies in the literature is also presented.     

一维均匀Morse晶格体系的热流棘齿效应

本文系统研究了系统两端无平均温差时一维均匀Morse晶格中的热流棘齿效应. Morse晶格的两端分别与两个热浴相接触, 其中一端热浴温度周期调制,另一端热浴温度保持不变, 两端热浴温度长时平均相等. 数值结果表明, 当对一端热浴温度进行周期调制时, 系统中便会有稳定的定向热流产生. 通过改变调制频率和强度, 可以控制热流的大小及方向. 在合适的频率范围内, 可观察到一种非常有趣的现象——非定态负热导现象, 即系统中产生的定向热流逆着系统温度梯度方向由低温端流向高温端. 通过热波动力学分析(分析热流及温度分     

利用Taylor展开法研究Lamb波在功能梯度材料中的传播特性

超声技术可用于对功能梯度材料(FGMs)的性质进行评估. 由于FGMs性质的非均匀性,采用分布函数来描述FGMs弹性常数和密度沿厚度方向的变化趋势,并提出利用Taylor展开的方法来解决分布函数为任意函数时的FGMs中Lamb波的传播问题. 利用本征函数展开法得到了铁基氧化铝FGMs中Lamb波的相速度色散曲线,讨论了材料性质分布对铁基氧化铝FGMs中Lamb波传播特性的影响. 为FGMs性质(沿板厚方向变化)的反演提供了理论依据. 关键词： 功能梯度材料 Lamb波 Legendre多项式 分布函数     

Equilibration and universal heat conduction in fermi-pasta-ulam chains

It is shown numerically that for Fermi-Pasta-Ulam (FPU) chains with alternating masses and heat baths at slightly different temperatures at the ends, the local temperature (LT) on small scales behaves paradoxically in steady state. This expands the long established problem of equilibration of FPU chains. A well-behaved LT appears to be achieved for equal mass chains; the thermal conductivity is shown to diverge with chain length N as N(1/3), relevant for the much debated question of the universality of one-dimensional heat conduction. The reason why earlier simulations have obtained systematically higher exponents is explained.     

Controlling Thermal Conduction by Graded Materials

Manipulating thermal conductivities are fundamentally important for controlling the conduction of heat at will. Thermal cloaks and concentrators, which have been extensively studied recently, are actually graded materials designed according to coordinate transformation approaches, and their effective thermal conductivity is equal to that of the host medium outside the cloak or concentrator. Here we attempt to investigate a more general problem: what is the effective thermal conductivity of graded materials? In particular, we perform a first-principles approach to the analytic exact results of effective thermal conductivities of materials possessing either power-law or linear gradation profiles. On the other hand, by solving Laplace's equation, we derive a differential equation for calculating the effective thermal conductivity of a material whose thermal conductivity varies along the radius with arbitrary gradation profiles. The two methods agree with each other for both external and internal heat sources, as confirmed by simulation and experiment. This work provides different methods for designing new thermal metamaterials (including thermal cloaks and concentrators), in order to control or manipulate the transfer of heat.     

Ergodicity and Energy Distributions for Some Boundary Driven Integrable Hamiltonian Chains

We consider systems of moving particles in 1-dimensional space interacting through energy storage sites. The ends of the systems are coupled to heat baths, and resulting steady states are studied. When the two heat baths are equal, an explicit formula for the (unique) equilibrium distribution is given. The bulk of the paper concerns nonequilibrium steady states, i.e., when the chain is coupled to two unequal heat baths. Rigorous results including ergodicity are proved. Numerical studies are carried out for two types of bath distributions. For chains driven by exponential baths, our main finding is that the system does not approach local thermodynamic equilibrium as system size tends to infinity. For bath distributions that are sharply peaked Gaussians, in spite of the near-integrable dynamics, transport properties are found to be more normal than expected.