Magneto–thermo–electro–elastic transient response in a piezoelectric hollow cylinder subjected to complex loadings

The article presents an analytical solution for magneto–thermo–electro–elastic problems of a piezoelectric hollow cylinder placed in an axial magnetic field subjected to arbitrary thermal shock, mechanical load and transient electric excitation. Using an interpolation method solves the Volterra integral equation of the second kind caused by interaction among magnetic, thermal, electric and mechanical fields, the electric displacement is determined. Thus, the exact expressions for the transient responses of displacement, stresses, electric displacement, electric potential and perturbation of the magnetic field vector in the piezoelectric hollow cylinder are obtained by means of Hankel transforms, Laplace transforms, and inverse Laplace transforms. From sample numerical calculations, it is seen that the present method is suitable for a piezoelectric hollow cylinder subjected to arbitrary thermal shock, mechanical load and transient electric excitation, and the result carried out may be used as a reference to solve other transient coupled problems of magneto–thermo–electro–elasticity.     

A face-based monolithic approach for the incompressible magnetohydrodynamics equations

A novel numerical algorithm has been developed to solve the incompressible resistive magnetohydrodynamics equations in a fully coupled form. The numerical method is based on the face-centered unstructured finite volume approximation, where the velocity and magnetic field vector components are defined at the center of edges/faces; meanwhile, the pressure term is defined at element centroid. In order to enforce a divergence-free magnetic field, the gradient of a scalar Lagrange multiplier is introduced into the induction equation. A special attention will be given to satisfy the continuity equation and the Gauss' law for magnetism within each element and the summation of the equations can be exactly reduced to the domain boundary. The first modification to the original algorithm involves the evaluation of the convective fluxes over the two neighboring elements, where the discrete continuity equations are exactly satisfied. The second modification is based on the neglecting electric field term from the Lorentz force in two dimensions. The resulting large-scale algebraic linear equations are solved in a fully coupled manner using the one- and two-level restricted additive Schwarz preconditioners to avoid any time step restrictions forced by stability requirements. The spatial convergence of the algorithm is confirmed by solving the Hartmann flow, and then the algorithm is applied to the classical lid-driven cavity and backward facing step benchmark problems in two and three dimensions. The lid-driven cavity flow calculations at relatively high Stuart numbers indicate the perfect braking effect of the magnetic field in two dimensions.     

Thermomagnetic viscoelastic responses in a functionally graded hollow structure

This paper presents an analytical solution for the interaction of electric potentials,electric displacements,elastic deformations,and thermoelasticity,and describes electromagnetoelastic responses and perturbation of the magnetic field vector in hollow structures(cylinder or sphere),subjected to mechanical load and electric potential.The material properties,thermal expansion coefficient and magnetic permeability of the structure are assumed to be graded in the radial direction by a power law distribution.In the present model we consider the solution for the case of a hollow structure made of viscoelastic isotropic material,reinforced by elastic isotropic fibers,this material is considered as structurally anisotropic material.The exact solutions for stresses and perturbations of the magnetic field vector in FGM hollow structures are determined using the infinitesimal theory of magnetothermoelasticity,and then the hollow structure model with viscoelastic material is solved using the correspondence principle and Illyushin’s approximation method.Finally,numerical results are carried out and discussed.     

考虑电场梯度的挠曲电纳米板弯曲性能分析

具有尺度依赖的挠曲电效应在器件的设计中扮演着越来越关键的角色, 研究人员在微纳米尺度多物理场分析中进行了大量工作. 基于考虑挠曲电和电场梯度效应的弹性介电材料非经典理论, 以二维纳米板为例, 通过理论建模, 分析纳米板在弯曲问题中的力?电耦合行为. 根据Mindlin假设给出板的位移场和电势场的一阶截断, 选取板的材料为立方晶体(m3m点群), 将广义三维本构方程代入到高阶应力、高阶偶应力、高阶电位移和高阶电四极矩的表达式中得到相应的二维本构方程, 利用弹性电介质变分原理得到板的控制方程和边界上的线积分等式, 分别将二维本构方程和边界上外法线的方向余弦代入, 得到板的高阶弯曲方程、高阶电势方程以及对应的四边简支边界条件. 利用四边简支矩形板的高阶弯曲方程、高阶电势方程和相应的边界条件, 根据Navier解理论, 求解纳米板的电势场, 重点分析电场梯度对板内一阶电势的影响. 数值计算结果表明: 电场梯度对纳米板中由挠曲电效应产生的一阶电势有削弱作用, 且材料参数g₁₁越大, 一阶电势受到的削弱越大; 同时电场梯度的存在消除了纳米板在受横向集中载荷作用时一阶电势的奇异性. 本文是对具有挠曲电效应和电场梯度效应的纳米板结构分析理论的一个扩展, 为微纳米尺度器件的结构设计提供参考.      

Exact Electromagnetothermoelastic Solution for a Transversely Isotropic Piezoelectric Hollow Sphere Subjected to Arbitrary Thermal Shock

This paper presents analytical study for electromagnetothermoelastic transient behavior of a transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Exact solutions for the transient responses of stresses, perturbation of magnetic field vector, electric displacement and electric potential in the transversely isotropic piezoelectric hollow sphere are obtained by means of the Hankel transform, the Laplace transform and their inverse transforms. An interpolation method is used to solve the Volterra integral equation of the second kind caused by interactions among electric, magnetic, thermal and elastic fields. From the sample numerical calculations, it is seen that the present method is suitable for the transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Finally, the result can be used as a reference to solve other transient coupling problems of electromagnetothermoelasticity.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

Magnetic and Gravitational Convection of Air in a Porous Cubic Enclosure with a Coil Inclined Around the Y Axis

The natural convection heat transfer of air in a porous media can be controlled by gradient magnetic field. Thermomagnetic convection of air in a porous cubic enclosure with an electric coil inclined around the $Y$ axis was numerically investigated. The Biot–Savart law was used to calculate the magnetic field. The governing equations in primitive variables were discretized by the finite-volume method and solved by the SIMPLE algorithm. The flow and temperature fields for the air natural convection were presented and the mean Nusselt number on the hot wall was calculated and compared. The results show that both the magnetic force and coil inclination have significant effect on the flow field and heat transfer in a porous cubic enclosure, the natural convection heat transfer of air can be enhanced or controlled by applying gradient magnetic field.     

功能梯度热释电材料矩形板的三维精确分析

对四边简支、接地、等温的功能梯度热释电材料矩形板进行精确三维分析．根据正交各向异性热释电材料基本方程，导出了功能梯度热释电材料的状态方程．假定材料的机械、电学和热学性质沿板厚方向按统一的指数函数形式梯度分布，获得了四边简支、接地和等温的矩形板，在上下表面作用任意的机械荷载、电荷载和热荷载情况下的三维精确解．通过算例，分析了在机械荷载、电荷载和热荷载分别作用下，材料性质的不同梯度变化对平板结构响应的影响．所获得的精确结果可作为评价其他近似方法的标准解答或者作为建立简化的功能梯度热释电材料平板理论的基础。     

Unsteady hydromagnetic Couette flow of dusty fluid with temperature dependent viscosity and thermal conductivity

In the present paper the unsteady Couette flow and heat transfer of a dusty conducting fluid between two parallel plates with temperature dependent viscosity and thermal conductivity are studied. A constant pressure gradient and an external uniform magnetic field are applied. The governing coupled momentum and energy equations are solved numerically using finite differences. The effect of the variable viscosity and thermal conductivity of the fluid and the uniform magnetic field on the velocity and temperature fields for both the fluid and dust particles is discussed.     

Entropy analysis in electrical magnetohydrodynamic(MHD) flow of nanofluid with effects of thermal radiation,viscous dissipation,and chemical reaction

The unsteady mixed convection flow of electrical conducting nanofluid and heat transfer due to a permeable linear stretching sheet with the combined effects of an electric field, magnetic field, thermal radiation, viscous dissipation, and chemical reaction have been investigated. A similarity transformation is used to transform the constitutive equations into a system of nonlinear ordinary differential equations.The resultant system of equations is then solved numerically using implicit finite difference method.The velocity, temperature, concentration, entropy generation, and Bejan number are obtained with the dependence of different emerging parameters examined. It is noticed that the velocity is more sensible with high values of electric field and diminished with a magnetic field. The radiative heat transfer and viscous dissipation enhance the heat conduction in the system. Moreover, the impact of mixed convection parameter and Buoyancy ratio parameter on Bejan number profile has reverse effects. A chemical reaction reduced the nanoparticle concentration for higher values.     

Unsteady hydromagnetic channel flow of dusty fluid with temperature dependent viscosity and thermal conductivity

In the present paper the unsteady flow and heat transfer of a dusty conducting fluid between two parallel plates with temperature dependent viscosity and thermal conductivity are studied. A constant pressure gradient and an external uniform magnetic field is applied. The governing coupled momentum and energy equations are solved numerically using finite differences. The effect of the variable viscosity and thermal conductivity of the fluid and the uniform magnetic field on the velocity and temperature fields for both the fluid and dust particles is discussed.On leave from Department of Mathematics and Physics, Faculty of Engineering, El-Fayoum University, Egypt     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

Heating of a medium as a result of Joule energy dissipation

It is well known that when a magnetic field is present and electric currents flow through a gas, terms over and above those present in the case of ordinary thermal conductivity appear in the heat flux density vector. If the gas is dense enough and the magnetic field not over large, then the anisotropy caused by the magnetic field may be neglected. However, for a sufficiently large electric current a term proportional to the temperature and to the current density vector remains in the heat flux density vector. This effect explains, for example, the asymmetry of heat fluxes in the electrodes of a continuously operated electromagnetic accelerator (1).We shall consider this situation in relation to the example of a fully ionized quasineutral gas with identical electron and ion temperatures.     

Effects of heat generation on g-Jitter induced natural convection flow in a channel with isothermal or isoflux walls

This paper discusses the behavior of g-jitter induced free convection in microgravity under the influence of a transverse magnetic field and in the presence of heat generation or absorption effects for a simple system consisting of two parallel impermeable infinite plates held at four different thermal boundary conditions. The governing equations for this problem are derived on the basis of the balance laws of mass, linear momentum, and energy modified to include the effects of thermal buoyancy, magnetic field and heat generation or absorption as well as Maxwell's equations. The fluid is assumed to be viscous, Newtonian and have constant properties except the density in the body force of the balance of linear momentum equation. The governing equations are solved analytically for the induced velocity and temperature distributions as well as for the electric field and total current for electrically-conducting and insulating walls. This is done for isothermal–isothermal, isoflux–isothermal, isothermal–isoflux and isoflux–isoflux thermal boundary conditions. Graphical results for the velocity amplitude and distribution are presented and discussed for various parametric physical conditions.     

Analytical solutions for plane problem of functionally graded magnetoelectric cantilever beam

In this paper, an exact analytical solution is presented for a transversely isotropic functionally graded magneto-electro-elastic (FGMEE) cantilever beam, which is subjected to a uniform load on its upper surface, as well as the concentrated force and moment at the free end. This solution can be applied for any form of gradient distribution. For the basic equations of plane problem, all the partial differential equations governing the stress field, electric, and magnetic potentials are derived. Then, the expressions of Airy stress, electric, and magnetic potential functions are assumed as quadratic polynomials of the longitudinal coordinate. Based on all the boundary conditions, the exact expressions of the three functions can be determined. As numerical examples, the material parameters are set as exponential and linear distributions in the thickness direction. The effects of the material parameters on the mechanical, electric, and magnetic fields of the cantilever beam are analyzed in detail.     

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

Spontaneous rotation in the exact solution of magnetohydrodynamic equations for flow between two stationary impermeable disks

Magnetohydrodynamic (MHD) flow of a viscous electrically conducting incompressible fluid between two stationary impermeable disks is considered. A homogeneous electric current density vector normal to the surface is specified on the upper disk, and the lower disk is nonconducting. The exact von Karman solution of the complete system of MHD equations is studied in which the axial velocity and the magnetic field depend only on the axial coordinate. The problem contains two dimensionless parameters: the electric current density on the upper plate Y and the Batchelor number (magnetic Prandtl number). It is assumed that there is no external source that produces an axial magnetic field. The problem is solved for a Batchelor number of 0–2. Fluid flow is caused by the electric current. It is shown that for small values of Y, the fluid velocity vector has only axial and radial components. The velocity of motion increases with increasing Y, and at a critical value of Y, there is a bifurcation of the new steady flow regime with fluid rotation, while the flow without rotation becomes unstable. A feature of the obtained new exact solution is the absence of an axial magnetic field necessary for the occurrence of an azimuthal component of the ponderomotive force, as is the case in the MHD dynamo. A new mechanism for the bifurcation of rotation in MHD flow is found.     

Steady MHD Couette flow with temperature-dependent physical properties

The influence of variation in physical variables on the steady magnetohydrodynamic (MHD) Couette flow with heat transfer is studied. An external uniform magnetic field is applied perpendicular to the parallel plates and the fluid is acted upon by a constant pressure gradient. The viscosity and the thermal as well as electric conductivities are assumed to be temperature dependent. The two plates are kept at two constant but different temperatures, and the viscous and Joule dissipations are considered in the energy equation. A numerical solution for the governing nonlinear coupled equations of motion and the energy equation is obtained. The effect of the temperature-dependent viscosity, thermal conductivity, and electrical conductivity on both the velocity and temperature distributions is examined. H.A. Attia - On leave from: Dept. of Eng. Mathematics and physics, El-Fayoum University, El-Fayoum, Egypt     

Electromagnetohydrodynamic (EMHD) flow of fractional viscoelastic fluids in a microchannel

This study investigates the electromagnetohydrodynamic (EMHD) flow of fractional viscoelastic fluids through a microchannel under the Navier slip boundary condition. The flow is driven by the pressure gradient and electromagnetic force where the electric field is applied horizontally, and the magnetic field is vertically (upward or downward). When the electric field direction is consistent with the pressure gradient direction, the changes of the steady flow rate and velocity with the Hartmann number Ha are irrelevant to the direction of the magnetic field (upward or downward). The steady flow rate decreases monotonically to zero with the increase in Ha. In contrast, when the direction of the electric field differs from the pressure gradient direction, the flow behavior depends on the direction of the magnetic field, i.e., symmetry breaking occurs. Specifically, when the magnetic field is vertically upward, the steady flow rate increases first and then decreases with Ha. When the magnetic field is reversed, the steady flow rate first reduces to zero as Ha increases from zero. As Ha continues to increase, the steady flow rate (velocity) increases in the opposite direction and then decreases, and finally drops to zero for larger Ha. The increase in the fractional calculus parameter α or Deborah number De makes it take longer for the flow rate (velocity) to reach the steady state. In addition, the increase in the strength of the magnetic field or electric field, or in the pressure gradient tends to accelerate the slip velocity at the walls. On the other hand, the increase in the thickness of the electric double-layer tends to reduce it.

     

无限长圆柱体中广义电磁热弹耦合问题研究

基于L-S广义热弹性理论,研究了处于磁场中无限长理想圆柱导体在边界受热冲击作用时的电磁热弹耦合问题.建立了广义电磁热弹耦合的有限元方程,为避免积分变换方法求解带来的精度丟失.采用将有限元方程直接在时间域求解的方法,得到了圆柱体中的温度、位移、应力、感应磁场和感应电场的分布规律,反映了热的波动性及电磁热弹的耦合效应.结果表明,将有限元方程直接在时间域求解,可以获得各物理量的准确分布.得到温度在热波波前处的阶跃,准确地反应热波的波动效应.