No abstract uploaded!

No abstract uploaded!

In this paper, we propose a new stabilized linear finite element method for solving reaction-convection-diffusion equations with arbitrary magnitudes of reaction and diffusion. The key feature of the new method is that the test function in the stabilization term is taken in the adjoint-operator-like form $-\varepsilon \Delta v- ({\mbf a}\cdot\nabla v)/\gamma+\sigma v$, where the parameter $\gamma$ is appropriately designed to adjust the convection strength to achieve high accuracy and stability. We derive the stability estimates for the finite element solutions and establish the explicit dependence of $L^2$ and $H^1$ error bounds on the diffusivity, modulus of the convection field, reaction coefficient and the mesh size. The analysis shows that the proposed method is suitable for a wide range of mesh P\'{e}clet numbers and mesh Damk\"{o}hler numbers. More specifically, if the diffusivity $\varepsilon$ is sufficiently small with $\varepsilon < \Vert{\mbf a}\Vert h$ and the reaction coefficient $\sigma$ is large enough such that $\Vert{\mbf a}\Vert < \sigma h$, then the method exhibits optimal convergence rates in both $L^2$ and $H^1$ norms. However, for a small reaction coefficient satisfying $\Vert{\mbf a}\Vert \ge \sigma h$, the method behaves like the well-known streamline upwind/Petrov-Galerkin formulation of Brooks and Hughes. Several numerical examples exhibiting boundary or interior layers are given to demonstrate the high performance of the proposed method. Moreover, we apply the developed method to time-dependent reaction-convection-diffusion problems and simulation results show the efficiency of the approach.

We deal with a stationary problem of a reaction–diffusion system with a conservation law under the Neu-mann boundary condition. It is shown that the stationary problem turns to be the Euler–Lagrange equation of an energy functional with a mass constraint. When the domain is the finite interval (0, 1), we investigate the asymptotic profile of a strictly monotone minimizer of the energy as d, the ratio of the diffusion coef-ficient of the system, tends to zero. In view of a logarithmic function in the leading term of the potential, we get to a scaling parameter κsatisfying the relation ε:=√d=√logκ/κ2. The main result shows that a sequence of minimizers converges to a Dirac mass multiplied by the total mass and that by a scaling with κthe asymptotic profile exhibits a parabola in the nonvanishing region. We also prove the existence of an unstable monotone solution when the mass is small.

In this article we show how the classical theory of Reeb and others extends to the case of codimension one singular foliations on closed oriented manifolds, provided that at the singular set sing(F), the foliation is locally defined by Bott-Morse functions which are transversely centers. We prove, in this setting, the equivalent of the local and the complete stability theorems of Reeb. We show that if F has a compact leaf with finite fundamental group, or if a component of sing(F) has codimension ≥ 3 and finite fundamental group, then all leaves of F are compact and diffeomorphic, sing(F) consists of two connected components, and there is a Bott-Morse function f : M → [0, 1] such that f : M \ sing(F) → (0, 1) is a fiber bundle defining F and sing(F) = f.1({0, 1}). This yields a topological description of the type of leaves that appear in these foliations, and also the type of manifolds admitting such foliations. These results unify and generalize well known results for cohomogeneity one isometric actions, and a theorem of Reeb for foliations with Morse singularities of center type.