Search

In this paper, we present port-Hamiltonian formulations of the incompressible Euler equations with a free surface governed by surface tension and gravity forces, modelling e.g. capillary and gravity waves and the evolution of droplets in air. Three sets of variables are considered, namely $(v,\Sigma)$, $(\eta,\phi_{\partial},\Sigma)$ and $(\omega,\phi_{\partial},\Sigma)$, with $v$ the velocity, $\eta$ the solenoidal velocity, $\phi_{\partial}$ a potential, $\omega$ the vorticity, and $\Sigma$ the free surface, resulting in the incompressible Euler equations in primitive variables and the vorticity equation. First, the Hamiltonian formulation for the incompressible Euler equations in a domain with a free surface combined with a fixed boundary surface with a homogeneous boundary condition will be derived in the proper Sobolev spaces of differential forms. Next, these results will be extended to port-Hamiltonian formulations allowing inhomogeneous boundary conditions and a non-zero energy flow through the boundaries. Our main results are the construction and proof of Dirac structures in suitable Sobolev spaces of differential forms for each variable set, which provides the core of any port-Hamiltonian formulation. Finally, it is proven that the state dependent Dirac structures are related to Poisson brackets that are linear, skew-symmetric and satisfy the Jacobi identity.

We develop entropy dissipative higher order accurate local discontinuous Galerkin (LDG) discretizations coupled with Diagonally Implicit Runge-Kutta (DIRK) methods for nonlinear degenerate parabolic equations with a gradient flow structure. Using the simple alternating numerical flux, we construct DIRK-LDG discretizations that combine the advantages of higher order accuracy, entropy dissipation and proper long-time behavior. The implicit time-discrete methods greatly alleviate the time-step restrictions needed for the stability of the numerical discretizations. Also, the larger time step significantly improves computational efficiency. We theoretically prove the unconditional entropy dissipation of the implicit Euler-LDG discretization. Next, in order to ensure the positivity of the numerical solution, we use the Karush-Kuhn-Tucker (KKT) limiter, which couples the positivity inequality constraint with higher order accurate DIRK-LDG discretizations using Lagrange multipliers. In addition, mass conservation of the positivity-limited solution is ensured by imposing a mass conservation equality constraint to the KKT equations. The unique solvability and unconditional entropy dissipation for an implicit first order accurate in time, but higher order accurate in space, KKT-LDG discretizations are proved, which provides a first theoretical analysis of the KKT limiter. Finally, numerical results demonstrate the higher order accuracy and entropy dissipation of the KKT-DIRK-LDG discretizations for problems requiring a positivity limiter.

A port-Hamiltonian (pH) system formulation is a geometrical notion used to formulate conservation laws for various physical systems. The distributed parameter port-Hamiltonian formulation models infinite dimensional Hamiltonian dynamical systems that have a non-zero energy flow through the boundaries. In this paper we propose a novel framework for discontinuous Galerkin (DG) discretizations of pH-systems. Linking DG methods with pH-systems gives rise to compatible structure preserving finite element discretizations along with flexibility in terms of geometry and function spaces of the variables involved. Moreover, the port-Hamiltonian formulation makes boundary ports explicit, which makes the choice of structure and power preserving numerical fluxes easier. We state the Discontinuous Finite Element Stokes-Dirac structure with a power preserving coupling between elements, which provides the mathematical framework for a large class of pH discontinuous Galerkin discretizations. We also provide an a priori error analysis for the port-Hamiltonian discontinuous Galerkin Finite Element Method (pH-DGFEM). The port-Hamiltonian discontinuous Galerkin finite element method is demonstrated for the scalar wave equation showing optimal rates of convergence.

Since thin-film silicon solar cells have limited optical absorption, we explore the effect of a nanostructured back reflector to recycle the unabsorbed light. As a back reflector we investigate a 3D photonic band gap crystal made from silicon that is readily integrated with the thin films. We numerically obtain the optical properties by solving the 3D time-harmonic Maxwell equations using the finite-element method, and model silicon with experimentally determined optical constants. The absorption enhancement relevant for photovoltaics is obtained by weighting the absorption spectra with the AM 1.5 standard solar spectrum. We study thin films either thicker ($L_{Si} = 2400$ nm) or much thinner ($L_{Si} = 80$ nm) than the wavelength of light. At $L_{Si} = 2400$ nm, the 3D photonic band gap crystal enhances the spectrally averaged ($\lambda = 680$ nm to $880$ nm) silicon absorption by $2.22$x (s-pol.) to $2.45$x (p-pol.), which exceeds the enhancement of a perfect metal back reflector ($1.47$ to $1.56$x). The absorption is enhanced by the (i) broadband angle and polarization-independent reflectivity in the 3D photonic band gap, and (ii) the excitation of many guided modes in the film by the crystal's surface diffraction leading to enhanced path lengths. At $L_{Si} = 80$ nm, the photonic crystal back reflector yields a striking average absorption enhancement of $9.15$x, much more than $0.83$x for a perfect metal, which is due to a remarkable guided mode confined within the combined thickness of the thin film and the photonic crystal's Bragg attenuation length. The broad bandwidth of the 3D photonic band gap leads to the back reflector's Bragg attenuation length being much shorter than the silicon absorption length. Consequently, light is confined inside the thin film and the absorption enhancements are not due to the additional thickness of the photonic crystal back reflector., Comment: 23 pages, 15 figures

An error analysis of a mixed discontinuous Galerkin (DG) method with Brezzi numerical flux for the time-harmonic Maxwell equations with minimal smoothness requirements is presented. The key difficulty in the error analysis for the DG method is that the tangential or normal trace of the exact solution is not well-defined on the mesh faces of the computational mesh. We overcome this difficulty by two steps. First, we employ a lifting operator to replace the integrals of the tangential/normal traces on mesh faces by volume integrals. Second, optimal convergence rates are proven by using smoothed interpolations that are well-defined for merely integrable functions. As a byproduct of our analysis, an explicit and easily computable stabilization parameter is given.

We present new and efficient quadrature rules for computing the stiffness matrices of mass-lumped tetrahedral elements for wave propagation modelling. These quadrature rules allow for a more efficient implementation of the mass-lumped finite element method and can handle materials that are heterogeneous within the element without loss of the convergence rate. The quadrature rules are designed for the specific function spaces of recently developed mass-lumped tetrahedra, which consist of standard polynomial function spaces enriched with higher-degree bubble functions. For the degree-2 mass-lumped tetrahedron, the most efficient quadrature rule seems to be an existing 14-point quadrature rule, but for tetrahedra of degrees 3 and 4, we construct new quadrature rules that require less integration points than those currently available in the literature. Several numerical examples confirm that this approach is more efficient than computing the stiffness matrix exactly and that an optimal order of convergence is maintained, even when material properties vary within the element.

Currently, nearly all positivity preserving discontinuous Galerkin (DG) discretizations of partial differential equations are coupled with explicit time integration methods. Unfortunately, for many problems this can result in severe time-step restrictions. The techniques used to develop explicit positivity preserving DG discretizations can, however, not easily be combined with implicit time integration methods. In this paper we therefore present a new approach. Using Lagrange multipliers the conditions imposed by the positivity preserving limiters are directly coupled to a DG discretization combined with a Diagonally Implicit Runge-Kutta time integration method. The positivity preserving DG discretization is then reformulated as a Karush-Kuhn-Tucker (KKT) problem, which is frequently encountered in constrained optimization. Since the limiter is only active in areas where positivity must be enforced it does not affect the higher order DG discretization elsewhere. The resulting non-smooth nonlinear algebraic equations have, however, a different structure compared to most constrained optimization problems. We therefore develop an efficient active set semi-smooth Newton method that is suitable for the KKT formulation of time-implicit positivity preserving DG discretizations. Convergence of this semi-smooth Newton method is proven using a specially designed quasi-directional derivative of the time-implicit positivity preserving DG discretization. The time-implicit positivity preserving DG discretization is demonstrated for several nonlinear scalar conservation laws, which include the advection, Burgers, Allen-Cahn, Barenblatt, and Buckley-Leverett equations., Comment: 23 Pages

We study numerically the transport and storage of light in a 3D photonic band gap crystal doped by a single embedded resonant cavity. The crystal has finite support since it is surrounded by vacuum, as in experiments. Therefore, we employ the finite element method to model the diamond-like inverse woodpile crystal that consists of two orthogonal arrays of pores in a high-index dielectric such as silicon and that has experimentally been realized by CMOS-compatible methods. A point defect that functions as the resonant cavity is formed in the proximal region of two selected orthogonal pores with a radius smaller than the ones in the bulk of the crystal. We present a field-field cross-correlation method to identify resonances in the finite-support crystal with defect states that appear in the 3D photonic band gap of the infinite crystal. Out of 5 observed angle-independent cavity resonances, one is s-polarized and 4 are p-polarized for light incident in the X or Z directions. It is remarkable that quality factors up to Q = 1000 appear in such thin structures (3 unit cells). We find that the optical energy density is remarkably enhanced at the cavity resonances by up to 2400x the incident energy density in vacuum or up to 1200x the energy density of the equivalent effective medium. We find that an inverse woodpile photonic band gap cavity with a suitably adapted lattice parameter reveals substantial absorption in the visible range. Below the 3D photonic band gap, Fano resonances arise due to interference between the discrete fundamental cavity mode and the continuum light scattered by the photonic crystal. We argue that the five eigenstates of our 3D photonic band gap cavity have quadrupolar symmetry, in analogy to d-like orbitals of transition metals. We conclude that inverse woodpile cavities offer interesting perspectives for applications in optical sensing and photovoltaics., Comment: 18 pages, 19 figures

We present a new accuracy condition for the construction of continuous mass-lumped elements. This condition is less restrictive than the one currently used and enabled us to construct new mass-lumped tetrahedral elements of degrees 2 to 4. The new degree-2 and degree-3 tetrahedral elements require 15 and 32 nodes per element, respectively, while currently, these elements require 23 and 50 nodes, respectively. The new degree-4 elements require 60, 61 or 65 nodes per element. Tetrahedral elements of this degree had not been found yet. We prove that our accuracy condition results in a mass-lumped finite element method that converges with optimal order in the $L^2$-norm and energy-norm. A dispersion analysis and several numerical tests confirm that our elements maintain the optimal order of accuracy and show that the new mass-lumped tetrahedral elements are more efficient than the current ones.

We analyse the dispersion properties of two types of explicit finite element methods for modelling acoustic and elastic wave propagation on tetrahedral meshes, namely mass-lumped finite element methods and symmetric interior penalty discontinuous Galerkin methods, both combined with a suitable Lax--Wendroff time integration scheme. The dispersion properties are obtained semi-analytically using standard Fourier analysis. Based on the dispersion analysis, we give an indication of which method is the most efficient for a given accuracy, how many elements per wavelength are required for a given accuracy, and how sensitive the accuracy of the method is to poorly shaped elements., Comment: J Sci Comput (2018)

We present sharp and sufficient bounds for the interior penalty term and time step size to ensure stability of the Symmetric Interior Penalty Discontinuous Galerkin (SIPDG) method combined with an explicit time-stepping scheme. These conditions hold for generic meshes, including unstructured non-conforming heterogeneous meshes of mixed element types, and apply to a large class of linear hyperbolic problems, including the acoustic wave equation, the (an)isotropic elastic wave equations and Maxwell's equations. The penalty term bounds are computed element-wise, while bounds for the time step size are computed at weighted submeshes requiring only a small number of elements and faces. Numerical results illustrate the sharpness of these bounds.

We present and analyze a new space-time finite element method for the solution of neural field equations with transmission delays. The numerical treatment of these systems is rare in the literature and currently has several restrictions on the spatial domain and the functions involved, such as connectivity and delay functions. The use of a space-time discretization, with basis functions that are discontinuous in time and continuous in space (dGcG-FEM), is a natural way to deal with space-dependent delays, which is important for many neural field applications. In this article we provide a detailed description of a space-time dGcG-FEM algorithm for neural delay equations, including an a-priori error analysis. We demonstrate the application of the dGcG-FEM algorithm on several neural field models, including problems with an inhomogeneous kernel.

We study numerically the reflectivity of three-dimensional (3D) photonic crystals with a complete 3D photonic band gap, with the aim to interpret recent experiments. We employ the finite element method to study crystals with the cubic diamond-like inverse woodpile structure. The high-index backbone has a dielectric function similar to silicon. We study crystals with a range of thicknesses up to ten unit cells ($L \leq 10 c$). The crystals are surrounded by vacuum, and have a finite support as in experiments. The polarization-resolved reflectivity spectra reveal Fabry-P{\'e}rot fringes related to standing waves in the finite crystal, as well as broad stop bands with nearly $100~\%$ reflectivity, even for thin crystals. From the strong reflectivity peaks, it is inferred that the maximum reflectivity observed in experiments is not limited by finite size. The frequency ranges of the stop bands are in excellent agreement with stop gaps in the photonic band structure, that pertain to infinite and perfect crystals. The frequency ranges of the observed stop bands hardly change with angle of incidence, which is plausible since the stop bands are part of the 3D band gap. Moreover, this result supports the previous assertion that intense reflection peaks measured with a large numerical aperture provide a faithful signature of the 3D photonic band gap. The Bragg attenuation lengths $L_{B}$ exceed the earlier estimates based on the width of the stop band by a factor $6$ to $9$. Hence crystals with a thickness of $12$ unit cells studied in experiments are in the thick crystal limit ($L >> L_{B}$). In our calculations for p-polarized waves, we also observe an intriguing hybridization of the zero reflection of Fabry-P{\'e}rot fringes and the Brewster angle, which has not yet been observed in experiments.

Two fourth order accurate Osher schemes are presented which maintain higher order accuracy on nonuniform grids. They use either a conservative finite difference or finite volume discretization. Both methods are successfully used for direct numerical simulations of flat plate boundary layer instability at different Mach numbers. Results of growth rates of Tollmien-Schlichting waves compare well with direct simulations of incompressible flow and for compressible flow with results obtained by solving the parabolic stability equations.