Scalable numerical approach for the steady-state ab initio laser theory

We present an efficient and flexible method for solving the non-linear lasing equations of the steady-state ab initio laser theory. Our strategy is to solve the underlying system of partial differential equations directly, without the need of setting up a parametrized basis of constant flux states. We validate this approach in one-dimensional as well as in cylindrical systems, and demonstrate its scalability to full-vector three-dimensional calculations in photonic-crystal slabs. Our method paves the way for efficient and accurate simulations of lasing structures which were previously inaccessible.Comment: 17 pages, 8 figure

Schnelle Löser für Partielle Differentialgleichungen

The workshop Schnelle Löser für partielle Differentialgleichungen, organised by Randolph E. Bank (La Jolla), Wolfgang Hackbusch (Leipzig), and Gabriel Wittum (Frankfurt am Main), was held May 22nd–May 28th, 2011. This meeting was well attended by 54 participants with broad geographic representation from 7 countries and 3 continents. This workshop was a nice blend of researchers with various backgrounds

Multigrid Discretization and Iterative Algorithm for Mixed Variational Formulation of the Eigenvalue Problem of Electric Field

This paper discusses highly finite element algorithms for the eigenvalue problem of electric field. Combining the mixed finite element method with the Rayleigh quotient iteration method, a new multi-grid discretization scheme and an adaptive algorithm are proposed and applied to the eigenvalue problem of electric field. Theoretical analysis and numerical results show that the computational schemes established in the paper have high efficiency

Recovery Techniques For Finite Element Methods And Their Applications

Recovery techniques are important post-processing methods to obtain improved approximate solutions from primary data with reasonable cost. The practical us- age of recovery techniques is not only to improve the quality of approximation, but also to provide an asymptotically exact posteriori error estimators for adaptive meth- ods. This dissertation presents recovery techniques for nonconforming finite element methods and high order derivative as well as applications of gradient recovery. Our first target is to develop a systematic gradient recovery technique for Crouzeix- Raviart element. The proposed method uses finite element solution to build a better approximation of the exact gradient based on local least square fittings. Due to poly- nomial preserving property of least square fitting, it is easy to show that the new proposed method preserves quadratic polynomials. In addition, the proposed gra- dient recovery is linearly bounded. Numerical tests indicate the recovered gradient is superconvergent to the exact gradient for both second order elliptic equation and Stokes equation. The gradient recovery technique can be used in a posteriori error estimates for Crouzeix-Raviart element, which is relatively simple to implement and problem independent. Our second target is to propose and analyze a new effective Hessian recovery for continuous finite element of arbitrary order. The proposed Hessian recovery is based on polynomial preserving recovery. The proposed method preserves polynomials of degree (k + 1) on general unstructured meshes and polynomials of degree (k + 2) on translation invariant meshes. Based on it polynomial preserving property, we can able to prove superconvergence of the proposed method on mildly structured meshes. In addition, we establish the ultraconvergence result for the new Hessian recovery technique on translation invariant finite element space of arbitrary order. Our third target is to demonstrate application of gradient recovery in eigenvalue computation. We propose two superconvergent two-grid methods for elliptic eigen- value problems by taking advantage of two-gird method, two-space method, shifted- inverse power method, and gradient recovery enhancement. Theoretical and numer- ical results reveal that the proposed methods provide superconvergent eigenfunction approximation and ultraconvergent eigenvalue approximation. In addition, two mul- tilevel adaptive methods based recovery type a posterior error estimate are proposed