Analytical solutions to general anti-plane shear problems in finite elasticity

This paper presents a pure complementary energy variational method for solving a general anti-plane shear problem in finite elasticity. Based on the canonical duality–triality theory developed by the author, the nonlinear/nonconvex partial differential equations for the large deformation problem are converted into an algebraic equation in dual space, which can, in principle, be solved to obtain a complete set of stress solutions. Therefore, a general analytical solution form of the deformation is obtained subjected to a compatibility condition. Applications are illustrated by examples with both convex and nonconvex stored strain energies governed by quadratic-exponential and power-law material models, respectively. Results show that the nonconvex variational problem could have multiple solutions at each material point, the complementary gap function and the triality theory can be used to identify both global and local extremal solutions, while the popular convexity conditions (including rank-one condition) provide mainly local minimal criteria and the Legendre–Hadamard condition (i.e., the so-called strong ellipticity condition) does not guarantee uniqueness of solutions. This paper demonstrates again that the pure complementary energy principle and the triality theory play important roles in finite deformation theory and nonconvex analysis.     

Traveling Waves in Spatial SIRS Models

We study traveling wavefront solutions for two reaction–diffusion systems, which are derived respectively as diffusion approximations to two nonlocal spatial SIRS models. These solutions characterize the propagating progress and speed of the spatial spread of underlying epidemic waves. For the first diffusion system, we find a lower bound for wave speeds and prove that the traveling waves exist for all speeds bigger than this bound. For the second diffusion system, we find the minimal wave speed and show that the traveling waves exist for all speeds bigger than or equal to the minimal speed. We further prove the uniqueness (up to translation) of these solutions for sufficiently large wave speeds. The existence of these solutions are proved by a shooting argument combining with LaSalle’s invariance principle, and their uniqueness by a geometric singular perturbation argument.     

板中热弹波传播: 一种改进的勒让德多项式方法

近年来, 超声导波因其衰减小, 传播距离远和信号覆盖范围广, 成为无损检测领域快速发展的方向之一. 然而, 基于超声导波的高温在线检测和激光超声技术却发展缓慢, 其关键在于热弹耦合波动方程求解难度大、传播与衰减特性研究困难. 作为一种有效的求解方法, 勒让德正交多项式方法已广泛应用于导波传播问题, 但该方法在求解热弹导波传播时存在两个不足, 限制其进一步的发展和应用. 这两个缺陷是: (1)求解过程中大量积分的存在, 致使计算效率低下; (2)仅能处理等热边界条件的热弹导波传播. 针对两项不足之处, 提出一种改进的勒让德正交多项式方法, 以求解分数阶热弹板中的导波传播. 推导求解方法中积分的解析表达式, 以提高计算效率; 引入温度梯度展开式, 发展适合勒让德多项式级数的绝热边界条件处理方法. 与已有文献结果对比表明改进方法的正确性; 与已有方法的计算时间对比说明改进方法的高效性. 最后将改进的方法用于求解分数阶热弹板中的导波传播, 研究分数阶次对频散、衰减曲线和应力、位移、温度分布等的影响.      

THERMOELASTIC WAVE PROPAGATION IN PLATES: AN IMPROVED LEGENDRE POLYNOMIAL APPROACH 1)

近年来, 超声导波因其衰减小, 传播距离远和信号覆盖范围广, 成为无损检测领域快速发展的方向之一. 然而, 基于超声导波的高温在线检测和激光超声技术却发展缓慢, 其关键在于热弹耦合波动方程求解难度大、传播与衰减特性研究困难. 作为一种有效的求解方法, 勒让德正交多项式方法已广泛应用于导波传播问题, 但该方法在求解热弹导波传播时存在两个不足, 限制其进一步的发展和应用. 这两个缺陷是: (1)求解过程中大量积分的存在, 致使计算效率低下; (2)仅能处理等热边界条件的热弹导波传播. 针对两项不足之处, 提出一种改进的勒让德正交多项式方法, 以求解分数阶热弹板中的导波传播. 推导求解方法中积分的解析表达式, 以提高计算效率; 引入温度梯度展开式, 发展适合勒让德多项式级数的绝热边界条件处理方法. 与已有文献结果对比表明改进方法的正确性; 与已有方法的计算时间对比说明改进方法的高效性. 最后将改进的方法用于求解分数阶热弹板中的导波传播, 研究分数阶次对频散、衰减曲线和应力、位移、温度分布等的影响.     

General stability of memory-type thermoelastic Timoshenko beam acting on shear force

In this paper, we consider a linear thermoelastic Timoshenko system with memory effects where the thermoelastic coupling is acting on shear force under Neumann–Dirichlet–Dirichlet boundary conditions. The same system with fully Dirichlet boundary conditions was considered by Messaoudi and Fareh (Nonlinear Anal TMA 74(18):6895–6906, 2011, Acta Math Sci 33(1):23–40, 2013), but they obtained a general stability result which depends on the speeds of wave propagation. In our case, we obtained a general stability result irrespective of the wave speeds of the system.     

Acoustic scattering by a sphere in the time domain

A sound pulse is scattered by a sphere leading to an initial–boundary value problem for the wave equation. A method for solving this problem is developed using integral representations involving Legendre polynomials in a similarity variable and Volterra integral equations. The method is compared and contrasted with the classical method, which uses Laplace transforms in time combined with separation of variables in spherical polar coordinates.     

On Stability Analyses of Three Classical Buckling Problems for the Elastic Strut

It is common practice in analyses of the configurations of an elastica to use Jacobi??s necessary condition to establish conditions for stability. Analyses of this type date to Born??s seminal work on the elastica in 1906 and continue to the present day. Legendre developed a treatment of the second variation which predates Jacobi??s. The purpose of this paper is to explore Legendre??s treatment with the aid of three classical buckling problems for elastic struts. Central to this treatment is the issue of existence of solutions to a Riccati differential equation. We present two different variational formulations for the buckling problems, both of which lead to the same Riccati equation, and we demonstrate that the conclusions from Legendre and Jacobi??s treatments are equivalent for some sets of boundary conditions. In addition, the failure of both treatments to classify stable configurations of a free-free strut are contrasted.     

Kinetic Relations for a Lattice Model of Phase Transitions

The aim of this article is to analyse travelling waves for a lattice model of phase transitions, specifically the Fermi–Pasta–Ulam chain with piecewise quadratic interaction potential. First, for fixed, sufficiently large subsonic wave speeds, we rigorously prove the existence of a family of travelling wave solutions. Second, it is shown that this family of solutions gives rise to a kinetic relation which depends on the jump in the oscillatory energy in the solution tails. Third, our constructive approach provides a very good approximate travelling wave solution.     

A numerical investigation of a spectral‐type nodal collocation discontinuous Galerkin approximation of the Euler and Navier–Stokes equations

In this paper, we focus on the applicability of spectral‐type collocation discontinuous Galerkin methods to the steady state numerical solution of the inviscid and viscous Navier–Stokes equations on meshes consisting of curved quadrilateral elements. The solution is approximated with piecewise Lagrange polynomials based on both Legendre–Gauss and Legendre–Gauss–Lobatto interpolation nodes. For the sake of computational efficiency, the interpolation nodes can be used also as quadrature points. In this case, however, the effect of the nonlinearities in the equations and/or curved elements leads to aliasing and/or commutation errors that may result in inaccurate or unstable computations. By a thorough numerical testing on a set of well known test cases available in the literature, it is here shown that the two sets of nodes behave very differently, with a clear advantage of the Legendre–Gauss nodes, which always displayed an accurate and robust behaviour in all the test cases considered.Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial Dynamics of a Discrete-Time Population Model in a Periodic Lattice Habitat

This paper is devoted to the study of spatial dynamics of a class of discrete-time population models in a periodic lattice habitat. In the general case of recruitment functions, we obtain the existence and computation formula of spreading speeds and show that they coincide with the minimal wave speeds for periodic traveling waves in the positive and negative directions.     

Investigation of guided waves propagation in orthotropic viscoelastic carbon–epoxy plate by Legendre polynomial method

The effect of viscoelasticity on the guided waves propagation in viscoelastic plate has been investigated according to multi-aspect. To this purpose, an extension of the Legendre polynomial method is proposed to formulate the guided waves equation in orthotropic viscoelastic plate composed of carbon–epoxy. The validity of the proposed Legendre polynomial method is illustrated by comparison with available data. The convergence of the method is discussed through a numerical example. The hysteretic and Kelvin–Voigt viscoelastic models are used to integrate the imaginary part of the complex stiffness matrix associated with the viscoelastic plate in this study. Accordingly, both viscoelastic models do not affect on the dispersion curves results. However, appreciable effects are seen in the attenuation curves. Also, the sensitivity of the guided waves propagation caused by variations of elastic and viscoelastic modulus has been studied in detail. Finally, the advantages of the Legendre polynomial method are described.     

Ill-Posedness of the Hydrostatic Euler and Singular Vlasov Equations

In this paper, we develop an abstract framework to establish ill-posedness, in the sense of Hadamard, for some nonlocal PDEs displaying unbounded unstable spectra. We apply this to prove the ill-posedness for the hydrostatic Euler equations as well as for the kinetic incompressible Euler equations and the Vlasov–Dirac–Benney system.     

Wave speeds in the macroscopic relativistic extended model with many moments

An exact macroscopic extended model for relativistic gases, with an arbitrary number of moments, is present in the literature. It is determined except for a numberable family $\mathcal {F}$ of single variable functions whose physical meaning remains an open problem; a possibility is that it allows to apply the theory to a wider set of materials. Other models appearing in literature are particular cases of this macroscopic one when all these arbitrary functions are zero. Here we exploit equations determining wave speeds for that model. We find interesting results; for example, the whole system for their determination can be divided into independent subsystems which are expressed by linear combinations, through scalar coefficients, of tensors all of the same order. As expected, these wave speeds for the macroscopic model depend on $\mathcal {F}$ . Moreover, some wave speeds (but not all of them) are expressed by square roots of rational numbers.     

Travelling Waves for Complete Discretizations of Reaction Diffusion Systems

In this paper we consider the impact that full spatial–temporal discretizations of reaction–diffusion systems have on the existence and uniqueness of travelling waves. In particular, we consider a standard second-difference spatial discretization of the Laplacian together with the six numerically stable backward differentiation formula methods for the temporal discretization. For small temporal time-steps and a fixed spatial grid-size, we establish some useful Fredholm properties for the operator that arises after linearizing the system around a travelling wave. In particular, we perform a singular perturbation argument to lift these properties from the natural limiting operator. This limiting operator is associated to a lattice differential equation, where space has been discretized but time remains continuous. For the backward-Euler temporal discretization, we also obtain travelling waves for arbitrary time-steps. In addition, we show that in the anti-continuum limit, in which the temporal time-step and the spatial grid-size are both very large, wave speeds are no longer unique. This is in contrast to the situation for the original continuous system and its spatial semi-discretization. This non-uniqueness is also explored numerically and discussed extensively away from the anti-continuum limit.     

Nonlinear waves in electromigration dispersion in a capillary

We construct exact solutions to an unusual nonlinear advection–diffusion equation arising in the study of Taylor–Aris (also known as shear) dispersion due to electroosmotic flow during electromigration in a capillary. An exact reduction to a Darboux equation is found under a traveling-wave ansatz. The equilibria of this ordinary differential equation are analyzed, showing that their stability is determined solely by the (dimensionless) wave speed without regard to any (dimensionless) physical parameters. Integral curves, connecting the appropriate equilibria of the Darboux equation that governs traveling waves, are constructed, which in turn are shown to be asymmetric kink solutions (i.e., non-Taylor shocks). Furthermore, it is shown that the governing Darboux equation exhibits bistability, which leads to two coexisting non-negative kink solutions for (dimensionless) wave speeds greater than unity. Finally, we give some remarks on other types of traveling-wave solutions and a discussion of some approximations of the governing partial differential equation of electromigration dispersion.     

Double-grid Chebyshev spectral elements for acoustic wave modeling

Highly accurate algorithms are needed for modeling wave propagation phenomena in realistic media. The spectral element methods, either based on a Chebyshev or a Legendre polynomial basis, have shown their excellent properties of high accuracy and flexibility in describing complex models outperforming other techniques. In contrast with standard grid methods, which use dense spatial meshes, spectral element methods discretize the computational domain in a very coarse mesh. With constant-property elements, this fact may in some cases reduce seriously the computational efficiency. For instance, if the medium is finely heterogeneous, it may need to be described in a much finer way than the acoustic wave field. The double-grid approach presented in this work is a viable way for overcoming this lack of the method and for handling problems where the medium changes continuously or even sharply on the small scale. The variation in the properties is taken into account by using an independent set of shape functions defined on a temporary local grid in such a way that either the small scale fluctuations are accurately handled, without the need of a global finer grid, and the macroscopic wave field propagation is solved with no loose of computational efficiency.     

On instabilities of single-integral constitutive equations for viscoelastic liquids

Various types of instabilities are exposed in this paper for time-strain separable single-integral viscoelastic constitutive equations (CE's). They were distinguished into two groups and defined as Hadamard and dissipative type of instabilities. As for the Hadamard-type, previously obtained criteria are found to be necessary only. They are necessary and sufficient only for thermodynamic stability. Improved, stricter Hadamard stability criteria are described briefly in this paper, and then applied to study of stability of several CE's. It is shown that the Currie potential with the K-BKZ equation and the model proposed by Papanastasiou et al. are Hadamard unstable. In the case of dissipative stability, the necessary and sufficient condition for stress boundedness in any regular flow with a given history, is proved. Then, this criterion was applied to the neoHookean, Mooney, and Yen and McIntire specifications of the general K-BKZ model, to exhibit unbounded solutions. In addition, Larson-Monroe potential which is later proved to be Hadamard unstable but satisfies the above criterion of boundedness, is shown to have unstable decreasing branch in steady simple shear flow. At present, to the authors' knowledge, there is no viscoelastic single-integral CE of factorable type proposed in the literature which can satisfy all the Hadamard and dissipative stability criteria.     

On well-posedness of the dynamic problem for an anisotropic fluid-saturated solid

It is known that a high degree of anisotropy in the constitutive behaviour of a solid may result in the loss of hyperbolicity of the dynamic equations in the form of either complex-conjugate or purely imaginary characteristic wave speeds (flutter ill-posedness and shear band formation, respectively). In the present paper we investigate the characteristic wave speeds in the dynamic problem for a transversely isotropic fluid-saturated porous solid. Three cases are considered: a dry solid and a saturated solid under locally undrained and drained conditions. It is shown that, for given constitutive parameters of the solid skeleton, the dynamic problem for a drained solid may become ill-posed due to the flutter-type loss of hyperbolicity, while the dynamic equations for a dry and an undrained solids remain hyperbolic. For a given solid skeleton, the characteristic wave speeds are strongly influenced by the pore fluid compressibility which, in turn, is extremely sensitive to the presence of a small amount of free gas.     

Elastic response of water-filled fiber composite tubes under shock wave loading

We experimentally and numerically investigate the response of fluid-filled filament-wound composite tubes subjected to axial shock wave loading in water. Our study focuses on the fluid–structure interaction occurring when the shock wave in the fluid propagates parallel to the axis of the tube, creating pressure waves in the fluid coupled to flexural waves in the shell. The in-house-developed computational scheme couples an Eulerian fluid solver with a Lagrangian shell solver, which includes a new and simple material model to capture the response of fiber composites in finite kinematics. In the experiments and simulations we examine tubes with fiber winding angles equal to 45° and 60°, and we measure the precursor and primary wave speeds, hoop and longitudinal strains, and pressure. The experimental and computational results are in agreement, showing the validity of the computational scheme in complex fluid–structure interaction problems involving fiber composite materials subjected to shock waves. The analyses of the measured quantities show the strong coupling of axial and hoop deformations and the significant effect of fiber winding angle on the composite tube response, which differs substantially from that of a metal tube in the same configuration.     

A Whitham–Boussinesq long-wave model for variable topography

We study the propagation of water waves in a channel of variable depth using the long-wave asymptotic regime. We use the Hamiltonian formulation of the problem in which the non-local Dirichlet–Neumann operator appears explicitly in the Hamiltonian, and propose a Hamiltonian model for bidirectional wave propagation in shallow water that involves pseudo-differential operators that simplify the variable-depth Dirichlet–Neumann operator. The model generalizes the Boussinesq system, as it includes the exact dispersion relation in the case of constant depth. Analogous models were proposed by Whitham for unidirectional wave propagation. We first present results for the normal modes and eigenfrequencies of the linearized problem. We see that variable depth introduces effects such as a steepening of the normal modes with the increase in depth variation, and a modulation of the normal mode amplitude. Numerical integration also suggests that the constant depth nonlocal Boussinesq model can capture qualitative features of the evolution obtained with higher order approximations of the Dirichlet–Neumann operator. In the case of variable depth we observe that wave-crests have variable speeds that depend on the depth. We also study the evolutions of Stokes waves initial conditions and observe certain oscillations in width of the crest and also some interesting textures and details in the evolution of wave-crests during the passage over obstacles.