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.     

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.     

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 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.     

Global Strichartz estimates for the wave equation with time-periodic potentials

We obtain global Strichartz estimates for the solutions u of the wave equation for time-periodic potentials V(t,x) with compact support with respect to x. Our analysis is based on the analytic properties of the cut-off resolvent Rχ(z)=χ(U⁻¹(T)−zI)ψ₁, where U(T)=U(T,0) is the monodromy operator and T>0 the period of V(t,x). We show that if Rχ(z) has no poles z∈C, |z|?1, then for n?3, odd, we have a exponential decal of local energy. For n?2, even, we obtain also an uniform decay of local energy assuming that Rχ(z) has no poles z∈C, |z|?1, and Rχ(z) remains bounded for z in a small neighborhood of 0.     

Recurrent points and non-wandering points of graph maps

Let G be a graph and f:G→G be a continuous map. Denote by P(f), R(f) and Ω(f) the sets of periodic points, recurrent points and non-wandering points of f, respectively. In this paper we show that: (1) If L=(x,y) is an open arc contained in an edge of G such that {fm(x),fk(y)}⊂(x,y) for some m,k∈N, then R(f)∩(x,y)≠∅; (2) Any isolated point of P(f) is also an isolated point of Ω(f); (3) If x∈Ω(f)−Ω(fn) for some n∈N, then x is an eventually periodic point. These generalize the corresponding results in W. Huang and X. Ye (2001) [9] and J. Xiong (1983, 1986) [17] and [19] on interval maps or tree maps.     

Heteroclinic solutions for a class of the second order Hamiltonian systems

We shall be concerned with the existence of heteroclinic orbits for the second order Hamiltonian system , where q∈Rn and V∈C¹(R×Rn,R), V?0. We will assume that V and a certain subset M⊂Rn satisfy the following conditions. M is a set of isolated points and #M?2. For every sufficiently small ε>0 there exists δ>0 such that for all (t,z)∈R×Rn, if d(z,M)?ε then −V(t,z)?δ. The integrals , z∈M, are equi-bounded and −V(t,z)→∞, as |t|→∞, uniformly on compact subsets of Rn?M. Our result states that each point in M is joined to another point in M by a solution of our system.     

A bilinear version of Orlicz-Pettis theorem

Given three Banach spaces X, Y and Z and a bounded bilinear map , a sequence x=n(xn)⊆X is called B-absolutely summable if is finite for any y∈Y. Connections of this space with are presented. A sequence x=n(xn)⊆X is called B-unconditionally summable if is finite for any y∈Y and z^∗∈Z^∗ and for any M⊆N there exists xM∈X for which _∑n∈M〈B(xn,y),z^∗〉=〈B(xM,y),z^∗〉 for all y∈Y and z^∗∈Z^∗. A bilinear version of Orlicz-Pettis theorem is given in this setting and some applications are presented.     

Generalized Jordan algebras

We study commutative algebras which are generalizations of Jordan algebras. The associator is defined as usual by (x, y, z) = (x y)z − x(y z). The Jordan identity is (x², y, x) = 0. In the three generalizations given below, t, β, and γare scalars. ((x x)y)x + t((x x)x)y = 0, ((x x)x)(y x) − (((x x)x)y)x = 0, β((x x)y)x + γ((x x)x)y − (β + γ)((y x)x)x = 0. We show that with the exception of a few values of the parameters, the first implies both the second and the third. The first is equivalent to the combination of ((x x)x)x = 0 and the third. We give examples to show that our results are in some reasonable sense, the best possible.     

An optimal filtering method for stable analytic continuation

In this paper, we consider the problem of numerical analytic continuation of an analytic function f(z)=f(x+iy) on a strip domain Ω⁺={z=x+iy∈C∣x∈R,0<y<y₀}, where the data is given approximately only on the real axis y=0. This problem is severely ill-posed: the solution does not depend continuously on the given data. A novel method (filtering) is used to solve this problem and an optimal error estimate with Hölder type is proved. Numerical examples show that this method works effectively.     

On solvability of two-dimensional Lp-Minkowski problem

Given p≠0 and a positive continuous function g, with g(x+T)=g(x), for some 0<T<1 and all real x, it is shown that for suitable choice of a constant C>0 the functional has a minimizer in the class of positive functions u∈C¹(R) for which u(x+T)=u(x) for all x∈R. This minimizer is used to prove the existence of a positive periodic solution y∈C²(R) of two-dimensional L_p-Minkowski problem y^1−p(x)(y″(x)+y(x))=g(x), where p∉{0,2}.     

One-step 5-stage Hermite-Birkhoff-Taylor ODE solver of order 12

A one-step 5-stage Hermite-Birkhoff-Taylor method, HBT(12)5, of order 12 is constructed for solving nonstiff systems of differential equations y^′=f(t,y), y(t₀)=y₀, where y∈Rⁿ. The method uses derivatives y^′ to y⁽⁹⁾ as in Taylor methods combined with a 5-stage Runge-Kutta method. Forcing an expansion of the numerical solution to agree with a Taylor expansion of the true solution to order 12 leads to Taylor- and Runge-Kutta-type order conditions which are reorganized into Vandermonde-type linear systems whose solutions are the coefficients of the method. HBT(12)5 has a larger interval of absolute stability than Dormand-Prince DP(8, 7)13M and Taylor method T12 of order 12. The new method has also a smaller norm of principal error term than T12. It is superior to DP(8, 7)13M and T12 on the basis the number of steps, CPU time and maximum global error on common test problems. The formulae of HBT(12)5 are listed in an appendix.     

Properties of the density for a three-dimensional stochastic wave equation

We consider a stochastic wave equation in space dimension three driven by a noise white in time and with an absolutely continuous correlation measure given by the product of a smooth function and a Riesz kernel. Let pt,x(y) be the density of the law of the solution u(t,x) of such an equation at points (t,x)∈]0,T]×R³. We prove that the mapping (t,x)?pt,x(y) owns the same regularity as the sample paths of the process {u(t,x),(t,x)∈]0,T]×R³} established in [R.C. Dalang, M. Sanz-Solé, Hölder-Sobolev regularity of the solution to the stochastic wave equation in dimension three, Mem. Amer. Math. Soc., in press]. The proof relies on Malliavin calculus and more explicitly, the integration by parts formula of [S. Watanabe, Lectures on Stochastic Differential Equations and Malliavin Calculus, Tata Inst. Fund. Res./Springer-Verlag, Bombay, 1984] and estimates derived from it.     

Polynomially bounded solutions to the Loewner differential equation in several complex variables

We determine the form of polynomially bounded solutions to the Loewner differential equation that is satisfied by univalent subordination chains of the form f(z,t)=etAz+?, where A∈L(Cn,Cn) has the property m(A)>0. Here m(A)=min{R〈A(z),z〉:‖z‖=1}. We also give sufficient conditions for g(z,t)=L(f(z,t)) to be polynomially bounded, where f(z,t) is an A-normalized polynomially bounded Loewner chain solution to the Loewner differential equation.     

A new inequality for symmetric functions

Let x_i ≥ 0, y_i ≥ 0 for i = 1,…, n; and let a_j(x) be the elementary symmetric function of n variables given by a_j(x) = ∑_{1 ≤ ii < … <ij ≤ n}x_ii … x_ij. Define the partical ordering x <y if a_j(x) ≤ a_j(y), j = 1,… n. We show that $\text{x}\text{\$}\text{?}\text{y ? x}{}^{\mathrm{\alpha }}\text{\$}\text{?}\text{y}{}^{\mathrm{\alpha }}\text{, 0}\text{\$}\text{?}\text{α ≤ 1}$, where {x^α}_i = x^α_i. We also give a necessary and sufficient condition on a function f(t) such that x <y ? f(x) <f(y). Both results depend crucially on the following: If x <y there exists a piecewise differentiable path z(t), with z_i(t) ≥ 0, such that z(0) = x, z(1) = y, and z(s) <z(t) if 0 ≤ s ≤ t ≤ 1.     

Control problems for parabolic equations on controlled domains

Given the one-dimensional heat equation v_t = v_xx on the controlled domain Q(y) = {(t, x); 0 < x < y(t), 0 < t < T} subject to some initial-boundary conditions, we study the problem of optimally selecting y(·) from some admissible class so as to maximize a given payoff of fixed duration. Q(y) is thus a controlled domain. We also study the problem in which the heat equation holds in Q(y, z) = {z(t) < x < y(t), 0 < t < T}; z minimizing, y maximizing, i.e., the differential game. The principle techniques involved are (i) transforming the controlled domain to an uncontrolled domain and then (ii) using the method of lines for parabolic equations to enable us to use known results for control systems governed by ordinary differential equations. Sufficient conditions for existence in an admissible class is given and the method of lines allows numerical techniques to be applied to determine the optimal control in our class.