A Fourth-Order Gas-Kinetic CPR Method for the Navier-Stokes Equations on Unstructured Meshes
|
|
- August Edwards
- 5 years ago
- Views:
Transcription
1 Tenth International Conference on Computational Fluid Dnamics (ICCFD), Barcelona,Spain, Jul 9-, 8 ICCFD-8- A Fourth-Order Gas-Kinetic CPR Method for the Navier-Stokes Equations on Unstructured Meshes Chao Zhang and Qibing Li Corresponding author: lqb@tsinghua.edu.cn AML, Department of Engineering Mechanics, Tsinghua Universit, Beijing 84, China. Abstract: A fourth-order accurate gas-kinetic CPR method is developed for the Navier-Stokes equations on triangular meshes. Different from the previous single-stage third-order gas-kinetic CPR method, the current method adopts the second-order gas-kinetic flu solver, which is simpler and less epensive, to reduce the computational cost of flu. Meanwhile, the temporal accurac is improved b using the two-stage temporal discretization, which is more efficient than the four-stage Runge-Kutta method usuall adopted in eisting fourth-order CPR. Numerical tests are presented to demonstrate the accurac and efficienc of the current method. Kewords: Gas-Kinetic Scheme, Correction Procedure via Reconstruction, Two-stage Temporal Discretization. Introduction During the past decades, high-order CFD methods have attracted man researchers due to their advantages on accurac and efficienc. There are man popular high-order methods, such as the Discontinuous Galerkin (DG) method, the P N P M and so on. In recent ears, the correction procedure via reconstruction (CPR) method has been developed rapidl, which provides a unified framework for man high-order methods []. The CPR method has been shown to be ver competitive in terms of accurac and efficienc, compared to the original DG method. Most of high-order methods mainl concentrate on the high-order reconstruction. As for the flu evolution, the Riemann solution of the Euler equations is the foundation of compressible flu solvers. However, developing a genuinel multidimensional Riemann solution is ver difficult. For solving viscous flows, the viscous terms in the Navier-Stokes equations require special treatment. Alternativel, the gas-kinetic scheme (GKS) offers a different wa to recover the N-S solutions, which has shown good performance in a wide range of flow problems. Based on the local integral solution of the BGK model, the time evolving kinetic flu function can be constructed, in which the inviscid and viscous terms are coupled and obtained simultaneousl. Through the high-order epansion of the gas distribution function, a high-order gas evolution model has been constructed successfull, which provides more abundant information to describe the flu evolution process []. Accordingl, a series of high-order gas-kinetic schemes have been developed. Recentl, a third-order gas-kinetic CPR method, denoted as CPR-GKS(P), has been developed for the Navier-Stokes equations on unstructured meshes []. It combines the high-efficient CPR framework with the third-order single-stage time stepping gas-kinetic flu solver, which shows high accurac and efficienc in man tpical flow problems. However, the third-order gas-kinetic flu solver is more complicated and time-consuming than its second-order counterpart. In this stud the fourth-order CPR framework is adopted to enhance the spatial accurac. The secondorder time-dependent gas-kinetic flu solver [4] is adopted for less computational cost, combined with the two-stage time-stepping method [5, 6] to improve the temporal accurac. In comparison with the four-stage Runge-Kutta time-stepping method, the temporal discretization is more efficient due to less intermediate stages. The developed scheme, denoted as CPR-GKS(P) is then validated through tpical numerical tests.
2 Fourth-Order Gas-Kinetic CPR Method. CPR framework B dividing the computational domain into N non-overlapping triangular cells Ω i, a solution polnomial Q i of degree k is constructed in each cell b using the Lagrange interpolation based on a set of solution points (SPs). For fourth-order CPR, i.e., k =, solution points are needed in each cell. Additionall, to compute the common flu, k + flu points (FPs) are set at each cell interface. For efficienc, the solution points are chosen the same as the flu points, as shown in Fig.. Figure. Solution points (circles) and flu points (squares) for k=. Then the CPR framework can be epressed as Q i,j t + Π j ( F ) (Q i ) + Ω i s Ω i l= 4 α j,s,l [F n ] s,l Γ s =, () where Q i,j is the conservative variable at each solution point, j =. The second term is the projection of the flu divergence onto P k. The last one is the correction term where [F n ] s,l is the normal flu difference at each flu point, i.e., [F n com(q i, Q i+, n) F n (Q i )] s,l. α j,s,l is the lifting coefficient, Ω i is the area of Ω i, Γ s is the length of triangular edge. More details of the flu divergence term and correction term can be found in [].. Gas-kinetic scheme In the present stud, the flu in the CPR framework is computed b the gas-kinetic flu solver. For completeness, the gas-kinetic scheme is briefl reviewed. For convenience, the summation convention is adopted in the following, where = (, ) = (, ), u = (u, u ) = (u, v), macroscopic velocit vector U = (U, U ) = (U, V ). The scheme is based on the BGK equation, f t + u f i = g f, () i τ where f is the gas distribution function, u i is the particle velocit, and τ = µ/p is the particle collision time related to the viscosit and p is the pressure. The equilibrium state g is the Mawellian distribution g = ρ ( ) K+ λ e λ( u U +ξ ), () π where ρ is the densit. λ is equal to /(RT ), where R is the gas constant and T is the temperature. Since mass, momentum, and energ are conserved during particle collisions, the collision term in Eq.() satisfies the compatibilit condition g f ψdξ =, (4) τ
3 where ψ is the the vector of moments ψ = (, u, ( u + ξ )) T, (5) and dξ = du du dξ dξ dξ K is the element of the phase space and ξ = ξ + ξ ξk. The total number of degrees of freedom K is equal to (5 γ)/(γ ) + for two-dimensional flow, in which γ is the specific heat ratio. The relations between macroscopic variables and the distribution function are Q = fψdξ, F = ufψdξ. (6) Taking moments of BGK equation, the N-S equations can be recovered based on the first-order Chapman- Enskog epansion f NS = g τ( g t + u g i ). (7) i In order to compute the flu at each solution point and flu point, the gas distribution function is constructed locall at these points. Since the flow distribution inside each cell is continuous, Eq.(7) can be directl adopted for computing the flu at each solution point (SP). The corresponding gas distribution function can be epressed as f SP (t, u, ξ) =g ( τa i u i + (t τ)a). (8) B introducing the notation = g ( )ψdξ, (9) the coefficients a i and A can be determined b the spatial derivatives of conservative variables and the compatibilit condition [] a i = Q i a i, a i u i + A = A. () To compute the common flu at each flu point (FP) where the discontinuit eists, the gas distribution function is constructed based on the local analtical solution of BGK equation, f(, t, u, ξ) = τ t g(, t, u, ξ)e (t t )/τ dt + e t/τ f ( ut,, u, ξ), () where f is the piece-wise discontinuous initial distribution function at the beginning of each time step (t = ), = u(t t ) is the particle trajector. Note that the construction is performed in a local coordinate sstem = (, ), in which the -ais is perpendicular to the cell interface. Based on Eq.(), the second-order time-dependent gas distribution function can be epressed as ( f F P (t, ũ, ) ) ( ξ) =g e t/τ + g ((t ) + τ)e t/τ τ a i ũ i + g t τ + τe t/τ A + e t/τ g R ( (τ + t)a R i ũ i τa R) H(u) () + e t/τ g L ( (τ + t)a L i ũ i τa L) ( H(u)), where H(u) is the Heaviside function. The equilibrium states g L, g R correspond to the flow states at the left and right sides of the cell interface. The coefficients a i, A, a L i, AL and a R i, AR come from derivatives of g, g L and g R respectivel, which can be determined b Eq.() as well. Besides, the equilibrium state g is
4 determined b the compatibilit condition g ψdξ = Q = u > g L ψdξ + u < g R ψdξ, () where Q are the conservative variables corresponding the equilibrium state g. In order to full determine g, a continuous flow distribution across the cell interface is constructed b using the least-square method, with the stencil shown in Fig.. Figure. The stencil for constructing the equilibrium state in the local coordinate.. Two-stage temporal discretization Based on the above gas-kinetic flu solver, the solution at each solution point inside a cell can be updated b adopting the following two-stage temporal discretization, ( t Q i,j = Q n i,j Π j ( ) Q n+ i,j = Q n i,j Π j ( F )) n i, t)dt Ω i Ω i s Ω i ) α j,s,l (ˆF l s Ω i s,l ( t α j,s,l l Γ s, ) ˆF (Q n i, t)dt Γ s, s,l where t = + t/ is the intermediate stage, ˆF (Q n i, t) indicates the normal flu difference in Eq.(), and F can be epressed as 8 t F = n i, t)dt tn+ n i, t)dt + 4 t + t i, t)dt 8 tn+ t i, t)dt, (5) t and ˆF has the same form as Eq.(5) b replacing F with ˆF. It is not eas to construct a fourth-order gas-kinetic flu solver with one-step time discretization as the corresponding gas distribution function is complicated and the computational cost is far more epensive. Fortunatel, with the above two-step time discretization, the second-order gas-kinetic flu solver can be used, which onl takes about one-sith the computational cost of its third-order counterpart. Therefore, although the current method computes the flu pointwisel and has two stages, the total cost of flu is at the same level as the single-stage third-order gas-kinetic CPR, however, both spatial and temporal accurac are improved. Furthermore, the current temporal discretization is more efficient than the fourth-order Runge-Kutta method due to less intermediate stages. As a result, although the gas-kinetic flu solver is more epensive than traditional flu solvers, the current scheme can achieve nearl the same efficienc as the original fourth-order CPR. What s more, when simulating high-speed compressible viscous flows, the slope limiting procedure is necessar to capture flow discontinuities, which takes a large part of computational time. Therefore, it can be epected that the current scheme can be more efficient than the original CPR in compressible viscous flows. (4) 4
5 Numerical Results. Compressible Couette flow In order to verif the accurac of the current scheme in viscous flows, the compressible Couette flow is simulated. The traditional CPR is also tested b using the HLLC and BR scheme with the fourth-stage Runge-Kutta method, denoted as CPR-BR(P). The lower wall with = is stationar and adiabatic. The upper wall at = is moving with a constant speed U and temperature T =. There is an analtical solution for this stead flow, i.e., ( + Pr γ ) ( Ma H = U + Pr γ Ma U ( ) ) U, U U U T V =, = + Pr γ ( ) ) (6) U Ma (, p = /γ, T U where the Prandtl number is Pr =, the Mach number is Ma =.5. The viscosit is determined b the linear law µ = µ T/T. The Renolds number is Re = 5. The CFL number to compute time steps is set as. in all test cases Figure. Computational mesh for compressible Couette flow. - CPR-GKS (P) CPR-BR (P) -5 4th order rd order - CPR-GKS (P) CPR-BR (P) -5 L error -7-9 L error Mesh size CPU time Figure 4. Temperature error vs. mesh size (left) and CPU time (right) for the compressible Couette flow. The coarsest mesh used in this case is shown in Fig.. The temperature error vs. mesh size and CPU time are presented in figure 4, in which the CPU time is the computational cost to achieve a stead state when the temperature residual less than 5. It shows that CPR-GKS(P) can achieve the fourth-order accurac. To achieve the same level of absolute error, CPR-GKS(P) is more efficient than CPR-GKS(P) and has nearl the same efficienc as CPR-BR(P), demonstrating the high accurac and efficienc of the current scheme. Nevertheless, further tests are necessar, especiall for supersonic flows, to evaluate the 5
6 robustness of the current scheme. It can be anticipated that, with less intermediate stages, CPR-GKS(P) can be more efficient than the CPR-BR(P) in supersonic flows when the slope limiting is included which has a considerable computational cost.. Lid-driven cavit flow The second case is the lid-driven cavit flow problem. The top boundar of the unit square [, ] [, ] is a moving plate at a speed U =. The Mach number is set as Ma =.5. The Renolds number is Re =. The non-slip and isothermal boundar conditions are imposed on all boundaries with temperature T w =. The initial flow is stationar with densit ρ = and tempture T = T w. As presented in Fig.5, the mesh contains elements with the minimum mesh size h = Figure 5. Computational mesh (left) and streamlines (right) in lid-driven cavit flow. The streamlines are also shown in Fig.5, in which the flow structure, including the primar and secondar vortices, are well captured. The results of U-velocities along the vertical centerline and V -velocities along the horizontal centerline are shown in Fig.6. It can be observed thaumerical results obtained b both CPR-GKS (P) and CPR-BR (P) match ver well with the benchmark data. CPR-GKS (P) is more accurate than with such a coarse mesh, which also demonstrates the high accurac and efficienc of the current scheme...8 CPR-GKS (P) CPR-BR (P) Ghia.6.4. CPR-GKS (P) CPR-BR (P) Ghia U.4 V Figure 6. U-velocities along the vertical centerline and V -velocities along the horizontal centerline with Re= in the lid-driven cavit flow. 6
7 4 Conclusion In this stud, based on the second-order gas-kinetic flu solver and the two-stage temporal discretization, a fourth-order gas-kinetic CPR method is developed for the Navier-Stokes equations on unstructured meshes. Compared with the previous third-order gas-kinetic CPR method, the second-order gas-kinetic flu solver is adopted to reduce the computational cost of flu evaluations, while the temporal accurac can be improved b using the two-stage temporal discretization. Numerical tests show that the current scheme is more efficient than its third-order counterpart and achieves nearl the same efficienc as the original CPR method. Further investigations is required to solve high-speed flows where a more complicated reconstruction with slope limiters is needed to capture flow discontinuities, which takes much more CPU time. Accordingl, it can be epected that, with less intermediate stages, the current scheme can be more efficient than the original CPR method with the four-stage Runge-Kutta method. 5 Acknowledgments This work is supported b the NSFC (6758), National Ke Basic Research and Development Program (4CB744) and Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase) (U55). References [] Z.J. Wang and H.Y. Gao. A unifing lifting collocation penalt formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mied grids. J. Comput. Phs., 8:86 886, 9. [] Q.B. Li, K. Xu, and S. Fu. A high-order gas-kinetic Navier-Stokes flow solver. J. Comput. Phs., 9:675 67,. [] C. Zhang, Q.B. Li, S. Fu, and Z.J. Wang. A third-order gas-kinetic CPR method for the Euler and NavierĺCStokes equations on triangular meshes. J. Comput. Phs., 6:9 5, 8. [4] K. Xu. A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method. J. Comput. Phs., 7:89 5,. [5] J. Li and Z. Du. A two-stage fourth order time-accurate discretization for La-Wendroff tpe flow solvers, I: hperbolic conservation laws. SIAM J. Sci. Comput., 8:A46 A69., 6. [6] L. Pan, K. Xu, Q.B. Li, and J. Li. An efficient and accurate two-stage fourth-order gas-kinetic scheme for the Euler and Navier-Stokes equations. J. Comput. Phs., 6:97, 6. 7
NUMERICAL PERFORMANCE OF COMPACT FOURTH ORDER FORMULATION OF THE NAVIER-STOKES EQUATIONS
Published in : Communications in Numerical Methods in Engineering (008 Commun.Numer.Meth.Engng. 008; Vol : pp 003-019 NUMERICAL PERFORMANCE OF COMPACT FOURTH ORDER FORMULATION OF THE NAVIER-STOKES EQUATIONS
More informationThe gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last
Parallel Implementation of Gas-Kinetic BGK Scheme on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer Middle East Technical
More informationDiscontinuous Galerkin Spectral Element Approximations for CFD
Discontinuous Galerkin Spectral Element Approimations for CFD D.A. Kopriva Florida State Universit Tallahassee, FL 3236 G.B. Jacobs San Diego State Universit San Diego, CA 92182 September 3, 211 1 Code
More informationThe use of gas-kinetic schemes for the simulation of compressible flows become widespread in the two last
A Gas-Kinetic BGK Scheme for Parallel Solution of 3-D Viscous Flows on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationApplication of a Coordinate Transformation and Discretization Method for Computational Fluid Dynamics
Application of a Coordinate Transformation and Discretization Method for Computational Fluid Dnamics 1 Douglas F. Hunsaker Utah State Universit An overview of the computational methods implemented in a
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More informationGrid and Mesh Generation. Introduction to its Concepts and Methods
Grid and Mesh Generation Introduction to its Concepts and Methods Elements in a CFD software sstem Introduction What is a grid? The arrangement of the discrete points throughout the flow field is simpl
More informationNATIONAL UNIVERSITY OF SINGAPORE. (Semester I: 1999/2000) EE4304/ME ROBOTICS. October/November Time Allowed: 2 Hours
NATIONAL UNIVERSITY OF SINGAPORE EXAMINATION FOR THE DEGREE OF B.ENG. (Semester I: 1999/000) EE4304/ME445 - ROBOTICS October/November 1999 - Time Allowed: Hours INSTRUCTIONS TO CANDIDATES: 1. This paper
More informationHierarchical Reconstruction for Spectral Volume Method on Unstructured Grids
Hierarchical Reconstruction for Spectral Volume Method on Unstructured Grids Zhiliang Xu, Yingjie Liu and Chi-Wang Shu April 4, 2009 Abstract The hierarchical reconstruction (HR) [, 24] is applied to a
More informationMarie Curie Actions. IEF Final report
Project No : 234782 Project Acronm : 3DZZI Project Full Name : Three-dimensional structure of stratified turbulence Marie Curie Actions IEF Final report Period covered : from 5/5/29 to 3/9/2 Start date
More informationParameteric study of the performance of two-dimensional Scramjet Intake
8 th Annual CFD Smposium, August -, 26, Bangalore. Parameteric stud of the performance of two-dimensional Scramjet Intake V. Jagadish Babu, Pratikkumar Raje, Rachit Singh, Subhajit Ro and Krishnendu Sinha
More informationPartial Semi-Coarsening Multigrid Method Based on the HOC Scheme on Nonuniform Grids for the Convection-diffusion Problems
International Journal of Computer Mathematics ISSN: 0020-760 (Print) 029-0265 (Online) Journal homepage: http://www.tandfonline.com/loi/gcom20 Partial Semi-Coarsening Multigrid Method Based on the HOC
More informationThe hybrid Grid Implemented DSMC Method used in 2D Triangular Micro Cavity Flows. Nevsan ŞENGİL*
The hbrid Grid Implemented DSMC Method used in 2D Triangular Micro Cavit Flows Nevsan ŞENGİL* * Corresponding author: Tel.: +90 (0)312 342 8458; Fax: +90 (0) 312 342 8460; Email: nsengil@thk.edu.tr Aeronautical
More informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
More informationHigh-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA )
High-Order Numerical Algorithms for Steady and Unsteady Simulation of Viscous Compressible Flow with Shocks (Grant FA9550-07-0195) Sachin Premasuthan, Kui Ou, Patrice Castonguay, Lala Li, Yves Allaneau,
More informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 Outline
More informationPartition Design and Optimization for High-Order Spectral Volume Schemes on Tetrahedral Grids
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 200, Orlando, Florida AIAA 200-720 Partition Design and Optimization for High-Order Spectral Volume
More informationA SCATTERED DATA INTERPOLANT FOR THE SOLUTION OF THREE-DIMENSIONAL PDES
European Conference on Computational Fluid Dnamics ECCOMAS CFD 26 P. Wesseling, E. Oñate and J. Périau (Eds) c TU Delft, The Netherlands, 26 A SCATTERED DATA INTERPOLANT FOR THE SOLUTION OF THREE-DIMENSIONAL
More informationALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving Objects
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona,Spain, July 9-13, 2018 ICCFD10-047 ALE Seamless Immersed Boundary Method with Overset Grid System for Multiple Moving
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationParallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes
Parallel Adaptive Tsunami Modelling with Triangular Discontinuous Galerkin Schemes Stefan Vater 1 Kaveh Rahnema 2 Jörn Behrens 1 Michael Bader 2 1 Universität Hamburg 2014 PDES Workshop 2 TU München Partial
More informationCFD SIMULATION OF FLOW CONTROL WITH ENERGY PROMOTERS IN S-SHAPED DIFFUSER
006-017 Asian Research Publishing Network (ARPN). All rights reserved. CFD SIMULATION OF FLOW CONTROL WITH ENERGY PROMOTERS IN S-SHAPED DIFFUSER Raed A. Jessam, Hussain H. Al-Kaiem and Mohammed Shakir
More informationECNS: Energy-Conserving Navier-Stokes Solver
ECNS: Energ-Conserving Navier-Stokes Solver Verification of stead laminar flows B. Sanderse ECN-E -42 June 2 Abstract This report describes the verification of a second order energ-conserving finite volume
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationComputational Complexity Study on Krylov Integration Factor WENO Method for High Spatial Dimension Convection Diffusion Problems
DOI.7/s9-7-398-7 Computational Complexit Stud on Krlov Integration Factor WENO Method for High Spatial Dimension Convection Diffusion Problems Dong Lu Yong-Tao Zhang Received: 8 Jul 6 / Revised: 3 Februar
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationNon-Linear Filtering and Limiting in High Order Methods for Ideal and Non-Ideal MHD
Universit of Nebraska - Lincoln DigitalCommons@Universit of Nebraska - Lincoln NASA Publications National Aeronautics and Space Administration - Non-Linear Filtering and Limiting in High Order Methods
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationNUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS
NUMERICAL ANALYSIS OF WIND EFFECT ON HIGH-DENSITY BUILDING AERAS Bin ZHAO, Ying LI, Xianting LI and Qisen YAN Department of Thermal Engineering, Tsinghua University Beijing, 100084, P.R. China ABSTRACT
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationMESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP
Vol. 12, Issue 1/2016, 63-68 DOI: 10.1515/cee-2016-0009 MESHLESS SOLUTION OF INCOMPRESSIBLE FLOW OVER BACKWARD-FACING STEP Juraj MUŽÍK 1,* 1 Department of Geotechnics, Faculty of Civil Engineering, University
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
More informationSensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method
2th AIAA Computational Fluid Dynamics Conference 27-3 June 2, Honolulu, Hawaii AIAA 2-348 Sensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method Li Wang
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationA high order moving boundary treatment for compressible inviscid flows 1. Abstract
A high order moving boundary treatment for compressible inviscid flows Sirui Tan and Chi-Wang Shu Abstract We develop a high order numerical boundary condition for compressible inviscid flows involving
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationDeveloping a Tracking Algorithm for Underwater ROV Using Fuzzy Logic Controller
5 th Iranian Conference on Fuzz Sstems Sept. 7-9, 2004, Tehran Developing a Tracking Algorithm for Underwater Using Fuzz Logic Controller M.H. Saghafi 1, H. Kashani 2, N. Mozaani 3, G. R. Vossoughi 4 mh_saghafi@ahoo.com
More informationIntegral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach
Integral Equation Methods for Vortex Dominated Flows, a High-order Conservative Eulerian Approach J. Bevan, UIUC ICERM/HKUST Fast Integral Equation Methods January 5, 2016 Vorticity and Circulation Γ =
More informationIJISET - International Journal of Innovative Science, Engineering & Technology, Vol. 2 Issue 11, November
P P P P IJISET - International Journal of Innovative Science, Engineering & Technology, Vol. 2 Issue, November 205. Numerical Simulation for Steady Incompressible Laminar Fluid Flow and Heat Transfer inside
More informationQuad-Tree Based Geometric-Adapted Cartesian Grid Generation
Quad-Tree Based Geometric-Adapted Cartesian Grid Generation EMRE KARA1, AHMET ĐHSAN KUTLAR1, MEHMET HALUK AKSEL 1Department of Mechanical Engineering Universit of Gaziantep 7310 Gaziantep TURKEY Mechanical
More informationSpatially Non-uniform Time-Step Adaptation in Unsteady Flow Problems
Spatiall Non-uniform Time-Step Adaptation in Unstead Flow Problems Karthik Mani and Dimitri J. Mavriplis Department of Mechanical Engineering, Universit of Woming, Laramie, Woming 827-3295 This paper presents
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationVerification of Moving Mesh Discretizations
Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions
More informationFOURTH ORDER COMPACT FORMULATION OF STEADY NAVIER-STOKES EQUATIONS ON NON-UNIFORM GRIDS
International Journal of Mechanical Engineering and Technology (IJMET Volume 9 Issue 10 October 2018 pp. 179 189 Article ID: IJMET_09_10_11 Available online at http://www.iaeme.com/ijmet/issues.asp?jtypeijmet&vtype9&itype10
More informationFormation of sprays from conical liquid sheets
Center for Turbulence Research Annual Research Briefs 999 249 Formation of spras from conical liquid sheets B Bill Peck AND N. N. Mansour. Motivation and objectives Our objective is to predict droplet
More informationDrop Impact Simulation with a Velocity-Dependent Contact Angle
ILASS Americas 19th Annual Conference on Liquid Atomization and Spray Systems, Toronto, Canada, May 2006 Drop Impact Simulation with a Velocity-Dependent Contact Angle S. Afkhami and M. Bussmann Department
More informationAvailable online at ScienceDirect. Procedia Engineering 99 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 99 (2015 ) 575 580 APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology, APISAT2014 A 3D Anisotropic
More informationNumerical Methods for Hyperbolic and Kinetic Equations
Numerical Methods for Hyperbolic and Kinetic Equations Organizer: G. Puppo Phenomena characterized by conservation (or balance laws) of physical quantities are modelled by hyperbolic and kinetic equations.
More informationAn Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid
An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid J. V. Lassaline Ryerson University 35 Victoria St, Toronto, ON, M5B 2K3, Canada D. W. Zingg University
More informationHigh-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson
More informationA New Concept on Automatic Parking of an Electric Vehicle
A New Concept on Automatic Parking of an Electric Vehicle C. CAMUS P. COELHO J.C. QUADRADO Instituto Superior de Engenharia de Lisboa Rua Conselheiro Emídio Navarro PORTUGAL Abstract: - A solution to perform
More informationNew LS-DYNA Fluids Solvers
7 th International LS-DYNA Users Conference Fluid/Structure New LS-DYNA Fluids Solvers Grant O. Cook, Jr.. Zeng-Chan Zhang Livermore Software Technology Corporation 7374 Las Positas Road Livermore, CA
More informationarxiv: v1 [cond-mat.dis-nn] 30 Dec 2018
A General Deep Learning Framework for Structure and Dynamics Reconstruction from Time Series Data arxiv:1812.11482v1 [cond-mat.dis-nn] 30 Dec 2018 Zhang Zhang, Jing Liu, Shuo Wang, Ruyue Xin, Jiang Zhang
More informationRunge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh 1
Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured mesh Jun Zhu, inghui Zhong 3, Chi-Wang Shu 4 and Jianxian Qiu 5 Abstract In this paper we generalize a new type
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More informationMultigrid Third-Order Least-Squares Solution of Cauchy-Riemann Equations on Unstructured Triangular Grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids ; : 6 Prepared using fldauth.cls [Version: /9/8 v.] Multigrid Third-Order Least-Squares Solution of Cauchy-Riemann Equations
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationThird Order WENO Scheme on Three Dimensional Tetrahedral Meshes
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 5, No. 2-4, pp. 86-848 Commun. Comput. Phys. February 29 Third Order WENO Scheme on Three Dimensional Tetrahedral Meshes Yong-Tao Zhang 1, and Chi-Wang Shu
More informationCase 3.1: Turbulent Flow over a 2D Multi-Element Airfoil. Summary of Results. Marco Ceze
Case 3.1: Turbulent Flow over a 2D Multi-Element Airfoil Summary of Results Marco Ceze (mceze@umich.edu) 2 nd International Workshop on High-Order CFD Methods, May 27-28, Cologne, Germany C3.1 1/14 Case
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationComparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils
Comparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils Moritz Kompenhans 1, Esteban Ferrer 2, Gonzalo Rubio, Eusebio Valero E.T.S.I.A. (School of Aeronautics)
More informationImage Metamorphosis By Affine Transformations
Image Metamorphosis B Affine Transformations Tim Mers and Peter Spiegel December 16, 2005 Abstract Among the man was to manipulate an image is a technique known as morphing. Image morphing is a special
More informationA Spectral Difference Method for the Euler and Navier-Stokes Equations on Unstructured Meshes
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 2006, Reno, Nevada AIAA 2006-304 A Spectral Difference Method for the Euler and Navier-Stokes Equations on Unstructured Meshes Georg May and
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationPARALLEL DNS USING A COMPRESSIBLE TURBULENT CHANNEL FLOW BENCHMARK
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 ECCOMAS PARALLEL DNS USING A COMPRESSIBLE
More informationAbstract. Euler Equations Solutions for Incompressible Fluid Flow. Copyright A. A. Frempong
Copright A. A. Frempong Euler Equations Solutions for Incompressible Fluid Flow Abstract The simplest solution is usuall the best solution---albert Einstein This paper covers the solutions of the Euler
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationA new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws
HYP 2012, Padova A new multidimensional-type reconstruction and limiting procedure for unstructured (cell-centered) FVs solving hyperbolic conservation laws Argiris I. Delis & Ioannis K. Nikolos (TUC)
More informationHigh Order Schemes for CFD: A Review. Juan Cheng 1. Institute of Applied Physics and Computational Mathematics, Beijing , China.
High Order Schemes for CFD: A Review Juan Cheng 1 Institute of Applied Physics and Computational Mathematics, Beijing 100088, China and Chi-Wang Shu 2 Division of Applied Mathematics, Brown University,
More informationRecent developments for the multigrid scheme of the DLR TAU-Code
www.dlr.de Chart 1 > 21st NIA CFD Seminar > Axel Schwöppe Recent development s for the multigrid scheme of the DLR TAU-Code > Apr 11, 2013 Recent developments for the multigrid scheme of the DLR TAU-Code
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationComputational Fluid Dynamics (CFD) using Graphics Processing Units
Computational Fluid Dynamics (CFD) using Graphics Processing Units Aaron F. Shinn Mechanical Science and Engineering Dept., UIUC Accelerators for Science and Engineering Applications: GPUs and Multicores
More informationObservations on the fifth-order WENO method with non-uniform meshes
Observations on the fifth-order WENO method with non-uniform meshes Rong Wang, a, Hui Feng, b, Raymond J. Spiteri a,, a Department of Computer Science, University of Saskatchewan, Saskatoon, SK, S7N 5C9,
More informationSensitivity analysis of a mesh refinement method using the numerical solutions of 2D lid-driven cavity flow
Author(s): Lal, R. ; Li, Z. Title: Sensitivity analysis of a mesh refinement method using the numerical solutions of D lid-driven cavity flow Journal: Journal of Mathematical Chemistry ISSN: 59-979 Year:
More informationA HYBRID SEMI-PRIMITIVE SHOCK CAPTURING SCHEME FOR CONSERVATION LAWS
Eighth Mississippi State - UAB Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conf. 9 (), pp. 65 73. ISSN: 7-669. URL: http://ejde.math.tstate.edu
More informationtwo-dimensional unsteady incompressible Navier-Stokes equations
*Manuscript revised Click here to view linked References A p-multigrid spectral difference method for two-dimensional unsteady incompressible Navier-Stokes equations Chunlei Liang a,, a Department of Mechanical
More informationDDFV Schemes for the Euler Equations
DDFV Schemes for the Euler Equations Christophe Berthon, Yves Coudière, Vivien Desveaux Journées du GDR Calcul, 5 juillet 2011 Introduction Hyperbolic system of conservation laws in 2D t W + x f(w)+ y
More informationComparison of h- and p- Adaptations for Spectral Difference Methods
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4435 Comparison of h- and p- Adaptations for Spectral Difference Methods Yi Li, Sachin Premasuthan and Antony
More informationWe present a high-order accurate space-time discontinuous Galerkin method for solving two-dimensional
A High-Order Discontinuous Galerkin Method with Unstructured Space-Time Meshes for Two-Dimensional Compressible Flows on Domains with Large Deformations Luming Wang a, Per-Olof Persson a, a Department
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationSpectral Analysis for Fourth Order Problem Using Three Different Basis Functions
Volume 119 o. 1, 09- ISS: 11-9 (on-line version url: http://www.ijpam.eu ijpam.eu Spectral Analsis for Fourth Order Problem Using Three Different Basis Functions 1 Sagitha and Rajeswari Seshadri 1 Department
More informationNumerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods
Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 13: The Lecture deals with:
The Lecture deals with: Some more Suggestions for Improvement of Discretization Schemes Some Non-Trivial Problems with Discretized Equations file:///d /chitra/nptel_phase2/mechanical/cfd/lecture13/13_1.htm[6/20/2012
More informationComputation of Fictitious Gas Flow with Euler Equations
1 Computation of Fictitious Gas Flow with Euler Equations Pei Li and Helmut Sobieczky DLR Göttingen, Germany Abstract The Fictitious Gas Concept supports some computational design methods to construct
More informationA Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids
A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids Ruo Li, Xin Wang and Weibo Zhao January 26, 27 Abstract We proposed an efficient and robust algorithm to solve the steady Euler
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationLiterature Report. Daniël Pols. 23 May 2018
Literature Report Daniël Pols 23 May 2018 Applications Two-phase flow model The evolution of the momentum field in a two phase flow problem is given by the Navier-Stokes equations: u t + u u = 1 ρ p +
More information3.6 Graphing Piecewise-Defined Functions and Shifting and Reflecting Graphs of Functions
76 CHAPTER Graphs and Functions Find the equation of each line. Write the equation in the form = a, = b, or = m + b. For Eercises through 7, write the equation in the form f = m + b.. Through (, 6) and
More informationAeroacoustic computations with a new CFD solver based on the Lattice Boltzmann Method
Aeroacoustic computations with a new CFD solver based on the Lattice Boltzmann Method D. Ricot 1, E. Foquet 2, H. Touil 3, E. Lévêque 3, H. Machrouki 4, F. Chevillotte 5, M. Meldi 6 1: Renault 2: CS 3:
More informationNonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangular Grids for Hyperbolic Systems of Conservation Laws Ivan Christov 1,* Bojan Popov 1 Peter Popov 2 1 Department of Mathematics, 2 Institute for Scientific
More information