This study deals with higher-ordered asymptotic equations for the water-waves problem. We considered the higher-order/extended Boussinesq equations over a flat bottom topography in the well-known long wave regime. Providing an existence and uniqueness of solution on a relevant time scale of order $1/\sqrt{\eps}$ and showing that the solution's behavior is close to the solution of the water waves equations with a better precision corresponding to initial data, the asymptotic model is well-posed in the sense of Hadamard. Then we compared several water waves solitary solutions with respect to the numerical solution of our model. At last, we solve explicitly this model and validate the results numerically.
We develop a mixed finite element method for the coupled problem arising in the interaction between a free fluid governed by the Stokes equations and flow in deformable porous medium modeled by the Biot system of poroelasticity. Mass conservation, balance of stress, and the Beavers-Joseph-Saffman condition are imposed on the interface. We consider a fully mixed Biot formulation based on a weakly symmetric stress-displacement-rotation elasticity system and Darcy velocity-pressure flow formulation. A velocity-pressure formulation is used for the Stokes equations. The interface conditions are incorporated through the introduction of the traces of the structure velocity and the Darcy pressure as Lagrange multipliers. Existence and uniqueness of a solution are established for the continuous weak formulation. Stability and error estimates are derived for the semi-discrete continuous-in-time mixed finite element approximation. Numerical experiments are presented to verify the theoretical results and illustrate the robustness of the method with respect to the physical parameters.
A new interpolation operator is defined, which preserves only P2 polynomials locally.
This work is devoted to the study of the approximation, using two nonlinear numerical methods, of a linear elliptic problem with measure data and heterogeneous anisotropic diffusion matrix. Both methods show convergence properties to a continuous solution of the problem in a weak sense, through the change of variable u = ψ(v), where ψ is a well chosen diffeomorphism between (−1, 1) and ℝ, and v is valued in (−1, 1). We first study a nonlinear finite element approximation on any simplicial grid. We prove the existence of a discrete solution, and, under standard regularity conditions, we prove its convergence to a weak solution of the problem by applying Hölder and Sobolev inequalities. Some numerical results, in 2D and 3D cases where the solution does not belong to H1 (Ω), show that this method can provide accurate results. We then construct a numerical scheme which presents a convergence property to the entropy weak solution of the problem in the case where the right-hand side belongs to L1. This is achieved owing to a nonlinear control volume finite element (CVFE) method, keeping the same nonlinear reformulation, and adding an upstream weighting evaluation and a nonlinear p-Laplace vanishing stabilisation term.
In this paper, we conduct a thorough mathematical analysis of the unified formulation advocated by Lauser et al. [Adv. Water Res. 34 (2011) 957-966] for compositional multiphase flows in porous media. The interest of this formulation lies in its potential to automatically handle the appearance and disappearance of phases. However, its practical implementation turned out to be not always robust for realistic fugacity laws associated with cubic equations of state, as shown by Ben Gharbia and Flauraud [Oil Gas Sci. Technol. 74 (2019) 43]. By focusing on the subproblem of phase equilibrium, we derive sufficient conditions for the existence of the corresponding system of equations. We trace back the difficulty of cubic laws to a deficiency of the Gibbs functions that comes into play due to the ``unifying'' feature of the new formulation. We propose a partial remedy for this problem by extending the domain of definition of these functions in a natural way. Besides, we highlight the crucial but seemingly unknown fact that the unified formulation encapsulates all the properties known to physicists on phase equilibrium, such as the tangent plane criterion and the minimization of the Gibbs energy of the mixture.
Electrical Impedance Tomography (EIT) is a well-known imaging technique for detecting the electrical properties of an object in order to detect anomalies, such as conductive or resistive targets. More specifically, EIT has many applications in medical imaging for the detection and location of bodily tumors since it is an affordable and non-invasive method, which aims to recover the internal conductivity of a body using voltage measurements resulting from applying low frequency current at electrodes placed at its surface. Mathematically, the reconstruction of the internal conductivity is a severely ill-posed inverse problem and yields a poor quality image reconstruction. To remedy this difficulty, at least in part, we regularize and solve the nonlinear minimization problem by the aid of a Krylov subspace-type method for the linear sub problem during each iteration. In EIT, a tumor or general anomaly can be modeled as a piecewise constant perturbation of a smooth background, hence, we solve the regularized problem on a subspace of relatively small dimension by the Flexible Golub-Kahan process that provides solutions that have sparse representation. For comparison, we use a well-known modified Gauss–Newton algorithm as a benchmark. Using simulations, we demonstrate the effectiveness of the proposed method. The obtained reconstructions indicate that the Krylov subspace method is better adapted to solve the ill-posed EIT problem and results in higher resolution images and faster convergence compared to reconstructions using the modified Gauss–Newton algorithm.
We present a structured coagulation-fragmentation model which describes the population dynamics of oceanic phytoplankton. This model is formulated on the space of Radon measures equipped with the bounded Lipschitz norm and unifies the study of the discrete and continuous coagulation-fragmentation models. We prove that the model is well-posed and show it can reduce down to the classic discrete and continuous coagulation-fragmentation models. To understand the interplay between the physical processes of coagulation and fragmentation and the biological processes of growth, reproduction, and death, we establish a regularity result for the solutions and use it to show that stationary solutions are absolutely continuous under some conditions on model parameters. We develop a semi-discrete approximation scheme which conserves mass and prove its convergence to the unique weak solution. We then use the scheme to perform numerical simulations for the model.
We analyze the approximation by mixed finite element methods of solutions of equations of the form −div (a∇u) = g, where the coefficient a = a(x) can degenerate going to zero or infinity. First, we extend the classic error analysis to this case provided that the coefficient a belongs to the Muckenhoupt class A2. The analysis developed applies to general mixed finite element spaces satisfying the standard commutative diagram property, whenever some stability and interpolation error estimates are valid in weighted norms. Next, we consider in detail the case of Raviart–Thomas spaces of arbitrary order, obtaining optimal order error estimates for simplicial elements in any dimension and for convex quadrilateral elements in the two dimensional case, in both cases under a regularity assumption on the family of meshes. For the lowest order case we show that the regularity assumption can be removed and prove anisotropic error estimates which are of interest in problems with boundary layers. Finally we apply the results to a problem arising in the solution of the fractional Laplace equation.
We deal with the numerical investigation of the local limit of nonlocal conservation laws. Previous numerical experiments suggest convergence in the local limit. However, recent analytic results state that (i) in general convergence does not hold because one can exhibit counterexamples; (ii) convergence can be recovered provided viscosity is added to both the local and the nonlocal equations. Motivated by these analytic results, we investigate the role of numerical viscosity in the numerical study of the local limit of nonlocal conservation laws. In particular, we show that the numerical viscosity of Lax-Friedrichs type schemes jeopardizes the reliability of the numerical scheme and erroneously detects convergence in cases where convergence is ruled out by analytic results. We also test Godunov type schemes, less affected by numerical viscosity, and show that in some cases they provide more reliable results.
This paper is devoted to the mathematical modeling and analysis of a hyperbolic Maxwell quasi-variational inequality (QVI) for the Bean-Kim superconductivity model with temperature and magnetic field dependence in the critical current. Our analysis relies on local (resp. global) boundedness and local (resp. global) Lipschitz continuity assumptions on the critical current with respect to the temperature (resp. magnetic field). Emerging from the Euler time discretization, we analyze the corresponding H(curl)-elliptic QVI and prove its existence using a fixed-point argument in combination with techniques from variational inequalities and Maxwell’s equations. Based on the existence result for the H(curl)-elliptic QVI, we examine the stability and convergence of the Euler scheme, which serve as our fundament for the global well-posedness of the governing hyperbolic Maxwell QVI.
We consider a perforated half-cylindrical thin shell and investigate the limit behavior when the period and the thickness simultaneously go to zero. By using the decomposition of shell displacements presented in Griso [JMPA 89 (2008) 199–223] we obtain a priori estimates. With the unfolding and rescaling operator we transform the problem to a reference configuration. In the end this yields a homogenized limit problem for the shell.
A recent paper introduced the network element method (NEM) where the usual mesh was replaced by a discretization network. Using the associated network geometric coefficients and following the virtual element framework, a consistent and stable numerical scheme was proposed. The aim of the present paper is to derive a convergence theory for the NEM under mild assumptions on the exact problem. We also derive basic error estimates, which are sub-optimal in the sense that we have to assume more regularity than usual.
We prove well-posedness and regularity of solutions to a fractional diffusion porous media equation with a variable fractional order that may depend on the unknown solution. We present a linearly implicit time-stepping method to linearize and discretize the equation in time, and present rigorous analysis for the convergence of numerical solutions based on proved regularity results.
In this work we introduce and study a nonlocal version of the PageRank. In our approach, the random walker explores the graph using longer excursions than just moving between neighboring nodes. As a result, the corresponding ranking of the nodes, which takes into account a long-range interaction between them, does not exhibit concentration phenomena typical of spectral rankings which take into account just local interactions. We show that the predictive value of the rankings obtained using our proposals is considerably improved on different real world problems.
In this paper the nonlinear multi-species Boltzmann equation with random uncertainty coming from the initial data and collision kernel is studied. Well-posedness and long-time behavior – exponential decay to the global equilibrium – of the analytical solution, and spectral gap estimate for the corresponding linearized gPC-based stochastic Galerkin system are obtained, by using and extending the analytical tools provided in [M. Briant and E.S. Daus, Arch. Ration. Mech. Anal. 3 (2016) 1367–1443] for the deterministic problem in the perturbative regime, and in [E.S. Daus, S. Jin and L. Liu, Kinet. Relat. Models 12 (2019) 909–922] for the single-species problem with uncertainty. The well-posedness result of the sensitivity system presented here has not been obtained so far neither in the single species case nor in the multi-species case.
In this paper, thea priorierror estimates of an embedded discontinuous Galerkin method for the Oseen equations are presented. It is proved that the velocity error in theL2(Ω) norm, has an optimal error bound with convergence orderk + 1, where the constants are dependent on the Reynolds number (orν−1), in the diffusion-dominated regime, and in the convection-dominated regime, it has a Reynolds-robust error bound with quasi-optimal convergence orderk + 1/2. Here,kis the polynomial order of the velocity space. In addition, we also prove an optimal error estimate for the pressure. Finally, we carry out some numerical experiments to corroborate our analytical results.
This paper is concerned with the Richards equation in a heterogeneous domain, each subdomain of which is homogeneous and represents a rocktype. Our first contribution is to rigorously prove convergence toward a weak solution of cell-centered finite-volume schemes with upstream mobility and without Kirchhoff's transform. Our second contribution is to numerically demonstrate the relevance of locally refining the grid at the interface between subregions, where discontinuities occur, in order to preserve an acceptable accuracy for the results computed with the schemes under consideration.
We address fundamental aspects in the approximation theory of vector-valued finite element methods, using finite element exterior calculus as a unifying framework. We generalize the Clément interpolant and the Scott-Zhang interpolant to finite element differential forms, and we derive a broken Bramble-Hilbert lemma. Our interpolants require only minimal smoothness assumptions and respect partial boundary conditions. This permits us to state local error estimates in terms of the mesh size. Our theoretical results apply to curl-conforming and divergence-conforming finite element methods over simplicial triangulations.
The context of this paper is the simulation of parameter-dependent partial differential equations (PDEs). When the aim is to solve such PDEs for a large number of parameter values, Reduced Basis Methods (RBM) are often used to reduce computational costs of a classical high fidelity code based on Finite Element Method (FEM), Finite Volume (FVM) or Spectral methods. The efficient implementation of most of these RBM requires to modify this high fidelity code, which cannot be done, for example in an industrial context if the high fidelity code is only accessible as a "black-box" solver. The Non-Intrusive Reduced Basis (NIRB) method has been introduced in the context of finite elements as a good alternative to reduce the implementation costs of these parameter-dependent problems. The method is efficient in other contexts than the FEM one, like with finite volume schemes, which are more often used in an industrial environment. In this case, some adaptations need to be done as the degrees of freedom in FV methods have different meanings. At this time, error estimates have only been studied with FEM solvers. In this paper, we present a generalisation of the NIRB method to Finite Volume schemes and we show that estimates established for FEM solvers also hold in the FVM setting. We first prove our results for the hybrid-Mimetic Finite Difference method (hMFD), which is part the Hybrid Mixed Mimetic methods (HMM) family. Then, we explain how these results apply more generally to other FV schemes. Some of them are specified, such as the Two Point Flux Approximation (TPFA).
Finite difference kind of schemes are popular in approximating wave propagation problems in finite dimensional spaces. While Yee’s original paper on the finite difference method is already from the sixties, mathematically there still remains questions which are not yet satisfactorily covered. In this paper, we address two issues of this kind. Firstly, in the literature Yee’s scheme is constructed separately for each particular type of wave problem. Here, we explicitly generalize the Yee scheme to a class of wave problems that covers at large physics field theories. For this we introduce Yee’s scheme for all problems of a class characterised on a Minkowski manifold by (i) a pair of first order partial differential equations and by (ii) a constitutive relation that couple the differential equations with a Hodge relation. In addition, we introduce a strategy to systematically exploit higher order Whitney elements in Yee-like approaches. This makes higher order interpolation possible both in time and space. For this, we show that Yee-like schemes preserve the local character of the Hodge relation, which is to say, the constitutive laws become imposed on a finite set of points instead of on all ordinary points of space. As a result, the usage of higher order Whitney forms does not compel to change the actual solution process at all. This is demonstrated with a simple example.
