Discontinuous Galerkin method for nonlinear convection-diffusion problems with mixed Dirichlet-Neumann boundary conditions

The paper is devoted to the analysis of the discontinuous Galerkin finite element method (DGFEM) applied to the space semidiscretization of a nonlinear nonstationary convection-diffusion problem with mixed Dirichlet-Neumann boundary conditions. General nonconforming meshes are used and the NIPG, IIPG and SIPG versions of the discretization of diffusion terms are considered. The main attention is paid to the impact of the Neumann boundary condition prescribed on a part of the boundary on the truncation error in the approximation of the nonlinear convective terms. The estimate of this error allows to analyse the error estimate of the method. The results obtained represent the completion and extension of the analysis from V. Dolej?í, M. Feistauer, Numer. Funct. Anal. Optim. 26 (2005), 349–383, where the truncation error in the approximation of the nonlinear convection terms was proved only in the case when the Dirichlet boundary condition on the whole boundary of the computational domain was considered.     

A weak Galerkin finite element method for the second order elliptic problems with mixed boundary conditions

In this paper, a weak Galerkin finite element method is proposed and analyzed for the second-order elliptic equation with mixed boundary conditions. Optimal order error estimates are established in both discrete $H^1$ norm and the standard $L^2$ norm for the corresponding WG approximations. The numerical experiments are presented to verify the efficiency of the method.     

Boundedness of solutions to a quasilinear parabolic–elliptic Keller–Segel system with logistic source

We study a quasilinear parabolic–elliptic Keller–Segel system involving a source term of logistic type u_t = ? ? (?(u) ? u) ? χ ? ? (u ? v) + g(u), ? Δv = ? v + u in Ω × (0,T), subject to nonnegative initial data and the homogeneous Neumann boundary condition in a bounded domain with smooth boundary, n ≥ 1, χ > 0, ? ≥ c₁s^p for s ≥ s₀ > 1, and g(s) ≤ as ? μs² for s > 0 with a,g(0) ≥ 0, μ > 0. There are three nonlinear mechanisms included in the chemotaxis model: the nonlinear diffusion, aggregation and logistic absorption. The interaction among the triple nonlinearities shows that together with the nonlinear diffusion, the logistic absorption will dominate the aggregation such that the unique classical solution of the system has to be global in time and bounded, regardless of the initial data, whenever , or, equivalently, , which enlarge the parameter range , or , required by globally bounded solutions of the quasilinear K‐S system without the logistic source. Copyright © 2013 John Wiley & Sons, Ltd.     

Galerkin methods for parabolic equations with nonlinear boundary conditions

A variety of Galerkin methods are studied for the parabolic equationu _t =(a(x) u),x ⁿ ,t (O,T], subject to the nonlinear boundary conditionu _v =g(x,t,u),x,t (O,T] and the usual initial condition. Optimal order error estimates are derived both inL ² () andH ¹ () norms for all methods treated, including several that produce linear computational procedures.The authors were partially supported by The National Science Foundation during the preparation of this paper.     

Boundedness in a parabolic‐elliptic chemotaxis system with nonlinear diffusion and sensitivity and logistic source

In this paper, we study the zero‐flux chemotaxis‐system where Ω is a bounded and smooth domain of , n≥1, and where , k,μ>0 and α≤1. For any v≥0, the chemotactic sensitivity function is assumed to behave as the prototype χ(v)=χ₀/(1+av)², with a≥0 and χ₀>0. We prove that for any nonnegative and sufficiently regular initial data u(x,0), the corresponding initial‐boundary value problem admits a unique global bounded classical solution if α<1; indeed, for α=1, the same conclusion is obtained provided μ is large enough. Finally, we illustrate the range of dynamics present within the chemotaxis system in 1, 2, and 3 dimensions by means of numerical simulations.     

Boundedness and stabilization in a two‐species chemotaxis‐competition system of parabolic‐parabolic–elliptic type

This paper is concerned with the two‐species chemotaxis‐competition system where Ω is a bounded domain in with smooth boundary ∂Ω, n≥2; χ_i and μ_i are constants satisfying some conditions. The above system was studied in the cases that a₁,a₂∈(0,1) and a₁>1>a₂, and it was proved that global existence and asymptotic stability hold when are small. However, the conditions in the above 2 cases strongly depend on a₁,a₂, and have not been obtained in the case that a₁,a₂≥1. Moreover, convergence rates in the cases that a₁,a₂∈(0,1) and a₁>1>a₂ have not been studied. The purpose of this work is to construct conditions which derive global existence of classical bounded solutions for all a₁,a₂>0 which covers the case that a₁,a₂≥1, and lead to convergence rates for solutions of the above system in the cases that a₁,a₂∈(0,1) and a₁≥1>a₂.     

Solutions for systems of nonlinear elliptic equations with nonlinear boundary conditions

We look for solutions of systems of nonlinear elliptic equations with nonlinear boundary conditions and values in some compact convex set M. If the nonlinear terms satisfy a sign condition on the boundary of M and the inhomogeneous terms assume their values in this set existence of solutions is proved. The proof is based on the homotopy invariance of the Leray-Schauder degree and Weinberger's strong maximum principle.     

Eigenproblem for p-Laplacian and nonlinear elliptic equation with nonlinear boundary conditions

In this article, we study the solvability of nonlinear problem for p-Laplacian with nonlinear boundary conditions. We give some characterization of the first eigenvalue of an intermediary eigenvalue problem as simplicity, isolation and its strict monotonicity. Afterward, we character also the second eigenvalue and its strictly partial monotony. On the other hand, in some sense, we establish the non-resonance below the first and furthermore between the first and second eigenvalues of nonlinear Steklov–Robin.     

Discrete maximum principles for finite element solutions of nonlinear elliptic problems with mixed boundary conditions

Summary. One of the most important problems in numerical simulations is the preservation of qualitative properties of solutions of the mathematical models by computed approximations. For problems of elliptic type, one of the basic properties is the (continuous) maximum principle. In our work, we present several variants of the maximum principles and their discrete counterparts for (scalar) second-order nonlinear elliptic problems with mixed boundary conditions. The problems considered are numerically solved by the continuous piecewise linear finite element approximations built on simplicial meshes. Sufficient conditions providing the validity of the corresponding discrete maximum principles are presented. Geometrically, they mean that the employed meshes have to be of acute or nonobtuse type, depending of the type of the problem. Finally some examples of real-life problems, where the preservation of maximum principles plays an important role, are presented.The first author was supported by the Hungarian Research Fund OTKA under grant no. F034840The second author was supported by the Agora Center under Grant InBCT of TEKES, Finland, and by the Academy Research Fellowship no. 208628 from the Academy of FinlandMathematics Subject Classification (2000): 35B50, 35J65, 65N30, 65N50     

Multiplicity of solutions to a nonlinear elliptic problem with nonlinear boundary conditions

We study the problem $$\left\{\begin{array}{ll}\Delta_p u = |u|^{q-2}u, & \quad x \in \Omega ,\\ |\nabla u|^{p-2} \frac{\partial u}{\partial \nu}= \lambda |u|^{p-2}u, &\quad x \in \partial \Omega, \end{array}\right.$$ where \({\Omega \subset \mathbb{R}^N}\) is a bounded smooth domain, \({\nu}\) is the outward unit normal at \({\partial \Omega}\) and \({\lambda > 0}\) is regarded as a bifurcation parameter. When p = 2 and in the superlinear regime q > 2, we show existence of n nontrivial solutions for all \({\lambda > \lambda_n}\) , \({\lambda_n}\) being the n-th Steklov eigenvalue. It is proved in addition that bifurcation from the trivial solution takes place at all \({\lambda_n}\) ’s. Similar results are obtained in the sublinear case 1 < q < 2. In this case, bifurcation from infinity takes place in those \({\lambda_n}\) with odd multiplicity. Partial extensions of these features are shown in the nonlinear diffusion case \({p \neq 2}\) and related problems under spatially heterogeneous reactions are also addressed.     

Semilinear elliptic problems with mixed Dirichlet-Neumann boundary conditions

This work deals with the analysis of problems

     

On a parabolic‐elliptic chemotaxis system with periodic asymptotic behavior

We study a parabolic‐elliptic chemotactic PDEs system, which describes the evolution of a biological population “u” and a chemical substance “v” in a bounded domain . We consider a growth term of logistic type in the equation of “u” in the form μu(1 ? u + f(t,x)). The function “f,” describing the resources of the systems, presents a periodic asymptotic behavior in the sense where f ^? is independent of x and periodic in time. We study the global existence of solutions and its asymptotic behavior. Under suitable assumptions on the initial data and f ^?, if the constant chemotactic sensitivity χ satisfies we obtain that the solution of the system converges to a homogeneous in space and periodic in time function.     

Steklov-Neumann eigenproblems and nonlinear elliptic equations with nonlinear boundary conditions

We study the solvability of nonlinear second order elliptic partial differential equations with nonlinear boundary conditions. We introduce the notion of “eigenvalue-lines” in the plane; these eigenvalue-lines join each Steklov eigenvalue to the first eigenvalue of the Neumann problem with homogeneous boundary condition. We prove existence results when the nonlinearities involved asymptotically stay, in some sense, below the first eigenvalue-lines or in a quadrilateral region (depicted in Fig. 1) enclosed by two consecutive eigenvalue-lines. As a special case we derive the so-called nonresonance results below the first Steklov eigenvalue as well as between two consecutive Steklov eigenvalues. The case in which the eigenvalue-lines join each Neumann eigenvalue to the first Steklov eigenvalue is also considered. Our method of proof is variational and relies mainly on minimax methods in critical point theory.     

Nontrivial solutions for a nonlinear discrete elliptic equation with periodic boundary conditions

In this paper, existence of nontrivial solutions for a nonlinear discrete elliptic equation with periodic boundary conditions is considered by using the monotone operator principle, Mountain pass lemma, Linking theorem and some results in Morse theory. The obtained results are also valid and new for the corresponding difference periodic boundary value problem.     

Some results for non-linear elliptic problems with mixed boundary conditions

We prove existence and uniqueness results for non-linear elliptic equations with lower order terms, L¹ data, and mixed boundary conditions that include as particular cases the Dirichlet and the Neumann problems. Mathematics Subject Classification (2000) 35J25, 35D05, 35J70, 35J60     

Regularity results for quasilinear elliptic systems with nonlinear boundary conditions

It is proved that a solution of the boundary-value problem for a second-order quasilinear system with controlled order of nonlinearity is partially smooth all the way to the boundary of a domain. The boundary condition is imposed by means of a second-order nonlinear operator which can be regarded as a generalization of the “directional derivative” to the case of quasilinear systems. Bibliography: 6 titles. Translated fromProblemy Matematicheskogo Analiza, No. 14, 1995, pp. 23–50.     

A quasilinearization method for elliptic problems with a nonlinear boundary condition

We study a nonlinear elliptic second order problem with a nonlinear boundary condition. Assuming the existence of an ordered couple of a supersolution and a subsolution, we develop a quasilinearization method in order to construct an iterative scheme that converges to a solution. Furthermore, under an extra assumption we prove that the convergence is quadratic.     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

Two-sided a posteriori error estimates for linear elliptic problems with mixed boundary conditions

The paper is devoted to verification of accuracy of approximate solutions obtained in computer simulations. This problem is strongly related to a posteriori error estimates, giving computable bounds for computational errors and detecting zones in the solution domain where such errors are too large and certain mesh refinements should be performed. A mathematical model consisting of a linear elliptic (reaction-diffusion) equation with a mixed Dirichlet/Neumann/Robin boundary condition is considered in this work. On the base of this model, we present simple technologies for straightforward constructing computable upper and lower bounds for the error, which is understood as the difference between the exact solution of the model and its approximation measured in the corresponding energy norm. The estimates obtained are completely independent of the numerical technique used to obtain approximate solutions and are “flexible” in the sense that they can be, in principle, made as close to the true error as the resources of the used computer allow. This work was supported by the Academy Research Fellowship No. 208628 from the Academy of Finland.