Complete classification of local conservation laws for generalized Cahn–Hilliard–Kuramoto–Sivashinsky equation

In this paper, we consider nonlinear multidimensional Cahn–Hilliard and Kuramoto–Sivashinsky equations that have many important applications in physics and chemistry, and a certain natural generalization of these two equations to which we refer to as the generalized Cahn–Hilliard–Kuramoto–Sivashinsky equation. For an arbitrary number of spatial independent variables, we present a complete list of cases when the latter equation admits nontrivial local conservation laws of any order, and for each of those cases, we give an explicit form of all the local conservation laws of all orders modulo trivial ones admitted by the equation under study. In particular, we show that the original Kuramoto–Sivashinsky equation admits no nontrivial local conservation laws, and find all nontrivial local conservation laws for the Cahn–Hilliard equation.     

Local linear convergence analysis of Primal–Dual splitting methods

In this paper, we study the local linear convergence properties of a versatile class of Primal–Dual splitting methods for minimizing composite non-smooth convex optimization problems. Under the assumption that the non-smooth components of the problem are partly smooth relative to smooth manifolds, we present a unified local convergence analysis framework for these methods. More precisely, in our framework, we first show that (i) the sequences generated by Primal–Dual splitting methods identify a pair of primal and dual smooth manifolds in a finite number of iterations, and then (ii) enter a local linear convergence regime, which is characterized based on the structure of the underlying active smooth manifolds. We also show how our results for Primal–Dual splitting can be specialized to cover existing ones on Forward–Backward splitting and Douglas–Rachford splitting/ADMM (alternating direction methods of multipliers). Moreover, based on these obtained local convergence analysis result, several practical acceleration techniques are discussed. To exemplify the usefulness of the obtained result, we consider several concrete numerical experiments arising from fields including signal/image processing, inverse problems and machine learning. The demonstration not only verifies the local linear convergence behaviour of Primal–Dual splitting methods, but also the insights on how to accelerate them in practice.     

Optimal control of bioprocess systems using hybrid numerical optimization algorithms

In the paper, we consider the bioprocess system optimal control problem. Generally speaking, it is very difficult to solve this problem analytically. To obtain the numerical solution, the problem is transformed into a parameter optimization problem with some variable bounds, which can be efficiently solved using any conventional optimization algorithms, e.g. the improved Broyden–Fletcher–Goldfarb–Shanno algorithm. However, in spite of the improved Broyden–Fletcher–Goldfarb–Shanno algorithm is very efficient for local search, the solution obtained is usually a local extremum for non-convex optimal control problems. In order to escape from the local extremum, we develop a novel stochastic search method. By performing a large amount of numerical experiments, we find that the novel stochastic search method is excellent in exploration, while bad in exploitation. In order to improve the exploitation, we propose a hybrid numerical optimization algorithm to solve the problem based on the novel stochastic search method and the improved Broyden–Fletcher–Goldfarb–Shanno algorithm. Convergence results indicate that any global optimal solution of the approximate problem is also a global optimal solution of the original problem. Finally, two bioprocess system optimal control problems illustrate that the hybrid numerical optimization algorithm proposed by us is low time-consuming and obtains a better cost function value than the existing approaches.     

A Vanka-based parameter-robust multigrid relaxation for the Stokes–Darcy Brinkman problems

We consider a block-structured multigrid method based on Braess–Sarazin relaxation for solving the Stokes–Darcy Brinkman equations discretized by the marker and cell scheme. In the relaxation scheme, an element-based additive Vanka operator is used to approximate the inverse of the corresponding shifted Laplacian operator involved in the discrete Stokes–Darcy Brinkman system. Using local Fourier analysis, we present the stencil for the additive Vanka smoother and derive an optimal smoothing factor for Vanka-based Braess–Sarazin relaxation for the Stokes–Darcy Brinkman equations. Although the optimal damping parameter is dependent on meshsize and physical parameter, it is very close to one. In practice, we find that using three sweeps of Jacobi relaxation on the Schur complement system is sufficient. Numerical results of two-grid and V(1,1)-cycle are presented, which show high efficiency of the proposed relaxation scheme and its robustness to physical parameters and the meshsize. Using a damping parameter equal to one gives almost the same convergence results as these for the optimal damping parameter.     

Spreading speeds for time heterogeneous prey–predator systems with diffusion

We investigate the large time behavior for two components reaction–diffusion systems of prey–predator type in a time varying environment. Here we assume that these variations in time exhibit an averaging property, which will be called mean value in this work. This framework includes in particular time periodicity, almost periodicity and unique ergodicity. We describe the spreading behavior of the prey and the predator, wherein the two populations are able to co-invade the empty space. Our analysis is based the parabolic strong maximum principle for scalar equation and on the derivation of local pointwise estimates that are used to compare the solutions of the prey–predator problem with those of a KPP scalar equation on suitable spatio-temporal domains.     

Dynamics analysis of a reaction–diffusion system with Beddington–DeAngelis functional response and strong Allee effect

In this paper, the dynamics of a diffusive predator–prey model with modified Leslie–Gower term and strong Allee effect on prey under homogeneous Neumann boundary condition is considered. Firstly, we obtain the qualitative properties of the system including the existence of the global positive solution and the local and global asymptotical stability of the constant equilibria. In addition, we investigate a priori estimate and the nonexistence of nonconstant positive steady state solutions. Finally, we establish the existence and local structure of steady state patterns and time-periodic patterns for the system.     

A special complementarity function revisited

ABSTRACT

Recently, a local framework of Newton-type methods for constrained systems of equations has been developed. Applied to the solution of Karush–Kuhn–Tucker (KKT) systems, the framework enables local quadratic convergence under conditions that allow nonisolated and degenerate KKT points. This result is based on a reformulation of the KKT conditions as a constrained piecewise smooth system of equations. It is an open question whether a comparable result can be achieved for other (not piecewise smooth) reformulations. It will be shown that this is possible if the KKT system is reformulated by means of the Fischer–Burmeister complementarity function under conditions that allow degenerate KKT points and nonisolated Lagrange multipliers. To this end, novel constrained Levenberg–Marquardt subproblems are introduced. They allow significantly longer steps for updating the multipliers. Based on this, a convergence rate of at least 1.5 is shown.     

Bifurcation and stability of a two-species reaction–diffusion–advection competition model

This paper is concerned with the dynamics of a two-species reaction–diffusion–advection competition model subject to the no-flux boundary condition in a bounded domain. By the signs of the associated principal eigenvalues, we derive the existence and local stability of the trivial and semi-trivial steady-state solutions. Moreover, the nonexistence and existence of the coexistence steady-state solutions stemming from the two boundary steady states are obtained as well. In particular, we describe the feature of the coincidence of bifurcating coexistence steady-state solution branches. At the same time, the effect of advection on the stability of the bifurcating solution is also investigated, and our results suggest that the advection term may change the stability. Finally, we point out that the methods we applied here are mainly based on spectral analysis, perturbation theory, comparison principle, monotone theory, Lyapunov–Schmidt reduction, and bifurcation theory.     

Initial Boundary Value Problems of the Camassa–Holm Equation

In this paper we study initial boundary value problems of the Camassa–Holm equation on the half line and on a compact interval. Using rigorously the conservation of symmetry, it is possible to convert these boundary value problems into Cauchy problems for the Camassa–Holm equation on the line and on the circle, respectively. Applying thus known results for the latter equations we first obtain the local well-posedness of the initial boundary value problems under consideration. Then we present some blow-up and global existence results for strong solutions. Finally we investigate global and local weak solutions for the equation on the half line and on a compact interval, respectively. An interesting result of our analysis shows that the Camassa–Holm equation on a compact interval possesses no nontrivial global classical solutions.     

Hilbert-valued self-intersection local times for planar Brownian motion

Dynkin's construction for self-intersection local time of a planar Wiener process is extended to Hilbert-valued weights. In Dynkin's construction, the weight is bounded and measurable. Since the weight function describes the properties of the medium in which the Brownian motion moves, relative to the external medium's properties, the weight function can be random and unbounded. In this article, we discuss the possibility to consider Hilbert-space-valued weights. It appears that the existence of Hilbert-valued Dynkin-renormalized self-intersection local time is equivalent to the embedding of the values of Hilbert-valued weight into a Hilbert–Schmidt brick. Using Dorogovtsev's sufficient condition for the embedding of compact sets into a Hilbert–Schmidt brick in terms of an isonormal process, we prove the existence of Hilbert-valued Dynkin-renormalized self-intersection local time. Also using Dynkin's construction we construct the self-intersection local time for the deterministic image of the planar Wiener process.     

Long-time existence for the semilinear Klein–Gordon equation on a compact boundary-less Riemannian manifold

We investigate the long-time existence of small and smooth solutions for the semilinear Klein–Gordon equation on a compact boundary-less Riemannian manifold. Without any spectral or geometric assumption, our first result improves the lifespan obtained by the local theory. The previous result is proved under a generic condition of the mass. As a by-product of the method, we examine the particular case, where the manifold is a multidimensional torus, and we give explicit examples of algebraic masses for which we can improve the local existence time. The analytic part of the proof relies on multilinear estimates of eigenfunctions and estimates of small divisors proved by Delort–Szeftel. The algebraic part of the proof relies on a multilinear version of the Roth theorem proved by Schmidt.     

Blow-up criteria for the three dimensional nonlinear dissipative system modeling electro-hydrodynamics

In this paper, we investigate some sufficient conditions for the breakdown of local smooth solutions to the three dimensional nonlinear nonlocal dissipative system modeling electro-hydrodynamics. This model is a strongly coupled system by the well-known incompressible Navier–Stokes equations and the classical Poisson–Nernst–Planck equations. We show that the maximum of the vorticity field alone controls the breakdown of smooth solutions, which reveals that the velocity field plays a more dominant role than the density functions of charged particles in the blow-up theory of the system. Moreover, some Prodi–Serrin type blow-up criteria are also established.     

Stationary Solutions of Compressible Navier–Stokes Equations with Slip Boundary Conditions

We consider the Navier–Stokes equations for a compressible, viscous fluid with heat–conduction in a bounded domain of IR² or IR³. Under the assumption that the external force field and the external heat supply are small we prove the existence and local uniqueness of a stationary solution satisfying a slip boundary condition. For the temperature we assume a Dirichlet or an oblique boundary condition.     

A discrete-time model with non-monotonic functional response and strong Allee effect in prey

In this paper, we investigate the impact of strong Allee effect on the stability of a discrete-time predator–prey model with a non-monotonic functional response. The dynamics of discrete-time predator–prey models with strong Allee effect is studied earlier. But, the mathematical investigations of predator–prey dynamics in discrete-time set up with Holling type-IV functional response and strong Allee effect in prey are lacking. The proposed model supports the coexistence of two steady states, and the mathematical features of the model are analyzed based on local stability and bifurcation theory. By considering the Allee parameter as the bifurcation parameter, we provide sufficient conditions for the flip and the Neimark–Sacker bifurcations. We observe that Allee parameter plays a significant role in the dynamics of the system.     

New blow-up criteria for local solutions of the 3D generalized MHD equations in Lei–Lin–Gevrey spaces

This paper determines the existence of a unique local solution for the 3D generalized magnetohydrodynamics equations. In order to be more precise, our solution is obtained by involving Lei–Lin–Gevrey spaces as well as Lei–Lin spaces. Furthermore, we present five new blow-up criteria for this same system when the maximal time of existence is finite. It is important to point out that one of these criteria is obtained by assuming fractional Laplacians with equal powers.     

Discrete May–Leonard Competition Models III

A discretization method attributed to Kahan is used to approximate the May–Leonard (M–L) competition model for three species. The local dynamics of this discrete-time M–L model are analyzed. This model differs from the discrete M–L models being studied previously. This discrete model shows dynamical consistency with the continuous M–L model. Numerically, we showed that the discrete M–L model has a degenerate Hopf-bifurcation, which is consistent with the continuous M–L model.     

Optimal control problems for the reaction–diffusion–convection equation with variable coefficients

The solvability of optimal control problems is proved on both weak and strong solutions of a boundary value problem for the nonlinear reaction–diffusion–convection equation with variable coefficients. In the second case, the requirements for smoothness of the multiplicative control are reduced. The study of extremal problems is based on the proof of the solvability of the corresponding boundary value problems and on the qualitative analysis of their solutions properties. The large data existence results for weak solutions, the maximum principle as well as the local existence and uniqueness of a strong solution are established. Moreover, an optimal feedback control problem is considered. Using methods of the theory of topological degree for set-valued perturbations (with aspheric image sets) of generalized monotone operators, we obtain sufficient conditions for the solvability of this problem in the class of weak solutions.     

Probability density function of SDEs with unbounded and path-dependent drift coefficient

In this paper, we first prove that the existence of a solution of SDEs under the assumptions that the drift coefficient is of linear growth and path-dependent, and diffusion coefficient is bounded, uniformly elliptic and Hölder continuous. We apply Gaussian upper bound for a probability density function of a solution of SDE without drift coefficient and local Novikov condition, in order to use Maruyama–Girsanov transformation. The aim of this paper is to prove the existence with explicit representations (under linear/super-linear growth condition), Gaussian two-sided bound and Hölder continuity (under sub-linear growth condition) of a probability density function of a solution of SDEs with path-dependent drift coefficient. As an application of explicit representation, we provide the rate of convergence for an Euler–Maruyama (type) approximation, and an unbiased simulation scheme.     

Dunkl-Schrödinger semigroups and applications

We obtain dispersive estimates for the linear Dunkl–Schrödinger equations with and without quadratic potential. As a consequence, we prove the local well-posedness for semilinear Dunkl–Schrödinger equations with polynomial nonlinearity in certain magnetic field. Furthermore, we study many applications: as the uncertainty principles for the Dunkl transform via the Dunkl–Schrödinger semigroups, the embedding theorems for the Sobolev spaces associated with the generalized Hermite semigroup. Finally, almost every where convergence of the solutions of the Dunkl–Schrödinger equation is also considered.