Mathematical Models and Methods in Applied Sciences最新文献

We focus on an epidemiological model (the archetypical SIR system) defined on graphs and study the asymptotic behavior of the solutions as the number of vertices in the graph diverges. By relying on the theory of graphons we provide a characterization of the limit and establish convergence results. We also provide approximation results for both deterministic and random discretizations.

This paper introduces a novel approach to approximate a broad range of reaction–convection–diffusion equations using conforming finite element methods while providing a discrete solution respecting the physical bounds given by the underlying differential equation. The main result of this work demonstrates that the numerical solution achieves an accuracy of $O(h^k)$ in the energy norm, where $k$ represents the underlying polynomial degree. To validate the approach, a series of numerical experiments had been conducted for various problem instances. Comparisons with the linear continuous interior penalty stabilised method, and the algebraic flux-correction scheme (for the piecewise linear finite element case) have been carried out, where we can observe the favorable performance of the current approach.

We study the global well-posedness and asymptotic behavior for a semilinear damped wave equation with Neumann boundary conditions, modeling a one-dimensional linearly elastic body interacting with a rigid substrate through an adhesive material. The key feature of of the problem is that the interplay between the nonlinear force and the boundary conditions allows for a continuous set of equilibrium points. We prove an exponential rate of convergence for the solution towards a (uniquely determined) equilibrium point.

In this paper, 3D–2D-dimensional reduction for hyperelastic thin films modeled through energies with point-dependent growth, assuming that the sample is clamped on the lateral boundary, is performed in the framework of $\mathrm{\Gamma }$-convergence. Integral representation results, with a more regular Lagrangian related to the original energy density, are provided for the lower dimensional limiting energy, in different contexts.

In this paper, we present a mathematical study for the development of Multiple Sclerosis in which a spatio-temporal kinetic theory model describes, at the mesoscopic level, the dynamics of a high number of interacting agents. We consider both interactions among different populations of human cells and the motion of immune cells, stimulated by cytokines. Moreover, we reproduce the consumption of myelin sheath due to anomalously activated lymphocytes and its restoration by oligodendrocytes. Successively, we fix a small time parameter and assume that the considered processes occur at different scales. This allows us to perform a formal limit, obtaining macroscopic reaction–diffusion equations for the number densities with a chemotaxis term. A natural step is then to study the system, inquiring about the formation of spatial patterns through a Turing instability analysis of the problem and basing the discussion on the microscopic parameters of the model. In particular, we get spatial patterns oscillating in time that may reproduce brain lesions characteristic of different phases of the pathology.

In this work, using the theory of first-order macroscopic crowd models, we introduce a compartmental advection–diffusion model, describing the spatio-temporal dynamics of a population in different human behaviors (alert, panic and control) during a catastrophic event. For this model, we prove the local existence, uniqueness and regularity of a solution, as well as the positivity and $L^1$-boundedness of this solution. Then, in order to study the spatio-temporal propagation of these behavioral reactions within a population during a catastrophic event, we present several numerical simulations for different evacuation scenarios.

In this paper, we consider the heat-conducting compressible self-gravitating fluids in time-dependent domains, which typically describe the motion of viscous gaseous stars. The flow is governed by the 3D Navier–Stokes–Fourier–Poisson equations where the velocity is supposed to fulfill the full-slip boundary condition and the temperature on the boundary is given by a non-homogeneous Dirichlet condition. We establish the global-in-time weak solution to the system. Our approach is based on the penalization of the boundary behavior, viscosity, and the pressure in the weak formulation. Moreover, to accommodate the non-homogeneous boundary heat flux, the concept of ballistic energy is utilized in this work.