Adaptive multiscale model reduction with Generalized Multiscale Finite Element Methods

In this paper, we discuss a general multiscale model reduction framework based on multiscale finite element methods. We give a brief overview of related multiscale methods. Due to page limitations, the overview focuses on a few related methods and is not intended to be comprehensive. We present a general adaptive multiscale model reduction framework, the Generalized Multiscale Finite Element Method. Besides the method's basic outline, we discuss some important ingredients needed for the method's success. We also discuss several applications. The proposed method allows performing local model reduction in the presence of high contrast and no scale separation

Polyharmonic homogenization, rough polyharmonic splines and sparse super-localization

We introduce a new variational method for the numerical homogenization of divergence form elliptic, parabolic and hyperbolic equations with arbitrary rough ($L^\infty$) coefficients. Our method does not rely on concepts of ergodicity or scale-separation but on compactness properties of the solution space and a new variational approach to homogenization. The approximation space is generated by an interpolation basis (over scattered points forming a mesh of resolution $H$) minimizing the $L^2$ norm of the source terms; its (pre-)computation involves minimizing $\mathcal{O}(H^{-d})$ quadratic (cell) problems on (super-)localized sub-domains of size $\mathcal{O}(H \ln (1/ H))$. The resulting localized linear systems remain sparse and banded. The resulting interpolation basis functions are biharmonic for $d\leq 3$, and polyharmonic for $d\geq 4$, for the operator -\diiv(a\nabla \cdot) and can be seen as a generalization of polyharmonic splines to differential operators with arbitrary rough coefficients. The accuracy of the method ($\mathcal{O}(H)$ in energy norm and independent from aspect ratios of the mesh formed by the scattered points) is established via the introduction of a new class of higher-order Poincar\'{e} inequalities. The method bypasses (pre-)computations on the full domain and naturally generalizes to time dependent problems, it also provides a natural solution to the inverse problem of recovering the solution of a divergence form elliptic equation from a finite number of point measurements.Comment: ESAIM: Mathematical Modelling and Numerical Analysis. Special issue (2013

A Stochastic Multi-scale Approach for Numerical Modeling of Complex Materials - Application to Uniaxial Cyclic Response of Concrete

In complex materials, numerous intertwined phenomena underlie the overall response at macroscale. These phenomena can pertain to different engineering fields (mechanical , chemical, electrical), occur at different scales, can appear as uncertain, and are nonlinear. Interacting with complex materials thus calls for developing nonlinear computational approaches where multi-scale techniques that grasp key phenomena at the relevant scale need to be mingled with stochastic methods accounting for uncertainties. In this chapter, we develop such a computational approach for modeling the mechanical response of a representative volume of concrete in uniaxial cyclic loading. A mesoscale is defined such that it represents an equivalent heterogeneous medium: nonlinear local response is modeled in the framework of Thermodynamics with Internal Variables; spatial variability of the local response is represented by correlated random vector fields generated with the Spectral Representation Method. Macroscale response is recovered through standard ho-mogenization procedure from Micromechanics and shows salient features of the uniaxial cyclic response of concrete that are not explicitly modeled at mesoscale.Comment: Computational Methods for Solids and Fluids, 41, Springer International Publishing, pp.123-160, 2016, Computational Methods in Applied Sciences, 978-3-319-27994-

Computational Multiscale Methods

Many physical processes in material sciences or geophysics are characterized by inherently complex interactions across a large range of non-separable scales in space and time. The resolution of all features on all scales in a computer simulation easily exceeds today's computing resources by multiple orders of magnitude. The observation and prediction of physical phenomena from multiscale models, hence, requires insightful numerical multiscale techniques to adaptively select relevant scales and effectively represent unresolved scales. This workshop enhanced the development of such methods and the mathematics behind them so that the reliable and efficient numerical simulation of some challenging multiscale problems eventually becomes feasible in high performance computing environments

Invisibility and Inverse Problems

This survey of recent developments in cloaking and transformation optics is an expanded version of the lecture by Gunther Uhlmann at the 2008 Annual Meeting of the American Mathematical Society.Comment: 68 pages, 12 figures. To appear in the Bulletin of the AM