Abstract

We present a method to blend local maximum entropy (LME) meshfree approximants and isogeometric analysis. The coupling strategy exploits the optimization program behind LME approximation, treats isogeometric and LME basis functions on an equal footing in the reproducibility constraints, but views the former as data in the constrained minimization. The resulting scheme exploits the best features and overcomes the main drawbacks of each of these approximants. Indeed, it preserves the high fidelity boundary representation (exact CAD geometry) of isogeometric analysis, out of reach for bare meshfree methods, and easily handles volume discretization and unstructured grids with possibly local refinement, while maintaining the smoothness and non-negativity of the basis functions. We implement the method with B-Splines in two dimensions, but the procedure carries over to higher spatial dimensions or to other non-negative approximants such as NURBS or subdivision schemes. The performance of the method is illustrated with the heat equation, and linear and nonlinear elasticity. The ability of the proposed method to impose directly essential boundary conditions in non-convex domains, and to deal with unstructured grids and local refinement in domains of complex geometry and topology is highlighted by the numerical examples.