The ability of cells to sense and respond to the mechanical properties of their environments is fundamental to cellular behaviour, with stiffness found to be a key control parameter. The physical mechanisms underpinning mechanosensing are, however, not well understood. I here consider the key physical cellular behaviours of active contractility of the internal cytoskeleton and cell growth, coupling these into mechanical models. These models suggest new distinct mechanisms of mechanotransduction in cells and tissues.

I will discuss a new theoretical approach to information and decisions in signalling systems and relate this to new experimental results about the NF-kappaB signalling system. NF-kappaB is an exemplar system that controls inflammation and in different contexts has varying effects on cell death and cell division. It is commonly claimed that it is information processing hub, taking in signals about the infection and stress status of the tissue environment and as a consequence of the oscillations, transmitting higher amounts of information to the hundreds of genes it controls. My aim is to develop a conceptual and mathematical framework to enable a rigorous quantifiable discussion of information in this context in order to follow Francis Crick's counsel that it is better in biology to follow the flow of information than those of matter or energy. In my approach the value of the information in the signalling system is defined by how well it can be used to make the "correct decisions" when those "decisions" are made by molecular networks. As part of this I will introduce a new mathematical method for the analysis and simulation of large stochastic non-linear oscillating systems. This allows an analytic analysis of the stochastic relationship between input and response and shows that for tightly-coupled systems like those based on current models for signalling systems, clocks, and the cell cycle this relationship is highly constrained and non-generic.

Multivariate cryptography is one of a handful of proposals for post-quantum cryptographic schemes, i.e. cryptographic schemes that are secure also against attacks carried on with a quantum computer. Their security relies on the assumption that solving a system of multivariate (quadratic) equations over a finite field is computationally hard. 

Groebner bases allow us to solve systems of polynomial equations. Therefore, one of the key questions in assessing the robustness of multivariate cryptosystems is estimating how long it takes to compute the Groebner basis of a given system of polynomial equations. 

After introducing multivariate cryptography and Groebner bases, I will present a rigorous method to estimate the complexity of computing a Groebner basis. This approach is based on techniques from commutative algebra and is joint work with Alessio Caminata (University of Barcelona).

In this talk I will review recent results on the analysis of shock-capturing-type methods applied to convection-dominated problems. The method of choice is a variant of the Algebraic Flux-Correction (AFC) scheme. This scheme has received some attention over the last two decades due to its very satisfactory numerical performance. Despite this attention, until very recently there was no stability and convergence analysis for it. Thus, the purpose of the works reviewed in this talk was to bridge that gap. The first step towards the full analysis of the method is a rewriting of it as a nonlinear edge-based diffusion method. This writing makes it possible to present a unified analysis of the different variants of it. So, minimal assumptions on the components of the method are stated in such a way that the resulting scheme satisfies the Discrete Maximum Principle (DMP) and is convergence. One property that will be discussed in detail is the linearity preservation. This property has been linked to the good performance of methods of this kind. We will discuss in detail its role and the impact of it in the overall convergence of the method. Time permitting, some results on a posteriori error estimation will also be presented. 
This talk will gather contributions with A. Allendes (UTFSM, Chile), E. Burman (UCL, UK), V. John (WIAS, Berlin), F. Karakatsani (Chester, UK), P. Knobloch (Prague, Czech Republic), and 
R. Rankin (U. of Nottingham, China).

The development of reduced order models for complex applications, offering the promise for rapid and accurate evaluation of the output of complex models under parameterized variation, remains a very active research area. Applications are found in problems which require many evaluations, sampled over a potentially large parameter space, such as in optimization, control, uncertainty quantification and applications where near real-time response is needed.

However, many challenges remain to secure the flexibility, robustness, and efficiency needed for general large-scale applications, in particular for nonlinear and/or time-dependent problems.

After giving a brief general introduction to reduced order models, we discuss developments in two different directions. In the first part, we discuss recent developments of reduced methods that conserve chosen invariants for nonlinear time-dependent problems. We pay particular attention to the development of reduced models for Hamiltonian problems and propose a greedy approach to build the basis. As we shall demonstrate, attention to the construction of the basis must be paid not only to ensure accuracy but also to ensure stability of the reduced model. Time permitting, we shall also briefly discuss how to extend the approach to include more general dissipative problems through the notion of port-Hamiltonians, resulting in reduced models that remain stable even in the limit of vanishing viscosity and also touch on extensions to Euler and Navier-Stokes equations.

The second part of the talk discusses the combination of reduced order modeling for nonlinear problems with the use of neural networks to overcome known problems of on-line efficiency for general nonlinear problems. We discuss the general idea in which training of the neural network becomes part of the offline part and demonstrate its potential through a number of examples, including for the incompressible Navier-Stokes equations with geometric variations.

This work has been done with in collaboration with B.F. Afkram (EPFL, CH), N. Ripamonti EPFL, CH) and S. Ubbiali (USI, CH).

The rapidly increasing number of cores in high-performance computing systems causes a multitude of challenges for developers of numerical methods. New parallel algorithms are required to unlock future growth in computing power for applications and energy efficiency and algorithm-based fault tolerance are becoming increasingly important. So far, most approaches to parallelise the numerical solution of partial differential equations focussed on spatial solvers, leaving time as a bottleneck. Recently, however, time stepping methods that offer some degree of concurrency, so-called parallel-in-time integration methods, have started to receive more attention.

I will introduce two different numerical algorithms, Parareal (by Lions et al., 2001) and PFASST (by Emmett and Minion, 2012), that allow to exploit concurrency along the time dimension in parallel computer simulations solving partial differential equations. Performance results for both methods on different architectures and for different equations will be presented. The PFASST algorithm is based on merging ideas from Parareal, spectral deferred corrections (SDC, an iterative approach to derive high-order time stepping methods by Dutt et al. 2000) and nonlinear multi-grid. Performance results for PFASST on close to half a million cores will illustrate the potential of the approach. Algorithmic modifications like IPFASST will be introduced that can further reduce solution times. Also, recent results showing how parallel-in-time integration can provide algorithm-based tolerance against hardware faults will be shown.

In micro-fluidics, at small scales where inertial effects become negligible, surface to volume ratios are large and the interfacial processes are extremely important for the overall dynamics. Integral
equation based methods are attractive for the simulations of e.g. droplet-based microfluidics, with tiny water drops dispersed in oil, stabilized by surfactants. In boundary integral formulations for
Stokes flow, jumps in pressure and velocity gradients are naturally taken care of, viscosity ratios enter only in coefficients of the equations, and only the drop surfaces must be discretized and not the volume inside nor in between.

We present numerical methods for drops with insoluble surfactants, both in two and three dimensions. We discretize the integral equations using Nyström methods, and special care is taken in the evaluation of singular and also nearly singular integrals that is needed in the case of close drop interactions. A spectral method is used to solve the advection-diffusion equation on each drop surface that describes the evolution of surfactant concentration. The drop velocity and surfactant concentration couple together through an equation of state for the surface tension coefficient. An adaptive time-stepping strategy is developed for the coupled problem, with the constraint to minimize the number of Stokes solves, since this is the computationally most expensive part.

For high quality discretization of the drops throughout the simulations, a hybrid method is used in two dimensions, offering an arc-length parameterization of the interface. In three dimensions, a
reparameterization procedure is developed to optimize the spherical harmonics representation of the drop, while conserving the drop volume and amount of surfactant.

We present results from some validation tests and illustrate the ability of the numerical methods in different challenging problems.