Search

This paper is concerned with a two-species predator-prey reaction-diffusion system with Beddington-DeAngelis functional response and subject to homogeneous Neumann boundary conditions. By linearizing the system at the positive constant steady-state solution and analyzing the associated characteristic equation in detail, the asymptotic stability of the positive constant steady-state solution and the existence of local Hopf bifurcations are investigated. Also, it is shown that the appearance of the diffusion and homogeneous Neumann boundary conditions can lead to the appearance of codimension two Bagdanov-Takens bifurcation. Moreover, by applying the normal form theory and the center manifold reduction for partial differential equations (PDEs), the explicit algorithm determining the direction of Hopf bifurcations and the stability of bifurcating periodic solutions is given. Finally, numerical simulations supporting the theoretical analysis are also included. [ABSTRACT FROM AUTHOR]

This paper presents an efficient multiscale stochastic framework for uncertainty quantification in modeling of flow through porous media with multiple rock types. The governing equations are based on Darcy's law with nonstationary stochastic permeability represented as a sum of local Karhunen-Loève expansions. The approximation uses stochastic collocation on either a tensor product or a sparse grid, coupled with a domain decomposition algorithm known as the multiscale mortar mixed finite element method. The latter method requires solving a coarse scale mortar interface problem via an iterative procedure. The traditional implementation requires the solution of local fine scale linear systems on each iteration. We employ a recently developed modification of this method that precomputes a multiscale flux basis to avoid the need for subdomain solves on each iteration. In the stochastic setting, the basis is further reused over multiple realizations, leading to collocation algorithms that are more efficient than the traditional implementation by orders of magnitude. Error analysis and numerical experiments are presented. [ABSTRACT FROM AUTHOR]

Linear advection of a scalar quantity by a specified velocity field arises in a number of different applications. Of particular interest here is the transport of species and energy in low Mach number models for combustion, atmospheric flows, and astrophysics, as well as contaminant transport in Darcy models of saturated subsurface flow. An important characteristic of these problems is that the velocity field is not known analytically. Instead, an auxiliary equation is solved to compute averages of the velocities over faces in a finite volume discretization. In this paper, we present a customized three-dimensional finite volume advection scheme for this class of problems that provides accurate resolution for smooth problems while avoiding undershoot and overshoot for nonsmooth profiles. The method is an extension of an algorithm by Bell, Dawson, and Shubin (BDS), which was developed for a class of scalar conservation laws arising in porous media flows in two dimensions. The original BDS algorithm is a variant of unsplit. higher-order Godunov methods based on construction of a limited bilinear profile within each computational cell. Here we present a three-dimensional extension of the original BDS algorithm that is based on a limited trilinear profile within each cell. We compare this new method to several other unsplit approaches, including piecewise linear methods, piecewise parabolic methods, and wave propagation schemes. [ABSTRACT FROM AUTHOR]

In this paper, we construct a group of Saul'yev type asymmetric difference formulas for the dispersive equation. Based on these formulas we derive a new alternating 6-point group algorithm to solve dispersive equations with periodic boundary conditions. The algorithm has a high-order accuracy in space and an unconditional stability. The theoretical results are conformed to the numerical simulation. A comparison of this algorithm with the previous Alternating Group Explicit method is presented. [ABSTRACT FROM AUTHOR]