Sudoku’s french ancestors

This column is a place for those bits of contagious mathematics that travel from person to person in the community, because they are so elegant, suprising, or appealing that one has an urge to pass them on. Contributions are most welcome.     

Estimating Invariant Probability Densities for Dynamical Systems

Knowing a probability density (ideally, an invariant density) for the trajectories of a dynamical system allows many significant estimates to be made, from the well-known dynamical invariants such as Lyapunov exponents and mutual information to conditional probabilities which are potentially more suitable for prediction than the single number produced by most predictors. Densities on typical attractors have properties, such as singularity with respect to Lebesgue measure, which make standard density estimators less useful than one would hope. In this paper we present a new method of estimating densities which can smooth in a way that tends to preserve fractal structure down to some level, and that also maintains invariance. We demonstrate with applications to real and artificial data.     

Discrete duality finite volume schemes for Leray−Lions−type elliptic problems on general 2D meshes

Discrete duality finite volume schemes on general meshes, introduced by Hermeline and Domelevo and Omnès for the Laplace equation, are proposed for nonlinear diffusion problems in 2D with nonhomogeneous Dirichlet boundary condition. This approach allows the discretization of non linear fluxes in such a way that the discrete operator inherits the key properties of the continuous one. Furthermore, it is well adapted to very general meshes including the case of nonconformal locally refined meshes. We show that the approximate solution exists and is unique, which is not obvious since the scheme is nonlinear. We prove that, for general W^?1,p′(Ω) source term and W^1‐(1/p),p(?Ω) boundary data, the approximate solution and its discrete gradient converge strongly towards the exact solution and its gradient, respectively, in appropriate Lebesgue spaces. Finally, error estimates are given in the case where the solution is assumed to be in W^2,p(Ω). Numerical examples are given, including those on locally refined meshes. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Geometric and Combinatorial Realizations of Crystal Graphs

For irreducible integrable highest weight modules of the finite and affine Lie algebras of type A and D, we define an isomorphism between the geometric realization of the crystal graphs in terms of irreducible components of Nakajima quiver varieties and the combinatorial realizations in terms of Young tableaux and Young walls. For type A_n⁽¹⁾, we extend the Young wall construction to arbitrary level, describing a combinatorial realization of the crystals in terms of new objects which we call Young pyramids. Presented by P. Littleman Mathematics Subject Classifications (2000) Primary 16G10, 17B37. Alistair Savage: This research was supported in part by the Natural Sciences and Engineering Research Council (NSERC) of Canada and was partially conducted by the author for the Clay Mathematics Institute.     

A New Differential Method Applied to the Study of Arbitrary Cross Section Microstructured Optical Fibers

The present work adapts a recent grating theory called “Fast Fourier factorization” to cylindrical coordinates in order to study microstructured optical fibers (MOFs). Compared with the classical differential method, this new differential method takes into account the truncation of Fourier series and the discontinuities of the fields across the diffracting surface with the help of new factorization rules. The main advantage of this method is that the directrix of the diffracting cylindrical surface is arbitrary and permits anisotropic and inhomogeneous media although its numerical application needs longer computation time, compared with other well-known numerical methods. The S-propagation algorithm is used to avoid numerical contaminations. The numerical results are validated and compared with the well-established Multipole method in the case of a MOF with six circular cylinders. Further, a new cross-sectional profile (with sectorial inclusions) that the Multipole method cannot consider is studied.