Dependence of Topological Conjugacies on Parameters

For impulsive differential equations, we establish the existence of invariant stable manifolds under sufficiently small perturbations of a linear equation. We consider the general case of nonautonomous equations for which the linear part has a nonuniform exponential dichotomy. One of the main advantages of our work is that our results are optimal, in the sense that for vector fields of class C ¹ outside the jumping times, we show that the invariant manifolds are also of class C ¹ outside these times. The novelty of our proof is the use of the fiber contraction principle to establish the smoothness of the invariant manifolds. In addition, using the same approach we can also consider linear perturbations.     

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.     

Gradient Systems with Wiggly Energies¶and Related Averaging Problems

Gradient systems with wiggly energies of the form

$$

and A:? ^d →? wereproposed by Abeyaratne, Chu &; James [2] to study the kinetics of martensitic phase transitions. Their model may be recast in the framework of the theory of averaging as a dynamical system on ? ^d ×? ^d , with the slow variable x∈? ^d and fast variable θ∈? ^d . However, this problem lies completely outside the classical theory of averaging, since the vertical flow on ? ^d is not ergodic for sets of positive measure, and we must interpret averages to mean weak limits.

We obtain rigorous averaging results for d= 2. We use Schwartz's generalization of the Poincaré-Bendixson theorem [37] to heuristically derive homogenized equations for the weak limits. These equations depend on the ω-limit sets for the vertical flow on fibres. When the vertical flow is structurally stable, we use the persistence of hyperbolic structures to prove that these are the correct equations. We combine these theorems with a study of two-parameter bifurcations of flows on ?² to characterize the weak limits. Our results may be interpreted as follows. The space ?² breaks into: (˙1) a bounded open set surrounding {?F ^?1 (0)} where there is only sticking, (˙2) a transition region outside this set, where the dynamics is a combination of sticking and slipping, and (˙3) the rest of the plane, which contains a countable number of resonance zones, with nonempty interior, and their nowhere dense complement. Inside a resonance zone the direction of the weak limits is given by the rotation number ρ∈?. The Cantor set structure of the resonance zones is described by well-known results of Arnol'd [7] and Herman [27] in the theory of circle diffeomorphisms. Consequently, the homogenized equations vary on all scales. We also study the linear transport equation associated with the wiggly gradient flow, and show that its homogenization limit is not well posed.Smyshlyaev has studied this problem independently, and some of our results are similar [39].     

An Eulerian–Lagrangian Form for the Euler Equations in Sobolev Spaces

In 2000 Constantin showed that the incompressible Euler equations can be written in an “Eulerian–Lagrangian” form which involves the back-to-labels map (the inverse of the trajectory map for each fixed time). In the same paper a local existence result is proved in certain Hölder spaces \({C^{1,\mu}}\). We review the Eulerian–Lagrangian formulation of the equations and prove that given initial data in H ^s for \({n \geq 2}\) and \({s > \frac{n}{2}+1}\), a unique local-in-time solution exists on the n-torus that is continuous into H ^s and C ¹ into H ^s-1. These solutions automatically have C ¹ trajectories. The proof here is direct and does not appeal to results already known about the classical formulation. Moreover, these solutions are regular enough that the classical and Eulerian–Lagrangian formulations are equivalent, therefore what we present amounts to an alternative approach to some of the standard theory.     

Stability of plane oscillations of a satellite in a circular orbit

We study motions of a rigid body (a satellite) about the center of mass in a central Newtonian gravitational field in a circular orbit. There is a known particular motion of the satellite in which one of its principal central axes of inertia is perpendicular to the orbit plane and the satellite itself exhibits plane pendulum-like oscillations about this axis. Under the assumption that the satellite principal central moments of inertia A, B, and C satisfy the relation B = A + C corresponding to the case of a thin plate, we perform rigorous nonlinear analysis of the orbital stability of this motion.In the plane of the problem parameters, namely, the oscillation amplitude ε and the inertial parameter, there exist countably many domains of orbital stability of the satellite oscillations in the linear approximation. Nonlinear orbital stability analysis was carried out in thirteen of these domains. Isoenergetic reduction of the system of equations of the perturbed motion is performed at the energy level corresponding to the unperturbed periodic motion. Further, using the algorithm developed in [1], we construct the symplectic mapping generated by the equations of the reduced system, normalize it, and analyze the stability. We consider resonance and nonresonance cases. For small values of the oscillation amplitude, we perform analytic investigations; for arbitrary values of ε, numerical analysis is used.Earlier, numerical analysis of stability of plane pendulum-like motions of a satellite in a circular orbit was performed in several special cases in [1–4].     

The Moments Sliding Vector Field on the Intersection of Two Manifolds

In this work, we consider a special choice of sliding vector field on the intersection of two co-dimension 1 manifolds. The proposed vector field, which belongs to the class of Filippov vector fields, will be called moments vector field and we will call moments trajectory the associated solution trajectory. Our main result is to show that the moments vector field is a well defined, and smoothly varying, Filippov sliding vector field on the intersection \(\Sigma \) of two discontinuity manifolds, under general attractivity conditions of \(\Sigma \). We also examine the behavior of the moments trajectory at first order exit points, and show that it exits smoothly at these points. Numerical experiments illustrate our results and contrast the present choice with other choices of Filippov sliding vector field.     

The stress state of a plate subjected to a piecewise homogeneous load

We consider the stress-strain state of a plate having a doubly connected domain S bounded from the outside by a circle of radius R and from the inside by an ellipse with two rectilinear cuts. The cuts lie symmetrically on the x-axis. The plate is subjected to various forces: the hole contour (the ellipse) is under the action of uniformly distributed forces of intensity q, and the cut shores are free of loads; at the points ±ib of the imaginary axis, the plate is under the action of a lumped force P.The solution of the problem is reduced to determining two analytic functions φ(z) and ψ(z) satisfying certain boundary conditions (depending on the type of the acting loads).We use the Kolosov-Muskhelishvili method to reduce the problem to a system of linear algebraic equations for the coefficients in the expansions of the functions φ(z) and ψ(z). The solution thus obtained is illustrated by numerical examples.     

Long-time Stability in Systems of Conservation Laws,Using Relative Entropy/Energy

We study the long-time stability of shock-free solutions of hyperbolic systems of conservation laws, under an arbitrarily large initial disturbance in L ²∩ L ^∞. We use the relative entropy method, a robust tool which allows us to consider rough and large disturbances. We display practical examples in several space dimensions, for scalar equations as well as isentropic gas dynamics. For full gas dynamics, we use a trick from Chen [1], in which the estimate is made in terms of the relative mechanical energy instead of the relative mathematical entropy.     

Entire Solutions to Equivariant Elliptic Systems with Variational Structure

We consider the system Δu ? W _u (u) = 0, where \({u : \mathbb{R}^n \to \mathbb{R}^n}\) , for a class of potentials \({W : \mathbb{R}^n \to \mathbb{R}}\) that possess several global minima and are invariant under a general finite reflection group G. We establish existence of nontrivial G-equivariant entire solutions connecting the global minima of W along certain directions at infinity.     

Entire Solutions of Superlinear Problems with Indefinite Weights and Hardy Potentials

We provide the structure of regular/singular fast/slow decay radially symmetric solutions for a class of superlinear elliptic equations with an indefinite weight. In particular we are interested in the case where such a weight is positive in a ball and negative outside, or in the reversed situation. We extend the approach to elliptic equations in presence of Hardy potentials, i.e. to

$$\begin{aligned} \varDelta u +\frac{h(|\text {x}|)}{|\text {x}|^2} u+ f(u, |\text {x}|)=0 \end{aligned}$$

where h is not necessarily constant. By the use of Fowler transformation we study the corresponding dynamical systems, presenting the construction of invariant manifolds when the global existence of solutions is not ensured.

     

Miscible Thermo-Viscous Fingering Instability in Porous Media. Part 1: Linear Stability Analysis

The development of the thermo-viscous fingering instability of miscible displacements in homogeneous porous media is examined. In this first part of the study dealing with stability analysis, the basic equations and the parameters governing the problem in a rectilinear geometry are developed. An exponential dependence of viscosity on temperature and concentration is represented by two parameters, thermal mobility ratio β _T and a solutal mobility ratio β _C , respectively. Other parameters involved are the Lewis number Le and a thermal-lag coefficient λ. The governing equations are linearized and solved to obtain instability characteristics using either a quasi-steady-state approximation (QSSA) or initial value calculations (IVC). Exact analytical solutions are also obtained for very weakly diffusing systems. Using the QSSA approach, it was found that an increase in thermal mobility ratio β _T is seen to enhance the instability for fixed β _C , Le and λ. For fixed β _C and β _T , a decrease in the thermal-lag coefficient and/or an increase in the Lewis number always decrease the instability. Moreover, strong thermal diffusion at large Le as well as enhanced redistribution of heat between the solid and fluid phases at small λ is seen to alleviate the destabilizing effects of positive β _T . Consequently, the instability gets strictly dominated by the solutal front. The linear stability analysis using IVC approach leads to conclusions similar to the QSSA approach except for the case of large Le and unity λ flow where the instability is seen to get even less pronounced than in the case of a reference isothermal flow of the same β _C , but β _T  = 0. At practically, small value of λ, however, the instability ultimately approaches that due to β _C only.     

Decay rates of higher-order norms of solutions to the Navier-Stokes-Landau-Lifshitz system

In this paper, we investigate a system of the incompressible Navier-Stokes equations coupled with Landau-Lifshitz equations in three spatial dimensions. Under the assumption of small initial data, we establish the global solutions with the help of an energy method. Furthermore, we obtain the time decay rates of the higher-order spatial derivatives of the solutions by applying a Fourier splitting method introduced by Schonbek(SCHONBEK, M. E. L~2 decay for weak solutions of the Navier-Stokes equations. Archive for Rational Mechanics and Analysis, 88, 209–222(1985)) under an additional assumption that the initial perturbation is bounded in L~1(R~3).     

Global Attractors of Sixth Order PDEs Describing the Faceting of Growing Surfaces

A spatially two-dimensional sixth order PDE describing the evolution of a growing crystalline surface h(x, y, t) that undergoes faceting is considered with periodic boundary conditions, as well as its reduced one-dimensional version. These equations are expressed in terms of the slopes \(u_1=h_{x}\) and \(u_2=h_y\) to establish the existence of global, connected attractors for both equations. Since unique solutions are guaranteed for initial conditions in \(\dot{H}^2_{per}\), we consider the solution operator \(S(t): \dot{H}^2_{per} \rightarrow \dot{H}^2_{per}\), to gain our results. We prove the necessary continuity, dissipation and compactness properties.     

Linearized Stability for Semiflows Generated by a Class of Neutral Equations,with Applications to State-Dependent Delays

We prove a principle of linearized stability for semiflows generated by neutral functional differential equations of the form x′(t) = g(? x _t , x _t ). The state space is a closed subset in a manifold of C ²-functions. Applications include equations with state-dependent delay, as for example x′(t) = a x′(t + d(x(t))) + f (x(t + r(x(t)))) with \({a\in\mathbb{R}, d:\mathbb{R}\to(-h,0), f:\mathbb{R}\to\mathbb{R}, r:\mathbb{R}\to[-h,0]}\).     

On the interplay between fast reaction and slow diffusion in the concrete carbonation process: a matched-asymptotics approach

A matched-asymptotics approach is proposed to show the occurrence of two distinct characteristic length scales in the carbonation process. The separation of these scales arises due to the strong competition between reaction and diffusion effects. We show that for sufficiently large times τ the width of the carbonated region is proportional to \(\sqrt{\tau}\), while the width of the reaction front is proportional to \(\tau^{\frac{p-1}{2(p+1)}}\) for carbonation-reaction rates with a power law structure like k[CO₂] ^p [Ca(OH)₂] ^q , where k>0 and p,q>1 and identify the proportionality coefficient asymptotically. We emphasize the occurrence of a water barrier in the reaction zone which may hinder the penetration of CO₂ by locally filling with water air parts of the pores. This non-linear effect may be one of the causes why a purely linear extrapolation of accelerated carbonation test results to natural carbonation settings is (even theoretically) not reasonable. Finally, we compare our asymptotic penetration law against measured penetration depths from Bune (Zum Karbonatisierungsbedingten Verlust der Dauerhaftigkeit von Außenbauteilen aus Stahlbeton, 1994). The novelty consists in the fact that the factor multiplying \(\sqrt{\tau}\) is now identified asymptotically by solving a non-linear system of ordinary differential equations, and hence, fitting arguments are not necessary to estimate its size. We offer an alternative to the (asymptotic) \(\sqrt{\tau}\) expression of the carbonation-front position obtained in Papadakis et al. (AIChE J. 35:1639, 1989).     

Beltrami Fields with a Nonconstant Proportionality Factor are Rare

We consider the existence of Beltrami fields with a nonconstant proportionality factor f in an open subset U of \({\mathbb{R}^3}\). By reformulating this problem as a constrained evolution equation on a surface, we find an explicit differential equation that f must satisfy whenever there is a nontrivial Beltrami field with this factor. This ensures that there are no nontrivial regular solutions for an open and dense set of factors f in the C^k topology, \({k\geqq 7}\). In particular, there are no nontrivial Beltrami fields whenever f has a regular level set diffeomorphic to the sphere. This provides an explanation of the helical flow paradox of Morgulis et al. (Commun Pure Appl Math 48:571–582, 1995).     

Enhanced Dissipation and Inviscid Damping in the Inviscid Limit of the Navier–Stokes Equations Near the Two Dimensional Couette Flow

In this work we study the long time inviscid limit of the two dimensional Navier–Stokes equations near the periodic Couette flow. In particular, we confirm at the nonlinear level the qualitative behavior predicted by Kelvin’s 1887 linear analysis. At high Reynolds number Re, we prove that the solution behaves qualitatively like two dimensional Euler for times \({{t \lesssim Re^{1/3}}}\), and in particular exhibits inviscid damping (for example the vorticity weakly approaches a shear flow). For times \({{t \gtrsim Re^{1/3}}}\), which is sooner than the natural dissipative time scale O(Re), the viscosity becomes dominant and the streamwise dependence of the vorticity is rapidly eliminated by an enhanced dissipation effect. Afterwards, the remaining shear flow decays on very long time scales \({{t \gtrsim Re}}\) back to the Couette flow. When properly defined, the dissipative length-scale in this setting is \({{\ell_D \sim Re^{-1/3}}}\), larger than the scale \({{\ell_D \sim Re^{-1/2}}}\) predicted in classical Batchelor–Kraichnan two dimensional turbulence theory. The class of initial data we study is the sum of a sufficiently smooth function and a small (with respect to Re^?1) L² function.     

The Failure of Rank-One Connections

This article is concerned with interface problems for Lipschitz mappings f ₊:? ⁿ ₊→? ⁿ and f _?:? ⁿ _?→? ⁿ in the half spaces, which agree on the common boundary ? ^{n ? 1}=?? ⁿ ₊=?? ⁿ _?. These naturally occur in mathematical models for material microstructures and crystals. The task is to determine the relationship between the sets of values of the differentials Df ₊ and Df _?. For some time it has been thought that the polyconvex hulls [Df ₊] ^pc and [Df _?] ^pc satisfy Hadamard's jump condition or are at least rank-one connected. Our examples here refute this idea.The estimates of the Jacobians we obtain in the course of solving the so-called Monge-Ampère inequalities seem also to be of independent interest. As an application, we construct uniformly elliptic systems of first order partial differential equations in the same homotopy class as the familiar Cauchy-Riemann equations, for which the unique continuation property fails.     

Decay of Almost Periodic Solutions¶of Conservation Laws

We consider the asymptotic behavior of solutions of systems of inviscid or viscous conservation laws in one or several space variables, which are almost periodic in the space variables in a generalized sense introduced by Stepanoff and Wiener, which extends the original one of H. Bohr. We prove that if u(x,t) is such a solution whose inclusion intervals at time t, with respect to ?>0, satisfy l _epsiv;(t)/t→0 as t→∞, and such that the scaling sequence u ^T (x,t)=u(T x,T t) is pre-compact as t→∞ in L _loc ¹(? ^{d +1} ₊, then u(x,t) decays to its mean value \(\), which is independent of t, as t→∞. The decay considered here is in L ¹ _loc of the variable ξ≡x/t, which implies, as we show, that \(\) as t→∞, where M _x denotes taking the mean value with respect to x. In many cases we show that, if the initial data are almost periodic in the generalized sense, then so also are the solutions. We also show, in these cases, how to reduce the condition on the growth of the inclusion intervals l _?(t) with t, as t→∞, for fixed ? > 0, to a condition on the growth of l _?(0) with ?, as ?→ 0, which amounts to imposing restrictions only on the initial data. We show with a simple example the existence of almost periodic (non-periodic) functions whose inclusion intervals satisfy any prescribed growth condition as ?→ 0. The applications given here include inviscid and viscous scalar conservation laws in several space variables, some inviscid systems in chromatography and isentropic gas dynamics, as well as many viscous 2 × 2 systems such as those of nonlinear elasticity and Eulerian isentropic gas dynamics, with artificial viscosity, among others. In the case of the inviscid scalar equations and chromatography systems, the class of initial data for which decay results are proved includes, in particular, the L ^∞ generalized limit periodic functions. Our procedures can be easily adapted to provide similar results for semilinear and kinetic relaxations of systems of conservation laws.     

Inviscid Models Generalizing the Two-dimensional Euler and the Surface Quasi-geostrophic Equations

Any classical solution of the two-dimensional incompressible Euler equation is global in time. However, it remains an outstanding open problem whether classical solutions of the surface quasi-geostrophic (SQG) equation preserve their regularity for all time. This paper studies solutions of a family of active scalar equations in which each component u _j of the velocity field u is determined by the scalar θ through \({u_j =\mathcal{R}\Lambda^{-1}P(\Lambda) \theta}\) , where \({\mathcal{R}}\) is a Riesz transform and Λ = (?Δ)^1/2. The two-dimensional Euler vorticity equation corresponds to the special case P(Λ) = I while the SQG equation corresponds to the case P(Λ) = Λ. We develop tools to bound \({\|\nabla u||_{L^\infty}}\) for a general class of operators P and establish the global regularity for the Loglog-Euler equation for which P(Λ) = (log(I + log(I ? Δ))) ^γ with 0 ≦ γ ≦ 1. In addition, a regularity criterion for the model corresponding to P(Λ) = Λ ^β with 0 ≦ β ≦ 1 is also obtained.