Towards Higher-Order Schemes for Compressible Flow
|
|
- Patrick Norton
- 5 years ago
- Views:
Transcription
1 WDS'6 Proceedings of Contributed Papers, Part I, 5 4, 6. ISBN MATFYZPRESS Towards Higher-Order Schemes for Compressible Flow K. Findejs Charles Universit, Facult of Mathematics and Phsics, Prague, Czech Republic. Abstract. Two different approaches to the construction of higher-order schemes for the stationar compressible Euler equations are presented. The former is a generalization of first order Finite Volume Method (FVM) based on a higher-order polnomial reconstruction. The latter one is the so-called Residual Distribution Method (RDM) which is basicall second order, but some generalizations to higher-order are alread known (e.g. [Abgrall, 5], [Caraeni, 5]). We made a comparative stud and present the numerical results for the Ringleb flow problem. Introduction The Finite Volume Method is formall of first order of accurac. This causes that discontinuities in a solution are often smeared and some details are not resolved sufficientl accuratel. To avoid this deficienc, one tries to generalize the FVM to higher-order. One wa is to use a higher-order polnomial reconstruction of piecewise constant approimation used in the FVM. This works onl if sufficientl large stencils are used and the special attention has to be paid to the reconstruction in the vicinit of discontinuities. Another ver popular approach is the Discontinuous Galerkin Method. Similarl to the FVM it uses discontinuous piecewise polnomial data representation and it turns out to be ver accurate and robust method. Its onl drawback is the fact, that the computational and memor costs are ver high. The Residual Distribution Method, which is considered in this paper, uses continuous piecewise polnomial data representation. Therefore, the number of degrees of freedom is lower than in the Discontinuous Galerkin Method. The results of man authors show that this is ver powerful method even for problems with discontinuous solution. Formulation of the problem Let us consider the flow of an inviscid gas in a bounded domain Ω IR N where N = or for D or D flow, respectivel. We suppose that Ω is polgonal in D or polhedral in D, respectivel. Further, we suppose that the flow is adiabatic and we neglect the outer volume force. Our goal is to find stationar solution of the Euler equations equipped with the initial condition w t + N s= with a given vector function w and boundar conditions f s (w) s = in Ω (, + ) () w(, ) = w (), Ω, () B(w(, t)) = for(, t) Ω (, + ). () Here B is suitable boundar operator. The specifications of the boundar conditions and their approimation can be found in [Feistauer et al., ]. The state vector w = (ρ, ρv,...,ρv N, E) T IR m, m = N+ (i.e. m = 4 or 5 for D or D flow, respectivel). Here ρ, v,..., v N and E are phsical quantities denoting the densit, the velocit components and the total energ, respectivel. The flues f s = (ρv s, ρv v s + δ s p,..., ρv N v s + δ Ns p, (E + p)v s ) are m-dimensional mappings of w. Here p denotes the pressure and δ is Kronecker s smbol. Further we denote b A s (w) = Df s (w)/dw the Jacobi matrices of the flues f s, which pla ver important role in the design of presented numerical schemes. Their forms and properties can be found e.g. in [Feistauer et al., ]. 5
2 Finite volume method The sstem () () can be discretized b an eplicit finite volume scheme w k+ i = w k i τ k H(w k i, w k D i j,n ij ) Γ ij, D i D h, t k [, T). (4) j S(i) This scheme is based on the construction of a finite volume mesh D h = {D i } i J, where J Z + = {,,...} is an inde set and h >. Here D i, i J, are closed polgons or polhedrons, if N = or, respectivel, with mutuall disjoint interiors such that Ω = i J D i. If two elements D i, D j D h, i j, contain a nonempt common open face, we call them neighbours. For two neighbours D i, D j D h we set Γ ij = D i D j = Γ ji. (5) We denote b D i or Γ ij the N-dimensional measure of D i or the (N )-dimensional measure of Γ ij, respectivel and b S(i) J the inde set of all neighbours of volume D i. We construct a partition = t < t <... of the time interval [, T] and denote b τ k = t k+ t k the time step between t k and t k+. Scheme (4) represents the algorithm for finding the approimations w k i of the integral averages D i w(, t k )d/ D i of the quantit w over the finite volume D i at time instant t k. Further, we approimate the flu N s= f s(w)(n ij ) s of the quantit w through the face Γ ij in the direction n ij with the aid of a numerical flu H(w k i, wk j, n ij), depending on the value of the approimate solution w k i on the finite volume D i, the value w k j on D j, and on the normal n ij : tk+ t k Γ ij s= n s f s (w) ds dt τ k H(w k i, wk j, n ij) Γ ij. (6) The strateg for finding the stationar solution will be following. We will iterate using the scheme (4) unless the following stopping criterion is satisfied: ρ k i ρk+ i τ k ρ k < TOL, (7) i i J where ρ k i is densit corresponding to the state wk i and TOL is the prescribed tolerance. Etension to higher-order at stead state To increase the order of the FVM, we use a higher-order polnomial reconstruction. For the scheme of order p we define vectors of (p )-degree polnomials ŵ k i P p satisfing ŵ k i D i () d = wk i, (8) D i w(, t k ) Di = ŵ k i () + O(h p ), if w C p (Ω). (9) In the higher-order approach the space integral in (6) is approimated using the quadrature rule with integration points µ and weights α µ tk+ τ k Γ ij t k Γ ij K n s f s (w) ds dt Ĥ(wk i, wk j, n µ ij) = α µ H(ŵ k i ( µ), ŵ k j ( µ), n ij ). () s= µ= The numerical flu H in scheme (4) is then replaced b higher-order flu scheme in space. Ĥ to obtain higher-order Residual distribution method In this section we present a different approach. We define the approimate solution w h of problem () () as a continuous piecewise linear mapping. We denote b P i the nodes of N-simplicial mesh and b w k i the approimate solution at the node P i at the time instant t k. First the cell residual Φ D on the element D will be computed as Φ D = n s f s (w h ) ds D s= s= w h A s ( w) D s d = D s= A s ( w) w h s, () D 6
3 where n s are components of outer unit normal to the boundar D and w is average state obtained using Roe s parameter vector z = ρ(, v...v N, H) T, where H = (E + p)/ρ is the enthalp. It is well known that both state vectors w = w(z) and flues f s (w) = f s ( w(z)) = f s (z) are quadratic functions P i D z i, P i are the vertices of N-simple D of the components of z. Then w = w( z), where z = N+ and z i are computed from data values w i. According to the scalar case, the cell residual of each component of w h D is distributed to the downstream nodes of each element D using the matrices B i which have to satisf the consistenc condition P i D BD i = I: Φ D i = B D i ΦD () The use of Euler eplicit time integration leads to the following scheme which is formall scheme of the tpe (4), but the difference is that w k i are now data values in nodes P i: w k+ i = w k i τ k D i D j,p i D j Φ Dj i. () Here Di = N+ D j,p i D j D j is the area of the median dual cell around the node P i and τ k = t k+ t k is the time step. Further, we will be concerned with the definition of the distribution matrices B D i. In what follows the superscript D will be omitted and the nodes of the element D will be indeed locall b numbers,..., N +. The construction of the matrices B i is direct generalization of the scalar case. For details see e.g. [Deconinck et al., 5]. First we define the matrices K i = N (n i ) s A s ( w), (4) s= where n i is the inward scaled normal to the side of N-simple D opposite to the node P i and its magnitude equals the length of the side. A ver useful propert is that these matrices are diagonalizable: K i = RΛ\R. (5) Using this, we can split the matrices K i into the positive and negative part: K + i = RΛ\ + R, K i = RΛ\ R, (6) where the matrices Λ\ +, Λ\ contain positive or negative eigenvalues, respectivel. Several was are known to define the suitable distribution matrices B i and the properties of the resulting scheme strongl depend on them. One possibilit is to construct the matrices B i in order to get desired second order. This leads to so-called LDA scheme (low diffusion advective), which uses following definition of nodal residuals: Φ LDA i = B LDA i Φ D, where B LDA i = K + i ( N+ i= K + i ). (7) On the other hand it is known from Godunov s theorem (see e.g. [Deconinck et al., 5]) that a linear scheme cannot be both positive and more than first order. This means that LDA scheme cannot be positive which causes a bad behaviour of the solution near discontinuities. Another useful choice is the so-called N scheme - the linear monotone scheme satisfing positivit and therefore (in view of Godunov s theorem) not reaching the second order. The distribution matrices B i are now dependent also on the solution w k i. The are not eplicitl defined, but we can easil epress the nodal residuals directl as Φ N i = K + i N+ j= K + j N+ K j (w i w j ). (8) Our goal is to construct a positive scheme which is satisfactor accurate. The idea is to combine previous two linear schemes into a non-linear LDA-N blended scheme whose nodal residuals are epressed as j= Φ i = LΦ N i + (I L)Φ LDA i, (9) 7
4 where L is the diagonal matri with the diagonal values l ii. This matri should be defined in such a wa that the scheme () behaves like N scheme near discontinuities and like LDA scheme elsewhere. Some heuristic was are known, but the optimal definition of such a matri L is still an open task. As a stopping criterion for the iteration process we again use (7) as in Finite Volume Method where J is the inde set of nodes. Numerical eperiment For a comparison of Residual Distribution Schemes and Finite Volume Schemes we have chosen the so-called Ringleb problem. Its eact stationar solution is smooth and can be epressed analticall. It is depicted in Fig.. It is widel used for testing higher-order methods, because the first order methods do not give satisfactor results. Four numerical schemes were compared: first order Finite Volume Scheme, second order Finite Volume Scheme with least squares linear reconstruction (FVM-LS), N scheme and LDA scheme on four different meshes. As numerical flu for FVM the HLLCF lu was used. The least squares reconstruction (LS) was constructed on stencil containing all volumes D j having at least one verte with D i in common. Let us denote Ŝ(i) the set of indees of such volumes. The reconstructed gradient ŵk i was then computed as the least squares solution of the following sstem of linear algebraic equations ŵ k i ( i j ) = w k i wk j, for j Ŝ(i), () where i and j are centers of gravit of the volumes D i and D j, respectivel. The reconstructed linear function is then ŵ k i () = w k i + ŵ k i ( i ). () All computations were done unless the condition (7) with tolerance TOL = 6 for stationar error was satisfied. For all methods were used the same constant time steps restricted with respect to used mesh to satisf the CFL stabilit condition of Finite Volume Method. The errors in L norm were computed using three point quadrature rule of second order Figure. Computational mesh with nodes (left) and densit isolines of the eact solution (right) 8
5 FVM scheme FVM - LS scheme N scheme LDA scheme Figure. Densit isolines in the domain Ω.95.9 FVM scheme FVM-LS scheme N scheme LDA scheme eact line FVM scheme FVM-LS scheme N scheme LDA scheme Figure. Densit distribution along the boundar (left) and histor of convergence to stationar solution (right) 9
6 Table. Comparison of errors in Densit Discussion method number L norm L norm L norm L norm of elements error order error order FVM scheme FVM-LS scheme N scheme LDA scheme The results in Tab. show that using Finite Volume Method with Least Squares reconstruction we obtained a more accurate solution than using second order LDA scheme. Further we see from Fig. and that the Residual Distribution Scheme suffers from numerical boundar conditions. This causes the loss of accurac close to inlet. Farther from the boundar the scheme gives quite satisfactor results. The convergence histor depicted in Fig. on the right-hand side shows, that the reconstruction used in the Finite Volume Scheme causes a slowdown of convergence to stationar solution. The second order LDA scheme converges much faster. Conclusion From these results we can conclude that RDM is probabl more sensitive to the treatment of boundar conditions near curved boundaries. However, in view of convergence to a stationar solution the RDM seems to be more suitable method than FVM with reconstruction. Acknowledgments. The work is a part of the research project MSM 689 financed b MŠMT and was partl supported b the grant No. 4/5/B-MAT/MFF of the Grant Agenc of the Charles Universit. The author acknowledges this support. References Abgrall R., Construction of high order residual distribution schemes for scalar problems preliminar results for sstems, 4th VKI Lecture Series CFD, von Karman Institute, 5. Abgrall R. and Barth T., Residual distribution schemes for conservation laws via adaptive quadrature, SIAM J. Sci. Comput., Vol. 4, No., pp ,. Caraeni D., A Third order residual distribution methods for stead/unstead simulations: formulation and benchmarking including LES, 4th VKI Lecture Series CFD, von Karman Institute, 5. Deconinck H., Abgrall R., Introduction to residual distribution methods, 4th VKI Lecture Series CFD, von Karman Institute, 5. De Palma P. Pascazio G., Rossiello G., Napolitano M., Implicit third-order-acurate residual distribution schemes for unstead hperbolic problems, 4th VKI Lecture Series CFD, von Karman Institute, 5. Feistauer M., Felcman J., Straškraba I. Mathematical and Computational Methods for Compressible Flow, Oford Universit Press, Oford,. Felcman J. and Findejs K., Higher-Order finite volume method for multidimensional compressible flow, Proceedings of ICNAAM 5 (ed. T. E. Simons), WILEY-VCH, Weinheim, 86-89, 5. Nishikawa H., and Roe P. L., Towards high-order fluctuation-splitting schemes for Navier-Stokes equations, Department of Aerospace engineering, The Universit of Michigan, MI
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM 1
RANGE CONTROL MPC APPROACH FOR TWO-DIMENSIONAL SYSTEM Jirka Roubal Vladimír Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Karlovo náměstí
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains: An Initial Investigation
Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains: An Initial Investigation Samira Nikkar and Jan Nordström Abstract A time-dependent coordinate
More informationDivergence Formulation of Source Term
Preprint accepted for publication in Journal of Computational Physics, 2012 http://dx.doi.org/10.1016/j.jcp.2012.05.032 Divergence Formulation of Source Term Hiroaki Nishikawa National Institute of Aerospace,
More informationFinite volume method on unstructured grids
Finite volume method on unstructured grids Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationBASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS
BASE VECTORS FOR SOLVING OF PARTIAL DIFFERENTIAL EQUATIONS J. Roubal, V. Havlena Department of Control Engineering, Facult of Electrical Engineering, Czech Technical Universit in Prague Abstract The distributed
More informationAnti-diffusive Methods for Hyperbolic Heat Transfer
Anti-diffusive Methods for Hperbolic Heat Transfer Wensheng Shen and Leigh Little Department of Computational Science, SUNY Brockport Brockport, NY 44, USA Liangjian Hu Department of Applied Mathematics,
More informationEQUADIFF 9. Mária Lukáčová-Medviďová Numerical solution of compressible flow. Terms of use:
EQUADIFF 9 Mária Lukáčová-Medviďová Numerical solution of compressible flow In: Zuzana Došlá and Jaromír Kuben and Jaromír Vosmanský (eds.): Proceedings of Equadiff 9, Conference on Differential Equations
More informationRoe Scheme in Generalized Coordinates; Part I- Formulations
Roe Scheme in Generalized Coordinates; Part I- Formulations M.. Kermani. G. Plett Department Of Mechanical & Aerospace ngineering Carleton Universit, Ottawa, Canada Abstract. Most of the computations b
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More informationHigh-order finite difference methods with subcell resolution for 2D detonation waves
Center for Turbulence Research Annual Research Briefs 7 High-order finite difference methods with subcell resolution for D detonation waves B W. Wang, C.-W. Shu, H. C. Yee AND B. Sjögreen. Motivation and
More informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationKrylov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations
DOI.7/s95-6-6-7 Krlov Integration Factor Method on Sparse Grids for High Spatial Dimension Convection Diffusion Equations Dong Lu Yong-Tao Zhang Received: September 5 / Revised: 9 March 6 / Accepted: 8
More informationOn the Extension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models
On the Etension of Goal-Oriented Error Estimation and Hierarchical Modeling to Discrete Lattice Models J. T. Oden, S. Prudhomme, and P. Bauman Institute for Computational Engineering and Sciences The Universit
More informationLecture 4.2 Finite Difference Approximation
Lecture 4. Finite Difference Approimation 1 Discretization As stated in Lecture 1.0, there are three steps in numerically solving the differential equations. They are: 1. Discretization of the domain by
More informationDATA STRUCTURES FOR THE NONLINEAR SUBGRID STABILIZATION METHOD
DATA STRUCTURES FOR TE NONLINEAR SUBGRID STABILIZATION METOD Isaac P. Santos Lucia Catabriga Departamento de Informática Universidade Federal do Espírito Santo - UFES Av. Fernando Ferrari, 54 - Vitória,
More informationMMJ1153 COMPUTATIONAL METHOD IN SOLID MECHANICS PRELIMINARIES TO FEM
B Course Content: A INTRODUCTION AND OVERVIEW Numerical method and Computer-Aided Engineering; Phsical problems; Mathematical models; Finite element method;. B Elements and nodes, natural coordinates,
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationAN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES *
Journal of Computational Mathematics, Vol.27, No.2-3, 29, 388 399. AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES * Yin Yang Hunan Ke Laborator for Computation and Simulation
More informationThe Control-Volume Finite-Difference Approximation to the Diffusion Equation
The Control-Volume Finite-Difference Approimation to the Diffusion Equation ME 448/548 Notes Gerald Recktenwald Portland State Universit Department of Mechanical Engineering gerr@mepdedu ME 448/548: D
More informationPHYSICAL REVIEW E 78, Lattice Boltzmann method for simulation of compressible flows on standard lattices
PHYSICAL REVIEW E 7, 016704 00 Lattice Boltzmann method for simulation of compressible flows on standard lattices Nikolaos I. Prasianakis 1, * and Ilia V. Karlin 1,, 1 Aerothermochemistr and Combustion
More informationVelocity Limit in DPD Simulations of Wall-Bounded Flows
Velocit Limit in DPD Simulations of Wall-Bounded Flows Dmitr A. Fedosov, Igor V. Pivkin and George Em Karniadakis Division of Applied Mathematics, Brown Universit, Providence, RI 2912 USA Abstract Dissipative
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT -23 Numerical Methods For Industrial Problems With Phase Changes Ibrahim, F. J. Vermolen, C. Vuik ISSN 389-652 Reports of the Department of Applied Mathematical Analsis
More informationQUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS
Journal of the Operations Research Societ of Japan 008, Vol. 51, No., 191-01 QUADRATIC AND CONVEX MINIMAX CLASSIFICATION PROBLEMS Tomonari Kitahara Shinji Mizuno Kazuhide Nakata Toko Institute of Technolog
More informationModule 1 : The equation of continuity. Lecture 4: Fourier s Law of Heat Conduction
1 Module 1 : The equation of continuit Lecture 4: Fourier s Law of Heat Conduction NPTEL, IIT Kharagpur, Prof. Saikat Chakrabort, Department of Chemical Engineering Fourier s Law of Heat Conduction According
More informationCOMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER PARABOLIC EQUATIONS. Lili Ju. Xinfeng Liu.
DISCRETE AND CONTINUOUS doi:13934/dcdsb214191667 DYNAMICAL SYSTEMS SERIES B Volume 19, Number 6, August 214 pp 1667 1687 COMPACT IMPLICIT INTEGRATION FACTOR METHODS FOR A FAMILY OF SEMILINEAR FOURTH-ORDER
More informationA FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS NUMBER FLUID FLOWS
Mathematical Models and Methods in Applied Sciences Vol. 12, No. 4 (2002) 525 539 c World Scientific Publishing Compan A FLOW-CONDITION-BASED INTERPOLATION MIXED FINITE ELEMENT PROCEDURE FOR HIGHER REYNOLDS
More informationSTABILIZED FEM SOLUTIONS OF MHD EQUATIONS AROUND A SOLID AND INSIDE A CONDUCTING MEDIUM
Available online: March 09, 2018 Commun. Fac. Sci. Univ. Ank. Ser. A1 Math. Stat. Volume 68, Number 1, Pages 197 208 (2019) DOI: 10.1501/Commua1_0000000901 ISSN 1303 5991 http://communications.science.ankara.edu.tr/index.php?series=a1
More information7. Choice of approximating for the FE-method scalar problems
7. Choice of approimating for the FE-method scalar problems Finite Element Method Differential Equation Weak Formulation Approimating Functions Weighted Residuals FEM - Formulation Weak form of heat flow
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 1. Recap: SDE. Euler Maruyama. Lecture 2, Part 1: Euler Maruyama
Numerical Simulation of Stochastic Differential Equations: Lecture, Part 1 Des Higham Department of Mathematics Universit of Strathclde Lecture, Part 1: Euler Maruama Definition of Euler Maruama Method
More information15th AIAA Computational Fluid Dynamics Conference June 2001 Anaheim, California
AIAA 2 2623 Multi-Dimensional Upwind Constrained Transport on Unstructured Grids for Shallow Water Magnetohdrodnamics Hans De Sterck,2,3 von Karman Institute for Fluid Dnamics, Sint-Genesius-Rode, Belgium
More informationMA8502 Numerical solution of partial differential equations. The Poisson problem: Mixed Dirichlet/Neumann boundary conditions along curved boundaries
MA85 Numerical solution of partial differential equations The Poisson problem: Mied Dirichlet/Neumann boundar conditions along curved boundaries Fall c Einar M. Rønquist Department of Mathematical Sciences
More informationApplications of Proper Orthogonal Decomposition for Inviscid Transonic Aerodynamics
Applications of Proper Orthogonal Decomposition for Inviscid Transonic Aerodnamics Bui-Thanh Tan, Karen Willco and Murali Damodaran Abstract Two etensions to the proper orthogonal decomposition (POD) technique
More informationAE/ME 339. K. M. Isaac Professor of Aerospace Engineering. December 21, 2001 topic13_grid_generation 1
AE/ME 339 Professor of Aerospace Engineering December 21, 2001 topic13_grid_generation 1 The basic idea behind grid generation is the creation of the transformation laws between the phsical space and the
More informationChapter 11 Three-Dimensional Stress Analysis. Chapter 11 Three-Dimensional Stress Analysis
CIVL 7/87 Chapter - /39 Chapter Learning Objectives To introduce concepts of three-dimensional stress and strain. To develop the tetrahedral solid-element stiffness matri. To describe how bod and surface
More informationChapter 2 Basic Conservation Equations for Laminar Convection
Chapter Basic Conservation Equations for Laminar Convection Abstract In this chapter, the basic conservation equations related to laminar fluid flow conservation equations are introduced. On this basis,
More informationWell-Balanced Discontinuous Galerkin Methods for the Euler Equations Under Gravitational Fields
J Sci Comput 6 67:49 5 DOI.7/s95-5-9-5 Well-Balanced Discontinuous Galerkin Methods for the Euler Equations Under Gravitational Fields Gang Li Yulong Xing Received: April 5 / Revised: Jul 5 / Accepted:
More informationNUMERICAL SIMULATION OF DISLOCATION DYNAMICS BY MEANS OF PARAMETRIC APPROACH
Proceedings of Czech Japanese Seminar in Applied Mathematics August -7,, Czech Technical Universit in Prague http://geraldine.fjfi.cvut.cz pp. 1 13 NUMERICAL SIMULATION OF DISLOCATION DYNAMICS BY MEANS
More informationTWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM
Proceedings of ALGORITMY 212 pp. 85 94 TWO DIMENSIONAL SIMULATION OF FLUID-STRUCTURE INTERACTION USING DGFEM JAROSLAVA HASNEDLOVÁ-PROKOPOVÁ,, MILOSLAV FEISTAUER, JAROMÍR HORÁČEK, ADAM KOSÍK, AND VÁCLAV
More informationLecture 5) Eigenstructure and Approximate Riemann Solvers for Hyperbolic Conservation Laws
ecture 5) Eigenstructure and Approimate iemann Solvers for Hperbolic Conservation aws B Prof. Dinshaw S. Balsara (dbalsara@nd.edu) es Houches Summer School in Computational Astrophsics http://www.nd.edu/~dbalsara/numerical-pde-course
More informationApplications of Gauss-Radau and Gauss-Lobatto Numerical Integrations Over a Four Node Quadrilateral Finite Element
Avaiable online at www.banglaol.info angladesh J. Sci. Ind. Res. (), 77-86, 008 ANGLADESH JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH CSIR E-mail: bsir07gmail.com Abstract Applications of Gauss-Radau
More informationarxiv: v1 [physics.comp-ph] 11 Dec 2017
A Compact Fourth-order Gas-kinetic Scheme for the Euler and Navier-Stokes Solutions Xing Ji a, Liang Pan c, Wei Sh b, Kun Xu a,b, arxiv:72.363v [phsics.comp-ph] Dec 27 a Department of Mathematics, Hong
More informationFunctions of Several Variables
Chapter 1 Functions of Several Variables 1.1 Introduction A real valued function of n variables is a function f : R, where the domain is a subset of R n. So: for each ( 1,,..., n ) in, the value of f is
More informationProper Orthogonal Decomposition Extensions For Parametric Applications in Transonic Aerodynamics
Proper Orthogonal Decomposition Etensions For Parametric Applications in Transonic Aerodnamics T. Bui-Thanh, M. Damodaran Singapore-Massachusetts Institute of Technolog Alliance (SMA) School of Mechanical
More informationSolving Variable-Coefficient Fourth-Order ODEs with Polynomial Nonlinearity by Symmetric Homotopy Method
Applied and Computational Mathematics 218; 7(2): 58-7 http://www.sciencepublishinggroup.com/j/acm doi: 1.11648/j.acm.21872.14 ISSN: 2328-565 (Print); ISSN: 2328-5613 (Online) Solving Variable-Coefficient
More informationA high order Padé ADI method for unsteady convection-diffusion equations
Center for Turbulence Research Annual Research Briefs 2005 85 A high order Padé ADI method for unstead convection-diffusion equations B D. You 1. Motivation and objectives The unstead convection-diffusion
More informationConvex ENO Schemes for Hamilton-Jacobi Equations
Convex ENO Schemes for Hamilton-Jacobi Equations Chi-Tien Lin Dedicated to our friend, Xu-Dong Liu, notre Xu-Dong. Abstract. In one dimension, viscosit solutions of Hamilton-Jacobi (HJ equations can be
More informationOn discontinuity capturing methods for convection diffusion equations
On discontinuity capturing methods for convection diffusion equations Volker John 1 and Petr Knobloch 2 1 Universität des Saarlandes, Fachbereich 6.1 Mathematik, Postfach 15 11 50, 66041 Saarbrücken, Germany,
More informationAn improved ADI-DQM based on Bernstein polynomial for solving two-dimensional convection-diffusion equations
Mathematical Theor and Modeling ISSN -50 (Paper) ISSN 5-05(Online) Vol 3, No.1, 013 An improved ADI-DQM based on Bernstein polnomial for solving two-dimensional convection-diffusion equations A.S.J. Al-
More informationHFVS: An Arbitrary High Order Flux Vector Splitting Method
HFVS: An Arbitrary High Order Flu Vector Splitting Method Yibing Chen, Song Jiang and Na Liu Institute of Applied Physics and Computational Mathematics, P.O. Bo 8009, Beijing 00088, P.R. China E-mail:
More informationQUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS
arxiv:1803.0591v1 [math.gm] QUALITATIVE ANALYSIS OF DIFFERENTIAL EQUATIONS Aleander Panfilov stable spiral det A 6 3 5 4 non stable spiral D=0 stable node center non stable node saddle 1 tr A QUALITATIVE
More informationUNIT 4 HEAT TRANSFER BY CONVECTION
UNIT 4 HEAT TRANSFER BY CONVECTION 4.1 Introduction to convection 4. Convection boundar laers 4..1 Hdrodnamic boundar laer over flat plate 4.. Thermal boundar laer 4..3 Concentration boundar laer 4.3 Dimensional
More informationA Meshless Simulations for 2D Nonlinear Reaction-diffusion Brusselator System
Copright 3 Tech Science Press CMES, vol.95, no.4, pp.59-8, 3 A Meshless Simulations for D Nonlinear Reaction-diffusion Brusselator Sstem Ahmad Shirzadi, Vladimir Sladek, Jan Sladek 3 Abstract: This paper
More informationAnalysis and Implementation of Recovery-Based Discontinuous Galerkin for Diffusion 1 Marcus Lo 2 and Bram van Leer 3
9th AIAA Computational Fluid Dnamics - 5 June 009, San Antonio, Texas AIAA 009-3786 Analsis and Implementation of Recover-Based Discontinuous Galerkin for Diffusion Marcus Lo and Bram van Leer 3 Department
More informationEstimation Of Linearised Fluid Film Coefficients In A Rotor Bearing System Subjected To Random Excitation
Estimation Of Linearised Fluid Film Coefficients In A Rotor Bearing Sstem Subjected To Random Ecitation Arshad. Khan and Ahmad A. Khan Department of Mechanical Engineering Z.. College of Engineering &
More informationCHAPTER 2: Partial Derivatives. 2.2 Increments and Differential
CHAPTER : Partial Derivatives.1 Definition of a Partial Derivative. Increments and Differential.3 Chain Rules.4 Local Etrema.5 Absolute Etrema 1 Chapter : Partial Derivatives.1 Definition of a Partial
More informationInteraction of Incompressible Fluid and Moving Bodies
WDS'06 Proceedings of Contributed Papers, Part I, 53 58, 2006. ISBN 80-86732-84-3 MATFYZPRESS Interaction of Incompressible Fluid and Moving Bodies M. Růžička Charles University, Faculty of Mathematics
More informationTHE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY. V. A. Sharifulin.
THE HEATED LAMINAR VERTICAL JET IN A LIQUID WITH POWER-LAW TEMPERATURE DEPENDENCE OF DENSITY 1. Introduction V. A. Sharifulin Perm State Technical Universit, Perm, Russia e-mail: sharifulin@perm.ru Water
More informationFirst, Second, and Third Order Finite-Volume Schemes for Diffusion
First, Second, and Third Order Finite-Volume Schemes for Diffusion Hiro Nishikawa 51st AIAA Aerospace Sciences Meeting, January 10, 2013 Supported by ARO (PM: Dr. Frederick Ferguson), NASA, Software Cradle.
More informationAn Error Minimized Pseudospectral Penalty Direct Poisson Solver
An Error Minimized Pseudospectral Penalt Direct Poisson Solver b Horng Tz-Leng Dept. of Applied Mathematics Feng Chia Universit Taichung Taiwan tlhorng123@gmail.com http://newton.math.fcu.edu.tw/~tlhorng
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationA NEW APPROACH TO TREAT THE RANS-LES
ETMM: th International ERCOFTAC Smposium on Engineering, Marbella, 7-9 September 4 A NEW APPROACH TO TREAT THE RANS-LES INTERFACE IN PANS Lars Davidson Div. of Fluid Dnamics, Dept. of Applied Mechanics
More informationON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
Conference Applications of Mathematics 212 in honor of the 6th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 212 ON USING ARTIFICIAL COMPRESSIBILITY METHOD FOR SOLVING TURBULENT FLOWS
More informationc 1999 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol., No. 6, pp. 978 994 c 999 Societ for Industrial and Applied Mathematics A STUDY OF MONITOR FUNCTIONS FOR TWO-DIMENSIONAL ADAPTIVE MESH GENERATION WEIMING CAO, WEIZHANG HUANG,
More informationDiscontinuous Galerkin method for a class of elliptic multi-scale problems
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 000; 00: 6 [Version: 00/09/8 v.0] Discontinuous Galerkin method for a class of elliptic multi-scale problems Ling Yuan
More informationS-I-R Model with Directed Spatial Diffusion
S-I-R Model with Directed Spatial Diffusion Fabio A. Milner Ruijun Zhao Department of Mathematics, Purdue Universit, West Lafaette, IN 4797-267 Abstract A S-I-R epidemic model is described in which susceptible
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 59 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Finite Volume Method These slides are partially based on the recommended textbook: Culbert B.
More informationNumerical Computation of Unsteady Compressible Flow past a Plunging Aerofoil
Journal of Applied Fluid Mechanics, Vol. 5, No. 4, pp.131-138, 2012. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. Numerical Computation of Unstead Compressible Flow past a Plunging
More informationA Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations
Int. Conference on Boundary and Interior Layers BAIL 6 G. Lube, G. Rapin Eds c University of Göttingen, Germany, 6 A Discontinuous Galerkin Moving Mesh Method for Hamilton-Jacobi Equations. Introduction
More informationFully Discrete Energy Stable High Order Finite Difference Methods for Hyperbolic Problems in Deforming Domains
Department of Mathematics Full Discrete Energ Stable High Order Finite Difference Methods for Hperbolic Problems in Deforming Domains Samira Nikkar and Jan Nordström LiTH-MAT-R--/--SE Department of Mathematics
More informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationKINEMATIC RELATIONS IN DEFORMATION OF SOLIDS
Chapter 8 KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS Figure 8.1: 195 196 CHAPTER 8. KINEMATIC RELATIONS IN DEFORMATION OF SOLIDS 8.1 Motivation In Chapter 3, the conservation of linear momentum for a
More informationSimultaneous Orthogonal Rotations Angle
ELEKTROTEHNIŠKI VESTNIK 8(1-2): -11, 2011 ENGLISH EDITION Simultaneous Orthogonal Rotations Angle Sašo Tomažič 1, Sara Stančin 2 Facult of Electrical Engineering, Universit of Ljubljana 1 E-mail: saso.tomaic@fe.uni-lj.si
More informationChapter 2 Overview of the Simulation Methodology
Chapter Overview of the Simulation Methodolog. Introduction This chapter presents an overview of the simulation methodolog that comprises both the art and science involved in simulating phsical phenomena.
More information520 Chapter 9. Nonlinear Differential Equations and Stability. dt =
5 Chapter 9. Nonlinear Differential Equations and Stabilit dt L dθ. g cos θ cos α Wh was the negative square root chosen in the last equation? (b) If T is the natural period of oscillation, derive the
More informationThird-Order Active-Flux Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and Newton Solver
Preprint accepted for publication in Computers and Fluids, 16 http://d.doi.org/1.116/j.compfluid.15.1. Third-Order Active-Flu Scheme for Advection Diffusion: Hyperbolic Diffusion, Boundary Condition, and
More informationResearch Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free Surface Flow
Journal of Applied Mathematics Volume 2, Article ID 793653, 6 pages doi:.55/2/793653 Research Article Development of a Particle Interaction Kernel Function in MPS Method for Simulating Incompressible Free
More informationEnforcing the Discrete Maximum Principle for Linear Finite Element Solutions of Second-Order Elliptic Problems
COMMUNICATIONS IN COMPUTATIONAL PHYSICS Vol. 3, No. 4, pp. 852-877 Commun. Comput. Phs. April 28 Enforcing the Discrete Maimum Principle for Linear Finite Element Solutions of Second-Order Elliptic Problems
More informationAdaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws
Adaptive backward difference formula discontinuous Galerkin finite element method for the solution of conservation laws V. Dolejší and P. Kůs Charles University Prague, Faculty of Mathematics and Physics,
More informationFirst Order Equations
10 1 Linear and Semilinear Equations Chapter First Order Equations Contents 1 Linear and Semilinear Equations 9 Quasilinear Equations 19 3 Wave Equation 6 4 Sstems of Equations 31 1 Linear and Semilinear
More informationResidual Distribution. basics, recents developments, relations with other techniques
Residual Distribution basics, recents developments, relations with other techniques MARIO RICCHIUTO July 19, 2012 INTRODUCTION with historical perspective t u+ F(u) = 0 In the 80 s numerical techniques
More informationDISCONTINUOUS GALERKIN DISCRETIZATION WITH EMBEDDED BOUNDARY CONDITIONS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol.3(23), No., pp.35 58 c 23 Editorial Board of the Journal Computational Methods in Applied Mathematics Joint Co. Ltd. DISCONTINUOUS GALERKIN DISCRETIZATION
More informationAnd similarly in the other directions, so the overall result is expressed compactly as,
SQEP Tutorial Session 5: T7S0 Relates to Knowledge & Skills.5,.8 Last Update: //3 Force on an element of area; Definition of principal stresses and strains; Definition of Tresca and Mises equivalent stresses;
More informationSome recent results on positivity-preserving discretisations for the convection-diffusion equation
Some recent results on positivity-preserving discretisations for the convection-diffusion equation Gabriel R. Barrenechea 1, Volker John 2 & Petr Knobloch 3 1 Department of Mathematics and Statistics,
More information16.5. Maclaurin and Taylor Series. Introduction. Prerequisites. Learning Outcomes
Maclaurin and Talor Series 6.5 Introduction In this Section we eamine how functions ma be epressed in terms of power series. This is an etremel useful wa of epressing a function since (as we shall see)
More informationAnalytic Solutions of Partial Differential Equations
ii Analtic Solutions of Partial Differential Equations 15 credits Taught Semester 1, Year running 3/4 MATH3414 School of Mathematics, Universit of Lee Pre-requisites MATH36 or MATH4 or equivalent. Co-requisites
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationA NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS
A NUMERICAL STUDY FOR THE PERFORMANCE OF THE RUNGE-KUTTA FINITE DIFFERENCE METHOD BASED ON DIFFERENT NUMERICAL HAMILTONIANS HASEENA AHMED AND HAILIANG LIU Abstract. High resolution finite difference methods
More informationUNCORRECTED. To recognise the rules of a number of common algebraic relations: y = x 1 y 2 = x
5A galler of graphs Objectives To recognise the rules of a number of common algebraic relations: = = = (rectangular hperbola) + = (circle). To be able to sketch the graphs of these relations. To be able
More informationConstruction of very high order Residual Distribution Schemes for steady problems
Construction of very high order Residual Distribution Schemes for steady problems Rémi Abgrall, Mario Ricchiuto, Cédric Tavé, Nadège Villedieu and Herman Deconinck Mathématiques Appliquées de Bordeaux,
More informationWavelet-Based Linearization Method in Nonlinear Systems
Wavelet-Based Linearization Method in Nonlinear Sstems Xiaomeng Ma A Thesis In The Department of Mechanical and Industrial Engineering Presented in Partial Fulfillment of the Requirements for the Degree
More informationOn the Missing Modes When Using the Exact Frequency Relationship between Kirchhoff and Mindlin Plates
On the Missing Modes When Using the Eact Frequenc Relationship between Kirchhoff and Mindlin Plates C.W. Lim 1,*, Z.R. Li 1, Y. Xiang, G.W. Wei 3 and C.M. Wang 4 1 Department of Building and Construction,
More informationH. L. Atkins* NASA Langley Research Center. Hampton, VA either limiters or added dissipation when applied to
Local Analysis of Shock Capturing Using Discontinuous Galerkin Methodology H. L. Atkins* NASA Langley Research Center Hampton, A 68- Abstract The compact form of the discontinuous Galerkin method allows
More informationA high order adaptive finite element method for solving nonlinear hyperbolic conservation laws
A high order adaptive finite element method for solving nonlinear hyperbolic conservation laws Zhengfu Xu, Jinchao Xu and Chi-Wang Shu 0th April 010 Abstract In this note, we apply the h-adaptive streamline
More informationA Method for Geometry Optimization in a Simple Model of Two-Dimensional Heat Transfer
A Method for Geometr Optimization in a Simple Model of Two-Dimensional Heat Transfer X. Peng, K. Niakhai B. Protas Jul, 3 arxiv:37.48v [math.na] 4 Jul 3 Abstract This investigation is motivated b the problem
More informationConservation of Linear Momentum for a Differential Control Volume
Conservation of Linear Momentum for a Differential Control Volume When we applied the rate-form of the conservation of mass equation to a differential control volume (open sstem in Cartesian coordinates,
More informationChapter #6 EEE State Space Analysis and Controller Design
hapter #6 EEE3-83 State Space nalsis and ontroller Design State feedack Linear Quadratic Regulator (LQR) Open and closed loop estimators Reduced Order Estimators (not assessed) racking Module Leader: Dr
More informationChapter 3: BASIC ELEMENTS. solid mechanics)
Chapter 3: BASIC ELEMENTS Section 3.: Preliminaries (review of solid mechanics) Outline Most structural analsis FE codes are displacement based In this chapter we discuss interpolation methods and elements
More informationAPPLICATION OF MFCAV RIEMANN SOLVER TO MAIRE S CELL-CENTERED LAGRANGIAN METHOD *
Journal of Computational Mathematics Vol.33, No.2, 25, 28 45. http://www.global-sci.org/jcm doi:.428/jcm.48-m44 APPLICATION OF MFCAV RIEMANN SOLVER TO MAIRE S CELL-CENTERED LAGRANGIAN METHOD * Yan Liu,
More information3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES
3.0 PROBABILITY, RANDOM VARIABLES AND RANDOM PROCESSES 3.1 Introduction In this chapter we will review the concepts of probabilit, rom variables rom processes. We begin b reviewing some of the definitions
More information