On the controllability for the Khokhlov–Zabolotskaya–Kuznetsov-like equation

Recalling the proprieties of the Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation, we prove the controllability of moments result for the linear part of the KZK equation and its non-linear perturbation.     

On the constraints problem for the Einstein–Yang–Mills–Higgs system

We provide proofs of some key propositions that were used in previous work by Dossa and Tadmon dealing with the characteristic initial value problem for the Einstein–Yang–Mills–Higgs (EYMH) system. The aforesaid proofs were missing, making the considered work difficult to understand. This work is presented with a view to have an almost self-contained paper. With this respect we completely recall the process of constructing initial data for the EYMH system on two intersecting smooth null hypersurfaces as done in the work of Dossa and Tadmon mentioned above. This is achieved by successfully adapting the hierarchical method set up by Rendall to solve the same problem for the Einstein equations in vacuum and with perfect fluid source. Many delicate calculations and expressions are given in details so as to address, in a forthcoming work, the issue of global resolution of the characteristic initial value problem for the EYMH system. The method obviously applies to the Einstein–Maxwell and the Einstein-scalar field models as well.     

On the Aleksandrov–Bakelman–Pucci estimate for the infinity Laplacian

We prove L ^∞ bounds and estimates of the modulus of continuity of solutions to the Poisson problem for the normalized infinity and p-Laplacian, namely $$\begin{array}{ll} -\Delta_p^N u=f\quad{\rm for } \; n < p \leq\infty.\end{array}$$ We are able to provide a stable family of results depending continuously on the parameter p. We also prove the failure of the classical Alexandrov–Bakelman–Pucci estimate for the normalized infinity Laplacian and propose alternate estimates.     

Dualities for the Domany–Kinzel Model

We study the Domany–Kinzel model, which is a class of discrete-time Markov processes in one-dimension with two parameters (p ₁,p ₂)[0,1]². When p ₁= and p ₂=(2– ²) with (,)[0,1]², the process can be identified with the mixed site-bond oriented percolation model on a square lattice with probabilities of a site being open and of a bond being open. This paper treats dualities for the Domany–Kinzel model _t ^A and the DKdual _t ^A starting from A. We prove that , as long as one of A,B is finite and p ₂p ₁.     

Tower-of-sets analysis for the Kise–Ibaraki–Mine algorithm

Kise et al. (Oper. Res. 26:121–126, 1978) give an O(n ²) time algorithm to find an optimal schedule for the single-machine number of late jobs problem with agreeable job release dates and due dates. Li et al. (Oper. Res. 58:508–509, 2010a) point out that their proof of optimality for their algorithm is incorrect by giving a counter-example. In this paper, using the concept of “tower-of-sets” from Lawler (Math. Comput. Model. 20:91–106, 1994), we construct the tower-of-sets of the early job set generated by the algorithm. Then we give a correct proof of optimality for the algorithm and show a new result that the early job set by the algorithm obtains not only the maximum number of jobs but also the smallest total processing time among all optimal schedules. The result can be applied to study the problems of the hard real-time systems.     

Riemann–Hilbert approach for the Camassa–Holm equation on the line

We present a Riemann–Hilbert problem formalism for the initial value problem for the Camassa–Holm equation $u_t?u_{txx}+2\omega u_x+3u{u}_x=2u_x{u}_{xx}+u{u}_{xxx}$ on the line (CH). We show that: (i) for all $\omega >0$, the solution of this problem can be obtained in a parametric form via the solution of some associated Riemann–Hilbert problem; (ii) for large time, it develops into a train of smooth solitons; (iii) for small ω, this soliton train is close to a train of peakons, which are piecewise smooth solutions of the CH equation for $\omega =0$. To cite this article: A. Boutet de Monvel, D. Shepelsky, C. R. Acad. Sci. Paris, Ser. I 343 (2006).     

Reorthogonalization for the Golub–Kahan–Lanczos bidiagonal reduction

The Golub–Kahan–Lanczos (GKL) bidiagonal reduction generates, by recurrence, the matrix factorization of $X \in \mathbb{R }^{m \times n}, m \ge n$ , given by $$\begin{aligned} X = UBV^T \end{aligned}$$ where $U \in \mathbb{R }^{m \times n}$ is left orthogonal, $V \in \mathbb{R }^{n \times n}$ is orthogonal, and $B \in \mathbb{R }^{n \times n}$ is bidiagonal. When the GKL recurrence is implemented in finite precision arithmetic, the columns of $U$ and $V$ tend to lose orthogonality, making a reorthogonalization strategy necessary to preserve convergence of the singular values. The use of an approach started by Simon and Zha (SIAM J Sci Stat Comput, 21:2257–2274, 2000) that reorthogonalizes only one of the two left orthogonal matrices $U$ and $V$ is shown to be very effective by the results presented here. Supposing that $V$ is the matrix reorthogonalized, the reorthogonalized GKL algorithm proposed here is modeled as the Householder Q–R factorization of $\left( \begin{array}{c} 0_{n \times k} \\ X V_k \end{array}\right) $ where $V_k = V(:,1:k)$ . That model is used to show that if $\varepsilon _M $ is the machine unit and $$\begin{aligned} \bar{\eta }= \Vert \mathbf{tril }(I-V^T\!~V)\Vert _F, \end{aligned}$$ where $\mathbf{tril }(\cdot )$ is the strictly lower triangular part of the contents, then: (1) the GKL recurrence produces Krylov spaces generated by a nearby matrix $X + \delta X$ , $\Vert \delta X\Vert _F = \mathcal O (\varepsilon _M + \bar{\eta }) \Vert X\Vert _F$ ; (2) singular values converge in the Lanczos process at the rate expected from the GKL algorithm in exact arithmetic on a nearby matrix; (3) a new proposed algorithm for recovering leading left singular vectors produces better bounds on loss of orthogonality and residual errors.     

Convergent discretizations for the Nernst–Planck–Poisson system

We propose and compare two classes of convergent finite element based approximations of the nonstationary Nernst–Planck–Poisson equations, whose constructions are motivated from energy versus entropy decay properties for the limiting system. Solutions of both schemes converge to weak solutions of the limiting problem for discretization parameters tending to zero. Our main focus is to study qualitative properties for the different approaches at finite discretization scales, like conservation of mass, non-negativity, discrete maximum principle, decay of discrete energies, and entropies to study long-time asymptotics.     

Criterion for the Solvability of the Weighted Cauchy Problem for an Abstract Euler–Poisson–Darboux Equation

In a Banach space E, we consider the abstract Euler–Poisson–Darboux equation u″(t) + kt^?1u′(t) = Au(t) on the half-line. (Here k ∈ ? is a parameter, and A is a closed linear operator with dense domain on E.) We obtain a necessary and sufficient condition for the solvability of the Cauchy problem u(0) = 0, lim _t→0+t ^k u′(t) = u₁, k < 0, for this equation. The condition is stated in terms of an estimate for the norms of the fractional power of the resolvent of A and its derivatives. We introduce the operator Bessel function with negative index and study its properties.     

On the Wave Operators for the Friedrichs–Faddeev Model

We provide new formulae for the wave operators in the context of the Friedrichs–Faddeev model. Continuity with respect to the energy of the scattering matrix and a few results on eigenfunctions corresponding to embedded eigenvalues are also derived.     

Algorithm for solving the Navier–Stokes equations for the modeling of creeping flows

We study a coupled algorithm for solving the two-dimensional Navier–Stokes equations in the stream function–vorticity variables. The algorithm is based on a finite-difference scheme in which the inertial terms in the vortex transport equation are taken from the lower time layer and the dissipative terms, from the upper time layer. In the linear approximation, we study the stability of this algorithm and use test computations to show its advantages when modeling flows at moderate Reynolds numbers.     

Second-order implicit–explicit scheme for the Gray–Scott model

A second-order scheme for the Gray–Scott (GS) model used to describe the pattern formation is studied. The linear part of the GS equation for the time derivative and the viscous terms is discretized implicitly, while the other (or nonlinear) part of the GS equation explicitly. Galerkin finite element approximation methods are presented and analyzed, as well as methods for solving the resulting system of equations. The optimal L²$L^2$-norm error estimates are derived. Numerical experiments are presented.     

Heat Flow for the Yang–Mills–Higgs Field and the Hermitian Yang–Mills–Higgs Metric

For a parameter > 0, we study a type of vortex equations, which generalize the well-known Hermitian–Einstein equation, for a connection A and a section of a holomorphic vector bundle E over a Kähler manifold X. We establish a global existence of smooth solutions to heat flow for a self-dual Yang–Mills–Higgs field on E. Assuming the -stability of (E, ), we prove the existence of the Hermitian Yang–Mills–Higgs metric on the holomorphic bundle E by studying the limiting behaviour of the gauge flow.     

On the wave-breaking phenomena for the periodic two-component Dullin–Gottwald–Holm system

Considered herein is the well-posedness problem of the periodic two-component Dullin–Gottwald–Holm (DGH) system on the circle, which can be derived from Euler?s equation with constant vorticity in shallow water waves moving over a linear shear flow. The result of blow-up solutions for certain initial profiles in a manner which corresponds to wave-breaking is established.