On {SL(2, \mathbb R)} Valued Smooth Proximal Cocycles and Cocycles with Positive Lyapunov Exponents Over Irrational Rotation Flows

Consider the class of C ^r -smooth SL(2, \mathbb R){SL(2, \mathbb R)} valued cocycles, based on the rotation flow on the two torus with irrational rotation number α. We show that in this class, (i) cocycles with positive Lyapunov exponents are dense and (ii) cocycles that are either uniformly hyperbolic or proximal are generic, if α satisfies the following Liouville type condition: |a-\fracp_nq_n| ￡ C exp (-q^r+1+k_n)\left|\alpha-\frac{p_n}{q_n}\right| \leq C {\rm exp} (-q^{r+1+\kappa}_{n}), where C >  0 and 0 < k < 1{0 < \kappa <1 } are some constants and \fracP_nq_n{\frac{P_n}{q_n}} is some sequence of irreducible fractions.     

Fully Developed Forced Convection in a Parallel Plate Channel with a Centered Porous Layer

In this study, fully developed heat and fluid flow in a parallel plate channel partially filled with porous layer is analyzed both analytically and numerically. The porous layer is located at the center of the channel and uniform heat flux is applied at the walls. The heat and fluid flow equations for clear fluid and porous regions are separately solved. Continues shear stress and heat flux conditions at the interface are used to determine the interface velocity and temperature. The velocity and temperature profiles in the channel for different values of Darcy number, thermal conductivity ratio, and porous layer thickness are plotted and discussed. The values of Nusselt number and friction factor of a fully clear fluid channel (Nu _cl = 4.12 and fRe _cl = 24) are used to define heat transfer increment ratio (e_th = Nu_p/Nu_cl)({\varepsilon _{\rm th} =Nu_{\rm p}/Nu_{\rm cl})} and pressure drop increment ratio (e_p = fRe_p/fRe_cl )({\varepsilon_{\rm p} =fRe_{\rm p}/fRe_{\rm cl} )} and observe the effects of an inserted porous layer on the increase of heat transfer and pressure drop. The heat transfer and pressure drop increment ratios are used to define an overall performance (e = e_th/e_p)({\varepsilon = \varepsilon_{\rm th}/\varepsilon_{\rm p})} to evaluate overall benefits of an inserted porous layer in a parallel plate channel. The obtained results showed that for a partially porous filled channel, the value of e{\varepsilon} is highly influenced from Darcy number, but it is not affected from thermal conductivity ratio (k _r) when k _r > 2. For a fully porous material filled channel, the value of e{\varepsilon} is considerably affected from thermal conductivity ratio as the porous medium is in contact with the channel walls.     

On the fundamental linear fractional order differential equation

This paper deals with the rational function approximation of the irrational transfer function G(s) = \fracX(s)E(s) = \frac1[(t₀s)^2m + 2z(t₀s)^m + 1]G(s) = \frac{X(s)}{E(s)} = \frac{1}{[(\tau _{0}s)^{2m} + 2\zeta (\tau _{0}s)^{m} + 1]} of the fundamental linear fractional order differential equation (t₀)^2m\fracd^2mx(t)dt^2m + 2z(t₀)^m\fracd^mx(t)dt^m + x(t) = e(t)(\tau_{0})^{2m}\frac{d^{2m}x(t)}{dt^{2m}} + 2\zeta(\tau_{0})^{m}\frac{d^{m}x(t)}{dt^{m}} + x(t) = e(t), for 0<m<1 and 0<ζ<1. An approximation method by a rational function, in a given frequency band, is presented and the impulse and the step responses of this fractional order system are derived. Illustrative examples are also presented to show the exactitude and the usefulness of the approximation method.     

The Critical Dimension for a Fourth Order Elliptic Problem with Singular Nonlinearity

We study the regularity of the extremal solution of the semilinear biharmonic equation ${{\Delta^2} u=\frac{\lambda}{(1-u)^2}}We study the regularity of the extremal solution of the semilinear biharmonic equation D² u=\fracl(1-u)²{{\Delta^2} u=\frac{\lambda}{(1-u)^2}}, which models a simple micro-electromechanical system (MEMS) device on a ball B ì \mathbbR^N{B\subset{\mathbb{R}}^N}, under Dirichlet boundary conditions u=?_n u=0{u=\partial_\nu u=0} on ?B{\partial B}. We complete here the results of Lin and Yang [14] regarding the identification of a “pull-in voltage” λ* > 0 such that a stable classical solution u _λ with 0 < u _λ < 1 exists for l ? (0,l^*){\lambda\in (0,\lambda^*)}, while there is none of any kind when λ > λ*. Our main result asserts that the extremal solution u_l^*{u_{\lambda^*}} is regular (sup_B u_l^* < 1 ){({\rm sup}_B u_{\lambda^*} <1 )} provided N \leqq 8{N \leqq 8} while u_l^*{u_{\lambda^*}} is singular (sup_B u_l^* = 1){({\rm sup}_B u_{\lambda^*} =1)} for N \geqq 9{N \geqq 9}, in which case 1-C₀|x|^4/3 \leqq u_l^* (x) \leqq 1-|x|^4/3{1-C_0|x|^{4/3} \leqq u_{\lambda^*} (x) \leqq 1-|x|^{4/3}} on the unit ball, where C₀:=(\fracl^*[`(l)])<sup>\frac</sup>13{C_0:=\left(\frac{\lambda^*}{\overline{\lambda}}\right)^\frac{1}{3}} and [`(l)]: = \frac89(N-\frac23)(N- \frac83){\bar{\lambda}:= \frac{8}{9}\left(N-\frac{2}{3}\right)\left(N- \frac{8}{3}\right)}.     

Metastable Periodic Patterns in Singularly Perturbed Delayed Equations

We consider the scalar delayed differential equation e[(x)\dot](t)=-x(t) +f(x(t-1)){\epsilon\dot x(t)=-x(t)\,+f(x(t-1))}, where ${\epsilon\,{>}\,0}${\epsilon\,{>}\,0} and f verifies either df/dx > 0 or df/dx < 0 and some other conditions. We present theorems indicating that a generic initial condition with sign changes generates a solution with a transient time of order exp(c/e){{\rm exp}(c{/}\epsilon)}, for some c > 0. We call it a metastable solution. During this transient a finite time span of the solution looks like that of a periodic function. It is remarkable that if df/dx > 0 then f must be odd or present some other very special symmetry in order to support metastable solutions, while this condition is absent in the case df/dx < 0. Explicit e{\epsilon}-asymptotics for the motion of zeroes of a solution and for the transient time regime are presented.     

Effect of thermal losses on the microscopic hyperbolic heat conduction model

The effects of radiative losses on the thermal behavior of thin metal films, as described by the microscopic two-step hyperbolic heat conduction model, are investigated. Different criteria, which determine the ranges within which thermal radiative losses are significant, are derived. It is found that radiative losses from the electron gas are significant in thin films having [(C_R e_e^4/₃ T_￥ ⁴ )/(k_e^1/₃ L^2/₃ G)] 3 4.6 ×10⁷{{C_R \epsilon _e^{{4 \over 3}} T_\infty ^4 } \over {k_e^{{1 \over 3}} L^{{2 \over 3}} G}}\geq 4.6 \times 10^7 for /_o > 4 and F_F < 1 and [(C_R e_e^3/₂ T_￥ ^9/₂)/(k_e^1/₂ L^1/₂ G)] 3 7.4 ×10¹⁰{{C_R \epsilon _e^{{3 \over 2}} T_\infty ^{{9 \over 2}}} \over {k_e^{{1 \over 2}} L^{{1 \over 2}} G}}\geq 7.4 \times 10^{10} for /_o < 4 and F_F > 1.     

On the Dynamics of Solitons in the Nonlinear Schr?dinger Equation

We study the behavior of the soliton solutions of the equation i\frac?y?t = - \frac12m Dy+ \frac12W_e^￠(y) + V(x)y,i\frac{\partial\psi}{{\partial}t} = - \frac{1}{2m} \Delta\psi + \frac{1}{2}W_{\varepsilon}^{\prime}(\psi) + V(x){\psi},     

Hydrodynamic Limit of Gradient Exclusion Processes with Conductances

Fix a strictly increasing right continuous with left limits function ${W: \mathbb{R} \to \mathbb{R}}Fix a strictly increasing right continuous with left limits function W: \mathbbR ? \mathbbR{W: \mathbb{R} \to \mathbb{R}} and a smooth function F: [l,r] ? \mathbb R{\Phi : [l,r] \to \mathbb R}, defined on some interval [l, r] of \mathbb R{\mathbb R}, such that 0 < b\leqq F￠\leqq b^-1{0 < b\leqq \Phi'\leqq b^{-1}}. On the diffusive time scale, the evolution of the empirical density of exclusion processes with conductances given by W is described by the unique weak solution of the non-linear differential equation ?_t r = (d/dx)(d/dW) F(r){\partial_t \rho = ({\rm d}/{\rm d}x)({\rm d}/{\rm d}W) \Phi(\rho)}. We also present some properties of the operator (d/dx)(d/dW).     

On the Existence and Conditional Energetic Stability of Solitary Gravity-Capillary Surface Waves on Deep Water

This paper presents an existence and stability theory for gravity-capillary solitary waves on the surface of a body of water of infinite depth. Exploiting a classical variational principle, we prove the existence of a minimiser of the wave energy E{{\mathcal E}} subject to the constraint I=?2m{{\mathcal I}=\sqrt{2}\mu}, where I{{\mathcal I}} is the wave momentum and 0 < m << 1{0 < \mu \ll 1} . Since E{{\mathcal E}} and I{{\mathcal I}} are both conserved quantities a standard argument asserts the stability of the set D _μ of minimisers: solutions starting near D _μ remain close to D _μ in a suitably defined energy space over their interval of existence. In the applied mathematics literature solitary water waves of the present kind are modelled as solutions of the nonlinear Schr?dinger equation with cubic focussing nonlinearity. We show that the waves detected by our variational method converge (after an appropriate rescaling) to solutions of this model equation as mˉ 0{\mu \downarrow 0} .     

Regional Blow-up for a Doubly Nonlinear Parabolic Equation with a Nonlinear Boundary Condition

In this work, positive solutions to a doubly nonlinear parabolic equation with a nonlinear boundary condition are considered. We study the problem where 0 < m, r, α < ∞ are parameters. It is known that for some values of the parameters there are solutions that blow up in finite time. We determine in terms of α ,m, r the blow-up sets for these solutions. We prove that single point blow-up occurs if max{m, r} < α, global blow-up appears for the range of parameters 0 < m < α < r and regional blow-up takes place if 0 < m < α = r and . In this case the blow-up set consists of the interval .     

A scaling law near the primary resonance of the quasiperiodic Mathieu equation

We analyze the quasiperiodic damped Mathieu equation

Positive Least Energy Solutions and Phase Separation for Coupled Schr?dinger Equations with Critical Exponent

In this paper we study the following coupled Schr?dinger system, which can be seen as a critically coupled perturbed Brezis–Nirenberg problem: {ll-Du +l₁ u = m₁ u³+buv²,     x ? W,-Dv +l₂ v = m₂ v³+bvu²,     x ? W,u\geqq 0, v\geqq 0 in W,    u=v=0     on ?W.\left\{\begin{array}{ll}-\Delta u +\lambda_1 u = \mu_1 u^3+\beta uv^2, \quad x\in \Omega,\\-\Delta v +\lambda_2 v =\mu_2 v^3+\beta vu^2, \quad x\in \Omega,\\u\geqq 0, v\geqq 0\, {\rm in}\, \Omega,\quad u=v=0 \quad {\rm on}\, \partial\Omega.\end{array}\right.     

Well-posedness and Regularity for the Elasticity Equation with Mixed Boundary Conditions on Polyhedral Domains and Domains with Cracks

We prove a regularity result for the anisotropic linear elasticity equation ${P u := {\rm div} \left( \boldmath\mathsf{C} \cdot \nabla u\right) = f}We prove a regularity result for the anisotropic linear elasticity equationP u : = div ( C ·?u) = f{P u := {\rm div} \left( \boldmath\mathsf{C} \cdot \nabla u\right) = f} , with mixed (displacement and traction) boundary conditions on a curved polyhedral domain W ì \mathbbR³{\Omega \subset \mathbb{R}^3} in weighted Sobolev spaces K^m+1_a+1(W){\mathcal {K}^{m+1}_{a+1}(\Omega)} , for which the weight is given by the distance to the set of edges. In particular, we show that there is no loss of K^m_a{\mathcal {K}^{m}_{a}} -regularity. Our curved polyhedral domains are allowed to have cracks. We establish a well-posedness result when there are no neighboring traction boundary conditions and |a| < η, for some small η > 0 that depends on P, on the boundary conditions, and on the domain Ω. Our results extend to other strongly elliptic systems and higher dimensions.     

Global Nonexistence Theorems for a Class of Evolution Equations with Dissipation

We study abstract evolution equations with nonlinear damping terms and source terms, including as a particular case a nonlinear wave equation of the type $ \ba{cl} u_{tt}-\Delta u+ b|u_t|^{m-2}u_t=c|u|^{p-2}u, &;(t,x)\in [0,T)\times\Omega,\\[6pt] u(t,x)=0, &;(t,x)\in [0,T)\times\partial \Omega,\\[6pt] u(0,\cdot)=u_0\in H_0^1(\Omega), \quad u_t(0,\cdot)=v_0\in L^2(\Omega),\es&; \ea $ \ba{cl} u_{tt}-\Delta u+ b|u_t|^{m-2}u_t=c|u|^{p-2}u, &;(t,x)\in [0,T)\times\Omega,\\[6pt] u(t,x)=0, &;(t,x)\in [0,T)\times\partial \Omega,\\[6pt] u(0,\cdot)=u_0\in H_0^1(\Omega), \quad u_t(0,\cdot)=v_0\in L^2(\Omega),\es&; \ea where 0 < T ￡ ￥0\Omega is a bounded regular open subset of \mathbbRⁿ\mathbb{R}^n, n 3 1n\ge 1, b,c > 0b,c>0, p > 2p>2, m > 1m>1. We prove a global nonexistence theorem for positive initial value of the energy when 1 < m < p,    2 < p ￡ \frac2nn-2. 1-Laplacian operator, q > 1q>1.     

Isolated Singularity for the Stationary Navier—Stokes System

In this paper the classical method to prove a removable singularity theorem for harmonic functions near an isolated singular point is extended to solutions to the stationary Stokes and Navier—Stokes system. Finding series expansion of solutions in terms of homogeneous harmonic polynomials, we establish some known results and new theorems concerning the behavior of solutions near an isolated singular point. In particular, we prove that if (u, p) is a solution to the Navier—Stokes system in B_R \{0} B_R \setminus \{0\} , n 3 3 n \geq 3 and |u(x)| = o (|x|^{-(n - 1)/2}) |u(x)| = o\,(|x|^{-(n - 1)/2}) as |x| ? 0 |x| \to 0 or u ? L^{2n/(n - 1)}(B_R) u \in L^{2n/(n - 1)}(B_R) , then (u, p) is a distribution solution and if in addition, u ? L^b(B_R) u \in L^{\beta}(B_R) for some b > n \beta > n then ( u, p) is smooth in BR.     

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.      17. Some Aspects of the Asymptotic Dynamics of Solutions of the Homogeneous Navier–Stokes Equations in General Domains      Zdeněk Skalák 《Journal of Mathematical Fluid Mechanics》2010,12(4):503-535 In the first part of the paper we study decays of solutions of the Navier–Stokes equations on short time intervals. We show, for example, that if w is a global strong nonzero solution of homogeneous Navier–Stokes equations in a sufficiently smooth (unbounded) domain Ω ⊆ R3 and β ∈[1/2, 1) , then there exist C0 > 1 and δ0 ∈ (0, 1) such that \frac |||w(t)|||b|||w(t + d)|||b ￡ C0{\frac {|||w(t)|||_\beta}{|||w(t + \delta)|||_{\beta}}} \leq C_0      18. The Erpenbeck High Frequency Instability Theorem for Zeldovitch–von Neumann–D?ring Detonations      Olivier Lafitte  Mark Williams  Kevin Zumbrun 《Archive for Rational Mechanics and Analysis》2012,204(1):141-187 The rigorous study of spectral stability for strong detonations was begun by Erpenbeck (Phys. Fluids 5:604–614 1962). Working with the Zeldovitch–von Neumann–D?ring (ZND) model (more precisely, Erpenbeck worked with an extension of ZND to general chemistry and thermodynamics), which assumes a finite reaction rate but ignores effects such as viscosity corresponding to second order derivatives, he used a normal mode analysis to define a stability function V(t,e){V(\tau,\epsilon)} whose zeros in ${\mathfrak{R}\tau > 0}${\mathfrak{R}\tau > 0} correspond to multidimensional perturbations of a steady detonation profile that grow exponentially in time. Later in a remarkable paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966; Stability of detonations for disturbances of small transverse wavelength, 1965) he provided strong evidence, by a combination of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of V exist for perturbations of sufficiently large transverse wavenumber e{\epsilon} , even when the von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly 20 years later) by Majda. In spite of a great deal of later numerical work devoted to computing the zeros of V(t,e){V(\tau,\epsilon)} , the paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966) remains one of the few works we know of [another is Erpenbeck (Phys. Fluids 7:684–696, 1964), which considers perturbations for which the ratio of longitudinal over transverse components approaches ∞] that presents a detailed and convincing theoretical argument for detecting them. The analysis in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) points the way toward, but does not constitute, a mathematical proof that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) and provide proofs, together with certain simplifications and extensions, of the main conclusions about stability and instability of detonations contained in that paper. The main mathematical problem, and our principal focus here, is to determine the precise asymptotic behavior as e?￥{\epsilon\to\infty} of solutions to a linear system of ODEs in x, depending on e{\epsilon} and a complex frequency τ as parameters, with turning points x * on the half-line [0,∞).      19. Spectral Stability of Small-Amplitude Viscous Shock Waves in Several Space Dimensions      Heinrich Freistühler  Peter Szmolyan 《Archive for Rational Mechanics and Analysis》2010,195(2):353-373 A planar viscous shock profile of a hyperbolic–parabolic system of conservation laws is a steady solution in a moving coordinate frame. The asymptotic stability of viscous profiles and the related vanishing-viscosity limit are delicate questions already in the well understood case of one space dimension and even more so in the case of several space dimensions. It is a natural idea to study the stability of viscous profiles by analyzing the spectrum of the linearization about the profile. The Evans function method provides a geometric dynamical-systems framework to study the eigenvalue problem. In this approach eigenvalues correspond to zeros of an essentially analytic function E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} which detects nontrivial intersections of the so-called stable and unstable spaces, that is, spaces of solutions that decay on one (“−∞”) or the other side (“ + ∞”) of the shock wave, respectively. In a series of pioneering papers, Kevin Zumbrun and collaborators have established in various contexts that spectral stability, that is, the non-vanishing of E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} and the non-vanishing of the Lopatinski–Kreiss–Majda function Δ(λ,ω), imply nonlinear stability of viscous shock profiles in several space dimensions. In this paper we show that these conditions hold true for small amplitude extreme shocks under natural assumptions. This is done by exploiting the slow-fast nature of the small-amplitude limit, which was used in a previous paper by the authors to prove spectral stability of small-amplitude shock waves in one space dimension. Geometric singular perturbation methods are applied to decompose the stable and unstable spaces into subbundles with good control over their limiting behavior. Three qualitatively different regimes are distinguished that relate the small strength e{\epsilon} of the shock wave to appropriate ranges of values of the spectral parameters (ρλ, ρ ω). Various rescalings are used to overcome apparent degeneracies in the problem caused by loss of hyperbolicity or lack of transversality.      20. Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall      M. Pavlidis  Y. Dimakopoulos  J. Tsamopoulos 《Rheologica Acta》2009,48(9):1031-1048 The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]s/([(h)\tilde]s +[(h)\tilde]p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]p\widetilde{\eta}_p and [(h)\tilde]s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]2/([(h)\tilde]p +[(h)\tilde]s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号

Some Aspects of the Asymptotic Dynamics of Solutions of the Homogeneous Navier–Stokes Equations in General Domains

In the first part of the paper we study decays of solutions of the Navier–Stokes equations on short time intervals. We show, for example, that if w is a global strong nonzero solution of homogeneous Navier–Stokes equations in a sufficiently smooth (unbounded) domain Ω ⊆ R³ and β ∈[1/2, 1) , then there exist C₀ > 1 and δ₀ ∈ (0, 1) such that

The Erpenbeck High Frequency Instability Theorem for Zeldovitch–von Neumann–D?ring Detonations

The rigorous study of spectral stability for strong detonations was begun by Erpenbeck (Phys. Fluids 5:604–614 1962). Working with the Zeldovitch–von Neumann–D?ring (ZND) model (more precisely, Erpenbeck worked with an extension of ZND to general chemistry and thermodynamics), which assumes a finite reaction rate but ignores effects such as viscosity corresponding to second order derivatives, he used a normal mode analysis to define a stability function V(t,e){V(\tau,\epsilon)} whose zeros in ${\mathfrak{R}\tau > 0}${\mathfrak{R}\tau > 0} correspond to multidimensional perturbations of a steady detonation profile that grow exponentially in time. Later in a remarkable paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966; Stability of detonations for disturbances of small transverse wavelength, 1965) he provided strong evidence, by a combination of formal and rigorous arguments, that for certain classes of steady ZND profiles, unstable zeros of V exist for perturbations of sufficiently large transverse wavenumber e{\epsilon} , even when the von Neumann shock, regarded as a gas dynamical shock, is uniformly stable in the sense defined (nearly 20 years later) by Majda. In spite of a great deal of later numerical work devoted to computing the zeros of V(t,e){V(\tau,\epsilon)} , the paper (Erpenbeck in Phys. Fluids 9:1293–1306, 1966) remains one of the few works we know of [another is Erpenbeck (Phys. Fluids 7:684–696, 1964), which considers perturbations for which the ratio of longitudinal over transverse components approaches ∞] that presents a detailed and convincing theoretical argument for detecting them. The analysis in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) points the way toward, but does not constitute, a mathematical proof that such unstable zeros exist. In this paper we identify the mathematical issues left unresolved in Erpenbeck (Phys. Fluids 9:1293–1306, 1966) and provide proofs, together with certain simplifications and extensions, of the main conclusions about stability and instability of detonations contained in that paper. The main mathematical problem, and our principal focus here, is to determine the precise asymptotic behavior as e?￥{\epsilon\to\infty} of solutions to a linear system of ODEs in x, depending on e{\epsilon} and a complex frequency τ as parameters, with turning points x _* on the half-line [0,∞).     

Spectral Stability of Small-Amplitude Viscous Shock Waves in Several Space Dimensions

A planar viscous shock profile of a hyperbolic–parabolic system of conservation laws is a steady solution in a moving coordinate frame. The asymptotic stability of viscous profiles and the related vanishing-viscosity limit are delicate questions already in the well understood case of one space dimension and even more so in the case of several space dimensions. It is a natural idea to study the stability of viscous profiles by analyzing the spectrum of the linearization about the profile. The Evans function method provides a geometric dynamical-systems framework to study the eigenvalue problem. In this approach eigenvalues correspond to zeros of an essentially analytic function E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} which detects nontrivial intersections of the so-called stable and unstable spaces, that is, spaces of solutions that decay on one (“−∞”) or the other side (“ + ∞”) of the shock wave, respectively. In a series of pioneering papers, Kevin Zumbrun and collaborators have established in various contexts that spectral stability, that is, the non-vanishing of E(rl,rw){\mathcal{E}(\rho\lambda,\rho\omega)} and the non-vanishing of the Lopatinski–Kreiss–Majda function Δ(λ,ω), imply nonlinear stability of viscous shock profiles in several space dimensions. In this paper we show that these conditions hold true for small amplitude extreme shocks under natural assumptions. This is done by exploiting the slow-fast nature of the small-amplitude limit, which was used in a previous paper by the authors to prove spectral stability of small-amplitude shock waves in one space dimension. Geometric singular perturbation methods are applied to decompose the stable and unstable spaces into subbundles with good control over their limiting behavior. Three qualitatively different regimes are distinguished that relate the small strength e{\epsilon} of the shock wave to appropriate ranges of values of the spectral parameters (ρλ, ρ ω). Various rescalings are used to overcome apparent degeneracies in the problem caused by loss of hyperbolicity or lack of transversality.     

Fully developed flow of a viscoelastic film down a vertical cylindrical or planar wall

The one-dimensional, gravity-driven film flow of a linear (l) or exponential (e) Phan-Thien and Tanner (PTT) liquid, flowing either on the outer or on the inner surface of a vertical cylinder or over a planar wall, is analyzed. Numerical solution of the governing equations is generally possible. Analytical solutions are derived only for: (1) l-PTT model in cylindrical and planar geometries in the absence of solvent, b o [(h)\tilde]_s/([(h)\tilde]_s +[(h)\tilde]_p)=0\beta\equiv {\tilde{\eta}_s}/\left({\tilde{\eta}_s +\tilde{\eta}_p}\right)=0, where [(h)\tilde]_p\widetilde{\eta}_p and [(h)\tilde]_s\widetilde{\eta}_s are the zero-shear polymer and solvent viscosities, respectively, and the affinity parameter set at ξ = 0; (2) l-PTT or e-PTT model in a planar geometry when β = 0 and x 1 0\xi \ne 0; (3) e-PTT model in planar geometry when β = 0 and ξ = 0. The effect of fluid properties, cylinder radius, [(R)\tilde]\tilde{R}, and flow rate on the velocity profile, the stress components, and the film thickness, [(H)\tilde]\tilde{H}, is determined. On the other hand, the relevant dimensionless numbers, which are the Deborah, De=[(l)\tilde][(U)\tilde]/[(H)\tilde]De={\tilde{\lambda}\tilde{U}}/{\tilde{H}}, and Stokes, St=[(r)\tilde][(g)\tilde][(H)\tilde]²/([(h)\tilde]_p +[(h)\tilde]_s )[(U)\tilde]St=\tilde{\rho}\tilde{g}\tilde{\rm H}^{2}/\left({\tilde{\eta}_p +\tilde{\eta}_s} \right)\tilde{U}, numbers, depend on [(H)\tilde]\tilde{H} and the average film velocity, [(U)\tilde]\widetilde{U}. This makes necessary a trial and error procedure to obtain [(H)\tilde]\tilde{H} a posteriori. We find that increasing De, ξ, or the extensibility parameter ε increases shear thinning resulting in a smaller St. The Stokes number decreases as [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to zero for a film on the outer cylindrical surface, while it asymptotes to very large values when [(R)\tilde]/[(H)\tilde]{\tilde{R}}/{\tilde{H}} decreases down to unity for a film on the inner surface. When x 1 0\xi \ne 0, an upper limit in De exists above which a solution cannot be computed. This critical value increases with ε and decreases with ξ.