An official website of the United States government
We consider the mathematical analysis and numerical approximation of a system of nonlinear partial differential equations that arises in models that have relevance to steady isochoric flows of colloidal suspensions. The symmetric velocity gradient is assumed to be a monotone nonlinear function of the deviatoric part of the Cauchy stress tensor. We prove the existence of a weak solution to the problem, and under the additional assumption that the nonlinearity involved in the constitutive relation is Lipschitz continuous we also prove uniqueness of the weak solution. We then construct mixed finite element approximations of the system using both conforming and nonconforming finite element spaces. For both of these we prove the convergence of the method to the unique weak solution of the problem, and in the case of the conforming method we provide a bound on the error between the analytical solution and its finite element approximation in terms of the best approximation error from the finite element spaces. We propose first a Lions–Mercier type iterative method and next a classical fixed-point algorithm to solve the finite-dimensional problems resulting from the finite element discretisation of the system of nonlinear partial differential equations under consideration and present numerical experiments that illustrate the practical performance of the proposed numerical method.
more »
« less
Uniqueness of discrete solutions of nonmonotone PDEs without a globally fine mesh condition
Uniqueness of the finite element solution for nonmonotone quasilinear problems of elliptic type is established in one and two dimensions. In each case, we prove a comparison theorem based on locally bounding the variation of the discrete solution over each element. The uniqueness follows from this result, and does not require a globally small meshsize.
more »
« less
- Award ID(s):
- 1719849
- PAR ID:
- 10061975
- Date Published:
- Journal Name:
- Numerische Mathematik
- ISSN:
- 0029-599X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We develop and analyze a model for the interaction of a quasi-Newtonian free fluid with a poroelastic medium. The flow in the fluid region is described by the nonlinear Stokes equations and in the poroelastic medium by the nonlinear quasi-static Biot model. Equilibrium and kinematic conditions are imposed on the interface. We establish existence and uniqueness of a solution to the weak formulation and its semidiscrete continuous-in-time finite element approximation. We present error analysis, complemented by numerical experiments.more » « less
-
Abstract In this article, we consider nondiffusive variational problems with mixed boundary conditions and (distributional and weak) gradient constraints. The upper bound in the constraint is either a function or a Borel measure, leading to the state space being a Sobolev one or the space of functions of bounded variation. We address existence and uniqueness of the model under low regularity assumptions, and rigorously identify its Fenchel pre-dual problem. The latter in some cases is posed on a nonstandard space of Borel measures with square integrable divergences. We also establish existence and uniqueness of solution to this pre-dual problem under some assumptions. We conclude the article by introducing a mixed finite-element method to solve the primal-dual system. The numerical examples illustrate the theoretical findings.more » « less
-
This work presents inGRASS, a novel algorithm designed for incremental spectral sparsification of large undirected graphs. The proposed inGRASS algorithm is highly scalable and parallel-friendly, having a nearly linear time complexity for the setup phase and the ability to update the spectral sparsifier in O(logN) time for each incremental change made to the original graph with N nodes. A key component in the setup phase of inGRASS is a multilevel resistance embedding framework introduced for efficiently identifying spectrally-critical edges and effectively detecting redundant ones, which is achieved by decomposing the initial sparsifier into many node clusters with bounded effective-resistance diameters leveraging a low-resistance-diameter decomposition (LRD) scheme. The update phase of inGRASS exploits low-dimensional node embedding vectors for efficiently estimating the importance and uniqueness of each newly added edge. As demonstrated through extensive experiments, inGRASS achieves up to over 200× speedups while retaining comparable solution quality in incremental spectral sparsification of graphs obtained from various datasets, such as circuit simulations, finite element analysis, and social networks.more » « less
-
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.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1007/s00211-018-0956-4
