Wed, 14 Mar 2012
10:15
OCCAM Common Room (RI2.28)

Meshless methods: from carbon nano-tubes to carbonate reservoir

Alexander Lukyanov
(Schlumberger)
Abstract

In many fields of science and engineering, such as fluid or structural mechanics and electric circuit design, large scale dynamical systems need to be simulated, optimized or controlled. They are often described by discretizations of systems of nonlinear partial differential equations yielding high-dimensional discrete phase spaces. For this reason, in recent decades, research has mainly focused on the development of sophisticated analytical and numerical tools to help understand the overall system behavior. During this time meshless methods have enjoyed significant interest in the research community and in some commercial simulators (e.g., LS-DYNA). In this talk I will describe a normalized-corrected meshless method which ensures linear completeness and improved accuracy. The resulting scheme not only provides first order consistency O(h) but also alleviates the particle deficiency (kernel support incompleteness) problem at the boundary. Furthermore, a number of improvements to the kernel derivative approximation are proposed.

To illustrate the performance of the meshless method, we present results for different problems from various fields of science and engineering (i.e. nano-tubes modelling, solid mechanics, damage mechanics, fluid mechanics, coupled interactions of solids and fluids). Special attention is paid to fluid flow in porous media. The primary attraction of the present approach is that it provides a weak formulation for Darcy's law which can be used in further development of meshless methods.

Tue, 13 Mar 2012

15:45 - 16:45
L3

A Uniqueness Theorem for Gluing Special Lagrangian Submanifolds

Yohsuke Imagi
(Kyoto)
Abstract

Special Lagrangian submanifolds are area minimizing Lagrangian submanifolds discovered by Harvey and Lawson. There is no obstruction to deforming compact special Lagrangian

submanifolds by a theorem of Mclean. It is however difficult to understand singularities of

special Lagrangian submanifolds (varifolds). Joyce has studied isolated singularities with multiplicity one smooth tangent cones. Suppose that there exists a compact special Lagrangian submanifold M of dimension three with one point singularity modelled on the Clliford torus cone. We may apply the gluing technique to M by a theorem of Joyce.

We obtain then a compact non-singular special Lagrangian submanifold sufficiently close to M as varifolds in Geometric Measure Theory. The main result of this talk is as follows: all special Lagrangian varifolds sufficiently close to M are obtained by the gluing technique.

The proof is similar to that of a theorem of Donaldson in the Yang-Mills theory.

One first proves an analogue of Uhlenbeck's removable singularities theorem in the Yang-Mills theory. One uses here an idea of a theorem of Simon, who proved the uniqueness of multiplicity one tangent cones of minimal surfaces. One proves next the uniqueness of local models for desingularizing M (see above) using symmetry of the Clifford torus cone.

These are the main part of the proof.

Tue, 13 Mar 2012

10:15 - 11:15
OCCAM Common Room (RI2.28)

Some symmetry results for the Ginzburg Landau equations

Adriano Pisante
(University of Rome Sapienza)
Abstract

We discuss new symmetry results for nonconstant entire local minimizers of the standard Ginzburg-Landau functional  for maps in ${H}^{1}_{\rm{loc}}(\mathbb{R}^3;\mathbb{R}^3)$ satisfying a natural energy bound.

Up to  translations and rotations, such solutions of the Ginzburg-Landau system are given by an explicit map equivariant under the action of the orthogonal group.

More generally, for any $N\geq 3$ we  characterize the $O(N)-$equivariant vortex solution for Ginzburg-Landau type equations in the $N-$dimensional Euclidean space and we prove its local energy minimality for the corresponding energy functional.

Mon, 12 Mar 2012
12:30
Gibson 1st Floor SR

Computational modeling in high-frequency MEMS resonator design

Sanjay Govindjee
(University of California)
Abstract

In the operation of high frequency resonators in micro electromechanical systems (MEMS)there is a strong need to be able to accurately determine the energy loss rates or alternativelythe quality of the resonance. The resonance quality is directly related to a designer’s abilityto assemble high fidelity system response for signal filtering, for example. This hasimplications on robustness and quality of electronic communication and also stronglyinfluences overall rates of power consumption in such devices – i.e. battery life. Pastdesign work was highly focused on the design of single resonators; this arena of work hasnow given way to active efforts at the design and construction of arrays of coupledresonators. The behavior of such systems in the laboratory shows un-necessarily largespread in operational characteristics, which are thought to be the result of manufacturingvariations. However, such statements are difficult to prove due to a lack of availablemethods for predicting resonator damping – even the single resonator problem is difficult.The physical problem requires the modeling of the behavior of a resonant structure (or setof structures) supported by an elastic half-space. The half-space (chip) serves as a physicalsupport for the structure but also as a path for energy loss. Other loss mechanisms can ofcourse be important but in the regime of interest for us, loss of energy through theanchoring support of the structure to the chip is the dominant effect.

The construction of a basic discretized model of such a system leads to a system ofequations with complex-symmetric (not Hermitian) structure. The complex-symmetryarises from the introduction of a radiation boundary conditions to handle the semi-infinitecharacter of the half-space region. Requirements of physical accuracy dictate rather finediscretization and, thusly, large systems of equations. The core to the extraction of relevantphysical performance parameters is dependent upon the underlying modeling framework.In three dimensional settings of practical interest, such systems are too large to be handleddirectly and must be solved iteratively. In this talk, I will cover the physical background ofthe problem class of interest, how such systems can be modeled, and then solved. Particularinterest will be paid to the radiation boundary conditions (perfectly matched layers versushigher order absorbing boundary conditions), issues associated with frequency domainversus time domain methods, and how these choices interact with iterative solvertechnologies in sometimes unexpected ways. Time permitting I will also touch upon the issue of harmonic inversion methods of this class of problems.

Fri, 09 Mar 2012

14:30 - 15:30
DH 3rd floor SR

Probability Forecasting: Looking Under the Hood and at the Road Ahead

Prof. Leonard A. Smith
(London School of Economics and Pembroke College)
Abstract

Probability does not exist. At least no more so than "mass" "spin" or "charm" exist. Yet probability forecasts are common, and there are fine reasons for deprecating point forecasts, as they require an unscientific certainty in exactly what the future holds. What roles do our physical understanding and laws of physics play in the construction of probability forecasts to support of decision making and science-based policy? Will probability forecasting more likely accelerate or retard the advancement of our scientific understanding?

Model-based probability forecasts can vary significantly with alterations in the method of data assimilation, ensemble formation, ensemble interpretation, and forecast evaluation, not to mention questions of model structure, parameter selection and the available forecast-outcome archive. The role of each of these aspects of forecasting, in the context of interpreting the forecast as a real-world probability, is considered and contrasted in the cases of weather forecasting, climate forecasting, and economic forecasting. The notion of what makes a probability forecast "good" will be discussed, including the goals of "sharpness given calibration" and "value".

For a probability forecast to be decision-relevant as such, it must be reasonably interpreted as a basis for rational action through the reflection of the probability of the outcomes forecast. This rather obvious sounding requirement proves to be the source of major discomfort as the distinct roles of uncertainty (imprecision) and error (structural mathematical "misspecification") are clarified. Probabilistic forecasts can be of value to decision makers even when it is irrational to interpret them as probability forecasts. A similar statement, of course, can be said for point forecasts, or for spin. In this context we explore the question: do decision-relevant probability forecasts exist?

Fri, 09 Mar 2012
14:15
DH 1st floor SR

G-Expectation for General Random Variables

Marcel Nutz
(Columbia)
Abstract

We provide a general construction of time-consistent sublinear expectations on the space of continuous paths. In particular, we construct the conditional G-expectation of a Borel-measurable (rather than quasi-continuous) random variable.

Fri, 09 Mar 2012

11:00 - 13:00
OCCAM Common Room (RI2.28)

OCCAM Group Meeting

Various
Abstract
  • Graham Morris - 'Topics in Voltammetry'
  • James Lottes - 'Algebraic Multigrid for Nonsymmetric Systems'
  • Amy Smith - 'Multi-scale modelling of blood flow in the coronary microcirculation'
Thu, 08 Mar 2012

17:00 - 18:00
L3
Thu, 08 Mar 2012

15:00 - 16:00
L3

Umbral Moonshine

John Duncan
Abstract

In April 2010 Eguchi--Ooguri--Tachikawa observed a fascinating connection between the elliptic genus of a K3 surface and the largest Mathieu group. We will report on joint work with Miranda Cheng and Jeff Harvey that identifies this connection as one component of a system of surprising relationships between a family of finite groups, their representation theory, and automorphic forms of various kinds Mock modular forms, and particularly their shadows, play a key role in the analysis, and we find several of Ramanujan's mock theta functions appearing as McKay--Thompson series arising from the umbral groups.

Thu, 08 Mar 2012

14:00 - 15:00
Gibson Grd floor SR

Solution of ill-posed inverse problems pertaining to signal restoration

Professor Rosie Renaut
(Arizona State University)
Abstract

In this talk I review the use of the spectral decomposition for understanding the solution of ill-posed inverse problems. It is immediate to see that regularization is needed in order to find stable solutions. These solutions, however, do not typically allow reconstruction of signal features such as edges. Generalized regularization assists but is still insufficient and methods of total variation are commonly suggested as an alternative. In the talk I consider application of standard approaches from Tikhonov regularization for finding appropriate regularization parameters in the total variation augmented Lagrangian implementations. Areas for future research will be considered.

Thu, 08 Mar 2012

13:00 - 14:00
L3

Twistor Geometry

Markus Röser
Abstract

Twistor theory is a technology that can be used to translate analytical problems on Euclidean space $\mathbb R^4$ into problems in complex algebraic geometry, where one can use the powerful methods of complex analysis to solve them. In the first half of the talk we will explain the geometry of the Twistor correspondence, which realises $\mathbb R^4$ , or $S^4$, as the space of certain "real" lines in the (projective) Twistor space $\mathbb{CP}^3$. Our discussion will start from scratch and will assume very little background knowledge. As an application, we will discuss the Twistor description of instantons on $S^4$ as certain holomorphic vector bundles on $\mathbb{CP}^3$ due to Ward.

Thu, 08 Mar 2012
13:00
DH 1st floor SR

Pertubative method for quadratic reflected backward stochastic differential equations

Arnaud Lionnet
Abstract

In this talk, I will present reflected backward stochastic differential equations (reflected BSDEs) and their connection with the pricing of American options. Then I will present a simple perturbative method for studying them. Under the appropriate assumptions on the coefficient, the terminal condition and the lower obstacle, similar to those used by Kobylankski, this method allows to prove the existence of a solution. I will also provide the usual comparison theorem and a new proof for a refined comparison theorem, specific to RBSDEs.

Thu, 08 Mar 2012

12:30 - 13:30
Gibson 1st Floor SR

Dynamics for an evolution equation describing micro phase separation

Yoshihito Oshita
(Okayama University)
Abstract

We study the mean-field models describing the evolution of distributions of particle radii obtained by taking the small volume fraction limit of the free boundary problem describing the micro phase separation of diblock copolymer melts, where micro phase separation consists of an ensemble of small balls of one component. In the dilute case, we identify all the steady states and show the convergence of solutions.

Next we study the dynamics for a free boundary problem in two dimension, obtained as a gradient flow of Ohta- Kawasaki free energy, in the case that one component is a distorted disk with a small volume fraction. We show the existence of solutions that a small, almost circular interface moves along a curve determined via a Green’s function of the domain. This talk is partly based on a joint work with Xiaofeng Ren.

Wed, 07 Mar 2012

16:00 - 17:00
SR2

TBA

John MacKay
Wed, 07 Mar 2012
12:30
Gibson 1st Floor SR

Chaos and its frequency in topological dynamical systems

Emma D'Aniello
(Seconda Università degli Studi di Napoli)
Abstract

Let $M$ be the Cantor space or an $n$-dimensional manifold with $C(M,M)$ the set of continuous self-maps of $M$. We analyse the behaviour of the generic $f$ in $C(M,M)$ in terms of attractors and some notions of chaos.

Wed, 07 Mar 2012

10:15 - 11:15
OCCAM Common Room (RI2.28)

The graph realization problem and eigenvector synchronization

Mihai Cucuringu
(Princeton)
Abstract

The graph realization problem has received a great deal of attention in recent years, due to its importance in applications such as wireless sensor networks and structural biology. We introduce the ASAP algorithm, for the graph realization problem in R^d, given a sparse and noisy set of distance measurements associated to the edges of a globally rigid graph. ASAP is a divide and conquer, non-incremental and non-iterative algorithm, which integrates local distance information into a global structure determination. Our approach starts with identifying, for every node, a subgraph of its 1-hop neighborhood graph, which can be accurately embedded in its own coordinate system. In the noise-free case, the computed coordinates of the sensors in each patch must agree with their global positioning up to some unknown rigid motion, that is, up to translation, rotation and possibly reflection. In other words, to every patch there corresponds an element of the Euclidean group Euc(3) of rigid transformations in R^3, and the goal is to estimate the group elements that will properly align all the patches in a globally consistent way. The reflections and rotations are estimated using a recently developed eigenvector synchronization algorithm, while the translations are estimated by solving an overdetermined linear system. Furthermore, the algorithm successfully incorporates information specific to the molecule problem in structural biology, in particular information on known substructures and their orientation. In addition, we also propose SP-ASAP, a faster version of ASAP, which uses a spectral partitioning algorithm as a preprocessing step for dividing the initial graph into smaller subgraphs. Our extensive numerical simulations show that ASAP and SP-ASAP are very robust to high levels of noise in the measured distances and to sparse connectivity in the measurement graph, and compare favorably to similar state-of-the art localization algorithms. Time permitting, we briefly discuss the analogy between the graph realization and the low-rank matrix completion problems, as well as an application of synchronization over Z_2 and its variations to bipartite multislice networks.

Tue, 06 Mar 2012
17:00
L2

Type theories and algebraic theories.

Dr Kobi Kremnitzer
(Oxford)
Abstract

By recent work of Voevodsky and others, type theories are now considered as a candidate

for a homotopical foundations of mathematics. I will explain what are type theories using the language

of (essentially) algebraic theories. This shows that type theories are in the same "family" of algebraic

concepts such as groups and categories. I will also explain what is homotopic in (intensional) type theories.

Tue, 06 Mar 2012

15:45 - 16:45
L3

(HoRSe seminar) Joyce-Song wall-crossing as an asymptotic expansion II

Jacopo Stoppa
(Cambridge)
Abstract

Joyce and Song expressed the wall-crossing behaviour of Donaldson-Thomas invariants using a sum over graphs. Joyce expected that these would have something to do with the Feynman diagrams of suitable physical theories. I will show how this can be achieved in the framework for wall-crossing proposed by Gaiotto, Moore and Neitzke. JS diagrams emerge from small corrections to a hyperkahler metric. The basics of GMN theory will be explained during the first talk.