ICOSAHOM 2020
International Conference on Spectral and High Order Methods
12th  16th July 2021  Vienna, Austria
Conference Agenda
Overview and details of the sessions of this conference. Please select a date or location to show only sessions at that day or location. Please select a single session for detailed view (with abstracts and downloads if available).
Please note that all times are shown in the time zone of the conference. The current conference time is: 1st Dec 2022, 08:14:24am CET

Session Overview 
Session  
Z5: adaptivity, meshing
 
Presentations  
4:00pm  4:20pm
Comparison of implicit and explicit pmultigrid strategies for compressible NavierStokes discontinuous Galerkin solvers ^{1}ETSIAEUPM (School of Aeronautics  Universidad Politécnica de Madrid); ^{2}Center for Computational Simulation  Universidad Politécnica de Madrid In the last decades, multigrid methods have become an important component of the convergence acceleration techniques for viscous compressible flows. In the context of the highorder methods, pmultigrid is a natural choice since the coarse levels are obtained just by reducing the approximation order. We analyze the computational cost and memory footprint of several popular relaxation strategies for multigrid: explicit RungeKutta methods with local time stepping, GMRES and different elementblock implicit smoothers (BlockJacobi, GaussSeidel) on representative test cases including twodimensional and threedimensional flow simulations. Moreover, we compare the performance of the smoothers in the context of standard DG discretization, where the Gauss nodal points are selected and the entropy stable DG scheme with GaussLobatto points employed. The former is of a main interest, as it has enhanced stability and is convenient for underresolved simulations. 4:20pm  4:40pm
A Direct Method of Generating Highorder Tetrahedral Meshes Using an Advancing Front Approach University of Kansas, United States of America The use of highorder methods has attracted the interest of the scientific computing community, thanks to the ability to deliver highly accurate results at a low computational cost while solving partial differential equations (PDEs). However, while working with curved boundaries, meshes used with highorder PDE solvers need to be highorder so as to accurately capture the curvature of the geometric domain. A highorder mesh is composed of both straightsided and curved elements depending on the curvature of the geometric domain. To date, there are a limited number of algorithms that can generate highorder meshes in a robust manner, i.e., without tangling the mesh. In this talk, we propose a novel method for generating highorder curvilinear tetrahedral meshes using an advancing front approach [1]. There are two categories of highorder mesh generation methods: a posteriori methods and direct methods. A posteriori methods are most commonly used; such methods deform an enriched linear mesh of the geometry into a highorder, curvilinear mesh. Since the mesh curving process can create tangled mesh elements, such methods often require a postprocessing step, such as mesh untangling. Our method uses a direct approach to generate meshes on geometries with curved boundaries. It can generate meshes based on various types of boundary representations, e.g., from computeraided design (CAD) files and patientspecific boundary meshes obtained from medical images. Since our method uses a direct approach instead of an a posteriori approach, it can generate highquality meshes of models obtained from medical images where no CAD representation is available. Our advancing front method generates high quality tetrahedral mesh elements on each iteration, thus omitting the need for postprocessing steps. We present several numerical examples of secondorder tetrahedral meshes generated using our method based on input triangular surface meshes. Our results show that our highorder tetrahedral meshes are of highquality according to the scaled Jacobian and equiangular skewness metrics. Examples are drawn from mechanical engineering and biomedical engineering. We will conclude the talk by discussing the possibility of extending the current algorithm to generate anisotropic meshes in 3D. References [1] F. Mohammadi and S.M. Shontz, A direct method of generating highorder tetrahedral meshes using an advancing front approach, Proceedings of the 29th International Meshing Roundtable, Under review, March 2021. 4:40pm  5:00pm
A secondorder optimizer for parallel distributed mesh curving on virtual geometry Barcelona Supercomputing Center  BSC, Spain We present a mesh curving method for virtual geometry ready for parallel distributed computing. Our technique generates geometrically accurate largescale meshes composed of highquality elements on complex geometry. To this end, we advocate for a penaltybased secondorder optimizer that uses global tight tolerances to converge the distortion residuals. To curve larger meshes without further increasing the computational resources, we aim to reduce the memory footprint, waiting time, and energy consumption. To this end, we consider four main ingredients. First, we propose a degree continuation technique that reduces the number of linear iterations performed on higher polynomial degrees. Second, we adaptively update the penalty parameter to reduce the number of nonlinear problems to solve. Third, we introduce a matrixfree GMRES solver preconditioned with a blockbased successive overrelaxation method. This preconditioner reduces the memory footprint three times. Fourth, we consider an indicator to set a dynamic tolerance for the linear solver and thus, reduce the number of iterations. Compared with our previous solver, we reduce the memory footprint, waiting time, and energy consumption. Using tight tolerances, we curve meshes composed of millions of quartic elements with high stretching for complex virtual geometries using thousands of processors. These capabilities are critical for highfidelity simulations on complex domains using unstructured highorder methods. 5:00pm  5:20pm
Sparse spectral methods for curvilinear domains in the Dedalus code ^{1}Massachusetts Institute of Technology; ^{2}University of Sydney; ^{3}Northwestern University; ^{4}Bates College; ^{5}University of Colorado Boulder Dedalus is an opensource Python framework for solving general partial differential equations at scale with modern spectral methods. Here we will describe recent additions to the code which support general tensor fields in curvilinear domains. This capability is based on recently developed polynomial bases that possess sparse operators for arbitrary tensor calculus in full disks and balls. We will detail our interface for the coordinatefree entry of systems of PDEs including prognostic equations and algebraic constraints. New fast direct solvers allow the inclusion of general boundary conditions, including nonlocal conditions, via generalized tau corrections that can also be specified symbolically. We will also describe ongoing extensions of this framework to support highorder spectral elements and direct couplings with other PDE and BIE solvers. 5:20pm  5:40pm
A Padaptive discontinuous Galerkin incompressible NavierStokes/CahnHilliard solver for multiphase flows ^{1}ETSIAEUPM  School of Aeronautics, Universidad Politécnica de Madrid, Spain; ^{2}Center for Computational Simulation, Universidad Politécnica de Madrid, Spain In this work we present an extension of the works presented in [1] and [2] for the simulation of two phase immiscible fluid flows on p–non–conforming meshes. The solver consists of an incompressible NavierStokes/CahnHilliard system of equations which uses the artificial compressibility method to enforce a divergence–free velocity field. The system is approximated with a high–order discontinuous Galerkin spectral element method. The interelement coupling with different polynomial orders is treated through the standard mortar method [3] in such a way that the new scheme mimics the entropy analysis of the continuous system. Therefore, the entropy remains bounded even for p–non–conforming meshes. A novelty of this work is that the p–adaptation method is applied as a means of reducing the computational cost of a multiphase simulation and we suggest it as an alternative or compliment to h–adaptation, which has been thoroughly studied. Since multiphase flows are typically of unsteady nature, we introduce a methodology to perform dynamic p–adaptation based on the location of the interface between the two phases. The aim is to retain the accuracy of the uniform solver while reducing the associated cost. The validation and assessment of the proposed scheme is done though numerical experiments such as that of a rising bubble or a dam–break, which show that the computational cost can be reduced by 40% to 50% compared to a uniform solver for the same level of accuracy. References [1] Juan Manzanero, Gonzalo Rubio, David A Kopriva, Esteban Ferrer, and Eusebio Valero. Entropy–stable discontinuous Galerkin approximation with summation–by–parts property for the incompressible Navier–Stokes/Cahn–hilliard system. Journal of Computational Physics, 408:109363, 2020. [2] Gerasimos Ntoukas, Juan Manzanero, Gonzalo Rubio, Eusebio Valero, and Esteban Ferrer. A freeenergy stable padaptive nodal discontinuous Galerkin for the CahnHilliard equation. Submitted toJournal of Computational Physics, 2020. [3] David A Kopriva, Stephen L Woodruff, and M Yousuff Hussaini. Computation of electromagnetic scattering with a nonconforming discontinuous spectral element method. International journal for numerical methods in engineering, 53(1):105–122, 2002 5:40pm  6:00pm
A functional oriented truncation error adaptation method UPM, Spain Mesh adaptation allows to localize the numerical degrees of freedom in mesh regions that require increased accuracy. This is particularly important in Computational Fluid Dynamics (CFD) where millions of degrees of freedom are required. Adaptation can be performed by moving the nodes (radaptation) subdividing or merging elements (hadaptation) or by enriching or reducing the polynomial order (padaptation). Mesh adaptation algorithms require a sensor to identify the flow regions that need refinement or coarsening. Three classic approaches can be found in the literature: featurebased, local errorbased and adjointbased [1]. Featurebased adaptation refines where high solution gradients are found. This lowcost approach does not have an explicit connection between the sensor and the error in the solution, and thus often results inefficiently adapted meshes. Adjointbased adaptations define a functional target (e.g., energy, drag or lift) to generate an optimal mesh to reduce the error associated to the selected functional. This approach is efficient in reducing the selected functional error but cannot ensure the reduction of other functionals. Additionally, adjointbased adaption has a high cost since an adjoint problem needs to be solved. Local errorbased include truncation error adaptations, which aims at equidistributing the local error in the mesh domain. For hyperbolic problems, truncation errorbased sensors localize the elements where the discretization error is generated such that the degrees of freedom can be efficiently located at the source, and not wasted in elements polluted by the advected error. The computation of this sensor can be done using the tauestimation approach, which does not require the solution of additional systems and is cheap (see [24] for the application of the tauestimation approach to discontinuous Galerkin (DG) and spectral element method). The advantages of truncation errors have encouraged the development of adaptation algorithms. In particular, padaptation methods have been developed for DG techniques [57]. However, truncation error adaptation techniques have a substantial drawback: they are not related to functional errors. Therefore, the user needs to find arbitrary truncation thresholds to improve the mesh. We propose here to solve this inconvenience by linking the truncation error and the functional error, which enables a cheap truncation error adaptation with the advantages of adjoints, i.e. targeting a specific functional error. [1] C. Roy, Strategies for driving mesh adaptation in CFD (invited), in: 47th AIAA aerospace sciences meeting including the new horizonsforum and aerospace exposition, 2009, p. 1302 [2] G. Rubio, F. Fraysse, J. de Vicente, E. Valero, The estimation of truncation error by tauestimation for Chebyshev spectral collocation method, Journal of Scientific Computing 57 (2013) 146–173 [3] G. Rubio, F. Fraysse, D. A. Kopriva, E. Valero, Quasia priori truncation error estimation in the DGSEM, Journal of Scientific Computing (2014) 1–31 [4] A. M. RuedaRamírez, G. Rubio, E. Ferrer, E. Valero, Truncation error estimation in the panisotropic discontinuous Galerkin spectral element method, Journal of Scientific Computing 78 (2019) 433–466 [5] M. Kompenhans, G. Rubio, E. Ferrer, and E. Valero, Adaptation strategies for high order discontinuous Galerkin methods based on tauestimation, Journal of Computational Physics 306 (2016) 216 – 236 [6] M. Kompenhans, G. Rubio, E. Ferrer, E. Valero, Comparisons of padaptation strategies based on truncation and discretisationerrors for high order discontinuous Galerkin methods, Computers and Fluids 139 (2016) 36–46 [7] A. M. RuedaRamírez, J. Manzanero, E. Ferrer, G. Rubio, E. Valero, A pmultigrid strategy with anisotropic padaptation based on truncation errors for highorder discontinuous Galerkin methods, Journal of Computational Physics 378 (2019) 209–233 
Contact and Legal Notice · Contact Address: Privacy Statement · Conference: ICOSAHOM2020 
Conference Software  ConfTool Pro 2.6.145+CC © 2001–2022 by Dr. H. Weinreich, Hamburg, Germany 