Numerical investigation of the deformation characteristics and heat generation in pneumatic aircraft tires : Part I. Mechanical modeling

A numerical methodology to study the temperature rise in aircraft tires under free rolling conditions is presented in this article. In the first part, we study the deformation characteristics of the tire to determine the heat generation due to the inelastic deformation (viscoelasticity), which is assumed to be the main source of heat generation under free rolling conditions. The heat generation is then used as an input to solve the heat transfer problem which is addressed in the second part of this article. A methodology which considers a 2-D formulation with the contribution of out-of-plane forces is presented. This methodology allows for a significant reduction of the computational requirements of 3-D analysis while capturing the main 3-D effects. The deformation pattern as well as the stress and strain fields are presented for a cross-section of the tire at several locations.     

Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical evaluation of the temperature field of steady-state rolling tires

Rubber is the main component of pneumatic tires. The tire heating is caused by the hysteresis effects due to the deformation of the rubber during operation. Tire temperatures can depend on many factors, including tire geometry, inflation pressure, vehicle load and speed, road type and temperature and environmental conditions. The focus of this study is to develop a finite element approach to computationally evaluate the temperature field of a steady-state rolling tire. For simplicity, the tire is assumed to be composed of rubber and body-ply. The nonlinear mechanical behavior of the rubber is characterized by a Mooney–Rivlin model while the body-ply is assumed to be linear elastic material. The coupled effects of the inflation pressure and vehicle loading are investigated. The influences of body-ply stiffness are studied as well. The simulation results show that loading is the main factor to determine the temperature field. The stiffer body-ply causes less deformation of rubber and consequently decreases the temperature.     

Thermal-fluid-structure coupling analysis for valve plate friction pair of axial piston pump in electrohydrostatic actuator (EHA) of aircraft

This paper deals with thermal-fluid-structure coupling analysis for valve plate friction pair of axial piston pump in electrohydrostatic actuator (EHA) of aircraft. The axial piston pump with high pressure and high rotational speed to be widely applied in EHA of more electric aircraft can increase the power density, but it also deteriorates thermal-fluid-structure coupling of the friction pairs. In order to reveal its interior multiphysics field coupling mechanism, taking the valve plate friction pair in three key friction pairs for example, this study carries out the research on multiphysics field coupling. Firstly, Navier–Stokes equations and energy equation of the incompressible fluid considering the influence of temperature and pressure on the oil properties, heat conduction governing equation with many boundary conditions including heat flux, heat convection, heat radiation and considering the influence of the structure deformation on the temperature and the influence of the temperature on the material properties, the elastic mechanics model of the structure exerted together by temperature, fluid pressure and mechanical load, are established. On this basis, a complete set of fast and effective thermal-fluid-structure coupling method is originally presented, and the numerical analysis is conducted using it for the valve plate friction pair. By the calculation results, the evolution laws with time and space are revealed regarding to the pressure and temperature of the fluid in the chambers, and the temperature, stress and deformation of the valve plate friction pair, the wedge-shaped clearance forms between them, even mixed friction occurs, and the corresponding improving measures aimed at the discovered problems are discussed. These results can provide the theoretical evidence for the design and development of the pump of EHA.     

Extensional and flexural modes of Rayleigh–Lamb wave in an orthotropic thermoelastic layer lying over a viscoelastic half-space

In this article the characteristics of the extensional and flexural modes, propagating in a thermoelastic orthotropic layer lying over a viscoelastic half-space, are analyzed. The complete analysis is carried out in the framework of a thermodynamically consistent hyperbolic type heat conduction model without energy dissipation. The normal-mode-analysis is adopted and a general form of dispersive equation is derived for an anisotropic thermoelastic layered medium. A prominent distinction with the isotropic elastic solids is observed in the symmetric as well as anti-symmetric modes of dispersion curves. In turn, such deformation reshapes the wave propagation while the deformation stiffening changes significantly the phase velocities of the wave till the acoustic radiation stresses are balanced by elastic stresses in the current configuration of the hyperelastic medium.     

Transient thermoelastic response in a cracked strip of functionally graded materials via generalized fractional heat conduction

This work is devoted to analyzing a thermal shock problem of an elastic strip made of functionally graded materials containing a crack parallel to the free surface based on a generalized fractional heat conduction theory. The embedded crack is assumed to be insulated. The Fourier transform and the Laplace transform are employed to solve a mixed initial-boundary value problem associated with a time-fractional partial differential equation. Temperature and thermal stresses in the Laplace transform domain are evaluated by solving a system of singular integral equations. Numerical results of the thermoelastic fields in the time domain are given by applying a numerical inversion of the Laplace transform. The temperature jump between the upper and lower crack faces and the thermal stress intensity factors at the crack tips are illustrated graphically, and phase lags of heat flux, fractional orders, and gradient index play different roles in controlling heat transfer process. A comparison of the temperature jump and thermal stress intensity factors between the non-Fourier model and the classical Fourier model is made. Numerical results show that wave-like behavior and memory effects are two significant features of the fractional Cattaneo heat conduction, which does not occur for the classical Fourier heat conduction.     

Optimal control of heat-pump/heat-storage systems with time-of-day energy price incentive

The maximum principle of optimal control theory is applied to the problem of optimizing the operation of a heat pump, when a storage capability is available and the electrical utility offers time-of-day price incentives in order to help level its diurnal load profile. The cost functional for optimal control is the monetary cost of purchased electrical energy. A bilinear model for the heat pump is assumed. When the ambient temperature is cyclic over the 24-hour period of the price pattern, periodic boundary conditions apply and the closed extremal trajectories are found to be unique and easily determined with a one-dimensional numerical search. These extremals have simple characteristics and reveal plausible strategies for minimizing the cost of purchased energy. They are potentially implementable with a simple, micro-processor-based controller.This work was supported in part by Energy, Mines, and Resources, Canada, Research Agreement No. 89.     

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.     

基于分数阶热传导方程激光加热瞬态温度场研究

基于分数阶Taylor(泰勒)级数展开原理,建立单相延迟一阶分数阶近似方程,获得分数阶热传导方程.针对短脉冲激光加热问题建立分数阶热传导方程组,并运用Laplace(拉普拉斯)变换方法进行求解,给出非Gauss(高斯)时间分布的激光内热源温度场解析解.针对具体算例数值研究温度波传播特性.结果表明热传播速度与分数阶阶次有关,分数阶阶次增加,热传播速度减小,温度变化幅度增加.分数阶方程可以用于描述介于扩散方程和热波方程间的热传输过程,且对热传播机制与分数阶热传导方程中分数阶项的关系做了深入剖析.     

Shape-preserving signal forms in heat conduction

Let us consider the heat conduction problem described by a parabolic equation. We study under which conditions is the time-dependence on the boundary preserved inside the solid. The question is how information entering on the boundary penetrates the solid. E.g. consider a heat conducting solid subject to sinusoidally varying boundary condition. After decay of the transients, the temperature at any inner point varies in time sinusoidally with the same circular frequency, with space dependent amplitude and phase delay. So, sinusoidal signals inserted on the boundary are preserved. Information is also preserved in case of linear signals. Farkas and Mudri [H. Farkas, I. Mudri, Shape-preserving time-dependences in heat conduction, Acta Phys. Hung. 55 (1984) 267–273] have formulated this phenomenon, defined the notion of the boundary following solution and the shape-preserving signal forms, determined necessary and heuristic sufficient conditions for the shape-preserving signal forms.

Their work is extended by rigorous proofs of some sufficient conditions in this paper, and the minimum of the phase delay, expected to be attained on the boundary for physical reasons, is examined.     

Similarity solution for a cylindrical shock wave in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux

The propagation of shock waves in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux, which has a variable azimuthally fluid velocity together with a variable axial fluid velocity, is investigated. The dusty gas is assumed to be a mixture of non-ideal (or perfect) gas and small solid particles, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-condition is maintained and variable energy input is continuously supplied by the piston (or inner expanding surface). The fluid velocities in the ambient medium are assume to be vary and obey power laws. The density of the ambient medium is assumed to be constant, the heat conduction is express in terms of Fourier’s law and the radiation is considered to be of the diffusion type for an optically thick grey gas model. The thermal conductivity K and the absorption coefficient α_R are assumed to vary with temperature and density. In order to obtain the similarity solutions the angular velocity of the ambient medium is assume to be decreasing as the distance from the axis increases. The effects of the variation of the heat transfer parameter and non-idealness of the gas in the mixture are investigated. The effects of an increase in (i) the mass concentration of solid particles in the mixture and (ii) the ratio of the density of solid particles to the initial density of the gas on the flow variables are also investigated.     

A finite difference scheme for two-dimensional semiconductor device of heat conduction on composite triangular grids

The momentary state of a semiconductor device of heat conduction is described by a system of four nonlinear partial differential equations. One elliptic equation is for the electrostatic, two parabolic equations are for the electron concentration and the hole concentration, and one heat exchange equation is for the temperature. According to the necessary of practical numerical simulations and based on the balance equation, finite difference schemes for two-dimensional transient behavior of a semiconductor device of heat conduction on composite triangular grids are constructed. Studying their stability and convergence properties, the error estimate in the energy norm is obtained. Finally, a numerical example is given.     

Coupled fluid-solid thermal interaction modeling for efficient transient simulation of biphasic water-steam energy systems

This paper proposes a fluid-solid coupled finite element formulation for the transient simulation of water-steam energy systems with phase change due to boiling and condensation. As it is commonly assumed in the study of thermal systems, the transient effects considered are exclusively originated by heat transfer processes. A homogeneous mixture model is adopted for the analysis of biphasic flow, resulting in a nonlinear transient advection-diffusion-reaction energy equation and an integral form for mass conservation in the fluid, coupled to the linear transient heat conduction equation for the solid. The conservation equations are approximated applying a stabilized Petrov-Galerkin FEM formulation, providing a set of coupled nonlinear equations for mass and energy conservation. This numerical model, combined with experimental heat transfer coefficients, provides a comprehensive simulation tool for the coupled analysis of boiling and condensation processes. For the treatment of enthalpy discontinuities traveling with the flow, a novel explicit-implicit time integration method based on Crank-Nicolson scheme is proposed, analyzing its accuracy and stability properties. To reduce problem size and enhance numerical efficiency, a modal superposition method with balanced truncation is applied to the solid equations. Finally, different example problems are solved to demonstrate the capabilities, flexibility and accuracy of the proposed formulation.     

Optimization of the design of recuperative heat exchangers in the exhaust nozzle of an aero engine

Today the needs for safer, cleaner and more affordable civil aero engines are found to be of great importance. Five years ago, the EU initiated an action for the design and the construction of efficient and environmentally friendly aero engines (EEFAE). One of the major European gas turbine industries, MTU, has presented a new technology for an advanced aero engine design, which uses an alternative thermodynamic cycle. The basis of this cycle is the adoption of a recuperation part with the use of a system of heat exchangers, installed in the exhaust nozzle of the aircraft engine. Thermal energy in the turbine exhaust is used in the recuperator to pre-heat the compressor outlet air before combustion. The benefits of this technique are focused on reduced pollutants and decreased fuel consumption. In this work, the procedure of the optimization of this installation, by means of the imposed pressure drop downstream the aircraft engine and the balanced mass inflow to the heat exchangers is presented. The optimization is based on experimental measurements in laboratory conditions and preliminary 2D CFD modeling for the flow inside the exhaust duct and through the heat exchangers. It is shown that with a careful approach, a better arrangement of the heat exchangers can be achieved in order to have a minimum pressure drop in the exhaust nozzle which can positively affect the engine’s performance.     

Estimation of dry road braking distance considering frictional energy of patterned tires

This paper is concerned with the braking distance estimation of tire controlled by anti-lock brake system (ABS) according to a numerical–analytical method. While the frictional heat dissipation at disc pad is derived analytically, the tire frictional energy loss is computed by the 3D dynamic rolling analysis of patterned tire. Since the tire rolling analysis to obtain the time history of the frictional energy rate for the entire braking period is impractical, we alternatively seek the tire frictional energy rate curve versus the lapse of time by interpolating the discretized frictional energy rates computed at intervals of 10 km/h. The effect of ABS is numerically implemented by specifying the corresponding tire angular velocity to the dynamic rolling analysis. Applying the energy conservation law to each speed interval determines the interval-wise braking times and distances from which the total braking time and distance are predicted. Illustrative numerical experiment is presented together with the comparison with the experimental estimation.     

Equal-Stress Reinforcement of Elastoplastic Momentless Shells with a Protective Coating under Thermal and Force Loadings

The problem of equal-stress reinforcement of three-layered momentless shells with fibers of constant cross section under thermoelastic and thermoplastic deformation of phase materials is formulated. A qualitative analysis of systems of resolving equations is carried out. Analytical solutions to the axisymmetric problems of equal-stress reinforcement of shells of revolution under elastic and inelastic deformation are constructed. It is shown that the problem can have several alternative solutions, which can be additionally controlled by varying the reinforcement intensities on the shell contour.     

Numerical modelling of EHD effects on heat transfer and bubble shapes of nucleate boiling

In this article, mathematical and numerical models are developed to study pure electrohydrodynamic (EHD) effects on heat transfer and bubble shapes when an initial bubble attached to a superheated horizontal wall in nucleate boiling. In the modelling of EHD effects on heat transfer, an undeformed bubble is considered; the electric body force and Joule heat are added to the momentum and energy equations; governing equations for heat, fluid flow and electric fields are coupled numerically and solved using a non-orthogonal body-fitted mesh system with necessary interfacial treatments at the gas–liquid boundary. While, to study the pure effect of EHD on the deformation of the bubble, the evaluation of a deformable bubble without heat transfer is simulated by volume of fluid (VOF) method based on an axial symmetric Cartesian coordinate system. The simulations indicate that EHD can effectively enhance heat transfer rate of nucleate boiling by influencing the motion of the ring vortex around the bubble and that bubble can be elongated due to the pull in axial direction and push in the negative radial direction by the electric field force.     

Estimation of temperature in rubber-like materials using non-linear finite element analysis based on strain history

A finite element procedure for hyper-elastic materials such as rubber has been developed to estimate the temperature rise during cyclic loading. The irreversible mechanical work developed in rubber has been used to determine the heat generation rate for carrying out thermal analysis. The evaluation of the heat energy is dependent on the strains. The finite element analysis assumes Green–Lagrangian strain displacement relations, Mooney–Rivlin strain energy density function for constitutive relationship, incremental equilibrium equations, and Total Lagrangian approach and the stress and strain of the rubber-like materials are evaluated using a degenerated shell element with assumed strain field technique, considering both material and geometric non-linearities. A transient heat conduction analysis has been carried out to estimate the temperature rise for different time steps in rubber-like materials using Galerkin's formulations. A numerical example is presented and the computed temperature values for various load steps agree closely with the experimental results reported in the literature.     

Spatial behaviour of solutions of the three-phase-lag heat equation

In this paper we study the spatial behaviour of solutions for the three-phase-lag heat equation on a semi-infinite cylinder. The theory of three-phase-lag heat conduction leads to a hyperbolic partial differential equation with a fourth-order derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary-value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial-time lines. A class of non-standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T₀ are assumed proportional to their initial values. These results are relevant in the context of the Saint-Venant Principle for heat conduction problems.