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.     

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     

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.