Search

We provide a quantitative lower bound to the Cheeger constant of a set $\Omega$ in both the Euclidean and the Gaussian settings in terms of suitable asymmetry indexes. We provide examples which show that these quantitative estimates are sharp., Comment: 18 pages, 3 figures

The definition of a Green’s function of a Cauchy–Dirichlet problem for the hyperbolic equation in a quarter plane is given. Its existence and uniqueness have been proven. Representation of the Green’s function is given. It is shown that the Green’s function can be represented by the Riemann–Green function.

We consider nonlinear reaction systems satisfying mass-action kinetics with slow and fast reactions. It is known that the fast-reaction-rate limit can be described by an ODE with Lagrange multipliers and a set of nonlinear constraints that ask the fast reactions to be in equilibrium. Our aim is to study the limiting gradient structure which is available if the reaction system satisfies the detailed-balance condition. The gradient structure on the set of concentration vectors is given in terms of the relative Boltzmann entropy and a cosh-type dissipation potential. We show that a limiting or effective gradient structure can be rigorously derived via EDP-convergence, i.e. convergence in the sense of the energy-dissipation principle for gradient flows. In general, the effective entropy will no longer be of Boltzmann type and the reactions will no longer satisfy mass-action kinetics. Deutsche Forschungsgemeinschafthttps://doi.org/10.13039/501100001659

In this paper, we study the inviscid and zero Froude number limits of the viscous shallow water system. We prove that the limit system is represented by the incompressible Euler equations on the whole space. Furthermore, the rate of convergence is also obtained.

Every non-convex pair $(C, D)$ ( C , D ) may not have proximal normal structure even in a Hilbert space. In this article, we use cyclic relatively nonexpansive maps with respect to orbits to show the presence of best proximity points in $C\cup D$ C ∪ D , where $C\cup D$ C ∪ D is a cyclic T-regular set and $(C, D)$ ( C , D ) is a non-empty, non-convex proximal pair in a real Hilbert space. Moreover, we show the presence of best proximity points and fixed points for non-cyclic relatively nonexpansive maps with respect to orbits defined on $C\cup D$ C ∪ D , where C and D are T-regular sets in a uniformly convex Banach space satisfying $T(C)\subseteq C$ T ( C ) ⊆ C , $T(D)\subseteq D$ T ( D ) ⊆ D wherein the convergence of Kranoselskii’s iteration process is also discussed.

In this paper, we propose a new iterative algorithm for solving the multiple-sets split feasibility problem (MSSFP for short) and the split equality fixed point problem (SEFPP for short) with firmly quasi-nonexpansive operators or nonexpansive operators in real Hilbert spaces. Under mild conditions, we prove strong convergence theorems for the algorithm by using the projection method and the properties of projection operators. The result improves and extends the corresponding ones announced by some others in the earlier and recent literature.

In this work, by introducing several parameters, a new kernel function including both the homogeneous and non-homogeneous cases is constructed, and a Hilbert-type inequality related to the newly constructed kernel function is established. By convention, the equivalent Hardy-type inequality is also considered. Furthermore, by introducing the partial fraction expansions of trigonometric functions, some special and interesting Hilbert-type inequalities with the constant factors represented by the higher derivatives of trigonometric functions, the Euler number and the Bernoulli number are presented at the end of the paper.

Apartial discrete Dirichlet boundary value problem involving mean curvature operator is concerned in this paper. Under proper assumptions on the nonlinear term, we obtain some feasible conditions on the existence of multiple solutions by the method of critical point theory. We further separately determine open intervals of the parameter to attain at least two positive solutions and an unbounded sequence of positive solutions with the help of the maximum principle.

Near-optimization is as sensible and important as optimization for both theory and applications. This paper concerns the near-optimal control of an avian influenza model with saturation on heterogeneous complex networks. Firstly, the basic reproduction number $ \mathcal{R}_{0} $ is defined for the model, which can be used to govern the threshold dynamics of influenza disease. Secondly, the near-optimal control problem was formulated by slaughtering poultry and treating infected humans while keeping the loss and cost to a minimum. Thanks to the maximum condition of the Hamiltonian function and the Ekeland's variational principle, we establish both necessary and sufficient conditions for the near-optimality by several delicate estimates for the state and adjoint processes. Finally, a number of examples presented to illustrate our theoretical results.

In this contribution, we investigate an initial-boundary value problem for a fractional diffusion equation with Caputo fractional derivative of space-dependent variable order where the coefficients are dependent on spatial and time variables. We consider a bounded Lipschitz domain and a homogeneous Dirichlet boundary condition. Variable-order fractional differential operators originate in anomalous diffusion modelling. Using the strongly positive definiteness of the governing kernel, we establish the existence of a unique weak solution in $u\in \operatorname{L}^{\infty } ((0,T),\operatorname{H}^{1}_{0}( \Omega ) )$ u ∈ L ∞ ( ( 0 , T ) , H 0 1 ( Ω ) ) to the problem if the initial data belongs to $\operatorname{H}^{1}_{0}(\Omega )$ H 0 1 ( Ω ) . We show that the solution belongs to $\operatorname{C} ([0,T],{\operatorname{H}^{1}_{0}(\Omega )}^{*} )$ C ( [ 0 , T ] , H 0 1 ( Ω ) ∗ ) in the case of a Caputo fractional derivative of constant order. We generalise a fundamental identity for integro-differential operators of the form $\frac{\mathrm{d}}{\mathrm{d}t} (k\ast v)(t)$ d d t ( k ∗ v ) ( t ) to a convolution kernel that is also space-dependent and employ this result when searching for more regular solutions. We also discuss the situation that the domain consists of separated subdomains.

In this paper, we study the fractional Schrödinger-Poisson system(−Δ)su+V(x)u+K(x)ϕ|u|q−2u=h(x)f(u)+|u|2s∗−2u,in R3,(−Δ)tϕ=K(x)|u|q,in R3,$$\begin{array}{} \displaystyle \left\{ \begin{array}{ll} (-{\it\Delta})^{s}u+V(x)u+ K(x) \phi|u|^{q-2}u=h(x)f(u)+|u|^{2^{\ast}_{s}-2}u,&\mbox{in}~ {\mathbb R^{3}},\\ (-{\it\Delta})^{t}\phi=K(x)|u|^{q},&\mbox{in}~ {\mathbb R^{3}}, \end{array}\right. \end{array}$$wheres,t∈ (0, 1), 3 < 4s< 3 + 2t,q∈ (1,2s∗$\begin{array}{} \displaystyle 2^*_s \end{array}$/2) are real numbers, (−Δ)sstands for the fractional Laplacian operator,2s∗:=63−2s$\begin{array}{} \displaystyle 2^{*}_{s}:=\frac{6}{3-2s} \end{array}$is the fractional critical Sobolev exponent,K,Vandhare non-negative potentials andV,hmay be vanish at infinity.fis aC1-function satisfying suitable growth assumptions. We show that the above fractional Schrödinger-Poisson system has a positive and a sign-changing least energy solution via variational methods.

This paper studies the mass-critical variable coefficient nonlinear Schrödinger equation. We first get the existence of the ground state by solving a minimization problem. Then we prove a compactness result by the variational characterization of the ground state solutions. In addition, we construct the blow-up solutions at the minimal mass threshold and further prove the uniqueness result on the minimal mass blow-up solutions which are pseudo-conformal transformation of the ground states.

In this paper, our aim is mathematical analysis and numerical simulation of a prey-predator model to describe the effect of predation between prey and predator with nonlinear functional response. First, we develop results concerning the boundedness, the existence and uniqueness of the solution. Furthermore, the Lyapunov principle and the Routh–Hurwitz criterion are applied to study respectively the local and global stability results. We also establish the Hopf-bifurcation to show the existence of a branch of nontrivial periodic solutions. Finally, numerical simulations have been accomplished to validate our analytical findings.

In this paper, we consider the maximal operator related to the Laplace-Bessel differential operator ( B B -maximal operator) on L p ( ⋅ ) , γ ( R k , + n ) {L}_{p\left(\cdot ),\gamma }\left({{\mathbb{R}}}_{k,+}^{n}) variable exponent Lebesgue spaces. We will give a necessary condition for the boundedness of the B B -maximal operator on variable exponent Lebesgue spaces. Moreover, we will obtain that the B B -maximal operator is not bounded on L p ( ⋅ ) , γ ( R k , + n ) {L}_{p\left(\cdot ),\gamma }\left({{\mathbb{R}}}_{k,+}^{n}) variable exponent Lebesgue spaces in the case of p − = 1 {p}_{-}=1 . We will also prove the boundedness of the fractional maximal function associated with the Laplace-Bessel differential operator (fractional B B -maximal function) on L p ( ⋅ ) , γ ( R k , + n ) {L}_{p\left(\cdot ),\gamma }\left({{\mathbb{R}}}_{k,+}^{n}) variable exponent Lebesgue spaces.

In this paper, we extend existing population growth models and propose a model based on a nonlinear cubic differential equation that reveals itself as a special subclass of Abel differential equations of first kind. We first summarize properties of the time-continuous problem formulation. We state the boundedness, global existence, and uniqueness of solutions for all times. Proofs of these properties are thoroughly given in the Appendix to this paper. Subsequently, we develop an explicit–implicit time-discrete numerical solution algorithm for our time-continuous population growth model and show that many properties of the time-continuous case transfer to our numerical explicit–implicit time-discrete solution scheme. We provide numerical examples to illustrate different behaviors of our proposed model. Furthermore, we compare our explicit–implicit discretization scheme to the classical Eulerian discretization. The latter violates the nonnegativity constraints on population sizes, whereas we prove and illustrate that our explicit–implicit discretization algorithm preserves this constraint. Finally, we describe a parameter estimation approach to apply our algorithm to two different real-world data sets.

Many mathematicians have studied degenerate versions of quite a few special polynomials and numbers since Carlitz’s work (Utilitas Math. 15 (1979), 51–88). Recently, Kim et al. studied the degenerate gamma random variables, discrete degenerate random variables and two-variable degenerate Bell polynomials associated with Poisson degenerate central moments, etc. This paper is divided into two parts. In the first part, we introduce a new type of degenerate Bell polynomials associated with degenerate Poisson random variables with parameter α > 0 \alpha \hspace{-0.15em}\gt \hspace{-0.15em}0 , called the fully degenerate Bell polynomials. We derive some combinatorial identities for the fully degenerate Bell polynomials related to the n n th moment of the degenerate Poisson random variable, special numbers and polynomials. In the second part, we consider the fully degenerate Bell polynomials associated with degenerate Poisson random variables with two parameters α > 0 \alpha \gt 0 and β > 0 \beta \hspace{-0.15em}\gt \hspace{-0.15em}0 , called the two-variable fully degenerate Bell polynomials. We show their connection with the degenerate Poisson central moments, special numbers and polynomials.

The nontrivial behavior of wave packets in the space fractional coupled nonlinear Schrödinger equation has received considerable theoretical attention. The difficulty comes from the fact that the Riesz fractional derivative is inherently a prehistorical operator. In contrast, nonlinear Schrödinger equation with both time and space nonlocal operators, which is the cornerstone in the modeling of a new type of fractional quantum couplers, is still in high demand of attention. This paper is devoted to numerically study the propagation of solitons through a new type of quantum couplers which can be called time-space fractional quantum couplers. The numerical methodology is based on the finite-difference/Galerkin Legendre spectral method with an easy to implement numerical algorithm. The time-fractional derivative is considered to describe the decay behavior and the nonlocal memory of the model. We conduct numerical simulations to observe the performance of the tunable decay and the sharpness behavior of the time-space fractional strongly coupled nonlinear Schrödinger model as well as the performance of the numerical algorithm. Numerical simulations show that the time and space fractional-order operators control the decay behavior or the memory and the sharpness of the interface and undergo a seamless transition of the fractional-order parameters. © 2021, The Author(s). This study was supported financially by RFBR Grant (19-01-00019), the National Research Centre of Egypt (NRC) and Ghent university.

Our aim in this paper is presenting an attractive numerical approach giving an accurate solution to the nonlinear fractional Abel differential equation based on a reproducing kernel algorithm with model endowed with a Caputo–Fabrizio fractional derivative. By means of such an approach, we utilize the Gram–Schmidt orthogonalization process to create an orthonormal set of bases that leads to an appropriate solution in the Hilbert space $\mathcal{H}^{2}[a,b]$ H 2 [ a , b ] . We investigate and discuss stability and convergence of the proposed method. The n-term series solution converges uniformly to the analytic solution. We present several numerical examples of potential interests to illustrate the reliability, efficacy, and performance of the method under the influence of the Caputo–Fabrizio derivative. The gained results have shown superiority of the reproducing kernel algorithm and its infinite accuracy with a least time and efforts in solving the fractional Abel-type model. Therefore, in this direction, the proposed algorithm is an alternative and systematic tool for analyzing the behavior of many nonlinear temporal fractional differential equations emerging in the fields of engineering, physics, and sciences.

We consider a class of semilinear equations with an absorption nonlinear zero order term of power type, where elliptic condition is given in terms of Gauss measure. In the case of the superlinear equation we introduce a suitable definition of solutions in order to prove the existence and uniqueness of a solution in ℝ N without growth restrictions at infinity. A comparison result in terms of the half-space Gaussian symmetrized problem is also proved. As an application, we give some estimates in measure of the growth of the solution near the boundary of its support for sublinear equations. Finally we generalize our results to problems with a nonlinear zero order term not necessary of power type.

Mass-media coverage is one of the most widely used government strategies on influencing public opinion in times of crisis. Awareness campaigns are highly influential tools to expand healthy behavior practices among individuals during epidemics and pandemics. Mathematical modeling has become an important tool in analyzing the effects of media awareness on the spread of infectious diseases. In this paper, a fractional-order epidemic model incorporating media coverage is presented and analyzed. The problem is formulated using susceptible, infectious and recovered compartmental model. A long-term memory effect modeled by a Caputo fractional derivative is included in each compartment to describe the evolution related to the individuals’ experiences. The well-posedness of the model is investigated in terms of global existence, positivity, and boundedness of solutions. Moreover, the disease-free equilibrium and the endemic equilibrium points are given alongside their local stabilities. By constructing suitable Lyapunov functions, the global stability of the disease-free and endemic equilibria is proven according to the basic reproduction number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}R0. Finally, numerical simulations are performed to support our analytical findings. It was found out that the long-term memory has no effect on the stability of the equilibrium points. However, for increased values of the fractional derivative order parameter, each solution reaches its equilibrium state more rapidly. Furthermore, it was observed that an increase of the media awareness parameter, decreases the magnitude of infected individuals, and consequently, the height of the epidemic peak.

Asymptotic stability analysis applied to stabilize unstable fixed points and to control chaotic motions in two and three-dimensional discrete dynamical systems. A new set of parameter values obtained which stabilizes an unstable fixed point and control the chaotic evolution to regularity. The output of the considered model and that of the adjustable system continuously compared by a typical feedback and the difference used by the adaptation mechanism to modify the parameters. Suitable numerical simulation which are used thoroughly discussed and parameter values are adjusted. The findings are significant and interesting. This strategy has some advantages over many other chaos control methods in discrete systems but, however it can be applied within some limitations.

A diffusive phytoplankton–zooplankton model with nonlinear harvesting is considered in this paper. Firstly, using the harvesting as the parameter, we get the existence and stability of the positive steady state, and also investigate the existence of spatially homogeneous and inhomogeneous periodic solutions. Then, by applying the normal form theory and center manifold theorem, we give the stability and direction of Hopf bifurcation from the positive steady state. In addition, we also prove the existence of the Bogdanov–Takens bifurcation. These results reveal that the harvesting and diffusion really affect the spatiotemporal complexity of the system. Finally, numerical simulations are also given to support our theoretical analysis.

This article proposes new strategies for solving two-point Fractional order Nonlinear Boundary Value Problems (FNBVPs) with Robin Boundary Conditions (RBCs). In the new numerical schemes, a two-point FNBVP is transformed into a system of Fractional order Initial Value Problems (FIVPs) with unknown Initial Conditions (ICs). To approximate ICs in the system of FIVPs, we develop nonlinear shooting methods based on Newton’s method and Halley’s method using the RBC at the right end point. To deal with FIVPs in a system, we mainly employ High-order Predictor–Corrector Methods (HPCMs) with linear interpolation and quadratic interpolation (Nguyen and Jang in Fract. Calc. Appl. Anal. 20(2):447–476, 2017) into Volterra integral equations which are equivalent to FIVPs. The advantage of the proposed schemes with HPCMs is that even though they are designed for solving two-point FNBVPs, they can handle both linear and nonlinear two-point Fractional order Boundary Value Problems (FBVPs) with RBCs and have uniform convergence rates of HPCMs, $\mathcal{O}(h^{2})$ O ( h 2 ) and $\mathcal{O}(h^{3})$ O ( h 3 ) for shooting techniques with Newton’s method and Halley’s method, respectively. A variety of numerical examples are demonstrated to confirm the effectiveness and performance of the proposed schemes. Also we compare the accuracy and performance of our schemes with another method.

In this work, we consider a fractional diffusion equation with nonlocal integral condition. We give a form of the mild solution under the expression of Fourier series which contains some Mittag-Leffler functions. We present two new results. Firstly, we show the well-posedness and regularity for our problem. Secondly, we show the ill-posedness of our problem in the sense of Hadamard. Using the Fourier truncation method, we construct a regularized solution and present the convergence rate between the regularized and exact solutions.

In this paper, a stochastic SIRV epidemic model with general nonlinear incidence and vaccination is investigated. The value of our study lies in two aspects. Mathematically, with the help of Lyapunov function method and stochastic analysis theory, we obtain a stochastic threshold of the model that completely determines the extinction and persistence of the epidemic. Epidemiologically, we find that random fluctuations can suppress disease outbreak, which can provide us some useful control strategies to regulate disease dynamics. In other words, neglecting random perturbations overestimates the ability of the disease to spread. The numerical simulations are given to illustrate the main theoretical results.

Algorithms for studying transitions and instabilities in incompressible flows typically require the solution of linear systems with the full Jacobian matrix. Other popular approaches, like gradient-based design optimization and fully implicit time integration, also require very robust solvers for this type of linear system. We present a parallel fully coupled multilevel incomplete factorization preconditioner for the 3D stationary incompressible Navier-Stokes equations on a structured grid. The algorithm and software are based on the robust two-level method developed by Wubs and Thies. In this article, we identify some of the weak spots of the two-level scheme and propose remedies such as a different domain partitioning and recursive application of the method. We apply the method to the well-known 3D lid-driven cavity benchmark problem, and demonstrate its superior robustness by comparing with a segregated SIMPLE-type preconditioner.

In this paper, exact solutions to the problem of acoustic scattering by elastic spherical symmetric scatterers are developed. The scatterer may consist of an arbitrary number of fluid and solid layers, and scattering with single Neumann conditions (replacing Neumann-to-Neumann conditions) is added. The solution is obtained by separation of variables, resulting in an infinite series which must be truncated for numerical evaluation. The implemented numerical solution is exact in the sense that numerical error is solely due to round-off errors, which will be shown using the symbolic toolbox in MATLAB. A system of benchmark problems is proposed for future reference. Numerical examples are presented, including comparisons with reference solutions, far-field patterns and near-field plots of the benchmark problems, and time-dependent solutions obtained by Fourier transformation.

In this article, the sine-Gordon expansion method is employed to find some new traveling wave solutions to the nonlinear Schrödinger equation with the coefficients of both group velocity dispersion and second-order spatiotemporal dispersion. The nonlinear model is reduced to an ordinary differential equation by introducing an intelligible wave transformation. A set of new exact solutions are observed corresponding to various parameters. These novel soliton solutions are depicted in figures, revealing the new physical behavior of the acquired solutions. The method proves its ability to provide good new approximate solutions with some applications in science. Moreover, the associated solution of the presented method can be extended to solve more complex models.

Evaluating the action of a matrix function on a vector, that is $$x=f({\mathcal {M}})v$$ x = f ( M ) v , is an ubiquitous task in applications. When $${\mathcal {M}}$$ M is large, one usually relies on Krylov projection methods. In this paper, we provide effective choices for the poles of the rational Krylov method for approximating x when f(z) is either Cauchy–Stieltjes or Laplace–Stieltjes (or, which is equivalent, completely monotonic) and $${\mathcal {M}}$$ M is a positive definite matrix. Relying on the same tools used to analyze the generic situation, we then focus on the case $${\mathcal {M}}=I \otimes A - B^T \otimes I$$ M = I ⊗ A - B T ⊗ I , and v obtained vectorizing a low-rank matrix; this finds application, for instance, in solving fractional diffusion equation on two-dimensional tensor grids. We see how to leverage tensorized Krylov subspaces to exploit the Kronecker structure and we introduce an error analysis for the numerical approximation of x. Pole selection strategies with explicit convergence bounds are given also in this case.

In a recently published theorem on the split common fixed point problem for strict pseudocontractive and asymptotically nonexpansive mappings, Tang et al. (J. Inequal. Appl. 2015:305, 2015) studied a uniformly convex and 2-uniformly smooth real Banach space with the Opial property and best smoothness constant κ satisfying the condition $0 0 < κ < 1 2 , as a real Banach space more general than Hilbert spaces. A well-known example of a uniformly convex and 2-uniformly smooth real Banach space with the Opial property is $E=l_{p}$ E = l p , $2\leq p 2 ≤ p < ∞ . It is shown in this paper that, if κ is the best smoothness constant of E and satisfies the condition $0 0 < κ ≤ 1 2 , then E is necessarily $l_{2}$ l 2 , a real Hilbert space. Furthermore, some important remarks concerning the proof of this theorem are presented.

In this paper, we develop an iterative algorithm whose architecture comprises a modified version of the forward–backward splitting algorithm and the hybrid shrinking projection algorithm. We provide theoretical results concerning weak and strong convergence of the proposed algorithm towards a common solution of the fixed point problem associated to a finite family of demicontractive operators, the split equilibrium problem and the monotone inclusion problem in Hilbert spaces. Moreover, we compute a numerical experiment to show the efficiency of the proposed algorithm. As a consequence, our results improve various existing results in the current literature.

In 2020, El-Deeb et al. proved several dynamic inequalities. It is our aim in this paper to give the retarded time scales case of these inequalities. We also give a new proof and formula of Leibniz integral rule on time scales. Beside that, we also apply our inequalities to discrete and continuous calculus to obtain some new inequalities as special cases. Furthermore, we study boundedness of some delay initial value problems by applying our results as application.

By introducing a kernel involving an exponent function with multiple parameters, we establish a new Hilbert-type inequality and its equivalent Hardy form. We also prove that the constant factors of the obtained inequalities are the best possible. Furthermore, by introducing the Bernoulli number, Euler number, and the partial fraction expansion of cotangent function and cosecant function, we get some special and interesting cases of the newly obtained inequality.

In this paper, we investigate the analysis of the proportional Caputo derivative that recently has been constructed. We create some useful relations between this new derivative and beta function. We discretize the new derivative. We investigate the stability and obtain a stability condition for the new derivative.

Various collective motions of camphor boats have been studied. Camphor boats are self-driven particles that interact with each other through the change of surface tension on water surface by camphor molecules. Consequently, even in a one-dimensional circuit, a congested state or jamming can be observed (Suematsu et al. in Phys Rev E 81:056210, 2010). In this phenomenon, the unidirectional motion of each particle is considered a traveling wave, and the concentration of camphor molecules forms a pulse shape. Hence, each pair of particles interacts with each other like a pulse–pulse interaction. Thus, we expect that the center manifold theories proposed in Ei (J Dyn Differ Equ 14:85–137, 2002) and Ei et al. (Physica D 165:176–198 2002) are applicable for the analysis of the collective motion of camphor boats. However, spatial discontinuity in our model, in particular the existence of Dirac delta functions in a linearized operator, does not fulfill the requirement in the reduction process because the authors developed their theories in $$L^2$$ -framework and for smooth nonlinearity in Ei (2002) and Ei et al. (2002). In this article, we extend the results obtained in Ei (2002) and Ei et al. (2002) and establish a new center manifold approach in $$(H^1)^*$$ -framework. Finally, we succeed to rigorously reduce a mathematical model (Nagayama et al. in Physica D: Nonlinear Phenom 194:151–165, 2004) coupled with a Newtonian equation and a reaction–diffusion equation including delta functions to an ordinary differential system that represents the motions of camphor boats.

This paper presents a numerical scheme of arbitrary order of accuracy in both space and time, based on the arbitrary high-order derivatives methodology, for transient acoustic simulations. The scheme combines the nodal discontinuous Galerkin method for the spatial discretization and the Taylor series integrator (TSI) for the time integration. The main idea of the TSI is a temporal Taylor series expansion of all unknown acoustic variables in which the time derivatives are replaced by spatial derivatives via the Cauchy-Kovalewski procedure. The computational cost for the time integration is linearly proportional to the order of accuracy. To increase the computational efficiency for simulations involving strongly varying mesh sizes or material properties, a local time-stepping (LTS) algorithm accompanying the arbitrary high-order derivatives discontinuous Galerkin (ADER-DG) scheme, which ensures correct communications between domains with different time step sizes, is proposed. A numerical stability analysis in terms of the maximum allowable time step sizes is performed. Based on numerical convergence analysis, we demonstrate that for nonuniform meshes, a consistent high-order accuracy in space and time is achieved using ADER-DG with LTS. An application to the sound propagation across a transmissive noise barrier validates the potential of the proposed method for practical problems demanding high accuracy.

This paper is concerned with a stochastic predator-prey model with Holling II increasing function in the predator. By applying the Lyapunov analysis method, we demonstrate the existence and uniqueness of the global positive solution. Then we show there is a stationary distribution which implies the stochastic persistence of the predator and prey in the model. Moreover, we obtain respectively sufficient conditions for weak persistence in the mean and extinction of the prey and extinction of the predator. Finally, some numerical simulations are given to illustrate our main results and the discussion and conclusion are presented.

In order to completely separate objects with large sections of occluded boundaries in an image, we devise a new variational level set model for image segmentation combining the Chan-Vese model with elastica and landmark constraints. For computational efficiency, we design its Augmented Lagrangian Method (ALM) or Alternating Direction Method of Multiplier (ADMM) method by introducing some auxiliary variables, Lagrange multipliers, and penalty parameters. In each loop of alternating iterative optimization, the sub-problems of minimization can be easily solved via the Gauss-Seidel iterative method and generalized soft thresholding formulas with projection, respectively. Numerical experiments show that the proposed model can not only recover larger broken boundaries but can also improve segmentation efficiency, as well as decrease the dependence of segmentation on parameter tuning and initialization., 17pages,9figure

In this paper, we investigate a delay reaction–diffusion–advection model with ideal free dispersal. The stability of positive steady-state solutions and the existence of the associated Hopf bifurcation are obtained by analyzing the principal eigenvalue of an elliptic operator. By the normal form theory and the center manifold reduction, the stability and bifurcation direction of Hopf bifurcating periodic solutions are obtained. Moreover, numerical simulations and a brief discussion are presented to illustrate our theoretical results.

We mainly consider the limit behaviors of the Riemann solutions to Chaplygin Euler equations for nonisentropic fluids. The formation of delta shock wave and the appearance of vacuum state are found as parameter ε \varepsilon tends to a certain value. Different from the isentropic fluids, the weight of delta shock wave is determined by variance density ρ \rho and internal energy H. Meanwhile, involving the entropy inequality, the uniqueness of delta shock wave is obtained.

In this paper, we give two Ramanujan-type circular summation formulas by applying the way of elliptic functions and the properties of theta functions. As applications, we obtain the corresponding imaginary transformation formulas for Ramanujan-type circular summations and some theta function identities.

We study weak solutions of the two-dimensional (2D) filtered Euler equations for measure-valued initial vorticity. The filtered Euler equations are a regularized model based on a spatial filtering to the Euler equations and have a unique global weak solution for initial vorticity in the space of a finite Radon measure. In this paper, we treat the vortex sheet problem and the point-vortex problem on the 2D filtered Euler equations. We prove that vortex sheet solutions of the 2D filtered Euler equations converge to those of the 2D Euler equations in the limit of the filtering parameter when initial vortex sheet has a distinguished sign. In addition, we make it clear what kind of condition should be imposed on the spatial filter to show the convergence and, according to the condition, our result is applicable to well-known regularized models such as the Euler- $$\alpha $$ model and the vortex blob model. We also show that filtered point vortices converge to a solution of the point-vortex system provided that collision of point vortices does not occur.

In this paper, we deduce the generalized $q$-difference equations for general Al-Salam--Carlitz polynomials and generalize Arjika's recently results [$q$-difference equation for homogeneous $q$-difference operators and their applications, J. Differ. Equ. Appl. {\bf 26}, 987--999 (2020)]. In addition, we obtain transformational identities by the method of $q$-difference equation. Moreover, we deduce $U(n+1)$ type generating functions and Ramanujan's integrals involving general Al-Salam--Carlitz polynomials by $q$-difference equation., 14 pages

Identities of symmetry in two variables for Bernoulli polynomials and power sums had been investigated by considering suitable symmetric identities. T. Kim used a completely different tool, namely thep-adic Volkenborn integrals, to find the same identities of symmetry in two variables. Not much later, it was observed that thisp-adic approach can be generalized to the case of three variables and shown that it gives some new identities of symmetry even in the case of two variables upon specializing one of the three variables. In this paper, we generalize the results in three variables to those in an arbitrary number of variables in a suitable setting and illustrate our results with some examples.

In this paper, we are concerned with the global asymptotic behavior of an infection age-structured SIR epidemic model with diffusion in a general n-dimensional bounded spatial domain under the homogeneous Dirichlet boundary condition. By using the method of characteristics, we reformulate the model into a system of a reaction-diffusion equation and a Volterra integral equation. We define the basic reproduction number R 0 by the spectral radius of a compact positive linear operator and show that if R 0 1 , then the disease-free steady state is globally attractive, whereas if R 0 > 1 , then a positive endemic steady state exists and the system is uniformly persistent. By numerical simulation for the 2-dimensional case, we show that R 0 depends on the shape of the spatial domain. This result is in contrast with the case of the homogeneous Neumann boundary condition, in which R 0 is independent of the spatial domain.

We prove the existence of tripled fixed points (TFPs) of a new generalized nonlinear contraction mapping in complete cone b-metric spaces (CCbMSs). Also, we present some exciting consequences as corollaries and three nontrivial examples. Finally, we find a solution for a tripled-system of integral equations (TSIE) and discussed a unique stationary distribution for the Markov process (SDMP).

In this paper we establish an invariant set bifurcation theory for the nonautonomous dynamical system $(\varphi_\lambda,\theta)_{X,\mathcal H}$ generated by the evolution equation \begin{equation}\label{e0}u_t+Au=\lambda u+p(t,u),\hspace{0.4cm} p\in \mathcal H=\mathcal H[f(\cdot,u)]\end{equation} on a Hilbert space $X$, where $A$ is a sectorial operator, $\lambda$ is the bifurcation parameter, $f(\cdot,u):\mathbb{R}\rightarrow X$ is translation compact, $f(t,0)\equiv0$ and $\mathcal H[f]$ is the hull of $f(\cdot,u)$. Denote by $\varphi_\lambda:=\varphi_\lambda(t,p)u$ the cocycle semiflow generated by the equation. Under some other assumptions on $f$, we show that as the parameter $\lambda$ crosses an eigenvalue $\lambda_0\in\mathbb{R}$ of $A$, the system bifurcates from $0$ to a nonautonomous invariant set $B_\lambda(\cdot)$ on one-sided neighborhood of $\lambda_0$. Moreover, $$\lim_{\lambda\rightarrow\lambda_0}H_{X^\alpha}\left(B_\lambda(p),0\right)=0,\hspace{0.4cm} p\in P,$$ where $H_{X^\alpha}(\cdot,\cdot)$ denotes the Hausdorff semidistance in $X^\alpha$ (here $X^\alpha$ ($\alpha\geq0$) defined below is the fractional power spaces associated with $A$). Our result is based on the pullback attractor bifurcation on the local central invariant manifolds $\mathcal {M}^\lambda_{\operatorname{loc}}(\cdot)$.

In this paper, we study a kind of generalized Riemann–Hilbert problems (R-HPs) with several unknown functions in strip domains. We mainly discuss methods of solving R-HPs with two unknown functions and obtain general solutions and conditions of solvability on function spaces $H_{\ast }$ H ∗ . At the end of this paper, we consider in detail the behavior of the solution at ∞ and in different domains. Thus the results in this paper generalize and improve the theory of the classical Riemann–Hilbert problems.

In this study, random linear difference equations obtained by transforming the components of deterministic difference equations to random variables are investigated. Uniform, Bernoulli, binomial, negative binomial (or Pascal), geometric, hypergeometric and Poisson distributions have been used for the random effects for obtaining the random behavior of linear difference equations. The random version of the Z-transform, the RZ-transform, has been used to obtain an approximation for the random linear difference equation. Approximate expected values and variances are calculated by using the RZ-transform. The results have been obtained with Maple and are shown in graphs. It is shown that the random Z-transform is an effective tool for the investigation of random linear difference equations. Gumu shane UniversityGumushane University We would like to thank Gumu shane University for its support in this study. WOS:000581693400003 2-s2.0-85092569388

In this paper, we have investigated thermally stratified MHD flow of an Oldroyd-B fluid over an inclined stretching surface in the presence of heat generation/absorption. Similarity solutions for the transformed governing equations are obtained. The reduced equations are solved numerically using the Runge–Kutta Fehlberg method with shooting technique. The influences of various involved parameters on velocity profiles, temperature profiles, local skin friction, and local Nusselt number are discussed. Numerical values of local skin friction and local Nusselt number are computed. The significant outcomes of the study are that the velocity decreases when the radiation parameter $$R_{d}$$ R d is increased while the temperature profile is increased for higher values of radiation parameter $$R_{d}$$ R d in case of opposing flow, moreover, growth in Deborah number $$\beta_{2}$$ β 2 enhance the velocity and momentum boundary layer. The heat transfer rate is decrease due to magnetic strength but increase with the increased values of Prandtl and Deborah numbers. The results of this model are closely matched with the outputs available in the literature.