Small block overlapping, and non-overlapping, Schwarz methods are theoretically highly attractive as multilevel smoothers for a wide variety of problems that are not amenable to point relaxation methods.  Examples include monolithic Vanka smoothers for Stokes, overlapping vertex-patch decompositions for $H(\text{div})$ and  $H(\text{curl})$ problems, along with nearly incompressible elasticity, and augmented Lagrangian schemes.

 While it is possible to manually program these different schemes,  their use in general purpose libraries has been held back by a lack   of generic, composable interfaces. We present a new approach to the   specification and development such additive Schwarz methods in PETSc  that cleanly separates the topological space decomposition from the  discretisation and assembly of the equations. Our preconditioner is  flexible enough to support overlapping and non-overlapping additive  Schwarz methods, and can be used to formulate line, and plane smoothers, Vanka iterations, amongst others. I will illustrate these new features with some examples utilising the Firedrake finite element library, in particular how the design of an approriate computational interface enables these schemes to be used as building blocks inside block preconditioners.

This is joint work with Patrick Farrell and Florian Wechsung (Oxford), and Matt Knepley (Buffalo).

Classical methods for X-ray computed tomography (CT) are based on the assumption that the X-ray source intensity is known. In practice, however, the intensity is measured and hence uncertain. Under normal circumstances, when the exposure time is sufficiently high, this kind of uncertainty typically has a negligible effect on the reconstruction quality. However, in time- or dose-limited applications such as dynamic CT, this uncertainty may cause severe and systematic artifacts known as ring artifacts.
By modeling the measurement process and by taking uncertainties into account, it is possible to derive a convex reconstruction model that leads to improved reconstructions when the signal-to-noise ratio is low. We discuss some computational challenges associated with the model and illustrate its merits with some numerical examples based on simulated and real data.

When a liquid drop is deposited over a solid surface whose temperature is sufficiently above the boiling point of the liquid, the drop does not experience nucleate boiling but rather levitates over a thin layer of its own vapor. This is known as the Leidenfrost effect. Whilst highly undesirable in certain cooling applications, because of a drastic decrease of the energy transferred between the solid and the evaporating liquid due to poor heat conductivity of the vapor, this effect can be of great interest in many other processes profiting from this absence of contact with the surface that considerably reduces the friction and confers an extreme mobility on the drop. During this presentation, I hope to provide a good vision of some of the knowledge on this subject through some recent studies that we have done. First, I will present a simple fitting-parameter-free theory of the Leidenfrost effect, successfully validated with experiments, covering the full range of stable shapes, i.e., from small quasi-spherical droplets to larger puddles floating on a pocketlike vapor film. Then, I will discuss the end of life of these drops that appear either to explode or to take-off. Finally, I will show that the Leidenfrost effect can also be observed over hot baths of non-volatile liquids. The understanding of the latter situation, compare to the classical Leidenfrost effect on solid substrate, provides new insights on the phenomenon, whether it concerns levitation or its threshold.

We discuss the key role that bespoke linear algebra plays in modelling PDEs with random coefficients using stochastic Galerkin approximation methods. As a specific example, we consider nearly incompressible linear elasticity problems with an uncertain spatially varying Young's modulus. The uncertainty is modelled with a finite set of parameters with prescribed probability distribution.  We introduce a novel three-field mixed variational formulation of the PDE model and and  assess the stability with respect to a weighted norm. The main focus will be  the efficient solution of the associated high-dimensional indefinite linear system of equations. Eigenvalue bounds for the preconditioned system can be  established and shown to be independent of the discretisation parameters and the Poisson ratio.  We also  discuss an associated a posteriori error estimation strategy and assess proxies for the error reduction associated with selected enrichments of the approximation spaces.  We will show by example that these proxies enable the design of efficient  adaptive solution algorithms that terminate when the estimated error falls below a user-prescribed tolerance.

This is joint work with Arbaz Khan and Catherine Powell

Just like optimization needs derivatives, shape optimization needs shape derivatives. Their definition and computation is a classical subject, at least concerning first order shape derivatives. Second derivatives have been studied as well, but some aspects of their theory still remains a bit mysterious for practitioners. As a result, most algorithms for shape optimization are first order methods.

To understand this situation better and in a general way, we consider first and second order shape sensitivities of integrals on smooth submanifolds using a variant of shape differentiation. Instead of computing the second derivative as the derivative of the first derivative, we choose a one-parameter family of perturbations  and compute first and second derivatives with respect to that parameter. The result is a  quadratic form in terms of a perturbation vector field that yields a second order quadratic model of the perturbed functional, which can be used as the basis of a second order shape optimization algorithm. We discuss the structure of this derivative, derive domain expressions and Hadamard forms in a general geometric framework, and give a detailed geometric interpretation of the arising terms.

Finally, we use our results to construct a second order SQP-algorithm for shape optimization that exhibits indeed local fast convergence.

Analytic continuation is ill-posed, but becomes merely ill-conditioned (though with an infinite condition number) if it is known that the function in question is bounded in a given region of the complex plane.
This classical, seemingly theoretical subject has many connections with numerical practice.  One argument indicates that if one tracks an analytic function from z=1 around a branch point at z=0 and back to z=1 again by a Weierstrass chain of disks, the number of accurate digits is divided by about exp(2 pi e) ~= 26,000,000.

Likely instabilities in stochastic hyperelastic solids

L. Angela Mihai

School of Mathematics, Cardiff University, Senghennydd Road, Cardiff, CF24 4AG, UK

E-mail: @email.uk

Nonlinear elasticity has been an active topic of fundamental and applied research for several decades. However, despite numerous developments and considerable attention it has received, there are important issues that remain unresolved, and many aspects still elude us. In particular, the quantification of uncertainties in material parameters and responses resulting from incomplete information remain largely unexplored. Nowadays, it is becoming increasingly apparent that deterministic approaches, which are based on average data values, can greatly underestimate, or overestimate, mechanical properties of many materials. Thus, stochastic representations, accounting for data dispersion, are needed to improve assessment and predictions. In this talk, I will consider stochastic hyperelastic material models described by a strain-energy density where the parameters are characterised by probability distributions. These models, which are constructed through a Bayesian identification procedure, rely on the maximum entropy principle and enable the propagation of uncertainties from input data to output quantities of interest. Similar modelling approaches can be developed for other mechanical systems. To demonstrate the effect of probabilistic model parameters on large strain elastic responses, specific case studies include the classic problem of the Rivlin cube, the radial oscillatory motion of cylindrical and spherical shells, and the cavitation and finite amplitude oscillations of spheres.