Decay estimates for some semilinear damped hyperbolic problems

Let Ω be a bounded open domain in R ⁿ , g ∶ R → R a non-decreasing continuous function such that g(0)=0 and h ε L _loc ¹ (R⁺; L ²(Ω)). Under suitable assumptions on g and h, the rate of decay of the difference of two solutions is studied for some abstract evolution equations of the general form u ^′′ + Lu + g(u ^′) = h(t,x) as t → + ∞. The results, obtained by use of differential inequalities, can be applied to the case of the semilinear wave equation $$u_u - \Delta u + g{\text{(}}u_t {\text{) = }}h{\text{ in }}R^ + \times \Omega ,{\text{ }}u = {\text{0 on }}R^ + \times \partial \Omega$$ in R ⁺×Ω, u=0 on R ⁺×?Ω. For instance if \(g(s) = c\left| s \right|^{p - 1} s + d\left| s \right|^{q - 1} s\) with c, d>0 and 1 < p≦q, (n?2)q≦n+2, then if \(h \in L^\infty (R + ;L^2 (\Omega ))\) , all solutions are bounded in the energy space for t≧0 and if u, v are two such solutions, the energy norm of u(t) ? v(t) decays like t ^?1/p?1 as t → + ∞.     

Asymptotic analysis of linearly elastic shells. III. Justification of Koiter's shell equations

We consider as in Parts I and II a family of linearly elastic shells of thickness 2?, all having the same middle surfaceS=?(?)?R ³, whereω?R ² is a bounded and connected open set with a Lipschitz-continuous boundary, and? ∈ ?³ (?;R ³). The shells are clamped on a portion of their lateral face, whose middle line is?(γ ₀), whereγ ₀ is a portion of?ω withlength γ ₀>0. For all?>0, let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $u_i^\varepsilon g^{i,\varepsilon }$ of the points of the shell, obtained by solving the three-dimensional problem; let $\zeta _i^\varepsilon$ denote the covariant components of the displacement $\zeta _i^\varepsilon$ a ⁱ of the points of the middle surfaceS, obtained by solving the two-dimensional model ofW.T. Koiter, which consists in finding $$\zeta ^\varepsilon = \left( {\zeta _i^\varepsilon } \right) \in V_K (\omega ) = \left\{ {\eta = (\eta _\iota ) \in {\rm H}^1 (\omega ) \times H^1 (\omega ) \times H^2 (\omega ); \eta _i = \partial _v \eta _3 = 0 on \gamma _0 } \right\}$$ such that $$\begin{gathered} \varepsilon \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \gamma _{\sigma \tau } (\zeta ^\varepsilon )\gamma _{\alpha \beta } (\eta )\sqrt a dy + \frac{{\varepsilon ^3 }}{3} \mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \rho _{\sigma \tau } (\zeta ^\varepsilon )\rho _{\alpha \beta } (\eta )\sqrt a dy \hfill \\ = \mathop \smallint \limits_\omega p^{i,\varepsilon } \eta _i \sqrt a dy for all \eta = (\eta _i ) \in V_K (\omega ), \hfill \\ \end{gathered}$$ where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor ofS, $\gamma _{\alpha \beta }$ (η) and $\rho _{\alpha \beta }$ (η) are the components of the linearized change of metric and change of curvature tensors ofS, and $p^{i,\varepsilon }$ are the components of the resultant of the applied forces. Under the same assumptions as in Part I, we show that the fields $\frac{1}{{2_\varepsilon }}\smallint _{ - \varepsilon }^\varepsilon u_i^\varepsilon g^{i,\varepsilon } dx_3^\varepsilon$ and $\zeta _i^\varepsilon$ a ⁱ , both defined on the surfaceS, have the same principal part as? → 0, inH ¹ (ω) for the tangential components, and inL ²(ω) for the normal component; under the same assumptions as in Part II, we show that the same fields again have the same principal part as? → 0, inH ¹ (ω) for all their components. For “membrane” and “flexural” shells, the two-dimensional model ofW.T. Koiter is therefore justified.     

Asymptotic analysis of linearly elastic shells. I. Justification of membrane shell equations

We consider a family of linearly elastic shells with thickness 2?, clamped along their entire lateral face, all having the same middle surfaceS=φ(?ω) ?R ³, whereω ?R ² is a bounded and connected open set with a Lipschitz-continuous boundaryγ, andφ ∈l ³ ( $\overline \omega$ ;R ³). We make an essential geometrical assumption on the middle surfaceS, which is satisfied ifγ andφ are smooth enough andS is “uniformly elliptic”, in the sense that the two principal radii of curvature are either both>0 at all points ofS, or both<0 at all points ofS. We show that, if the applied body force density isO(1) with respect to?, the fieldtu(?)=(u _i(?)), whereu _i (?) denote the three covariant components of the displacement of the points of the shell given by the equations of three-dimensional elasticity, one “scaled” so as to be defined over the fixed domain Ω=ω×]?1, 1[, converges inH ¹(Ω)×H ¹(Ω)×L ²(Ω) as?→0 to a limitu, which is independent of the transverse variable. Furthermore, the averageξ=1/2ε _?1 ¹ u dx ₃, which belongs to the space $$V_M (\omega ) = H_0^1 (\omega ) \times H_0^1 (\omega ) \times L^2 (\omega ),$$ satisfies the (scaled) two-dimensional equations of a “membrane shell” viz., $$\mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \gamma _{\sigma \tau } (\zeta )\gamma _{\alpha \beta } (\eta ) \sqrt \alpha dy = \mathop \smallint \limits_\omega \left\{ {\mathop \smallint \limits_{ - 1}^1 f^i dx_3 } \right\}\eta _i \sqrt a dy$$ for allη=(η _i) εV _M(ω), where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor of the surfaceS, $$\gamma _{\alpha \beta } (\eta ) = \frac{1}{2}\left( {\partial _{\alpha \eta \beta } + \partial _{\beta \eta \alpha } } \right) - \Gamma _{\alpha \beta }^\sigma \eta _\sigma - b_{\alpha \beta \eta 3} $$ are the components of the linearized change of metric tensor ofS, $\Gamma _{\alpha \beta }^\sigma$ are the Christoffel symbols ofS, $b_{\alpha \beta }$ are the components of the curvature tensor ofS, andf ⁱ are the scaled components of the applied body force. Under the above assumptions, the two-dimensional equations of a “membrane shell” are therefore justified.     

Asymptotic analysis of linearly elastic shells. II. Justification of flexural shell equations

We consider as in Part I a family of linearly elastic shells of thickness 2?, all having the same middle surfaceS=?(?)?R ³, whereω?R ² is a bounded and connected open set with a Lipschitz-continuous boundary, and?∈l ³ (?;R ³). The shells are clamped on a portion of their lateral face, whose middle line is?(γ ₀), whereγ ₀ is any portion of?ω withlength γ ₀>0. We make an essential geometrical assumption on the middle surfaceS and on the setγ ₀, which states that the space of inextensional displacements $$\begin{gathered} V_F (\omega ) = \{ \eta = (\eta _i ) \in H^1 (\omega ) \times H^1 (\omega ) \times H^2 (\omega ); \hfill \\ \eta _i = \partial _v \eta _3 = 0 on \gamma _0 ,\gamma _{\alpha \beta } (\eta ) = 0 in \omega \} , \hfill \\ \end{gathered}$$ where $\gamma _{\alpha \beta }$ (η) are the components of the linearized change is metric tensor ofS, contains non-zero functions. This assumption is satisfied in particular ifS is a portion of cylinder and?(γ ₀) is contained in a generatrix ofS. We show that, if the applied body force density isO(? ²) with respect to?, the fieldu(?)=(u _i (?)), whereu _i (?) denote the three covariant components of the displacement of the points of the shell given by the equations of three-dimensional elasticity, once “scaled” so as to be defined over the fixed domain Ω=ω×]?1, 1[, converges as?→0 inH ¹(Ω) to a limitu, which is independent of the transverse variable. Furthermore, the averageζ=1/2ts _?1 ¹ u dx ₃, which belongs to the spaceV _F (ω), satisfies the (scaled) two-dimensional equations of a “flexural shell”, viz., $$\frac{1}{3}\mathop \smallint \limits_\omega a^{\alpha \beta \sigma \tau } \rho _{\sigma \tau } (\zeta )\rho _{\alpha \beta } (\eta )\sqrt {a } dy = \mathop \smallint \limits_\omega \left\{ {\mathop \smallint \limits_{ - 1}^1 f^i dx_3 } \right\} \eta _i \sqrt {a } dy$$ for allη=(η _i ) ∈V _F (ω), where $a^{\alpha \beta \sigma \tau }$ are the components of the two-dimensional elasticity tensor of the surfaceS, $$\begin{gathered} \rho _{\alpha \beta } (\eta ) = \partial _{\alpha \beta } \eta _3 - \Gamma _{\alpha \beta }^\sigma \partial _\sigma \eta _3 + b_\beta ^\sigma \left( {\partial _\alpha \eta _\sigma - \Gamma _{\alpha \sigma }^\tau \eta _\tau } \right) \hfill \\ + b_\alpha ^\sigma \left( {\partial _\beta \eta _\sigma - \Gamma _{\beta \sigma }^\tau \eta _\tau } \right) + b_\alpha ^\sigma {\text{|}}_\beta \eta _\sigma - c_{\alpha \beta } \eta _3 \hfill \\ \end{gathered} $$ are the components of the linearized change of curvature tensor ofS, $\Gamma _{\alpha \beta }^\sigma$ are the Christoffel symbols ofS, $b_\alpha ^\beta$ are the mixed components of the curvature tensor ofS, andf ⁱ are the scaled components of the applied body force. Under the above assumptions, the two-dimensional equations of a “flexural shell” are therefore justified.     

In the Chapman-Enskog Expansion¶the Burnett Coefficients Satisfy¶the Universal Relation ω3+ω4+θ3= 0

The Chapman–Enskog expansion when applied to a gas of spherical molecules yields formal expressions for the stress deviator P and energy-flux vector q, P=μP ⁽¹⁾+μ² P ⁽²⁾+…, q=μq ⁽¹⁾+μ² q ⁽²⁾+…. The Burnett terms P ⁽²⁾, q ⁽²⁾ depend on 11 coefficients ω _i , 1≦i≦6, θ&; _i , 1≦i≦ 5. This paper shows that ω₃+ω₄+θ₃= 0.     

Differential Systems with Feedback: Time Discretizations and Lyapunov Functions

In this paper we examine a class of Eulerian time discretizations for a monotone cyclic feedback system with a time delay; see Mallet-Paret and Sell (1996a, 1996b) for background information. We construct an integer-valued function V for the discrete-time problem. The Main Theorem shows that V is a Lyapunov function, that is, V(x ⁿ⁺¹)≤V(x ⁿ ) along a solution {x ⁿ }^∞ _n=0, where the time steps can be relatively large.     

Operator method for solving the plane elasticity problem for a strip with periodic cuts

In the present paper, we use the conformal mapping z/c = ζ?2a sin ζ (a, c?const, ζ = u + iv) of the strip {|v| ≤ v ₀, |u| < ∞} onto the domain D, which is a strip with symmetric periodic cuts. For the domain D, in the orthogonal system of isometric coordinates u, v, we solve the plane elasticity problem. We seek the biharmonic function in the form F = C ψ ₀ + S ψ*₀ + x(C ψ ₁ ? S ψ ₂) + y(C ψ ₂ + S ψ ₁), where C(v) and S(v) are the operator functions described in [1] and ψ ₀(u), …, ψ ₂(u) are the desired functions. The boundary conditions for the function F posed for v = ±v ₀ are equivalent to two operator equations for ψ ₁(u) and ψ ₂(u) and to two ordinary differential equations of first order for ψ ₀(u) and ψ*₀(u) [2]. By finding the functions ψ _j (u) in the form of trigonometric series with indeterminate coefficients and by solving the operator equations, we obtain infinite systems of linear equations for the unknown coefficients. We present an efficient method for solving these systems, which is based on studying stable recursive relations. In the present paper, we give an example of analysis of a specific strip (a = 1/4, v ₀ = 1) loaded on the boundary v = v ₀ by a normal load of intensity p. We find the particular solutions corresponding to the extension of the strip by the longitudinal force X ^∞ and to the transverse and pure bending of the strip due to the transverse force Y ^∞ and the constant moment M ^∞, respectively. We also present the graphs of normal and tangential stresses in the transverse cross-section x = 0 and study the stress concentration effect near the cut bottom.     

A Refinement of the Local Serrin-Type Regularity Criterion for a Suitable Weak Solution to the Navier–Stokes Equations

We formulate a new criterion for regularity of a suitable weak solution v to the Navier–Stokes equations at the space-time point (x ₀, t ₀). The criterion imposes a Serrin-type integrability condition on v only in a backward neighbourhood of (x ₀, t ₀), intersected with the exterior of a certain space-time paraboloid with vertex at point (x ₀, t ₀). We make no special assumptions on the solution in the interior of the paraboloid.     

Regularity for Energy-Minimizing Area-Preserving Deformations

In this paper we establish the square integrability of the nonnegative hydrostatic pressure p, that emerges in the minimization problem $$\inf_{\mathcal{K}}\int_{\varOmega}|\nabla \textbf {v}|^2, \quad\varOmega\subset \mathbb {R}^2 $$ as the Lagrange multiplier corresponding to the incompressibility constraint det?v=1 a.e. in Ω. Our method employs the Euler-Lagrange equation for the mollified Cauchy stress C satisfied in the image domain Ω ^?=u(Ω). This allows to construct a convex function ψ, defined in the image domain, such that the measure of the normal mapping of ψ controls the L ² norm of the pressure. As a by-product we conclude that $\textbf {u}\in C^{\frac{1}{2}}_{\textrm {loc}}(\varOmega)$ if the dual pressure (introduced in Karakhanyan, Manuscr. Math. 138:463, 2012) is nonnegative.     

Dynamics as a mechanism preventing the formation of finer and finer microstructure

We study the dynamics of pattern formation in the one-dimensional partial differential equation $$u_u - (W'(u_x ))_x - u_{xxt} + u = 0{\text{ (}}u = u(x,t),{\text{ }}x \in (0,1),{\text{ }}t > 0)$$ proposed recently by Ball, Holmes, James, Pego & Swart [BHJPS] as a mathematical “cartoon” for the dynamic formation of microstructures observed in various crystalline solids. Here W is a double-well potential like 1/4((u _x )² ?1)². What makes this equation interesting and unusual is that it possesses as a Lyapunov function a free energy (consisting of kinetic energy plus a nonconvex “elastic” energy, but no interfacial energy contribution) which does not attain a minimum but favours the formation of finer and finer phase mixtures: $$E[u,u_t ] = \int\limits_0^1 {(\frac{{u_t^2 }}{2} + W(u_x ) + \frac{{u^2 }}{2})dx.}$$ Our analysis of the dynamics confirms the following surprising and striking difference between statics and dynamics, conjectured in [BHJPS] on the basis of numerical simulations of Swart & Holmes [SH]:

?While minimizing the above energy predicts infinitely fine patterns (mathematically: weak but not strong convergence of all minimizing sequences (u _nv_n) of E[u,v] in the Sobolev space W ^{1 p}(0, 1)×L²(0,1)), solutions to the evolution equation of ball et al. typically develop patterns of small but finite length scale (mathematically: strong convergence in W ₁ ^p(0,1)×L²(0,1) of all solutions (u(t),u_t(t)) with low initial energy as time t → ∞).

Moreover, in order to understand the finer details of why the dynamics fails to mimic the behaviour of minimizing sequences and how solutions select their limiting pattern, we present a detailed analysis of the evolution of a restricted class of initial data — those where the strain field u _x has a transition layer structure; our analysis includes proofs that

?at low energy, the number of phases is in fact exactly preserved, that is, there is no nucleation or coarsening

?transition layers lock in and steepen exponentially fast, converging to discontinuous stationary sharp interfaces as time t → ∞

?the limiting patterns — while not minimizing energy globally — are ‘relative minimizers’ in the weak sense of the calculus of variations, that is, minimizers among all patterns which share the same strain interface positions.

     

Minimizers of Caffarelli–Kohn–Nirenberg Inequalities with the Singularity on the Boundary

Let Ω be a bounded smooth domain in ${{\bf R}^N, N\geqq 3}Let Ω be a bounded smooth domain in R^N, N\geqq 3{{\bf R}^N, N\geqq 3}, and D_a^1,2(W){D_a^{1,2}(\Omega)} be the completion of C₀^￥(W){C_0^\infty(\Omega)} with respect to the norm:

||u||a2=òW |x|-2a|?u|2dx.||u||_a^2=\int_\Omega |x|^{-2a}|\nabla u|^2{d}x.      12. Stress and finite-strain analysis of a circular ring under diametral compression      A. J. Durelli  V. J. Parks  T. L. Chen 《Experimental Mechanics》1969,9(5):210-214 Using a moiré, large-strain analysis method, a complete solution is shown in this paper of the fields of strain and stress for a circular ring subjected to diametral compression between two flat platens. The isotheticsu andv, obtained using 1000-lines-per-inch gratings, were differentiated photographically by the shifting technique (moiré-of-moiré) to determine ?u/?x, ?v/?y, ?u/?y and ?v/?x. Using the exact finite strain-displacement relationship, the Eulerian strains ? x E , ? y E and γ xy E were computed. From these, the principal Eulerian strains were obtained. These results were verified with the isochromatics obtained from a large-deformation photoelasticity analysis. The ring was made of a polyurethane rubber which exhibits a linear relationship between natural strain and a newly introduced concept of “natural stress”. The Eulerian strains were converted to natural strains, and from these natural stresses were computed using the newly developed concept. Results are presented graphically for the whole field of the ring.      13. Asymptotic behavior of Toeplitz matrices and determinants      Robert E. Hartwig  Michael E. Fisher 《Archive for Rational Mechanics and Analysis》1969,32(3):190-225 We consider the inverse X N and determinant DN(c) of an N×N Toeplitz matrix CN=[ci?j] 0 N?1 as N ar∞. Under the condition that there exists a monotonic decreasing summable bound b n ≧|c n |+|c ?n |, and that the generating function \(c(\theta ) = \sum\limits_{n = - \infty }^\infty {c_n e^{i{\text{ }}n{\text{ }}\theta } }\) does not vanish, we construct a matrix iterative process which yields (i) explicit asymptotic formulae for the elements of XN when v(c) = (2π)?1 [arg{c(2π)}?arg{c(0)}] is zero. Thence we obtain (ii) expressions for the constants, and bounds on the remainder, in the asymptotic formula $$\ln D_N (c) = N{\text{ }}k_0 (c) + E_0 (c) + E_{1,N} (c) + \mathcal{R}_N (c),$$ and (iii) the extension of this formula to the case of general integral v(c). Under certain further conditions the monotonicity of E1,N+?N is proved. We discuss various identities for DN which apply when c(θ) is a rational function of eiθ and mention a conjecture for D N when c(θ) has zeros, and is discontinuous with arbitrary v(c).      14. A continuum description of partially coherent interfaces      P. H. Leo  J. Hu 《Continuum Mechanics and Thermodynamics》1995,7(1):39-56 A continuum model of a two-phase crystal-crystal system is constructed in which the structure of the interface between the phases is determined by energy minimization, rather than by being specified a priori. The interfacial structure is parameterized by a variable? corresponding to the jump in the surface deformation gradient (or strain) at the interface, so that coherence is defined locally by the condition? = 0. The energy of the system is taken to be the sum of the bulk and interfacial energies, where the interfacial energy densityf xs depends on?. In order to explore how the equilibrium interfacial structure depends on the functionf xs (?), a model system consisting of an elastic film on a rigid substrate is studied, and the interfacial energy density is taken to be nonconvex with a sharp minimum associated with coherence. In this case, it can be shown that the energy of the system is driven to its infimum by separating the interface into coherent and incoherent regions, which may be viewed as a continuum analog to a partially coherent interface. Further, this solution only appears above a certain critical thickness of the film, in agreement with misfit dislocation models of partially coherent interfaces.      15. Regularity of Optimal Maps on the Sphere: the Quadratic Cost and the Reflector Antenna      Grégoire Loeper 《Archive for Rational Mechanics and Analysis》2011,199(1):269-289 Building on the results of Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)), and of the author Loeper (in Acta Math., to appear), we study two problems of optimal transportation on the sphere: the first corresponds to the cost function d 2(x, y), where d(·, ·) is the Riemannian distance of the round sphere; the second corresponds to the cost function ?log |x ? y|, known as the reflector antenna problem. We show that in both cases, the cost-sectional curvature is uniformly positive, and establish the geometrical properties so that the results of Loeper (in Acta Math., to appear) and Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)) can apply: global smooth solutions exist for arbitrary smooth positive data and optimal maps are Hölder continuous under weak assumptions on the data.      16. An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow      Jing Li  Hanfeng Wang  Zhaohui Liu  Sheng Chen  Chuguang Zheng 《Experiments in fluids》2012,53(5):1385-1403 Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y +?<?250 at Re τ (based on the wall friction velocity u τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d p + ?=?1.71, normalized by the fluid kinematic viscosity ν and u τ) are used as dispersed phase, and three low mass loading ratios (Φ m ) ranging from 10?4 to 10?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y +?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y +?>?20).      17. Eulerian formulation of transport equations for three-dimensional shock waves in simple elastic solids      T. C. T. Ting  Yongchi Li 《Journal of Elasticity》1983,13(3):295-310 The propagation of a three-dimensional shock wave in an elastic solid is studied. The material is assumed to be a simple elastic solid in which the Cauchy stress depends on the deformation gradient only. It is shown that the growth or decay of a discontinuity ψ depends on (i) an unknown quantity φ? behind the shock wave, (ii) the two principal curvatures of the shock surface, (iii) the gradient on the shock surface of the shock wave speeds and (iv) the inhomogeneous term which depends on the motion ahead of the shock surface and vanishes when the motion ahead of the shock surface is uniform. If a proper choice is made of the propagation vectorb along which the growth or decay of the discontinuity is measured, the dependence on item (iii) can be avoided. However,b assumes different directions depending on the choice of discontinuity ψ with which one is concerned and the unknown quantity φ? behind the shock wave on which one chooses to depend. As in the case of one-dimensional shock waves, the growth (or decay) of one discontinuity may not be accompanied by the growth (or decay) of other discontinuities. A universal equation relating the growth or decay of discontinuities in the normal stress, normal velocity and specific volume is also presented.      18. Non-piezoactivity in piezoacoustics: basic properties and topological features      V.I. Alshits  V.N. Lyubimov  A. Radowicz 《Archive of Applied Mechanics (Ingenieur Archiv)》2005,74(11-12):739-745 The existence conditions of zero electric fields E and zero electric displacements D are studied for bulk acoustic waves in piezoelectric crystals. General equations are derived for lines of zero electric fields, E(m)=0, and for specific points m 0 of vanishing electric displacements, D(m 0)=0, on the unit sphere of propagation directions m 2=1. The obtained equations are solved for a series of examples of particular crystal symmetry. It is shown that the vectors D α (m) being generally orthogonal to the wave normal m are characterized by definite orientational singularities in the vicinity of m 0 and can be described by the Poincaré indices n=0, ±1 or ±2. The algebraic expressions for the indices n are found both for unrestricted anisotropy and for a series of particular cases.      19. Compatibility Equations of Nonlinear Elasticity for Non-Simply-Connected Bodies      Arash Yavari 《Archive for Rational Mechanics and Analysis》2013,209(1):237-253 Compatibility equations of elasticity are almost 150 years old. Interestingly, they do not seem to have been rigorously studied, to date, for non-simply-connected bodies. In this paper we derive necessary and sufficient compatibility equations of nonlinear elasticity for arbitrary non-simply-connected bodies when the ambient space is Euclidean. For a non-simply-connected body, a measure of strain may not be compatible, even if the standard compatibility equations (“bulk” compatibility equations) are satisfied. It turns out that there may be topological obstructions to compatibility; this paper aims to understand them for both deformation gradient F and the right Cauchy-Green strain C = F T F. We show that the necessary and sufficient conditions for compatibility of deformation gradient F are the vanishing of its exterior derivative and all its periods, that is, its integral over generators of the first homology group of the material manifold. We will show that not every non-null-homotopic path requires supplementary compatibility equations for F and linearized strain e. We then find both necessary and sufficient compatibility conditions for the right Cauchy-Green strain tensor C for arbitrary non-simply-connected bodies when the material and ambient space manifolds have the same dimensions. We discuss the well-known necessary compatibility equations in the linearized setting and the Cesàro-Volterra path integral. We then obtain the sufficient conditions of compatibility for the linearized strain when the body is not simply-connected. To summarize, the question of compatibility reduces to two issues: i) an integrability condition, which is d(F dX) = 0 for the deformation gradient and a curvature vanishing condition for C, and ii) a topological condition. For F dx this is a homological condition because the equation one is trying to solve takes the form dφ = F dX. For C, however, parallel transport is involved, which means that one needs to solve an equation of the form dR/ ds = RK, where R takes values in the orthogonal group. This is, therefore, a question about an orthogonal representation of the fundamental group, which, as the orthogonal group is not commutative, cannot, in general, be reduced to a homological question.      20. Finite Amplitude Transverse Waves in Special Incompressible Viscoelastic Solids      Michael A. Hayes  Giuseppe Saccomandi 《Journal of Elasticity》2000,59(1-3):213-225 We consider the propagation of finite amplitude plane transverse waves in a class of homogeneous isotropic incompressible viscoelastic solids. It is assumed that the Cauchy stress may be written as the sum of an elastic part and a dissipative viscoelastic part. The elastic part is of the form of the stress corresponding to a Mooney–Rivlin material, whereas the dissipative part is a linear combination of A 1, A 1 2 and A 2, where A 1, A 2 are the first and second Rivlin–Ericksen tensors. The body is first subject to a homogeneous static deformation. It is seen that two finite amplitude transverse plane waves may propagate in every direction in the deformed body. It is also seen that a finite amplitude circularly polarized wave may propagate along either n + or n ?, where n +, n ? are the normals to the planes of the central circular section of the ellipsoid x?B ?1 x=1. Here B is the left Cauchy–Green strain tensor corresponding to the finite static homogeneous deformation.

Stress and finite-strain analysis of a circular ring under diametral compression

Using a moiré, large-strain analysis method, a complete solution is shown in this paper of the fields of strain and stress for a circular ring subjected to diametral compression between two flat platens. The isotheticsu andv, obtained using 1000-lines-per-inch gratings, were differentiated photographically by the shifting technique (moiré-of-moiré) to determine ?u/?x, ?v/?y, ?u/?y and ?v/?x. Using the exact finite strain-displacement relationship, the Eulerian strains ? _x ^E , ? _y ^E and γ _xy ^E were computed. From these, the principal Eulerian strains were obtained. These results were verified with the isochromatics obtained from a large-deformation photoelasticity analysis. The ring was made of a polyurethane rubber which exhibits a linear relationship between natural strain and a newly introduced concept of “natural stress”. The Eulerian strains were converted to natural strains, and from these natural stresses were computed using the newly developed concept. Results are presented graphically for the whole field of the ring.     

Asymptotic behavior of Toeplitz matrices and determinants

We consider the inverse X _N and determinant D_N(c) of an N×N Toeplitz matrix C_N=[c_i?j] ₀ ^N?1 as N ar∞. Under the condition that there exists a monotonic decreasing summable bound b _n ≧|c _n |+|c _?n |, and that the generating function \(c(\theta ) = \sum\limits_{n = - \infty }^\infty {c_n e^{i{\text{ }}n{\text{ }}\theta } }\) does not vanish, we construct a matrix iterative process which yields (i) explicit asymptotic formulae for the elements of X_N when v(c) = (2π)^?1 [arg{c(2π)}?arg{c(0)}] is zero. Thence we obtain (ii) expressions for the constants, and bounds on the remainder, in the asymptotic formula $$\ln D_N (c) = N{\text{ }}k_0 (c) + E_0 (c) + E_{1,N} (c) + \mathcal{R}_N (c),$$ and (iii) the extension of this formula to the case of general integral v(c). Under certain further conditions the monotonicity of E_1,N+?_N is proved. We discuss various identities for D_N which apply when c(θ) is a rational function of e^iθ and mention a conjecture for D _N when c(θ) has zeros, and is discontinuous with arbitrary v(c).     

A continuum description of partially coherent interfaces

A continuum model of a two-phase crystal-crystal system is constructed in which the structure of the interface between the phases is determined by energy minimization, rather than by being specified a priori. The interfacial structure is parameterized by a variable? corresponding to the jump in the surface deformation gradient (or strain) at the interface, so that coherence is defined locally by the condition? = 0. The energy of the system is taken to be the sum of the bulk and interfacial energies, where the interfacial energy densityf ^xs depends on?. In order to explore how the equilibrium interfacial structure depends on the functionf ^xs (?), a model system consisting of an elastic film on a rigid substrate is studied, and the interfacial energy density is taken to be nonconvex with a sharp minimum associated with coherence. In this case, it can be shown that the energy of the system is driven to its infimum by separating the interface into coherent and incoherent regions, which may be viewed as a continuum analog to a partially coherent interface. Further, this solution only appears above a certain critical thickness of the film, in agreement with misfit dislocation models of partially coherent interfaces.     

Regularity of Optimal Maps on the Sphere: the Quadratic Cost and the Reflector Antenna

Building on the results of Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)), and of the author Loeper (in Acta Math., to appear), we study two problems of optimal transportation on the sphere: the first corresponds to the cost function d ²(x, y), where d(·, ·) is the Riemannian distance of the round sphere; the second corresponds to the cost function ?log |x ? y|, known as the reflector antenna problem. We show that in both cases, the cost-sectional curvature is uniformly positive, and establish the geometrical properties so that the results of Loeper (in Acta Math., to appear) and Ma et al. (in Arch. Rational Mech. Anal. 177(2), 151–183 (2005)) can apply: global smooth solutions exist for arbitrary smooth positive data and optimal maps are Hölder continuous under weak assumptions on the data.     

An experimental study on turbulence modification in the near-wall boundary layer of a dilute gas-particle channel flow

Turbulence modifications of a dilute gas-particle flow are experimentally investigated in the lower boundary layer of a horizontal channel by means of a simultaneous two-phase PIV measurement technique. The measurements are conducted in the near-wall region with y ⁺?<?250 at Re _τ (based on the wall friction velocity u _τ and half channel height h)?=?430. High spatial resolution and small interrogation window are used to minimize the PIV measurement uncertainty due to the velocity gradient near the wall. Polythene beads with the diameter of 60?μm (d _p ⁺ ?=?1.71, normalized by the fluid kinematic viscosity ν and u _τ) are used as dispersed phase, and three low mass loading ratios (Φ _m ) ranging from 10^?4 to 10^?3 are tested. It is found that the addition of the particles noticeably modifies the mean velocity and turbulent intensities of the gas-phase, as well as the turbulence coherent structures, even at Φ _m ?=?0.025?%. Particle inertia changes the viscous sublayer of the gas turbulence with a smaller thickness and a larger streamwise velocity gradient, which increases the peak value of the streamwise fluctuation velocity ( $ u_{\text{rms}}^{ + } $ ) of the gas-phase with its location shifting to the wall. Particle sedimentation increases the roughness of the bottom wall, which significantly increases the wall-normal fluctuation velocity ( $ v_{\text{rms}}^{ + } $ ) and Reynolds shear stress ( $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ ) of the gas-phase in the inner region of the boundary layer (y ⁺?<?10). Under effect of particle–wall collision, the Q2 events (ejections) of the gas-phase are slightly increased by particles, while the Q4 events (sweeps) are obviously decreased. The spatial scale of the coherent structures near the wall shrinks remarkably with the presence of the particles, which may be attributed to the intensified crossing-trajectory effects due to particle saltation near the bottom wall. Meanwhile, the $ v_{\text{rms}}^{ + } $ and $ - \langle u^{ \prime } v^{\prime } \rangle^{ + } $ of the gas-phase are significantly reduced in the outer region of the boundary layer (y ⁺?>?20).     

Eulerian formulation of transport equations for three-dimensional shock waves in simple elastic solids

The propagation of a three-dimensional shock wave in an elastic solid is studied. The material is assumed to be a simple elastic solid in which the Cauchy stress depends on the deformation gradient only. It is shown that the growth or decay of a discontinuity ψ depends on (i) an unknown quantity φ^? behind the shock wave, (ii) the two principal curvatures of the shock surface, (iii) the gradient on the shock surface of the shock wave speeds and (iv) the inhomogeneous term which depends on the motion ahead of the shock surface and vanishes when the motion ahead of the shock surface is uniform. If a proper choice is made of the propagation vectorb along which the growth or decay of the discontinuity is measured, the dependence on item (iii) can be avoided. However,b assumes different directions depending on the choice of discontinuity ψ with which one is concerned and the unknown quantity φ^? behind the shock wave on which one chooses to depend. As in the case of one-dimensional shock waves, the growth (or decay) of one discontinuity may not be accompanied by the growth (or decay) of other discontinuities. A universal equation relating the growth or decay of discontinuities in the normal stress, normal velocity and specific volume is also presented.     

Non-piezoactivity in piezoacoustics: basic properties and topological features

The existence conditions of zero electric fields E and zero electric displacements D are studied for bulk acoustic waves in piezoelectric crystals. General equations are derived for lines of zero electric fields, E(m)=0, and for specific points m ₀ of vanishing electric displacements, D(m ₀)=0, on the unit sphere of propagation directions m ²=1. The obtained equations are solved for a series of examples of particular crystal symmetry. It is shown that the vectors D _α (m) being generally orthogonal to the wave normal m are characterized by definite orientational singularities in the vicinity of m ₀ and can be described by the Poincaré indices n=0, ±1 or ±2. The algebraic expressions for the indices n are found both for unrestricted anisotropy and for a series of particular cases.     

Compatibility Equations of Nonlinear Elasticity for Non-Simply-Connected Bodies

Compatibility equations of elasticity are almost 150 years old. Interestingly, they do not seem to have been rigorously studied, to date, for non-simply-connected bodies. In this paper we derive necessary and sufficient compatibility equations of nonlinear elasticity for arbitrary non-simply-connected bodies when the ambient space is Euclidean. For a non-simply-connected body, a measure of strain may not be compatible, even if the standard compatibility equations (“bulk” compatibility equations) are satisfied. It turns out that there may be topological obstructions to compatibility; this paper aims to understand them for both deformation gradient F and the right Cauchy-Green strain C = F ^T F. We show that the necessary and sufficient conditions for compatibility of deformation gradient F are the vanishing of its exterior derivative and all its periods, that is, its integral over generators of the first homology group of the material manifold. We will show that not every non-null-homotopic path requires supplementary compatibility equations for F and linearized strain e. We then find both necessary and sufficient compatibility conditions for the right Cauchy-Green strain tensor C for arbitrary non-simply-connected bodies when the material and ambient space manifolds have the same dimensions. We discuss the well-known necessary compatibility equations in the linearized setting and the Cesàro-Volterra path integral. We then obtain the sufficient conditions of compatibility for the linearized strain when the body is not simply-connected. To summarize, the question of compatibility reduces to two issues: i) an integrability condition, which is d(F dX) = 0 for the deformation gradient and a curvature vanishing condition for C, and ii) a topological condition. For F dx this is a homological condition because the equation one is trying to solve takes the form dφ = F dX. For C, however, parallel transport is involved, which means that one needs to solve an equation of the form dR/ ds = RK, where R takes values in the orthogonal group. This is, therefore, a question about an orthogonal representation of the fundamental group, which, as the orthogonal group is not commutative, cannot, in general, be reduced to a homological question.     

Finite Amplitude Transverse Waves in Special Incompressible Viscoelastic Solids

We consider the propagation of finite amplitude plane transverse waves in a class of homogeneous isotropic incompressible viscoelastic solids. It is assumed that the Cauchy stress may be written as the sum of an elastic part and a dissipative viscoelastic part. The elastic part is of the form of the stress corresponding to a Mooney–Rivlin material, whereas the dissipative part is a linear combination of A ₁, A ₁ ² and A ₂, where A ₁, A ₂ are the first and second Rivlin–Ericksen tensors. The body is first subject to a homogeneous static deformation. It is seen that two finite amplitude transverse plane waves may propagate in every direction in the deformed body. It is also seen that a finite amplitude circularly polarized wave may propagate along either n ⁺ or n ^?, where n ⁺, n ^? are the normals to the planes of the central circular section of the ellipsoid x?B ^?1 x=1. Here B is the left Cauchy–Green strain tensor corresponding to the finite static homogeneous deformation.