Weslley da Silva Pereira

Numerical exceptions, which may be caused by overflow, operations like division by 0 or sqrt(−1), or convergence failures, are unavoidable in many cases, in particular when software is used on unforeseen and difficult inputs. As more aspects of society become automated e.g., self-driving cars, health monitors, and cyber-physical systems more generally, it is becoming increasingly important to design software that is resilient to exceptions, and that responds to them in a consistent way. Consistency is needed to allow users to build higher-level software that is also resilient and consistent (and so on recursively). In this paper we explore the design space of consistent exception handling for the widely used BLAS and LAPACK linear algebra libraries, pointing out a variety of instances of inconsistent exception handling in the current versions, and propose a new design that balances consistency, complexity, ease of use, and performance. Some compromises are needed, because there are preexisting inconsistencies that are outside our control, including in or between existing vendor BLAS implementations, different programming languages, and even compilers for the same programming language. And user requests from our surveys are quite diverse. We also propose our design as a possible model for other numerical software, and welcome comments on our design choices.

This work proposes a family of multiscale hybrid-mixed methods for the two-dimensional linear elasticity problem on general polygonal meshes. The new methods approximate displacement, stress, and rotation using two-scale discretizations. The first scale level setting consists of approximating the traction variable (Lagrange multiplier) in discontinuous polynomial spaces, and of computing elementwise rigid body modes. In the second level, the methods are made effective by solving completely independent local boundary Neumann elasticity problems written in a mixed form with weak symmetry enforced via the rotation multiplier. Since the finite-dimensional space for the traction variable constraints the local stress approximations, the discrete stress field lies in the H(div) space globally and stays in local equilibrium with external forces. We propose different choices to approximate local problems based on pairs of finite element spaces defined on affine second-level meshes. Those choices generate the family of multiscale finite element methods for which stability and convergence are proved in a unified framework. Notably, we prove that the methods are optimal and high-order convergent in the natural norms. Also, it emerges that the approximate displacement and stress divergence are super-convergent in the L2-norm. Numerical verifications assess theoretical results and highlight the high precision of the new methods on coarse meshes for multilayered heterogeneous material problems.

The multiscale hybrid-mixed (MHM) method consists of a multi-level strategy to approximate the solution of boundary value problems with heterogeneous coefficients. In this context, we propose a new family of finite elements for the linear elasticity equation defined on coarse polytopal partitions of the domain. The finite elements rely on face degrees of freedom associated with multiscale bases obtained from local Neumann problems with piecewise polynomial interpolations on faces. We establish sufficient conditions on the fine-scale interpolations such that the MHM method is well-posed and optimally convergent under local regularity conditions. Also, a multi-level error analysis demonstrates that the MHM method achieves convergence without refining the coarse partition. The upshot is that the Poincaré and Korn’s inequalities do not degenerate, and then the convergence arises on general meshes. Two- and three-dimensional numerical tests assess theoretical results and verify the robustness of the method on a multi-layer media case.

The injection and in situ generation of foam in porous media successfully control gas mobility and improve the fluids’sweep efficiency inside porous media. Mathematical models describing this problem use two phases, foamed gas, and fluid, and usually have a term for foam generation and destruction. Moreover, the non-Newtonian foam behavior is frequently modeled using Hirasaki and Lawson’s formula for foamed gas viscosity. In this paper, we detail how the traveling wave analysis can be used to estimate the propagation profiles and velocity for a range of non-Newtonian foam models in porous media at constant total superficial flow velocity. We reformulate Hirasaki and Lawson’s formula in an explicit form allowing us to find traveling wave solutions for a foam model with non-Newtonian gas viscosity and a foam generation linearly dependent on the foam texture. Comparing the solution with the one for the Newtonian version allows us to analyze qualitatively and quantitatively the rheology of the foam flow in porous media.