Stability interval for explicit difference schemes for multi-dimensional second-order hyperbolic equations with significant first-order space derivative terms

In this piece of work, we introduce a new idea and obtain stability interval for explicit difference schemes of O(k²+h²) for one, two and three space dimensional second-order hyperbolic equations u_tt=a(x,t)u_xx+α(x,t)u_x-2η²(x,t)u,u_tt=a(x,y,t)u_xx+b(x,y,t)u_yy+α(x,y,t)u_x+β(x,y,t)u_y-2η²(x,y,t)u, and u_tt=a(x,y,z,t)u_xx+b(x,y,z,t)u_yy+c(x,y,z,t)u_zz+α(x,y,z,t)u_x+β(x,y,z,t)u_y+γ(x,y,z,t)u_z-2η²(x,y,z,t)u,0<x,y,z<1,t>0 subject to appropriate initial and Dirichlet boundary conditions, where h>0 and k>0 are grid sizes in space and time coordinates, respectively. A new idea is also introduced to obtain explicit difference schemes of O(k²) in order to obtain numerical solution of u at first time step in a different manner.     

Multiple positive solutions of nonhomogeneous elliptic systems with strongly indefinite structure and critical Sobolev exponents

In this paper, we study the existence of multiple positive solutions to some Hamiltonian elliptic systems −Δv=λu+u^p+εf(x), −Δu=μv+v^q+δg(x) in Ω;u,v>0 in Ω; u=v=0 on ∂Ω, where Ω is a bounded domain in R^N (N?3); 0?f, g∈L∞(Ω); 1/(p+1)+1/(q+1)=(N−2)/N, p,q>1; λ,μ>0. Using sub- and supersolution method and based on an adaptation of the dual variational approach, we prove the existence of at least two nontrivial positive solutions for all λ,μ∈(0,λ₁) and ε,δ∈(0,δ₀), where λ₁ is the first eigenvalue of the Laplace operator −Δ with zero Dirichlet boundary conditions and δ₀ is a positive number.     

Dynamics of a family of piecewise-linear area-preserving plane maps III. Cantor set spectra

This paper studies the behavior under iteration of the maps T _ab (x,y) = (F _ab (x) ? y, x) of the plane ?², in which F _ab (x) = ax if x ≥ 0 and bx if x < 0. These maps are area-preserving homeomorphisms of ?² that map rays from the origin to rays from the origin. Orbits of the map correspond to solutions of the nonlinear difference equation x _n+2 = 1/2(a ? b)|x _n+1|+1/2(a+b)x _n+1 ? x _n . This difference equation can be rewritten in an eigenvalue form for a nonlinear difference operator of Schrödinger type ? x _n+2+2x _n+1 ? x _n +V _μ(x _n+1)x _n+1 = Ex _n+1, in which μ = (1/2)(a ? b) is fixed, and V _μ(x) = μ(sgn(x)) is an antisymmetric step function potential, and the energy E = 2 ? 1/2(a+b). We study the set Ω _SB of parameter values where the map T _ab has at least one bounded orbit, which correspond to l _∞-eigenfunctions of this difference operator. The paper shows that for transcendental μ the set Spec_∞[μ] of energy values E having a bounded solution is a Cantor set. Numerical simulations suggest the possibility that these Cantor sets have positive (one-dimensional) measure for all real values of μ.     

An implicit function theorem

Suppose thatF:DR ⁿ×R^mRⁿ, withF(x ⁰,y ⁰)=0. The classical implicit function theorem requires thatF is differentiable with respect tox and moreover that ₁ F(x ⁰,y ⁰) is nonsingular. We strengthen this theorem by removing the nonsingularity and differentiability requirements and by replacing them with a one-to-one condition onF as a function ofx.     

The Instability of Nonmonotonic Wave Solutions of Parabolic Equations

We consider the class of equations u_t=f(u_xx, u_x, u) under the restriction that for all a,b,c. We first consider this equation over the unbounded domain ? ∞ < x < + ∞, and we show that very nearly every bounded nonmonotonic solution of the form u(t, x)=?(x?ct) is unstable to all nonnegative and all nonpositive perturbations. We then extend these results to nonmonotonic plane wave solutions u(t, x, y)=?(x?ct) of u_t = F(u_xx, u_xy, u_x, u_y, u). Finally, we consider the class of equations u_t=f(u_xx, u_x, u) over the bounded domain 0 < x < 1 with the boundary conditions u(t, x)=A at x=0 and u(t, x)=B at x=1, and we find the stability of all steady solutions u(t, x)=?(x).     

Circuits including a given set of vertices

Let G = (V, E) be a digraph of order n, satisfying Woodall's condition ? x, y ∈ V, if (x, y) ? E, then d⁺(x) + d^?(y) ≥ n. Let S be a subset of V of cardinality s. Then there exists a circuit including S and of length at most Min(n, 2s). In the case of oriented graphs we obtain the same result under the weaker condition d⁺(x) + d^?(y) ≥ n – 2 (which implies hamiltonism).     

A finite difference solution to a two-dimensional parabolic inverse problem

In this paper, we consider the problem of finding u = u(x, y, t) and p = p(t) which satisfy u_t = u_xx + u_yy + p(t)u + ? in R × [0, T], u(x, y, 0) = f(x, y), (x, y) ∈ R = [0, 1] × [0, 1], u is known on the boundary of R and u(x^∗, y^∗, t) = E(t), 0 < t ? T, where E(t) is known and (x^∗, y^∗) is a given point of R. Through a function transformation, the nonlinear two-dimensional diffusion problem is transformed into a linear problem, and a backward Euler scheme is constructed. It is proved by the maximum principle that the scheme is uniquely solvable, unconditionally stable and convergent in L_∞ norm. The convergence orders of u and p are of O(τ + h²). The impact of initial data errors on the numerical solution is also considered. Numerical experiments are presented to illustrate the validity of the theoretical results.     

Nonlinear equations with mixed derivatives

The nonlinear hyperbolic equation ∂²u(x, y)/∂x ∂y + g(x, y)f(u(x, y)) = 0 with u(x, 0) = φ(x) and u(0, y) = Ψ(y), considered by [1.], 31–45) under appropriate smoothness conditions, is solvable by the author's decomposition method (“Stochastic Systems,” Academic Press, 1983 and “Nonlinear Stochastic Operator Equations,” Academic Press, 1986).     

A new continuation method for complementarity problems with uniformP-functions

The complementarity problem with a nonlinear continuous mappingf from the nonnegative orthantR ₊ ⁿ ofR ⁿ intoR ⁿ can be written as the system of equationsF(x, y) = 0 and(x, y) R ₊ ²ⁿ , whereF denotes the mapping from the nonnegative orthantR ₊ ²ⁿ ofR ²ⁿ intoR ₊ ⁿ × Rⁿ defined byF(x, y) = (x ₁y₁,,x_ny_n, f₁(x) – y₁,, f_n(x) – y_n) for every(x, y) R ₊ ²ⁿ . Under the assumption thatf is a uniformP-function, this paper establishes that the mappingF is a homeomorphism ofR ₊ ²ⁿ ontoR ₊ ⁿ × Rⁿ. This result provides a theoretical basis for a new continuation method of tracing the solution curve of the one parameter family of systems of equationsF(x, y) = tF(x ⁰, y⁰) and(x, y) R ₊ ²ⁿ from an arbitrary initial point(x ⁰, y⁰) R ₊ ²ⁿ witht = 1 until the parametert attains 0. This approach is an extension of the one used in the polynomially bounded algorithm recently given by Kojima, Mizuno and Yoshise for solving linear complementarity problems with positive semi-definite matrices.     

Nonlocal initial-boundary-value problems for a degenerate hyperbolic equation

We consider the equation y ^m u _xx − u _yy − b ² y ^m u = 0 in the rectangular area {(x, y) | 0 < x < 1, 0 < y < T}, where m < 0, b ≥ 0, T > 0 are given real numbers. For this equation we study problems with initial conditions u(x, 0) = τ(x), u _y (x, 0) = ν(x), 0 ≤ x ≤ 1, and nonlocal boundary conditions u(0, y) = u(1, y), u _x (0, y) = 0 or u _x (0, y) = u _x (1, y), u(1, y) = 0 with 0≤y≤T. Using the method of spectral analysis, we prove the uniqueness and existence theorems for solutions to these problems     

Differential Invariants of Second Order ODEs,I

This paper is devoted to differential invariants of equations

Existence and Uniqueness of Strong Solutions for a Class of Quasi‐Linear Hyperbolic Equations with Order Degeneration

Let k(y) > 0, 𝓁(y) > 0 for y > 0, k(0) = 𝓁(0) = 0 and lim_{y → 0}k(y)/𝓁(y) exists; then the equation L(u) ≔ k(y)u_xx – ∂_y(𝓁(y)u_y) + a(x, y)u_x = f(x, y, u) is strictly hyperbolic for y > 0 and its order degenerates on the line y = 0. Consider the boundary value problem Lu = f(x, y, u) in G, u|_AC = 0, where G is a simply connected domain in ℝ² with piecewise smooth boundary ∂G = AB∪AC∪BC; AB = {(x, 0) : 0 ≤ x ≤ 1}, AC : x = F(y) = ∫^y₀(k(t)/𝓁(t))^1/2dt and BC : x = 1 – F(y) are characteristic curves. Existence of generalized solution is obtained by a finite element method, provided f(x, y, u) satisfies Carathéodory condition and |f(x, y, u)| ≤ Q(x, y) + b|u| with Q ∈ L²(G), b = const > 0. It is shown also that each generalized solution is a strong solution, and that fact is used to prove uniqueness under the additional assumption |f(x, y, u₁) – f(x, y, u₂| ≤ C|u₁ – u₂|, where C = const > 0.     

Asymptotic behaviour of a Brownian motion on exterior domains

We study the asymptotic behaviour of the transition density of a Brownian motion in ?, killed at ∂?, where ? ^c is a compact non polar set. Our main result concern dimension d = 2, where we show that the transition density p ^? _t (x, y) behaves, for large t, as u(x)u(y)(t(log t)²)⁻¹ for x, y∈?, where u is the unique positive harmonic function vanishing on (∂?) ^r , such that u(x) ∼ log ∣x∣. Received: 29 January 1999 / Revised version: 11 May 1999     

On Cesáro summability of fourier series with respect to double complete orthonormal systems

Let {ϕn(x), n = 1, 2,...} be an arbitrary complete orthonormal system on the interval I:= [0, 1]which consists of a.e. bounded functions. Suppose that E ₀ ⊂ I ² is any Lebesgue measurable set such that μ₂ E ₀ > 0, and φ, φ(0) = 0, is an increasing continuous function on [0, ∞) with φ(u) = o(u ln u) as u → ∞. Then there exist a function f ∈ L¹(I ²) and a set E ₀ ^′ , ⊂ E ₀, μ₂ E ₀ ^′ > 0, such that

Correlation structure of intermittency in the parabolic Anderson model

Consider the Cauchy problem ∂u(x, t)/∂t = ℋu(x, t) (x∈ℤ^d, t≥ 0) with initial condition u(x, 0) ≡ 1 and with ℋ the Anderson Hamiltonian ℋ = κΔ + ξ. Here Δ is the discrete Laplacian, κ∈ (0, ∞) is a diffusion constant, and ξ = {ξ(x): x∈ℤ ^d } is an i.i.d.random field taking values in ℝ. G?rtner and Molchanov (1990) have shown that if the law of ξ(0) is nondegenerate, then the solution u is asymptotically intermittent. In the present paper we study the structure of the intermittent peaks for the special case where the law of ξ(0) is (in the vicinity of) the double exponential Prob(ξ(0) > s) = exp[−e ^{s /θ}] (s∈ℝ). Here θ∈ (0, ∞) is a parameter that can be thought of as measuring the degree of disorder in the ξ-field. Our main result is that, for fixed x, y∈ℤ ^d and t→∈, the correlation coefficient of u(x, t) and u(y, t) converges to ∥w _ρ∥⁻² _ℓ2Σ_{z ∈ℤd} w _ρ(x+z)w _ρ(y+z). In this expression, ρ = θ/κ while w _ρ:ℤ^d→ℝ⁺ is given by w _ρ = (v _ρ)^{⊗ d} with v _ρ: ℤ→ℝ⁺ the unique centered ground state (i.e., the solution in ℓ²(ℤ) with minimal l ²-norm) of the 1-dimensional nonlinear equation Δv + 2ρv log v = 0. The uniqueness of the ground state is actually proved only for large ρ, but is conjectured to hold for any ρ∈ (0, ∞). empty It turns out that if the right tail of the law of ξ(0) is thicker (or thinner) than the double exponential, then the correlation coefficient of u(x, t) and u(y, t) converges to δ _{x, y} (resp.the constant function 1). Thus, the double exponential family is the critical class exhibiting a nondegenerate correlation structure. Received: 5 March 1997 / Revised version: 21 September 1998     

Solution of a restricted Hadamard problem on Minkowski spaces

We give a complete solution of Hadamard's problem concerning Huygens' principle on even-dimensional Minkowski spaces 𝕄ⁿ⁺¹ for a restricted class of linear, second-order, normal hyperbolic operators ℒ︁ = □_n+1 + u(x¹, x²) with real (locally) analytic potentials u = u(x¹, x²) depending on two spatial variables and homogeneous of degree −2. © 1997 John Wiley & Sons, Inc.     

Eigenvalue Intervals for 2mth Order Sturm—Liouville Boundary Value Problems

Values of λ are determined for which there exist positive solutions of the 2mth order differential equation on a measure chain, (-1)^m x ^?2m (t)=λa(t)f(u(σ(t))), y? [0,1], satisfying α _i+1 u ^?21(0)+0, γ _i+1 u ^?21(σ(1))=0, 0≤i≤m?1 with αi,βi,γi,δi≥0, where a and f are positive valued, and both lim _x-0+ (f(x)/x) and lim _x-0+ (f(x)/x) exist.     

Spreading Properties and Complex Dynamics for Monostable Reaction–Diffusion Equations

This paper is concerned with the study of the large-time behavior of the solutions u of a class of one-dimensional reaction–diffusion equations with monostable reaction terms f, including in particular the classical Fisher-KPP nonlinearities. The nonnegative initial data u ₀(x) are chiefly assumed to be exponentially bounded as x tends to + ∞ and separated away from the unstable steady state 0 as x tends to ? ∞. On the one hand, we give some conditions on u ₀ which guarantee that, for some λ > 0, the quantity c _λ = λ +f′(0)/λ is the asymptotic spreading speed, in the sense that lim  _t→+∞ u(t, ct) = 1 (the stable steady state) if c < c _λ and lim  _t→+∞ u(t, ct) = 0 if c > c _λ. These conditions are fulfilled in particular when u ₀(x) e ^λx is asymptotically periodic as x → + ∞. On the other hand, we also construct examples where the spreading speed is not uniquely determined. Namely, we show the existence of classes of initial conditions u ₀ for which the ω-limit set of u(t, ct + x) as t tends to + ∞ is equal to the whole interval [0, 1] for all x ∈ ? and for all speeds c belonging to a given interval (γ₁, γ₂) with large enough γ₁ < γ₂ ≤ + ∞.     

Infinitely many solutions for p‐Kirchhoff equation with concave–convex nonlinearities in

In this paper, we study the existence of infinitely many solutions to p‐Kirchhoff‐type equation (0.1) where f(x,u) = λh₁(x)|u|^{m ? 2}u + h₂(x)|u|^{q ? 2}u,a≥0,μ > 0,τ > 0,λ≥0 and . The potential function verifies , and h₁(x),h₂(x) satisfy suitable conditions. Using variational methods and some special techniques, we prove that there exists λ₀>0 such that problem 0.1 admits infinitely many nonnegative high‐energy solutions provided that λ∈[0,λ₀) and . Also, we prove that problem 0.1 has at least a nontrivial solution under the assumption f(x,u) = h₂|u|^{q ? 2}u,p < q< min{p*,p(τ + 1)} and has infinitely many nonnegative solutions for f(x,u) = h₁|u|^{m ? 2}u,1 < m < p. Copyright © 2015 John Wiley & Sons, Ltd.     

On the number of y-smooth natural numbers ≤x representable as a sum of two integer squares

LetB(x,y) be the sum taken over alln, 1≤n≤x, such that n can be represented as a sum of two squares of integers andn has no prime factors exceedingy. It is shown foru smaller than about .5log logx/log log logx thatB(x,x ^1/u)≈cxlog^-1/2 xσ(u), where σ(u satisfies a delay differential equation similar to the one satisfied by the Dickman function andc is a positive constant.