A new off-step high order approximation for the solution of three-space dimensional nonlinear wave equations

In this paper, we propose a new high accuracy numerical method of O(k² + k²h² + h⁴) based on off-step discretization for the solution of 3-space dimensional non-linear wave equation of the form u_tt = A(x,y,z,t)u_xx + B(x,y,z,t)u_yy + C(x,y,z,t)u_zz + g(x,y,z,t,u,u_x,u_y,u_z,u_t), 0 < x,y,z < 1,t > 0 subject to given appropriate initial and Dirichlet boundary conditions, where k > 0 and h > 0 are mesh sizes in time and space directions respectively. We use only seven evaluations of the function g as compared to nine evaluations of the same function discussed in  and . We describe the derivation procedure in details of the algorithm. The proposed numerical algorithm is directly applicable to wave equation in polar coordinates and we do not require any fictitious points to discretize the differential equation. The proposed method when applied to a telegraphic equation is also shown to be unconditionally stable. Comparative numerical results are provided to justify the usefulness of the proposed method.     

An inverse problem for computing a leading coefficient in the Sturm-Liouville operator by using the boundary data

We consider an inverse problem for identifying a leading coefficient α(x) in −(α(x)y′(x))′ + q(x)y(x) = H(x), which is known as an inverse coefficient problem for the Sturm-Liouville operator. We transform y(x) to u(x, t) =  (1 + t)y(x) and derive a parabolic type PDE in a fictitious time domain of t. Then we develop a Lie-group adaptive method (LGAM) to find the coefficient function α(x). When α(x) is a continuous function of x, we can identify it very well, by giving boundary data of y, y′ and α. The efficiency of LGAM is confirmed by comparing the numerical results with exact solutions. Although the data used in the identification are limited, we can provide a rather accurate solution of α(x).     

Duality for a Class of Multiobjective Control Problems

In this paper, a class of multiobjective control problems is considered, where the objective and constraint functions involved are f(t, x(t), ?(t), y(t), z(t)) with x(t) ∈ Rⁿ, y(t) ∈ Rⁿ, and z(t) ∈ R^m, where x(t) and z(t) are the control variables and y(t) is the state variable. Under the assumption of invexity and its generalization, duality theorems are proved through a parametric approach to related properly efficient solutions of the primal and dual problems.     

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.     

Distributive equations of implications based on nilpotent triangular norms

In this paper, we explore the distributive equations of implications, both independently and along with other equations. In detail, we consider three classes of equations. (1) By means of the section of I, we give out the sufficient and necessary conditions of solutions for the distributive equation of implication I(x, T(y, z)) = T(I(x, y), (x, z)) based on a nilpotent triangular norm T and an unknown function I, which indicates that there are no continuous solutions satisfying the boundary conditions of implications. Under the assumptions that I is continuous except the vertical section I(0, y), y ∈ [0, 1), we get its complete characterizations. (2) We prove that there are no solutions for the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, I(y, z)) = I(T(x, y), z). (3) We obtain the sufficient and necessary conditions on T and I to be solutions of the functional equations I(x, T(y, z)) = T(I(x, y), I(x, z)), I(x, y) = I(N(y), N(x)).     

Automatic Controls for Nonlinear Dynamical Systems with Lipschitzian Trajectories

We would like to investigate on the solution to the automatic control problem given by the differential equation y′(t) = f(t, y(t), w(t)) for a given initial function x in the initial domain D(x, ω, Y) for almost all t in the interval I, with controls given by w(t) = g(t, y(t), T(y)(t)), where T is a nonanticipating and Lipschitzian operator. The result will be generalized for a dynamical system y′(t) = f(t, y(t), T(y), u(t)).     

A nonlinear evolution system with two subdifferentials and monotone differential games

It is shown that the first order multivalued equation for V = V(t, x, y, z) involving the sum of two subdifferentials composed with the partials of V (V_t +f(t, x, y, z) · ▽_xV + β(V_y) + γ(V_z) + h(t, x, y, z) ? 0 a.e.) has a Lipschitz solution. This solution is shown to be the value of a differential game in which the players are restricted to choosing monotone nondecreasing functions of time. Accordingly, the multivalued equation is interpreted as the corresponding Hamilton-Jacobi equation of the game.     

Fractal oscillations of self-adjoint and damped linear differential equations of second-order

For a prescribed real number s ∈ [1, 2), we give some sufficient conditions on the coefficients p(x) and q(x) such that every solution y = y(x), y ∈ C²((0, T]) of the linear differential equation (p(x)y′)′ + q(x)y = 0 on (0, T], is bounded and fractal oscillatory near x = 0 with the fractal dimension equal to s. This means that y oscillates near x = 0 and the fractal (box-counting) dimension of the graph Γ(y) of y is equal to s as well as the s dimensional upper Minkowski content (generalized length) of Γ(y) is finite and strictly positive. It verifies that y admits similar kind of the fractal geometric asymptotic behaviour near x = 0 like the chirp function y_ch(x) = a(x)S(φ(x)), which often occurs in the time-frequency analysis and its various applications. Furthermore, this kind of oscillations is established for the Bessel, chirp and other types of damped linear differential equations given in the form y″ + (μ/x)y′ + g(x)y = 0, x ∈ (0, T]. In order to prove the main results, we state a new criterion for fractal oscillations near x = 0 of real continuous functions which essentially improves related one presented in [1].     

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.     

Meromorphic Functions That Share Fixed-Points

In this paper, we prove the following result: Let f(z) and g(z) be two nonconstant meromorphic(entire) functions, n ≥ 11(n ≥ 6) a positive integer. If fⁿ(z)f′(z) and gⁿ(z)g′(z) have the same fixed-points, then either f(z) = c₁e^cz2, g(z) = c₂e^− cz2, where c₁, c₂, and c are three constants satisfying 4(c₁c₂)^n + 1c² = −1, or f(z) ≡ tg(z) for a constant t such that t^n + 1 = 1.     

Multiple positive solutions to a singular boundary value problem for a superlinear Emden-Fowler equation

A multiplicity result for the singular ordinary differential equation y^″+λx⁻²yσ=0, posed in the interval (0,1), with the boundary conditions y(0)=0 and y(1)=γ, where σ>1, λ>0 and γ?0 are real parameters, is presented. Using a logarithmic transformation and an integral equation method, we show that there exists Σ^?∈(0,σ/2] such that a solution to the above problem is possible if and only if λγσ−1?Σ^?. For 0<λγσ−1<Σ^?, there are multiple positive solutions, while if γ=(λ⁻¹Σ^?)^1/(σ−1) the problem has a unique positive solution which is monotonic increasing. The asymptotic behavior of y(x) as x→⁰⁺ is also given, which allows us to establish the absence of positive solution to the singular Dirichlet elliptic problem −Δu=d⁻²(x)uσ in Ω, where Ω⊂RN, N?2, is a smooth bounded domain and d(x)=dist(x,∂Ω).     

A variational approach to dislocation problems for periodic Schrödinger operators

As a simple model for lattice defects like grain boundaries in solid state physics we consider potentials which are obtained from a periodic potential V=V(x,y) on R² with period lattice Z² by setting Wt(x,y)=V(x+t,y) for x<0 and Wt(x,y)=V(x,y) for x?0, for t∈[0,1]. For Lipschitz-continuous V it is shown that the Schrödinger operators Ht=−Δ+Wt have spectrum (surface states) in the spectral gaps of H₀, for suitable t∈(0,1). We also discuss the density of these surface states as compared to the density of the bulk. Our approach is variational and it is first applied to the well-known dislocation problem (Korotyaev (2000, 2005) [15] and [16]) on the real line. We then proceed to the dislocation problem for an infinite strip and for the plane. In Appendix A, we discuss regularity properties of the eigenvalue branches in the one-dimensional dislocation problem for suitable classes of potentials.     

A nonoscillation theorem for superlinear Emden-Fowler equations with near-critical coefficients

We are interested in the oscillatory behavior of solutions of the Emden-Fowler equation y^″+a(x)|y|^γ−1y=0, γ>1, where a(x) is a positive continuous function on (0,∞). In the special case when the coefficient a(x) is a power of x, i.e. a(x)=xα for some constant α, the value α^∗=−(γ+3)/2 plays a critical role: The equation has both oscillatory and nonoscillatory solutions if α>α^∗, while all solutions are nonoscillatory if α<α^∗. When a(x) is close to the critical exponent, one of the known results is that if a(x)=x−(γ+3)/2log−σ(x), where σ>0, then all solutions are nonoscillatory. In this paper, this result is further extended to include a class of coefficients in which the above condition with log(x) can be replaced by loglog(x), or logloglog(x) and so on.     

Generalized stability of the Cauchy and the Jensen functional equations on spheres

We study a generalized stability problem for Cauchy and Jensen functional equations satisfied for all pairs of vectors x,y from a linear space such that γ(x)=γ(y) or γ(x+y)=γ(x−y) with a given function γ.     

Stable periodic traveling waves for a predator-prey model with non-constant death rate and delay

In this paper we will consider a predator-prey model with a non-constant death rate and distributed delay, described by a partial integro-differential system. The main goal of this work is to prove that the partial integro-differential system has periodic orbitally asymptotically stable solutions in the form of periodic traveling waves; i.e. N(x, t) = N(σt − μ · x), P(x, t) = P(σt − μ · x), where σ > 0 is the angular frequency and μ is the vector number of the plane wave, which propagates in the direction of the vector μ with speed c = σ/∥μ∥; and N(x, t) and P(x, t) are the spatial population densities of the prey and the predator species, respectively. In order to achieve our goal we will use singular perturbation’s techniques.     

On the chromaticity of complete multipartite graphs with certain edges added

Let P(G,λ) be the chromatic polynomial of a graph G. A graph G is chromatically unique if for any graph H, P(H,λ)=P(G,λ) implies H is isomorphic to G. For integers k≥0, t≥2, denote by K((t−1)×p,p+k) the complete t-partite graph that has t−1 partite sets of size p and one partite set of size p+k. Let K(s,t,p,k) be the set of graphs obtained from K((t−1)×p,p+k) by adding a set S of s edges to the partite set of size p+k such that 〈S〉 is bipartite. If s=1, denote the only graph in K(s,t,p,k) by K⁺((t−1)×p,p+k). In this paper, we shall prove that for k=0,1 and p+k≥s+2, each graph G∈K(s,t,p,k) is chromatically unique if and only if 〈S〉 is a chromatically unique graph that has no cut-vertex. As a direct consequence, the graph K⁺((t−1)×p,p+k) is chromatically unique for k=0,1 and p+k≥3.     

Asymptotic behaviour of three-dimensional singularly perturbed convection-diffusion problems with discontinuous data

We consider three singularly perturbed convection-diffusion problems defined in three-dimensional domains: (i) a parabolic problem −?(uxx+uyy)+ut+v₁ux+v₂uy=0 in an octant, (ii) an elliptic problem −?(uxx+uyy+uzz)+v₁ux+v₂uy+v₃uz=0 in an octant and (iii) the same elliptic problem in a half-space. We consider for all of these problems discontinuous boundary conditions at certain regions of the boundaries of the domains. For each problem, an asymptotic approximation of the solution is obtained from an integral representation when the singular parameter ?→⁰⁺. The solution is approximated by a product of two error functions, and this approximation characterizes the effect of the discontinuities on the small ?− behaviour of the solution and its derivatives in the boundary layers or the internal layers.     

Periodic solutions to one-dimensional wave equation with x-dependent coefficients

In this paper, we study the nonlinear one-dimensional periodic wave equation with x-dependent coefficients u(x)ytt−(ux(x)yx)+g(x,t,y)=f(x,t) on (0,π)×R under the boundary conditions a₁y(0,t)+b₁yx(0,t)=0, a₂y(π,t)+b₂yx(π,t)=0 ( for i=1,2) and the periodic conditions y(x,t+T)=y(x,t), yt(x,t+T)=yt(x,t). Such a model arises from the forced vibrations of a nonhomogeneous string and the propagation of seismic waves in nonisotropic media. A main concept is the notion “weak solution” to be given in Section 2. For T is the rational multiple of π, we prove some important properties of the weak solution operator. Based on these properties, the existence and regularity of weak solutions are obtained.