High-order discontinuous Galerkin discretization for flows with strong moving shocks
Supersonic flows over both simple and complex geometries involve features over a wide spectrum of spatial and temporal scales, whose resolution in a numerical solution is of significant importance for accurate predictions in engineering applications. While CFD has been greatly developed in the la...
| Κύριος συγγραφέας: | |
|---|---|
| Άλλοι συγγραφείς: | |
| Μορφή: | Thesis |
| Γλώσσα: | English |
| Έκδοση: |
2013
|
| Θέματα: | |
| Διαθέσιμο Online: | http://hdl.handle.net/10889/5839 |
| id |
nemertes-10889-5839 |
|---|---|
| record_format |
dspace |
| institution |
UPatras |
| collection |
Nemertes |
| language |
English |
| topic |
High-order Discontinous galerkin Shock capturing Strong moving shocks Ασυνεχή πεπερασμένα στοιχεία Περιοριστές 629.132 305 |
| spellingShingle |
High-order Discontinous galerkin Shock capturing Strong moving shocks Ασυνεχή πεπερασμένα στοιχεία Περιοριστές 629.132 305 Κοντζιάλης, Κωνσταντίνος High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| description |
Supersonic flows over both simple and complex geometries involve features over a
wide spectrum of spatial and temporal scales, whose resolution in a numerical
solution is of significant importance for accurate predictions in engineering
applications. While CFD has been greatly developed in the last 30 years, the
desire and necessity to perform more complex, high fidelity simulations still
remains.
The present thesis has introduced two major innovations regarding the fidelity
of numerical solutions of the compressible \ns equations. The first one is the
development of new a priori mesh quality measures for the Finite
Volume (FV) method on mixed-type (quadrilateral/triangular) element meshes.
Elementary types of mesh distortion were identified expressing grid distortion
in terms of stretching, skewness, shearing and non-alignment of the mesh.
Through a rigorous truncation error analysis, novel grid quality measures were
derived by emphasizing on the direct relation between mesh distortion and the
quality indicators. They were applied over several meshes and their ability was
observed to identify faithfully irregularly-shaped small or large distortions in
any direction. It was concluded that accuracy degradation occurs even for small
mesh distortions and especially at mixed-type element mesh interfaces the formal
order of the FV method is degraded no matter of the mesh geometry and local mesh
size.
Therefore, in the present work, the high-order Discontinuous Galerkin (DG)
discretization of the compressible flow equations was adopted as a means of
achieving and attaining high resolution of flow features on irregular mixed-type
meshes for flows with strong moving shocks. During the course of the
thesis a code was developed and named HoAc (standing for High Order Accuracy),
which can perform via the domain decomposition method parallel $p$-adaptive
computations for flows with strong shocks on mixed-type element meshes over
arbitrary geometries at a predefined arbitrary order of accuracy. In HoAc in
contrast to other DG developments, all the numerical operations are performed in
the computational space, for all element types. This choice constitutes the key
element for the ability to perform $p$-adaptive computations along with modal
hierarchical basis for the solution expansion. The time marching of the DG
discretized Navier-Stokes system is performed with the aid of explicit Runge-Kutta methods or with a matrix-free implicit approach.
The second innovation of the present thesis, which is also based on the choice
of implementing the DG method on the regular computational space, is the
development of a new $p$-adaptive limiting procedure for shock capturing of the
implemented DG discretization. The new limiting approach along with positivity
preserving limiters is suitable for computations of high speed flows with strong
shocks around complex geometries. The unified approach for $p$-adaptive limiting
on mixed-type meshes is achieved by applying the limiters on the transformed
canonical elements, and it is fully automated without the need of
ad hoc specification of parameters as it has been done with
standard limiting approaches and in the artificial dissipation method for shock
capturing.
Verification and validation studies have been performed, which prove the
correctness of the implemented discretization method in cases where the linear
elements are adequate for the tessellation of the computational domain both for
subsonic and supersonic flows. At present HoAc can handle only linear elements
since most grid generators do not provide meshes with curved elements.
Furthermore, p-adaptive computations with the implemented DG method were
performed for a number of standard test cases for shock capturing schemes to
illustrate the outstanding performance of the proposed $p$-adaptive limiting
approach. The obtained results are in excellent agreement with analytical
solutions and with experimental data, proving the excellent efficiency of the
developed shock capturing method for the DG discretization of the equations of
gas dynamics. |
| author2 |
Αικατερινάρης, Ιωάννης |
| author_facet |
Αικατερινάρης, Ιωάννης Κοντζιάλης, Κωνσταντίνος |
| format |
Thesis |
| author |
Κοντζιάλης, Κωνσταντίνος |
| author_sort |
Κοντζιάλης, Κωνσταντίνος |
| title |
High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| title_short |
High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| title_full |
High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| title_fullStr |
High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| title_full_unstemmed |
High-order discontinuous Galerkin discretization for flows with strong moving shocks |
| title_sort |
high-order discontinuous galerkin discretization for flows with strong moving shocks |
| publishDate |
2013 |
| url |
http://hdl.handle.net/10889/5839 |
| work_keys_str_mv |
AT kontzialēskōnstantinos highorderdiscontinuousgalerkindiscretizationforflowswithstrongmovingshocks |
| _version_ |
1771297296062873600 |
| spelling |
nemertes-10889-58392022-09-05T20:35:52Z High-order discontinuous Galerkin discretization for flows with strong moving shocks Κοντζιάλης, Κωνσταντίνος Αικατερινάρης, Ιωάννης Καλλιντέρης, Ιωάννης Γιαννάκογλου, Κυριάκος Κωστόπουλος, Βασίλειος Σαραβάνος, Δημήτριος Καλλιντέρης, Ιωάννης Γιαννάκογλου, Κυριάκος Γαλόπουλος, Στρατής Τσαμόπουλος, Ιωάννης Kontzialis, Konstantinos High-order Discontinous galerkin Shock capturing Strong moving shocks Ασυνεχή πεπερασμένα στοιχεία Περιοριστές 629.132 305 Supersonic flows over both simple and complex geometries involve features over a wide spectrum of spatial and temporal scales, whose resolution in a numerical solution is of significant importance for accurate predictions in engineering applications. While CFD has been greatly developed in the last 30 years, the desire and necessity to perform more complex, high fidelity simulations still remains. The present thesis has introduced two major innovations regarding the fidelity of numerical solutions of the compressible \ns equations. The first one is the development of new a priori mesh quality measures for the Finite Volume (FV) method on mixed-type (quadrilateral/triangular) element meshes. Elementary types of mesh distortion were identified expressing grid distortion in terms of stretching, skewness, shearing and non-alignment of the mesh. Through a rigorous truncation error analysis, novel grid quality measures were derived by emphasizing on the direct relation between mesh distortion and the quality indicators. They were applied over several meshes and their ability was observed to identify faithfully irregularly-shaped small or large distortions in any direction. It was concluded that accuracy degradation occurs even for small mesh distortions and especially at mixed-type element mesh interfaces the formal order of the FV method is degraded no matter of the mesh geometry and local mesh size. Therefore, in the present work, the high-order Discontinuous Galerkin (DG) discretization of the compressible flow equations was adopted as a means of achieving and attaining high resolution of flow features on irregular mixed-type meshes for flows with strong moving shocks. During the course of the thesis a code was developed and named HoAc (standing for High Order Accuracy), which can perform via the domain decomposition method parallel $p$-adaptive computations for flows with strong shocks on mixed-type element meshes over arbitrary geometries at a predefined arbitrary order of accuracy. In HoAc in contrast to other DG developments, all the numerical operations are performed in the computational space, for all element types. This choice constitutes the key element for the ability to perform $p$-adaptive computations along with modal hierarchical basis for the solution expansion. The time marching of the DG discretized Navier-Stokes system is performed with the aid of explicit Runge-Kutta methods or with a matrix-free implicit approach. The second innovation of the present thesis, which is also based on the choice of implementing the DG method on the regular computational space, is the development of a new $p$-adaptive limiting procedure for shock capturing of the implemented DG discretization. The new limiting approach along with positivity preserving limiters is suitable for computations of high speed flows with strong shocks around complex geometries. The unified approach for $p$-adaptive limiting on mixed-type meshes is achieved by applying the limiters on the transformed canonical elements, and it is fully automated without the need of ad hoc specification of parameters as it has been done with standard limiting approaches and in the artificial dissipation method for shock capturing. Verification and validation studies have been performed, which prove the correctness of the implemented discretization method in cases where the linear elements are adequate for the tessellation of the computational domain both for subsonic and supersonic flows. At present HoAc can handle only linear elements since most grid generators do not provide meshes with curved elements. Furthermore, p-adaptive computations with the implemented DG method were performed for a number of standard test cases for shock capturing schemes to illustrate the outstanding performance of the proposed $p$-adaptive limiting approach. The obtained results are in excellent agreement with analytical solutions and with experimental data, proving the excellent efficiency of the developed shock capturing method for the DG discretization of the equations of gas dynamics. - 2013-02-04T11:59:19Z 2013-02-04T11:59:19Z 2012 2013-02-04 Thesis http://hdl.handle.net/10889/5839 en_US Η ΒΚΠ διαθέτει αντίτυπο της διατριβής σε έντυπη μορφή στο βιβλιοστάσιο διδακτορικών διατριβών που βρίσκεται στο ισόγειο του κτιρίου της. 0 application/pdf |
