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.     

A gradient based iterative algorithm for solving structural dynamics model updating problems

The procedure of updating an existing but inaccurate model is an essential step toward establishing an effective model. Updating damping and stiffness matrices simultaneously with measured modal data can be mathematically formulated as following two problems. Problem 1: Let M _a ∈SR ^n×n be the analytical mass matrix, and Λ=diag{λ ₁,…,λ _p }∈C ^p×p , X=[x ₁,…,x _p ]∈C ^n×p be the measured eigenvalue and eigenvector matrices, where rank(X)=p, p<n and both Λ and X are closed under complex conjugation in the sense that $\lambda_{2j} = \bar{\lambda}_{2j-1} \in\nobreak{\mathbf{C}} $ , $x_{2j} = \bar{x}_{2j-1} \in{\mathbf{C}}^{n} $ for j=1,…,l, and λ _k ∈R, x _k ∈R ⁿ for k=2l+1,…,p. Find real-valued symmetric matrices D and K such that M _a XΛ ²+DXΛ+KX=0. Problem 2: Let D _a ,K _a ∈SR ^n×n be the analytical damping and stiffness matrices. Find $(\hat{D}, \hat{K}) \in\mathbf{S}_{\mathbf{E}}$ such that $\| \hat{D}-D_{a} \|^{2}+\| \hat{K}-K_{a} \|^{2}= \min_{(D,K) \in \mathbf{S}_{\mathbf{E}}}(\| D-D_{a} \|^{2} +\|K-K_{a} \|^{2})$ , where S _E is the solution set of Problem 1 and ∥?∥ is the Frobenius norm. In this paper, a gradient based iterative (GI) algorithm is constructed to solve Problems 1 and 2. A sufficient condition for the convergence of the iterative method is derived and the range of the convergence factor is given to guarantee that the iterative solutions consistently converge to the unique minimum Frobenius norm symmetric solution of Problem 2 when a suitable initial symmetric matrix pair is chosen. The algorithm proposed requires less storage capacity than the existing numerical ones and is numerically reliable as only matrix manipulation is required. Two numerical examples show that the introduced iterative algorithm is quite efficient.     

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. 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.     

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).     

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.     

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.     

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 → + ∞.     

A symmetry analysis of some classes of evolutionary nonlinear (2+1)-diffusion equations with variable diffusivity

In this paper the (2+1)-nonlinear diffusion equation u _t ?div(f(u)grad u)=0 with variable diffusivity is considered. Using the Lie method, a complete symmetry classification of the equation is presented. Reductions, via two-dimensional Lie subalgebras of the equation, to first- or second-order ordinary differential equations are given. In a few interesting cases exact solutions are presented.     

On the Nonhomogeneous Boundary Value Problem for the Navier–Stokes System in a Class of Unbounded Domains

The stationary Navier–Stokes system with nonhomogeneous boundary conditions is studied in a class of domains Ω having “paraboloidal” outlets to infinity. The boundary ${\partial\Omega}$ is multiply connected and consists of M infinite connected components S _m , which form the outer boundary, and I compact connected components Γ _i forming the inner boundary Γ. The boundary value a is assumed to have a compact support and it is supposed that the fluxes of a over the components Γ _i of the inner boundary are sufficiently small. We do not pose any restrictions on fluxes of a over the infinite components S _m . The existence of at least one weak solution to the Navier–Stokes problem is proved. The solution may have finite or infinite Dirichlet integral depending on geometrical properties of outlets to infinity.     

On the use of interaction tensors to describe and predict rod interactions in rod suspensions

Recently, Férec et al. (2009a) proposed a model for nondilute rod-like suspensions, where particle interactions are taking into account via a micromechanical approach. The derived governing equation used the well-known second- and fourth-order orientation tensors (a ₂ and a ₄ ) and novel second- and fourth-order interaction tensors (b ₂ and b ₄ ). To completely close the model, it is necessary to express a ₄ , b ₂ , and b ₄ in terms of a ₂ . This paper gives the general framework to elaborate these new relations. Firstly, approximations for b ₂ are developed based on linear combinations of a ₂ and a ₄ . Moreover, a new closure approximation is also derived for b ₄ , based on the orthotropic fitted closure approach. Unknown parameters are determined by a least-square fitting technique with assumed exact solutions constructed from the probability distribution function (PDF). As numerical solutions for the PDF are difficult to obtain given the nonlinearity of the problem, a combination of steady state solutions is used to generate PDF designed to cover uniformly the entire domain of possible orientations. All these proposed approximations are tested against the particle-based simulations in a variety of flow fields. Improvements of the different approximations are observed, and the couple iORW-CO4P3 gives efficient results.     

A thermodynamic approach to isotropic-nematic phase transitions in liquid crystals

We propose a dynamical model for (non-isothermal) phase transitions in liquid crystals. Macroscopic motions of the liquid crystal (LC) are neglected, while the coupling with the electromagnetic field is considered. The LC is described in terms of the classical order tensor Q, which is split as Q=s N, where N is a normalized tensor; two independent evolution laws are given for s and N. The model includes an evolutive equation for the temperature field obtained from an appropriate form of the energy balance, in which the internal powers associated to the equations for s and N are accounted for. The thermodynamic restrictions in the constitutive relations which ensure the Clausius–Duhem inequality have been pointed out.     

Singular PDFs of a Dispersing Scalar in Turbulence

When a finite quantity of scalar of uniform concentration C _m is injected at time t= 0 into a turbulent flow containing no scalar, the initial PDF (probability density function) of the concentration of the scalar is a weighted sum of two delta functions. Because of molecular mixing, p _C must eventually become a single delta function. In an earlier paper (Chatwin, Eur. J. Appl. Math. 13 (2002) 95-108). it was suggested that the evolution from the initial state could take place with p _C having modes at its end-points that merged into one another as time increased. This is different from what has normally been supposed. The present paper extends an examination of two families of PDFs that first appeared in Kowe and Chatwin (J. Eng. Math. 19 (1985) 217-231); one of these was obtained again in two separate investigations (Chatwin, Eur. J. Appl. Math. 13 (2002) 95-108; Chatwin and Zimmerman, Environmetrics 9 (1998) 131-138). Here the two families are extended, but in a preliminary and idealized way, by allowing an isotropic spatial dependence and by including the diffusion and advection terms, present in the evolution equation for p _C, but assumed to vanish identically in Chatwin (Eur. J. Appl Math. 13 (2002) 95-108). It is shown that there is consistency, thereby supporting a more detailed investigation using DNS.     

Finite-dimensional asymptotic behavior of some semilinear damped hyperbolic problems

We prove that the solution semigroup $$S_t \left[ {u_0 ,v_0 } \right] = \left[ {u(t),u_t (t)} \right]$$ generated by the evolutionary problem $$\left\{ P \right\}\left\{ \begin{gathered} u_{tt} + g(u_t ) + Lu + f(u) = 0, t \geqslant 0 \hfill \\ u(0) = u_0 , u_t (0) = \upsilon _0 \hfill \\ \end{gathered} \right.$$ possesses a global attractorA in the energy spaceE _o=V×L ²(Ω). Moreover,A is contained in a finite-dimensional inertial setA attracting bounded subsets ofE ₁=D(L)×V exponentially with growing time.     

Existence of Time-Periodic Solutions for a Magnetoelastic System in Bounded Domains

We investigate the existence of periodic solutions for a semilinear (nonlinearly coupled) magnetoelastic system in bounded, simply connected, three-dimensional domains with boundaries of class C ². The mathematical model includes a nonlinear mechanical dissipation like ρ(u′)=|u′| ^p u′ and a periodic forcing function of period T. We prove the existence of T-periodic weak solutions when p∈[3,4] (p=0 being a simpler case). In the corresponding two-dimensional case, the existence result holds under the assumption that p≥2.     

On the Brézis–Nirenberg Problem

We study the following Brézis–Nirenberg problem (Comm Pure Appl Math 36:437–477, 1983): $-\Delta u=\lambda u+ |u|^{2^\ast-2}u, \quad u\in H_0^1(\Omega),$ where Ω is a bounded smooth domain of R ^N (N ≧ 7) and 2* is the critical Sobolev exponent. We show that, for each fixed λ > 0, this problem has infinitely many sign-changing solutions. In particular, if λ ≧ λ₁, the Brézis–Nirenberg problem has and only has infinitely many sign-changing solutions except zero. The main tool is the estimates of Morse indices of nodal solutions.     

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.     

A Concentration Phenomenon for Semilinear Elliptic Equations

For a domain ${\Omega \subset \mathbb{R}^{N}}$ we consider the equation $$-\Delta{u} + V(x)u = Q_n(x)|{u}|^{p-2}u$$ with zero Dirichlet boundary conditions and ${p\in(2, 2^*)}$ . Here ${V \geqq 0}$ and Q _n are bounded functions that are positive in a region contained in ${\Omega}$ and negative outside, and such that the sets {Q _n  > 0} shrink to a point ${x_0 \in \Omega}$ as ${n \to \infty}$ . We show that if u _n is a nontrivial solution corresponding to Q _n , then the sequence (u _n ) concentrates at x ₀ with respect to the H ¹ and certain L ^q -norms. We also show that if the sets {Q _n  > 0} shrink to two points and u _n are ground state solutions, then they concentrate at one of these points.