| 1. |
|
|
| 2. |
- Birken, Philipp, et al.
(författare)
-
A Note on Terminology in Multigrid Methods
- 2016
-
Ingår i: PAMM - Proceedings in Applied Mathematics and Mechanics. - : Wiley. - 1617-7061. ; 16, s. 721-722
-
Tidskriftsartikel (refereegranskat)abstract
- We compare terminology used in the literature on multigrid methods for compressible computational fluid dynamics to that used in linear multigrid theory. Several popular iterative and direct smoothers are presented side-by-side using the same terminology. We argue for greater analysis of these methods in order to place them into a more rigorous framework and to identify the most promising candidates for future development.
|
|
| 3. |
- Birken, Philipp, et al.
(författare)
-
A STUDY of multigrid smoothers used in compressible CFD based on the convection diffusion equation
- 2016
-
Ingår i: ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering. - 9786188284401 ; 2, s. 2648-2663
-
Konferensbidrag (refereegranskat)abstract
- We look at multigrid methods for unsteady viscous compressible flows. We specifically target smoothers that can be used in parallel and without computation of a Jacobian, which are particlarly attractive candidates in the context of Discontinuous Galerkin discretizations. In CFD, a plethora of nonlinear smoothers have been suggested which are hard to analyze. Our methodology is to use a linear model problem, here the convection diffusion equation, to be able to classify and compare smoothers better. Specifically, we consider explicit and implicit pseudo time iterations, GMRES as a smoother, SGS and implicit line smoothers. We relate GMRES to explicit Runge-Kutta smoothers, identify implicit line smoothers as Block Jacobi and analyze the potential of implicit pseudo time iterations. Finally, we discuss the relation between methods for steady and unsteady flows. Numerical results show that GMRES is a very attractive smoother in this context.
|
|
| 4. |
|
|
| 5. |
- Birken, Philipp, et al.
(författare)
-
Conservation Properties of Iterative Methods for Implicit Discretizations of Conservation Laws
- 2022
-
Ingår i: Journal of Scientific Computing. - : Springer Science and Business Media LLC. - 0885-7474 .- 1573-7691. ; 92:2
-
Tidskriftsartikel (refereegranskat)abstract
- Conservation properties of iterative methods applied to implicit finite volume discretizations of nonlinear conservation laws are analyzed. It is shown that any consistent multistep or Runge-Kutta method is globally conservative. Further, it is shown that Newton’s method, Krylov subspace methods and pseudo-time iterations are globally conservative while the Jacobi and Gauss-Seidel methods are not in general. If pseudo-time iterations using an explicit Runge-Kutta method are applied to a locally conservative discretization, then the resulting scheme is also locally conservative. However, the corresponding numerical flux can be inconsistent with the conservation law. We prove an extension of the Lax-Wendroff theorem, which reveals that numerical solutions based on these methods converge to weak solutions of a modified conservation law where the flux function is multiplied by a particular constant. This constant depends on the choice of Runge-Kutta method but is independent of both the conservation law and the discretization. Consistency is maintained by ensuring that this constant equals unity and a strategy for achieving this is presented. Simulations show that this strategy improves the convergence rate of the pseudo-time iterations. Experiments with GMRES suggest that it also suffers from inconsistency but that this is automatically accounted for after some finite number of iterations. Similar experiments with coarse grid corrections based on agglomeration indicate no inconsistency.
|
|
| 6. |
- Birken, Philipp, et al.
(författare)
-
Extrapolation in Time in Thermal Fluid Structure Interaction
- 2015
-
Ingår i: Recent Trends in Computational Engineering - CE2014 : Optimization, Uncertainty, Parallel Algorithms, Coupled and Complex Problems (Lecture Notes in Computational Science and Engineering). - Cham : Springer International Publishing. ; 105, s. 215-231
-
Konferensbidrag (refereegranskat)abstract
- We consider time dependent thermal fluid structure interaction. The respective models are the compressible Navier-Stokes equations and the nonlinear heat equation. A partitioned coupling approach via a Dirichlet-Neumann method and a fixed point iteration is employed. As a reference solver a previously developed efficient time adaptive higher order time integration scheme is used. To improve upon this, we work on reducing the number of fixed point coupling iterations. Thus, we explore the idea of extrapolation based on data given from the time integration and derive such methods for SDIRK2. This allows to reduce the number of fixed point iterations further by up to a factor of two with linear extrapolation performing better than quadratic.
|
|
| 7. |
- Birken, Philipp, et al.
(författare)
-
Fast Solvers for Unsteady Thermal Fluid Structure Interaction
- 2015
-
Ingår i: International Journal for Numerical Methods in Fluids. - : Wiley. - 1097-0363 .- 0271-2091. ; 79:1, s. 16-29
-
Tidskriftsartikel (refereegranskat)abstract
- We consider time-dependent thermal fluid structure interaction. The respective models are the compressible Navier–Stokes equations and the nonlinear heat equation. A partitioned coupling approach via a Dirichlet–Neumann method and a fixed point iteration is employed. As a reference solver, a previously developed efficient time-adaptive higher-order time integration scheme is used. To improve on this, we work on reducing the number of fixed point coupling iterations. Using the idea of extrapolation based on data given from the time integration by deriving such methods for SDIRK2, it is possible to reduce the number of fixed point iterations further by up to a factor of two with linear extrapolation performing better than quadratic. This leads to schemes that can use less than two iterations per time step. Furthermore, widely used vector extrapolation methods for convergence acceleration of the fixed point iteration are tested, namely Aitken relaxation, minimal polynomial extrapolation and reduced rank extrapolation. These have no beneficial effects.
|
|
| 8. |
|
|
| 9. |
- Birken, Philipp, et al.
(författare)
-
Numerical methods for unsteady thermal fluid structure interaction
- 2017
-
Ingår i: Fluid-Structure Interaction Modeling, Adaptive Discretisations and Solvers. - : De Gruyter. - 9783110494259 ; , s. 129-168
-
Bokkapitel (refereegranskat)abstract
- We discuss thermal fluid-structure interaction processes, where a simulation of the time-dependent temperature field is of interest. Thereby, we consider partitioned coupling schemes with a Dirichlet-Neumann method. We present an analysis of the method on a model problem of discretized coupled linear heat equations. This shows that for large quotients in the heat conductivities, the convergence rate will be very small. The time dependencymakes the use of time-adaptive implicitmethods imperative. This gives rise to the question as to how accurately the appearing nonlinear systems should be solved, which is discussed in detail for both the nonlinear and linear case. The efficiency of the resulting method is demonstrated using realistic test cases.
|
|
| 10. |
- Birken, Philipp, et al.
(författare)
-
On stability and conservation properties of (S)epirk integrators in the context of discretized pdes
- 2018
-
Ingår i: Theory, Numerics and Applications of Hyperbolic Problems II. - Cham : Springer International Publishing. - 9783319915470 ; 237, s. 617-629
-
Konferensbidrag (refereegranskat)abstract
- Exponential integrators are becoming increasingly popular for stiff problems of high dimension due to their attractive property of solving the linear part of the system exactly and hence being A-stable. In practice, however, exponential integrators are implemented using approximation techniques to matrix-vector products involving functions of the matrix exponential (the so-called ϕ-functions) to make them efficient and competitive to other state-of-the-art schemes. We will examine linear stability and provide a Courant–Friedrichs–Lewy (CFL) condition of special classes of exponential integrator schemes called EPIRK and sEPIRK and demonstrate their dependence on the parameters of the embedded approximation technique. Furthermore, a conservation property of the EPIRK schemes is proven.
|
|
| 11. |
- Birken, Philipp, et al.
(författare)
-
Preconditioned Smoothers for the Full Approximation Scheme for the RANS Equations
- 2019
-
Ingår i: Journal of Scientific Computing. - : Springer Science and Business Media LLC. - 0885-7474 .- 1573-7691. ; 78:2, s. 995-1022
-
Tidskriftsartikel (refereegranskat)abstract
- We consider multigrid methods for finite volume discretizations of the Reynolds averaged Navier–Stokes equations for both steady and unsteady flows. We analyze the effect of different smoothers based on pseudo time iterations, such as explicit and additive Runge–Kutta (AERK) methods. Furthermore, by deriving them from Rosenbrock smoothers, we identify some existing schemes as a class of additive W (AW) methods. This gives rise to two classes of preconditioned smoothers, preconditioned AERK and AW, which are implemented the exact same way, but have different parameters and properties. This derivation allows to choose some of these based on results for time integration methods. As preconditioners, we consider SGS preconditioners based on flux vector splitting discretizations with a cutoff function for small eigenvalues. We compare these methods based on a discrete Fourier analysis. Numerical results on pitching and plunging airfoils identify AW3 as the best smoother regarding overall efficiency. Specifically, for the NACA 64A010 airfoil steady-state convergence rates of as low as 0.85 were achieved, or a reduction of 6 orders of magnitude in approximately 25 pseudo-time iterations. Unsteady convergence rates of as low as 0.77 were achieved, or a reduction of 11 orders of magnitude in approximately 70 pseudo-time iterations.
|
|
| 12. |
- Birken, Philipp, et al.
(författare)
-
Subcell finite volume multigrid preconditioning for high-order discontinuous Galerkin methods
- 2019
-
Ingår i: International Journal of Computational Fluid Dynamics. - : Informa UK Limited. - 1061-8562 .- 1029-0257. ; 33:9, s. 353-361
-
Tidskriftsartikel (refereegranskat)abstract
- We suggest a new multigrid preconditioning strategy for use in Jacobian-free Newton–Krylov (JFNK) methods for the solution of algebraic equation systems arising from implicit Discontinuous Galerkin (DG) discretisations. To define the new preconditioner, use is made of an auxiliary first-order finite volume discretisation that refines the original DG mesh, but can still be implemented algebraically. As smoother, we consider the pseudo-time iteration W3 with a symmetric Gauss–Seidel-type approximation of the Jacobian. As a proof of concept numerical tests are presented for the one-dimensional Euler equations, demonstrating the potential of the new approach.
|
|
| 13. |
- Birken, Philipp
(författare)
-
Termination criteria for inexact fixed-point schemes
- 2015
-
Ingår i: Numerical Linear Algebra with Applications. - : Wiley. - 1070-5325. ; 22:4, s. 702-716
-
Tidskriftsartikel (refereegranskat)abstract
- We analyze inexact fixed-point iterations where the generating function contains an inexact solve of an equation system to answer the question of how tolerances for the inner solves influence the iteration error of the outer fixed-point iteration. Important applications are the Picard iteration and partitioned fluid-structure interaction. For the analysis, the iteration is modeled as a perturbed fixed-point iteration, and existing analysis is extended to the nested case x=F(S(x)). We prove that if the iteration converges, it converges to the exact solution irrespective of the tolerance in the inner systems, provided that a nonstandard relative termination criterion is employed, whereas standard relative and absolute criteria do not have this property. Numerical results demonstrate the effectiveness of the approach with the nonstandard termination criterion.
|
|
| 14. |
- Blom, David S., et al.
(författare)
-
A comparison of Rosenbrock and ESDIRK methods combined with iterative solvers for unsteady compressible flows
- 2016
-
Ingår i: Advances in Computational Mathematics. - : Springer Science and Business Media LLC. - 1572-9044 .- 1019-7168. ; 42:6, s. 1401-1426
-
Tidskriftsartikel (refereegranskat)abstract
- In this article, we endeavour to find a fast solver for finite volume discretizations for compressible unsteady viscous flows. Thereby, we concentrate on comparing the efficiency of important classes of time integration schemes, namely time adaptive Rosenbrock, singly diagonally implicit (SDIRK) and explicit first stage singly diagonally implicit Runge-Kutta (ESDIRK) methods. To make the comparison fair, efficient equation system solvers need to be chosen and a smart choice of tolerances is needed. This is determined from the tolerance TOL that steers time adaptivity. For implicit Runge-Kutta methods, the solver is given by preconditioned inexact Jacobian-free Newton-Krylov (JFNK) and for Rosenbrock, it is preconditioned Jacobian-free GMRES. To specify the tolerances in there, we suggest a simple strategy of using TOL/100 that is a good compromise between stability and computational effort. Numerical experiments for different test cases show that the fourth order Rosenbrock method RODASP and the fourth order ESDIRK method ESDIRK4 are best for fine tolerances, with RODASP being the most robust scheme.
|
|
| 15. |
- Frenander, Hannes
(författare)
-
High-order finite difference approximations for hyperbolic problems : multiple penalties and non-reflecting boundary conditions
- 2017
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- In this thesis, we use finite difference operators with the Summation-By-Partsproperty (SBP) and a weak boundary treatment, known as SimultaneousApproximation Terms (SAT), to construct high-order accurate numerical schemes.The SBP property and the SAT’s makes the schemes provably stable. The numerical procedure is general, and can be applied to most problems, but we focus on hyperbolic problems such as the shallow water, Euler and wave equations.For a well-posed problem and a stable numerical scheme, data must be available at the boundaries of the domain. However, there are many scenarios where additional information is available inside the computational domain. In termsof well-posedness and stability, the additional information is redundant, but it can still be used to improve the performance of the numerical scheme. As a first contribution, we introduce a procedure for implementing additional data using SAT’s; we call the procedure the Multiple Penalty Technique (MPT).A stable and accurate scheme augmented with the MPT remains stable and accurate. Moreover, the MPT introduces free parameters that can be used to increase the accuracy, construct absorbing boundary layers, increase the rate of convergence and control the error growth in time.To model infinite physical domains, one need transparent artificial boundary conditions, often referred to as Non-Reflecting Boundary Conditions (NRBC). In general, constructing and implementing such boundary conditions is a difficult task that often requires various approximations of the frequency and range of incident angles of the incoming waves. In the second contribution of this thesis,we show how to construct NRBC’s by using SBP operators in time.In the final contribution of this thesis, we investigate long time error bounds for the wave equation on second order form. Upper bounds for the spatial and temporal derivatives of the error can be obtained, but not for the actual error. The theoretical results indicate that the error grows linearly in time. However, the numerical experiments show that the error is in fact bounded, and consequently that the derived error bounds are probably suboptimal.
|
|
| 16. |
- Gleim, Tobias, et al.
(författare)
-
Experimental and Numerical Aspects of a Thermal Fluid-Structure Phenomenon
- 2017
-
Ingår i: International Conference of Numerical Analysis and Applied Mathematics (ICNAAM 2016). - : Author(s). - 0094-243X .- 1551-7616. - 9780735415386 ; 1863
-
Konferensbidrag (refereegranskat)abstract
- A fundamental research experiment for thermal fluid-structure-interaction for the verification of a partitioned approach with non-linear material properties is examined. In the following, a specimen is heated as well as cooled within a wind tunnel. The thermal fluid-structure-interaction is first experimentally investigated and subsequently numerically validated. For the numerical simulation, two existing programs (a fluid and a structure code) are coupled using a partitioned approach.
|
|
| 17. |
- Gleim, Tobias, et al.
(författare)
-
Thermal Fluid-Structure-Interaction - Experimental and Numerical Analysis
- 2015
-
Ingår i: Proceedings of the International Conference of Numerical Analysis and Applied Mathematics 2014 (ICNAAM-2014). - : AIP Publishing LLC. - 0094-243X .- 1551-7616. ; 1648
-
Konferensbidrag (refereegranskat)abstract
- In the present paper the thermal fluid-structure-interaction is experimentally and numerically investigated. Therefore, the interaction phenomena is modeled by the Reynolds-averaged Navier-Stokes equations and the nonlinear Fourier heat conduction equation are used for the fluid and the solid phase, respectively. The simulation is performed using a partitioned approach using the finite volume method for the fluid domain, the finite element method for the solid domain and Runge-Kutta integration schemes for the time domain. Furthermore, as a basis for the understanding of thermal fluid-structure-interaction and also for the veri cation and validation of the applied continuum mechanical models and numerical methods, respectively, a fundamental wind tunnel experiment is presented.
|
|
| 18. |
- Görtz, Morgan, et al.
(författare)
-
On the convergence rate of the dirichlet-neumann iteration for coupled poisson problems on unstructured grids
- 2020
-
Ingår i: Finite Volumes for Complex Applications IX - Methods, Theoretical Aspects, Examples, FVCA 2020. - Cham : Springer International Publishing. - 2194-1009 .- 2194-1017. - 9783030436506 ; 323, s. 355-363
-
Konferensbidrag (refereegranskat)abstract
- We consider thermal fluid structure interaction with a partitioned approach, where typically, a finite volume and a finite element code would be coupled. As a model problem, we consider two coupled Poisson problems with heat conductivities $$\lambda _1$$, $$\lambda _2$$ in one dimension on intervals of length $$l:1$$ and $$l:2$$. Hereby, we consider linear discretizations on arbitrary meshes, such as finite volumes, finite differences, finite elements. For these, we prove that the convergence rate of the Dirichlet-Neumann iteration is given by $$\lambda _1l_2/\lambda _2l_1$$ and is thus independent of discretization and mesh.
|
|
| 19. |
- Kasimir, Johannes, et al.
(författare)
-
An finite volume based multigrid preconditioner for dg-sem for convection-diffusion
- 2021
-
Ingår i: Fluid Dynamics and Transport Phenomena. - : CIMNE. ; 600, s. 1-12
-
Konferensbidrag (refereegranskat)abstract
- The goal of our research is the construction of efficient Jacobian-free preconditioners for high order Discontinuous Galerkin (DG) discretizations with implicit time integration. We are motivated by three-dimensional unsteady compressible flow applications, which often result in large stiff systems. Implicit time integrators overcome the impact upon restrictive CFL conditions on explicit ones but leave the problem to solve huge nonlinear systems. In this paper we consider a multigrid preconditioning strategy for Jacobian-free Newton-Krylov (JFNK) methods for the solution of algebraic equation systems arising from implicit Discontinuous Galerkin (DG) discretizations. The preconditioner is defined by an auxiliary first order Finite Volume (FV) discretization that refines the original DG mesh, but can still be implemented algebraically. Different options exist to define the grid transfer between DG and FV. We suggest an ad hoc assignment of the unknowns as well as L2 projections. We present new numerical results for the two-dimensional convection-diffusion equation in combination with the different transfer options, which demonstrate the quality and efficiency of the suggested preconditioner with regards to convergence speed up and CPU time. The suggested L2 projection from this paper result in the best convergence speed up.
|
|
| 20. |
- Linders, Viktor, et al.
(författare)
-
Locally conservative and flux consistent iterative methods
- 2023
-
Annan publikation (övrigt vetenskapligt/konstnärligt)abstract
- Conservation and consistency are fundamental properties of discretizations of systems of hyperbolic conservation laws. Here, these concepts are extended to the realm of iterative methods by formally defining locally conservative and flux consistent iterations. These concepts are of both theoretical and practical importance: Based on recent work by the authors, it is shown that pseudo-time iterations using explicit Runge-Kutta methods are locally conservative but not necessarily flux consistent. An extension of the Lax-Wendroff theorem is presented, revealing convergence towards weak solutions of a temporally retarded system of conservation laws. Each equation is modified in the same way, namely by a particular scalar factor multiplying the spatial flux terms. A technique for enforcing flux consistency, and thereby recovering convergence, is presented. Further, local conservation is established for all Krylov subspace methods, with and without restarts, and for Newton's method under certain assumptions on the discretization. Thus it is shown that Newton-Krylov methods are locally conservative, although not necessarily flux consistent. Numerical experiments with the 2D compressible Euler equations corroborate the theoretical results. Further numerical investigations of the impact of flux consistency on Newton-Krylov methods indicate that its effect is case dependent, and diminishes as the number of iterations grow.
|
|
| 21. |
- Linders, Viktor, et al.
(författare)
-
LOCALLY CONSERVATIVE AND FLUX CONSISTENT ITERATIVE METHODS
- 2024
-
Ingår i: SIAM Journal on Scientific Computing. - 1064-8275. ; 46:2, s. 424-444
-
Tidskriftsartikel (refereegranskat)abstract
- Conservation and consistency are fundamental properties of discretizations of conservation laws, necessary to ensure physically meaningful solutions. In the context of systems of nonlinear hyperbolic conservation laws, conservation and consistency additionally play an important role in convergence theory via the Lax-Wendroff theorem. Here, these concepts are extended to the realm of iterative methods by formally defining locally conservative and flux consistent iterations. These concepts are used to prove an extension of the Lax-Wendroff theorem incorporating pseudotime iterations with explicit Runge-Kutta methods. This result reveals that lack of flux consistency implies convergence towards weak solutions of a time dilated system of conservation laws, where each equation is modified by a particular scalar factor multiplying the spatial flux terms. Local conservation is further established for Krylov subspace methods with and without restarts, and for Newton's method under certain assumptions on the discretization. It is thus shown that Newton-Krylov methods are locally conservative, although not necessarily flux consistent. Numerical experiments with the 2 dimensional compressible Euler equations corroborate the theoretical results. A simple technique for enforcing flux consistency of Newton-Krylov methods is presented. Experiments indicate that its efficacy is case dependent, and diminishes as the number of iterations grow.
|
|
| 22. |
- Linders, Viktor, et al.
(författare)
-
On the Consistency of Arnoldi-Based Krylov Methods for Conservation Laws
- 2023
-
Ingår i: PAMM - Proceedings in Applied Mathematics and Mechanics. - 1617-7061. ; 23:1
-
Tidskriftsartikel (refereegranskat)abstract
- Conservation and consistency are fundamental properties of discretizations of systems of hyperbolic conservation laws. Re- cently, these concepts have been extended to the realm of iterative methods by defining locally conservative and flux consistent iterations. In this note, the current status of such iterative methods is summarized. In particular, it has been shown that Krylov subspace methods are locally conservative, but that they are not flux consistent. Here, we approach the problem of quantifying the flux inconsistency of Krylov subspace methods. Krylov methods introduce a time retardation factor into discretizations of linear conservation laws. It has thusfar been unknown how to compute the precise value of this factor. This issue is resolved herein for Arnoldi-based Krylov subspace methods.
|
|
| 23. |
- Linders, Viktor, et al.
(författare)
-
Resolving entropy growth from iterative methods
- 2023
-
Ingår i: BIT Numerical Mathematics. - 0006-3835. ; 63:4
-
Tidskriftsartikel (refereegranskat)abstract
- We consider entropy conservative and dissipative discretizations of nonlinear conservation laws with implicit time discretizations and investigate the influence of iterative methods used to solve the arising nonlinear equations. We show that Newton’s method can turn an entropy dissipative scheme into an anti-dissipative one, even when the iteration error is smaller than the time integration error. We explore several remedies, of which the most performant is a relaxation technique, originally designed to fix entropy errors in time integration methods. Thus, relaxation works well in consort with iterative solvers, provided that the iteration errors are on the order of the time integration method. To corroborate our findings, we consider Burgers’ equation and nonlinear dispersive wave equations. We find that entropy conservation results in more accurate numerical solutions than non-conservative schemes, even when the tolerance is an order of magnitude larger.
|
|
| 24. |
- Linders, Viktor, et al.
(författare)
-
Resolving Entropy Growth from Iterative Methods
- 2023
-
Annan publikation (övrigt vetenskapligt/konstnärligt)abstract
- We consider entropy conservative and dissipative discretizations of nonlinear conservation laws with implicit time discretizations and investigate the influence of iterative methods used to solve the arising nonlinear equations. We show that Newton's method can turn an entropy dissipative scheme into an anti-dissipative one, even when the iteration error is smaller than the time integration error. We explore several remedies, of which the most performant is a relaxation technique, originally designed to fix entropy errors in time integration methods. Thus, relaxation works well in consort with iterative solvers, provided that the iteration errors are on the order of the time integration method. To corroborate our findings, we consider Burgers' equation and nonlinear dispersive wave equations. We find that entropy conservation results in more accurate numerical solutions than non-conservative schemes, even when the tolerance is an order of magnitude larger.
|
|
| 25. |
- Meisrimel, Peter, et al.
(författare)
-
A time adaptive multirate Dirichlet–Neumann waveform relaxation method for heterogeneous coupled heat equations
- 2023
-
Ingår i: ZAMM Zeitschrift fur Angewandte Mathematik und Mechanik. - 0044-2267. ; 103:11
-
Tidskriftsartikel (refereegranskat)abstract
- We consider partitioned time integration for heterogeneous coupled heat equations. First and second order multirate, as well as time-adaptive Dirichlet-Neumann Waveform relaxation (DNWR) methods are derived. In 1D and for implicit Euler time integration, we analytically determine optimal relaxation parameters for the fully discrete scheme. We test the robustness of the relaxation parameters on the second order multirate method in 2D. DNWR is shown to be very robust and consistently yielding fast convergence rates, whereas the closely related Neumann-Neumann Waveform relaxtion (NNWR) method is slower or even diverges. The waveform approach naturally allows for different timesteps in the subproblems. In a performance comparison for DNWR, the time-adaptive method dominates the multirate method due to automatically finding suitable stepsize ratios. Overall, we obtain a fast, robust, multirate and time adaptive partitioned solver for unsteady conjugate heat transfer.
|
|
| 26. |
- Meisrimel, Peter, et al.
(författare)
-
Goal Oriented Time Adaptivity Using Local Error Estimates
- 2020
-
Ingår i: Algorithms. - : MDPI AG. - 1999-4893. ; 13.5:113
-
Tidskriftsartikel (refereegranskat)abstract
- We consider initial value problems (IVPs) where we are interested in a quantity of interest(QoI) that is the integral in time of a functional of the solution. For these, we analyze goal orientedtime adaptive methods that use only local error estimates. A local error estimate and timestepcontroller for step-wise contributions to the QoI are derived. We prove convergence of the error in theQoI for tolerance to zero under a controllability assumption. By analyzing global error propagationwith respect to the QoI, we can identify possible issues and make performance predictions. Numericaltests verify these results. We compare performance with classical local error based time-adaptivityand a posteriori based adaptivity using the dual-weighted residual (DWR) method. For dissipativeproblems, local error based methods show better performance than DWR and the goal orientedmethod shows good results in most examples, with significant speedups in some cases.
|
|
| 27. |
- Meisrimel, Peter, et al.
(författare)
-
On Goal Oriented Time Adaptivity using Local Error Estimates
- 2017. - 1
-
Ingår i: PAMM : Proceedings in applied mathematics and mechanics - Proceedings in applied mathematics and mechanics. - : Wiley. - 1617-7061. ; 17, s. 849-850
-
Konferensbidrag (refereegranskat)abstract
- We consider adaptive time discretization methods for ordinary differential equations where one aims to control the error in a quantity of interest of the form J(u) = ∫ j(u(t))dt with j : Rd -> R. In this setting we propose a new timestep controller based on local error estimates of the quantity of interest. The new method converges when the tolerance goes to zero.We experimentally compare the new scheme with the classic norm-based time-adaptivity based on local error estimates as well as the dual-weighted residual (DWR) method. The results show significantly lower efficiency for the DWR method. The local error based schemes are similarly efficient, with the new scheme showing significant improvement in some cases.
|
|
| 28. |
- Meisrimel, Peter, et al.
(författare)
-
Waveform Iteration with asynchronous communication
- 2019
-
Ingår i: Special Issue: 90th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM). - 1617-7061. ; 19:1
-
Konferensbidrag (refereegranskat)abstract
- Consider the coupling of two multi-physics systems of time-dependent ODEs. We propose a new Waveform iteration type method for coupling which uses asynchronous communication to be both parallel in time and fast in convergence. Analytical results show convergence in the continuous setting and numerical results for two coupled head equations show good performance of this method, which is both parallel and converging fast.
|
|
| 29. |
- Meisrimel, Peter, et al.
(författare)
-
Waveform Relaxation with Asynchronous Time-integration
- 2022
-
Ingår i: ACM Transactions on Mathematical Software. - : Association for Computing Machinery (ACM). - 0098-3500 .- 1557-7295. ; 48:4
-
Tidskriftsartikel (refereegranskat)abstract
- We consider Waveform Relaxation (WR) methods for parallel and partitioned time-integration of surface-coupled multiphysics problems. WR allows independent time-discretizations on independent and adaptive time-grids, while maintaining high time-integration orders. Classical WR methods such as Jacobi or Gauss-Seidel WR are typically either parallel or converge quickly.We present a novel parallel WR method utilizing asynchronous communication techniques to get both properties. Classical WR methods exchange discrete functions after time-integration of a subproblem. We instead asynchronously exchange time-point solutions during time-integration and directly incorporate all new information in the interpolants. We show both continuous and time-discrete convergence in a framework that generalizes existing linear WR convergence theory. An algorithm for choosing optimal relaxation in our new WR method is presented. Convergence is demonstrated in two conjugate heat transfer examples. Our new method shows an improved performance over classical WR methods. In one example, we show a partitioned coupling of the compressible Euler equations with a nonlinear heat equation, with subproblems implemented using the open source libraries DUNE and FEniCS.
|
|
| 30. |
- Monge, Azahar, et al.
(författare)
-
A multirate Neumann-Neumann waveform relaxation method for heterogeneous coupled heat equations
- 2019
-
Ingår i: SIAM Journal on Scientific Computing. - 1064-8275. ; 41:5, s. 86-105
-
Tidskriftsartikel (refereegranskat)abstract
- An important challenge when coupling two different time dependent problems is to increase parallelization in time. We suggest a multirate Neumann-Neumann waveform relaxation algorithm to solve two heterogeneous coupled heat equations. In order to fix the mismatch produced by the multirate feature at the space-time interface a linear interpolation is constructed. The heat equations are discretized using a finite element method in space, whereas two alternative time integration methods are used: implicit Euler and SDIRK2. We perform a one-dimensional convergence analysis for the nonmultirate fully discretized heat equations using implicit Euler to find the optimal relaxation parameter in terms of the material coefficients, the step size, and the mesh resolution. This gives a very efficient method which needs only two iterations. Numerical results confirm the analysis and show that the one-dimensional nonmultirate optimal relaxation parameter is a very good estimator for the multirate one-dimensional case and even for multirate and nonmultirate two-dimensional examples using both implicit Euler and SDIRK2.
|
|
| 31. |
- Monge, Azahar, et al.
(författare)
-
A time‐adaptive Dirichlet‐Neumann waveform relaxation method for coupled heterogeneous heat equations
- 2019
-
Ingår i: Special Issue: 90th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM). - : Wiley. - 1617-7061. ; 19:1
-
Konferensbidrag (refereegranskat)abstract
- We introduce a time adaptive multirate method based on the Dirichlet-Neumann waveform relaxation (DNWR) algorithm for the simulation of two coupled linear heat equations with strong jumps in the material coefficients across the interface. Numerical results are included to illustrate the advantages of the time adaptive approach over the multirate approach and the robustness of the multirate DNWR method with respect to its sibling, the multirate Neumann-Neumann waveform relaxation (NNWR) method introduced in a previous work [3].
|
|
| 32. |
- Monge, Azahar, et al.
(författare)
-
Convergence analysis of coupling iterations for the unsteady transmission problem with mixed discretizations
- 2016
-
Ingår i: ECCOMAS Congress 2016 - Proceedings of the 7th European Congress on Computational Methods in Applied Sciences and Engineering. - 9786188284401 ; 1, s. 1530-1544
-
Konferensbidrag (refereegranskat)abstract
- We analyze the convergence rate of the Dirichlet-Neumann iteration for the fully discretized one dimensional unsteady transmission problem. Specifically, we consider the coupling of two linear heat equations on two identical non overlapping intervals. The Laplacian is discretized using finite differences on one interval and finite elements on the other and the implicit Euler method is used for the time discretization. Following previous analysis where finite elements where used on both subdomains, we provide an exact formula for the spectral radius of the iteration matrix for this specific mixed discretizations. We then show that these tend to the ratio of heat conductivities in the semidiscrete spatial limit, but to a factor of the ratio of the products of density and specific heat capacity in the semidiscrete temporal one. In the previous finite element analysis, the same result was obtained in the semidiscrete spatial limit but the factor in the temporal limit was lower. This explains the fast convergence previously observed for cases with strong jumps in the material coefficients. Numerical results confirm the analysis.
|
|
| 33. |
- Monge, Azahar, et al.
(författare)
-
Convergence Analysis of the Dirichlet-Neumann Iteration for Finite Element Discretizations
- 2016
-
Ingår i: PAMM - Proceedings in Applied Mathematics and Mechanics. - : Wiley. - 1617-7061. ; 16, s. 733-734
-
Tidskriftsartikel (refereegranskat)abstract
- We analyze the convergence rate of the Dirichlet-Neumann iteration for the fully discretized unsteady transmission problem. Specifically, we consider the coupling of two linear heat equations on two identical non overlapping domains with jumps in the material coefficients across these. In this context, we derive the iteration matrix of the coupled problem. In the 1D case, the spectral radius of the iteration matrix tends to the ratio of heat conductivities in the semidiscrete spatial limit, but to the ratio of the products of density and specific heat capacity in the semidiscrete temporal one. This explains the fast convergence previously observed for cases with strong jumps in the material coefficients.
|
|
| 34. |
- Monge, Azahar, et al.
(författare)
-
Convergence Speed of Coupling Iterations for the Unsteady Transmission Problem
- 2015
-
Ingår i: VI International Conference on Computational Methods for Coupled Problems in Science and Engineering (COUPLED PROBLEMS 2015), Proceedings of the. - 9788494392832 ; , s. 452-463
-
Konferensbidrag (refereegranskat)abstract
- We present an estimate for the convergence rate of the Dirichlet-Neumann iteration for the discretized unsteady transmission problem. Specifically, we consider the coupling of two heat equations on two identical squared domains. The Laplacian is discretized by second order central finite differences and the implicit Euler method is used for the time discretization. For the semidiscrete case, Henshaw and Chad provided in 2009 a method to analyse stability and convergence speed based on applying the continuous Fourier transform to the semi-discretized equations. Numerical results for the fully discrete case show differences, which is why we propose a complementary analysis based on approximating the spectral radius of the iteration matrix. Numerical results are presented to illustrate the analysis.
|
|
| 35. |
- Monge, Azahar, et al.
(författare)
-
On the convergence rate of the Dirichlet-Neumann iteration for unsteady thermal fluid structure interaction
- 2018
-
Ingår i: Computational Mechanics. - : Springer Science and Business Media LLC. - 0178-7675 .- 1432-0924. ; 62:3, s. 525-541
-
Tidskriftsartikel (refereegranskat)abstract
- We consider the Dirichlet-Neumann iteration for partitioned simulation of thermal fluid-structure interaction, also called conjugate heat transfer. We analyze its convergence rate for two coupled fully discretized 1D linear heat equations with jumps in the material coefficients across these. These are discretized using implicit Euler in time, a finite element method on one domain, a finite volume method on the other one and variable aspect ratio. We provide an exact formula for the spectral radius of the iteration matrix. This shows that for large time steps, the convergence rate is the aspect ratio times the quotient of heat conductivities and that decreasing the time step will improve the convergence rate. Numerical results confirmthe analysis and show that the 1D formula is a good estimator in 2D and even for nonlinear thermal FSI applications.
|
|
| 36. |
- Monge, Azahar, et al.
(författare)
-
Towards a Time Adaptive Neumann-Neumann Waveform Relaxation Method for Thermal Fluid-Structure Interaction
- 2020
-
Ingår i: Domain Decomposition Methods in Science and Engineering XXV, DD 2018. - Cham : Springer International Publishing. - 2197-7100 .- 1439-7358. - 9783030567491 ; 138, s. 466-473
-
Konferensbidrag (refereegranskat)abstract
- Our prime motivation is thermal fluid-structure interaction (FSI) where two domains with jumps in the material coefficients are connected through an interface. There exist two main strategies to simulate FSI models: the monolithic approach where a new code is tailored for the coupled equations and the partitioned approach that allows to reuse existing software for each sub-problem. Here we want to develop multirate methods that contribute to the time parallelization of the sub-problems for the partitioned simulation of FSI problems.
|
|
| 37. |
- Osswald, Kai, et al.
(författare)
-
L2Roe: A low-dissipation version of Roe's approximate Riemann solver for low Mach numbers
- 2015
-
Ingår i: International Journal for Numerical Methods in Fluids. - : Wiley. - 1097-0363 .- 0271-2091.
-
Tidskriftsartikel (refereegranskat)abstract
- A modification of the Roe scheme called L2Roe for low dissipation low Mach Roe is presented. It reduces the dissipation of kinetic energy at the highest resolved wave numbers in a low Mach number test case of decaying isotropic turbulence. This is achieved by scaling the jumps in all discrete velocity components within the numerical flux function. An asymptotic analysis is used to show the correct pressure scaling at low Mach numbers and to identify the reduced numerical dissipation in that regime. Furthermore, the analysis allows a comparison with two other schemes that employ different scaling of discrete velocity jumps, namely, LMRoe and a method of Thornber et al. To this end, we present for the first time an asymptotic analysis of the last method. Numerical tests on cases ranging from low Mach number (M∞=0.001) to hypersonic (M∞=5) viscous flows are used to illustrate the differences between the methods and to show the correct behavior of L2Roe. No conflict is observed between the reduced numerical dissipation and the accuracy or stability of the scheme in any of the investigated test cases.
|
|
| 38. |
- Rüth, Benjamin, et al.
(författare)
-
Quasi-Newton waveform iteration for partitioned surface-coupled multiphysics applications
- 2021
-
Ingår i: International Journal for Numerical Methods in Engineering. - : Wiley. - 0029-5981 .- 1097-0207. ; 122:19, s. 5236-5257
-
Tidskriftsartikel (refereegranskat)abstract
- We present novel coupling schemes for partitioned multiphysics simulation that combine four important aspects for strongly coupled problems: implicit coupling per time step, fast and robust acceleration of the corresponding iterative coupling, support for multirate time stepping, and higher-order convergence in time. To achieve this, we combine waveform relaxation—a known method to achieve higher-order in applications with split time stepping based on continuous representations of coupling variables in time— with interface quasi-Newton coupling, which has been developed throughout the last decade and is generally accepted as a very robust iterative coupling method even for gluing together black-box simulation codes. We show convergence results (in terms of convergence of the iterative solver and in terms of approximation order in time) for two academic testcases—a heat transfer scenario and a fluid-structure interaction simulation. We show that we achieve the expected approximation order and that our iterative method is competitive in terms of iteration counts with those designed for simpler first-order-in-time coupling.
|
|
| 39. |
- Schlutow, Mark, et al.
(författare)
-
Spectral stability of nonlinear gravity waves in the atmosphere
- 2019
-
Ingår i: Mathematics of Climate and Weather Forecasting. - : Portico. - 2353-6438. ; 5:1, s. 12-33
-
Tidskriftsartikel (refereegranskat)abstract
- We apply spectral stability theory to investigate nonlinear gravity waves in the atmosphere. These waves are determined by modulation equations that result from Wentzel-Kramers-Brillouin theory. First, we establish that plane waves, which represent exact solutions to the inviscid Boussinesq equations, are spectrally stable with respect to their nonlinear modulation equations under the same conditions as what is known as modulational stability from weakly nonlinear theory. In contrast to Boussinesq, the pseudo-incompressible regime does fully account for the altitudinal varying background density. Second, we show for the first time that upward-traveling non-plane wave fronts solving the inviscid nonlinear modulation equations, that compare to pseudo-incompressible theory, are unconditionally unstable. Both inviscid regimes turn out to be ill-posed as the spectra allow for arbitrarily large instability growth rates. Third, a regularization is found by including dissipative effects. The corresponding nonlinear traveling wave solutions have localized amplitude. As a consequence of the nonlinearity, envelope and linear group velocity, as given by the derivative of the frequency with respect to wavenumber, do not coincide anymore. These waves blow up unconditionally by embedded eigenvalue instabilities but the instability growth rate is bounded from above and can be computed analytically. Additionally, all three types of nonlinear modulation equations are solved numerically to further investigate and illustrate the nature of the analytic stability results.
|
|
| 40. |
|
|
| 41. |
- Schäfer, Jonas, et al.
(författare)
-
A Memory-Efficient Finite Volume Method for Advection-Diffusion-Reaction Systems with Non-Smooth Sources
- 2015
-
Ingår i: Numerical Methods for Partial Differential Equations. - : Wiley. - 1098-2426 .- 0749-159X. ; 31:1, s. 143-167
-
Tidskriftsartikel (refereegranskat)abstract
- We present a parallel matrix-free implicit nite volume scheme for the solution of unsteady three-dimensional advection-diusion-reaction equations with smooth and Dirac-Delta source terms. The scheme is formally second order in space and a Newton-Krylov method is employed for the appearing nonlinear systems in the implicit time integration. The matrix-vector product required is hardcoded without any approximations, obtaining a matrix-free method that needs little storage and is well suited for parallel implementation. We describe the matrix-free implementation of the method in detail and give numerical evidence of its second order convergence in the presence of smooth source terms. For non-smooth source terms the convergence order drops to one half. Furthermore, we demonstrate the method's applicability for the long time simulation of calcium ow in heart cells and show its parallel scaling.
|
|
| 42. |
|
|
| 43. |
- Straub, Veronika, et al.
(författare)
-
A new domain‐based implicit‐explicit time stepping scheme based on the class of exponential integrators called sEPIRK
- 2019. - 1
-
Ingår i: Special Issue: 90th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM). - : Wiley. - 1617-7061. ; 19:1
-
Konferensbidrag (refereegranskat)abstract
- Reliable simulations of flows in real applications involve the task of discretizing both space and time in an accurate and efficient way. To cope with the large semidiscrete systems resulting from a space discretization on appropriate grids, which often include comparatively few very small cells near solid walls for boundary layer resolution, an implicit-explicit time stepping scheme can be the most efficient variant. We present such a type of scheme utilizing recent exponential integrators called sEPIRK and show its computational advantages opposed to IMEX-Runge-Kutta schemes.
|
|
| 44. |
- Straub, Veronika, et al.
(författare)
-
Adopting (s)EPIRK schemes in a domain-based IMEX setting
- 2017
-
Ingår i: Proceedings of ICNAAM 2016, AIP Conference Proceedings 1863. - : Author(s). - 9780735415386 ; 410008
-
Konferensbidrag (refereegranskat)abstract
- The simulation of viscous, compressible flows around complex geometries or similar applications often inherit the task of solving large, stiff systems of ODEs. Domain-based implicit-explicit (IMEX) type schemes offer the possibility to apply two different schemes to different parts of the computational domain. The goal hereby is to decrease the computational cost by increasing the admissible step sizes with no loss of stability and by reducing the system sizes of the linear solver within the implicit integrator. But which combination of methods reaches the largest gain in efficiency? Coupling of Runge-Kutta methods or different multistep methods has been investigated so far by other authors. Here, we inspect the adoption of the recently introduced exponential integrators called EPIRK and sEPIRK in the IMEX setting, since they are perfectly suited for large, stiff systems of ODEs.
|
|
| 45. |
- Veronika, Straub, et al.
(författare)
-
Efficient Time Integration of IMEX Type using Exponential Integrators for Compressible, Viscous Flow Simulation
- 2016
-
Ingår i: PAMM - Proceedings in Applied Mathematics and Mechanics. - : Wiley. - 1617-7061. ; 16, s. 867-868
-
Tidskriftsartikel (refereegranskat)abstract
- We investigate the adaption of the recently developed exponential integrators called EPIRK in the so-called domain-based implicit-explicit (IMEX) setting of spatially discretized PDE's. The EPIRK schemes were shown to be efficient for sufficiently stiff problems and offer high precision and good stability properties like A- and L-stability in theory. In practice, however, we can show that these stability properties are dependent on the parameters of the interior approximation techniques.Here, we introduce the IMEX-EPIRK method, which consists of coupling an explicit Runge-Kutta scheme with an EPIRK scheme. We briefly analyze its linear stability, show its conservation property and set up a CFL condition. Though the method is convergent of only first order, it demonstrates the advantages of this novel type of schemes for stiff problems very well.
|
|
| 46. |
- Versbach, Lea Miko, et al.
(författare)
-
Finite volume based multigrid preconditioners for discontinuous Galerkin methods
- 2018
-
Ingår i: Special Issue: 89th Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM). - : Wiley. - 1617-7061. ; 18:1
-
Konferensbidrag (övrigt vetenskapligt/konstnärligt)abstract
- Our aim is to construct efficient preconditioners for high order discontinuous Galerkin (DG) methods. We consider the DG spectral element method with Gauss‐Lobatto‐Legendre nodes (DGSEM‐GL) for the 1D linear advection equation. It has been shown in [4] that DGSEM‐GL has the summation‐by‐parts (SBP) property and an equivalent finite volume (FV) discretization is presented in [3]. Thus we present a multigrid (MG) preconditioner based on a simplified FV discretization.
|
|
| 47. |
- Versbach, Lea Miko, et al.
(författare)
-
Theoretical and Practical Aspects of Space-Time DG-SEM Implementations
- 2023
-
Ingår i: The SMAI Journal of computational mathematics. - 2426-8399. ; 9, s. 61-93
-
Tidskriftsartikel (refereegranskat)abstract
- We discuss two approaches for the formulation and implementation of space-time discontinuous Galerkin spectral element methods (DG-SEM). In one, time is treated as an additional coordinate direction and a Galerkin procedure is applied to the entire problem. In the other, the method of lines is used with DG-SEM in space and the fully implicit Runge–Kutta method Lobatto IIIC in time. The two approaches are mathematically equivalent in the sense that they lead to the same discrete solution. However, in practice they differ in several important respects, including the terminology used to describe them, the structure of the resulting software, and the interaction with nonlinear solvers. Challenges and merits of the two approaches are discussed with the goal of providing the practitioner with sufficient consideration to choose which path to follow. Additionally, implementations of the two methods are provided as a starting point for further development. Numerical experiments validate the theoretical accuracy of these codes and demonstrate their utility, even for 4D problems.
|
|
