Modeling of electrodynamic-mechanical coupling phenomena in smart magnetoelectroelastic materials

The coupling of electric, magnetic and mechanical phenomena may have various reasons. The famous Maxwell equations of electrodynamics describe the interaction of transient magnetic and electric fields. On the constitutive level of dielectric materials, coupling mechanisms are manyfold comprising piezoelectric, magnetostrictive or magnetoelectric effects. Electromagnetically induced specific forces acting at the boundary and within the domain of a dielectric body are, within a continuum mechanics framework, commonly denoted as Maxwell stresses. In transient electromagnetic fields, the Poynting vector gives another contribution to mechanical stresses. First, a system of transient partial differential equations is presented. Introducing scalar and vector potentials for the electromagnetic fields and representing the mechanical strain by displacement fields, seven coupled differential equations govern the boundary value problem, accounting for linear constitutive equations of magnetoelectroelasticity. To reduce the effort of numerical solution, the system of equations is partly decoupled applying generalized forms of Coulomb and Lorenz gauge transformations [1,2]. A weak formulation is given to establish a basis for a finite element solution. The influence of constitutive magnetoelectric coupling on electromagnetic wave propagation is finally demonstrated with a simple one-dimensional example. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

离散外微分与计算电磁学

计算电磁学的核心之一是数值求解Maxwell方程组.适当的离散方式是保证结果能真实反映物理现象的关键.为了在离散的过程中保持该方程组的几何性质,我们建立了基于棱柱网格的系数为R的格点规范理论,其离散曲率满足相应的Bianchi恒等式.通过适当定义离散微分形式之间的内积和棱柱网格上的Hodge星算子,我们由离散变分导出源方程和连续性方程,和Bianchi恒等式一起称为真空中的离散Maxwell方程组.这组方程是内蕴的,并具有规范不变性.     

The integral equation method in problems of diffraction by a finite impedance section of a medium interface

A 3D problem of reflection of a plane electromagnetic wave by a local impedance section of a wavy surface is considered. The boundary value problem for the system of Maxwell’s equations in a region with an irregular boundary is reduced to solution of systems of hypersingular integral equations. A numerical algorithm is proposed for solution of these systems. Results of numerical computations are presented.     

Solitary waves in the nolinear wave equation and in gauge theories

Roughly speaking a solitary wave is a solution of a field equation whose energy travels as a localized packet and which preserves this localization in time. This paper is an introduction to the study of solitary waves relative to the nonlinear wave equation and to the Abelian gauge theories. Abelian gauge theories consist of a class of field equations obtained by coupling in a suitable way the nonlinear wave equation with the Maxwell equations. They provide a model for the interaction of matter with the electromagnetic field. One of the motivations of this study lies in the fact that the nonlinear wave equation and the Abelian gauge theories are the simplest equations which satisfy the basic principles of modern physics. Dedicated to the memory of Jean Leray     

Method of Integral Equations for the Three-Dimensional Problem of Wave Reflection from an Irregular Surface

The problem on the reflection of the field of a plane H-polarized three-dimensional electromagnetic wave from a perfectly conducting interface between media which contains a local perfectly conducting inhomogeneity is considered. To construct a numerical algorithm, the boundary value problem for the system of Maxwell equations in an infinite domain with irregular boundary is reduced to a system of singular integral equations, which is solved by the approximation–collocation method. The elements of the resulting complex matrix are calculated by a specially developed algorithm. The solution of the system of singular integral equations is used to obtain an integral representation for the reflected electromagnetic field and computational formulas for the directional diagram of the reflected electromagnetic field in the far region.     

Light scattering by a homogeneous sphere near a plane boundary

The scattering of electromagnetic waves by a homogeneous sphere near a plane boundary is presented in this paper. The vector wave equations derived from Maxwell’s equations are solved by means of the two orthogonal solutions to the scalar wave equation. Hankel transformation and Erdélyi’s formula are used to satisfy the planar boundary conditions and the determination of the unknown coefficients in the scattered field and internal fields is achieved by matching the electromagnetic boundary conditions on the surface of the sphere. Existence and uniqueness of the solution of the series involving these unknown coefficients are shown.     

Scattering of the electromagnetic field on a finite impedance section of an interface

The scattering of a three-dimensional electromagnetic field on a local impedance section of a wavy surface is considered. The boundary-value problem for the system of Maxwell equations is reduced to solving a system of hypersingular integral equations. A numerical algorithm is developed and a program is designed for computing the electrodynamic characteristics of the given scattering problem. The accuracy of the simulation results is investigated. __________ Translated from Prikladnaya Matematika i Informatika, No. 25, pp. 56–69, 2007.     

Numerical Model for Electro-Mechanical Analysis Based on Finite Element Method

In this paper, we present a numerical procedure that can be used to model the electro-mechanical coupled behavior of the dielectric actuator domain. The equation describing the electrostatical part is given by the reduced form of the Maxwell equation and the electrostatic potential [1]. The mechanical problem is described by the constitutive equations and equilibrium equations. Using the finite element method, this technique is to divide a whole problem into sub-problems. The complexity of the original problem is therefore reduced by focusing only on the most relevant areas. A finite element analysis is then performed by applying the electrostatic Maxwell pressure as Neumann boundary conditions to compute the displacements. Once the displacement is computed, the electrostatic domain or the conductor is updated. Electrostatic analysis is performed on the updated geometry and the finite element method is then used to determine the change in potential due to geometric perturbations. Once the surface charge densities are known, the new electrostatic Maxwell pressure is computed. The mechanical and electrostatic analysis is repeated until an equilibrium state is computed. The procedure is demonstrated in the paper by the solution of some two-dimensional and three-dimensional problems. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Formalism of two potentials for the numerical solution of Maxwell’s equations

A new formulation of Maxwell’s equations based on the introduction of two vector and two scalar potentials is proposed. As a result, the electromagnetic field equations are written as a hyperbolic system that contains, in contrast to the original Maxwell system, only evolution equations and does not involve equations in the form of differential constraints. This makes the new equations especially convenient for the numerical simulation of electromagnetic processes. Specifically, they can be solved by applying powerful modern shock-capturing methods based on the approximation of spatial derivatives by upwind differences. The cases of an electromagnetic field in a vacuum and an inhomogeneous material are considered. Examples are given in which electromagnetic wave propagation is simulated by solving the formulated system of equations with the help of modern high-order accurate schemes.     

On a model of electromagnetic field propagation in ferroelectric media

The Maxwell system in an anisotropic, inhomogeneous medium with non-linear memory effect produced by a Maxwell type system for the polarization is investigated under low regularity assumptions on data and domain. The particular form of memory in the system is motivated by a model for electromagnetic wave propagation in ferromagnetic materials suggested by Greenberg, MacCamy and Coffman [J.M. Greenberg, R.C. MacCamy, C.V. Coffman, On the long-time behavior of ferroelectric systems, Phys. D 134 (1999) 362-383]. To avoid unnecessary regularity requirements the problem is approached as a system of space-time operator equation in the framework of extrapolation spaces (Sobolev lattices), a theoretical framework developed in [R. Picard, Evolution equations as space-time operator equations, Math. Anal. Appl. 173 (2) (1993) 436-458; R. Picard, Evolution equations as operator equations in lattices of Hilbert spaces, Glasnik Mat. 35 (2000) 111-136]. A solution theory for a large class of ferromagnetic materials confined to an arbitrary open set (with suitably generalized boundary conditions) is obtained.     

Maxwell’s equations in inhomogeneous bi-anisotropic materials: Existence,uniqueness and stability for the initial value problem

In the present paper inhomogeneous bi-anisotropic materials characterized by matrices of electric permittivity, magnetic permeability and magnetoelectric characteristics are considered. All elements of these matrices are functions of the position in three dimensional space. The time-dependent Maxwell’s equations describe the electromagnetic wave propagation in these materials. Maxwell’s equations together with zero initial data are analyzed and a statement of the initial value problem (IVP) is formulated. This IVP is reduced to the IVP for a symmetric hyperbolic system of partial differential equations of the first order. Applying the theory of a symmetric hyperbolic system, new existence, uniqueness and stability estimate theorems have been obtained for the IVP of Maxwell’s equations in inhomogeneous bi-anisotropic materials.     

Problem of determining the permittivity in the stationary system of Maxwell equations

The stationary system of Maxwell equations for a unmagnetized nonconducting medium is considered. For this system, the problem of determining the permittivity ε from given electric or magnetic fields is studied. It is assumed that the electromagnetic field is induced by a plane wave coming from infinity in the direction ν. It is also assumed that the permittivity is different from a given positive constant ε₀ only inside a compact domain Ω ? R ³ with a smooth boundary S. To find ε inside Ω, the solution of the corresponding direct problem for the system of electrodynamic equations on the shadow portion of the boundary of Ω is specified for all frequencies starting at some fixed ω₀ and for all ν. The high-frequency asymptotics of the solution to the direct problem is studied. It is shown that the information specified makes it possible to reduce the original problem to the well-known inverse kinematic problem of determining the refraction coefficient inside Ω from the traveling times of an electromagnetic wave. This leads to a uniqueness theorem for the solution of the problem under consideration and opens up the opportunity of its constructive solution.     

The Well-Posedness of Dynamical Equations of Magneto-electro-elasticity

A mathematical model of wave propagation in magneto-electro-elastic materials is obtained in the form of a symmetric hyperbolic system of the first-order partial differential equations. This model is a result of the qualitative analysis of the coupled time-dependent Maxwell’s equations and equations of elastodynamics which are considered together with constitutive relations in non-homogeneous anisotropic magneto-electro-elastic materials. Applying the theory and methods of symmetric hyperbolic systems, we have proved that the reported model of wave propagation in magneto-electro-elastic materials satisfies the Hadamards correctness requirements: solvability, uniqueness and stability with respect to perturbation of data.     

On soliton and periodic solutions of Maxwell–Bloch system for two-level medium with degeneracy

Maxwell–Bloch (MB) system describing the resonant propagation of electromagnetic pulses in either two-level media with degeneracy in angle moment projection or a three-level media with equal oscillator forces is considered. An inhomogeneous broadening of energy levels and a polarization of the wave are accounted. The equations are integrated by the binary Darboux transformations technique. Pulses corresponding to a transition between levels with the largest population difference are shown to be stable. The solution describing the propagation of pulses in the medium excited by a periodic wave is obtained. The hierarchy of infinitesimal symmetries is obtained by means of Darboux transformation.     

On the Fredholm property of the electric field equation in the vector diffraction problem for a partially screened solid

We consider a vector problem of diffraction of an electromagnetic wave on a partially screened anisotropic inhomogeneous dielectric body. The boundary conditions and the matching conditions are posed on the boundary of the inhomogeneity domain, and under passage through it, the medium parameters have jump changes. A boundary value problem for the system of Maxwell equations in unbounded space is studied in a semiclassical statement and is reduced to a system of integro-differential equations on the body domain and the screen surfaces. We show that the quadratic form of the problem operator is coercive and the operator itself is Fredholm with zero index.     

Flow and heat transfer of double fractional Maxwell fluids over a stretching sheet with variable thickness

This paper presents research on the fractional boundary layer flow and heat transfer over a stretching sheet with variable thickness. Based on the Caputo operators, the double fractional Maxwell model and generalized Fourier's law are introduced to the constitutive relationships. The governing equations are solved numerically by utilizing the finite difference method. The effects of fractional parameters on the velocity and temperature field are analyzed. The results indicate that the larger is the fractional stress parameter, the stronger is the elastic characteristic. However, fluids show viscous fluid-like behavior for a larger value of fractional strain parameter. Moreover, the numerical solutions are in good agreement with the exact solution and the convergence order can achieve the expected first order. The numerical method in this study is reliable and can be extended to other fractional boundary layer problems over a variable thickness sheet.     

电动力学电磁场边值问题的广义变分原理

给出了线性各项异性电磁场边值问题的广义虚功原理表达式,运用钱伟长教授提出的方法建立了该问题的广义变分原理,可直接反映该问题的全部特征,即4个Maxwell方程、2个场强-位势方程、2个本构方程和8个边界条件．继而导出了一族有先决条件的广义变分原理．作为例证,导出了两个退化形式的广义变分原理,和已知的广义变分原理等价．此外还导出了两个修正的广义变分原理,可为该问题提供杂交有限元模型．建立的各广义变分原理可为电磁场边值问题的有限元应用提供更为完善的理论基础．     

Symmetry properties of two-dimensional Ciarlet–Mooney–Rivlin constitutive models in nonlinear elastodynamics

Nonlinear dynamic equations for isotropic homogeneous hyperelastic materials are considered in the Lagrangian formulation. An explicit criterion of existence of a natural state for a given constitutive law is presented, and is used to derive natural state conditions for some common constitutive relations.For two-dimensional planar motions of Ciarlet–Mooney–Rivlin solids, equivalence transformations are computed that lead to a reduction of the number parameters in the constitutive law. Point symmetries are classified in a general dynamical setting and in traveling wave coordinates. A special value of traveling wave speed is found for which the nonlinear Ciarlet–Mooney–Rivlin equations admit an additional infinite set of point symmetries. A family of essentially two-dimensional traveling wave solutions is derived for that case.     

Einstein Equations on a 5-Manifold with a Causal Structure in the Absence of Matter Fields

A 5-manifold with a restricted smooth structure and an appropriate group of coordinate transformations including general relativity and gauge transformations is considered. An explicit expression for the Riemannian curvature of a 4-dimensional vector distribution is obtained, which implies the classical Einstein and Maxwell equations. Bibliography: 14 titles.