{"id":1656,"date":"2015-03-11T13:43:40","date_gmt":"2015-03-11T12:43:40","guid":{"rendered":"http:\/\/www.mi.uni-koeln.de\/NumSim\/news\/"},"modified":"2020-12-16T10:15:00","modified_gmt":"2020-12-16T09:15:00","slug":"news","status":"publish","type":"page","link":"https:\/\/www.mi.uni-koeln.de\/NumSim\/news\/","title":{"rendered":"News"},"content":{"rendered":"<ul class=\"display-posts-listing\"><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2026\/03\/16\/workshop-trixi-user-and-developer-interaction-trudi-2026\/\">Workshop: Trixi User and Developer Interaction (TRUDI) 2026<\/a> <span class=\"date\">March 16, 2026<\/span><div class=\"content\"><p>From 16 to 20 March, we are hosting the TRUDI workshop at the Division of Mathematics at the University of Cologne.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024-258x300.jpg\" alt=\"\" width=\"258\" height=\"300\" class=\"aligncenter size-medium wp-image-6693\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024-258x300.jpg 258w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024-882x1024.jpg 882w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024-768x892.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/03\/TRUDI2026_1024.jpg 1024w\" sizes=\"(max-width: 258px) 100vw, 258px\" \/><\/a><\/p>\n<p>TRUDI 2026 is a workshop bringing together users and developers of the Trixi Framework, which includes the Julia packages <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a>, <a href=\"https:\/\/github.com\/trixi-framework\/TrixiShallowWater.jl\">TrixiShallowWater.jl<\/a>, <a href=\"https:\/\/github.com\/trixi-framework\/TrixiAtmo.jl\">TrixiAtmo.jl<\/a>, and <a href=\"https:\/\/github.com\/trixi-framework\/TrixiParticles.jl\">TrixiParticles.jl<\/a>. The program will feature talks, tutorials, and hackathons aimed at consolidating the existing community, aligning long-term strategy, and advancing key features in high-performance scientific computing.<\/p>\n<p>In addition to the main Trixi Framework packages, we rely on &#8211; and actively collaborate with &#8211; a broader ecosystem of upstream dependencies and partner projects maintained by members of our community, including (but not limited to):<\/p>\n<ul>\n<li><a href=\"https:\/\/github.com\/trixi-framework\/HOHQMesh\">HOHQMesh<\/a> \u00b7 <a href=\"https:\/\/github.com\/trixi-framework\/HOHQMesh.jl\">HOHQMesh.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/trixi-framework\/libtrixi\">libtrixi\u00a0<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/cburstedde\/p4est\">p4est<\/a> \u00b7 <a href=\"https:\/\/github.com\/trixi-framework\/P4est.jl\">P4est.jl\u00a0<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/trixi-framework\/PointNeighbors.jl\">PointNeighbors.jl <\/a><\/li>\n<li><a href=\"https:\/\/github.com\/jlchan\/StartUpDG.jl\">StartUpDG.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/ranocha\/SummationByPartsOperators.jl\">SummationByPartsOperators.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/DLR-AMR\/t8code\">t8code<\/a> \u00b7 <a href=\"https:\/\/github.com\/DLR-AMR\/T8code.jl\">T8code.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/trixi-framework\/Trixi2Vtk.jl\">Trixi2Vtk.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/trixi-framework\/TrixiBottomTopography.jl\">TrixiBottomTopography.jl<\/a><\/li>\n<li><a href=\"https:\/\/github.com\/Arpit-Babbar\/TrixiLW.jl\">TrixiLW.jl<\/a><\/li>\n<\/ul>\n<p>The event is co-hosted by <a href=\"https:\/\/trixi-framework.github.io\/trudi-2026\/snubic.html\">SNuBIC<\/a>, the DFG research unit on \u201cStructure-preserving Numerical Methods for Bulk and Interface-Coupling of Heterogeneous Models\u201d (FOR5409). We gratefully acknowledge further support by the key profile area <a href=\"https:\/\/trixi-framework.github.io\/trudi-2026\/universe.html\">Dynamics of the Universe<\/a> of the University of Cologne.<\/p>\n<p>Further information can be found here: <a href=\"https:\/\/trixi-framework.github.io\/trudi-2026\/\">https:\/\/trixi-framework.github.io\/trudi-2026\/<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2026\/03\/09\/new-preprint-published-sbp-fdec-summation-by-parts-finite-difference-exterior-calculus\/\">New preprint published: SBP-FDEC: Summation-by-Parts Finite Difference Exterior Calculus<\/a> <span class=\"date\">March 9, 2026<\/span><div class=\"content\"><p>We extend the framework of Finite Element Exterior Calculus (FEEC) to Summation-by-Parts (SBP) Finite Difference (FD) operators to derive an energy-conservative and divergence-free scheme for the homogeneous Maxwell\u2019s equations.<\/p>\n<p><a href=\"https:\/\/arxiv.org\/abs\/2511.20529v1\">https:\/\/arxiv.org\/abs\/2511.20529v1<\/a><\/p>\n<p><strong>Abstract<\/strong><br \/>\nWe demonstrate that we can carry over the strategy of Finite Element Exterior Calculus (FEEC) to Summation-by-Parts (SBP) Finite Difference (FD) methods to achieve divergence- and curl-free discretizations. This is not obvious at first sight, as for SBP-FD no basis functions are known, but only values and derivatives at points. The key is a remarkable analytic relationship that enables us to construct compatible operators using integral and nodal degrees of freedom. Pre-existing SBP-FD matrix operators can then be used to obtain nodal values from the integral degrees of freedom to derive a scheme with the desired properties.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2026\/01\/30\/snapshot-outreach-visit-to-kreuzgasse-gymnasium-cologne\/\">Snapshot: Outreach visit to Kreuzgasse Gymnasium Cologne<\/a> <span class=\"date\">January 30, 2026<\/span><div class=\"content\"><p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small-1024x768.jpg\" alt=\"\" width=\"584\" height=\"438\" class=\"aligncenter size-large wp-image-6679\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small-1024x768.jpg 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small-300x225.jpg 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small-768x576.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small-400x300.jpg 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/kreuzgasse_visit_small.jpg 1152w\" sizes=\"(max-width: 584px) 100vw, 584px\" \/><\/a><\/p>\n<p>On Thursday, 15 January 2026, Simone Chiocchetti visited students of the 11. Klasse at Kreuzgasse Gymnasium in Cologne for an interactive meeting on computational physics, engineering, and applied mathematics, with a particular focus on the numerical simulation of fluids.<br \/>\nThe discussion covered a wide range of topics: from the mathematical foundations of partial differential equations and their discretization, to real-world applications in civil engineering, astrophysics, and industrial design. Students learned how researchers model phenomena ranging from river floods and debris flows to turbulence in aerodynamics, and why such simulations matter when experiments are expensive, dangerous, or simply impossible. The conversation also touched on the practical challenges of supercomputing, multiphase flows one can observe in a cup of coffee, and the surprising behavior of steel structures under large deformations.<br \/>\nThe students participated actively throughout, contributing thoughtful questions on subjects including chaos in planetary orbits and its connection to fluid turbulence, the multiscale nature of turbulent flows, and the computational tricks needed to make large simulations feasible.<br \/>\nMany thanks to the students for their engagement and curiosity, and to Markus Klein, their physics teacher, for contributing to organizing the event.<br \/>\nThis activity is funded by the European Union\u2019s Horizon Europe Research and Innovation Programme under the Marie Sk\u0142odowska-Curie Postdoctoral Fellowship MoMeNTUM (grant agreement No. 101109532).<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-300x63.png\" alt=\"\" width=\"300\" height=\"63\" class=\"aligncenter size-medium wp-image-6678\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-300x63.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-1024x215.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-768x161.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-500x105.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image.png 1100w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/12\/05\/talk-alex-bercik-with-volume-dissipation-for-finite-difference-and-spectral-element-summation-by-parts-methods\/\">Talk: Alex Bercik (University of Toronto) with Volume Dissipation for Finite-Difference and Spectral-Element Summation-By-Parts Methods<\/a> <span class=\"date\">December 5, 2025<\/span><div class=\"content\"><p>We construct provably stable, conservative, and accurate volume dissipation within a tensor-product framework. The dissipation operators can be applied to any scheme that uses the SBP framework, including high-order entropy-stable schemes. We clarify the incorporation of a variable coefficient within the operator structure, and connect the presented volume dissipation to the use of upwind SBP operators, both in a finite-difference and spectral-element context. Numerical examples featuring the linear convection, Burgers, and Euler equations verify the properties of the constructed dissipation operators, and assess their performance compared to existing upwind schemes, including local linear stability behaviour. When applied to entropy-stable schemes, the presented approach results in accurate and robust methods that can solve a broader range of problems where comparable existing methods fail.<\/p>\n<p>Link to the paper: <a href=\"https:\/\/arxiv.org\/abs\/2503.12670\">arXiv:2503.12670<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/11\/05\/snapshot-de-mach-2000-astrophysical-jet-with-dg-and-subcell-positivity-limiters-using-amr\/\">Snapshot: DE Mach 2000 astrophysical Jet with DG and subcell positivity limiters using AMR<\/a> <span class=\"date\">November 5, 2025<\/span><div class=\"content\"><p>Numerical simulation of an astrophysical Jet with a Mach number of about 2000. We solve the compressible Euler equations of gas dynamics with a FCT-type subcell limiting approach combining a (forth-order accurate) discontinuous Galerkin (DG) method with a first-order accurate finite volume (FV) method at the node level to impose positivity of density and pressure. Additionally, a new mortar approach also ensures positivity of density and pressure.<\/p>\n<p>At end time T=0.0015, we have 79,708 elements which result in 1,275,328 degrees of freedom. We use the entropy-conserving and kinetic energy preserving flux of Chandrashekar for the volume fluxes and the local Lax-Friedrichs flux for the surface fluxes of the DG and FV methods.<\/p>\n<p>The results were obtained using <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a>.<\/p>\n<p><iframe loading=\"lazy\" title=\"Mach 2000 astrophysical Jet with DG and subcell positivity limiters using AMR\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/RDcyKNc_OXA?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe><\/p>\n<p>Reference:<br \/>\n[1] A. M. Rueda-Ram\u00edrez, B. Bolm, D. Kuzmin, G. J. Gassner: Monolithic Convex Limiting for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods, <a href=\"https:\/\/link.springer.com\/article\/10.1007\/s42967-023-00321-6\">https:\/\/link.springer.com\/article\/10.1007\/s42967-023-00321-6<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/10\/06\/snapshot-modern-transfer-learning-architectures-improve-archaeological-artifact-classification\/\">Snapshot: Modern transfer learning architectures improve archaeological artifact classification<\/a> <span class=\"date\">October 6, 2025<\/span><div class=\"content\"><p>With the increasing complexity of deep neural networks and continual architectural improvements, AI models have achieved remarkable success in image classification, surpassing 90% TOP-1 accuracy on ImageNet. Such high performance highlights their effectiveness across diverse domains and supports reliable transfer learning for smaller, specialized datasets. Artifact classification is one example where knowledge transfer from large-scale datasets proves highly beneficial. In (a), a pretrained CNN from ImageNet provides powerful feature extraction for archaeological image classification. In (b), this transferred knowledge is further extended to multi-image classification tasks through optimized feature fusion, enhancing overall model performance.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure.png\" alt=\"\" width=\"1024\" height=\"778\" class=\"aligncenter size-full wp-image-6609\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure-300x228.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure-768x584.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/10\/Figure-395x300.png 395w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/09\/08\/snapshot-subcell-idp-fct-limiting-simulation-of-the-mach-2000-astrophysical-jet\/\">Snapshot: Subcell IDP\/FCT Limiting Simulation of the Mach 2000 Astrophysical Jet<\/a> <span class=\"date\">September 8, 2025<\/span><div class=\"content\"><p>The videos show a result of simulations with DGSEM and a polynomial degree of 3. At the end of the simulation, the solution has about 150000 degrees of freedom. Subcell IDP\/FCT-type Limiting provides for stability. Moreover, the simulation uses adaptive mesh refinement, combined with a new IDP mortar type that ensures positivity of density and pressure.<br \/>\nOn the left side, local (shock-capturing) limiters are used in the volume integral, while on the right side only limiters that ensure physics admissibility (positivity of density and pressure) are enabled.<br \/>\n<div style=\"width: 584px;\" class=\"wp-video\"><!--[if lt IE 9]><script>document.createElement('video');<\/script><![endif]-->\n<video class=\"wp-video-shortcode\" id=\"video-6599-1\" width=\"584\" height=\"296\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/09\/video.mp4?_=1\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/09\/video.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/09\/video.mp4<\/a><\/video><\/div><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/07\/21\/snapshot-sph-simulation-of-a-freediving-bifin\/\">Snapshot: SPH Simulation of a freediving bifin<\/a> <span class=\"date\">July 21, 2025<\/span><div class=\"content\"><p>Simulation of a freediving bifin, modelling the fluid with \u03b4+-SPH and periodic boundaries and the fin with Total Lagrangian SPH.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6584-2\" width=\"584\" height=\"329\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/07\/out_fin_adami_consistent_40_compressed.mp4?_=2\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/07\/out_fin_adami_consistent_40_compressed.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/07\/out_fin_adami_consistent_40_compressed.mp4<\/a><\/video><\/div>\n<p>The software <a href=\"https:\/\/github.com\/trixi-framework\/TrixiParticles.jl\">TrixiParticles.jl<\/a> was used for this simulation.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/05\/12\/snapshot-investigating-atmospheric-dynamics-the-held-suarez-test-case-with-trixiatmo-jl\/\">Snapshot: Investigating Atmospheric Dynamics: The Held-Suarez Test Case with TrixiAtmo.jl<\/a> <span class=\"date\">May 12, 2025<\/span><div class=\"content\"><p>Our research group is currently utilizing <a href=\"https:\/\/github.com\/trixi-framework\/TrixiAtmo.jl\">TrixiAtmo.jl<\/a> to explore the Held-Suarez test case, a widely recognized benchmark for atmospheric general circulation and climate models. First proposed by Held and Suarez in their 1994 <a href=\"https:\/\/doi.org\/10.1175\/1520-0477(1994)075%3C1825:APFTIO%3E2.0.CO;2\">paper<\/a>, this idealized setup is designed to capture fundamental large-scale flow features of Earth&#8217;s atmosphere.<\/p>\n<p>The model incorporates a simplified atmospheric forcing that establishes a decreasing temperature gradient from the equator to the poles and maintains a vertically balanced state. It also accounts for idealized boundary layer friction affecting the wind field. This configuration drives a characteristic atmospheric circulation: warm air rises at the equator, while colder air flows towards it at lower levels. The Coriolis force deflects this moving air, leading to the development of prominent jet streams in both the Northern and Southern Hemispheres.<\/p>\n<p>For our simulations, TrixiAtmo.jl employs a Discontinuous Galerkin Spectral Element Method (DGSEM) with flux differencing. We discretize the Earth&#8217;s atmosphere using a cubed sphere grid with 6 patches, each consisting of 10 x 10 x 8 cells. Starting from a state with constant temperature and without any motion, the video visualizes the simulated near-surface air temperature and the developing flow patterns.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6574-3\" width=\"584\" height=\"584\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/05\/hs_T_surface.mp4?_=3\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/05\/hs_T_surface.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/05\/hs_T_surface.mp4<\/a><\/video><\/div>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/04\/07\/new-preprint-published-mimetic-metrics-for-the-dgsem\/\">New preprint published:  Mimetic Metrics for the DGSEM<\/a> <span class=\"date\">April 7, 2025<\/span><div class=\"content\"><p>Our new preprint is available on arXiv. We explore an alternative approach to the Kopriva metric terms, which is based on finite element exterior calculus.<\/p>\n<p><a href=\"https:\/\/doi.org\/10.48550\/arXiv.2410.14502\">https:\/\/doi.org\/10.48550\/arXiv.2410.14502<\/a><\/p>\n<p><strong>Abstract<\/strong><br \/>\nFree-stream preservation is an essential property for numerical solvers on curvilinear grids. Key to this property is that the metric terms of the curvilinear mapping satisfy discrete metric identities, i.e., have zero divergence. Divergence-free metric terms are furthermore essential for entropy stability on curvilinear grids. We present a new way to compute the metric terms for discontinuous Galerkin spectral element methods (DGSEMs) that guarantees they are divergence-free. Our proposed mimetic approach uses projections that fit within the de Rham Cohomology.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/02\/10\/snapshot-barotropic-instability-with-and-without-orography\/\">Snapshot: Barotropic Instability with and without Orography<\/a> <span class=\"date\">February 10, 2025<\/span><div class=\"content\"><div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6506-4\" width=\"584\" height=\"343\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/02\/Barotropic_instability_topography_compressed.mp4?_=4\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/02\/Barotropic_instability_topography_compressed.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2025\/02\/Barotropic_instability_topography_compressed.mp4<\/a><\/video><\/div>\n<p>The video shows the evolution of a barotropic instability in Earth&#8217;s polar jet stream. Initially uniform, the jet stream undergoes perturbations, leading to vortex formation driven by the Coriolis force due to Earth&#8217;s rotation. To simulate this phenomenon, we discretize the Shallow Water Equations on the sphere using a discontinuous Galerkin method. We compare two scenarios: an ocean-covered Earth (without orography) and a more realistic representation that includes Earth&#8217;s orography. Orography is incorporated into the equations as a non-conservative term, with values sourced from the <a href=\"https:\/\/www.ncei.noaa.gov\/products\/etopo-global-relief-model\">ETOPO<\/a> dataset provided by the National Oceanic and Atmospheric Administration (NOAA). The simulations are performed using <a href=\"https:\/\/github.com\/trixi-framework\/TrixiAtmo.jl\">TrixiAtmo.jl<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2025\/01\/28\/new-paper-published-trixiparticles-jl-particle-based-multiphysics-simulation-in-julia\/\">New paper published: TrixiParticles.jl: Particle-based multiphysics simulation in Julia<\/a> <span class=\"date\">January 28, 2025<\/span><div class=\"content\"><p>Our new paper &#8220;TrixiParticles.jl: Particle-based multiphysics simulation in Julia&#8221; has been published in the Journal of Open Source Software.<br \/>\nWe are happy we were able to contribute to the publication and thank all our collaborators for the great experience.<\/p>\n<p><a href=\"https:\/\/doi.org\/10.21105\/joss.07044\"><br \/>\n    <img decoding=\"async\" src=\"https:\/\/joss.theoj.org\/papers\/10.21105\/joss.07044\/status.svg\" alt=\"DOI\"><br \/>\n<\/a><\/p>\n<h1>Summary<\/h1>\n<p>TrixiParticles.jl is a Julia-based open-source package for particle-based multiphysics simulations and part of the <a href=\"https:\/\/trixi-framework.github.io\/\">Trixi Framework<\/a>. It handles complex geometries and specialized applications, such as computational fluid dynamics (CFD) and structural dynamics, by providing a versatile platform for particle-based methods. TrixiParticles.jl allows for the straightforward addition of new particle systems and their interactions, facilitating the setup of coupled multiphysics simulations such as fluid-structure interaction (FSI). Furthermore, simulations are set up directly with Julia code, simplifying the integration of custom functionalities and promoting rapid prototyping<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/12\/20\/snapshot-hyperbolic-viscous-flow-three-dimensional-cubic-lid-driven-cavity\/\">Snapshot: Hyperbolic viscous flow &#8211; Three-dimensional cubic Lid-Driven Cavity<\/a> <span class=\"date\">December 20, 2024<\/span><div class=\"content\"><div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6483-5\" width=\"584\" height=\"329\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/12\/re1000_xz_short_compressed-1.mp4?_=5\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/12\/re1000_xz_short_compressed-1.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/12\/re1000_xz_short_compressed-1.mp4<\/a><\/video><\/div>\n<p>x-z central cross section of 3d Lid-Driven Cavity flow, solving the unified model of Godunov, Peshkov, and Romenski for hyperbolic viscous flow. The upper row shows the macroscopic flow state, while the bottom row shows the distortion field, tracking fluid flow deformations as if it were a solid. The numerical solver is a simple explicit second order Finite Volume method (MUSCL-Hancock) using an HLL-type Riemann solver based on Toro-Vazquez flux splitting. The Reynolds number is 1000 and the mesh resolution is 384^3.<\/p>\n<p>This research has been carried out by Simone Chiocchetti, funded by the European Union\u2019s Horizon Europe Research and Innovation Programme under the Marie Sk\u0142odowska-Curie Postdoctoral Fellowship MoMeNTUM (grant agreement No. 101109532).<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-300x63.png\" alt=\"\" width=\"300\" height=\"63\" class=\"aligncenter size-medium wp-image-6678\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-300x63.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-1024x215.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-768x161.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image-500x105.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2026\/01\/image.png 1100w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/11\/07\/snapshot-postdoctoral-researcher-dr-boqiang-huang-joins-our-research-group\/\">Snapshot: Postdoctoral researcher Dr. Boqiang Huang joins our research group<\/a> <span class=\"date\">November 7, 2024<\/span><div class=\"content\"><p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/11\/BQ_Huang.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/11\/BQ_Huang.jpg\" alt=\"\" width=\"517\" height=\"664\" class=\"aligncenter size-full wp-image-6472\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/11\/BQ_Huang.jpg 689w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/11\/BQ_Huang-233x300.jpg 233w\" sizes=\"(max-width: 517px) 100vw, 517px\" \/><\/a><\/p>\n<p>Boqiang Huang received his Ph.D. in Biomedical Engineering from the Department of Electronic Engineering at Fudan University, Shanghai, China, in 2010. Following his doctoral studies, he was awarded the Alexander von Humboldt Postdoctoral Fellowship and worked under the mentorship of Prof. Angela Kunoth. During this time, he contributed as a research scientist in Prof. Kunoth&#8217;s group (AG-Kunoth) at the Institute of Mathematics, University of Paderborn, and later at the University of Cologne. In 2021, Dr. Huang joined the group of Prof. Dorit Merhof (AG-Merhof) at the Institute of Imaging and Computer Vision at RWTH Aachen University, and subsequently at the University of Regensburg. His research expertise spans multiple disciplines in the field of data science, including biomedical engineering, electronic engineering, applied mathematics and applied physics.<\/p>\n<p>In July 2024, Dr. Huang became a research scientist in the group of Prof. Gregor Gassner (AG-Gassner). He is also supported by the HESCOR research project (&#8220;Human &#038; Earth System Coupled Research&#8221; <a href=\"https:\/\/hescor.uni-koeln.de\/\">https:\/\/hescor.uni-koeln.de\/<\/a>). His current research focuses on the &#8220;Machine learning &#038; Culture Clusters&#8221; over large timescales, ranging from the Paleolithic era (1.4 million to 24,000 years ago) to the present. He collaborates with researchers from archaeology, geophysics, and the humanities on this interdisciplinary project.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/08\/28\/snapshot-postdoctoral-researcher-dr-tristan-montoya-joins-our-research-group\/\">Snapshot: Postdoctoral researcher Dr. Tristan Montoya joins our research group<\/a> <span class=\"date\">August 28, 2024<\/span><div class=\"content\"><p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/01\/headshot_blackwhite-e1705574100657.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/01\/headshot_blackwhite-e1705574100657.jpg\" alt=\"\" width=\"400\" height=\"400\" class=\"aligncenter size-full wp-image-6233\" \/><\/a><br \/>\nTristan Montoya completed his PhD at the University of Toronto Institute for Aerospace Studies under the supervision of <a href=\"http:\/\/goldfinger.utias.utoronto.ca\/dwz\/\">Prof. David W. Zingg<\/a>. His doctoral work focused on the development and analysis of provably stable discontinuous spectral-element methods (DSEMs) with the summation-by-parts property for systems of conservation laws, particularly on triangular and tetrahedral meshes.<\/p>\n<p>Tristan joined the Numerical Simulation research group at the University of Cologne in January 2024, and is part of the <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a>  development team. His current research involves the development of split-form and entropy-stable DSEMs for atmospheric models and their application to global weather and climate prediction. You can keep up-to-date with Tristan\u2019s research and software contributions on his <a href=\"https:\/\/scholar.google.com\/citations?hl=en&#038;user=1nEXVtwAAAAJ\">Google Scholar<\/a> and <a href=\"https:\/\/github.com\/tristanmontoya\">GitHub<\/a> pages, as well as on his <a href=\"https:\/\/tjbmontoya.com\/\">personal website<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/07\/04\/snapshot-multiphysics-sph-simulation-of-the-liquefaction-of-the-julia-logo\/\">Snapshot: Multiphysics SPH simulation of the liquefaction of the Julia logo<\/a> <span class=\"date\">July 4, 2024<\/span><div class=\"content\"><p>2D single-phase simulation of the liquefaction of the Julia logo with <a href=\"https:\/\/github.com\/trixi-framework\/TrixiParticles.jl\">TrixiParticles.jl<\/a>. The letters are simulated as a fluid with the <a href=\"https:\/\/trixi-framework.github.io\/TrixiParticles.jl\/stable\/systems\/weakly_compressible_sph\/\">Weakly Compressible SPH  (WCSPH) method<\/a> and the dots are simulated as an elastic solid with the <a href=\"https:\/\/trixi-framework.github.io\/TrixiParticles.jl\/stable\/systems\/total_lagrangian_sph\/\">Total Lagrangian SPH (TLSPH) method<\/a> for structure mechanics.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/07\/julia_logo_2d.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/07\/julia_logo_2d.gif\" alt=\"\" width=\"880\" height=\"596\" class=\"aligncenter size-full wp-image-6389\" \/><\/a><\/p>\n<p><em>This Snapshot was created by Niklas Neher (<a href=\"https:\/\/www.hlrs.de\/de\/team\/niklas-neher\">https:\/\/www.hlrs.de\/de\/team\/niklas-neher<\/a>).<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/05\/31\/snapshot-discontinuous-galerkin-solver-for-the-spherical-shallow-water-equations\/\">Snapshot: Discontinuous Galerkin Solver for the Spherical Shallow Water Equations<\/a> <span class=\"date\">May 31, 2024<\/span><div class=\"content\"><p>A preliminary step in the development of a global atmospheric model is the construction of a suitable horizontal discretization for partial differential equations in spherical geometry. The spherical shallow water equations, which consist of a scalar equation governing mass conservation coupled with a vector equation for the momentum balance under gravitational and Coriolis forces, serve as a simplified model for the horizontal dynamics of the Earth&#8217;s atmopheric circulation. Since the spherical shallow water equations exhibit many of the characteristic features and numerical challenges associated with atmospheric fluid flow, they provide a useful testbed for the development and assessment of numerical schemes for weather prediction and climate modelling (see, for example, Williamson et al. [1]). <\/p>\n<p>By extending <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a> to solve hyperbolic partial differential equations on curved manifolds, we are able to simulate shallow-water flows on the surface of the sphere using discontinuous Galerkin methods formulated with respect to the two-dimensional tangent space associated with a cubed-sphere grid [2, 3]. The video below depicts the relative vorticity field for the numerical solution to the spherical shallow water equations in flux form, which we discretize similarly to Bao et al. [4] using a discontinuous Galerkin method employing 5400 curved quadrilateral elements of polynomial degree seven. The initial condition corresponds to a mid-latitude jet with a small perturbation added to initiate a barotropic instability, which was proposed by Galewsky et al. [5] as a test case exhibiting complex nonlinear dynamics characteristic of those present in numerical weather prediction and climate models.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6304-6\" width=\"584\" height=\"467\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/05\/spherical_shallow_water_barotropic_instability.mp4?_=6\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/05\/spherical_shallow_water_barotropic_instability.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/05\/spherical_shallow_water_barotropic_instability.mp4<\/a><\/video><\/div>\n<p><strong><br \/>\nReferences<\/strong><\/p>\n<p>[1] D. L. Williamson, J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber. A standard test set for numerical approximations to the shallow water equations in spherical geometry. Journal of Computational Physics, 102(1):211-224, 1992.<\/p>\n<p>[2] R. Sadourny. Conservative \ufb01nite-di\ufb00erence approximations of the primitive equations on quasi-uniform spherical grids. Monthly Weather Review, 100(2):136-144, 1972.<\/p>\n<p>[3] C. Ronchi, R. Iacono, and P. S. Paolucci. The &#8220;cubed sphere&#8221;: A new method for the solution of partial differential equations in spherical geometry. Journal of Computational Physics, 124(1):93-114, 1996.<\/p>\n<p>[4] L. Bao, R. D. Nair, and H. M. Tufo. A mass and momentum \ufb02ux-form high-order discontinuous Galerkin shallow water model on the cubed-sphere. Journal of Computational Physics, 271:224-243, 2014.<\/p>\n<p>[5] J. Galewsky, R. K. Scott, and L. M. Polvani. An initial-value problem for testing numerical models of the global shallow-water equations. Tellus A, 56(5):429\u2013440, 2004.<\/p>\n<p><em>This snapshot was created by <a href=\"https:\/\/tjbmontoya.com\/\">Tristan Montoya<\/a>.<\/em><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/04\/11\/snapshot-postdoctoral-researcher-dr-simone-chiocchetti-joins-our-research-group\/\">Snapshot: Postdoctoral researcher Dr. Simone Chiocchetti joins our research group<\/a> <span class=\"date\">April 11, 2024<\/span><div class=\"content\"><p><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/04\/SIMONE_F.png\" alt=\"\" width=\"512\" height=\"831\" class=\"aligncenter size-full wp-image-6289\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/04\/SIMONE_F.png 512w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/04\/SIMONE_F-185x300.png 185w\" sizes=\"(max-width: 512px) 100vw, 512px\" \/><\/p>\n<p>Simone Chiocchetti studied Civil, Environmental, and Mechanical Engineering at the  University of Trento, Italy, under the supervision of Prof. Dr.-Ing. Michael Dumbser.<\/p>\n<p>He spent a year in Stuttgart at the IAG (Institute of Aerodynamics and Gasdynamics), funded by the DFG project &#8220;DROPIT startup grant&#8221; and he is currently working at the University of Cologne in the Numerical Simulation Research group led by Prof. Dr.-Ing. Gregor Gassner, under a Marie Sk\u0142odowska-Curie Postdoctoral Fellowship (European Union&#8217;s Horizon Europe Research and Innovation Programme, grant agreement No. 101109532).<\/p>\n<p>Simone&#8217;s research interests include: high order numerical methods, of the Finite Volume and Discontinuous Galerkin families, in the context of first order hyperbolic partial differential equations; algorithms for generation and optimization of high quality unstructured meshes; treatment of stiff algebraic source terms; modeling of fluid and solid mechanics as well as multiphase flow;  high-performance computing; semi-implicit numerical methods.<\/p>\n<p>Simone&#8217;s current MSCA project (MoMeNTUM: <em>Modern high order numerical Methods based on No-compromise moving Voronoi Tessellations, a Unified solver for continuum Mechanics<\/em>) concerns the development of high order numerical methods on Voronoi grids, together with new meshing algorithms that allow to fully take advantage of the versatility provided by unstructured moving polygonal meshes, with special attention to the computational efficiency of the implementation.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/03\/28\/snapshot-postdoctoral-researcher-dr-aleksey-sikstel-joins-our-research-group\/\">Snapshot: Postdoctoral researcher Dr. Aleksey Sikstel joins our research group<\/a> <span class=\"date\">March 28, 2024<\/span><div class=\"content\"><p><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/03\/photo_bw-5-e1678111188648.jpeg\" alt=\"\" width=\"400\" height=\"400\" class=\"aligncenter size-full wp-image-5945\" \/><br \/>\nAleksey defended his doctoral thesis on coupling of hyperbolic conservation laws in 2020 at the RWTH Aachen University while working with Prof. Siegfried M\u00fcller and Prof. Michael Herty. After that he continued working with Prof. Jan Giesselmann at the TU Darmstadt on a posteriori error estimators for systems of hyperbolic conservation laws before switching to Cologne in 2022.<\/p>\n<p>Alekseys research focus lies on adaptivity for high-order numerical methods for hyperbolic conservation laws, coupling of conservation laws via sharp interfaces, multiphase flows and stochastic-Galerkin PDE. In the Numerical Simulation Research Group Aleksey is working on multi-ion MHD systems and novel coupling and temporal schemes within the Snubic-project.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/02\/21\/snapshot-postdoctoral-researcher-dr-benedict-geihe-joins-our-research-group\/\">Snapshot: Postdoctoral researcher Dr. Benedict Geihe joins our research group<\/a> <span class=\"date\">February 21, 2024<\/span><div class=\"content\"><p><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/09\/DSC_1337_square.jpg\" alt=\"\" width=\"400\" height=\"400\" class=\"aligncenter size-full wp-image-6093\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/09\/DSC_1337_square.jpg 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/09\/DSC_1337_square-300x300.jpg 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/09\/DSC_1337_square-150x150.jpg 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/09\/DSC_1337_square-45x45.jpg 45w\" sizes=\"(max-width: 400px) 100vw, 400px\" \/><\/p>\n<p>Benedict Geihe studied Mathematics and Computer Science at the University of Bonn, Germany. After graduation he joined the research group of Professor Martin Rumpf at the Institute for Numerical Simulation as a research and teaching assistant. In his PhD studies he investigated and implemented finite element based numerical simulation and shape optimization of multi-scale microstructures in elastic composite materials.<\/p>\n<p>Benedict then moved to the German Aerospace Center (DLR) and worked as a researcher at the Numerical Methods department of the Institute of Propulsion Technology. He was part of the development team of the flow solver TRACE, used for academic and industrial simulation and analysis of flows in turbomachines. His research covered resonance phenomena induced by fluid structure interaction and frequency based numerical methods.<\/p>\n<p>Benedict is now part of the Numerical Simulation Research Group at the University of Cologne, where he joined the <a href=\"https:\/\/trixi-framework.github.io\/\">Trixi.jl<\/a> development team. He is currently working on project <a href=\"https:\/\/dlr-amr.github.io\/adaptex\/\">ADAPTEX<\/a>, funded by the German Federal Ministry of Education and Research within the SCALEXA initiative. The objective is to deploy Trixi.jl in exascale-capable flow simulations on CPUs and GPUs for applications in earth system modeling on dynamic adaptive meshes.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2024\/01\/26\/snapshot-2d-sloshing-simulation\/\">Snapshot: 2D Sloshing simulation<\/a> <span class=\"date\">January 26, 2024<\/span><div class=\"content\"><p>2D Sloshing simulation in a round tank oscillating from left to right. The simulation was done with the CFD code <a href=\"https:\/\/trixi-framework.github.io\/\">Trixi.jl<\/a> [1] which was extended with a three-equations model [2]. The model consists of a two-phase flow simplification of the Baer\u2013Nunziato equations of compressible multi-phase flows [3]. The shown simulation is part of combined research efforts of the German Aerospace Center (DLR) and the University of Cologne to create computational models of sloshing dynamics in liquid hydrogen storage systems on mobile platforms such as cars, ships, and planes.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6240-7\" width=\"584\" height=\"389\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/01\/sloshing_round_tank_publish1.mp4?_=7\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/01\/sloshing_round_tank_publish1.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2024\/01\/sloshing_round_tank_publish1.mp4<\/a><\/video><\/div>\n<p><strong><br \/>\nReferences:<\/strong><br \/>\n[1] <a href=\"https:\/\/trixi-framework.github.io\/\">https:\/\/trixi-framework.github.io\/<\/a><br \/>\n[2] Dumbser, Michael. &#8220;A simple two-phase method for the simulation of complex free surface flows.&#8221; Computer methods in applied mechanics and engineering 200.9-12 (2011): 1204-1219.<br \/>\n[3] M.R. Baer, J.W. Nunziato: &#8220;A two-phase mixture theory for the deflagration-todetonation transition (DDT) in reactive granular materials&#8221;, J. Multiphase Flow 12 (1986) 861\u2013889.<br \/>\n[4] Project HyTaZer: <a href=\"https:\/\/elib.dlr.de\/201347\/\">https:\/\/elib.dlr.de\/201347\/<\/a> and <a href=\"https:\/\/www.dlr.de\/sy\/PortalData\/17\/Resources\/dokumente\/wissenschaftstag\/2022\/6_Wissenschaftstag_2022_Hytazer_Freund.pdf\">https:\/\/www.dlr.de\/sy\/PortalData\/17\/Resources\/dokumente\/wissenschaftstag\/2022\/6_Wissenschaftstag_2022_Hytazer_Freund.pdf<\/a><\/p>\n<p><em>This Snapshot was created by Johannes Markert (<a href=\"http:\/\/jmark.de\/\">http:\/\/jmark.de\/<\/a>).<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/12\/01\/snapshot-discontinuous-galerkin-simulation-on-a-spherical-shell\/\">Snapshot: Discontinuous Galerkin simulation on a spherical shell<\/a> <span class=\"date\">December 1, 2023<\/span><div class=\"content\"><p>Utilizing a mapping $(\\xi, \\eta) \\in \\mathbb{R}^2 \\rightarrow (x, y, z) \\in \\mathbb{R}^3$ and specifically tailored tensor-product Legendre&#8211;Gauss&#8211;Lobatto basis functions [1,2], discontinuous Galerkin (DG) simulations can be performed on a curved surface. In particular, we are interested in running DG simulations on a spherical shell.<\/p>\n<p>To achieve this, we create a two-dimensional cubed-sphere mesh for tessellating the sphere&#8217;s surface. This mesh proves advantageous as it avoids singularities at the poles present in latitude-longitude grids, while still facilitating a highly regular tessellation of the simulation domain through the use of quadrilaterals.<\/p>\n<p>Within this mesh, we address the linear advection equations incorporating position-dependent advection velocity:<br \/>\n\\[ \\frac{\\partial \\rho}{\\partial t} + \\nabla \\cdot \\left( \\vec{v} (x, y, z) \\rho \\right) = 0. \\]<br \/>\nAdditionally, a solid-body rotation velocity tangential to the spherical surface is imposed.<\/p>\n<p>The presented video showcases the advection of a Gaussian pulse across the globe using Trixi.jl (<a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">https:\/\/github.com\/trixi-framework\/Trixi.jl<\/a>) with the two-dimensional p4est solver. The simulation also incorporates a custom implementation of the three-dimensional linear advection equations.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6175-8\" width=\"584\" height=\"359\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/12\/earth.mp4?_=8\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/12\/earth.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/12\/earth.mp4<\/a><\/video><\/div>\n<p><strong><br \/>\nReferences<\/strong><\/p>\n<p>[1] Song, C. &#038; Wolf, J. P. (1999). The scaled boundary finite element method\u2014alias consistent infinitesimal finite element cell method\u2014for diffusion. International Journal for Numerical Methods in Engineering, 45(10), 1403-1431.<br \/>\n[2] Giraldo, F. X. (2001). A spectral element shallow water model on spherical geodesic grids. International Journal for Numerical Methods in Fluids, 35(8), 869-901.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/10\/27\/snapshot-blending-finite-volume-fluxes-with-reinforcement-learning\/\">Snapshot: Blending Finite Volume Fluxes with Reinforcement Learning<\/a> <span class=\"date\">October 27, 2023<\/span><div class=\"content\"><p>When solving PDEs with the finite volume method, one must choose a numerical flux function. We use a convex combination of the central flux $F^C$ and local Lax Friedrich flux $F^{LLF}$ to solve the Burgers equation. <\/p>\n<p>\\begin{align*}<br \/>\n&#038;u_t + \\left(\\frac{u^2}{2}\\right)_x =0\\\\<br \/>\n&#038;\\frac{du_j}{dt} =\\frac{1}{\\Delta x} \\left[ \\alpha_j \\left(F_{j-1\/2}^{LLF} -F_{j+1\/2}^{LLF}\\right) + (1-\\alpha_j)\\left(F_{j-1\/2}^{C} &#8211; F_{j+1\/2}^{C}\\right) \\right]<br \/>\n\\end{align*}<\/p>\n<p>The central flux is highly accurate but can lead to oscillating and thus unstable solutions. The local Lax Friedrich flux is stable but also dissipative. Therefore, we want to choose alpha so that the solution is stable and at the same time as accurate as possible. We train a Reinforcement Learning Agent to choose alpha. One agent performs locally in one grid cell. The action of agent j is $\\alpha_j$ and the state is given by<\/p>\n<p>\\begin{equation*}<br \/>\nr_j=\\frac{u_j &#8211; u_{j-1}}{u_{j+1} &#8211; u_{j}}.<br \/>\n\\end{equation*}<\/p>\n<p>The policy is a Neural Network with one hidden layer of size 10, a relu activation function in the first layer, and a hard sigmoid function in the output layer. The deep deterministic policy gradient (DDPG) approach is used to train the agent.<\/p>\n<p>The video below shows the chosen amount of local Lax Friedrich flux and the resulting solution of the Burgers equation. The initial solution is a sine wave with two discontinuities. Randomly changed variants of this were used for training.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/10\/sinus_square.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/10\/sinus_square.gif\" alt=\"\" width=\"600\" height=\"400\" class=\"aligncenter size-full wp-image-6168\" \/><\/a><\/p>\n<p>The agent performs well in unseen situations, e.g. negative sine wave, as shown below.<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/10\/minus_sinus.gif\" alt=\"\" width=\"600\" height=\"400\" class=\"aligncenter size-full wp-image-6170\" \/><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/09\/29\/open-postdoctoral-positions-in-scientific-computing-and-machine-learning\/\">Open Postdoctoral Positions in Scientific Computing and Machine Learning<\/a> <span class=\"date\">September 29, 2023<\/span><div class=\"content\"><p><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">In the framework of the new project on \u201cCultural Evolution in Changing Climate: Human and Earth System Coupled Research (HESCOR)\u201d, we invite applications for postdoctoral research associate positions in the area of Earth system modeling, Human system modeling, applied numerical mathematics, scientific computing and\/or machine learning, among many others. <\/span><\/span><\/p>\n<p><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">The job description<\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">s<\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">, duty statement and qualification requirements for each of these positions can be found<\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\"> on <\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">the website<\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">s<\/span><\/span><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">:<br \/>\n<a href=\"https:\/\/imfess.uni-koeln.de\/hescor\/positions\">https:\/\/imfess.uni-koeln.de\/hescor\/positions<\/a> and<br \/>\n<a href=\"https:\/\/ufg.phil-fak.uni-koeln.de\/hescor\">https:\/\/ufg.phil-fak.uni-koeln.de\/hescor<\/a><\/span><\/span><\/p>\n<p><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">The research associate <\/span><\/span><strong><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">RA6 (RA in Scientific Computing) <\/span><\/span><\/strong><strong><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">will be hosted at the Numerical Simulation Research Group of the University of Cologne<\/span><\/span><\/strong><strong><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">.<\/span><\/span><\/strong><\/p>\n<p><span style=\"font-family: Sansserif, sans-serif;\"><span style=\"font-size: small;\">We provide a collaborative and friendly environment for cutting edge research and ensure measures are in place to promote early-career researchers.<\/span><\/span><\/p>\n<p>&nbsp;<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/09\/04\/studentische-wissenschaftliche-hilfskraft-w-m-d-zur-betreuung-und-weiterentwicklung-einer-automatisierten-testumgebung-fur-den-stromungsloser-trixi-gesucht\/\">Studentische\/Wissenschaftliche Hilfskraft (w\/m\/d) zur Betreuung und Weiterentwicklung einer automatisierten Testumgebung f\u00fcr den Str\u00f6mungsl\u00f6ser Trixi gesucht<\/a> <span class=\"date\">September 4, 2023<\/span><div class=\"content\"><p>In der Arbeitsgruppe von Prof. Gassner besch\u00e4ftigen wir uns mit mehrskaligen Str\u00f6mungsproblemen in verschiedenen physikalischen Disziplinen. Um numerisch effizient simulieren zu k\u00f6nnen, wurde in den vergangenen Jahren der Str\u00f6mungsl\u00f6ser Trixi entwickelt, basierend auf der modernen und nutzerfreundlichen Programmiersprache Julia. Die Arbeit der zahlreichen Entwickler wird auf github koordiniert und erm\u00f6glicht automatisierte Prozesse im Rahmen von continuous integration, testing und benchmarking. Des Weiteren wurde k\u00fcrzlich mit der Entwicklung einer C\/Fortran-Schnittstelle begonnen, die k\u00fcnftig interessierten Anwender m\u00f6glichst einfach zur Verf\u00fcgung gestellt werden soll.<\/p>\n<p>F\u00fcr die Wartung und Weiterentwicklung der automatisierten Prozesse sowie der Schnittstelle suchen wir ab sofort interessierte Hilfskr\u00e4fte. Technische Vorkenntnisse im Bereich des Software-Engineering unter Linux sind vorteilhaft, aber nicht zwingend n\u00f6tig. Die Arbeit bietet Gelegenheit, sich in die Thematik einzuarbeiten. Der Starttermin sowie die Wochenstundenzahl k\u00f6nnen flexibel vereinbart werden. Bei Interesse wenden Sie sich gerne an\u00a0<a class=\"c-link\" href=\"mailto:bgeihe@uni-koeln.de\" target=\"_blank\" rel=\"noopener noreferrer\" data-stringify-link=\"mailto:bgeihe@uni-koeln.de\" data-sk=\"tooltip_parent\" aria-haspopup=\"menu\" aria-expanded=\"false\">bgeihe@uni-koeln.de<\/a>, so dass wir Ihnen vor Ort weitere Details vorstellen k\u00f6nnen.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/08\/25\/snapshot-second-order-fv-method-on-hybrid-mesh-using-t8code\/\">Snapshot: Second-order FV method on hybrid mesh using t8code<\/a> <span class=\"date\">August 25, 2023<\/span><div class=\"content\"><p>This FV simulation uses linear reconstruction to reach second order. The minmod slope limiter is used. The mesh contains triangles and quads and is a high-resolution version (with 327680 cells) of this one:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3.png\" alt=\"\" width=\"896\" height=\"811\" class=\"aligncenter size-full wp-image-6085\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3.png 896w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3-300x272.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3-768x695.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/mesh_level3-331x300.png 331w\" sizes=\"(max-width: 896px) 100vw, 896px\" \/><\/a><\/p>\n<p>It is created and organized using <a href=\"https:\/\/github.com\/DLR-AMR\/t8code\" rel=\"noopener\" target=\"_blank\">t8code<\/a>.<\/p>\n<p>The simulation uses an initial condition with a pressure blast wave in the center and periodic boundaries.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6084-9\" width=\"584\" height=\"469\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/FV2_level8_hybrid_blast.mp4?_=9\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/FV2_level8_hybrid_blast.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/08\/FV2_level8_hybrid_blast.mp4<\/a><\/video><\/div>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/07\/28\/snapshot-sph-simulation-of-an-elastic-ball-falling-into-a-tank-of-water\/\">Snapshot: SPH Simulation of an Elastic Ball Falling Into a Tank of Water<\/a> <span class=\"date\">July 28, 2023<\/span><div class=\"content\"><p>This simulation uses weakly compressible smoothed particle hydrodynamics for the fluid and a total Lagrangian SPH formulation for the elastic solid.<br \/>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6080-10\" width=\"584\" height=\"584\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/07\/falling_sphere_3d_negative_pressure.mp4?_=10\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/07\/falling_sphere_3d_negative_pressure.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/07\/falling_sphere_3d_negative_pressure.mp4<\/a><\/video><\/div><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/06\/02\/snapshot-subcell-limiting-results-for-dgsem-simulations\/\">Snapshot: Subcell limiting results for DGSEM simulations<\/a> <span class=\"date\">June 2, 2023<\/span><div class=\"content\"><p>DGSEM simulations of the two-dimensional Euler equations considering the astrophysical jet with Mach number $\\textrm{Ma} \\approx 2000$ [1].<br \/>\nThe computational domain, $\\Omega = [-0.5,0.5]^2$, is filled with a mono-atomic gas ($\\gamma = 5\/3$) at rest with<\/p>\n<p>\\begin{equation}<br \/>\n\\rho(x,y) = 0.5, \\qquad<br \/>\np(x,y) = 0.4127, \\qquad<br \/>\nv_1(x,y) = 0, \\qquad<br \/>\nv_2(x,y) = 0,<br \/>\n\\end{equation}<\/p>\n<p>and on the left boundary there is a hypersonic inflow with<\/p>\n<p>\\begin{equation}<br \/>\n\\rho(-0.5,y_B) = 5, ~<br \/>\np(-0.5,y_B) = 0.4127, ~<br \/>\nv_1(-0.5,y_B) = 800, ~<br \/>\nv_2(-0.5,y_B) = 0,<br \/>\n\\end{equation}<\/p>\n<p>for $y_B \\in [-0.05, 0.05]$, which corresponds to a Mach number of $\\textrm{Ma}=2156.91$ with respect to the speed of sound of the jet gas, and $\\textrm{Ma}=682.08$ with respect to the speed of sound of the ambient gas.<\/p>\n<p>A spatial resolution of $256\\times 256$ quadrilateral elements, a polynomial degree of $N=3$, and four different CFL numbers are used.<\/p>\n<p>The simulations are stabilized by using two different subcell limiting approaches, the monolithic convex limiting (MCL) and the FCT\/IDP limiting.<\/p>\n<p>The resulting density contours at $t=10^{-3}$:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-570x1024.jpg\" alt=\"\" width=\"570\" height=\"1024\" class=\"aligncenter size-large wp-image-6072\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-570x1024.jpg 570w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-167x300.jpg 167w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet-768x1381.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/AstroJet.jpg 850w\" sizes=\"(max-width: 570px) 100vw, 570px\" \/><\/a><\/p>\n<p>The dependence of the spatial discretization on the time-step size for FCT\/IDP methods causes the number of vortical structures to be highly dependent on the CFL number and the total number of time steps to be not inversely proportional to the CFL number.<\/p>\n<p>In fact, the amount of dissipation is reduced for small CFL numbers, which leads to lower minimum densities and higher maximum pressures in FCT\/IDP, as can be seen in this Table.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg\" alt=\"\" width=\"962\" height=\"378\" class=\"aligncenter size-full wp-image-6074\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table.jpg 962w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-300x118.jpg 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-768x302.jpg 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Table-500x196.jpg 500w\" sizes=\"(max-width: 962px) 100vw, 962px\" \/><\/a><\/p>\n<p><strong><br \/>\nReferences<\/strong><br \/>\n[1] Rueda-Ram\u00edrez, A., Bolm, B., Kuzmin, D. &#038; Gassner, G. (2023). Monolithic Convex Limiting for Legendre-Gauss-Lobatto Discontinuous Galerkin Spectral Element Methods. arXiv preprint arXiv:2303.00374<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/06\/02\/thesis-snapshot-a-wet-dry-treatment-to-numerically-solve-the-shallow-water-equations-with-a-high-order-discontinuous-galerkin-method-with-trixi-jl-masters-thesis-by-sven-goldberg\/\">Thesis Snapshot: A wet-dry treatment to numerically solve the shallow water equations with a high-order discontinuous Galerkin method with Trixi.jl (Master&#8217;s Thesis by Sven Goldberg)<\/a> <span class=\"date\">June 2, 2023<\/span><div class=\"content\"><p>The shallow water equations are a well-known and often used physical model to simulate shallow water flow within given domains such as oceans. They are a system of hyperbolic partial differential equations and consist of conservation equations.<\/p>\n<p>The numerical solver <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a> is used to find approximate solutions to the equations. The shallow water equations are already embedded in this framework, but Trixi.jl does not yet allow for the appearance of dry sub-regions. Some strategies are combined and implemented to resolve the wet-dry problem resulting in a stable and well-balanced scheme: We use the positivity-preserving limiter by Zhang and Shu [1], the hydrostatic reconstruction method by Chen and Noelle [2] and a strategy to identify and mark dry cells.<\/p>\n<p>A standard test exists to check for accuracy, namely the parabolic bowl test by Niklas Wintermeyer [3], first presented by William Thacker [4]. We use a StructuredMesh with 75625 elements and polynomial degree five. Since it has periodic analytical solutions, we can show the precision by realizing the initial state and its linear structure are recovered after one period. In the left picture, the whole domain is shown. In the right image, a slice at $y=0$ is visualized. Both are snapshots after one period:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-1024x629.png\" alt=\"\" width=\"584\" height=\"359\" class=\"aligncenter size-large wp-image-6055\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-1024x629.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-300x184.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-768x472.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-1536x943.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic-488x300.png 488w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/005_2D_parabolic_bowl_periodic.png 2032w\" sizes=\"(max-width: 584px) 100vw, 584px\" \/><\/a><\/p>\n<p>We also use the <a href=\"https:\/\/github.com\/trixi-framework\/TrixiBottomTopography.jl\">TrixiBottomTopography.jl<\/a> package by Maximilian Bertrand to account for real-world bottom topography data from the\u00a0<a href=\"https:\/\/www.opengeodata.nrw.de\/produkte\/geobasis\/hm\/dgm1_xyz\/dgm1_xyz\/\">DGM1<\/a>\u00a0data set. We give the simulation of flooding on the Rhine River valley using a TreeMesh with 1024 elements polynomial degree six and outflow boundary conditions:<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6048-11\" width=\"584\" height=\"355\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Rhine_animation.mp4?_=11\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Rhine_animation.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/06\/Rhine_animation.mp4<\/a><\/video><\/div>\n<p><strong><br \/>\nReferences:<\/strong><br \/>\n[1] Zhang, X. and Chi-Wang, S. (2011). Maximum-principle-satisfying and positivity-preserving high-order schemes for conservation laws: survey and new developments. Proc. R. Soc. A. 467:2752\u20132776. <a href=\"http:\/\/doi.org\/10.1098\/rspa.2011.0153\">http:\/\/doi.org\/10.1098\/rspa.2011.0153<\/a><br \/>\n[2] Chen, G. and Noelle, S. (2017). A New Hydrostatic Reconstruction Scheme Based on Subcell Reconstructions. SIAM Journal on Numerical Analysis, 55, 758-784. <a href=\"https:\/\/doi.org\/10.1137\/15M1053074\">https:\/\/doi.org\/10.1137\/15M1053074<\/a><br \/>\n[3] Wintermeyer, N.(2018). A novel entropy stable discontinuous Galerkin spectral element method for the shallow water equations on GPUs. Dissertation, Universit\u00e4t zu K\u00f6ln. <a href=\"https:\/\/kups.ub.uni-koeln.de\/9234\/\">https:\/\/kups.ub.uni-koeln.de\/9234\/<\/a><br \/>\n[4] Thacker, W. (1981). Some exact solutions to the nonlinear shallow-water wave equations. Journal of Fluid Mechanics, 107, 499-508. <a href=\"https:\/\/doi.org\/10.1017\/S0022112081001882\">doi:10.1017\/S0022112081001882<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/05\/04\/snapshot-subcell-limiting-strategies-for-the-dgsem-with-legendre-gauss-nodes\/\">Snapshot: Subcell limiting strategies for the DGSEM with Legendre&#8211;Gauss nodes<\/a> <span class=\"date\">May 4, 2023<\/span><div class=\"content\"><p>The discontinuous Galerkin spectral element method (DGSEM) is a nodal discontinuous Galerkin collocation scheme that uses a set of Legendre-Gauss or Legendre-Gauss-Lobatto nodes to represent the solution with Lagrange interpolating polynomials and to compute integrals with a quadrature rule. <\/p>\n<p>It was recently shown [1] that the DGSEM semi-discretization of a conservation law using Legendre-Gauss nodes can be written in finite-volume form as<br \/>\n\\begin{equation*}<br \/>\nJ_{j} \\dot{\\mathbf{u}}^{DG}_{j} =<br \/>\n\\frac{1}{\\omega_j}<br \/>\n\\left(<br \/>\n  \\hat{\\mathbf{f}}^{DG}_{(j-1,j)}<br \/>\n&#8211; \\hat{\\mathbf{f}}^{DG}_{(j,j+1)}<br \/>\n\\right),<br \/>\n\\end{equation*}<br \/>\nwhere $J_j$ is the mapping Jacobian, $\\omega_j$ are the Legendre-Gauss quadrature weights, and $\\hat{\\mathbf{f}}^{DG}_{(\\cdot,\\cdot)}$ are high-order DGSEM fluxes between adjacent nodes.<\/p>\n<p>The existence of a FV form for the Legendre&#8211;Gauss DGSEM enables the use of the subcell limiting strategies described in [2]. In these methods, a hybrid DGSEM\/FV method is used, where the semi-discretization at each node reads<br \/>\n\\begin{equation*}<br \/>\nJ_{j} \\dot{\\mathbf{u}}_{j} =<br \/>\n\\frac{1}{\\omega_j}<br \/>\n\\left(<br \/>\n  \\hat{\\mathbf{f}}_{(j-1,j)}<br \/>\n&#8211; \\hat{\\mathbf{f}}_{(j,j+1)}<br \/>\n\\right),<br \/>\n\\end{equation*}<br \/>\nand the local fluxes are obtained as a convex combination of high-order DGSEM and low-order FV fluxes,<br \/>\n\\begin{equation*}<br \/>\n    \\hat{\\mathbf{f}}_{(i,j)} = (1-\\alpha_{(i,j)}) \\hat{\\mathbf{f}}^{DG}_{(i,j)}<br \/>\n    +<br \/>\n    \\alpha_{(i,j)}<br \/>\n    \\hat{\\mathbf{f}}^{DG}_{(i,j)},<br \/>\n    ~~~~<br \/>\n    \\alpha_{(i,j)} \\in [0,1].<br \/>\n\\end{equation*}<\/p>\n<p>To illustrate the shock-capturing capacity of the hybrid DGSEM\/FV method, we simulate a Sedov blast problem describing the evolution of a blast wave expanding from an initial concentration of density and pressure, and adjust the blending coefficient $\\alpha_{(i,j)}$ to avoid non-physical density oscillations in the vicinity of shocks.<\/p>\n<p>The video shows the distribution of density and the blending coefficient obtained for entropy-stable variants of the Legendre-Gauss and Legendre-Gauss-Lobatto DGSEM methods.<\/p>\n<div style=\"width: 584px;\" class=\"wp-video\"><video class=\"wp-video-shortcode\" id=\"video-6004-12\" width=\"584\" height=\"362\" preload=\"metadata\" controls=\"controls\"><source type=\"video\/mp4\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/subcell_wise_gauss_vs_lgl.mp4?_=12\" \/><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/subcell_wise_gauss_vs_lgl.mp4\">https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/subcell_wise_gauss_vs_lgl.mp4<\/a><\/video><\/div>\n<p>The simulation runs stably without any spurious oscillations for both variants of the DGSEM method.<br \/>\nThe figures below show the mean blending coefficient and the number of time steps for the simulation of the blast wave.<br \/>\nAlthough the Legendre&#8211;Gauss DGSEM requires more limiting than the Legendre-Gauss-Lobatto DGSEM to avoid spurious density oscillations, it completes the simulation in fewer time steps.<br \/>\nThe Legendre-Gauss subcell distribution allows for longer time-step sizes than the Legendre-Gauss-Lobatto subcell distribution for the same CFL number.<\/p>\n<p><div id=\"attachment_6023\" style=\"width: 275px\" class=\"wp-caption alignleft\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/alpha_bar.png\"><img aria-describedby=\"caption-attachment-6023\" decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/alpha_bar-300x225.png\" alt=\"\" width=\"300\" height=\"225\" class=\"size-medium wp-image-6023\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/alpha_bar-300x225.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/alpha_bar-400x300.png 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/alpha_bar.png 640w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><p id=\"caption-attachment-6023\" class=\"wp-caption-text\">Mean blending coefficient<\/p><\/div><div id=\"attachment_6022\" style=\"width: 275px\" class=\"wp-caption alignright\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/time_steps.png\"><img aria-describedby=\"caption-attachment-6022\" decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/time_steps-300x225.png\" alt=\"\" width=\"300\" height=\"225\" class=\"size-medium wp-image-6022\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/time_steps-300x225.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/time_steps-400x300.png 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/05\/time_steps.png 640w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><p id=\"caption-attachment-6022\" class=\"wp-caption-text\">Number of time steps<\/p><\/div><\/p>\n<p>Figure 1: Evolution of the mean blending coefficient and number of time steps taken for the simulation of the blast wave with CFL=0.9.<\/p>\n<p><strong><br \/>\nReferences<\/strong><br \/>\n[1] Mateo-Gab\u00edn, A., Rueda-Ram\u00edrez, A. M., Valero, E., &#038; Rubio, G. (2022). Entropy-stable flux-differencing formulation with Gauss nodes for the DGSEM. arXiv preprint arXiv:2211.05066.<br \/>\n[2] Rueda-Ram\u00edrez, A. M., Pazner, W., &#038; Gassner, G. J. (2022). Subcell limiting strategies for discontinuous Galerkin spectral element methods. Computers &#038; Fluids, 247, 105627. arXiv preprint arXiv:2202.00576.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/04\/03\/snapshot-sph-simulation-of-spheres-falling-into-a-tank-of-water\/\">Snapshot: SPH simulation of spheres falling into a tank of water<\/a> <span class=\"date\">April 3, 2023<\/span><div class=\"content\"><p>SPH simulation of two elastic spheres with different densities falling into a tank of water:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/04\/out-1.gif\"><img decoding=\"async\" loading=\"lazy\" class=\"aligncenter size-full wp-image-5992\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/04\/out-1.gif\" alt=\"\" width=\"600\" height=\"226\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/04\/03\/new-paper-published-reliable-genotyping-of-recombinant-genomes-using-a-robust-hidden-markov-model\/\">New paper published:  Reliable genotyping of recombinant genomes using a robust hidden Markov model<\/a> <span class=\"date\">April 3, 2023<\/span><div class=\"content\"><p>Meiotic recombination is an essential mechanism during sexual reproduction and includes the exchange of chromosome segments between homologous chromosomes. New allelic combinations are transmitted to the new generation, introducing novel genetic variation in the offspring genomes. With the improvement of high-throughput whole-genome sequencing technologies, large numbers of recombinant individuals can now be sequenced with low sequencing depth at low costs, necessitating computational methods for reconstructing their haplotypes. The main challenge is the uncertainty in haplotype calling that arises from the low information content of a single genomic position. Straightforward sliding window-based approaches are difficult to tune and fail to place recombination breakpoints precisely. Hidden Markov Model-based approaches, on the other hand, tend to over-segment the genome. Here, we present RTIGER, an HMM-based model that exploits in a mathematically precise way the fact that true chromosome segments typically have a certain minimum length. We further separate the task of identifying the correct haplotype sequence from the accurate placement of haplotype borders, thereby maximizing the accuracy of border positions. By comparing segmentations based on simulated data with known underlying haplotypes, we highlight the reasons for RTIGER outperforming traditional segmentation approaches. We then analyze the meiotic recombination pattern of segregants of two Arabidopsis (Arabidopsis thaliana) accessions and a previously described hyper-recombining mutant. RTIGER is available as an R-package with an efficient Julia implementation of the core algorithm.<\/p>\n<p>Plants Physiology: <a href=\"https:\/\/doi.org\/10.1093\/plphys\/kiad191\">https:\/\/doi.org\/10.1093\/plphys\/kiad191<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/02\/22\/snapshot-test-simulation-with-novel-wet-dry-dg-implementation\/\">Snapshot: Test Simulation with novel wet\/dry DG implementation<\/a> <span class=\"date\">February 22, 2023<\/span><div class=\"content\"><p>We solve the shallow water equations with DGSEM on a mesh with 103 elements and polynomial degree 12 in combination with subcell shock capturing and limiting, positivity preservation and handling of wet\/dry regions.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/better_flood.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/better_flood.gif\" alt=\"\" width=\"640\" height=\"480\" class=\"aligncenter size-full wp-image-5895\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/02\/03\/thesis-snapshot-ameisenalgorithmen-metaheuristische-verfahren-zur-losung-des-travelling-salesman-problems-bachelorarbeit-von-marie-becker\/\">Thesis Snapshot: Ameisenalgorithmen \u2013 Metaheuristische Verfahren zur L\u00f6sung des Travelling Salesman Problems (Bachelorarbeit von Marie Becker)<\/a> <span class=\"date\">February 3, 2023<\/span><div class=\"content\">\r\n<div class=\"wp-block-columns is-layout-flex wp-container-2 wp-block-columns-is-layout-flex\">\r\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis: 100%;\">\r\n<p class=\"wp-block-preformatted\">In den neunziger Jahren wurden erstmals Ameisenalgorithmen von Marco Dorigo vorgestellt. Der erste Algorithmus war das sogenannte Ant System (AS), welches in der Bachelorarbeit implementiert wurde, um das Travelling Salesman Problem zu l\u00f6sen. Dabei wird das Verhalten von Ameisen in einen Algorithmus \u00fcbertragen, um den k\u00fcrzesten Hamiltonkreis in einem Graphen zu finden.\r\n\r\nDie Ameisenalgorithmen bauen auf dem Prinzip auf, dass jede virtuelle Ameise, ausgehend von einem zuf\u00e4lligen Startknoten, einen Hamiltonkreis in dem gegebenen Graphen abl\u00e4uft und auf dieser Tour dann Duftstoffe, sogenannte Pheromone, hinterl\u00e4sst. Die hinterlassene Menge der Pheromone h\u00e4ngt dabei davon ab, wie lang die Tour einer Ameise ist. Durch die Duftstoffe werden andere Ameisen dazu angeregt, auch diesen Weg zu w\u00e4hlen. Je mehr Pheromone sich auf einer Kante befinden, desto wahrscheinlicher ist es, dass eine Ameise diese Kante in ihrer Tour w\u00e4hlt. Nach einigen Iterationen befindet sich also die gr\u00f6\u00dfte Menge der Pheromone auf dem k\u00fcrzesten Hamiltonkreis.\r\n\r\nEine Erweiterung dieses Algorithmus ist das Ant Colony System (ACS), welches ebenfalls auf das Travelling Salesman Problem angewandt wurde. Einer der Unterschiede zum Ant System ist es, dass die Ameisen alle die gleiche Menge an Pheromonen auf ihrer Tour hinterlassen und au\u00dferdem die beste Ameise ausgew\u00e4hlt wird, welche eine zus\u00e4tzliche Menge an Pheromonen auf ihrem Hamiltonkreis abgibt.\r\n\r\nIm Laufe der Jahre wurden weitere Modifikationen entwickelt, die sich auch auf andere Optimierungsprobleme anwenden lassen.<\/p>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n<figure class=\"wp-block-gallery columns-2 wp-block-gallery-3 is-layout-flex wp-block-gallery-is-layout-flex\">\r\n<ul class=\"blocks-gallery-grid\">\r\n \t<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild1ThesisSnapshot.png\"><img decoding=\"async\" loading=\"lazy\" class=\"size-full wp-image-5910\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild1ThesisSnapshot-300x200.png\" alt=\"\" width=\"300\" height=\"200\" data-id=\"5910\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5910\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild1ThesisSnapshot-300x200.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild1ThesisSnapshot-450x300.png 450w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild1ThesisSnapshot.png 600w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/figure>\r\n<\/li>\r\n \t<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild2ThesisSnapshot.png\"><img decoding=\"async\" loading=\"lazy\" class=\"size-full wp-image-5913\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild2ThesisSnapshot-300x200.png\" alt=\"\" width=\"300\" height=\"200\" data-id=\"5913\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5913\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild2ThesisSnapshot-300x200.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild2ThesisSnapshot-450x300.png 450w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/02\/Bild2ThesisSnapshot.png 600w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/figure>\r\n<\/li>\r\n \t<li><\/li>\r\n<\/ul>\r\n<\/figure>\r\n\r\n\r\n\r\n<div class=\"wp-block-columns is-layout-flex wp-container-6 wp-block-columns-is-layout-flex\">\r\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis: 100%;\">\r\n<p class=\"wp-block-preformatted\">Abbildung: Beispiel f\u00fcr die L\u00f6sung des Travelling Salesman Problems mit dem ACS\r\n\r\n\r\n[1] M. Dorigo und T. Stuetzle, Ant Colony Optimization. MIT Press Ltd, 2004, isbn: 9780262042192<\/p>\r\n<\/div>\r\n<\/div>\r\n\r\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/01\/31\/new-open-position-postdoc-in-scientific-computing-2\/\">New Open Position: Postdoc in Scientific Computing<\/a> <span class=\"date\">January 31, 2023<\/span><div class=\"content\"><p>The Numerical Simulation group of Professor Gassner invites applications for a 2.5 year (possibility for extensions afterwards available) postdoc position in scientific computing (pay grade 100% TVL-13).<\/p>\n<p>Deadline for applications is 10.2.2023.<\/p>\n<p>For details, please have a look at the job ad: <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/01\/postdoc_uoc_hpc.pdf\">click here<\/a><\/p>\n<p>For questions, please contact Professor Gassner via Email (<a href=\"mailto:ggassner@uni-koeln.de\">ggassner@uni-koeln.de<\/a>).<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2023\/01\/30\/snapshot-prediction-of-annual-earth-surface-temperature-with-fourier-neural-operator\/\">Snapshot: Prediction of annual earth surface temperature with Fourier Neural Operator<\/a> <span class=\"date\">January 30, 2023<\/span><div class=\"content\"><p>Recently, a lot of research has been done on the development of Neural Operators. The aim is to map infinite-dimensional functions to infinite-dimensional functions with neural networks. One approach is the Fourier Neural Operator [1]. We use this model to estimate the annual temperature based on the global land-sea-ice distribution. The input of the neural network is the earth&#8217;s geography, the output is the global temperature over one year. The training and testing data set is generated with the simulation framework \u201cKlimakoffer\u201d (for more information see <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/09\/30\/snapshot-numerical-simulations-of-earths-climate\/\">Snapshot: Numerical simulations of earth\u2019s climate<\/a> and [2]). <\/p>\n<p>The best results are achieved when using 4 Fourier integral operator layers, 16 Fourier modes, and a network width of 20. The model thus has 52 433 557 trainable parameters. The remaining specifications are as chosen by Li et al. [1]. It took 4 hours and 26 minutes to train the model on an NVIDIA GeForce RTX 3090 GPU and 0.77 seconds to evaluate it for one given example on a CPU. A relative L2 loss of approximately 0.0025 was achieved during training. <\/p>\n<p>The first video shows one exemplary output. The second animation visualizes the difference between the temperature calculated by the FNO model and the \u201ctrue\u201d temperature determined by the framework \u201cKlimakoffer\u201d. <\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/01\/result_temp_no_change.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/01\/result_temp_no_change.gif\" alt=\"\" width=\"640\" height=\"480\" class=\"aligncenter size-full wp-image-5895\" \/><\/a><\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/01\/diff_temp_no_change.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2023\/01\/diff_temp_no_change.gif\" alt=\"\" width=\"640\" height=\"480\" class=\"aligncenter size-full wp-image-5896\" \/><\/a><\/p>\n<p>References:<br \/>\n[1] Li, Z., Kovachki, N., Azizzadenesheli, K., Liu, B., Bhattacharya, K., Stuart, A., Anandkumar, A. (2020). \u201cFourier Neural Operator for Parametric Partial Differential Equations\u201d, arXiv, <a href=\"https:\/\/arxiv.org\/abs\/2010.08895\">https:\/\/arxiv.org\/abs\/2010.08895<\/a><br \/>\n[2] Zhuang, K., North, G.R., Stevens, Mark J. (2017). \u201cA NetCDF version of the two-dimensional energy balance model based on the full multigrid algorithm\u201d, SoftwareX<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/12\/22\/thesis-snapshot-bottom-topography-modelling-for-the-shallow-water-equations-in-trixi-jl-masters-thesis-by-maximilian-bertrand\/\">Thesis Snapshot: Bottom topography modelling for the shallow water equations in Trixi.jl (Master\u2019s thesis by Maximilian Bertrand)<\/a> <span class=\"date\">December 22, 2022<\/span><div class=\"content\"><p>Back in the summer of 2021, severe floods struck Western Germany, resulting in several cities and villages being destroyed and hundreds of people losing their lives. Modelling and predicting those floods can be crucial for investigating such natural disasters. The shallow water equations are a system of physical conservation equations that can model water flow over a given domain. The two dimensional shallow water equations are given as follows:<\/p>\n<p>\\begin{equation*}<br \/>\n    \\left\\{\\begin{aligned}<br \/>\n         \\frac{\\partial}{\\partial t} h + \\left( h v_1\\right)_x + \\left( h v_2\\right)_y =0&#038;\\\\<br \/>\n        \\frac{\\partial}{\\partial t} hv_1 + \\left( hv_1^2 + \\frac{1}{2} gh^2\\right)_x + \\left(h v_1v_2\\right)_y &#038;= -ghb_x\\\\<br \/>\n        \\frac{\\partial}{\\partial t}h v_2 + \\left( hv_1v_2\\right)_x + \\left(hv_2^2 +\\frac{1}{2} gh^2\\right)_y &#038;= -ghb_y<br \/>\n    \\end{aligned}\\right.<br \/>\n\\end{equation*}<\/p>\n<p>Here the variable $b$ describes the underlying bottom topography function of the domain. As this system is set up of hyperbolic partial differential equations, we use the numerical solver <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">Trixi.jl<\/a> to approximate the solutions. This solver is part of the <a href=\"https:\/\/github.com\/trixi-framework\">trixi-framework<\/a>. Trixi.jl already implemented the shallow water equations but did not yet offer the functionality to include real-world topography data into the calculations. Therefore the add-on <a href=\"https:\/\/github.com\/trixi-framework\/TrixiBottomTopography.jl\">TrixiBottomTopography.jl<\/a> has been included in the trixi-framework. This tool uses B-spline interpolation approaches to remodel the bottom topography function  from real-world data, which is often given as elevation data at specific points. One example is the <a href=\"https:\/\/www.bezreg-koeln.nrw.de\/brk_internet\/geobasis\/hoehenmodelle\/digitale_gelaendemodelle\/gelaendemodell\/index.html\">DGM1<\/a> data set which provides uniformly spaced elevation data of the whole German state of North Rhine-Westphalia. Using this data set and the functionalities implemented in Trixi.jl as well as TrixiBottomTopography.jl, we can simulate a dam break problem over a section of the Rhine river valley.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break.gif\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break-1024x552.gif\" alt=\"\" width=\"584\" height=\"315\" class=\"aligncenter size-large wp-image-5880\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break-1024x552.gif 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break-300x162.gif 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break-768x414.gif 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/Dam_break-500x270.gif 500w\" sizes=\"(max-width: 584px) 100vw, 584px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/12\/13\/talk-gregor-gassner-uoc-with-robust-split-form-dgsem-for-hydro-and-magnetohydrodynamics\/\">Talk: Gregor Gassner (UoC) with Robust Split-Form DGSEM for Hydro- and Magnetohydrodynamics<\/a> <span class=\"date\">December 13, 2022<\/span><div class=\"content\"><p>In this talk, we present the class of split form discontinuous Galerkin methods. The notion of split form refers to different interpretations of the non-linear terms in the fluid dynamics equations. For instance the advective part of the momentum flux can be cast in analytically equivalent forms such as the advective form, or the conservative form, or convex combinations of both forms. While these forms are analytically equivalent for smooth solutions, it is interesting to understand that their discrete forms might have strongly different properties. It turns out that specific underlying split forms of the fluid equations give discontinuous Galerkin (DG) approximations with special favourable properties such as, kinetic energy preservation, energy consistency, pressure equilibrium preservation and even entropy conservation\/stability. The most important improvement that can be observed is drastically increased non-linear robustness of the DG discretization, in particular when simulating under-resolved turbulence. A necessary ingredient to retain fully discrete conservation when using split formulations is the so-called summation-by-parts (SBP) property of the discrete derivative and integral operator. It turns out hat specific DG discretizations, such as the Legendre-Gauss-Lobatto (LGL) spectral element method, satisfy the SBP property. We will further discuss in this talk that it is possible to construct compatible low-order finite volume discretizations on the LGL subcell grid, such that a convex blending of the high-order DG method with the robust low-order method is feasible. This allows to construct provably positive hybrid discretizations that enable simulations of problems with strong shock waves, such as e.g. an astrophysical jet with Ma=2000.<\/p>\n<p>Link to the Talk: <a href=\"https:\/\/cassyni.com\/events\/VrFL7HmR9YjHyFJvFcZNjP\" rel=\"noopener\" target=\"_blank\">https:\/\/cassyni.com\/events\/VrFL7HmR9YjHyFJvFcZNjP<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/12\/10\/snapshot-dam-break-against-a-flexible-structure-with-smoothed-particle-hydrodynamics\/\">Snapshot: Dam break against a flexible structure with Smoothed Particle Hydrodynamics<\/a> <span class=\"date\">December 10, 2022<\/span><div class=\"content\"><p>To simulate elastic structures, we use a total Lagrangian SPH formulation [1].<br \/>\nWe simulate a dam break against a flexible structure based on the setup given by [1].<br \/>\nFor the fluid dynamics, a weakly-compressible Lagrangian SPH formulation with an artificial viscosity by Monaghan [2] and boundary force particles to model the tank [3] was used.<br \/>\nThis animation shows the simulation in 0.1x real time with 8.7k fluid particles and 455 particles to model the flexible structure.<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-full wp-image-5855\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/12\/fsi_0.1x.gif\" alt=\"\" width=\"1100\" height=\"646\" \/><\/p>\n<p>References:<br \/>\n[1] O\u2019Connor, J., Rogers, B. D. (2021). \u201cA fluid\u2013structure interaction model for free-surface flows and flexible structures using smoothed particle hydrodynamics on a GPU.\u201d In: <em>Journal of Fluids and Structures<\/em>, 104. <a href=\"https:\/\/doi.org\/10.1016\/J.JFLUIDSTRUCTS.2021.103312\">https:\/\/doi.org\/10.1016\/J.JFLUIDSTRUCTS.2021.103312<\/a><br \/>\n[2] Monaghan, J. J. (1989). \u201cOn the Problem of Penetration in Particle Methods\u201d. In: <em>Journal of Computational Physics<\/em>, 82(1), 1\u201315. <a href=\"https:\/\/doi.org\/10.1016\/0021-9991(89)90032-6\">https:\/\/doi.org\/10.1016\/0021-9991(89)90032-6<\/a><br \/>\n[3] Monaghan, J. J., &amp; Kajtar, J. B. (2009). \u201cSPH particle boundary forces for arbitrary boundaries\u201d. In: <em>Computer Physics Communications<\/em>, 180(10), 1811\u20131820. <a href=\"https:\/\/doi.org\/10.1016\/j.cpc.2009.05.008\">https:\/\/doi.org\/10.1016\/j.cpc.2009.05.008<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/11\/11\/workshop-novel-adaptive-discontinuous-galerkin-approaches-for-the-simulation-of-atmospheric-flows\/\">Workshop: Novel Adaptive Discontinuous Galerkin Approaches for the Simulation of Atmospheric Flows<\/a> <span class=\"date\">November 11, 2022<\/span><div class=\"content\"><p>On November 30th and December 1st, the workshop on &#8220;Novel Adaptive Discontinuous Galerkin Approaches for the Simulation of Atmospheric Flows&#8221; will take place at the Department of Mathematics and Computer Science with researchers from the German Aerospace Center (DLR) Cologne, Deutscher Wetterdienst (DWD), and University of Cologne (UoC). This workshop is supported by and embedded into the <a href=\"https:\/\/cesoc.net\/\">Center for Earth System Observation and Computational Analysis (CESOC)<\/a>.<\/p>\n<p>The idea of the workshop is to discuss novel developments for discontinuous Galerkin methods applied to atmospheric flow simulations, with a focus on implicit time integration, meshing of the sphere, adaptivity, efficient implementations, benchmark test cases for atmospheric flows, entropy-stable split formulations of the nonlinear partial differential equations, and bounds-preserving numerical schemes.<\/p>\n<p>A detailed program with abstracts is available <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/11\/20221130-Workshop-Program-3.pdf\">here<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/11\/08\/snapshot-low-pressure-kelvin-helmholtz-instability-with-dg-and-subcell-positivity-limiter\/\">Snapshot: High Mach number Kelvin-Helmholtz instability with DG and subcell positivity limiter<\/a> <span class=\"date\">November 8, 2022<\/span><div class=\"content\"><p>We solve the compressible Euler equations of gas dynamics with a fourth-order accurate entropy-stable discontinuous Galerkin (DG) method and combine it with a first-order accurate finite volume (FV) method at the node level to impose positivity of density and pressure [1,2].<br \/>\nThe initial condition of this problem is given by:<br \/>\n\\[<br \/>\n\\begin{array}{rlrl}<br \/>\n\\rho (t=0) &#038;= \\frac{1}{2}<br \/>\n+ \\frac{3}{4} B,<br \/>\n&#038;<br \/>\np (t=0) &#038;= 0.1, \\\\<br \/>\nv_1 (t=0) &#038;= \\frac{1}{2} \\left( B-1 \\right),<br \/>\n&#038;<br \/>\nv_2 (t=0) &#038;= \\frac{1}{10} \\sin(2 \\pi x),<br \/>\n\\end{array}<br \/>\n\\]<br \/>\nwith $B=\\tanh \\left( 15 y + 7.5 \\right) &#8211; \\tanh(15y-7.5)$, where $\\rho$ is the density, $\\vec{v}=(v_1,v_2)$ is the velocity, and $p$ is the pressure.<\/p>\n<p><iframe loading=\"lazy\" title=\"Low-pressure Kelvin-Helmholtz instability with DG and subcell positivity limiter\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/sbgMK9YkusI?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<p>The video shows the evolution of the density in time using $1024^2$ and degrees of freedom. We use the entropy-conservative and kinetic energy preserving flux of Ranocha [3] for the volume fluxes and the Rusanov solver for the surface fluxes of the DG and FV methods. During this simulation, the FV method acts on average in 0.000086% of the computational domain, and never more than 0.002% of the computational domain at a specific time.<\/p>\n<p>References<br \/>\n[1] A. M. Rueda-Ram\u00edrez, G. J. Gassner, A Subcell Finite Volume Positivity-Preserving Limiter for DGSEM Discretizations of the Euler Equations, WCCM-ECCOMAS2020, pp. 1\u201312.<br \/>\n[2] A. M. Rueda-Ram\u00edrez, W. Pazner, G. J. Gassner, Subcell limiting strategies for discontinuous galerkin spectral element methods, Computers &#038; Fluids 247 (2022) 105627.<br \/>\n[3] H. Ranocha, Generalised summation-by-parts operators and entropy stability of numerical methods for hyperbolic balance laws, Cuvillier Verlag, 2018.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/10\/04\/snapshot-multi-component-magneto-hydrodynamics-simulations-with-a-bounds-preserving-discontinuous-galerkin-method\/\">Snapshot: Multi-component magneto-hydrodynamics simulations with a bounds-preserving discontinuous Galerkin method<\/a> <span class=\"date\">October 4, 2022<\/span><div class=\"content\"><p>We solve the multi-component magneto-hydrodynamics (MHD) equations using a hybrid finite volume (FV)\/discontinuous Galerkin (DG) method [1], which combines the two discretization approaches at the node level. To test the robustness of the hybrid FV\/DG method, we compute a modification of the Orszag-Tang vortex problem [2], where we use two ion species and initialize the flow with a sharp interface at $y=0$:<br \/>\n\\[<br \/>\n\\begin{array}{rlrl}<br \/>\n\\rho^1(x,y,t=0) &amp;= \\begin{cases}<br \/>\n{25}\/{36 \\pi} &amp; y \\ge 0 \\\\<br \/>\n10^{-8} &amp; y &lt; 0<br \/>\n\\end{cases} ,<br \/>\n&amp;\\rho^2(x,y,t=0) &amp;= \\begin{cases}<br \/>\n10^{-8} &amp; y &lt; 0\\\\<br \/>\n{25}\/{36 \\pi} &amp; y \\ge 0<br \/>\n\\end{cases} , \\\\<br \/>\nv_1 (x,y,t=0) &amp;= &#8211; \\sin (2 \\pi y),<br \/>\n&amp;v_2 (x,y,t=0) &amp;= \\sin (2 \\pi x), \\\\<br \/>\nB_1 (x,y,t=0) &amp;= -\\frac{1}{\\sqrt{4 \\pi}} \\sin (2 \\pi y),<br \/>\n&amp;B_2 (x,y,t=0) &amp;= -\\frac{1}{\\sqrt{4 \\pi}} \\sin (4 \\pi x). \\\\<br \/>\np (x,y,t=0) &amp;= \\frac{5}{12 \\pi},<br \/>\n&amp; &amp; \\\\<br \/>\n\\end{array}<br \/>\n\\]<br \/>\nwhere $\\rho^1$ is the density of the first ion species, $\\rho^2$ is the density of the second ion species,$\\vec{v}=(v_1,v_2)$ is the plasma velocity, $\\vec{B}=(B_1,B_2)$ is the magnetic field, and $p$ is the plasma pressure (without magnetic pressure).<br \/>\n&nbsp;<br \/>\nThe fourth-order DG method delivers high accuracy, but fails to describe shocks and material interfaces. Therefore, we use a modal indicator [3] to detect the elements where the solution changes abruptly. In those elements, we combine the DG method with a robust lower order FV scheme locally (at the node level) to impose TVD-like conditions on the density of the ion species,<br \/>\n\\[<br \/>\n\\min_{k \\in \\mathcal{N} (j)} {\\rho}^{i,\\mathrm{FV}}_{k}<br \/>\n\\le \\rho^i_{j} \\le<br \/>\n\\max_{k \\in \\mathcal{N} (j)} {\\rho}^{i,\\mathrm{FV}}_{k},<br \/>\n\\]<br \/>\nand a local minimum principle on the specific entropy,<br \/>\n\\[<br \/>\ns({\\mathbf{u}}_{j}) \\ge \\min_{k \\in \\mathcal{N} (j)} s(\\mathbf{u}^{\\mathrm{FV}}_{k}),.<br \/>\n\\]<br \/>\nwhere $\\mathcal{N} (j)$ denotes the collection of nodes in the low-order stencil of each node $j$.<br \/>\n&nbsp;<br \/>\nThe video shows the total density ($\\rho=\\rho^1+\\rho^2$), the density of the first ion species ($\\rho^1$), and the so-called blending coefficient ($\\alpha$) for the combination of a fourth-order accurate entropy stable DG method with first- and second-order accurate FV methods (Rusanov solver) and $1024^2$ degrees of freedom. A blending coefficient $\\alpha=0$ means that only the DG is being used, and a blending coefficient $\\alpha=1$ means that only the FV method is being used. The example shows that the hybrid FV\/DG scheme is able to capture the shocks of the simulation correctly, deal with near-vacuum conditions and sharp interfaces.<\/p>\n<p>[1] Rueda-Ram\u00edrez, A. M.; Pazner, W., &amp; Gassner, G. J. (2022). Subcell limiting strategies for discontinuous Galerkin spectral element methods. Computers &#038; Fluids. <a href=\"https:\/\/arxiv.org\/abs\/2202.00576v2\">arXiv:2202.00576<\/a>.<br \/>\n[2] Orszag, S.; Tang, C (1979). Small-scale structure of two-dimensional magnetohydrodynamic turbulence. Journal of Fluid Mechanics.<br \/>\n[3] Persson, P. O.; Peraire, J. (2006). Sub-cell shock capturing for discontinuous Galerkin methods. In 44th AIAA Aerospace Sciences Meeting and Exhibit (p. 112).<\/p>\n<p><iframe loading=\"lazy\" title=\"Multi-component Orszag-Tang vortex simulation (bounds-preserving discontinuous Galerkin method)\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/JeLKOtWXSiI?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/08\/29\/dam-break-simulation-with-smoothed-particle-hydrodynamics\/\">Snapshot: Dam break simulation with Smoothed Particle Hydrodynamics<\/a> <span class=\"date\">August 29, 2022<\/span><div class=\"content\"><p>To study the Smoothed Particle Hydrodynamics method, we simulate a dam break based on the setup given by [1].<br \/>\nWe use a weakly-compressible Lagrangian SPH formulation with an artificial viscosity by Monaghan [2] and boundary force particles to model the tank [3].<\/p>\n<p>This animation shows the simulation in 0.5x real time with 125k fluid particles and 7k boundary particles.<\/p>\n<p><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-full wp-image-5735\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/08\/out_half_speed.gif\" alt=\"\" width=\"1280\" height=\"778\" \/><\/p>\n<p>References:<br \/>\n[1] Marrone, S., Antuono, M., Colagrossi, A., Colicchio, G., le Touz\u00e9, D., &amp; Graziani, G. (2011). \u201c\u03b4-SPH model for simulating violent impact flows\u201d. In: <em>Computer Methods in Applied Mechanics and Engineering<\/em>, 200(13\u201316), 1526\u20131542. <a href=\"https:\/\/doi.org\/10.1016\/J.CMA.2010.12.016\">https:\/\/doi.org\/10.1016\/J.CMA.2010.12.016<\/a><br \/>\n[2] Monaghan, J. J. (1989). \u201cOn the Problem of Penetration in Particle Methods\u201d. In: <em>Journal of Computational Physics<\/em>, 82(1), 1\u201315. <a href=\"https:\/\/doi.org\/10.1016\/0021-9991(89)90032-6\">https:\/\/doi.org\/10.1016\/0021-9991(89)90032-6<\/a><br \/>\n[3] Monaghan, J. J., &amp; Kajtar, J. B. (2009). \u201cSPH particle boundary forces for arbitrary boundaries\u201d. In: <em>Computer Physics Communications<\/em>, 180(10), 1811\u20131820. <a href=\"https:\/\/doi.org\/10.1016\/j.cpc.2009.05.008\">https:\/\/doi.org\/10.1016\/j.cpc.2009.05.008<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/08\/23\/talk-deniz-bezgin-tum-and-aaron-buhendwa-tum-about-jax-fluids-friday-26th-august-2022-10am-cest\/\">Talk: Deniz Bezgin (TUM) and Aaron Buhendwa (TUM) about Differentiable Fluid Dynamics in JAX: Challenges and Perspectives, Friday, 26th August 2022, 10am CEST<\/a> <span class=\"date\">August 23, 2022<\/span><div class=\"content\"><p>JAX-FLUIDS is a CFD solver written in Python, which uses the JAX framework to enable automatic differentiation (AD). This allows one to easily create applications for data-driven simulations or other optimization problems.The talk is based on the recent preprint &#8220;JAX-FLUIDS: A fully-differentiable high-order computational fluid dynamics solver for compressible two-phase flows&#8221; (<a href=\"https:\/\/arxiv.org\/abs\/2203.13760\">arXiv:2203.13760<\/a>).<\/p>\n<p>To obtain the Zoom link for this online talk, please get in touch with <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/gregor-gassner\/\">Gregor Gassner<\/a> or <a href=\"https:\/\/www.hlrs.de\/people\/michael-schlottke-lakemper\">Michael Schlottke-Lakemper<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/06\/10\/snapshot-chaotic-behavior-of-the-kelvin-helmholtz-instability\/\">Snapshot: Chaotic behavior of the Kelvin-Helmholtz Instability<\/a> <span class=\"date\">June 10, 2022<\/span><div class=\"content\"><p>We investigate the Kelvin-Helmholtz instability setup of Fjordholm et al. [1], where the sharp interface in the initial condition is randomly perturbed. The simulation is done with an entropy stable DGSEM with polynomial degree N=3 and 32 x 32 grid cells. The final time is t = 2.0. Trixi.jl [2] is used with the setup found in examples\/tree_2d_dgsem\/elixir_euler_kelvin_helmholtz_instability_fjordholm_etal.jl.<\/p>\n<p>The following gif is showing 100 simulation results of the density distribution for different random interface perturbations:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/06\/10\/snapshot-chaotic-behavior-of-the-kelvin-helmholtz-instability\/density-2\/\" rel=\"attachment wp-att-5710\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/06\/density.gif\" alt=\"\" width=\"600\" height=\"400\" class=\"aligncenter size-full wp-image-5710\" \/><\/a><\/p>\n<p>Here is the average density distribution of all 100 results:<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/06\/10\/snapshot-chaotic-behavior-of-the-kelvin-helmholtz-instability\/average\/\" rel=\"attachment wp-att-5711\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/06\/Average.png\" alt=\"\" width=\"600\" height=\"400\" class=\"aligncenter size-full wp-image-5711\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/06\/Average.png 600w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/06\/Average-300x200.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/06\/Average-450x300.png 450w\" sizes=\"(max-width: 600px) 100vw, 600px\" \/><\/a><\/p>\n<p>[1] Ulrik S. Fjordholm, Roger K\u00e4ppeli, Siddhartha Mishra, Eitan Tadmor. Construction of approximate entropy measure valued solutions for hyperbolic systems of conservation laws, 2014. (<a href=\"https:\/\/arxiv.org\/abs\/1402.0909\">https:\/\/arxiv.org\/abs\/1402.0909<\/a>)<br \/>\n[2] Trixi.jl, a numerical simulation framework for hyperbolic conservation laws written in Julia. (<a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">https:\/\/github.com\/trixi-framework\/Trixi.jl<\/a>)<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/05\/10\/thesis-snapshot-planning-of-collision-free-trajectories-for-uavs-and-optimization-using-uniform-b-splines-bachelor-thesis-by-alexander-bast\/\">Thesis Snapshot: Planning of collision-free trajectories for UAVs and optimization using uniform B-splines. (Bachelor thesis by Alexander Bast)<\/a> <span class=\"date\">May 10, 2022<\/span><div class=\"content\"><p>Autonomous vehicles will be an integral part of our environment and will have a significant impact on transportation and other industries. In this context, drones, so-called UAVs will be used for aerial transportation. Due to the high information density, motion planning is a fundamental problem in the autonomous movement of vehicles.<br \/>\nIn the future, drones should transport goods or perhaps people in a fast, safe and thus collision-free manner. In doing so, the drone should be able to perform all calculations quickly, accurately and autonomously. Therefore, it makes sense to develop algorithms that meet these requirements.<br \/>\nAs part of a bachelor thesis, a method [1] was investigated that uses uniform B-splines to replan a drone&#8217;s trajectory in the presence of a potential collision. This method formulates the task of replanning the trajectory as the solution to a non-linear optimization problem.<br \/>\n\\begin{align*}<br \/>\nE=\\lambda_{ep}(p(t_{ep})-p_{ep})^{2}+\\lambda_{ep}(p'(t_{ep})-p&#8217;_{ep})^{2} + \\lambda_{c} \\int_{t_{min}}^{t_{max}}c(p(t))||p'(t)||dt + \\sum_{i=2}^{4}\\int_{t_{min}}^{t_{max}}\\lambda_{q}(p^{(i)}(t))^{2}dt<br \/>\n\\end{align*}<br \/>\nWith this simplified equation, three cost terms are considered: Endpoint error, collision cost and quadratic derivative cost. With the goal of computing an efficient and collision-free trajectory, an iterative process is used to construct a B-spline curve such that it follows a pre-planned and optimal trajectory. The focus of the analysis was on the selection of the control parameters \u03bbep, \u03bbc, \u03bbq, as well as the number of control points to be optimized and the resulting properties of the optimized trajectory. Using a uniform knot vector, the B-spline is expanded to include control points whose spatial coordinates are determined using nonlinear optimization so that the objective function is minimized. Due to the properties of B-splines, upcoming segments can be locally optimized so that a collision is prevented. The use of a uniform knot vector allows an efficient replanning of the trajectory in real time, which can be performed autonomously by a drone. The following figure shows one result.<\/p>\n<figure class=\"wp-block-gallery aligncenter columns-1 is-cropped\">\n<ul class=\"blocks-gallery-grid\">\n<li class=\"blocks-gallery-item\">\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/SAMPLE.gif\"><img decoding=\"async\" class=\"wp-image-5368 aligncenter\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/SAMPLE.gif\" alt=\"\" data-id=\"5368\" data-full-url=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/SAMPLE.gif\"><\/a><\/figure>\n<\/li>\n<\/ul>\n<\/figure>\n<p>[1] Vladyslav Usenko, Lukas von Stumberg, Andrej Pangercic, Daniel Cremers, Real-Time Trajectory Replanning for MAVs using Uniform B-splines and a 3D Circular Buffer, Technical University of Munich, 2017, https:\/\/arxiv.org\/abs\/1703.01416.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/05\/03\/thesis-snapshot-the-effect-of-sea-ice-seasonality-on-the-two-dimensional-energy-balance-model-klimakoffer-bachelor-thesis-by-xenia-ibach\/\">Thesis Snapshot: The effect of sea ice seasonality on the two-dimensional energy balance model &#8216;Klimakoffer&#8217; (Bachelor thesis by Xenia Ibach)<\/a> <span class=\"date\">May 3, 2022<\/span><div class=\"content\"><p>In a previous snapshot, the two-dimensional energy balance model has been already presented:\\begin{align*}<br \/>\nC(\\hat{r})\\frac{\\partial T(\\hat{r},t)}{\\partial t} +A+BT(\\hat{r},t)= \\nabla\\cdot(D(\\hat{r})\\nabla T(\\hat{r},t)) + QS(\\hat{r},t)(1-a(\\hat{r}))<br \/>\n\\end{align*}<br \/>\nThe model is based on the assumption of a static surface distribution of the surface types land, ocean, sea ice and snow. However, climate change at an increasing rate is also affecting surface types seasonally as well as permanently over the years. Melting ice and snow masses worldwide are a consequence of this, affecting not only snow and sea ice distribution, but also that of oceans and land surfaces through rising sea ice levels.<br \/>\nTherefore, for more accurate climate modeling with this EBM, it is interesting to modify the model with respect to a time-dependent surface distribution. This modification directly affects the model parameters albedo, solar forcing, and heat capacity, requiring the following reformulation of the model:<br \/>\n\\begin{align*}<br \/>\nC(\\hat{r},t)\\frac{\\partial T(\\hat{r},t)}{\\partial t} +A+BT(\\hat{r},t)= \\nabla\\cdot(D(\\hat{r})\\nabla T(\\hat{r},t)) + QS(\\hat{r},t)(1-a(\\hat{r},t))<br \/>\n\\end{align*}<br \/>\nAs part of a bachelor&#8217;s thesis, two options were presented for introducing a seasonal sea ice distribution into the EBM by using and modifying simulations from the &#8216;<a href=\"https:\/\/github.com\/klimakoffer\/Klimakoffer.jl\">Klimakoffer<\/a>&#8216;. Similar approaches are also conceivable for the other surface types.<br \/>\nOne possibility is to process documented measurements of sea ice distribution, available for example from institutions such as the NSIDC (National Snow and Sea Ice Data Center) or NASA, so that a monthly change in sea ice distribution at the edge of the sea ice extent at the two poles can be mapped using a specially developed algorithm.<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_march.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_march.png\" alt=\"\" width=\"451\" height=\"282\" class=\"aligncentral size-full wp-image-5666\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_march.png 451w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_march-300x188.png 300w\" sizes=\"(max-width: 451px) 100vw, 451px\" \/><\/a><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_september.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_september.png\" alt=\"\" width=\"456\" height=\"285\" class=\"aligncentral size-full wp-image-5665\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_september.png 456w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_september-300x188.png 300w\" sizes=\"(max-width: 456px) 100vw, 456px\" \/><\/a><br \/>\nA second option is to directly read in satellite images and categorize the color values of the image&#8217;s pixels and map them to the surface types required in the EBM. Here, as an example, a satellite image of the Blue Marble series from 2004 of NASA was read in and categorized on the basis of the color values to the different surface types.<br \/>\n<div id=\"attachment_5670\" style=\"width: 724px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_satellite.png\"><img aria-describedby=\"caption-attachment-5670\" decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_satellite.png\" alt=\"\" width=\"714\" height=\"357\" class=\"size-full wp-image-5670\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_satellite.png 714w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_satellite-300x150.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/marble_satellite-500x250.png 500w\" sizes=\"(max-width: 714px) 100vw, 714px\" \/><\/a><p id=\"caption-attachment-5670\" class=\"wp-caption-text\">Image from NASA&#8217;s Blue Marble Next Generation [3]<\/p><\/div><br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_november.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_november.png\" alt=\"\" width=\"742\" height=\"464\" class=\"aligncenter size-full wp-image-5669\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_november.png 742w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_november-300x188.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/05\/Marble_november-480x300.png 480w\" sizes=\"(max-width: 742px) 100vw, 742px\" \/><\/a><\/p>\n<p>[1] North, G.R.; Mengel, J.G.; Short, D.A. (1983). Simple Energy Balance Model Resolving the Seasons and the Continents\u2019 Application to the Astronomical Theory of the Ice Ages. Journal of Geophysical Research.<br \/>\n[2] Zhuang, K.; North, G.R.; Stevens, Mark J. (2017). A NetCDF version of the two-dimensional energy balance model based on the full multigrid algorithm. SoftwareX.<br \/>\n[3] EOS Project Science Office at NASA Goddard Space Flight Center, Blue Marble Next Generation, 2004, https:\/\/visibleearth.nasa.gov\/collection\/1484\/blue-marble.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/04\/15\/start-of-new-research-project-dfg-funded-research-unit-snubic-structure-preserving-numerical-methods-for-bulk-and-interface-coupling-of-heterogeneous-models-2022-2026\/\">Start of new Research Project: DFG funded Research Unit \u201cSNUBIC: Structure-Preserving Numerical Methods for Bulk- and Interface Coupling of Heterogeneous Models\u201d (2022 &#8211; 2026)<\/a> <span class=\"date\">April 15, 2022<\/span><div class=\"content\"><p>Funded by the German Research Foundation under the grant number DFG-FOR5409 to investigate the modeling and simulation of coupled systems described by partial differential equations (PDEs).<\/p>\n<p>Please visit the research unit webpage for all the details: <a href=\"https:\/\/www.snubic.io\/\">SNUBIC<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/04\/01\/start-of-new-research-project-klaus-tschira-stiftung-funded-project-hifilab-a-high-fidelity-laboratory-for-the-simulation-of-celestial-bodies-with-their-space-environment-2022\/\">Start of new Research Project: Klaus-Tschira-Stiftung funded Project \u201cHiFiLab: A High-Fidelity Laboratory for the Simulation of Celestial Bodies with their Space Environment\u201d (2022 &#8211; 2025)<\/a> <span class=\"date\">April 1, 2022<\/span><div class=\"content\"><p>In this project, we focus on generating a novel computational simulation framework to describe the interaction of plasma with celestial bodies. Understanding the interaction of celestial bodies with their space environment is very important, as it often reveals information about their inner structure and the existence\/composition of their atmospheres. Of fundamental importance is the question about liquid water under the icy surface of some moons of the solar system, as water is considered to be one of the essential ingredients for life as we know it.In the last years, we have successfully designed a high-order accurate 3D unstructured discontinuous Galerkin (DG) open source solver with fully parallel adaptive mesh refinement for single-fluid magnetohydrodynamics. DG methods are famous for their high accuracy, their high flexibility and extreme parallel scaling capabilities and are thus perfectly suited for complex plasma interaction simulations. We plan a major step forward regarding the physical modeling fidelity of our computational plasma framework, by extending our high-order DG solver to multi-ion MHD models that account for the interaction of electrons, ions, and neutrals. We will further apply the resulting novel computational plasma framework to simulate the Jovian moon Europa and compare our results with data taken by space missions during flybys of the moon and observations from the Hubble Space Telescope to gain insight and better understanding of the complex plasma interactions.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/03\/23\/snapshot-a-high-order-detonation-diffraction-simulation-with-reactive-euler-equations\/\">Snapshot: A high order detonation diffraction simulation with reactive Euler equations<\/a> <span class=\"date\">March 23, 2022<\/span><div class=\"content\"><p>We show a detonation diffraction simulation with multiple obstacles in a reflective box using reactive multicomponent Euler equations with Trixi.jl. These kind of detonation waves may result in pressure and density drops close to zero which is numerically difficult not only but especially for high order schemes.<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/03\/23\/snapshot-a-high-order-detonation-diffraction-simulation-with-reactive-euler-equations\/pressure_double_diff_high_res\/\" rel=\"attachment wp-att-5609\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/pressure_double_diff_high_res.gif\" alt=\"\" width=\"2132\" height=\"1352\" class=\"aligncenter size-full wp-image-5609\" \/><\/a><\/p>\n<p>In this simulation we use a reacting model which consists of one reaction and two species. This example has been calculated using the high order DG-FV blending method in Trixi.jl with a HOHQMesh generated mesh whereas the chemical network is solved with KROME.<\/p>\n<p><em>(Note: Not all of these features are currently available in the main code of Trixi.jl.)<\/em><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/03\/22\/new-paper-submitted-on-the-theoretical-foundation-of-overset-grid-methods-for-hyperbolic-problems-ii-entropy-bounded-formulations-for-nonlinear-conservation-laws\/\">New paper submitted: On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems II: Entropy Bounded Formulations for Nonlinear Conservation Laws<\/a> <span class=\"date\">March 22, 2022<\/span><div class=\"content\"><p>We derive entropy conserving and entropy dissipative overlapping domain formulations for systems of nonlinear hyperbolic equations in conservation form, such as would be approximated by overset mesh methods. The entropy conserving formulation imposes two-way coupling at the artificial interface boundaries through nonlinear penalty functions that vanish when the solutions coincide. The penalty functions are expressed in terms of entropy conserving fluxes originally introduced for finite volume schemes. Entropy dissipation and additional coupling in the overlap region are added through the use of linear penalties.<\/p>\n<p>arXiv: <a href=\"https:\/\/arxiv.org\/abs\/2203.11149\">https:\/\/arxiv.org\/abs\/2203.11149<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/03\/22\/new-paper-submitted-on-the-entropy-projection-and-the-robustness-of-high-order-entropy-stable-discontinuous-galerkin-schemes-for-under-resolved-flows\/\">New paper submitted: On the entropy projection and the robustness of high order entropy stable discontinuous Galerkin schemes for under-resolved flows<\/a> <span class=\"date\">March 22, 2022<\/span><div class=\"content\"><p>High order entropy stable schemes provide improved robustness for computational simulations of fluid flows. However, additional stabilization and positivity preserving limiting can still be required for variable-density flows with under-resolved features. We demonstrate numerically that entropy stable DG methods which incorporate an &#8220;entropy projection&#8221; are less likely to require additional limiting to retain positivity for certain types of flows. We conclude by investigating potential explanations for this observed improvement in robustness.<\/p>\n<p>arXiv: <a href=\"https:\/\/arxiv.org\/abs\/2203.10238\">https:\/\/arxiv.org\/abs\/2203.10238<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/03\/14\/new-paper-submitted-entropy-stable-gauss-collocation-methods-for-ideal-magneto-hydrodynamics\/\">New paper submitted: Entropy-Stable Gauss Collocation Methods for Ideal Magneto-Hydrodynamics<\/a> <span class=\"date\">March 14, 2022<\/span><div class=\"content\"><p>In this paper, we present an entropy-stable Gauss collocation discontinuous Galerkin (DG) method on 3D curvilinear meshes for the GLM-MHD equations: the single-fluid magneto-hydrodynamics (MHD) equations with a generalized Lagrange multiplier (GLM) divergence cleaning mechanism. For the continuous entropy analysis to hold and to ensure Galilean invariance in the divergence cleaning technique, the GLM-MHD system requires the use of non-conservative terms.<br \/>\nTraditionally, entropy-stable DG discretizations have used a collocated nodal variant of the DG method, also known as the discontinuous Galerkin spectral element method (DGSEM) on Legendre-Gauss-Lobatto (LGL) points. Recently, Chan et al. (&#8220;Efficient Entropy Stable Gauss Collocation Methods&#8221;. SIAM -2019) presented an entropy-stable DGSEM scheme that uses Legendre-Gauss points (instead of LGL points) for conservation laws. Our main contribution is to extend the discretization technique of Chan et al. to the non-conservative GLM-MHD system.<br \/>\nWe provide a numerical verification of the entropy behavior and convergence properties of our novel scheme on 3D curvilinear meshes. Moreover, we test the robustness and accuracy of our scheme with a magneto-hydrodynamic Kelvin-Helmholtz instability problem. The numerical experiments suggest that the entropy-stable DGSEM on Gauss points for the GLM-MHD system is more accurate than the LGL counterpart. <\/p>\n<p><a href=\"https:\/\/arxiv.org\/abs\/2203.06062\">arXiv:2203.06062<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/02\/25\/snapshot-supersonic-turbulence-simulations-with-fluxo\/\">Snapshot: Supersonic turbulence simulations with FLUXO<\/a> <span class=\"date\">February 25, 2022<\/span><div class=\"content\"><p>We show a supersonic turbulence simulation in a periodic box fueled by perpetual kinetic energy injection on the largest three modes. Such setups are subject of active research in astrophysics and are believed to play an important role in the dynamics of galactic clouds and star formation, e.g., [1,2,3]. An initially uniform density distribution (left plot in the video) is gradually driven to a turbulent state till an average velocity of Mach 3.5 is reached (about t=0.19). Locally, speeds can spike above Mach 35 (bottom right plot in the video). The simulation was computed on FLUXO (<a href=\"https:\/\/github.com\/project-fluxo\/fluxo\">https:\/\/github.com\/project-fluxo\/fluxo<\/a>) with a robust invariant domain preserving Discontinuous Galerkin scheme (DG) [4] solving the compressible, quasi-isothermal ($\\gamma =1.001$) Euler equations on a uniform grid with $128^3$ degrees-of-freedom.<\/p>\n<p>Smooth blending with a monotone finite volume scheme stabilizes the high-order DG at strongly shocked regions (filaments of high density concentration) displayed by the blending parameter plot (top right in the video). Blue encodes regions of 100% high order DG while red dots indicate focused blending with FV. Hypersonic turbulence regimes are very challenging for any numerical scheme considering the stark gradients in density, bubbles of near-vacuum and the extremely high speeds involved in such simulations.<\/p>\n<p>[1] https:\/\/academic.oup.com\/mnras\/article\/436\/2\/1245\/1126116<br \/>\n[2] https:\/\/adsabs.harvard.edu\/full\/1981MNRAS.194..809L]<br \/>\n[3] https:\/\/academic.oup.com\/mnras\/article\/480\/3\/3916\/5060766<br \/>\n[4] Rueda-Ram\u00edrez, A. M., Pazner, W., &amp; Gassner, G. J. (2022). Subcell limiting strategies for discontinuous Galerkin spectral element methods. <a href=\"https:\/\/arxiv.org\/abs\/2202.00576v2\">arXiv preprint arXiv:2202.00576<\/a>.<\/p>\n<p><iframe loading=\"lazy\" title=\"Supersonic turbulent box simulation using a Discontinuous Galerkin \/ Finite Volume Methods\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/Kd5aNS5Cm6k?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/01\/31\/snapshot-hypersonic-astrophysical-jet-mach2000-simulation-using-a-hybrid-dg-fv-solver\/\">Snapshot: Hypersonic astrophysical jet (Mach~2000) simulation using a hybrid DG\/FV solver<\/a> <span class=\"date\">January 31, 2022<\/span><div class=\"content\"><p>Simulation of a hypersonic astrophysical jet with Mach=2156.91 moving through a medium at rest.<\/p>\n<p>The simulation was performed on a 256 \u00d7 256 quadrilateral grid with a nodal fourth-order entropy-stable Discontinuous Galerkin (DG) scheme. Local bounds on the density and specific entropy are imposed by blending the DG scheme locally (at the node level) with a second-order Finite Volume (FV) method. Both the DG and the FV methods use the Harten-Lax-van Leer-Einfeldt (HLLE) Riemann solver.<\/p>\n<p>The simulation results were obtained with the open-source code FLUXO (https:\/\/github.com\/project-fluxo\/fluxo).<\/p>\n<p>The numerical methods used to perform the simulation are described in our preprint:<br \/>\nRueda-Ram\u00edrez, A. M., Pazner, W., Gassner, G.J. (2022). <em>Subcell Limiting Strategies for Discontinuous Galerkin Spectral Element Methods<\/em>. Submitted. <a href=\"https:\/\/arxiv.org\/abs\/2202.00576\">https:\/\/arxiv.org\/abs\/2202.00576<\/a>.<\/p>\n<p><iframe loading=\"lazy\" title=\"Hypersonic astrophysical jet (Mach~2000) simulation using a hybrid DG\/FV solver\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/s4AdclDe1vs?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/01\/31\/new-paper-submitted-subcell-limiting-strategies-for-discontinuous-galerkin-spectral-element-methods\/\">New paper submitted: Subcell limiting strategies for discontinuous Galerkin spectral element methods<\/a> <span class=\"date\">January 31, 2022<\/span><div class=\"content\"><p>We present a general family of subcell limiting strategies to construct robust high-order accurate nodal discontinuous Galerkin (DG) schemes. The main strategy is to construct compatible low order finite volume (FV) type discretizations that allow for convex blending with the high-order variant with the goal of guaranteeing additional properties, such as bounds on physical quantities and\/or guaranteed entropy dissipation. For an implementation of this main strategy, four main ingredients are identified that may be combined in a flexible manner: (i) a nodal high-order DG method on Legendre-Gauss-Lobatto nodes, (ii) a compatible robust subcell FV scheme, (iii) a convex combination strategy for the two schemes, which can be element-wise or subcell-wise, and (iv) a strategy to compute the convex blending factors, which can be either based on heuristic troubled-cell indicators, or using ideas from flux-corrected transport methods.<br \/>\nBy carefully designing the metric terms of the subcell FV method, the resulting methods can be used on unstructured curvilinear meshes, are locally conservative, can handle strong shocks efficiently while directly guaranteeing physical bounds on quantities such as density, pressure or entropy. We further show that it is possible to choose the four ingredients to recover existing methods such as a provably entropy dissipative subcell shock-capturing approach or a sparse invariant domain preserving approach.<br \/>\nWe test the versatility of the presented strategies and mix and match the four ingredients to solve challenging simulation setups, such as the KPP problem (a hyperbolic conservation law with non-convex flux function), turbulent and hypersonic Euler simulations, and MHD problems featuring shocks and turbulence. <\/p>\n<p><a href=\"https:\/\/arxiv.org\/abs\/2202.00576\">arXiv:2202.00576<\/a><\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-300x167.png\" alt=\"\" width=\"300\" height=\"167\" class=\"alignnone size-medium wp-image-5594\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-300x167.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-1024x569.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-768x427.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-1536x854.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-2048x1139.png 2048w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/03\/OT_2nd_LLF_densL-500x278.png 500w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2022\/01\/27\/new-open-position-postdoc-in-scientific-computing\/\">New Open Position: PostDoc in Scientific Computing<\/a> <span class=\"date\">January 27, 2022<\/span><div class=\"content\"><p>The Numerical Simulation group of Professor Gassner invites applications for a 2 year (possibility for extensions afterwards available) postdoc position in scientific computing (pay grade 100% TVL-13).<\/p>\n<p>Deadline for applications is 28.2.2022.<\/p>\n<p>For details, please have a look at the job ad:<\/p>\n<p>German version (<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/01\/german_ad.pdf\" target=\"_blank\" rel=\"noopener\">click here<\/a>)<br \/>\nEnglish version (<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2022\/01\/english_ad.pdf\" target=\"_blank\" rel=\"noopener\">click here<\/a>)<\/p>\n<p>For questions, please directly contact Professor Gassner via Email (<a href=\"mailto:ggassner@uni-koeln.de\">ggassner@uni-koeln.de<\/a>).<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/12\/22\/new-paper-submitted-a-discontinuous-galerkin-solver-in-the-flash-multi-physics-framework\/\">New paper submitted: A Discontinuous Galerkin Solver in the FLASH Multi-Physics Framework<\/a> <span class=\"date\">December 22, 2021<\/span><div class=\"content\"><p>In this paper, we present a discontinuous Galerkin solver based on previous work by Markert et al. (2021) for magneto-hydrodynamics in form of a new fluid solver module integrated into the established and well-known multi-physics simulation code FLASH. Our goal is to enable future research on the capabilities and potential advantages of discontinuous Galerkin methods for complex multi-physics simulations in astrophysical settings. We give specific details and adjustments of our implementation within the FLASH framework and present extensive validations and test cases, specifically its interaction with several other physics modules such as (self-)gravity and radiative transfer. We conclude that the new DG solver module in FLASH is ready for use in astrophysics simulations and thus ready for assessments and investigations.<\/p>\n<p><a href=\"https:\/\/arxiv.org\/abs\/2112.11318\">arXiv:2112.11318<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/12\/22\/new-paper-submitted-efficient-implementation-of-modern-entropy-stable-and-kinetic-energy-preserving-discontinuous-galerkin-methods-for-conservation-laws\/\">New paper submitted: Efficient implementation of modern entropy stable and kinetic energy preserving discontinuous Galerkin methods for conservation laws<\/a> <span class=\"date\">December 22, 2021<\/span><div class=\"content\"><p>Many modern discontinuous Galerkin (DG) methods for conservation laws make use of summation by parts operators and flux differencing to achieve kinetic energy preservation or entropy stability. While these techniques increase the robustness of DG methods significantly, they are also computationally more demanding than standard weak form nodal DG methods. We present several implementation techniques to improve the efficiency of flux differencing DG methods that use tensor product quadrilateral or hexahedral elements, in 2D or 3D respectively. Focus is mostly given to CPUs and DG methods for the compressible Euler equations, although these techniques are generally also useful for GPU computing and other physical systems including the compressible Navier-Stokes and magnetohydrodynamics equations. We present results using two open source codes, Trixi.jl written in Julia and FLUXO written in Fortran, to demonstrate that our proposed implementation techniques are applicable to different code bases and programming languages.<\/p>\n<p><a href=\"https:\/\/arxiv.org\/abs\/2112.10517\">arXiv:2112.10517<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/11\/30\/snapshot-immersed-boundary-discontinuous-galerkin-simulation-of-flow-past-a-cylinder-at-re100-using\/\">Snapshot: Immersed Boundary Discontinuous Galerkin Simulation of Flow Past a Cylinder at Re=100<\/a> <span class=\"date\">November 30, 2021<\/span><div class=\"content\"><p>Simulation of the viscous flow past a cylinder with a discontinuous Galerkin spectral element method on a hierarchical Cartesian mesh with adaptive mesh refinement (AMR). The polynomial degree is N = 3 and the mesh uses five refinement levels. The cylinder boundary is imposed with an immersed boundary method (IBM) based on volume penalization [1], which adds a source term to the compressible Navier-Stokes equations,<br \/>\n$$<br \/>\n\\partial_t \\mathbf{u} + \\nabla \\cdot \\vec{\\mathbf{f}} (\\mathbf{u},\\nabla \\mathbf{u}) = \\underbrace{\\beta (\\mathbf{u}_s &#8211; \\mathbf{u})}_{\\text{IBM source}},<br \/>\n$$<br \/>\nwhere $\\mathbf{u}$ is the vector of conservative quantities, $\\vec{\\mathbf{f}}$ is the Navier-Stokes flux, $\\beta$ is the penalty parameter, and $\\mathbf{u}_s$ is an imposed state. For this simulation, $\\mathbf{u}_s$ mimics an isothermal no-slip wall boundary condition.<\/p>\n<p>The simulation was obtained with the open-source code FLUXO (https:\/\/github.com\/project-fluxo\/fluxo).<\/p>\n<p><strong>References<\/strong><br \/>\n[1] Kou, J., Joshi, S., Hurtado-de-Mendoza, A., Puri, K., Hirsch, C., &#038; Ferrer, E. (2022). Immersed boundary method for high-order flux reconstruction based on volume penalization. Journal of Computational Physics, 448, 110721.<\/p>\n<p><iframe loading=\"lazy\" title=\"Immersed Boundary Discontinuous Galerkin Simulation of Flow Past a Cylinder at Re=100 using\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/F3KnEpcv2UY?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/10\/27\/new-paper-published-on-the-theoretical-foundation-of-overset-grid-methods-for-hyperbolic-problems-well-posedness-and-conservation\/\">New paper published: On the theoretical foundation of overset grid methods for hyperbolic problems: Well-posedness and conservation<\/a> <span class=\"date\">October 27, 2021<\/span><div class=\"content\"><p>We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is usually the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem.<\/p>\n<p>Published in Journal of Computational Physics (<a href=\"https:\/\/www.sciencedirect.com\/science\/article\/pii\/S0021999121006276?via%3Dihub\">ScienceDirect<\/a>)<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/10\/19\/snapshot-simulation-of-a-stiff-multi-species-reacting-flow\/\">Snapshot: Simulation of a stiff multi-species reacting flow<\/a> <span class=\"date\">October 19, 2021<\/span><div class=\"content\"><p>Simulation of an asymmetric two-dimensional stiff multi-species reacting flow for the compressible Euler equations. The computation domain <code>[0,6]x[0,2]<\/code> of this detonation problem is split into three zones. Zone A <code>[0,0.5]x[0,2]<\/code> as well as zone B <code>(0.5,6]x[1.2,2]<\/code> is filled with burnt gas, namely <code>O2,OH<\/code> and <code>H2O<\/code>. The last zone C <code>(0.5,6]x[0,1.2)<\/code> contains unburnt gas, namely <code>H2<\/code> and <code>O2<\/code>. <\/p>\n<p>The reacting model consists of two reactions <code>H2+O2-&gt;2OH<\/code>, <code>2OH+H2-&gt;2H2O<\/code> and five species <code>H2,O2,OH,H2O,N2<\/code> whereas reaction 1 has a smaller ignition temperature as well as a much faster chemical rate than reaction 2. Species <code>N2<\/code> acts as a dilute catalyst here. <\/p>\n<p>This example has been calculated using Trixi.jl with a P4est mesh using <code>3072x1024=3.145.728<\/code> elements, Dirichlet boundary conditions in x-direction as well as solid wall boundaries in y-direction. The final time is T = 0.1. The convection part is solved with FV whereas the chemical network is solved with KROME. <em>(Note: Not all of these features are currently available in the main code of Trixi.jl.)<\/em><\/p>\n<p>Pressure distribution at final time T=0.1:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-300x98.png\" alt=\"\" width=\"600\" height=\"196\" class=\"alignnone size-medium wp-image-5500\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-300x98.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-1024x336.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-768x252.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-1536x503.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure-500x164.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/pressure.png 1770w\" sizes=\"(max-width: 600px) 100vw, 600px\" \/><\/a><\/p>\n<p>Temperature distribution at final time T=0.1:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-300x98.png\" alt=\"\" width=\"600\" height=\"196\" class=\"alignnone size-medium wp-image-5510\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-300x98.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-1024x334.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-768x250.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-1536x501.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature-500x163.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/temperature.png 1766w\" sizes=\"(max-width: 600px) 100vw, 600px\" \/><\/a><\/p>\n<p>Density distribution at final time T=0.1:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/density.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/density-300x97.png\" alt=\"\" width=\"600\" height=\"196\" class=\"alignnone size-medium wp-image-5511\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/density-300x97.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/density-500x162.png 500w\" sizes=\"(max-width: 600px) 100vw, 600px\" \/><\/a><\/p>\n<p>OH-fraction distribution at final time T=0.1:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-300x98.png\" alt=\"\" width=\"600\" height=\"196\" class=\"alignnone size-medium wp-image-5512\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-300x98.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-1024x333.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-768x250.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-1536x500.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction-500x163.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/10\/OH_fraction.png 1782w\" sizes=\"(max-width: 600px) 100vw, 600px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/09\/30\/snapshot-numerical-simulations-of-earths-climate\/\">Snapshot: Numerical simulations of earth&#8217;s climate<\/a> <span class=\"date\">September 30, 2021<\/span><div class=\"content\"><p>Earth&#8217;s climate is changing due to human activity. Therefore, the ability to predict future climate change is of utmost importance if we want to reduce the human interference with the climate system (climate change mitigation) and our vulnerability to the harmful effects of climate change (climate change adaptation). <\/p>\n<p>It is possible to model the time-dependent local temperatures on the globe using a simple two-dimensional energy balance model [1,2],<br \/>\n$$\u00a0C(\\vec{x}) \\frac{\\partial T}{\\partial t} + A + BT = \\nabla \\cdot \\left( D(\\vec{x}) \\nabla T \\right) + Q S(\\vec{x},t) (1-a(\\vec{x})), $$<br \/>\nwhere $T$ is the local temperature; $C(\\vec{x})$ is the local heat capacity for a column of air on land, water or ice; the term $A+BT$ models the outgoing long-wave radiation, where $A$ depends on the amount of green-house gases in the atmosphere and $B$ models several feedback effects; $D$ is a latitude-dependent diffusion coefficient; and the last term includes the heat input in the system, where $Q$ is the solar radiation at earth&#8217;s position in the solar system, $S$ is the seasonal and latitudinal distribution of insolation, which depends on orbital parameters, such as the eccentricity of earth&#8217;s orbit, and the precession and tilt of earth&#8217;s rotation axis, and $a(\\vec{x})$ is the local albedo, which measures the reflection of solar radiation and depends on the surface cover.<\/p>\n<p>We developed a <a href=\"https:\/\/julialang.org\/\">julia<\/a> code that implements the EBM model on the sphere surface and used it to describe the mean temperatures on the globe assuming a constant atmospheric CO$_2$ concentration (of year 1950). The results are presented in the animation.<\/p>\n<figure class=\"wp-block-gallery aligncenter columns-1 is-cropped\">\n<ul class=\"blocks-gallery-grid\">\n<li class=\"blocks-gallery-item\">\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/09\/Earth_temp.gif\"><img decoding=\"async\" class=\"wp-image-5368 aligncenter\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/09\/Earth_temp.gif\" alt=\"\" data-id=\"5368\" data-full-url=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/09\/Earth_temp.gif\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5368\" \/><\/a><\/figure>\n<\/li>\n<\/ul>\n<\/figure>\n<p>[1] North, G.R.; Mengel, J.G.; Short, D.A. (1983). Simple Energy Balance Model Resolving the Seasons and the Continents&#8217; Application to the Astronomical Theory of the Ice Ages. Journal of Geophysical Research.<br \/>\n[2] Zhuang, K.; North, G.R.; Stevens, Mark J. (2017). A NetCDF version of the two-dimensional energy balance model based on the full multigrid algorithm. SoftwareX.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/08\/15\/snapshot-simulation-of-a-colliding-flow-setup-with-dgsem\/\">Snapshot: Simulation of a colliding flow setup with DGSEM<\/a> <span class=\"date\">August 15, 2021<\/span><div class=\"content\"><div>Simulation of a colliding flow setup with DGSEM (polynomial degree 3,<\/div>\n<div>compatible sub-cell LGL type FV for shock capturing) of the compressible<\/div>\n<div>Euler equations, where at x = +\/- 64 inflow boundary conditions with<\/div>\n<div>Mach number Ma = 70 is set (with periodic boundary conditions in y<\/div>\n<div>direction).<\/p>\n<div>Density distribution, where the color bar is plotted in logarithmic scale, at<\/div>\n<div>final simulation time t=25:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF1.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5417\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF1-300x245.png\" alt=\"\" width=\"300\" height=\"245\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF1-300x245.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF1-368x300.png 368w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF1.png 718w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<div>Simulation with AMR (level 3 up to level 9), showing the AMR indicator at<\/div>\n<div>final time t=25:<\/div>\n<\/div>\n<div><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF2.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5418\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF2-300x235.png\" alt=\"\" width=\"300\" height=\"235\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF2-300x235.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF2-383x300.png 383w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF2.png 661w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<div>Plot of the DGSEM\/FV blending factor for shock capturing at final time<\/div>\n<div>t=25:<\/div>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5419\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4-299x300.png 299w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF4.png 514w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<section class=\"pdf_page\" aria-label=\"Page 2\">\n<div class=\"textlayer\">\n<div>The combination of the AMR indicator and the shock capturing blending<\/div>\n<div>factor triggers the final mesh refinement, as can be seen here at the final time t=25:<\/div>\n<\/div>\n<\/section>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF3.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5420\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF3-300x241.png\" alt=\"\" width=\"300\" height=\"241\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF3-300x241.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF3-373x300.png 373w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/08\/CF3.png 666w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/div>\n<\/div>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/08\/09\/juliacon-2021-adaptive-and-extendable-numerical-simulations-with-trixi-jl\/\">JuliaCon 2021: Adaptive and extendable numerical simulations with Trixi.jl<\/a> <span class=\"date\">August 9, 2021<\/span><div class=\"content\">\n<p><a rel=\"noreferrer noopener\" href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" data-type=\"URL\" data-id=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" target=\"_blank\">Trixi.jl<\/a> is a numerical simulation framework for adaptive, high-order discretizations of conservation laws. It has a modular architecture that allows users to easily extend its functionality and was designed to be useful to experienced researchers and new users alike. In this talk, we give an overview of Trixi\u2019s current features, present a typical workflow for creating and running a simulation, and show how to add new capabilities for your own research projects.<\/p>\n\n\n\n<p>This talk was given on July 30th, 2021 by <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\" data-type=\"page\" data-id=\"3741\">Michael Schlottke-Lakemper<\/a> and <a rel=\"noreferrer noopener\" href=\"https:\/\/ranocha.de\" data-type=\"URL\" data-id=\"https:\/\/ranocha.de\" target=\"_blank\">Hendrik Ranocha<\/a> as part of <a rel=\"noreferrer noopener\" href=\"https:\/\/juliacon.org\/2021\/\" data-type=\"URL\" data-id=\"https:\/\/juliacon.org\/2021\/\" target=\"_blank\">JuliaCon 2021<\/a>.<\/p>\n\n\n\n<p>Talk on YouTube: <a rel=\"noreferrer noopener\" href=\"https:\/\/www.youtube.com\/watch?v=hoViWRAhCBE\" target=\"_blank\">https:\/\/www.youtube.com\/watch?v=hoViWRAhCBE<\/a><br>Repository: <a rel=\"noreferrer noopener\" href=\"https:\/\/github.com\/trixi-framework\/talk-2021-juliacon\" target=\"_blank\">https:\/\/github.com\/trixi-framework\/talk-2021-juliacon<\/a><br>Conference agenda entry: <a rel=\"noreferrer noopener\" href=\"https:\/\/live.juliacon.org\/talk\/VAGFD7\" target=\"_blank\">https:\/\/live.juliacon.org\/talk\/VAGFD7<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/08\/03\/new-paper-published-an-entropy-stable-nodal-discontinuous-galerkin-method-for-the-resistive-mhd-equations-part-ii-subcell-finite-volume-shock-capturing\/\">New paper published: An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part II: Subcell finite volume shock capturing<\/a> <span class=\"date\">August 3, 2021<\/span><div class=\"content\"><p id=\"sp0300\">The second paper of this series presents two robust entropy stable shock-capturing methods for discontinuous Galerkin spectral element (DGSEM) discretizations of the compressible magneto-hydrodynamics (MHD) equations. Specifically, we use the resistive GLM-MHD equations, which include a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM). For the continuous entropy analysis to hold, and due to the divergence-free constraint on the magnetic field, the GLM-MHD system requires the use of non-conservative terms, which need special treatment.<\/p>\n<p id=\"sp0310\">Hennemann et al. (2020) <a class=\"workspace-trigger\" href=\"https:\/\/www.sciencedirect.com\/science\/article\/pii\/S0021999121004757?via%3Dihub#br0250\" name=\"bbr0250\">[25]<\/a> recently presented an entropy stable shock-capturing strategy for DGSEM discretizations of the <a class=\"topic-link\" title=\"Learn more about Euler equations from ScienceDirect's AI-generated Topic Pages\" href=\"https:\/\/www.sciencedirect.com\/topics\/computer-science\/euler-equation\">Euler equations<\/a> that blends the DGSEM scheme with a subcell first-order finite volume (FV) method. Our first contribution is the extension of the method of Hennemann et al. to systems with non-conservative terms, such as the GLM-MHD equations. In our approach, the advective and non-conservative terms of the equations are discretized with a hybrid FV\/DGSEM scheme, whereas the visco-resistive terms are discretized only with the high-order DGSEM method. We prove that the extended method is semi-discretely entropy stable on three-dimensional unstructured curvilinear meshes. Our second contribution is the derivation and analysis of a second entropy stable shock-capturing method that provides enhanced resolution by using a subcell reconstruction procedure that is carefully built to ensure entropy stability.<\/p>\n<p>We provide a numerical verification of the properties of the hybrid FV\/DGSEM schemes on curvilinear meshes and show their robustness and accuracy with common benchmark cases, such as the Orszag-Tang <a class=\"topic-link\" title=\"Learn more about vortex from ScienceDirect's AI-generated Topic Pages\" href=\"https:\/\/www.sciencedirect.com\/topics\/physics-and-astronomy\/vortices\">vortex<\/a> and the GEM (Geospace Environmental Modeling) reconnection challenge. Finally, we simulate a space physics application: the interaction of <a class=\"topic-link\" title=\"Learn more about Jupiter's from ScienceDirect's AI-generated Topic Pages\" href=\"https:\/\/www.sciencedirect.com\/topics\/physics-and-astronomy\/jupiter\">Jupiter&#8217;s<\/a> magnetic field with the plasma <a class=\"topic-link\" title=\"Learn more about torus from ScienceDirect's AI-generated Topic Pages\" href=\"https:\/\/www.sciencedirect.com\/topics\/physics-and-astronomy\/toruses\">torus<\/a> generated by the moon Io.<\/p>\n<p>Published in Journal of Computational Physics (<a href=\"https:\/\/www.sciencedirect.com\/science\/article\/pii\/S0021999121004757?via%3Dihub\">ScienceDirect<\/a>)<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/07\/31\/snapshot-interaction-of-jupiters-moon-io-with-the-jovian-magnetosphere\/\">Snapshot: Interaction of Jupiter\u2019s moon Io with the Jovian magnetosphere<\/a> <span class=\"date\">July 31, 2021<\/span><div class=\"content\"><p>Io is the most volcanically active body of the solar system. Furthermore, it is embedded in Jupiter&#8217;s magnetic field, the largest and most powerful planetary magnetosphere of the solar system. Due to its strong volcanic activity, Io expels ions and neutrals, which are in turn ionized by ultraviolet and electron impact ionization, forming a plasma torus around Jupiter [1,2]. As Io moves inside the plasma torus, elastic collisions of ions and neutrals inside its atmosphere generate a magnetospheric disturbance that propagates away from Io along the background magnetic field lines at the Alfv\u00e9n wave speed. This phenomenon creates a pair of Alfv\u00e9n current tubes that are commonly called Alfv\u00e9n wings, which have been observed by several flybys [1].<\/p>\n<p>The figure shows the momentum (non-dimensional) and magnetic field of the plasma that surrounds Io, obtained with the magneto-hydrodynamic (MHD) plasma model of FLUXO*. The Alfv\u00e9n wings can be observed as disturbances of both the background magnetic and momentum fields. In yellow is the trajectory of the I31 flyby of the Galileo spacecraft, which visited Io in 2001 [1]. The Galileo spacecraft is also depicted (not to scale). A 3D model of the moon\u2019s surface developed by NASA [3] was superposed on the simulation results.<\/p>\n<p>*FLUXO (<a href=\"http:\/\/www.github.com\/project-fluxo\/fluxo\">www.github.com\/project-fluxo\/fluxo<\/a>) is an MPI parallel high-order Discontinuous Galerkin code, which supports unstructured curvilinear hexahedral grids, and is able to perform Adaptive Mesh Refinement (AMR).<\/p>\n<p style=\"text-align: center;\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5410\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-300x234.png\" alt=\"\" width=\"300\" height=\"234\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-300x234.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-1024x799.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-768x599.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-1536x1198.png 1536w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1-385x300.png 385w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/07\/20210731_IoSimulation_Galileo-1.png 1962w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p style=\"text-align: left;\">\n<p>[1] M. Kivelson, K. Khurana, C. Russell, R. Walker, S. Joy, J. Mafi, GALILEO ORBITER AT JUPITER CALIBRATED MAG HIGH RES V1.0, GO-J-MAG-3-RDR-HIGHRES-V1.0, Technical Report, NASA Planetary Data System, 1997.<br \/>\n[2] J. Saur, F. M. Neubauer, J. E. P. Connerney, Plasma interaction of Io with its plasma torus, 2004.<br \/>\n[3] https:\/\/solarsystem.nasa.gov\/resources\/2379\/io-3d-model\/<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/06\/28\/talk-on-2021-07-02-on-wave-propagation-characteristics-upwind-sbp-properties-and-energy-stability-of-dg-viscous-flux-discretizations\/\">Talk on 2021-07-02: On Wave Propagation Characteristics, Upwind SBP Properties and Energy Stability of DG Viscous Flux Discretizations<\/a> <span class=\"date\">June 28, 2021<\/span><div class=\"content\">\n<p><em>Speaker<\/em>: <strong>Dr. Sigrun Ortleb<\/strong>, <a rel=\"noreferrer noopener\" href=\"https:\/\/www.mathematik.uni-kassel.de\/~ortleb\/\" target=\"_blank\">University of Kassel<\/a>, Germany<br><em>Date &amp; time<\/em>: <strong>Friday, 2nd July 2021, 10 am (CEST)<\/strong><br><em>Venue<\/em>: Please request the Zoom meeting link from <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\" data-type=\"page\" data-id=\"3741\">Michael Schlottke-Lakemper<\/a><\/p>\n\n\n\n<p><em>Abstract<\/em>:<br>Regarding accuracy and stability of numerical schemes for computational fluid dynamics, the investigation of diffusion\/dispersion errors depending on the wave number is of utmost importance. Especially for high order methods, a desired small numerical dissipation competes with robustness and thus has to be carefully analyzed. This wave propagation analysis is often based on pure advection problems. In the literature, various approaches to discretize diffusion terms within a DG scheme have been introduced since the discretization of higher order spatial derivatives within the DG framework is less natural than in case of first order derivatives. In this talk, we will address significant differences in the disspation\/dispersion properties for linear advection-diffusion, depending on the specific DG viscous flux discretization which is employed. In addition, results on energy stability of DG viscous flux formulations are dealt with and we show how to formulate the well-known LDG and BR1 fluxes in terms of global upwind SBP operators which complements the derivation and analysis regarding element level SBP properties of the DG scheme.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/06\/24\/snapshot-evolution-of-a-magnetized-torus-differentially-rotating-around-a-static-point-mass\/\">Snapshot: Evolution of a magnetized torus differentially rotating around a static point mass<\/a> <span class=\"date\">June 24, 2021<\/span><div class=\"content\">\r\n<div class=\"wp-block-columns is-layout-flex wp-container-8 wp-block-columns-is-layout-flex\">\r\n<div class=\"wp-block-column is-layout-flow wp-block-column-is-layout-flow\" style=\"flex-basis: 100%;\">\r\n<p class=\"wp-block-preformatted\">Evolution of a magnetized torus differentially rotating around a static point mass. [1] The torus gradually forms into an accretion disk with expanding filaments into the surrounding space due to magneto-rotational instabilities similar to the rich structure of the corona of our sun. The simulation was done with a new Discontinuous-Galerkin based MHD solver [2] in FLASH [3]. Fig. 1 and Fig. 2 show the log. scale density cross section in the x-y and x-z plane. Fig 3 highlights the rich structure of the magnetic field (white stream lines) overlayed on top of the log. scale magnetic pressure in the x-z plane.<\/p>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n<figure class=\"wp-block-gallery columns-3 is-cropped wp-block-gallery-9 is-layout-flex wp-block-gallery-is-layout-flex\">\r\n<ul class=\"blocks-gallery-grid\">\r\n<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7.png\"><img decoding=\"async\" loading=\"lazy\" width=\"784\" height=\"732\" class=\"wp-image-5364\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7.png\" alt=\"\" data-id=\"5364\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5364\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7.png 784w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7-300x280.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7-768x717.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-7-321x300.png 321w\" sizes=\"(max-width: 784px) 100vw, 784px\" \/><\/a><\/figure>\r\n<\/li>\r\n<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-8.png\"><img decoding=\"async\" loading=\"lazy\" width=\"769\" height=\"732\" class=\"wp-image-5365\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-8.png\" alt=\"\" data-id=\"5365\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5365\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-8.png 769w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-8-300x286.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-8-315x300.png 315w\" sizes=\"(max-width: 769px) 100vw, 769px\" \/><\/a><\/figure>\r\n<\/li>\r\n<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6.png\"><img decoding=\"async\" loading=\"lazy\" width=\"779\" height=\"732\" class=\"wp-image-5366\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6.png\" alt=\"\" data-id=\"5366\" data-full-url=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6.png\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5366\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6.png 779w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6-300x282.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6-768x722.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/image-6-319x300.png 319w\" sizes=\"(max-width: 779px) 100vw, 779px\" \/><\/a><\/figure>\r\n<\/li>\r\n<li>\u00a0<\/li>\r\n<\/ul>\r\n<\/figure>\r\n\r\n\r\n\r\n<figure class=\"wp-block-gallery aligncenter columns-1 is-cropped wp-block-gallery-11 is-layout-flex wp-block-gallery-is-layout-flex\">\r\n<ul class=\"blocks-gallery-grid\">\r\n<li class=\"blocks-gallery-item\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/Torus_optimized.gif\"><img decoding=\"async\" loading=\"lazy\" width=\"500\" height=\"500\" class=\"wp-image-5368 aligncenter\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/Torus_optimized.gif\" alt=\"\" data-id=\"5368\" data-full-url=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/06\/Torus_optimized.gif\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5368\" \/><\/a><\/figure>\r\n<\/li>\r\n<\/ul>\r\n<\/figure>\r\n\r\n\r\n\r\n<p>&nbsp;<\/p>\r\n<p class=\"wp-block-preformatted\">[1] Machida, Mami, Mitsuru R. Hayashi, and R. Matsumoto. <br \/>&#8220;Global simulations of differentially rotating magnetized disks:<br \/>Formation of low-beta filaments and structured coronae.&#8221;<br \/>The Astrophysical Journal Letters 532.1 (2000): L67.<br \/>[2] Markert et al. &#8220;A Discontinuous Galerkin Method with Sub\u2011Cell Adaptive Shock Capturing for<br \/>FLASH&#8221; (in preparation)<br \/>[3] <a href=\"http:\/\/flash.uchicago.edu\/\" target=\"_blank\" rel=\"noreferrer noopener\">http:\/\/flash.uchicago.edu\/<\/a><\/p>\r\n\r\n\r\n\r\n<p>&nbsp;<\/p>\r\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/05\/25\/snapshot-gingerbread-man-simulation-with-unstructured-quadrilateral-elements\/\">Snapshot: Gingerbread man simulation with unstructured quadrilateral elements<\/a> <span class=\"date\">May 25, 2021<\/span><div class=\"content\">\n<p>Here we are simulating a gingerbread man, whereby the gingerbread man itself has&nbsp;<code>861<\/code>&nbsp;unstructured quadrilateral elements. For the run parameters, we used&nbsp;<code>polydeg = 10<\/code>&nbsp;with the compressible Euler equations and the manufactured (convergence test) solution. The visualization is done on&nbsp;<code>nvisnodes=12<\/code>&nbsp;uniform points in each direction, on each element.<br><\/p>\n\n\n\n<figure class=\"wp-block-image size-large\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/05\/ginger_optimize.gif\"><img decoding=\"async\" loading=\"lazy\" width=\"600\" height=\"393\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/05\/ginger_optimize.gif\" alt=\"\" class=\"wp-image-5357\"\/><\/a><\/figure>\n\n\n\n<p>We thank Prof. Dr. Andrew Winters (<a href=\"https:\/\/liu.se\/en\/employee\/andwi94\" data-type=\"URL\" data-id=\"https:\/\/liu.se\/en\/employee\/andwi94\">https:\/\/liu.se\/en\/employee\/andwi94<\/a>) for providing this video\/simulation. <\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/05\/19\/talk-on-2021-05-21-on-a-linear-stability-issue-of-split-form-schemes-for-compressible-flows\/\">Talk on 2021-05-21: On a linear stability issue of split form schemes for compressible flows<\/a> <span class=\"date\">May 19, 2021<\/span><div class=\"content\">\n<p><em>Speaker:<\/em> <strong>Dr. Vikram Singh<\/strong>, <a rel=\"noreferrer noopener\" href=\"https:\/\/mpimet.mpg.de\/en\/staff\/vikram-singh\" target=\"_blank\">Max-Planck-Institute for Meteorology<\/a><br><em>Date &amp; time:<\/em> <strong>Friday, 21st May 2021, 10 am (CEST)<\/strong><br><em>Meeting link:<\/em> Please request via email from <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\" data-type=\"page\" data-id=\"3741\">Michael Schlottke-Lakemper<\/a><\/p>\n\n\n\n<p><em>Abstract:<\/em><br>Split form schemes for Euler and Navier-Stokes equations are useful for computation of turbulent flows due to their better robustness. This is because they satisfy additional conservation properties of the governing equations like kinetic energy preservation leading to a reduction in aliasing errors at high orders. Recently, linear stability issues have been pointed out for these schemes for a density wave problem and we investigate this behaviour for some standard split forms. By deriving linearized equations of split form schemes, we show that most existing schemes do not satisfy a perturbation energy equation that holds at the continuous level. A simple modification to the energy flux of some existing schemes is shown to yield a scheme that is consistent with the energy perturbation equation. Numerical tests are given using a discontinuous Galerkin method to demonstrate these results.<\/p>\n\n\n\n<p><em>Note<\/em>: This meeting will be recorded.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/05\/12\/new-paper-submitted-on-the-theoretical-foundation-of-overset-grid-methods-for-hyperbolic-problems-well-posedness-and-conservation\/\">New paper submitted: On the Theoretical Foundation of Overset Grid Methods for Hyperbolic Problems: Well-Posedness and Conservation<\/a> <span class=\"date\">May 12, 2021<\/span><div class=\"content\"><p>We use the energy method to study the well-posedness of initial-boundary value problems approximated by overset mesh methods in one and two space dimensions for linear constant-coefficient hyperbolic systems. We show that in one space dimension, for both scalar equations and systems of equations, the problem where one domain partially oversets another is well-posed when characteristic coupling conditions are used. If a system cannot be diagonalized, as is ususally the case in multiple space dimensions, then the energy method does not give proper bounds in terms of initial and boundary data. For those problems, we propose a novel penalty approach. We show, by using a global energy that accounts for the energy in the overlap region of the domains, that under well-defined conditions on the coupling matrices the penalized overset domain problems are energy bounded, conservative, well-posed and have solutions equivalent to the original single domain problem.<\/p>\n<p>Preprint available at: https:\/\/arxiv.org\/abs\/2105.04664<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/05\/08\/talk-on-2021-05-10-a-two-dimensional-stabilized-discontinuous-galerkin-method-on-curvilinear-embedded-boundary-grids\/\">Talk on 2021-05-10: A two-dimensional stabilized discontinuous Galerkin method on curvilinear embedded boundary grids<\/a> <span class=\"date\">May 8, 2021<\/span><div class=\"content\">\n<p><em>Speaker:<\/em>\u00a0<strong>Dr. Andrew Giuliani<\/strong>, <a href=\"https:\/\/cims.nyu.edu\/~ag6961\/\" target=\"_blank\" rel=\"noreferrer noopener\">New York University<\/a><br><em>Date &amp; time:<\/em>\u00a0<strong>Monday, 10th May 2021, 4pm (CEST)\/10am (EDT)<\/strong><br><em>Meeting link: <\/em>Please request via email from <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\" data-type=\"page\" data-id=\"3741\">Michael Schlottke-Lakemper<\/a><\/p>\n\n\n\n<p><em>Abstract:<\/em><br>In this talk, we present a state redistribution method for high order discontinuous Galerkin methods on curvilinear embedded boundary grids. State redistribution relaxes the overly restrictive CFL condition that results from arbitrarily small cut cells when explicit time steppers are used. Thus, the scheme can take time steps that are proportional to the size of cells in the background grid. The discontinuous Galerkin scheme is stabilized by postprocessing the numerical solution after each stage or step of an explicit time stepping method. The advantage of this approach is that it uses only basic mesh information that is already available in many cut cell codes and does not require complex geometric manipulations. We prove that state redistribution is conservative and p-exact. Finally, we solve a number of test problems that demonstrate the encouraging potential of this technique for applications on curvilinear embedded geometries.&nbsp;Numerical experiments reveal that our scheme converges with order $p+1$ in $L_1$ and between $p$ and $p+1$ in $L_\\infty$.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/05\/04\/new-paper-published-stability-of-discontinuous-galerkin-spectral-element-schemes-for-wave-propagation-when-the-coefficient-matrices-have-jumps\/\">New paper published: Stability of Discontinuous Galerkin Spectral Element Schemes for Wave Propagation when the Coefficient Matrices have Jumps<\/a> <span class=\"date\">May 4, 2021<\/span><div class=\"content\"><p>Abstract:<\/p>\n<p>We use the behavior of the <span id=\"MathJax-Element-1-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-1\" class=\"math\"><span id=\"MathJax-Span-2\" class=\"mrow\"><span id=\"MathJax-Span-3\" class=\"msubsup\"><span id=\"MathJax-Span-4\" class=\"mi\">L<\/span><span id=\"MathJax-Span-5\" class=\"texatom\"><span id=\"MathJax-Span-6\" class=\"mrow\"><span id=\"MathJax-Span-7\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm of the solutions of linear hyperbolic equations with discontinuous coefficient matrices as a surrogate to infer stability of discontinuous Galerkin spectral element methods (DGSEM). Although the <span id=\"MathJax-Element-2-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-8\" class=\"math\"><span id=\"MathJax-Span-9\" class=\"mrow\"><span id=\"MathJax-Span-10\" class=\"msubsup\"><span id=\"MathJax-Span-11\" class=\"mi\">L<\/span><span id=\"MathJax-Span-12\" class=\"texatom\"><span id=\"MathJax-Span-13\" class=\"mrow\"><span id=\"MathJax-Span-14\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm is not bounded by the initial data for homogeneous and dissipative boundary conditions for such systems, the <span id=\"MathJax-Element-3-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-15\" class=\"math\"><span id=\"MathJax-Span-16\" class=\"mrow\"><span id=\"MathJax-Span-17\" class=\"msubsup\"><span id=\"MathJax-Span-18\" class=\"mi\">L<\/span><span id=\"MathJax-Span-19\" class=\"texatom\"><span id=\"MathJax-Span-20\" class=\"mrow\"><span id=\"MathJax-Span-21\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm is easier to work with than a norm that discounts growth due to the discontinuities. We show that the DGSEM with an upwind numerical flux that satisfies the Rankine-Hugoniot (or conservation) condition has the same energy bound as the partial differential equation does in the <span id=\"MathJax-Element-4-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-22\" class=\"math\"><span id=\"MathJax-Span-23\" class=\"mrow\"><span id=\"MathJax-Span-24\" class=\"msubsup\"><span id=\"MathJax-Span-25\" class=\"mi\">L<\/span><span id=\"MathJax-Span-26\" class=\"texatom\"><span id=\"MathJax-Span-27\" class=\"mrow\"><span id=\"MathJax-Span-28\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm, plus an added dissipation that depends on how much the approximate solution fails to satisfy the Rankine-Hugoniot jump.<\/p>\n<p>Accepted in Journal of Scientific Computing.<\/p>\n<p>arXiv link:\u00a0<a href=\"https:\/\/arxiv.org\/abs\/2011.11746\">https:\/\/arxiv.org\/abs\/2011.11746<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/04\/29\/snapshot-astrophysical-colliding-flow-simulation-run-by-a-novel-discontinuous-galerkin-finite-volume-dgfv-blending-scheme-in-flash\/\">Snapshot: Astrophysical colliding flow simulation run by a novel Discontinuous Galerkin\/Finite Volume (DGFV) blending scheme in FLASH<\/a> <span class=\"date\">April 29, 2021<\/span><div class=\"content\">\r\n<div class=\"wp-block-group\">\r\n<div class=\"wp-block-group__inner-container is-layout-flow wp-block-group-is-layout-flow\">\r\n<h6 class=\"wp-block-preformatted\">Simulation of an astrophysical colliding flow [1] run by a novel Discontinuous Galerkin\/Finite Volume (DGFV) blending scheme [2] which has been implemented in the FLASH code [3].<\/h6>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n<figure class=\"wp-block-gallery columns-1 is-cropped wp-block-gallery-14 is-layout-flex wp-block-gallery-is-layout-flex\">\r\n<figure><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image.png\"><img decoding=\"async\" loading=\"lazy\" width=\"1024\" height=\"288\" class=\"wp-image-5333\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image-1024x288.png\" alt=\"\" data-id=\"5333\" data-full-url=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image.png\" data-link=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/?attachment_id=5333\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image-1024x288.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image-300x84.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image-768x216.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image-500x141.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/image.png 1276w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/a><\/figure>\r\n<\/figure>\r\n\r\n\r\n\r\n<div class=\"wp-block-group\">\r\n<div class=\"wp-block-group__inner-container is-layout-flow wp-block-group-is-layout-flow\">\r\n<p class=\"wp-block-preformatted\">The code has following capabilities:<\/p>\r\n<ul>\r\n<li class=\"wp-block-preformatted\">fourth-order accurate ideal magneto-hydrodynamics with hyperbolic divergence cleaning [4]<\/li>\r\n<li class=\"wp-block-preformatted\">octree-based adaptive mesh refinement<\/li>\r\n<li class=\"wp-block-preformatted\">distributive computing and load balancing<\/li>\r\n<li class=\"wp-block-preformatted\">multi-species fluid dynamics (N_species &gt; 10)<\/li>\r\n<li class=\"wp-block-preformatted\">turbulent driving<\/li>\r\n<li class=\"wp-block-preformatted\">octree-based Poisson solver for self-gravity [5]<\/li>\r\n<li class=\"wp-block-preformatted\">octree-based radiation physics [6]<\/li>\r\n<li class=\"wp-block-preformatted\">external gravitional fields<\/li>\r\n<li class=\"wp-block-preformatted\">sink particles [7]<\/li>\r\n<li class=\"wp-block-preformatted\">chemical reaction networks [8]<\/li>\r\n<\/ul>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n<div class=\"wp-block-group\">\r\n<div class=\"wp-block-group__inner-container is-layout-flow wp-block-group-is-layout-flow\">\r\n<p class=\"wp-block-preformatted\">[1] Weis, Micheal et al. &#8220;The Virial Balance of CO-Substructures in Colliding Magnetised Flows&#8221; (in preparation)<br \/>[2] Markert, Johannes et al. &#8220;A Sub-Element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods&#8221; (submitted)<br \/>[3] Fryxell, Bruce, et al. &#8220;FLASH: An adaptive mesh hydrodynamics code for modeling astrophysical thermonuclear flashes.&#8221; The Astrophysical Journal Supplement Series 131.1 (2000): 273.<br \/>[4] Markert, Johannes et al. &#8220;Flash goes DG&#8221; (working title, in preparation)<br \/>[5] W\u00fcnsch, Richard, et al. &#8220;Tree-based solvers for adaptive mesh refinement code FLASH\u2013I: gravity and optical depths.&#8221; Monthly Notices of the Royal Astronomical Society 475.3 (2018): 3393-3418.<br \/>[6] W\u00fcnsch, Richard et al. &#8220;Tree-based solvers for adaptive mesh refinement code FLASH &#8211; II: radiation transport module TreeRay&#8221; (submitted)<br \/>[7] Federrath, Christoph, et al. &#8220;Modeling collapse and accretion in turbulent gas clouds: implementation and comparison of sink particles in AMR and SPH.&#8221; The Astrophysical Journal 713.1 (2010): 269.<br \/>[8] Seifried, D., and S. Walch. &#8220;Modelling the chemistry of star-forming filaments\u2013I. H2 and CO chemistry.&#8221; Monthly Notices of the Royal Astronomical Society: Letters 459.1 (2016): L11-L15.<\/p>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n<p>&nbsp;<\/p>\r\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/04\/01\/snapshot-simulation-of-a-kelvin-helmholtz-instability-using-second-order-finite-volume-schemes-and-fourth-order-discontinuous-galerkin-methods\/\">Snapshot: Simulation of a Kelvin-Helmholtz instability using second order Finite Volume schemes and fourth order Discontinuous Galerkin methods<\/a> <span class=\"date\">April 1, 2021<\/span><div class=\"content\"><p>We present in-viscid and viscous simulations of a Kelvin-Helmholtz instability using second a order accurate monotoniced-central finite volume (FV) method and a fourth order accurate discontinuous Galerkin (DG) method. The initial condition is given by [1]:<\/p>\n<p>$$\\rho (t=0) = \\frac{1}{2}<br \/>\n+ \\frac{3}{4} B,<br \/>\n~~~~~~~~~<br \/>\np (t=0) = 1,~~~~~~~~~~<br \/>\n$$<\/p>\n<p>$$<br \/>\nv_1 (t=0) = \\frac{1}{2} \\left( B-1 \\right),<br \/>\n~~~~~~~<br \/>\nv_2 (t=0) = \\frac{1}{10} \\sin(2 \\pi x),<br \/>\n$$<\/p>\n<p>with $$B=\\tanh \\left( 15 y + 7.5 \\right) &#8211; \\tanh(15y-7.5).$$<\/p>\n<p>We first present the FV results at end time $t=3.7$, which are computed using a monotoniced-central second order discretization of the Euler equations of gas dynamics on uniform grids.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5307\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-1022x1024.png 1022w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-768x769.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_2048_Density.png 1363w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5308\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-768x767.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_NEMO_4096_Density.png 1367w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>The next results use a fourth order DG discretization of the Navier-Stokes equations on uniform grids using $Re=320.000$ at end time $t=3.7$. The highest resolution (4096\u00b2 DOFs) is a direct numerical simulation (DNS) of the problem, where all scales are resolved.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5309\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-768x768.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_2048_Density_JohannesScale.png 1368w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5310\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-300x300.png\" alt=\"\" width=\"300\" height=\"300\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-300x300.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-1024x1024.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-150x150.png 150w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-768x769.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/04\/cropLegend_FLUXO_4096_Density.png 1367w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>It is remarkable that the numerical dissipation of the second order FV scheme causes the in-viscid simulation with 2048\u00b2 DOFs to look very similar to the viscous DNS solution at $Re=320.000$.<\/p>\n<p>[1] A.M. Rueda-Ram\u00edrez, G.J Gassner (2021). <em>A Subcell Finite Volume Positivity-Preserving Limiter for DGSEM Discretizations of the Euler Equations<\/em>. <a href=\"https:\/\/arxiv.org\/pdf\/2102.06017.pdf\">https:\/\/arxiv.org\/pdf\/2102.06017.pdf<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/03\/12\/talk-on-2021-03-18-a-tour-of-bifurcationkit-and-some-results-on-mean-fields-of-spiking-neurons\/\">Talk on 2021-03-18: A tour of BifurcationKit and some results on mean fields of spiking neurons<\/a> <span class=\"date\">March 12, 2021<\/span><div class=\"content\">\n<p><em>Speaker<\/em>: <strong>Dr. Romain Veltz<\/strong>, <a rel=\"noreferrer noopener\" aria-label=\"INRIA (opens in a new tab)\" href=\"http:\/\/romainveltz.pythonanywhere.com\/\" target=\"_blank\">INRIA<\/a>, France<br><em>Date: <\/em><strong>Thursday, 18th March 2021, 10am (CET)<\/strong><br><em>Zoom-Link:<\/em> Please request via email from <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\">Michael Schlottke-Lakemper<\/a><\/p>\n\n\n\n<h2 class=\"wp-block-heading\">Abstract<\/h2>\n\n\n\n<p>In this talk, I will first present the basics of bifurcation theory. Then, I will give a panorama of\u00a0<code><a href=\"https:\/\/github.com\/rveltz\/BifurcationKit.jl\" target=\"_blank\" rel=\"noreferrer noopener\" aria-label=\"BifurcationKit.jl (opens in a new tab)\">BifurcationKit.jl<\/a><\/code>, a Julia package to perform numerical bifurcation analysis of large dimensional equations (PDE, nonlocal equations, etc) possibly on GPUs using Matrix-Free \/ Sparse Matrix formulations of the problem. Julia programming language gives access to a rich ecosystem (PDE, GPU, AD, cluster&#8230;). Notably, numerical bifurcation analysis can be done\u00a0<strong>entirely<\/strong>\u00a0on GPU as will be shown in an example.<\/p>\n\n\n\n<p><code>BifurcationKit<\/code>\u00a0incorporates continuation algorithms (PALC, deflated continuation, &#8230;) which can be used to perform fully automatic bifurcation diagram computation of stationary states. I will showcase this with the 2d Bratu problem.<\/p>\n\n\n\n<p>Additionally, by leveraging on the above methods, the package can also seek for periodic orbits of Cauchy problems by casting them into an equation of high dimension. It is by now, one of the only software which provides parallel (Standard \/ Poincar\u00e9) shooting methods and finite differences based methods to compute periodic orbits in high dimensions. I will present an application highlighting the ability to fine tune\u00a0<code>BifurcationKit<\/code>\u00a0to get performance.<\/p>\n\n\n\n<p>In a last part, I will describe a mean field model of stochastic spiking neurons described with a 2d measure valued equation. I will present a numerical scheme based on an implicit Finite Volume method. I will then provide some mathematical properties of the mean field concerning well posedness and stationary solutions. Additionally, I will show how\u00a0<code>BifurcationKit.jl<\/code>\u00a0can be used to study numerically the model. Finally, I will conclude on open problems, some of which could hopefully be tackled numerically with\u00a0<code><a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" target=\"_blank\" rel=\"noreferrer noopener\" aria-label=\"Trixi.jl (opens in a new tab)\">Trixi.jl<\/a><\/code>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/02\/25\/new-proceedings-paper-submitted-a-subcell-finite-volume-positivity-preserving-limiter-for-dgsem-discretizations-of-the-euler-equations\/\">New proceedings paper submitted: A Subcell Finite Volume Positivity-Preserving Limiter for DGSEM Discretizations of the Euler Equations<\/a> <span class=\"date\">February 25, 2021<\/span><div class=\"content\"><p>In this paper, we present a positivity-preserving limiter for nodal Discontinuous Galerkin disctretizations of the compressible Euler equations. We use a Legendre-Gauss-Lobatto (LGL) Discontinuous Galerkin Spectral Element Method (DGSEM) and blend it locally with a consistent LGL-subcell Finite Volume (FV) discretization using a hybrid FV\/DGSEM scheme that was recently proposed for entropy stable shock capturing. We show that our strategy is able to ensure robust simulations with positive density and pressure when using the standard and the split-form DGSEM. Furthermore, we show the applicability of our FV positivity limiter in extremely under-resolved vortex dominated simulations and in problems with shocks.<\/p>\n<p>Preprint available at: <a href=\"https:\/\/arxiv.org\/pdf\/2102.06017.pdf\">https:\/\/arxiv.org\/pdf\/2102.06017.pdf<\/a><\/p>\n<p><div id=\"attachment_5269\" style=\"width: 310px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1.png\"><img aria-describedby=\"caption-attachment-5269\" decoding=\"async\" loading=\"lazy\" class=\"size-medium wp-image-5269\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1-300x185.png\" alt=\"Evolution of the density for a Sedov blast simulation with periodic boundaries\" width=\"300\" height=\"185\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1-300x185.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1-768x474.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1-1024x632.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1-486x300.png 486w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/PictureNewsWebsite-1.png 1160w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><p id=\"caption-attachment-5269\" class=\"wp-caption-text\">Evolution of the density for a Sedov blast simulation with periodic boundaries<\/p><\/div><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/02\/24\/snapshot-accuracy-of-the-lgl-subcell-fv-scheme-with-and-without-inner-element-reconstruction\/\">Snapshot: Accuracy of the LGL-subcell FV scheme with and without inner-element reconstruction<\/a> <span class=\"date\">February 24, 2021<\/span><div class=\"content\">\n<p>We are interested in the accuracy of the finite volume scheme, on a LGL subcell grid induced by a LGL-DGSEM discretization. For comparison, we look at a Kelvin-Helmholtz-Instability test problem, simulated with 4\u00b2 LGL nodes per element and 32\u00b2 elements, i.e. 128\u00b2 spatial degrees of freedom. The high-order DGSEM uses the entropy-conserving split-form powered by the flux of Chandrashekar in the volume, and the HLLC flux at the element interfaces. Positivity is controlled by adding subcell FV where the positivity is not fulfilled with the amount needed. All FV discretizations use the HLLC flux at the element interfaces and at the subcell interfaces as well.\u00a0\u00a0<br><br>The first results show the high-order DGSEM result, which is formally 4<sup>th<\/sup>\u00a0order accurate for smooth problems:<br><br><\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" loading=\"lazy\" width=\"570\" height=\"392\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/snap1.png\" alt=\"\" class=\"wp-image-5264\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/snap1.png 570w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/snap1-300x206.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/snap1-436x300.png 436w\" sizes=\"(max-width: 570px) 100vw, 570px\" \/><\/figure>\n\n\n\n<p>The next result uses a subcell finite volume approximation on the LGL-subcell grid, directly, i.e. without spatial reconstruction (piecewise constant approximation). The result is very dissipative:<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" loading=\"lazy\" width=\"592\" height=\"418\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap2.png\" alt=\"\" class=\"wp-image-5263\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap2.png 592w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap2-300x212.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap2-425x300.png 425w\" sizes=\"(max-width: 592px) 100vw, 592px\" \/><\/figure>\n\n\n\n<p>The last results show the improvement when using a piecewise linear reconstruction with monotonized-central slope limiter on the LGL-subcell grid. The reconstruction is inner-element only, as it does not use element neighbor information. Thus, at the very first and very last LGL subcell, no reconstruction is used. Still, the accuracy is recovered nicely by this simple modification of the subcell FV scheme:\u00a0<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" loading=\"lazy\" width=\"689\" height=\"474\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap3.png\" alt=\"\" class=\"wp-image-5262\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap3.png 689w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap3-300x206.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/Snap3-436x300.png 436w\" sizes=\"(max-width: 689px) 100vw, 689px\" \/><\/figure>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/01\/30\/snapshot-single-node-performance-comparison-of-four-different-magnetohydrodynamics-codes\/\">Snapshot: Single-Node Performance Comparison of Four Different Magnetohydrodynamics Codes<\/a> <span class=\"date\">January 30, 2021<\/span><div class=\"content\"><p>We compare the runtime performance of four different magnetohydrodynamics codes\u00a0on a single compute node on the in-house high performance cluster ODIN. A\u00a0compute node on ODIN consists of 16 cores. We run the &#8216;3D Alfven wave&#8217; test\u00a0case up to a fixed simulation time and measure the elapsed wall clock time of\u00a0each code minus initialization time and input\/output operations. For each run\u00a0the number of cores is successively increased. This allows us to get insights\u00a0into the scaling behavior (speedup gain wih increasing number of cores) on a\u00a0single compute node. Furthermore we plot the performance index (PID) over the\u00a0number of nodes which is a measure of throughput, i.e. how many millions of\u00a0data points (degrees-of-freedom\/DOF) per second are processed by the each\u00a0code.<\/p>\n<p>The four contestants are:<\/p>\n<ul>\n<li>Flash [1] with Paramesh 4.0 and the Five-wave Bouchut Finite-Volume solver written in Fortran<\/li>\n<li>Fluxo [2] an MHD Discontinuous Galerkin Spectral Element code written in Fortran<\/li>\n<li>Trixi [3] an MHD Discontinuous Galerkin Spectral Element code written in Julia<\/li>\n<li>Nemo an in-house prototype code providing a 2nd order monotonized-central MHD-FV scheme (MCFV) and a\u00a0hybrid MHD Discontinuous Galerkin Spectral Element \/ Finite Volume scheme (DGFV) written in Fortran.<\/li>\n<\/ul>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5242\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2-300x166.png\" alt=\"\" width=\"300\" height=\"166\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2-300x166.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2-768x426.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2-500x277.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-2.png 942w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a> <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5243\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1-300x166.png\" alt=\"\" width=\"300\" height=\"166\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1-300x166.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1-768x426.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1-1024x567.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1-500x277.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2021\/02\/image-1.png 1884w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>Trixi uses a thread-based parallelization approach while Flash, Fluxo and Nemo<br \/>\nuse the standard MPI method. Furthermore, Nemo also provides OpenMP-based<br \/>\nparallelization.<\/p>\n<p>[1] https:\/\/flash.uchicago.edu\/site\/<br \/>\n[2] https:\/\/github.com\/project-fluxo\/fluxo<br \/>\n[3] https:\/\/github.com\/trixi-framework\/Trixi.jl<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/01\/25\/talk-julia-for-adaptive-high-order-multi-physics-simulations\/\">Talk on 2021-01-21: Julia for adaptive high-order multi-physics simulations<\/a> <span class=\"date\">January 25, 2021<\/span><div class=\"content\">\n<p>On <strong>Wednesday, 27th January 2021, 3:15pm CET<\/strong>, <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\/\">Michael Schlottke-Lakemper<\/a> will give an online talk on<\/p>\n\n\n\n<p><strong>Julia for adaptive high-order simulations<\/strong><\/p>\n\n\n\n<p>To obtain the Zoom link, please contact the organizers via the <a rel=\"noreferrer noopener\" aria-label=\"official meeting announcement (opens in a new tab)\" href=\"http:\/\/www.maths.lu.se\/kalendarium\/?evenemang=numerical-analysis-seminar-0\" target=\"_blank\">official meeting announcement<\/a>.<\/p>\n\n\n\n<p><strong>Abstract<\/strong><\/p>\n\n\n\n<p>Julia has been touted as a programming language especially \nwell-suited for numerical analysis and scientific computing. However, \nwhile its prevalence is steadily increasing, it has not yet seen \nwidespread adoption in the computational science or high-performance \ncomputing communities. One of the hurdles is a (perceived) lack of \nreal-world examples that show how Julia can be used to conduct numerical\n simulations and what its advantages and drawbacks are for scientific \napplications.<br>\nTo remediate this, in this talk we discuss the development of a \npurely hyperbolic method for self-gravitating gas dynamics within our \nJulia-based open source simulation framework Trixi.jl (<a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\">https:\/\/github.com\/trixi-framework\/Trixi.jl<\/a>).\n In this approach, we reformulate the elliptic gravity problem into a \nhyperbolic diffusion problem, which is solved in pseudotime using the \nsame explicit high-order discontinuous Galerkin method we use for the \nflow solution. A key benefit is that in the resulting multi-physics \nsimulation problem, we can reuse existing hyperbolic solvers while \nretaining advanced features such as non-conforming and solution-adaptive\n meshes.<br>\nNext to presenting numerical results, we will critically examine \nour experience with building a multi-physics simulation framework with \nJulia. We will discuss its strengths and weaknesses as a programming \nlanguage for research software engineering, including an assessment of \nhow Julia\u2019s claimed benefits hold up against scientific reality, and \ngive a live demonstration of Julia and Trixi.jl in action.<br>\n<br>\nTo make the shown examples reproducible by the audience, the Jupyter notebook used for the live demonstration is available at&nbsp;<a href=\"https:\/\/github.com\/trixi-framework\/talk-2021-julia-adaptive-multi-physics-simulations\">https:\/\/github.com\/trixi-framework\/talk-2021-julia-adaptive-multi-physics-simulations<\/a>.\n It can be either run from a local Julia\/Jupyter installation or in the \ncloud via Binder (without having to install Julia locally).<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2021\/01\/04\/new-paper-submitted-an-entropy-stable-nodal-discontinuous-galerkin-method-for-the-resistive-mhd-equations-part-ii-subcell-finite-volume-shock-capturing\/\">New paper submitted: An Entropy Stable Nodal Discontinuous Galerkin Method for the resistive MHD Equations. Part II: Subcell Finite Volume Shock Capturing<\/a> <span class=\"date\">January 4, 2021<\/span><div class=\"content\"><p>The second paper of this series presents two robust entropy stable shock-capturing methods for discontinuous Galerkin spectral element (DGSEM) discretizations of the compressible magneto-hydrodynamics (MHD) equations. Specifically, we use the resistive GLM-MHD equations, which include a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM). For the continuous entropy analysis to hold, and due to the divergence-free constraint on the magnetic field, the GLM-MHD system requires the use of non-conservative terms, which need special treatment.<\/p>\n<p>Hennemann et al. [DOI:<a class=\"link-https link-external\" href=\"https:\/\/arxiv.org\/ct?url=https%3A%2F%2Fdx.doi.org%2F10.1016%2Fj.jcp.2020.109935&amp;v=31c3f2c5\" rel=\"external noopener nofollow\" data-doi=\"10.1016\/j.jcp.2020.109935\">10.1016\/j.jcp.2020.109935<\/a>] recently presented an entropy stable shock-capturing strategy for DGSEM discretizations of the Euler equations that blends the DGSEM scheme with a subcell first-order finite volume (FV) method. Our first contribution is the extension of the method of Hennemann et al. to systems with non-conservative terms, such as the GLM-MHD equations. In our approach, the advective and non-conservative terms of the equations are discretized with a hybrid FV\/DGSEM scheme, whereas the visco-resistive terms are discretized only with the high-order DGSEM method. We prove that the extended method is entropy stable on three-dimensional unstructured curvilinear meshes. Our second contribution is the derivation and analysis of a second entropy stable shock-capturing method that provides enhanced resolution by using a subcell reconstruction procedure that is carefully built to ensure entropy stability.<\/p>\n<p>We provide a numerical verification of the properties of the hybrid FV\/DGSEM schemes on curvilinear meshes and show their robustness and accuracy with common benchmark cases, such as the Orszag-Tang vortex and the GEM reconnection challenge. Finally, we simulate a space physics application: the interaction of Jupiter&#8217;s magnetic field with the plasma torus generated by the moon Io.<\/p>\n<p>Preprint available at: <a href=\"https:\/\/arxiv.org\/abs\/2012.12040\">arXiv:2012.12040<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/12\/24\/snapshot-hybrid-discontinuous-galerkin-finite-volume-dg-fv-simulation-of-the-multicomponent-shock-bubble-interaction-problem\/\">Snapshot: Hybrid Discontinuous Galerkin \/ Finite Volume (DG\/FV) simulation of the multicomponent Shock-Bubble Interaction Problem<\/a> <span class=\"date\">December 24, 2020<\/span><div class=\"content\"><p><span class=\"style-scope yt-formatted-string\" dir=\"auto\">Simulation of the Shock-Bubble Interaction Problem with a hybrid entropy-stable DG\/FV method and Adaptive Mesh Refinement (AMR) using Trixi (<\/span><a class=\"yt-simple-endpoint style-scope yt-formatted-string\" dir=\"auto\" spellcheck=\"false\" href=\"https:\/\/www.youtube.com\/redirect?redir_token=QUFFLUhqbThmNlpCcENDVmlwaEwzVFlUaDlxcDk4RnlhUXxBQ3Jtc0trSnc3OXdMY2FDNjVFSG9CSklkbFlETjRaeGgxSl9HNHpISmVlREVVMktZQmYyVVUwQlBUQUFXcF9aamNwUmstTGdTM2pVX2FuM3IxUjdTN205T0JBY1d5di1zaEtTOXdqdnd1VWtxU0hUeUtoeWdRcw%3D%3D&amp;q=https%3A%2F%2Fgithub.com%2Ftrixi-framework%2FTrixi.jl&amp;v=qaBlIU9VgXE&amp;event=video_description\" target=\"_blank\" rel=\"nofollow noopener noreferrer\">https:\/\/github.com\/trixi-framework\/Tr&#8230;<\/a><span class=\"style-scope yt-formatted-string\" dir=\"auto\">). <\/span><\/p>\n<p><span class=\"style-scope yt-formatted-string\" dir=\"auto\">The simulation uses a multicomponent model and computes the spatial operator as a blend of a first-order subcell FV method and a fourth-order DG scheme on a squared domain.<\/span><\/p>\n<p><iframe loading=\"lazy\" title=\"Hybrid DG\/FV simulation of the multicomponent Shock-Bubble-Problem\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/qaBlIU9VgXE?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/12\/15\/new-paper-published-a-novel-robust-strategy-for-discontinuous-galerkin-methods-in-computational-fluid-mechanics-why-when-what-where\/\">New paper published: A Novel Robust Strategy for Discontinuous Galerkin Methods in Computational Fluid Mechanics: Why? When? What? Where?<\/a> <span class=\"date\">December 15, 2020<\/span><div class=\"content\"><p>In this paper we will review a recent emerging paradigm shift in the construction and analysis of high order Discontinuous Galerkin (DG) methods applied to approximate solutions of hyperbolic or mixed hyperbolic-parabolic partial differential equations (PDEs) in computational physics. There is a long history using DG methods to approximate the solution of PDEs in computational physics with successful applications in linear wave propagation, like those governed by Maxwell&#8217;s equations, incompressible and compressible fluid and plasma dynamics governed by the Navier-Stokes and the Magnetohydrodynamics equations, or as a solver for ordinary differential equations (ODEs), e.g., in structural mechanics. The DG method amalgamates ideas from several existing methods such as the Finite Element Galerkin method (FEM) and the Finite Volume method (FVM) and is specifically applied to problems with advection dominated properties, such as fast moving fluids or wave propagation. In the numerics community, DG methods are infamous for being computationally complex and, due to their high order nature, as having issues with robustness, i.e., these methods are sometimes prone to crashing easily. In this article we will focus on efficient nodal versions of the DG scheme and present recent ideas to restore its robustness, its connections to and influence by other sectors of the numerical community, such as the finite difference community, and further discuss this young, but rapidly developing research topic by highlighting the main contributions and a closing discussion about possible next lines of research.<\/p>\n<p>Published in: <a href=\"https:\/\/www.frontiersin.org\/articles\/10.3389\/fphy.2020.500690\/abstract\">https:\/\/www.frontiersin.org\/articles\/10.3389\/fphy.2020.500690\/abstract<\/a><\/p>\n<p>Front. Phys. | doi: 10.3389\/fphy.2020.500690<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/12\/15\/new-paper-submitted-a-split-form-stable-cg-dg-sem-for-wave-propagation-modeled-by-linear-hyperbolic-systems\/\">New paper submitted: A Split-Form, Stable CG\/DG-SEM for Wave Propagation Modeled by Linear Hyperbolic Systems<\/a> <span class=\"date\">December 15, 2020<\/span><div class=\"content\"><p>We present a hybrid continuous and discontinuous Galerkin spectral element approximation that leverages the advantages of each approach. The continuous Galerkin approximation is used on interior element faces where the equation properties are continuous. A discontinuous Galerkin approximation is used at physical boundaries and if there is a jump in properties at a face. The approximation uses a split form of the equations and two-point fluxes to ensure stability for unstructured quadrilateral\/hexahedral meshes with curved elements. The approximation is also conservative and constant state preserving on such meshes. Spectral accuracy is obtained for all examples, which include wave scattering at a discontinuous medium boundary.<\/p>\n<p>Preprint available at: <span class=\"arxivid\"><a href=\"https:\/\/arxiv.org\/abs\/2012.06510\">arXiv:2012.06510<\/a><\/span><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/11\/27\/hybrid-discontinuous-galerkin-finite-volume-dg-fv-simulation-of-the-orszag-tang-vortex\/\">Snapshot: Hybrid Discontinuous Galerkin \/ Finite Volume (DG\/FV) simulation of the Orszag-Tang vortex<\/a> <span class=\"date\">November 27, 2020<\/span><div class=\"content\"><p>Simulation of the Orszag-Tang vortex test with a hybrid entropy-stable DG\/FV method and Adaptive Mesh Refinement (AMR) using FLUXO (<a href=\"https:\/\/github.com\/project-fluxo\/fluxo\">https:\/\/github.com\/project-fluxo\/fluxo<\/a>).<\/p>\n<p>The simulation uses the GLM-MHD model to control div(B) errors and computes the spatial operator as a blend of a first-order subcell FV method and a fourth-order DG scheme.<\/p>\n<p>The initial grid has 16&#215;16 elements (64\u00b2 DOFs), and the maximum resolution is achieved with three refinement levels (equivalent to 512\u00b2 DOFs).<\/p>\n<p><iframe loading=\"lazy\" title=\"Hybrid Discontinuous Galerkin  \/ Finite Volume (DG\/FV) simulation of the Orszag-Tang vortex\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/Tr-3h40T8PU?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/11\/25\/new-paper-submitted-stability-of-discontinuous-galerkin-spectral-element-schemes-for-wave-propagation-when-the-coefficient-matrices-have-jumps\/\">New paper submitted: Stability of Discontinuous Galerkin Spectral Element Schemes for Wave Propagation when the Coefficient Matrices have Jumps<\/a> <span class=\"date\">November 25, 2020<\/span><div class=\"content\"><p>We use the behavior of the <span id=\"MathJax-Element-1-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-1\" class=\"math\"><span id=\"MathJax-Span-2\" class=\"mrow\"><span id=\"MathJax-Span-3\" class=\"msubsup\"><span id=\"MathJax-Span-4\" class=\"mi\">L<\/span><span id=\"MathJax-Span-5\" class=\"texatom\"><span id=\"MathJax-Span-6\" class=\"mrow\"><span id=\"MathJax-Span-7\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm of the solutions of linear hyperbolic equations with discontinuous coefficient matrices as a surrogate to infer stability of discontinuous Galerkin spectral element methods (DGSEM). Although the <span id=\"MathJax-Element-2-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-8\" class=\"math\"><span id=\"MathJax-Span-9\" class=\"mrow\"><span id=\"MathJax-Span-10\" class=\"msubsup\"><span id=\"MathJax-Span-11\" class=\"mi\">L<\/span><span id=\"MathJax-Span-12\" class=\"texatom\"><span id=\"MathJax-Span-13\" class=\"mrow\"><span id=\"MathJax-Span-14\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm is not bounded by the initial data for homogeneous and dissipative boundary conditions for such systems, the <span id=\"MathJax-Element-3-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-15\" class=\"math\"><span id=\"MathJax-Span-16\" class=\"mrow\"><span id=\"MathJax-Span-17\" class=\"msubsup\"><span id=\"MathJax-Span-18\" class=\"mi\">L<\/span><span id=\"MathJax-Span-19\" class=\"texatom\"><span id=\"MathJax-Span-20\" class=\"mrow\"><span id=\"MathJax-Span-21\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm is easier to work with than a norm that discounts growth due to the discontinuities. We show that the DGSEM with an upwind numerical flux that satisfies the Rankine-Hugoniot (or conservation) condition has the same energy bound as the partial differential equation does in the <span id=\"MathJax-Element-4-Frame\" class=\"MathJax\" tabindex=\"0\"><span id=\"MathJax-Span-22\" class=\"math\"><span id=\"MathJax-Span-23\" class=\"mrow\"><span id=\"MathJax-Span-24\" class=\"msubsup\"><span id=\"MathJax-Span-25\" class=\"mi\">L<\/span><span id=\"MathJax-Span-26\" class=\"texatom\"><span id=\"MathJax-Span-27\" class=\"mrow\"><span id=\"MathJax-Span-28\" class=\"mn\">2<\/span><\/span><\/span><\/span><\/span><\/span><\/span> norm, plus an added dissipation that depends on how much the approximate solution fails to satisfy the Rankine-Hugoniot jump.<\/p>\n<p>Preprint available at\u00a0<a href=\"https:\/\/arxiv.org\/abs\/2011.11746\">arXiv:2011.11746<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/11\/23\/new-paper-submitted-a-sub-element-adaptive-shock-capturing-approach-for-discontinuous-galerkin-methods\/\">New paper submitted: A Sub-Element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods<\/a> <span class=\"date\">November 23, 2020<\/span><div class=\"content\">\n<p>In this paper, a new strategy for a sub-element based shock capturing for discontinuous Galerkin (DG) approximations is presented. The idea is to interpret a DG element as a collection of data and construct a hierarchy of low to high order discretizations on this set of data, including a first order finite volume scheme up to the full order DG scheme. The different DG discretizations are then blended according to sub-element troubled cell indicators, resulting in a final discretization that adaptively blends from low to high order within a single DG element. The goal is to retain as much high order accuracy as possible, even in simulations with very strong shocks, as e.g. presented in the Sedov test. The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing. The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.<br><\/p>\n\n\n\n<p>Preprint available at\u00a0<strong><a href=\"https:\/\/arxiv.org\/abs\/2011.03338\">arXiv:2011.03338<\/a><\/strong>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/10\/23\/snapshot-2d-riemann-problem-computed-with-a-hybrid-discontinuous-galerkin-finite-volume-scheme\/\">Snapshot: 2D Riemann problem computed with a hybrid Discontinuous Galerkin \/ Finite Volume scheme<\/a> <span class=\"date\">October 23, 2020<\/span><div class=\"content\">\n<p>We present a 2D Riemann problem (configuration 3) computed with a hybrid Discontinuous Galerkin \/ Finite Volume scheme on a uniform grid with 1024^2 degrees of freedom.<\/p>\n\n\n\n<figure class=\"wp-block-embed-youtube wp-block-embed is-type-video is-provider-youtube wp-embed-aspect-4-3 wp-has-aspect-ratio\"><div class=\"wp-block-embed__wrapper\">\n<iframe loading=\"lazy\" title=\"2D Riemann problem computed with a hybrid Discontinuous Galerkin \/ Finite Volume scheme\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/T4HUW66kOio?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe>\n<\/div><\/figure>\n\n\n\n<p>In this hybrid approach we are blending an 8th order, 4th order and 2nd order DG as well as a 1st order FV scheme. <\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/09\/29\/new-paper-submitted-preventing-pressure-oscillations-does-not-fix-local-linear-stability-issues-of-entropy-based-split-form-high-order-schemes\/\">New paper submitted: Preventing pressure oscillations does not fix local linear stability issues of entropy-based split-form high-order schemes<\/a> <span class=\"date\">September 29, 2020<\/span><div class=\"content\">\n<p>Recently, it was discovered that the entropy-conserving\/dissipative high-order split-form discontinuous Galerkin discretizations have robustness issues when trying to solve the simple density wave propagation example for the compressible Euler equations. The issue is related to missing local linear stability, i.e. the stability of the discretization towards perturbations added to a stable base flow. This is strongly related to an anti-diffusion mechanism, that is inherent in entropy-conserving two-point fluxes, which are a key ingredient for the high-order discontinuous Galerkin extension. In this paper, we investigate if pressure equilibrium preservation is a remedy to these recently found local linear stability issues of entropy-conservative\/dissipative high-order split-form discontinuous Galerkin methods for the compressible Euler equations. Pressure equilibrium preservation describes the property of a discretization to keep pressure and velocity constant for pure density wave propagation. We present the full theoretical derivation, analysis, and show corresponding numerical results to underline our findings. The source code to reproduce all numerical experiments presented in this article is available online (DOI: <a rel=\"noreferrer noopener\" aria-label=\"10.5281\/zenodo.4054366 (opens in a new tab)\" href=\"https:\/\/doi.org\/10.5281\/zenodo.4054366\" target=\"_blank\">10.5281\/zenodo.4054366<\/a>).<\/p>\n\n\n\n<p>Preprint available at <a rel=\"noreferrer noopener\" aria-label=\"arXiv:2009.13139 (opens in a new tab)\" href=\"https:\/\/arxiv.org\/abs\/2009.13139\" target=\"_blank\">arXiv:2009.13139<\/a>. Numerical results were obtained with <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" target=\"_blank\" rel=\"noreferrer noopener\" aria-label=\"Trixi.jl (opens in a new tab)\">Trixi.jl<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/09\/28\/snapshot-purely-hyperbolic-gravity-simulations-with-trixi-jl\/\">Snapshot: Purely hyperbolic gravity simulations with Trixi.jl<\/a> <span class=\"date\">September 28, 2020<\/span><div class=\"content\">\n<p>One of the challenges when simulating astrophysical flows with self-gravity is to compute the gravitational forces. In contrast to the hyperbolic compressible Euler equations, the gravity field is described by an elliptic Poisson equation. In <a href=\"#ref-1\">[1]<\/a>, we present a purely hyperbolic approach by reformulating the elliptic problem into a hyperbolic diffusion problem, which is solved in pseudotime using the same explicit high-order discontinuous Galerkin method we use for the flow solution.<\/p>\n\n\n\n<p>We start with the Poisson equation for the Newtonian gravitational potential $\\phi$,<\/p>\n\n\n\n<p>$$-\\vec{\\nabla}^2\\phi = -4\\pi G \\rho\\label{grav}$$<\/p>\n\n\n\n<p>where G is the universal gravitational constant and $\\rho$ is the mass density. Following <a href=\"https:\/\/doi.org\/10.1016\/j.jcp.2007.07.029\" target=\"_blank\" rel=\"noreferrer noopener\" aria-label=\"Nishikawa's work (opens in a new tab)\">Nishikawa&#8217;s work<\/a>, we convert the Poisson equation for the gravitational potential into the <em>hyperbolic gravity equations<\/em>,<\/p>\n\n\n\n<p>$$\\frac{\\partial}{\\partial t}\\begin{bmatrix}                                                                    <br> \\phi\\\\[0.1cm]                                                                                       <br> q_1\\\\[0.1cm]                                                                                      <br> q_2\\\\[0.1cm]                                                                                      <br> \\end{bmatrix}<br>+<br> \\frac{\\partial}{\\partial x}\\begin{bmatrix}<br> -q_1\\\\[0.1cm]<br> -\\phi\/T_r\\\\[0.1cm]<br> 0\\\\[0.1cm]<br> \\end{bmatrix}<br>+<br> \\frac{\\partial}{\\partial y}\\begin{bmatrix}<br> -q_2\\\\[0.1cm]<br> 0\\\\[0.1cm]<br> -\\phi\/T_r\\\\[0.1cm]<br>\\end{bmatrix}<br>=<br>\\begin{bmatrix}                                                                                     <br> -4\\pi G \\rho\\\\[0.1cm]                                                                               <br> -q_1\/T_r\\\\[0.1cm]                                                                                 <br> -q_2\/T_r\\\\[0.1cm]                                                                                 <br> \\end{bmatrix},\\label{hyp}$$<\/p>\n\n\n\n<p>where the auxiliary variables $(q_1,q_2)^\\intercal\\approx\\vec{\\nabla}\\phi$ and $T_r$ is the relaxation time. The steady-state solution of the hyperbolic system is, in fact, the desired solution of the original Poisson problem.<\/p>\n\n\n\n<p>The flow and the gravity solvers operate on a joint hierarchical Cartesian mesh and are two-way coupled in each Runge-Kutta stage via the source terms. A key benefit of our strategy is that it allows the reuse of existing explicit hyperbolic solvers without modifications, while retaining their advanced features such as non-conforming and solution-adaptive grids.<\/p>\n\n\n\n<p>We implemented this purely hyperbolic approach for self-gravitating gas dynamics in <a rel=\"noreferrer noopener\" aria-label=\"Trixi.jl (opens in a new tab)\" href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" target=\"_blank\">Trixi.jl<\/a>, a tree-based numerical simulation framework for hyperbolic PDEs, which is written in Julia <a href=\"#ref-2\">[2]<\/a>. It was first validated by computing the Jeans gravitational instability, which demonstrates excellent agreement of the numerical results with the analytical solution:<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" loading=\"lazy\" width=\"976\" height=\"775\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-1.png\" alt=\"\" class=\"wp-image-5117\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-1.png 976w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-1-300x238.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-1-768x610.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-1-378x300.png 378w\" sizes=\"(max-width: 976px) 100vw, 976px\" \/><figcaption>Evolution of the computed kinetic, internal, and potential energies for the Jeans<br> instability using polynomial order N = 3 on a uniform 16 \u00d7 16 mesh. The analytical energies are included for reference.<\/figcaption><\/figure>\n\n\n\n<p>In a second example, we consider a modification of the Sedov blast wave problem that incorporates the effects of gravitational acceleration and which involves strong shocks and complex fluid interactions. We include it to demonstrate the shock capturing and AMR capabilities of Trixi.jl to resolve the cylindrical Sedov blast wave:<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" loading=\"lazy\" width=\"718\" height=\"1024\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-2-718x1024.png\" alt=\"\" class=\"wp-image-5118\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-2-718x1024.png 718w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-2-210x300.png 210w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-2-768x1096.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/09\/image-2.png 865w\" sizes=\"(max-width: 718px) 100vw, 718px\" \/><figcaption>Self-gravitating Sedov explosion approximated with polynomial order $N = 3$ on an AMR grid with four levels of possible refinement at $T = 0.5$ (left) and $T = 1.0$ (right). (first and third row) Two dimensional plots of the density and gravitational potential at two times. The white overlay of squares in the density plots shows the AMR grids used by both solvers. (second and last row) One dimensional slices of the solutions along a line from the origin to the edge of the domain in the positive x-direction.<\/figcaption><\/figure>\n\n\n\n<p>Finally, we show the evolution of the solution together with the adaptive mesh, which is dynamically refined and coarsened to maintain high accuracy while reducing the overall computation time:<\/p>\n\n\n\n<figure class=\"wp-block-image\"><img decoding=\"async\" src=\"https:\/\/raw.githubusercontent.com\/trixi-framework\/paper-self-gravitating-gas-dynamics\/master\/assets\/sedov-rho-phi-mesh.gif\" alt=\"https:\/\/raw.githubusercontent.com\/trixi-framework\/paper-self-gravitating-gas-dynamics\/master\/assets\/sedov-rho-phi-mesh.gif\"\/><figcaption>Simulation of self-gravitating Sedov blast wave with adaptive mesh refinement, where adaptation performed after every time step of the compressible Euler solver. The total run time compared to a uniformly refined mesh is reduced by a factor of 19, while the solutions are virtually unchanged.<\/figcaption><\/figure>\n\n\n\n<p>All results in the paper can be easily reproduced by using the resources provided in <a href=\"#ref-3\">[3]<\/a>, including step-by-step instructions for how to obtain and install Trixi.jl.<\/p>\n\n\n\n<h2 class=\"wp-block-heading\">References<\/h2>\n\n\n\n<p><a name=\"ref-1\"><\/a>[1] Schlottke-Lakemper, Michael and Winters, Andrew R and Ranocha, Hendrik and Gassner, Gregor J, <em>A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics<\/em>, submitted to J Comp Phys, 2020. <a href=\"https:\/\/arxiv.org\/abs\/2008.10593\">arXiv:2008.10593<\/a><\/p>\n\n\n\n<p><a name=\"ref-2\"><\/a>[2] Schlottke-Lakemper, Michael and Gassner, Gregor J and Ranocha, Hendrik and Winters, Andrew R, <em>Trixi.jl: A tree-based numerical simulation framework for hyperbolic PDEs written in Julia<\/em>, 2020. <a rel=\"noreferrer noopener\" aria-label=\"doi:10.5281\/zenodo.3996439 (opens in a new tab)\" href=\"https:\/\/doi.org\/10.5281\/zenodo.3996439\" target=\"_blank\">doi:10.5281\/zenodo.3996439<\/a><\/p>\n\n\n\n<p><a name=\"ref-3\"><\/a>[3] Schlottke-Lakemper, Michael and Winters, Andrew R and Ranocha, Hendrik and Gassner, Gregor J, <em>Self-gravitating gas dynamics simulations with Trixi.jl<\/em>, 2020. <a rel=\"noreferrer noopener\" aria-label=\"doi:10.5281\/zenodo.3996575 (opens in a new tab)\" href=\"https:\/\/doi.org\/10.5281\/zenodo.3996575\" target=\"_blank\">doi:10.5281\/zenodo.3996575<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/09\/01\/trixi-jl-a-tree-based-numerical-simulation-framework-for-hyperbolic-pdes-written-in-julia\/\">Trixi.jl: A tree-based numerical simulation framework for hyperbolic PDEs written in Julia<\/a> <span class=\"date\">September 1, 2020<\/span><div class=\"content\">\n<div class=\"wp-block-image\"><figure class=\"aligncenter\"><a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\/blob\/master\/docs\/src\/assets\/logo.png\" target=\"_blank\" rel=\"noreferrer noopener\"><img decoding=\"async\" src=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\/raw\/master\/docs\/src\/assets\/logo.png\" alt=\"\"\/><\/a><\/figure><\/div>\n\n\n\n<p><strong><a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\" target=\"_blank\" rel=\"noreferrer noopener\" aria-label=\"Trixi.jl (opens in a new tab)\">Trixi.jl<\/a><\/strong> is a numerical simulation framework for hyperbolic conservation laws written in <a href=\"https:\/\/julialang.org\">Julia<\/a>. A key objective for the framework is to be useful to both scientists and students. Therefore, next to having an extensible design with a fast implementation, Trixi is focused on being easy to use for new or inexperienced users, including the installation and postprocessing procedures. Its features include:<\/p>\n\n\n\n<ul><li>Hierarchical quadtree\/octree grid with adaptive mesh refinement<\/li><li>Native support for 2D and 3D simulations<\/li><li>High-order accuracy in space in time<\/li><li>Nodal discontinuous Galerkin spectral element methods <ul><li>Kinetic energy-preserving and entropy-stable split forms<\/li><li>Entropy-stable shock capturing <\/li><\/ul><\/li><li>Explicit low-storage Runge-Kutta time integration<\/li><li>Square\/cubic domains with periodic and Dirichlet boundary conditions<\/li><li>Multiple governing equations: <ul><li>Compressible Euler equations<\/li><li>Magnetohydrodynamics equations<\/li><li>Hyperbolic diffusion equations for elliptic problems<\/li><li>Scalar advection <\/li><\/ul><\/li><li>Multi-physics simulations <ul><li><a href=\"https:\/\/github.com\/trixi-framework\/paper-self-gravitating-gas-dynamics\">Self-gravitating gas dynamics<\/a> <\/li><\/ul><\/li><li>Shared-memory parallelization via multithreading<\/li><li>Visualization of results with Julia-only tools (2D) or ParaView\/VisIt (2D\/3D)<\/li><\/ul>\n\n\n\n<p>Trixi.jl was initiated by <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/schlottke-lakemper\">Michael Schlottke-Lakemper<\/a> and <a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/gregor-gassner\">Gregor Gassner<\/a> (both University of Cologne, Germany). Together with <a href=\"https:\/\/ranocha.de\">Hendrik Ranocha<\/a> (KAUST, Saudi Arabia) and <a href=\"https:\/\/liu.se\/en\/employee\/andwi94\">Andrew Winters<\/a> (Link\u00f6ping University, Sweden), they are the principal developers of Trixi.<\/p>\n\n\n\n<p>In case of questions, please feel free to <a href=\"https:\/\/github.com\/trixi-framework\/Trixi.jl\/issues\/new\">create an issue<\/a>. We are looking forward to feedback and\/or potential scientific collaboration.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/08\/26\/new-paper-submitted-a-purely-hyperbolic-discontinuous-galerkin-approach-for-self-gravitating-gas-dynamics\/\">New paper submitted: A purely hyperbolic discontinuous Galerkin approach for self-gravitating gas dynamics<\/a> <span class=\"date\">August 26, 2020<\/span><div class=\"content\">\n<p>One of the challenges when simulating astrophysical flows with self-gravity is to compute the gravitational forces. In contrast to the hyperbolic hydrodynamic equations, the gravity field is described by an elliptic Poisson equation. We present a purely hyperbolic approach by reformulating the elliptic problem into a hyperbolic diffusion problem, which is solved in pseudotime using the same explicit high-order discontinuous Galerkin method we use for the flow solution. The flow and the gravity solvers operate on a joint hierarchical Cartesian mesh and are two-way coupled via the source terms. A key benefit of our approach is that it allows the reuse of existing explicit hyperbolic solvers without modifications, while retaining their advanced features such as non-conforming and solution-adaptive grids. By updating the gravitational field in each Runge-Kutta stage of the hydrodynamics solver, high-order convergence is achieved even in coupled multi-physics simulations. After verifying the expected order of convergence for single-physics and multi-physics setups, we validate our approach by a simulation of the Jeans gravitational instability. Furthermore, we demonstrate the full capabilities of our numerical framework by computing a self-gravitating Sedov blast with shock capturing in the flow solver and adaptive mesh refinement for the entire coupled system.<\/p>\n\n\n\n<p>Available at <a href=\"https:\/\/arxiv.org\/abs\/2008.10593\">arXiv:2008.10593<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/08\/24\/snapshot-a-new-approach-for-approximating-conservation-laws\/\">Snapshot: A new approach for approximating conservation laws<\/a> <span class=\"date\">August 24, 2020<\/span><div class=\"content\"><p>A new approach for approximating conservation laws has been tested: instead of monitoring the changes of the means of a quantity within a certain volume over time like a finite volume procedure, this method mimics the behavior of simple solutions.<\/p>\n<p>First the current state is split into waves w, each aligned along one of the eigenvectors of the Fluxjacobian and living on its own grid. Then the corresponding grids are being moved with their eigenvelocity and finally all waves w are being overlayed and normalized yielding the solution after one timestep.<\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/initialcond.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/initialcond-300x225.png\" alt=\"\" width=\"300\" height=\"225\" class=\"alignnone size-medium wp-image-5061\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/initialcond-300x225.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/initialcond-400x300.png 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/initialcond.png 640w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/MyPlot.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/MyPlot-300x225.png\" alt=\"\" width=\"300\" height=\"225\" class=\"alignnone size-medium wp-image-5062\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/MyPlot-300x225.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/MyPlot-400x300.png 400w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/MyPlot.png 640w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>The plots show the solution of isothermal Euler equations at Time T=0.0000 (initial condition) and T=0.0025 with N=40 gridpoints for each mesh.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/07\/31\/snapshot-higher-order-schemes-for-the-mhd-equations\/\">Snapshot: Higher-order schemes for the MHD equations<\/a> <span class=\"date\">July 31, 2020<\/span><div class=\"content\"><p>A robust and easy way of simulating a hyperbolic test case with discontinuities is to use a 1st order finite volume scheme. Using such a method for magnetohydrodynamics (MHD) problems like the Orszag-Tang vortex leads to following results:<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5046\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1-300x292.png\" alt=\"\" width=\"300\" height=\"292\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1-300x292.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1-768x749.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1-308x300.png 308w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1-45x45.png 45w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/fv1.png 950w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><br \/>\nFor this and the following examples we used a 4th order time integration scheme and 256 degrees of freedom in each spatial direction.<\/p>\n<p>A way to generate more accurate results is to increase the order of the scheme, which has to be treated with caution near discontinuities because oscillations may occur. To overcome this issue one could use higher order schemes in smooth regions and lower order schemes in regions with discontinuities. An example for such an approach is our so called DGFV scheme, which blends e.g. a 4th order Discontinuous Galerkin scheme with a 1st order Finite Volume scheme.<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5049\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv-300x289.png\" alt=\"\" width=\"300\" height=\"289\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv-300x289.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv-768x740.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv-311x300.png 311w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/dgfv.png 946w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>Another way of doing it, is to use a suitable 4th order finite-volume scheme for MHD with a fitting limiter.<br \/>\n<a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25.png\"><img decoding=\"async\" loading=\"lazy\" class=\"alignnone size-medium wp-image-5051\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25-300x272.png\" alt=\"\" width=\"300\" height=\"272\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25-300x272.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25-768x697.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25-1024x929.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25-331x300.png 331w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/08\/Bildschirmfoto-2020-08-01-um-13.17.25.png 1170w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/07\/16\/new-paper-published-entropy-stable-p-nonconforming-discretizations-with-the-summation-by-parts-property-for-the-compressible-navier-stokes-equations\/\">New paper published: Entropy-Stable p-Nonconforming Discretizations with the Summation-by-Parts Property for the Compressible Navier\u2013Stokes Equations<\/a> <span class=\"date\">July 16, 2020<\/span><div class=\"content\"><p>The entropy-conservative\/stable, curvilinear, nonconforming, p-refinement algorithm for hyperbolic conservation laws of Del Rey Fern\u00e1ndez  et al. (2019) is extended from the compressible Euler equations to the compressible Navier\u2013Stokes equations. A simple and flexible coupling procedure with planar interpolation operators between adjoining nonconforming elements is used. Curvilinear volume metric terms are numerically approximated via a minimization procedure and satisfy the discrete geometric conservation law conditions. Distinct curvilinear surface metrics are used on the adjoining interfaces to construct the interface coupling terms, thereby localizing the discrete geometric conservation law constraints to each individual element. The resulting scheme is entropy conservative\/stable, element-wise conservative, and freestream preserving. Viscous interface dissipation operators that retain the entropy stability of the base scheme are developed. The accuracy and stability of the resulting numerical scheme are shown to be comparable to those of the original conforming scheme in Carpenter et al. (2014) and Parsani et al. (2016), i.e., this scheme achieves ~p+1\/2 convergence on geometrically high-order distorted element grids; this is demonstrated in the context of the viscous shock problem, the Taylor\u2013Green vortex problem at a Reynolds number of Re = 1,600 and a subsonic turbulent flow past a sphere at Re = 2,000.<\/p>\n<p><a href=\"https:\/\/www.sciencedirect.com\/science\/article\/abs\/pii\/S0045793020302036?via%3Dihub\">sciencedirect<\/a><br \/>\n<a href=\"https:\/\/arxiv.org\/abs\/1909.12546\">arxiv<\/a><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/07\/16\/new-article-published-in-eccomas-newsletter-p-16-20-split-form-discontinuous-galerkin-methods-for-implicit-large-eddy-simulation-of-compressible-turbulence\/\">New article published in ECCOMAS Newsletter (p.16-20): Split Form Discontinuous Galerkin Methods For Implicit Large Eddy Simulation Of Compressible Turbulence<\/a> <span class=\"date\">July 16, 2020<\/span><div class=\"content\"><p>As a teaser we show a numerical demonstration of the capabilities of the split form DG approach. We consider the flow past a plunging SD7003 airfoil at Mach number and cord length based Reynolds number Rec = 40,000. We subdivide the computational domain into 58,490 unstructured curvilinear hexahedral elements and use a polynomial degree N = 7 resulting in a total of about 150 million degrees of freedom. <\/p>\n<p><a href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA.png\"><img decoding=\"async\" loading=\"lazy\" src=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA-300x169.png\" alt=\"\" width=\"300\" height=\"169\" class=\"alignnone size-medium wp-image-5028\" srcset=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA-300x169.png 300w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA-768x432.png 768w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA-1024x576.png 1024w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA-500x281.png 500w, https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-content\/uploads\/2020\/07\/renderA.png 1035w\" sizes=\"(max-width: 300px) 100vw, 300px\" \/><\/a><\/p>\n<p>The Newsletter is available via: <a href=\"https:\/\/www.eccomas.org\/wp-content\/uploads\/sites\/15\/2020\/07\/Eccomas_Newsletter_July_2020_online.pdf\">eccomas Newsletter<\/a>.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/06\/26\/snapshot-variant-of-the-boss-bodenheimer-test-with-a-self-gravitating-high-order-hydro-code\/\">Snapshot: Variant of the Boss-Bodenheimer test with a self-gravitating high order hydro code<\/a> <span class=\"date\">June 26, 2020<\/span><div class=\"content\"><p>Variant of the classical Boss-Bodenheimer Test: We simulate the collapse and fragmentation of a rotating cloud core using a Discontinous Galerkin \/ Finite Volume hybrid high-order code. The self-gravity of the gas is calculated via the Fast-Multipole-Method.<\/p>\n<p><iframe loading=\"lazy\" title=\"Variant of the Boss-Bodenheimer test with a self-gravitating high order hydro code\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/sB2I1mKVh1k?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/05\/27\/snapshot-2d-pyopengl-real-time-simulation-of-the-forward-facing-step-test-case\/\">Snapshot: 2D PyOpenGL real-time simulation of the forward-facing step test case<\/a> <span class=\"date\">May 27, 2020<\/span><div class=\"content\"><p>WENO methods refers to a class of nonlinear finite volume or finite difference methods which are well suited to approximate solutions of hyperbolic conservation laws with high order accuracy in smooth regions and essentially non-oscillatory transition for solution discontinuities. In the following video we used a 5th order accurate WENO-FD scheme with immersed boundary conditions. <\/p>\n<p><iframe loading=\"lazy\" title=\"2D PyOpenGL real-time simulation of the forward-facing step test case\" width=\"584\" height=\"329\" src=\"https:\/\/www.youtube.com\/embed\/RQ3jDFw3MG0?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<p>Here we solve the forward-facing step test case in real-time using Python3\/PyOpenGL\/GLFW with an Intel Corporation HD Graphics 620 (rev 02). As one can see, it is also possible to add pressure peaks to the simulation by simply clicking with the mouse on the desired field.<\/p>\n<\/div><\/li><li class=\"listing-item\"><a class=\"title\" href=\"https:\/\/www.mi.uni-koeln.de\/NumSim\/2020\/04\/28\/snapshot-smoothed-particle-hydrodynamics-sph-simulation-of-the-numerical-simulation-research-group\/\">Snapshot: Smoothed-Particle Hydrodynamics (SPH) simulation of the Numerical Simulation Research Group<\/a> <span class=\"date\">April 28, 2020<\/span><div class=\"content\"><p>Smoothed-particle hydrodynamics (SPH) is a computational method used for simulating the mechanics of continuum media, such as solid mechanics and is being increasingly used to model fluid motion as well.<br \/>\nThis method is well-suited for problems dominated by complex boundary dynamics, since SPH is a mesh-free method, as well as for mass conservation problems since the particles themselves represent mass.<\/p>\n<p><iframe loading=\"lazy\" title=\"Smoothed-Particle Hydrodynamics (SPH) simulation of the Numerical Simulation Research Group\" width=\"584\" height=\"438\" src=\"https:\/\/www.youtube.com\/embed\/L3qoXcMQZms?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; encrypted-media; gyroscope; picture-in-picture\" allowfullscreen><\/iframe><\/p>\n<p>We used a 2D Python\/OpenCL SPH code that solves the incompressible Navier-Stokes equations in real time. In this simulation, 52,700 particles were used to smash the people of the Numerical Simulation Research Group against each other.<\/p>\n<\/div><\/li><\/ul>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":9,"featured_media":0,"parent":0,"menu_order":8,"comment_status":"closed","ping_status":"closed","template":"","meta":{"footnotes":""},"post_mailing_queue_ids":[],"_links":{"self":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/pages\/1656"}],"collection":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/users\/9"}],"replies":[{"embeddable":true,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/comments?post=1656"}],"version-history":[{"count":104,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/pages\/1656\/revisions"}],"predecessor-version":[{"id":6112,"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/pages\/1656\/revisions\/6112"}],"wp:attachment":[{"href":"https:\/\/www.mi.uni-koeln.de\/NumSim\/wp-json\/wp\/v2\/media?parent=1656"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}