65N30 Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Refine
Year of publication
Document Type
- ZIB-Report (27)
- Doctoral Thesis (1)
- Habilitation (1)
Language
- English (29)
Is part of the Bibliography
- no (29)
Keywords
- finite elements (6)
- Maxwell's equations (4)
- preconditioned conjugate gradient methods (3)
- Helmholtz equation (2)
- adaptive mesh refinement (2)
- discontinuous solution (2)
- edge elements (2)
- embedded interface (2)
- finite element method (2)
- finite volume methods (2)
- generic programming (2)
- hierarchical basis (2)
- optimal control (2)
- partial differential equations (2)
- perfectly matched layer (2)
- transparent boundary conditions (2)
- Adaptive Finite Elements (1)
- Algebraic multigrid (1)
- Currents (1)
- Diffeomorphic Registration (1)
- Domain Decomposition (1)
- Elliptic equations (1)
- Fini (1)
- Finite Elements (1)
- Finite Volumen (1)
- Large Deformation (1)
- Matching (1)
- Multigrid Methods (1)
- N\'{e}d\'{e}lec spaces (1)
- Newton method (1)
- PML (1)
- Pelvic heating (1)
- Poisson equation (1)
- Pole condition (1)
- Rayleigh-Ritz and Galerkin Methods (1)
- Schroedinger equation (1)
- Transparent boundary conditions (1)
- Unbounded domains (1)
- Wide angle one-way equation (1)
- Wiener-Hopf techniques (1)
- adaptive (1)
- adaptive finite elements (1)
- adaptive grid orientation (1)
- adaptive grid refinement (1)
- adaptive methods (1)
- affine conjugate Newton methods (1)
- ariable and discontinuous coefficients (1)
- collocation method (1)
- convection diffusion equation (1)
- countable (1)
- cranio-maxillofacial surgery (1)
- data compression (1)
- differential equations: (1)
- directed refinement (1)
- discontinuous finite element (1)
- diskontinuierliche Lösung (1)
- eigenvalue problem (1)
- eingebettete Ränder (1)
- electromagnetic scattering (1)
- elliptische Gleichungen (1)
- finite element methods: (1)
- floating point compression (1)
- grid generation (1)
- hamiltonian (1)
- hyperthermia (1)
- implant design (1)
- incomplete factorization (1)
- internal and boundary layers (1)
- lossy compression (1)
- medical t (1)
- mesh coarsening (1)
- mesh independence (1)
- multigrid (1)
- multigrid methods (1)
- multilev (1)
- multilevel (1)
- multilevel methods (1)
- neutron transport equation (1)
- non-symmetric linear systems (1)
- nonconvex minimization (1)
- nonlinear elastomechnics (1)
- nonlinear heat transfer (1)
- nonlinear partial differential equations (1)
- nonuniform meshes (1)
- optical chip design (1)
- optimal convergence rates (1)
- ordinary (1)
- partial (1)
- partial differential equation (1)
- perfectly matched layers (1)
- phased array optimization (1)
- pole condition (1)
- polyconvex elasticity (1)
- soft tissue simulation (1)
- sparse secant method (1)
- symmetric Gauss-Seidel (1)
- variable and discontinuous coefficients (1)
- variable und diskontinuierliche Koeffizienten (1)
- waveguide (1)
Institute
Subspace decompositions of finite element spaces based on $L2$-like orthogonal projections play an important role for the construction and analysis of multigrid like iterative methods. Recently several authors proved the equivalence of the associated discrete norms with the $H^1$-norm. The present report gives an elementary, self-contained derivation of this result which is based on the use of $ K$-functionals known from the theory of interpolation spaces. {\bf Keywords:} multilevel methods, nonuniform meshes, optimal convergence rates. {\bf AMS(MOS) Subject classifications:} 65N55, 65N30, 65N50.
We present a finite volume method for the solution of the two-dimensional Poisson equation $ \nabla\cdot( \beta( {\mbox{\boldmath $x$}}) \nabla u({\mbox{\boldmath $x$}})) = f(\mbox{\boldmath $x$}) $ with variable, discontinuous coefficients and solution discontinuities on irregular domains. The method uses bilinear ansatz functions on Cartesian grids for the solution $u({\mbox{\boldmath $x$})$ resulting in a compact nine-point stencil. The resulting linear problem has been solved with a standard multigrid solver. Singularities associated with vanishing partial volumes of intersected grid cells or the dual bilinear ansatz itself are removed by a two-step asymptotic approach. The method achieves second order of accuracy in the $L^\infty$ and $L^2$ norm.
The work presents a new approach to the numerical solution of time-harmonic and time-dependent scattering problems. We replace Sommerfeld's radiation condition valid for the Helmholtz equation by a more general concept called pole condition. The pole condition is based on the Laplace transform of the exterior solution and allows a characterization of outgoing waves. Both new insight into the analysis of scattering problems as well as new numerical algorithms are obtained.
In this paper we propose a new finite element realization of the Perfectly Matched Layer method (PML-method). Our approach allows to deal with arbitrary shaped polygonal domains and with certain types of inhomogeneous exterior domains. Among the covered inhomogeneities are open waveguide structures playing an essential role in integrated optics. We give a detailed insight to implementation aspects. Numerical examples show exponential convergence behavior to the exact solution with the thickness of the PML sponge layer.
We present a second order sharp interface finite volume method for the solution of the three-dimensional poisson equation with variable coefficients on Cartesian grids. In particular, we focus on interface problems with discontinuities in the coefficient, the source term, the solution, and the fluxes across the interface. The method uses standard piecewiese trilinear finite elements for normal cells and a double piecewise trilinear ansatz for the solution on cells intersected by the interface resulting always in a compact 27-point stencil. Singularities associated with vanishing partial volumes of intersected grid cells are removed by a two-term asymptotic approach. In contrast to the 2D method presented by two of the authors in [M.~Oevermann, R.~Klein: A Cartesian grid finite volume method for elliptic equations with variable coefficients and embedded interfaces, J.~Comp.~Phys.~219 (2006)] we use a minimization technique to determine the unknown coefficients of the double trilinear ansatz. This simplifies the treatment of the different cut-cell types and avoids additional special operations for degenerated interface topologies. The resulting set of linear equations has been solved with a BiCGSTAB solver preconditioned with an algebraic multigrid. In various testcases -- including large coefficient ratios and non-smooth interfaces -- the method achieves second order of accuracy in the L_inf and L_2 norm.
In this paper it is shown that for highly nonuniformly refined triangulations the condition number of the BPX preconditioner for elliptic finite element problems grows at most linearly in the depth of refinement. This is achieved by viewing the computational available version of the BPX preconditioner as an abstract additive Schwarz method with exact solvers. {\bf AMS CLASSIFICATION:} 65F10, 65F35, 65N20, 65N30.
A cascadic multigrid (CMG) method for elliptic problems with strong material jumps is proposed and analyzed. Non--matching grids at interfaces between subdomains are allowed and treated by mortar elements. The arising saddle point problems are solved by a subspace confined conjugate gradient method as smoother for the CMG. Details of algorithmic realization including adaptivity are elaborated. Numerical results illustrate the efficiency of this CMG algorithm.
Adaptive Multigrid Methods for the Vectorial Maxwell Eigenvalue Problem for Optical Waveguide Design
(2000)
This paper has been motivated by the need for a fast robust adaptive multigrid method to solve the vectorial Maxwell eigenvalue problem arising from the design of optical chips. Our nonlinear multigrid methods are based on a previous method for the scalar Helmholtz equation, which must be modified to cope with the null space of the Maxwell operator due to the divergence condition. We present two different approaches. First, we present a multigrid algorithm based on an edge element discretization of time-harmonic Maxwell's equations, including the divergence condition. Second, an explicit elimination of longitudinal magnetic components leads to a nodal discretization known to avoid discrete \emph{spurious modes} also and a vectorial eigenvalue problem, for which we present a multigrid solver. Numerical examples show that the edge element discretization clearly outperforms the nodal element approach.
The paper surveys recent progress in a joint mathematical-medical project on cancer therapy planning. Within so-called regional hyperthermia the computational task is to tune a set of coupled radiofrequency antennas such that a carefully measured tumor is locally heated, but any outside hot spots are avoided. A mathematical model of the whole clinical system -- air, applicator with antennas, water bolus, individual patient body -- involves Maxwell's equations in inhomogeneous media and a parabolic bioheat transfer equation, which represents a simplified model of heat transfer in the human body (ignoring strong blood vessel heat transport). Both PDEs need to be computed fast and to medical reliability (!) on a workstation within a clinical environment. This requirement triggered a series of new algorithmic developments to be reported here, among which is an adaptive multilevel FEM for Maxwell's equations, which dominates the numerical simulation time. In total, however, the main bulk of computation time (see Table 3 in Section 4 below) still goes into segmentation -- a necessary preprocessing step in the construction a 3D virtual patient from the input of a stack of 2D computed tomograms (left out here).
The paper extends affine conjugate Newton methods from convex to nonconvex minimization, with particular emphasis on PDE problems originating from compressible hyperelasticity. Based on well-known schemes from finite dimensional nonlinear optimization, three different algorithmic variants are worked out in a function space setting, which permits an adaptive multilevel finite element implementation. These algorithms are tested on two well-known 3D test problems and a real-life example from surgical operation planning.