4,495 research outputs found

    Review of Summation-by-parts schemes for initial-boundary-value problems

    Full text link
    High-order finite difference methods are efficient, easy to program, scales well in multiple dimensions and can be modified locally for various reasons (such as shock treatment for example). The main drawback have been the complicated and sometimes even mysterious stability treatment at boundaries and interfaces required for a stable scheme. The research on summation-by-parts operators and weak boundary conditions during the last 20 years have removed this drawback and now reached a mature state. It is now possible to construct stable and high order accurate multi-block finite difference schemes in a systematic building-block-like manner. In this paper we will review this development, point out the main contributions and speculate about the next lines of research in this area

    A high order compact scheme for hypersonic aerothermodynamics

    Get PDF
    A novel high order compact scheme for solving the compressible Navier-Stokes equations has been developed. The scheme is an extension of a method originally proposed for solving the Euler equations, and combines several techniques for the solution of compressible flowfields, such as upwinding, limiting and flux vector splitting, with the excellent properties of high order compact schemes. Extending the method to the Navier-Stokes equations is achieved via a Kinetic Flux Vector Splitting technique, which represents an unusual and attractive way to include viscous effects. This approach offers a more accurate and less computationally expensive technique than discretizations based on more conventional operator splitting. The Euler solver has been validated against several inviscid test cases, and results for several viscous test cases are also presented. The results confirm that the method is stable, accurate and has excellent shock-capturing capabilities for both viscous and inviscid flows

    Improved shock-capturing of Jameson's scheme for the Euler equations

    Get PDF
    It is known that Jameson's scheme is a pseudo-second-order-accurate scheme for solving discrete conservation laws. The scheme contains a non-linear artificial dissipative flux which is designed to capture shocks. In this paper, it is shown that the, shock-capturing of Jameson's scheme for the Euler equations can be improved by replacing the Lax-Friedrichs' type of dissipative flux with Roe's dissipative flux. This replacement is at a moderate expense of the calculation time

    Implicit High-Order Flux Reconstruction Solver for High-Speed Compressible Flows

    Full text link
    The present paper addresses the development and implementation of the first high-order Flux Reconstruction (FR) solver for high-speed flows within the open-source COOLFluiD (Computational Object-Oriented Libraries for Fluid Dynamics) platform. The resulting solver is fully implicit and able to simulate compressible flow problems governed by either the Euler or the Navier-Stokes equations in two and three dimensions. Furthermore, it can run in parallel on multiple CPU-cores and is designed to handle unstructured grids consisting of both straight and curved edged quadrilateral or hexahedral elements. While most of the implementation relies on state-of-the-art FR algorithms, an improved and more case-independent shock capturing scheme has been developed in order to tackle the first viscous hypersonic simulations using the FR method. Extensive verification of the FR solver has been performed through the use of reproducible benchmark test cases with flow speeds ranging from subsonic to hypersonic, up to Mach 17.6. The obtained results have been favorably compared to those available in literature. Furthermore, so-called super-accuracy is retrieved for certain cases when solving the Euler equations. The strengths of the FR solver in terms of computational accuracy per degree of freedom are also illustrated. Finally, the influence of the characterizing parameters of the FR method as well as the the influence of the novel shock capturing scheme on the accuracy of the developed solver is discussed

    A Study of Flow Separation in Transonic Flow Using Inviscid and Viscous Computational Fluid Dynamics (CFD) Schemes

    Get PDF
    A comparison of flow separation in transonic flows is made using various computational schemes which solve the Euler and the Navier-Stokes equations of fluid mechanics. The flows examined are computed using several simple two-dimensional configurations including a backward facing step and a bump in a channel. Comparison of the results obtained using shock fitting and flux vector splitting methods are presented and the results obtained using the Euler codes are compared to results on the same configurations using a code which solves the Navier-Stokes equations

    Entropy stable boundary conditions for the Euler equations

    Get PDF
    We consider the initial-boundary value Euler equations with the aim to derive boundary conditions that yield an entropy bound for the physical (Navier-Stokes) entropy. We begin by reviewing the entropy bound obtained for standard no-penetration wall boundary conditions and propose a numerical implementation. The main results are the derivation of full-state boundary conditions (far-field, inlet, outlet) and the accompanying entropy stable implementations. We also show that boundary conditions obtained from linear theory are unable to bound the entropy and that non-linear bounds require additional boundary conditions. We corroborate our theoretical findings with numerical experiments.publishedVersio
    • ā€¦
    corecore