A remark on free boundary problem of 1‐D compressible Navier–Stokes equations with density‐dependent viscosity

In this paper, we prove the existence and uniqueness of the weak solution of the one‐dimensional compressible Navier–Stokes equations with density‐dependent viscosity µ(ρ)=ρ^θ with θ∈(0, γ?2], γ>1. The initial data are a perturbation of a corresponding steady solution and continuously contact with vacuum on the free boundary. The obtained results apply for the one‐dimensional Siant–Venant model of shallow water and generalize ones in (Arch. Rational Mech. Anal. 2006; 182: 223–253). Copyright © 2009 John Wiley & Sons, Ltd.     

Robust a-posteriori error control of Cahn-Hilliard type equations with elasticity

Phase separation of an initially homogeneous mixture into two different phases can be modeled on a mesoscopic scale by the Cahn-Hilliard equation. The interface thickness between the pure phases enters as a small parameter γ into this mass conserving fourth order semilinear parabolic equation. Numerical analysis is well established for a fixed parameter size, but error estimates depend exponentially on γ^–1 and thus become useless if γ → 0. We consider the case, that elastic stresses due to a lattice misfit become important and the equation has to be coupled to a system of linear elasticity. Applications include e. g. the simulation of Sn-Cu alloys for the production of lead free solder or Ni-Al alloys used for rotor blade surfaces. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Canonical Diffraction Problem With Two Media

A Wiener–Hopf equation in L² being equivalent [5] to a boundary value problem (of the first kind) for a wave-scattering Sommerfeld half-plane Σ=ℝ₊×{0} which faces two different media Ω_-: x₂<0, Ω₊: x₂>0, as a special configuration in [3], is solved by canonical Weiner–Hopf factorization of its L²-regular scalar symbol γ_o=γ_o- γ_o+. The factors are calculated by solving a Riemann–Hilbert boundary value problem on the semi-infinite branch cuts of t_j(ξ):=(ξ²−k²_j)^1/2, k_j∈ℂ⁺⁺ for j=1,2: taken parallel to the imaginary axis. The procedure following this idea is known as the Wiener–Hopf–Hilbert(–Hurd) method [2] and requires the evaluation of elliptic-type integrals. Formula (3.7) seems not to be contained in tables of integrals.     

The degree sequence of a scale‐free random graph process

Recently, Barabási and Albert [2] suggested modeling complex real‐world networks such as the worldwide web as follows: consider a random graph process in which vertices are added to the graph one at a time and joined to a fixed number of earlier vertices, selected with probabilities proportional to their degrees. In [2] and, with Jeong, in [3], Barabási and Albert suggested that after many steps the proportion P(d) of vertices with degree d should obey a power law P(d)αd^?γ. They obtained γ=2.9±0.1 by experiment and gave a simple heuristic argument suggesting that γ=3. Here we obtain P(d) asymptotically for all d≤n^1/15, where n is the number of vertices, proving as a consequence that γ=3. © 2001 John Wiley & Sons, Inc. Random Struct. Alg., 18, 279–290, 2001     

Upper bound of the constant in the strengthened C.B.S. inequality for FEM 2D elasticity equations

The basic theory of the strengthened Cauchy–Buniakowskii–Schwarz (C.B.S.) inequality is the main tool in the convergence analysis of the recently proposed algebraic multilevel iterative methods. An upper bound of the constant γ in the strengthened C.B.S. inequality for the case of the finite element solution of 2D elasticity problems is obtained. It is assumed that linear triangle finite elements are used, the initial mesh consisting of right isosceles triangles and the mesh refinement procedure being uniform. For the resulting linear algebraic systems we have proved that γ²<0.75 uniformly on the mesh parameter and on Poisson's ratio ν ? (0, 1/2). Furthermore, the presented numerical tests show that the same relation holds for arbitrary initial right triangulations, even in the case of degeneracy of triangles. The theoretical results obtained are practically important for successful implementation of the finite element method to large-scale modeling of complicated structures. They allow us to construct optimal order algebraic multilevel iterative solvers for a wide class of real–life elasticity problems.     

Generalized impedance boundary condition at high frequency for a domain with thin layer: the circular case

Consider a conducting disk surrounded by a thin dielectric layer submitted to an electric field at the pulsation ω. The conductivity of the layer grows like ω^1?γ, γ∈[0,1], when the pulsation ω?tends to infinity. Using a pseudodifferential approach on the torus, we build an equivalent boundary condition with the help of an appropriate factorization of Helmholtz operator in the layer. This generalized impedance condition approximates the thin membrane in the high frequency limit for small thickness of the layer. L ²-error estimates are given and we illustrate our results with numerical simulations. This work extends, in the circular geometry, previous works of Lafitte and Lebeau (Lafitte O. Lebeau G. 1993, Équations de Maxwell et opérateur d’impédance sur le bord d’un obstacle convexe absorbant. Comptes Rendus de l ' Académic dis Science, Paris, Série I, Mathématiques, 316(11), 1177–1182); (Lafitte O.D., 1999, Diffraction in the high frequency regime by a thin layer of dielectric material. I. The equivalent impedance boundary condition. SIAM Journal on Applied Mathematics, 59(3), 1028–1052 (electronic)) in which γ?identically equals zero.     

Generalized conic curves and their applications in curve approximation

The appriximation properties of generalized conic curves are studied in this paper. A generalized conic curve is defined as one of the following curves or their affine and translation equivalent curves:

1. conic curves, including parabolas, hyperbolas and ellipses;
2. generalized monomial curves, including curves of the form x=y^γ, γ∈R, γ≠0,1, in the x?y Cartesian coordinate system;
3. exponential spiral curves of the form ρ(?)=Ae^γ?, A>0, γ≠0, in the ρ-? polar coordinate system.

This type of curves has many important properties such as convexity, approximation property, effective numerical computation property and the subdivision property etc. Applications of these curves in both interpolation and approximations using piecewise generalized conic segment are also developed. It is shown that these generalized conic splines are very similar to the cubic polynomial splines and the best error of approximation isO(h ⁵) or at leastO(h ⁴) in general provided appropriate procedures are used. Finally some numerical examples of interpolation and approximations with generalized conic splines are given.     

Periodic solutions for the transport of a pollutant concentration in a porous medium

We study the existence and uniqueness of the periodic solution for a model which describes the transport of a pollutant in a porous medium, with periodic boundary conditions. The prove of existence, is based upon a topological tool as the Schauder fixed point theorem, applied to a PDE. The uniqueness is instead of obtained, assuming thats→d _i,j(x,? ^–1(s)) is Hölder continuous of order γ∈ ∈ (1/2, 1).     

Note on the heights of random recursive trees and random m-ary search trees

A process of growing a random recursive tree T_n is studied. The sequence {T_n} is shown to be a sequence of “snapshots” of a Crump–Mode branching process. This connection and a theorem by Kingman are used to show quickly that the height of T_n is asymptotic, with probability one, to c log n. In particular, c = e = 2.718 … for the uniform recursive tree, and c = (2γ)^?1, where γe^1+γ = 1, for the ordered recursive tree. An analogous reduction provides a short proof of Devroye's limit law for the height of a random m-ary search tree. We show finally a close connection between another Devroye's result, on the height of a random union-find tree, and our theorem on the height of the uniform recursive tree. © 1994 John Wiley & Sons, Inc.     

Norm estimates related to self-commulators

For a given self-adjoint matrix H with ∥H∥ = 1 and tr(H) = 0, we consider the number γ(H) which is defined to be the minimum of ∥T∥² for those T satisfying [T¹,T] = H. We show that 1 ⩽ γ(H) ⩽ 2 and that γ(H) is close to 2 if H is suitably chosen.     

Asymptotic expansions for derivatives of stable laws

For the derivativesp ^(k)(x; α, γ) of the stable density of index α asymptotic formulae (of Plancherel Rotach type) are computed ask→∞ thereby exhibiting the detailed analytic structure for large orders of derivatives. Generalizing known results for the special case of the one-sided stable laws (O<α<1, γ=-α) the whole range for the index of stability and the asymmetry parameter γ is covered.     

A note on similarity solutions for boundary layer flows with prescribed heat flux

We study the third‐order nonlinear equation: f′′′ + (m + 2)ff′′ ? (2m + 1)f′² = 0 on (0, ∞), subject to the boundary conditions f(0) = ? γ ∈ ?, f′(∞) = 0 and f′′(0) = ?1. The problem arises in the study of similarity solutions for boundary layer flows with prescribed heat flux. We will address the following two open questions which were proposed by Brighi and Hoernel (Math. Methods Appl. Sci. 2005; 28 : 479–503): The first one is the uniqueness of bounded solutions for and γ>0, and the second one is the structure of solutions for and γ?0. Copyright © 2007 John Wiley & Sons, Ltd.     

Interior C 1, γ-Regularity for Weak Solutions of Nonlinear Second Order Parabolic Systems

The interior C ^{γ, γ/2}-regularity for the first gradient with respect to the space variables of a weak solution to a class of nonlinear second order parabolic systems is proven under the assumption that oscillations of coefficients are controlled by the ellipticity constant.     

A note on stable exponents

Given an element γ in a group γ, the stable exponent p+(γ) of γ is defined as p+(γ) =lim sup_n→∞P(γⁿ) denotes the exponent of P(γⁿ) = sup{k/ ?γ_o ∈ γ such that γⁿ = γ^k _o We prove that if γ acts properly discontinuously by isometrics on a proper geodesic Gromov-hyperbolic metric space and γ ∈ γ is of hyperbolic type, then P+(γ) is an integer. This implies that the stable exponent of every element of infinite order in a word hyperbolic group is an integer. We also show that, in a translation discrete group, the stable exponent of every element of infinite order is finite.     

Randomness-optimal oblivious sampling

We present the first efficient oblivious sampler that uses an optimal number of random bits, up to an arbitrary constant factor bigger than 1. Specifically, for any α>0, it uses (1+α)(m+log γ⁻¹) random bits to output d=poly(ϵ⁻¹, log γ⁻¹, m) sample points z₁,…,z_d∈{0, 1}^m such that for any function f: {0, 1}^m→[0, 1], Pr [|(1/d)∑_i=1^df(z_i)− E f|≤ϵ]≥1−γ. Our proof is based on an improved extractor construction. An extractor is a procedure which takes as input the output of a defective random source and a small number of truly random bits, and outputs a nearly random string. We present the first optimal extractor, up to constant factors, for defective random sources with constant entropy rate. We give applications to constructive leader election and reducing randomness in interactive proofs. © 1997 John Wiley & Sons, Inc. Random Struct. Alg., 11 , 345–367 (1997)     

Interior C1, γ‐regularity for weak solutions of nonlinear second orderelliptic systems

The interior C^{0, γ}‐regularity for the first gradient of a weak solution to a class of nonlinear second order elliptic systems is proved under the assumption that oscillations of coefficients are controlled by the ellipticity constant. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Vertex criticality for upper domination and irredundance

Let π be any of the domination parameters ir γ, i, β, Γ or IR. The graph G is π‐critical (π⁺‐critical) if the removal of any vertex of G causes π(G) to decrease (increase). We show that the classes of IR‐critical and Γ‐critical graphs coincide, and exhibit a class of Γ⁺‐critical graphs. © 2001 John Wiley & Sons, Inc. J Graph Theory 37: 205–212, 2001     

Infinite families of biembedding numbers

Let N(γ, γ′) denote the size of the smallest complete graph that cannot be edge-partitioned into two parts embeddable in closed orientable sufaces of genera γ, γ′, respectively. Well-known embedding theorems are used to obtain several infinite families of values of N(γ, γ′). Some related small values of N(γ, γ′) are also discussed.     

A characterization of domination perfect graphs

Let γ(G) and i(G) be the domination number and independent domination number of a graph G, respectively. Sumner and Moore [8] define a graph G to be domination perfect if γ(H) = i(H), for every induced subgraph H of G. In this article, we give a finite forbidden induced subgraph characterization of domination perfect graphs. Bollobás and Cockayne [4] proved an inequality relating γ(G) and i(G) for the class of K^1,k -free graphs. It is shown that the same inequality holds for a wider class of graphs.     

A numerical method for computing radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation

In this article we develop a finite‐difference scheme to approximate radially symmetric solutions of a dissipative nonlinear modified Klein‐Gordon equation subject to smooth initial conditions ? and ψ in an open sphere D around the origin, with constant internal and external damping coefficients—β and γ, respectively—, and nonlinear term of the form G′(w) = w^p, with p > 1 an odd number. The functions ? and ψ are radially symmetric in D, and ?, ψ, r?, and rψ are assumed to be small at infinity. We prove that our scheme is consistent order ??(Δt²) + ??(Δr²) for G′ identically equal to zero and provide a necessary condition for it to be stable order n. Part of our study will be devoted to compare the physical effects of β and γ. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005