TY - JOUR A1 - Schenkl, Sebastian A1 - Muggenthaler, Holger A1 - Hubig, Michael A1 - Erdmann, Bodo A1 - Weiser, Martin A1 - Zachow, Stefan A1 - Heinrich, Andreas A1 - Güttler, Felix Victor A1 - Teichgräber, Ulf A1 - Mall, Gita T1 - Automatic CT-based finite element model generation for temperature-based death time estimation: feasibility study and sensitivity analysis JF - International Journal of Legal Medicine N2 - Temperature based death time estimation is based either on simple phenomenological models of corpse cooling or on detailed physical heat transfer models. The latter are much more complex, but allow a higher accuracy of death time estimation as in principle all relevant cooling mechanisms can be taken into account. Here, a complete work flow for finite element based cooling simulation models is presented. The following steps are demonstrated on CT-phantoms: • CT-scan • Segmentation of the CT images for thermodynamically relevant features of individual geometries • Conversion of the segmentation result into a Finite Element (FE) simulation model • Computation of the model cooling curve • Calculation of the cooling time For the first time in FE-based cooling time estimation the steps from the CT image over segmentation to FE model generation are semi-automatically performed. The cooling time calculation results are compared to cooling measurements performed on the phantoms under controlled conditions. In this context, the method is validated using different CTphantoms. Some of the CT phantoms thermodynamic material parameters had to be experimentally determined via independent experiments. Moreover the impact of geometry and material parameter uncertainties on the estimated cooling time is investigated by a sensitivity analysis. KW - temperature based death time estimation KW - finite element method KW - CT segmentation KW - sensitivity analysis Y1 - 2017 U6 - https://doi.org/doi:10.1007/s00414-016-1523-0 VL - 131 IS - 3 SP - 699 EP - 712 ER - TY - GEN A1 - Hammerschmidt, Martin A1 - Barth, Carlo A1 - Pomplun, Jan A1 - Burger, Sven A1 - Becker, Christiane A1 - Schmidt, Frank T1 - Reconstruction of photonic crystal geometries using a reduced basis method for nonlinear outputs N2 - Maxwell solvers based on the hp-adaptive finite element method allow for accurate geometrical modeling and high numerical accuracy. These features are indispensable for the optimization of optical properties or reconstruction of parameters through inverse processes. High computational complexity prohibits the evaluation of the solution for many parameters. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem allowing to compute solutions for a parameter configuration orders of magnitude faster. The RBM allows to evaluate linear and nonlinear outputs of interest like Fourier transform or the enhancement of the electromagnetic field in milliseconds. We apply the RBM to compute light-scattering off two dimensional photonic crystal structures made of silicon and reconstruct geometrical parameters. T3 - ZIB-Report - 16-06 KW - finite element method KW - rigorous optical modeling KW - photonic crystals KW - reduced basis method KW - parameter estimation KW - optical metrology Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57249 SN - 1438-0064 ER - TY - GEN A1 - Hammerschmidt, Martin A1 - Herrmann, Sven A1 - Pomplun, Jan A1 - Burger, Sven A1 - Schmidt, Frank T1 - Model order reduction for the time-harmonic Maxwell equation applied to complex nanostructures T2 - Proc. SPIE N2 - Fields such as optical metrology and computational lithography require fast and efficient methods for solving the time-harmonic Maxwell’s equation. Highly accurate geometrical modeling and numerical accuracy atcomputational costs are a prerequisite for any simulation study of complex nano-structured photonic devices. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem based on the hp-adaptive finite element solver JCMsuite capable of handling geometric and non-geometric parameter dependencies allowing for online evaluations in milliseconds. We apply the RBM to compute light-scatteringoptical wavelengths off periodic arrays of fin field-effect transistors (FinFETs) where geometrical properties such as the width and height of the fin and gate can vary in a large range. T3 - ZIB-Report - 16-05 KW - reduced basis method KW - finite element method KW - rigorous optical modeling KW - reduced order models KW - electromagnetic field solver Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-57239 SN - 1438-0064 VL - 9742 SP - 97420M ER - TY - GEN A1 - Hammerschmidt, Martin A1 - Herrmann, Sven A1 - Pomplun, Jan A1 - Zschiedrich, Lin A1 - Burger, Sven A1 - Schmidt, Frank T1 - Reduced basis method for Maxwell's equations with resonance phenomena T2 - Proc. SPIE N2 - Rigorous optical simulations of 3-dimensional nano-photonic structures are an important tool in the analysis and optimization of scattering properties of nano-photonic devices or parameter reconstruction. To construct geometrically accurate models of complex structured nano-photonic devices the finite element method (FEM) is ideally suited due to its flexibility in the geometrical modeling and superior convergence properties. Reduced order models such as the reduced basis method (RBM) allow to construct self-adaptive, error-controlled, very low dimensional approximations for input-output relationships which can be evaluated orders of magnitude faster than the full model. This is advantageous in applications requiring the solution of Maxwell's equations for multiple parameters or a single parameter but in real time. We present a reduced basis method for 3D Maxwell's equations based on the finite element method which allows variations of geometric as well as material and frequency parameters. We demonstrate accuracy and efficiency of the method for a light scattering problem exhibiting a resonance in the electric field. T3 - ZIB-Report - 15-37 KW - reduced basis method KW - finite element method KW - maxwell equation KW - photonic crystal KW - nano-photonics Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-55687 SN - 1438-0064 VL - 9630 SP - 96300R ER - TY - GEN A1 - Hammerschmidt, Martin A1 - Lockau, Daniel A1 - Zschiedrich, Lin A1 - Schmidt, Frank T1 - Optical modelling of incoherent substrate light-trapping in silicon thin film multi-junction solar cells with finite elements and domain decomposition N2 - In many experimentally realized applications, e.g. photonic crystals, solar cells and light-emitting diodes, nano-photonic systems are coupled to a thick substrate layer, which in certain cases has to be included as a part of the optical system. The finite element method (FEM) yields rigorous, high accuracy solutions of full 3D vectorial Maxwell's equations [1] and allows for great flexibility and accuracy in the geometrical modelling. Time-harmonic FEM solvers have been combined with Fourier methods in domain decomposition algorithms to compute coherent solutions of these coupled system. [2,3] The basic idea of a domain decomposition approach lies in a decomposition of the domain into smaller subdomains, separate calculations of the solutions and coupling of these solutions on adjacent subdomains. In experiments light sources are often not perfectly monochromatic and hence a comparision to simulation results might only be justified if the simulation results, which include interference patterns in the substrate, are spectrally averaged. In this contribution we present a scattering matrix domain decomposition algorithm for Maxwell's equations based on FEM. We study its convergence and advantages in the context of optical simulations of silicon thin film multi-junction solar cells. This allows for substrate light-trapping to be included in optical simulations and leads to a more realistic estimation of light path enhancement factors in thin-film devices near the band edge. T3 - ZIB-Report - 14-21 KW - finite element method KW - rigorous optical modeling KW - domain decomposition KW - multi-junction solar cells KW - thin-film silicon solar cells KW - incoherent layers KW - incoherent light-trapping Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-50410 SN - 1438-0064 ER - TY - GEN A1 - Hammerschmidt, Martin A1 - Pomplun, Jan A1 - Burger, Sven A1 - Schmidt, Frank T1 - Adaptive sampling strategies for efficient parameter scans in nano-photonic device simulations N2 - Rigorous optical simulations are an important tool in optimizing scattering properties of nano-photonic devices and are used, for example, in solar cell optimization. The finite element method (FEM) yields rigorous, time-harmonic, high accuracy solutions of the full 3D vectorial Maxwell's equations [1] and furthermore allows for great flexibility and accuracy in the geometrical modeling of these often complex shaped 3D nano-structures. A major drawback of frequency domain methods is the limitation of single frequency evaluations. For example the accurate computation of the short circuit current density of an amorphous silicon / micro-crystalline multi-junction thin film solar cell may require the solution of Maxwell's equations for over a hundred different wavelengths if an equidistant sampling strategy is employed. Also in optical metrology, wavelength scans are frequently used to reconstruct unknown geometrical and material properties of optical systems numerically from measured scatterometric data. In our contribution we present several adaptive numerical integration and sampling routines and study their efficiency in the context of the determination of generation rate profiles of solar cells. We show that these strategies lead to a reduction in the computational effort without loss of accuracy. We discuss the employment of tangential information in a Hermite interpolation scheme to achieve similar accuracy on coarser grids. We explore the usability of these strategies for scatterometry and solar cell simulations. T3 - ZIB-Report - 14-20 KW - finite element method KW - optical simulations KW - adaptive sampling KW - optical metrology KW - parameter scans KW - solar cells Y1 - 2014 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-50395 SN - 1438-0064 ER - TY - GEN A1 - Burger, Sven A1 - Klose, Roland A1 - Schädle, Achim A1 - Zschiedrich, Lin T1 - HelmPole - A finite element solver for scattering problems on unbounded domains: Implementation based on PML N2 - The solution of scattering problems described by the Helmholtz equation on unbounded domains is of importance for a wide variety of applications, for example in electromagnetics and acoustics. An implementation of a solver for scattering problems based on the programming language Matlab is introduced. The solver relies on the finite-element-method and on the perfectly-matched-layer-method, which allows for the simulation of scattering problems on complex geometries surrounded by inhomogeneous exterior domains. This report gives a number of detailed examples and can be understood as a user manual to the freely accessible code of the solver HelmPole. T3 - ZIB-Report - 03-38 KW - Helmholtz equation KW - finite element method KW - perfectly matched layer Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7609 ER - TY - GEN A1 - Kober, Cornelia A1 - Erdmann, Bodo A1 - Sader, Robert A1 - Zeilhofer, Hans-Florian T1 - Simulation of the Human Mandible: Comparison of Bone Mineral Density and Stress/Strain Profiles due to Masticatory Muscles' Traction N2 - The correlation of the inner architecture of bone and its functional loading was already stated by Wolff in 1892. Our objective is to demonstrate this interdependence in the case of the human mandible. For this purpose, stress/strain profiles occuring at a human lateral bite were simulated. Additionally, by a combination of computer graphics modules, a three--dimensional volumetric visualization of bone mineral density could be given. Qualitative correspondences between the density profile of the jaw and the simulated stress/strain profiles could be pointed out. In the long run, this might enable the use of the simulation for diagnosis and prognosis. The solution of the underlying partial differential equations describing linear elastic material behaviour was provided by an adaptive finite element method. Estimates of the discretization errors, local grid refinement, and multilevel techniques guaranteed the reliability and efficiency of the method. T3 - ZIB-Report - 03-23 KW - linear elasticity KW - finite element method KW - adaptive grid refinement KW - stress/strain analysis KW - human mandible Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7458 ER - TY - GEN A1 - Weiser, Martin A1 - Schiela, Anton A1 - Deuflhard, Peter T1 - Asymptotic Mesh Independence of Newton's Method Revisited N2 - The paper presents a new affine invariant theory on asymptotic mesh independence of Newton's method in nonlinear PDEs. Compared to earlier attempts, the new approach is both much simpler and more natural from the algorithmic point of view. The theory is exemplified at collocation methods for ODE boundary value problems and at finite element methods for elliptic PDE problems. T3 - ZIB-Report - 03-13 KW - mesh independence KW - nonlinear partial differential equations KW - Newton method KW - finite element method KW - collocation method Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7352 ER - TY - GEN A1 - Erdmann, Bodo A1 - Kober, Cornelia A1 - Lang, Jens A1 - Sader, Robert A1 - Zeilhofer, Hans-Florian A1 - Deuflhard, Peter T1 - Efficient and Reliable Finite Element Methods for Simulation of the Human Mandible N2 - By computed tomography data (CT), the individual geometry of the mandible is quite well reproduced, also the separation between cortical and trabecular bone. Using anatomical knowledge about the architecture and the functional potential of the masticatory muscles, realistic situations were approximated. The solution of the underlying partial differential equations describing linear elastic material behaviour is provided by an adaptive finite element method. Estimations of the discretization error, local grid refinement, and multilevel techniques guarantee the reliability and efficiency of the method. T3 - ZIB-Report - 01-14 KW - mandible KW - sensitivity analysis KW - finite element method KW - adaptive grid refinement Y1 - 2001 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-6403 ER -