For the purpose of uncertainty quantification with collocation, a method is proposed for generating families of one-dimensional nested quadrature rules with positive weights and symmetric nodes. This is achieved through a reduction procedure: we start with a high-degree quadrature rule with positive weights and remove nodes while preserving symmetry and positivity. This is shown to be always possible, by a lemma depending primarily on Carathéodory's theorem. The resulting one-dimensional rules can be used within a Smolyak procedure to produce sparse multi-dimensional rules, but weight positivity is lost then. As a remedy, the reduction procedure is directly applied to multi-dimensional tensor-product cubature rules. This allows to produce a family of sparse cubature rules with positive weights, competitive with Smolyak rules. Finally the positivity constraint is relaxed to allow more flexibility in the removal of nodes. This gives a second family of sparse cubature rules, in which iteratively as many nodes as possible are removed. The new quadrature and cubature rules are applied to test problems from mathematics and fluid dynamics. Their performance is compared with that of the tensor-product and standard Clenshaw–Curtis Smolyak cubature rule.

In this work we present a robust and accurate arbitrary order solver for the fixed-boundary plasma equilibria in toroidally axisymmetric geometries. To achieve this we apply the mimetic spectral element formulation presented in [56] to the solution of the Grad-Shafranov equation. This approach combines a finite volume discretization with the mixed finite element method. In this way the discrete differential operators (∇, ∇×, ∇.) can be represented exactly and metric and all approximation errors are present in the constitutive relations. The result of this formulation is an arbitrary order method even on highly curved meshes. Additionally, the integral of the toroidal current J_φ is exactly equal to the boundary integral of the poloidal field over the plasma boundary. This property can play an important role in the coupling between equilibrium and transport solvers. The proposed solver is tested on a varied set of plasma cross sections (smooth and with an X-point) and also for a wide range of pressure and toroidal magnetic flux profiles. Equilibria accurate up to machine precision are obtained. Optimal algebraic convergence rates of order p+1 and geometric convergence rates are shown for Soloviev solutions (including high Shafranov shifts), field-reversed configuration (FRC) solutions and spheromak analytical solutions. The robustness of the method is demonstrated for non-linear test cases, in particular on an equilibrium solution with a pressure pedestal.