Uniqueness Of Discrete Solution Research Articles

A Convergent Entropy-Dissipating BDF2 Finite-Volume Scheme for a Population Cross-Diffusion System

Abstract A second-order backward differentiation formula (BDF2) finite-volume discretization for a nonlinear cross-diffusion system arising in population dynamics is studied. The numerical scheme preserves the Rao entropy structure and conserves the mass. The existence and uniqueness of discrete solutions and their large-time behavior as well as the convergence of the scheme are proved. The proofs are based on the G-stability of the BDF2 scheme, which provides an inequality for the quadratic Rao entropy and hence suitable a priori estimates. The novelty is the extension of this inequality to the system case. Some numerical experiments in one and two space dimensions underline the theoretical results.

Analysis and simulation of numerical schemes for nonlinear hyperbolic predator–prey models with spatial diffusion

Departing from a two-dimensional parabolic system that describes the spatial dynamics in a predator–prey system with Michaelis–Menten-type functional response, we investigate a general form of that model using a finite-difference approach. The model under investigation is a hyperbolic nonlinear system consisting of two coupled partial differential equations with generalized reaction terms. We impose initial conditions on a closed and bounded rectangle, and a fully discrete finite-difference methodology is proposed. Among the most important results of this work, we establish analytically the existence and uniqueness of discrete solutions, along with the second-order consistency of our scheme. Moreover, a discrete form of the energy method is employed to prove the stability and the quadratic convergence of the technique. Some numerical simulations obtained through our method show the appearance of Turing patterns in the parabolic case, in agreement with some reports found in the literature. Moreover, our simulations also show that Turing patterns are present in the hyperbolic scenario.

Simple efficient simulation of the complex dynamics of some nonlinear hyperbolic predator–prey models with spatial diffusion

In this manuscript, we will consider a parabolic partial differential equation that models the spatial-temporal dynamics in predator–prey systems with Michaelis–Menten-type functional response. We study numerically an extension of those systems by employing finite differences. The extension includes the presence of inertial times with extended reaction components in a diffusive system consisting of two populations. The spatial domain is a closed and bounded rectangle, and initial-boundary conditions are imposed on the model. An easy-to-implement fully discrete finite-difference methodology is designed here. We establish analytically the existence and uniqueness of discrete solutions as a consequence of the explicit character of the scheme, along with the second-order consistency. To prove the stability and the second-order convergence of the technique, the energy method will be applied conveniently in a discrete manner. From the point of view of the applications, we provide various numerical simulations that show the appearance of complex patterns in the parabolic case. The results will resemble qualitatively those of some papers already published in the specialized literature. Our simulations also show that different patterns are present in the hyperbolic scenario. We also perform a numerical analysis of the convergence of the method.

Augmented Lagrangian finite element methods for contact problems

We propose two different Lagrange multiplier methods for contact problems derived from the augmented Lagrangian variational formulation. Both the obstacle problem, where a constraint on the solution is imposed in the bulk domain and the Signorini problem, where a lateral contact condition is imposed are considered. We consider both continuous and discontinuous approximation spaces for the Lagrange multiplier. In the latter case the method is unstable and a penalty on the jump of the multiplier must be applied for stability. We prove the existence and uniqueness of discrete solutions, best approximation estimates and convergence estimates that are optimal compared to the regularity of the solution.

Galerkin least squares finite element method for the obstacle problem

We construct a consistent multiplier free method for the finite element solution of the obstacle problem. The method is based on an augmented Lagrangian formulation in which we eliminate the multiplier by use of its definition in a discrete setting. We prove existence and uniqueness of discrete solutions and optimal order a priori error estimates for smooth exact solutions. Using a saturation assumption we also prove an a posteriori error estimate. Numerical examples show the performance of the method and of an adaptive algorithm for the control of the discretization error.

Qualitative analysis of 3D elastostatic contact problems with orthotropic Coulomb friction and solution-dependent coefficients of friction

This paper analyzes a discrete form of 3D contact problems with local orthotropic Coulomb friction and coefficients of friction which may depend on the solution itself. The analysis is based on the fixed-point reformulation of the original problem. Conditions guaranteeing the existence and uniqueness of discrete solutions are established. Finally, numerical results of a model example are presented.

An Equal-Order DG Method for the Incompressible Navier-Stokes Equations

We introduce and analyze a discontinuous Galerkin method for the incompressible Navier-Stokes equations that is based on finite element spaces of the same polynomial order for the approximation of the velocity and the pressure. Stability of this equal-order approach is ensured by a pressure stabilization term. A simple element-by-element post-processing procedure is used to provide globally divergence-free velocity approximations. For small data, we prove the existence and uniqueness of discrete solutions and carry out an error analysis of the method. A series of numerical results are presented that validate our theoretical findings.

Finite-difference scheme for two-scale homogenized equations of one-dimensional motion of a thermoviscoelastic Voigt-type body

The Bakhvalov-Eglit two-scale homogenized equations are used to describe the motion of layered periodic compressible media with rapidly oscillating data. A new finite-difference scheme for a system of such equations is proposed and analyzed in the case of a thermoviscoelastic Voigt-type body. A priori estimates of solutions are derived for nonsmooth data. The existence and uniqueness of discrete solutions are established. A theorem is proved on the convergence of a subsequence of discrete solutions to a weak solution of the problem under study. Simultaneously, a new theorem on the existence of global weak solutions is deduced.