The Damped-Driven 2D Navier–Stokes System on Large Elongated Domains

The Navier–Stokes system with damping, which is motivated by Stommel–Charney model of ocean circulation, is considered in a large elongated periodic rectangular domain with area of the order α⁻¹, as α → 0. We obtain estimates for the dimension of the global attractor that are sharp as both α → 0 and ν → 0, where ν is the viscosity coefficient. This work was supported in part by the US Civilian Research and Development Foundation, grant no. RUM1-2654-MO-05 (A.A.I. and E.S.T.). The work of A.A.I. was supported in part by the Russian Foundation for Fundamental Research, grants no. 06-001-0096 and no. 05-01-429, and by the RAS Programme no. 1 ‘Modern problems of theoretical mathematics’. The work of E.S.T. was supported in part by the NSF, grant no. DMS-0204794, the MAOF Fellowship of the Israeli Council of Higher Education, and by the BSF, grant no. 200423.     

Vorticities in a LES Model for 3D Periodic Turbulent Flows

This paper is devoted to the study of a LES model to simulate turbulent 3D periodic flow. We focus our attention on the vorticity equation derived from this LES model for small values of the numerical grid size δ. We obtain entropy inequalities for the sequence of corresponding vorticities and corresponding pressures independent of δ, provided the initial velocity u₀ is in L_x² while the initial vorticity ω₀ = ∇ × u₀ is in L_x¹. When δ tends to zero, we show convergence, in a distributional sense, of the corresponding equations for the vorticities to the classical 3D equation for the vorticity.     

Stability of Discrete Stokes Operators in Fractional Sobolev Spaces

Using a general approximation setting having the generic properties of finite-elements, we prove uniform boundedness and stability estimates on the discrete Stokes operator in Sobolev spaces with fractional exponents. As an application, we construct approximations for the time-dependent Stokes equations with a source term in L ^p (0, T; L ^q (Ω)) and prove uniform estimates on the time derivative and discrete Laplacian of the discrete velocity that are similar to those in Sohr and von Wahl [20]. On long leave from LIMSI (CNRS-UPR 3251), BP 133, 91403, Orsay, France.     

Temporal and Spatial Decays for the Stokes Flow

We estimate the time decay rates in L ¹, in the Hardy space and in L ^∞ of the gradient of solutions for the Stokes equations on the half spaces. For the estimates in the Hardy space we adopt the ideas in [7], and also use the heat kernel and the solution formula for the Stokes equations. We also estimate the temporal-spatial asymptotic estimates in L ^q , 1 < q < ∞, for the Stokes solutions. This work was supported by grant No. (R05-2002-000-00002-0(2002)) from the Basic Research Program of the Korea Science & Engineering Foundation.     

Regularity of Flow of Anisotropic Fluid

We study the steady flow of an anisotropic generalised Newtonian fluid under Dirichlet boundary conditions in a bounded domain . The fluid is characterised by a nonlinear dependence of the stress tensor on the symmetric gradient of the velocity vector field. We prove the existence of a C ^1,α-solution of this problem under certain assumptions on the growth of the elliptic term. The result is global: we prove the regularity up to the boundary of the domain Ω. This research was supported by the grant GACR 201/03/0934; partially also by the grants MSM 0021620839 and GAUK 262/2002/B-MAT/MFF.     

Relaxation patterns of long,linear, flexible,monodisperse polymers: BSW spectrum revisited

Theoretical predictions for the dynamic moduli of long, linear, flexible, monodisperse polymers are summarized and compared with experimental observations. Surprisingly, the predicted 1/2 power scaling of the long-time modes of the relaxation spectrum is not found in the experiments. Instead, scaling with a power of about 1/4 extends all the way up to the longest relaxation times near τ/τ _max = 1. This is expressed in the empirical relaxation time spectrum of Baumgaertel-Schausberger-Winter, denoted as “BSW spectrum,” and justifies a closer look at the properties of the BSW spectrum. Working with the BSW spectrum, however, is made difficult by the fact that hypergeometric functions occur naturally in BSW-based rheological material functions. BSW provides no explicit solutions for the dynamic moduli, G ^′(ω), G ^″(ω), or the relaxation modulus G(t). To overcome this problem, close approximations of simple analytical form are shown for these moduli. With these approximations, analysis of linear viscoelastic data allows the direct determination of BSW parameters.     

Transonic Shock Problem for the Euler System in a Nozzle

In this paper, we study the well-posedness problem on transonic shocks for steady ideal compressible flows through a two-dimensional slowly varying nozzle with an appropriately given pressure at the exit of the nozzle. This is motivated by the following transonic phenomena in a de Laval nozzle. Given an appropriately large receiver pressure P _r , if the upstream flow remains supersonic behind the throat of the nozzle, then at a certain place in the diverging part of the nozzle, a shock front intervenes and the flow is compressed and slowed down to subsonic speed, and the position and the strength of the shock front are automatically adjusted so that the end pressure at exit becomes P _r , as clearly stated by Courant and Friedrichs [Supersonic flow and shock waves, Interscience Publishers, New York, 1948 (see section 143 and 147)]. The transonic shock front is a free boundary dividing two regions of C ^2,α flow in the nozzle. The full Euler system is hyperbolic upstream where the flow is supersonic, and coupled hyperbolic-elliptic in the downstream region Ω₊ of the nozzle where the flow is subsonic. Based on Bernoulli’s law, we can reformulate the problem by decomposing the 3 × 3 Euler system into a weakly coupled second order elliptic equation for the density ρ with mixed boundary conditions, a 2 × 2 first order system on u ₂ with a value given at a point, and an algebraic equation on (ρ, u ₁, u ₂) along a streamline. In terms of this reformulation, we can show the uniqueness of such a transonic shock solution if it exists and the shock front goes through a fixed point. Furthermore, we prove that there is no such transonic shock solution for a class of nozzles with some large pressure given at the exit. This research was supported in part by the Zheng Ge Ru Foundation when Yin Huicheng was visiting The Institute of Mathematical Sciences, The Chinese University of Hong Kong. Xin is supported in part by Hong Kong RGC Earmarked Research Grants CUHK-4028/04P, CUHK-4040/06P, and Central Allocation Grant CA05-06.SC01. Yin is supported in part by NNSF of China and Doctoral Program of NEM of China.     

Flux Norm Approach to Finite Dimensional Homogenization Approximations with Non-Separated Scales and High Contrast

We consider linear divergence-form scalar elliptic equations and vectorial equations for elasticity with rough (L ^∞(Ω), W ì \mathbb R^d{\Omega \subset \mathbb R^d}) coefficients a(x) that, in particular, model media with non-separated scales and high contrast in material properties. While the homogenization of PDEs with periodic or ergodic coefficients and well separated scales is now well understood, we consider here the most general case of arbitrary bounded coefficients. For such problems, we introduce explicit and optimal finite dimensional approximations of solutions that can be viewed as a theoretical Galerkin method with controlled error estimates, analogous to classical homogenization approximations. In particular, this approach allows one to analyze a given medium directly without introducing the mathematical concept of an e{\epsilon} family of media as in classical homogenization. We define the flux norm as the L ² norm of the potential part of the fluxes of solutions, which is equivalent to the usual H ¹-norm. We show that in the flux norm, the error associated with approximating, in a properly defined finite-dimensional space, the set of solutions of the aforementioned PDEs with rough coefficients is equal to the error associated with approximating the set of solutions of the same type of PDEs with smooth coefficients in a standard space (for example, piecewise polynomial). We refer to this property as the transfer property. A simple application of this property is the construction of finite dimensional approximation spaces with errors independent of the regularity and contrast of the coefficients and with optimal and explicit convergence rates. This transfer property also provides an alternative to the global harmonic change of coordinates for the homogenization of elliptic operators that can be extended to elasticity equations. The proofs of these homogenization results are based on a new class of elliptic inequalities. These inequalities play the same role in our approach as the div-curl lemma in classical homogenization.     

Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates

In this paper we consider the numerical approximation of steady and unsteady generalized Newtonian fluid flows using divergence free finite elements generated by the Powell–Sabin–Heindl elements. We derive a priori and a posteriori finite element error estimates and prove convergence of the method of successive approximations for the steady flow case. A priori error estimates of unsteady flows are also considered. These results provide a theoretical foundation and supporting numerical studies are to be provided in Part II. Copyright © 2003 John Wiley & Sons, Ltd.     

Existence of a Solution “in the Large” for Ocean Dynamics Equations

For the system of equations describing the large-scale ocean dynamics, an existence and uniqueness theorem is proved “in the large”. This system is obtained from the 3D Navier–Stokes equations by changing the equation for the vertical velocity component u ₃ under the assumption of smallness of a domain in z-direction, and a nonlinear equation for the density function ρ is added. More precisely, it is proved that for an arbitrary time interval [0, T], any viscosity coefficients and any initial conditions

The effect of non-steady wall shear on pressure surges in conduits carrying viscous liquids

Summary A study is made of the attenuation of pressure surges in a two-dimension a channel carrying a viscous liquid when a valve at the downstream end is suddenly closed. The analysis differs from previous work in this area by taking into account the transient nature of the wall shear, which in the past has been assumed as equivalent to that existing in steady flow. For large values of the frictional resistance parameter the transient wall shear analysis results in less attenuation than given by the steady wall shear assumption.Nomenclature c /, ft/sec - e base of natural logarithms - F(x, y) integration function, equation (38) - (x) mean value of F(x, y) - g local acceleration of gravity, ft/sec² - h width of conduit, ft - k (2m–1)^{2 2} L/h ² c, equation (35) - k* 12L/h ² c, frictional resistance parameter, equation (46) - L length of conduit, ft - m positive integer - n positive integer - p pressure, lb/ft² - p ₀ constant pressure at inlet of conduit, lb/ft² - P pressure plus elevation head, p+gz, equation (4) - mean value of P over the conduit width h - P ₀ p ₀+gz ₀, lbs/ft² - R frictional resistance coefficient for steady state wall shear, lb sec/ft⁴ - s positive integer; also, condensation, equation (6) - t time, sec - t ct/L, dimensionless time - u x component of fluid velocity, ft/sec - u _m mean velocity in conduit, equation (12), ft/sec - u ₀(y) velocity profile in Poiseuille flow, equation (19), ft/sec - transformed velocity - U initial mean velocity in conduit - x distance along conduit, measured from valve (fig. 1), ft - x x/L, dimensionless distance - y distance normal to conduit wall (fig. 1), ft - y y/h, equation (25) - z elevation, measured from arbitrary datum, ft - z ₀ elevation of constant pressure source, ft - isothermal bulk compression modulus, lbs/ft² - _n , equation (37) - _n (2n–1)/2, equation (36) - viscosity, slugs/ft sec - / = kinematic viscosity, ft²/sec - density of fluid, slugs/ft³ - ₀ density of undisturbed fluid, slugs/ft³ - ø angle between conduit and vertical (fig. 1) The research upon which this paper is based was supported by a grant from the National Science Foundation.     

THE SOLUTION TO THE DESTABILIZING CRITICAL LOAD OF CIRCULAR DOUBLE ARTICULATED ARCH UNDER GOING VERTICAL DISTRIBUTIVE LOAD g0/cos2θ

In this paper, after taking the effect of axis force on bending into consideration, the general potential energy for the circular double articulated arch is established undergoing vertical distributive load g₀/cos²θ. With sufficient engineering precision, the fourth approximations to the destabilizing critical load of the arch under t his load are obtained by Ritz method. The approximations to the critical load ta ble are listed for various center angles of arch, and are contrasted with the critical load circular arch undergoing radial uniform load. Some reference results have been obtained.     

Global stability of a virus dynamics model with Beddington–DeAngelis incidence rate and CTL immune response

In this paper, the global stability of virus dynamics model with Beddington–DeAngelis infection rate and CTL immune response is studied by constructing Lyapunov functions. We derive the basic reproduction number R ₀ and the immune response reproduction number R ⁰ for the virus infection model, and establish that the global dynamics are completely determined by the values of R ₀. We obtain the global stabilities of the disease-free equilibrium E ₀, immune-free equilibrium E ₁ and endemic equilibrium E ^∗ when R ₀≤1, R ₀>1, R ⁰>1, respectively.     

Regularity of the Inverse of a Planar Sobolev Homeomorphism

Let be a domain. Suppose that f ∈ W^1,1_loc(Ω,R²) is a homeomorphism such that Df(x) vanishes almost everywhere in the zero set of J_f. We show that f^-1 ∈ W^1,1_loc(f(Ω),R²) and that Df⁻¹(y) vanishes almost everywhere in the zero set of Sharp conditions to quarantee that f⁻¹ ∈ W^1,q(f(Ω),R²) for some 1<q≤2 are also given.     

The Dirichlet Problem for H-Systems with Small Boundary Data: BlowUp Phenomena and Nonexistence Results

Given H:ℝ³→ℝ of class C¹ and bounded, we consider a sequence (uⁿ) of solutions of the H-system in the unit open disc satisfying the boundary condition uⁿ=γⁿ on ∂. In the first part of this paper, assuming that (uⁿ) is bounded in H¹(,ℝ³) we study the behavior of (uⁿ) when the boundary data γⁿ shrink to zero. We show that either uⁿ→0 strongly in H¹(,ℝ³) or uⁿ blows up at least one H-bubble ω, namely a nonconstant, conformal solution of the H-system on ℝ². Under additional assumptions on H, we can obtain more precise information on the blow up. In the second part of this paper we investigate the multiplicity of solutions for the Dirichlet problem on the disc with small boundary datum. We detect a family of nonconstant functions H (even close to a nonzero constant in any reasonable topology) for which the Dirichlet problem cannot admit a ``large' solution at a mountain pass level when the boundary datum is small.     

Higher approximations to the homogeneous solutions for toroid

To analyze the strength of the elastic, isotropic toroid with axisymmetrical loads, it is necessary to find the solutions for a nonhomogeneous second-order ordinary-differential equation in complex form. When the parameter = 12(1-v²) a² / (r ₀ h), which occurs in the equation, is large, the asymptotic approach is often used. Up to now, only the first approximations to its homogeneous solutions are found. In this paper the author has obtained higher approximations to the homogeneous solutions, thus, reaching the precision of linear thin shell theory.This work was supported by the National Natural Science Foundation of China.     

Stability, Instability, and Error of the Force-based Quasicontinuum Approximation

Due to their algorithmic simplicity and high accuracy, force-based model coupling techniques are popular tools in computational physics. For example, the force-based quasicontinuum (QCF) approximation is the only known pointwise consistent quasicontinuum approximation for coupling a general atomistic model with a finite element continuum model. In this paper, we present a detailed stability and error analysis of this method. Our optimal order error estimates provide a theoretical justification for the high accuracy of the QCF approximation: they clearly demonstrate that the computational efficiency of continuum modeling can be utilized without a significant loss of accuracy if defects are captured in the atomistic region. The main challenge we need to overcome is the fact that the linearized QCF operator is typically not positive definite. Moreover, we prove that no uniform inf-sup stability condition holds for discrete versions of the W ^1,p -W ^1,q “duality pairing” with 1/p + 1/q = 1, if 1 ≤ p < ∞. However, we were able to establish an inf-sup stability condition for a discrete version of the W ^1,∞-W ^1,1 “duality pairing” which leads to optimal order error estimates in a discrete W ^1,∞-norm.     

Drag of a strongly heated sphere at small Reynolds numbers

Using an asymptotic small-perturbation method, the flow around a strongly heated sphere at small Reynolds numbers Re ≪ 1 is considered with account for thermal stresses in the gas in the higher-order approximations, beyond the Stokes one. It is assumed that the value of the Prandtl number Pr is arbitrary and the temperature dependence of the viscosity is described by a power law with an arbitrary exponent. In the O(Re²) and O(Re³ ln(Re)) approximations, the drag force of a heated sphere is found over a wide range of the ratios of sphere’s temperature to the gas free-stream temperature T _W /T _∞. The limits of applicability of the first (in Re) approximation are investigated, including the negative-drag effect, attributable to the action of the thermal stresses. The results are compared with numerical calculations of the flow around a hot sphere. The limits of applicability of the approximations found are examined. Similar results are obtained for the standard Navier-Stokes equations in which the thermal stresses are neglected.     

Nonlinear vibrations of rectangular plates with linearly varying thickness

Simplified nonlinear governing differential equations proposed by Berger for static cases and extended by Nash and Modeer for dynamic cases are used to analyse the title problem. Steady-state harmonic oscillations are assumed and the time variable is eliminated by a Kantorovich averaging method. The enclosure or comparison theorem of Collatz is then applied to the reduced equations to obtain the upper and lower bounds for the fundamental nonlinear frequency of simply-supported rectangular plates with linearly varying thickness. The fundamental eigenvalues are given for several taper and aspect ratios.Nomenclature a, b dimensions of plates - A _i series coefficients - D Eh ³/12(1– ²) flexural rigidity - D ₀ Eh ₀ ³ /12(1– ²) - E Young's modulus - h thickness, h ₀(1+x) - h ₀ thickness parameter - N _x , N _y stress resultants in the X and Y directions - N (N _x +N _y )/(1+) - P ₁, P ₂, ... parameters - Q ₁, Q ₂, ... parameters - R[X, (A/h ₀)²] bounding function - t time - u, v in-plane displacements - lateral deflections of plate - X=x/a dimensionless co-ordinate - x, y rectangular co-ordinates - y _n (X) series related to - thickness taper ratio - parameter in the neighbourhood of - error-function associated with differential equation - eigenvalue relating to frequency - Poisson's ra-tio - plate material specific weight - (X) function related to plate deflection - (X) admissible functions - circular frequency     

Singular Limits of a Two-Dimensional Boundary Value Problem Arising in Corrosion Modelling

We consider the boundary value problem where Ω is a smooth and bounded domain in ℝ² and λ > 0. We prove that for any integer k ≧ 1 there exist at least two solutions u _λ with the property that the boundary flux satisfies up to subsequences λ → 0, where the ξ _j are points of ∂Ω ordered clockwise in j.