Uniform asymptotics of the boundary values of the solution in a linear problem on the run-up of waves on a shallow beach

We consider the Cauchy problem with spatially localized initial data for a two-dimensional wave equation with variable velocity in a domain Ω. The velocity is assumed to degenerate on the boundary ?Ω of the domain as the square root of the distance to ?Ω. In particular, this problems describes the run-up of tsunami waves on a shallow beach in the linear approximation. Further, the problem contains a natural small parameter (the typical source-to-basin size ratio) and hence admits analysis by asymptotic methods. It was shown in the paper “Characteristics with singularities and the boundary values of the asymptotic solution of the Cauchy problem for a degenerate wave equation” [1] that the boundary values of the asymptotic solution of this problem given by a modified Maslov canonical operator on the Lagrangian manifold formed by the nonstandard characteristics associatedwith the problemcan be expressed via the canonical operator on a Lagrangian submanifold of the cotangent bundle of the boundary. However, the problem as to how this restriction is related to the boundary values of the exact solution of the problem remained open. In the present paper, we show that if the initial perturbation is specified by a function rapidly decaying at infinity, then the restriction of such an asymptotic solution to the boundary gives the asymptotics of the boundary values of the exact solution in the uniform norm. To this end, we in particular prove a trace theorem for nonstandard Sobolev type spaces with degeneration at the boundary.     

A non symmetric FVM-BEM coupling method

We consider the prototype model for flow and transport of a concentration in porous media in an interior domain and couple it with a diffusion process in the corresponding unbounded exterior domain. To guarantee mass conservation and stability with respect to dominating convection also for a discrete solution we introduce a non symmetric coupling of the vertex-centered finite volume method (FVM) and the boundary element method (BEM). BEM approximates the unbounded exterior problem which avoids truncation of the domain. One can also interpret that the (unbounded) exterior problem “replaces” the boundary conditions of the interior problem. We aim to provide a first rigorous analysis of the discrete system for different model parameters; existence and uniqueness, convergence, and a priori estimates. Numerical examples illustrate the strength of the chosen method which is computational cheaper than the previous three field FVM-BEM couplings. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A nonlinearity lagging for the solution of nonlinear steady state reaction diffusion problems

This paper concerns with the analysis of the iterative procedure for the solution of a nonlinear reaction diffusion equation at the steady state in a two dimensional bounded domain supplemented by suitable boundary conditions. This procedure, called Lagged Diffusivity Functional Iteration (LDFI)-procedure, computes the solution by “lagging” the diffusion term. A model problem is considered and a finite difference discretization for that model problem is described. Furthermore, properties of the finite difference operator are proved. Then, sufficient conditions for the convergence of the LDFI-procedure are given. At each stage of the LDFI-procedure a weakly nonlinear algebraic system has to be solved and the simplified Newton–Arithmetic Mean (Newton–AM) method is used. This method is particularly well suited for implementation on parallel computers. Numerical studies show the efficiency, for different test functions, of the LDFI-procedure combined with the simplified Newton–AM method. Better results are obtained when in the reaction diffusion equation also a convection term is present.     

在部分区域上的奇摄动反应扩散方程初始边值问题解的渐近性态

讨论一类在部分区域上的奇摄动反应扩散方程初始边值问题,利用算子理论,得到了相应问题解的渐近性态。     

Global asymptotic stability of Lotka–Volterra competition reaction–diffusion systems with time delays

This paper is concerned with a time-delayed Lotka–Volterra competition reaction–diffusion system with homogeneous Neumann boundary conditions. Some explicit and easily verifiable conditions are obtained for the global asymptotic stability of all forms of nonnegative semitrivial constant steady-state solutions. These conditions involve only the competing rate constants and are independent of the diffusion–convection and time delays. The result of global asymptotic stability implies the nonexistence of positive steady-state solutions, and gives some extinction results of the competing species in the ecological sense. The instability of the trivial steady-state solution is also shown.     

Convergence analysis of an hp finite element method for singularly perturbed transmission problems in smooth domains

We consider a two‐dimensional singularly perturbed transmission problem with two different diffusion coefficients, in a domain with smooth (analytic) boundary. The solution will contain boundary layers only in the part of the domain where the diffusion coefficient is high and interface layers along the interface. Utilizing existing and newly derived regularity results for the exact solution, we prove the robustness of an hp finite element method for its approximation. Under the assumption of analytic input data, we show that the method converges at an “exponential” rate, provided the mesh and polynomial degree distribution are chosen appropriately. Numerical results illustrating our theoretical findings are also included. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Singularly perturbed boundary value problems for a class of second order turning point on infinite interval

This article solved the asymptotic solution of a singularly perturbed boundary value problem with second order turning point, encountered in the dissipative equilibrium vector field of the coupled convection disturbance kinetic equations under the constrained filed and the gravity. Using the matching of asymptotic expansions, the formal asymptotic solution is constructed. By using the theory of differential inequality the uniform validity of the asymptotic expansion for the solution is proved.     

On one model problem for the reaction–diffusion–advection equation

The asymptotic behavior of the solution with boundary layers in the time-independent mathematical model of reaction–diffusion–advection arising when describing the distribution of greenhouse gases in the surface atmospheric layer is studied. On the basis of the asymptotic method of differential inequalities, the existence of a boundary-layer solution and its asymptotic Lyapunov stability as a steady-state solution of the corresponding parabolic problem is proven. One of the results of this work is the determination of the local domain of the attraction of a boundary-layer solution.     

Convergence of an operator splitting method on a bounded domain for a convection-diffusion-reaction system

We solve a convection-diffusion-sorption (reaction) system on a bounded domain with dominant convection using an operator splitting method. The model arises in contaminant transport in groundwater induced by a dual-well, or in controlled laboratory experiments. The operator splitting transforms the original problem to three subproblems: nonlinear convection, nonlinear diffusion, and a reaction problem, each with its own boundary conditions. The transport equation is solved by a Riemann solver, the diffusion one by a finite volume method, and the reaction equation by an approximation of an integral equation. This approach has proved to be very successful in solving the problem, but the convergence properties where not fully known. We show how the boundary conditions must be taken into account, and prove convergence in L1,loc of the fully discrete splitting procedure to the very weak solution of the original system based on compactness arguments via total variation estimates. Generally, this is the best convergence obtained for this type of approximation. The derivation indicates limitations of the approach, being able to consider only some types of boundary conditions. A sample numerical experiment of a problem with an analytical solution is given, showing the stated efficiency of the method.     

Averaging a transport equation with small diffusion and oscillating velocity

A complete asymptotic expansion is constructed for the transport equation with diffusion term small with respect to the convection. Error estimates are obtained by using matched asymptotic expansion technique and building all the boundary layer terms in time and in space, necessary for obtaining an accurate error estimate. Copyright © 2003 John Wiley & Sons, Ltd.     

Boundary layers for cellular flows at high Péclet numbers

We analyze the behavior of solutions of steady advection‐diffusion problems in bounded domains with prescribed Dirichlet data when the Péclet number Pe ? 1 is large. We show that the solution converges to a constant in each flow cell outside a boundary layer of width O(?^1/2), ? = Pe^?1, around the flow separatrices. We construct an ?‐dependent approximate “water pipe problem” purely inside the boundary layer that provides a good approximation of the solution of the full problem but has ?‐independent computational cost. We also define an asymptotic problem on the graph of streamline separatrices and show that solution of the water pipe problem itself may be approximated by an asymptotic, ?‐independent problem on this graph. Finally, we show that the Dirichlet‐to‐Neumann map of the water pipe problem approximates the Dirichlet‐to‐Neumann map of the separatrix problem with an error independent of the flow outside the boundary layers. © 2004 Wiley Periodicals, Inc.     

Feynman formulas and functional integrals for diffusion with drift in a domain on a manifold

We obtain representations for the solution of the Cauchy-Dirichlet problem for the diffusion equation with drift in a domain on a compact Riemannian manifold as limits of integrals over the Cartesian powers of the domain; the integrands are elementary functions depending on the geometric characteristics of the manifold, the coefficients of the equation, and the initial data. It is natural to call such representations Feynman formulas. Besides, we obtain representations for the solution of the Cauchy-Dirichlet problem for the diffusion equation with drift in a domain on a compact Riemannian manifold as functional integrals with respect to Weizsäcker-Smolyanov surface measures and the restriction of the Wiener measure to the set of trajectories in the domain; such a restriction of the measure corresponds to Brownian motion in a domain with absorbing boundary. In the proof, we use Chernoff’s theorem and asymptotic estimates obtained in the papers of Smolyanov, Weizsäcker, and their coauthors.     

A mean field equation as limit of nonlinear diffusions with fractional Laplacian operators

In the limit of a nonlinear diffusion model involving the fractional Laplacian we get a “mean field” equation arising in superconductivity and superfluidity. For this equation, we obtain uniqueness, universal bounds and regularity results. We also show that solutions with finite second moment and radial solutions admit an asymptotic large time limiting profile which is a special self-similar solution: the “elementary vortex patch”.     

Asymptotic properties of a nonlinear diffusion process arising in articular cartilage

Stress relaxation experiments in articular cartilage studies give rise to a nonlinear diffusion process. The paper demonstrates that the nonlinear model admits a solution which possesses certain salient features that are observed experimentally. A relaxation constant is defined for the nonlinear model in order to examine the asymptotic behavior of solutions as time approaches infinity. An explicit formula relating the relaxation constant to the nonlinear permeability function, and the displacement at the “motor cutoff time” is derived. Finally, upper and lower bounds to the solution and certain of its derivatives are provided.     

Direct and inverse boundary value problems for models of stationary reaction–convection–diffusion

Direct and inverse boundary value problems for models of stationary reaction–convection–diffusion are investigated. The direct problem consists in finding a solution of the corresponding boundary value problem for given data on the boundary of the domain of the independent variable. The peculiarity of the direct problem consists in the inhomogeneity and irregularity of mixed boundary data. Solvability and stability conditions are specified for the direct problem. The inverse boundary value problem consists in finding some traces of the solution of the corresponding boundary value problem for given standard and additional data on a certain part of the boundary of the domain of the independent variable. The peculiarity of the inverse problem consists in its ill-posedness. Regularizing methods and solution algorithms are developed for the inverse problem.     

On heat conductors with a stationary hot spot

We consider a heat conductor having initial constant temperature and zero boundary temperature at every time.The hot spot is the point at which temperature attains its maximum at each given time. For convex conductors, if the hot spot does not move in time, we prove symmetry results for planar triangular and quadrangular conductors.Then, we examine the case of a general conductor and, by an asymptotic formula, we prove that, if there is a stationary critical point, not necessarily the hot spot, then the conductor must satisfy a geometric condition. In particular, we show that there is no stationary critical point inside planar non-convex quadrangular conductors. Mathematics Subject Classification (2000) Primary 35K05, 35K20, 35J05; Secondary 35J25, 35B38, 35B40     

Stability and bifurcation in a two-species reaction–diffusion–advection competition model with time delay

In this paper, we are concerned with the dynamics of a class of two-species reaction–diffusion–advection competition models with time delay subject to the homogeneous Dirichlet boundary condition or no-flux boundary condition in a bounded domain. The existence of steady state solution is investigated by means of the Lyapunov–Schmidt reduction method. The stability and Hopf bifurcation at the spatially nonhomogeneous steady-state are obtained by analyzing the distribution of the associated eigenvalues. Finally, the effect of advection on Hopf bifurcation is explored, which shows that with the increase of convection rate, the Hopf bifurcation phenomenon is more likely to emerge.     

On the asymptotics of a solution to an equation with a small parameter at some of the highest derivatives

We study the asymptotic behavior of a solution of the first boundary value problem for a second-order elliptic equation in a nonconvex domain with smooth boundary in the case where a small parameter is a factor at only some of the highest derivatives and the limit equation is an ordinary differential equation. Although the limit equation has the same order as the initial equation, the problem is singulary perturbed. The asymptotic behavior of its solution is studied by the method of matched asymptotic expansions.     

Iterative solution of a singular convection-diffusion perturbation problem

An iterative domain decomposition method is developed to solve a singular perturbation problem. The problem consists of a convection-diffusion equation with a discontinuous (piecewise-constant) diffusion coefficient, and the problem domain is decomposed into two subdomains, on each of which the coefficient is constant. After showing that the boundary value problem is well posed, we indicate a specific numerical implementation of the iterative technique that combines the finite element method on one subdomain with the method of matched asymptotic expansions on the other subdomain. This procedure extends work by Carlenzoli and Quarteroni, which was originally intended for a boundary layer problem with an outer region and an inner region. Our extension carries over to a problem where the domain consists of the outer and inner boundary layer regions plus a region in which the diffusion coefficient is constant and significant in magnitude. An unexpected benefit of our new implementation is its efficiency, which is due to the fact that at each iteration the problem needs to be solved explicitly only on one subdomain. It is only when the final approximation on the entire domain is desired that the matched asymptotic expansions approximation need be computed on the second subdomain. Two-dimensional convergence results and numerical results illustrating the method for a two-dimensional test problem are given.     

Global asymptotic stability of positive equilibrium of three-species Lotka–Volterra mutualism models with diffusion and delay effects

In the mutualism system with three species if the effects of dispersion and time delays are both taken into consideration, then the densities of the cooperating species are governed by a coupled system of reaction–diffusion equations with time delays. The aim of this paper is to investigate the asymptotic behavior of the time-dependent solution in relation to a positive uniform solution of the corresponding steady-state problem in a bounded domain with Neumann boundary condition, including the existence and uniqueness of a positive steady-state solution. A simple and easily verifiable condition is given to ensure the global asymptotic stability of the positive steady-state solution. This result leads to the permanence of the mutualism system, the instability of the trivial and all forms of semitrivial solutions, and the nonexistence of nonuniform steady-state solution. The condition for the global asymptotic stability is independent of diffusion and time-delays as well as the net birth rate of species, and the conclusions for the reaction–diffusion system are directly applicable to the corresponding ordinary differential system and 2-species cooperating reaction–diffusion systems. Our approach to the problem is based on inequality skill and the method of upper and lower solutions for a more general reaction–diffusion system. Finally, the numerical simulation is given to illustrate our results.