{"id":{"repo_id":"uiuc","oai_identifier":"oai:www.ideals.illinois.edu:2142/124325"},"canonical_url":"https://search.dev.ndltd.org/etd/uiuc/oai:www.ideals.illinois.edu:2142/124325","repository":{"repo_id":"uiuc","name":"University of Illinois - Urbana-Champaign","base_url":"https://www.ideals.illinois.edu/oai-pmh"},"display":{"title":"Dispersion analysis, time parallelization, and GPU autotuning for finite element methods","abstract":"Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2024-09-16 without embargo terms","abstract_html":"Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2024-09-16 without embargo terms","abstract_has_math":false,"creators":["Christensen, Nicholas J"],"institution":"University of Illinois at Urbana-Champaign","degree_name":"Ph.D.","degree_level":"Dissertation","degree_discipline":"Computer Science","degree_department":null,"school":null,"contributors":["Fischer, Paul F","Klöckner, Andreas","Olson, Luke N","Parish, Eric J"],"advisors":[],"committee_chairs":[],"committee_members":[],"year":2024,"date_issued":"2024-05","date_published":"2024-05","updated_at":"2026-07-22T22:25:00Z","subjects":["Finite Element","Spectral Element","Dispersion","Pde","Autotuning","Parallel-in-time","Euler-lagrange","Least-squares","Gpu","Eigenvalue Avoidance","Einstein Summation"],"languages":["en","eng"],"rights":["Copyright 2024 Nicholas Christensen. CC0: This work has been marked as dedicated to the public domain."],"rights_urls":[],"identifier_entries":[]},"links":{"outbound_url":"https://hdl.handle.net/2142/124325","outbound_label":"Handle","outbound_source":"dc:identifier"},"metadata_groups":[{"id":"people","label":"People","entries":[{"key":"dc:contributor","label":"Contributor","values":["Fischer, Paul F","Klöckner, Andreas","Olson, Luke N","Parish, Eric J"]},{"key":"dc:creator","label":"Author","values":["Christensen, Nicholas J"]}]},{"id":"academic_context","label":"Academic Context","entries":[{"key":"dc:date","label":"Dc Date","values":["2024-05","2024-04-23"]},{"key":"dc:type","label":"Dc Type","values":["text"]},{"key":"thesis:degree_discipline","label":"Discipline","values":["Computer Science"]},{"key":"thesis:degree_level","label":"Degree Level","values":["Dissertation"]},{"key":"thesis:degree_name","label":"Degree Name","values":["Ph.D."]},{"key":"thesis:institution_name","label":"Thesis Institution Name","values":["University of Illinois at Urbana-Champaign"]}]},{"id":"subjects_keywords","label":"Subjects and Keywords","entries":[{"key":"dc:subject","label":"Dc Subject","values":["Finite Element","Spectral Element","Dispersion","Pde","Autotuning","Parallel-in-time","Euler-lagrange","Least-squares","Gpu","Eigenvalue Avoidance","Einstein Summation"]}]},{"id":"language_rights","label":"Language and Rights","entries":[{"key":"dc:language","label":"Dc Language","values":["en","eng"]},{"key":"dc:rights","label":"Dc Rights","values":["Copyright 2024 Nicholas Christensen. CC0: This work has been marked as dedicated to the public domain."]}]},{"id":"identifiers","label":"Identifiers","entries":[{"key":"dc:identifier","label":"Identifier","values":["https://hdl.handle.net/2142/124325"]}]},{"id":"additional","label":"Additional Metadata","entries":[{"key":"dc:description","label":"Description","values":["Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2024-09-16 without embargo terms","The student, Nicholas Christensen, accepted the attached license on 2024-04-22 at 11:01.","The student, Nicholas Christensen, submitted this Dissertation for approval on 2024-04-22 at 11:29.","This Dissertation was approved for publication on 2024-04-23 at 11:33.","DSpace SAF Submission Ingestion Package generated from Vireo submission #20471 on 2024-09-16 at 00:35:21","This dissertation explores strategies to enhance the accuracy and computational performance of finite element methods. Specifically, we analyze the dispersive error of the spectral element method, investigate a novel least-squares parallel-in-time formulation, and propose an autotuning approach for a discontinuous Galerkin finite element solver. We investigate the dispersion properties of the spectral element method (SEM) when applied to advection and advection-diffusion problems on a 1D periodic domain. Our analysis spans both the well-resolved (asymptotic) limit and the marginally resolved (pre-asymptotic) limit. To achieve this, we systematically explore a wide range of parameters, including wave numbers, element counts, and local polynomial orders. We observe that high-order methods demand fewer points-per-wavelength (PPW) than low-order methods to meet engineering tolerances during long time-integration. Remarkably, Gottlieb’s observation that polynomial-based spectral methods require approximately PPW > 5 for engineering tolerances holds true across various polynomial orders and element counts. Comparing use of exact quadrature on the Gauss-Legendre points and inexact quadrature (employing a diagonal mass matrix) on the Gauss-Lobotto-Legendre points, we find inexact quadrature does not significantly compromise solution accuracies at high polynomial orders. At high polynomial orders (N > 4), we observe error spikes near specific values of PPW disrupt the convergence behavior, regardless of the quadrature method used. These spikes arise from previously identified gaps in the eigenvalue spectrum of the discrete operators, leading to unrepresentable phase velocities. We demonstrate that diffusive mechanisms – whether introduced numerically (via time-relaxation or an upwind discontinuous Galerkin formulation) or arising naturally from the physics of the problem – can largely mitigate these error spikes. A two-dimensional model problem further illustrates the effectiveness of the mitigation strategy. We propose MG-HLS-PinT, a novel parallel-in-time method based on multigrid principles. Derived from a normal-equations formulation of a semi-discrete partial differential equation (PDE), this approach shows potential in accelerating the solution of hyperbolic PDEs. However, it currently demands high processor counts and problems with stringent accuracy demands. We explore several methods of speeding up the approach and identify promising avenues for future exploration. In the context of the MIRGE-Com DGFEM simulation library, we introduce sub-batching as a method to enhance the computation of large fused batched Einstein summation (einsum) GPU kernels. By limiting the number of concurrent einsum computations and minimizing contention for local memory and cache, sub-batching significantly improves performance compared to the baseline."]},{"key":"dc:format","label":"Dc Format","values":["application/pdf"]},{"key":"dc:title","label":"Title","values":["Dispersion analysis, time parallelization, and GPU autotuning for finite element methods"]}]}],"canonical_facts":{"dc:contributor":["Fischer, Paul F","Klöckner, Andreas","Olson, Luke N","Parish, Eric J"],"dc:creator":["Christensen, Nicholas J"],"dc:date":["2024-05","2024-04-23"],"dc:description":["Submission original under an indefinite embargo labeled 'Open Access'. The submission was exported from vireo on 2024-09-16 without embargo terms","The student, Nicholas Christensen, accepted the attached license on 2024-04-22 at 11:01.","The student, Nicholas Christensen, submitted this Dissertation for approval on 2024-04-22 at 11:29.","This Dissertation was approved for publication on 2024-04-23 at 11:33.","DSpace SAF Submission Ingestion Package generated from Vireo submission #20471 on 2024-09-16 at 00:35:21","This dissertation explores strategies to enhance the accuracy and computational performance of finite element methods. Specifically, we analyze the dispersive error of the spectral element method, investigate a novel least-squares parallel-in-time formulation, and propose an autotuning approach for a discontinuous Galerkin finite element solver. We investigate the dispersion properties of the spectral element method (SEM) when applied to advection and advection-diffusion problems on a 1D periodic domain. Our analysis spans both the well-resolved (asymptotic) limit and the marginally resolved (pre-asymptotic) limit. To achieve this, we systematically explore a wide range of parameters, including wave numbers, element counts, and local polynomial orders. We observe that high-order methods demand fewer points-per-wavelength (PPW) than low-order methods to meet engineering tolerances during long time-integration. Remarkably, Gottlieb’s observation that polynomial-based spectral methods require approximately PPW > 5 for engineering tolerances holds true across various polynomial orders and element counts. Comparing use of exact quadrature on the Gauss-Legendre points and inexact quadrature (employing a diagonal mass matrix) on the Gauss-Lobotto-Legendre points, we find inexact quadrature does not significantly compromise solution accuracies at high polynomial orders. At high polynomial orders (N > 4), we observe error spikes near specific values of PPW disrupt the convergence behavior, regardless of the quadrature method used. These spikes arise from previously identified gaps in the eigenvalue spectrum of the discrete operators, leading to unrepresentable phase velocities. We demonstrate that diffusive mechanisms – whether introduced numerically (via time-relaxation or an upwind discontinuous Galerkin formulation) or arising naturally from the physics of the problem – can largely mitigate these error spikes. A two-dimensional model problem further illustrates the effectiveness of the mitigation strategy. We propose MG-HLS-PinT, a novel parallel-in-time method based on multigrid principles. Derived from a normal-equations formulation of a semi-discrete partial differential equation (PDE), this approach shows potential in accelerating the solution of hyperbolic PDEs. However, it currently demands high processor counts and problems with stringent accuracy demands. We explore several methods of speeding up the approach and identify promising avenues for future exploration. In the context of the MIRGE-Com DGFEM simulation library, we introduce sub-batching as a method to enhance the computation of large fused batched Einstein summation (einsum) GPU kernels. By limiting the number of concurrent einsum computations and minimizing contention for local memory and cache, sub-batching significantly improves performance compared to the baseline."],"dc:format":["application/pdf"],"dc:identifier":["https://hdl.handle.net/2142/124325"],"dc:language":["en","eng"],"dc:rights":["Copyright 2024 Nicholas Christensen. CC0: This work has been marked as dedicated to the public domain."],"dc:subject":["Finite Element","Spectral Element","Dispersion","Pde","Autotuning","Parallel-in-time","Euler-lagrange","Least-squares","Gpu","Eigenvalue Avoidance","Einstein Summation"],"dc:title":["Dispersion analysis, time parallelization, and GPU autotuning for finite element methods"],"dc:type":["text"],"thesis:degree_discipline":["Computer Science"],"thesis:degree_level":["Dissertation"],"thesis:degree_name":["Ph.D."],"thesis:institution_name":["University of Illinois at Urbana-Champaign"]},"updated_at":"2026-07-22T22:25:00Z"}