GENUINELY MULTIDIMENSIONAL NON-DISSIPATIVE FINITE VOLUME SCHEMES FOR TRANSPORT
|
|
- Clemence Webb
- 5 years ago
- Views:
Transcription
1 GENUINELY MULTIDIMENSIONAL NON-DISSIPATIVE FINITE VOLUME SCHEMES FOR TRANSPORT BRUNO DESPRÉS, FRÉDÉRIC LAGOUTIÈRE Commissariat à l Énergie Atomique, Bruyères-le-Châtel. despres@cmpax.polytechnique.fr despres@cmpax.polytechnique.fr Université Paris-VII and Laboratoire Jacques-Louis Lions, UMR 7598, 75, rue du Chevaleret, 753 Paris lagoutie@math.jussieu.fr Keywords: grids. multidimensional transport, finite volume schemes, anti-dissipative schemes, triangular Abstract We develop a new multidimensional finite volume algorithm for transport equations. This algorithm is both stable and non-dissipative. It is based on a reconstruction of the discrete solution inside each cell, at every time step. The proposed reconstruction, genuinely multidimensional, allows to recover sharp profiles both in the direction of the transport velocity and in the transverse direction. It is an extension of the one-dimensional reconstructions analyzed in Lagoutière 27 and Lagoutière 27. Introduction The present study concerns reconstruction schemes for transport equations. We are especially interested in schemes that are not dissipative, in particular for initial conditions with discontinuities. This leads to consider reconstructions that are not smooth, on the contrary to usual reconstructions. We indeed develop a discontinuous reconstruction scheme which consists in reconstructing the constant-in-cell datum as a discontinuous inside each cell function. Typically, the goal is to develop transport schemes for mass or volume fractions in multi-fluids. These fractions can be discontinuous at interfaces between components or continuous in mixing zones. A former algorithm was already developed in Després and Lagoutière 27. It was based in dimension on the limited downwind scheme equivalent to the Ultra-bee limiter for advection with constant velocity and the multidimensional algorithms where obtained via an alternate direction splitting strategy. Results are satisfying, in particular for inter- This work was partially supported by CEA/DIF Bruyères-le-Châtel faces see results for interface instabilities in dimensions 2 and 3 in the cited reference. The main drawback of this method is precisely the dimensional splitting, which prevents the algorithm form being used on non-cartesian grids. We here present a way to generalize the limited downwind scheme to the case of a general triangular grid. It is based on the geometrical approach followed in Lagoutière 27, Lagoutière 27, that provided a new interpretation of the limited downwind algorithm, in terms of reconstruction scheme. The paper is organized as follows. Sections 2 and 3 present the mathematical and numerical problems and the notations. Then, we recall the reconstruction procedure leading to the limited downwind scheme in dimension section 4. In section 5, the main subject of the paper is addressed: the extension to upper dimensions of the preceding procedure. The presentation deals with dimension 2. Finally, section 6 exhibits numerical results. For the time-being, the most efficient algorithms for pure transport rest upon interface reconstruction: see for example SLIC and Youngs method, Noh and Woodward 976 and Youngs 984, and the method of Mosso Mosso and Cleancy 995 which is a recent promising enhancement. These methods are essentially limited to Cartesian grids. We here try to derive a truly multidimensional reconstruction algorithm.
2 2 Model problem The considered model is the linear transport equation with constant given velocity t ut, x + divaut, x = t ut, x + a x ut, x = for t >, u, = u L R 2, where a is a smooth divergence-free velocity field: divat, x = for every t, x R + R 2. The main issue for the numerical treatment of this PDE problem is the numerical diffusion. This phenomenon, easily understandable in dimension, is much more complex in higher dimension. It shall be decomposed in two different types of diffusion. The diffusion of the first type, which will be called longitudinal diffusion, is the one that occurs in the direction of the velocity. It is the diffusion which is present in one-dimensional classical algorithms. The second type diffusion, transverse diffusion, is typically multidimensional and is due to the fact that the mesh may not be aligned with the velocity. This distinction between to phenomena could appear arbitrary, but is in accordance with the numerical tests. It is illustrated on FIG. 2 and 3 representing numerical solutions obtained with the classical upwind scheme on a square mesh, with, as initial condition, the characteristic function of a square FIG.. The transport velocity is a =, aligned with the mesh for FIG. 2 and a =, diagonal, not aligned with the mesh for FIG. 3 and the boundary condition on [, ] 2 are periodic. The results are displayed for time t = after one period Figure. Initial condition. Figure 2. Upwind scheme. a =,, aligned with the mesh. Time t =. The longitudinal diffusion applies in the direction of transport. There is no transverse diffusion Figure 3. Upwind scheme. a =,, not aligned with the mesh. Time t =. The diffusion which applies orthogonally to the transport direction is brought to light. The difference between longitudinal and transverse diffusions appears clearly in Després and Lagoutière 2 which was a previous attempt to elaborate nondissipative schemes on non-cartesian grids. Based on mono-dimensional techniques, the schemes derived from this work where longitudinally anti-dissipative, but not transversely. The new method here described is based on an interpretation of the anti-dissipative limited downwind scheme which was proposed in Lagoutière 27 previous works are Després and Lagoutière 22 and Després and Lagoutière 27. This paper shows that the limited downwind scheme in dimension can be understood as a reconstruction scheme, decomposed in 3 stages at each time step: a reconstruction stage which, starting from a constant-in-cell datum, constructs a new datum, presenting in each cell one discontinuity that separates two constant values, a transport stage that solves the transport operation with the new datum, a projection stage that computes the mean value of the transported datum in each cell 2
3 see section 4 for precision on the reconstruction stage. The natural extension to the multidimensional transport problem consists in reconstructing the datum in two steps: in a first step, performing a transverse reconstruction, and performing a longitudinal reconstruction in a second step. This is the technique which is proposed in this work. The following describes precisely these two operations. 3 General numerical formalism It is the one of finite volume methods for problem. We consider a mesh of R 2 composed of non-empty open triangles T j j Z such that j Z T j = R 2 and T i T j = for every i, every j i. For each cell T j for j Z, one denotes Kj the set of indices of neighboring cells of T j that is to say cells having a common edge with T j, Kj = { k Z \ {j} s.t. meas T j T k > } where meas denotes the Lebesgue measure in dimension. For j Z and for k Kj T j and T k have an edge in common, one denotes l j,k the length of the common edge, l j,k = meas Tj T k and n j,k the unit normal vector to the common edge outward to T j. We thus have l j,k = l k,j and n j,k = n k,j for every j Z and every k Kj. We then denote by K + j the set of indices of the downwind neighbors of T j and by K j the set of indices of the upwind neighbors of T j : K + j = {k Kj s.t. a,n j,k > }, Let s j denote the surface of cell T j : s j = meas 2 T j where meas 2 denotes the Lebesgue measure in dimension 2. The general form of the here considered schemes is obtained by choosing a time step t > and by integrating the transport equation of on [n t, n + t] T j : u n+ j u n j t + l j,k a n j,k s,n j,ku n j,k + j K j K + j l j,k a n j,k,n j,ku n j,k =. 2 In this formula, the quantities a n j,k are approximate values of the given velocity at, x on the edges and the u n j,k are approximate values of the solution on the edges between times n t and n + t. The upwind scheme is obtained taking u n j,k = un j for k K+ j and u n j,k = un k for k K j. We will propose another definition of these numerical fluxes, intended to provide more precise numerical results. 4 Discontinuous reconstructions in dimension In the following, for every a, b R, a, b denotes the interval [a, b] if a b and the interval [b, a] else, that is to say that we adopt the convention a, b = [mina, b, maxa, b] a K + j = {k Kj s.t. a,n j,k < }. T m T k l j,k n j,k T j Figure 4. Mesh and notations k, l, m Kj k, l K + j m K j T l and one has a, b = b, a. We here recall the principle of discontinuous reconstruction schemes in dimension. Details and proofs lie in Lagoutière 27 and Lagoutière 27. We consider a mesh on R with constant space step x > whose cells are the intervals T j = j /2 x, j + /2 x for j Z. It is shown in Lagoutière 27 that for the transport equation t u + a x u =, the limited downwind scheme out of Després and Lagoutière 22 is equivalent to the 3-stages following scheme. Let u n j be a discrete datum associated j Z with a constant-in-cell function. Reconstruct in each cell the discrete datum in a form with more detail not constant following the algorithm detailed above. 2 Perform the exact transport of this reconstructed datum at velocity a for a time t. 3
4 3 Project the computed solution to obtain a new constant-in-cell datum and to be able to iterate these 3 operations. It now remains to describe the reconstruction operation, the two others being trivial. Let u n j be given. We define the discrete function ux = j Z un j χ [j /2 x,j+/2 xx. The al- j Z gorithm is defined by if u n j / u n j, un j+, that is to say if un j is a local extremum, we do not reconstruct the datum in cell T j, if u n j un j, un j+, we define u n j,l = u n j, u n j,r = un j+, d n j = un j+ un j u n j+, un j and the reconstructed datum is defined inside T j by ûx = u n j,l if x [j /2 x, j /2 + d n j x, ûx = u n j,r if x [j /2 + d n j x, j + /2 x. Notice that when u n j un j, un j+, formula 3 is welldefined. In particular, the denominator in the definition of d n j is not. More precisely, one always has dn j,, so that j /2 + d n j x T j. Furthermore, another essential property of this reconstruction is that it is conservative: d n j un j,l + dn j un j,r = un j. Of course, other values for u n j,l and un j,r are admissible, but this choice leads to the most anti-dissipative scheme under constraints of L stability and decreasing of the total variation, cf. Lagoutière 27. At last, notice that not to reconstruct the datum is equivalent to take u n j,l = un j,r = un j and gives the Godunov upwind scheme. The discontinuous reconstruction is illustrated on FIG Discontinuous reconstructions in dimension 2 The two-dimensional extension is not obvious. The genuine two-dimensional nature of the problem, due to the choice of a non-cartesian grid, leads to the fact that the finite volume method 2 cannot be seen as an exact Godunov one which would be composed of an exact transport step preceded or not by a reconstruction followed by a projection on the mesh since only edge fluxes are taken into account, with the exception of node fluxes. In order to separate the problems of longitudinal and transverse diffusions see the introduction of this paper, we propose to perform the reconstruction in 2 steps, each one being a one-dimensional reconstruction. In order to simplify the presentation, let us assume that a is a time and space constant. Let u n j be a discrete datum at time j Z step n. The first reconstruction, transverse, boils down to cut some cells in two in the direction of the velocity and to modify the value of the datum in each of these two subcells. The cell T j has at least one downwind neighbor and at most two thanks to the divergence-free hypothesis. If it has only one downwind neighbor, we do not perform the transverse reconstruction we do not cut the cell. Let us now assume that T j has two downwind neighbors, T k and T l. It has then one upwind neighbor, T m. We consider the intersection point of the two edges relative to the downwind neighbors and cut T j along the line passing on this intersection point and parallel to a. The two induced subcells are denoted T j,k and T j,l : T j,k has T k as downwind neighbor, and T j,l has T l as downwind neighbor. This cutting is illustrated on FIG. 6. a T k T j,k T j T j,l T l u n j+ u n j u n j u n r j u n l j j j j + dn j x d n j x Figure 5. Discontinuous reconstruction. T m Figure 6. Transverse reconstruction. We denote s j,k and s j,l the surfaces of cells T j,k and T j,l. One of course has s j,k + s j,l = s j and s j,k > and s j,l >. 4
5 We now have to assign to each of the sub-cells a value of the reconstructed solution. The adopted principle is quite similar to the one in dimension, except that the locus of the discontinuity is determined by the geometrical aspect and not the local values of the datum: it is the line parallel to a defined just above. The aim is to define a value u n j,k in the cell T j,k and a value u n j,l in the cell T j,l by maximizing u n j,l un j,k to guarantee the antidissipativity and making certain that s j,k u n j,k +s j,lu n j,l = s j u n j to guarantee the conservativity. On the other side, following the same rules as in dimension, we impose that the triplet {u n l, un j,l, un j,k } has the same monotony as the pair {u n l, un j } and that the triplet {un j,l, un j,k, un k } has the same monotony as the pair {u n j, un k }. These constraints imply in particular that the datum in T j would not be reconstructed if u n j was a local extremum in the transverse direction. The algorithm is if u n j / u n k, un l if un j is a local extremum in the transverse direction, we do not reconstruct u in the cell T j that is to say, u n j,l = un j,k = un j, if u n j un k, un l, if s j u n j s j,lu n l /s j,k [u n j, un k ], we define { u n j,l = un l, u n j,k = s ju n j s j,lu n l /s j,k ; if s j u n j s j,ku n k /s j,l [u n j, un j ], we define { u n j,l = s j u n j s j,k u n k/s j,l, u n j,k = un k. Lemma. The Courant-Friedrichs-Lewy CFL condition on the time step is not degraded by the transverse reconstruction. Proof. The standard condition for the upwind scheme for the cell T j is k K t + j l j,ka, n j,k, s j i.e., in the particular studied case, t l j,ka, n j,k + l j,l a, n j,l s j. 4 One has s j,k = l j l j,k a, n j,k 2 and and s j,l = l j l j,l a, n j,l, 2 s j = s j,k + s j,l = l j l j,k a, n j,k + l j,l a, n j,l 2 The two inequalities of 5 and inequality 4 thus rewrite they are equivalent. 2 t l j, Remark. Once the solution is reconstructed, the transport of it is related to one-dimensional transport since T j,l and T j,k have only one upwind and one downwind neighbor: respectively, T m, T l and T m, T k. When T j has only one downwind neighbor, we do not perform the transverse reconstruction and can consider the problem as monodimensional, performing virtually the cutting but assigning the value u n j to each sub-cell. After this transverse reconstruction, the algorithm is more classical, lying on the fact that the transport is mono-dimensional, as noticed in remark. We thus can use the algorithm of our own choice. In the following, we pitch on the limited downwind one. This can also be understood as a second reconstruction, longitudinal. 6 Numerical results We present a few results obtained with the algorithm described herein. Different velocity fields are used: translation and rotation fields. For all the test-cases, the spatial domain is [, ] 2. The triangular mesh has been generated by the software FreeFem++: see FIG. 7. The CFL condition for sub-cells T j,k and T j,l are t l j,ka, n j,k s j,k and t l j,la, n j,l s j,l. 5 Let us denote l j the length of the segment separating T j,k and T j,l : l j = mes Tj,k T j,l. Figure 7. Example of Freefem++ mesh, used for the numerical tests here, 3766 triangles. 5
6 For all P the test-cases, the Courant number k K max j Z t + j l j,ka,n j,k s j takes the value.. 6. Translation field Here we consider the veloc- translation of a square ity at, x, y = This first test is the translation of the characteristic function of a square, u x, y = χ [.3,.7] 2x, y. The boundary conditions are periodic in x and y. Initial condition and results at time t = after one revolution are reported in FIG. 8 and FIG. 9 reports the result by the upwind scheme Rotation field Rotation of a square at, x, y = The velocity field is here 2πy 2πx and the initial condition is the same as in the preceding test. The final time is t =. Figure presents the result Figure. Rotation of a square with 5874 cells. Zalesak s test-case the preceding test, at, x, y = The velocity field is the same as in 2πy 2πx and the initial condition is taken from the original paper Zalesak 979: u x, y = χ Z x, y Figure 8. Translation of a square with 5874 cells Figure 9. Translation of a square with 5874 cells, upwind scheme. where with Z = C \ R C = { x, y; x y }, D = {x, y; x.5 /4 and y 7/2}. Initial condition and result computed with 5874 cells at time t = lie in FIG.. Then, FIG. 2 represents the same for a mesh made of 2368 cells, and FIG. 3 for a mesh composed of cells. 6
7 which means that at time t = the field is reversed. The initial condition is the characteristic function of a rectangle: u x, y = χ[.2,.7] [.45,.55]x, y The exact solution at time 2 coincides with the initial condition. We observe the initial condition FIG. 4 and the results at time t = FIG. 5 with 5874 cells and finally the result at time t = 2 with 5874 cells FIG. 6, 2368 cells FIG. 7 and cells FIG. 8. The result given by the upwind scheme with cells is reported on FIG Figure. Zalesak s test-case with 5874 cells Figure 4. Initial condition with 5874 cells Figure 2. Zalesak s test-case with 2368 cells Figure 5. Numerical solution with 5874 cells at time t = Figure 3. Zalesak s test-case with cells Deformation of a rectangle The velocity field is now not constant in time:! 3πy.52 a = sgn t, 3πx Figure 6. Numerical solution with 5874 cells at time t = 2. 7
8 References Després B. and Lagoutière F. 2: Generalized Harten formalism and longitudinal variation diminishing schemes for linear advection on arbitrary grids. M2AN Math. Model. Numer. Anal. 35, no. 6, pp Després B. and Lagoutière F. 22: Contact discontinuity capturing schemes for linear advection and compressible gas dynamics. J. Sci. Comput,. 6 2, no. 4: pp Figure 7. Numerical solution with 2368 cells at time t = 2. Després B., Lagoutière F. 27: Numerical resolution of a two-component compressible fluid model with interfaces, to appear in Progress in Computational Fluid Dynamics Lagoutière F. 27: Stability of reconstruction schemes for scalar hyperbolic conservation laws, in revision Lagoutière F. 27: Non-dissipative reconstruction schemes satisfying entropy inequalities, in revision. Mosso S. and Cleancy S. 995: A geometrical derived priority system for Young s interface reconstruction. LA-CP-95-8, LANL report. Figure 8. Numerical solution with cells at time t = 2. Noh W. F. and Woodward P. R. 976: SLIC Simple Line Interface Calculation. Springer Lecture Notes in Physics, 25: pp D. L. Youngs, An interface tracking method for a 3D eulerian hydrodynamics code, Technical Report 44/92/35, AWRE Zalesak S. T. 979: Fully multidimensional fluxcorrected transport algorithms for fluids, Journal of Computational Physics, 3, pp Figure 9. Numerical solution with cells at time t = 2 with the upwind scheme. 7 Final comments We developed a new method for the numerical transport in dimension 2. The approach is truly multidimensional in the sense that is does not resume to a mono-dimensional reconstruction of interfaces. The numerical results show the anti-dissipative behavior of the algorithm. The results are nevertheless not perfect. Indeed, the discontinuity lines may be degraded in long time see FIG., 2, 3. Next study will concern the extension on general meshes non-triangular and in dimension 3, and the application to gas dynamics equations. 8
Development 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 informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationVector Image Polygon (VIP) limiters in ALE Hydrodynamics
New Models and Hydrocodes, Paris, France,24-28 May 2010 Vector Image Polygon (VIP) limiters in ALE Hydrodynamics 6 4 8 3 5 2 79 1 2 1 1 1 Gabi Luttwak 1 and Joseph Falcovitz 2 1 Rafael, P.O. Box 2250,
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationBACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION
BACK AND FORTH ERROR COMPENSATION AND CORRECTION METHODS FOR REMOVING ERRORS INDUCED BY UNEVEN GRADIENTS OF THE LEVEL SET FUNCTION TODD F. DUPONT AND YINGJIE LIU Abstract. We propose a method that significantly
More informationWeno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach
Weno Scheme for Transport Equation on Unstructured Grids with a DDFV Approach Florence Hubert and Rémi Tesson Abstract In this paper we develop a DDFV approach for WENO scheme on unstructred grids for
More informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
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 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 informationIntroduction. D t x = u. (2)
ESAI: PROCEEDINGS, August 28, Vol. 24, p. 46-59 C. Dobrzynsi, P. Frey, Ph. Pebay, Editors LAGRANGIAN ETHOD ENHANCED WITH EDGE SWAPPING FOR THE FREE FALL AND CONTACT PROBLE Étienne Bernard 1, Stéphane Del
More informationSkåne University Hospital Lund, Lund, Sweden 2 Deparment of Numerical Analysis, Centre for Mathematical Sciences, Lund University, Lund, Sweden
Volume Tracking: A New Method for Visualization of Intracardiac Blood Flow from Three-Dimensional, Time-Resolved, Three-Component Magnetic Resonance Velocity Mapping Appendix: Theory and Numerical Implementation
More informationLax-Wendroff and McCormack Schemes for Numerical Simulation of Unsteady Gradually and Rapidly Varied Open Channel Flow
Archives of Hydro-Engineering and Environmental Mechanics Vol. 60 (2013), No. 1 4, pp. 51 62 DOI: 10.2478/heem-2013-0008 IBW PAN, ISSN 1231 3726 Lax-Wendroff and McCormack Schemes for Numerical Simulation
More informationTraffic flow optimization on roundabouts
Available online at www.sciencedirect.com ScienceDirect Procedia - Social and Behavioral Sciences 111 ( 2014 ) 127 136 EWGT2013 16 th Meeting of the EURO Working Group on Transportation Traffic flow optimization
More informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
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 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 informationConservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow. Abstract
Conservative high order semi-lagrangian finite difference WENO methods for advection in incompressible flow Jing-Mei Qiu 1 and Chi-Wang Shu Abstract In this paper, we propose a semi-lagrangian finite difference
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationAdvective and conservative semi-lagrangian schemes on uniform and non-uniform grids
Advective and conservative semi-lagrangian schemes on uniform and non-uniform grids M. Mehrenberger Université de Strasbourg and Max-Planck Institut für Plasmaphysik 5 September 2013 M. Mehrenberger (UDS
More informationA Toolbox of Level Set Methods
A Toolbox of Level Set Methods Ian Mitchell Department of Computer Science University of British Columbia http://www.cs.ubc.ca/~mitchell mitchell@cs.ubc.ca research supported by the Natural Science and
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
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 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 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 informationA-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain
2th Annual CFD Symposium, August 9-1, 218, Bangalore A-posteriori Diffusion Analysis of Numerical Schemes in Wavenumber Domain S. M. Joshi & A. Chatterjee Department of Aerospace Engineering Indian Institute
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 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 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 informationc 2005 Society for Industrial and Applied Mathematics
SIAM J. SCI. COMPUT. Vol. 6, No. 3, pp. 785 89 c 5 Society for Industrial and Applied Mathematics A HIGH-RESOLUTION ROTATED GRID METHOD FOR CONSERVATION LAWS WITH EMBEDDED GEOMETRIES CHRISTIANE HELZEL,
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 informationarxiv: v1 [math.na] 20 Sep 2016
arxiv:1609.06236v1 [math.na] 20 Sep 2016 A Local Mesh Modification Strategy for Interface Problems with Application to Shape and Topology Optimization P. Gangl 1,2 and U. Langer 3 1 Doctoral Program Comp.
More informationAn efficient filtered scheme for some first order Hamilton-Jacobi-Bellman equat
An efficient filtered scheme for some first order Hamilton-Jacobi-Bellman equations joint work with Olivier Bokanowski 2 and Maurizio Falcone SAPIENZA - Università di Roma 2 Laboratoire Jacques-Louis Lions,
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 informationStability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow
Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow F. Haider 1, B. Courbet 1, J.P. Croisille 2 1 Département de Simulation Numérique des Ecoulements
More informationHomework 4A Due November 7th IN CLASS
CS207, Fall 2014 Systems Development for Computational Science Cris Cecka, Ray Jones Homework 4A Due November 7th IN CLASS Previously, we ve developed a quite robust Graph class to let us use Node and
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationSIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS
TASK QUARTERLY 12 No 3, 255 265 SIMPLIFIED VOLUME OF FLUID METHOD (SVOF) FOR TWO-PHASE FLOWS MACIEJ MAREK, WOJCIECH ANISZEWSKI AND ANDRZEJ BOGUSŁAWSKI Institute of Thermal Machinery, Czestochowa University
More informationOn the thickness of discontinuities computed by THINC and RK schemes
The 9th Computational Fluid Dynamics Symposium B7- On the thickness of discontinuities computed by THINC and RK schemes Taku Nonomura, ISAS, JAXA, Sagamihara, Kanagawa, Japan, E-mail:nonomura@flab.isas.jaxa.jp
More informationOn the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes
Journal of Computational Physics 8, 87 09 (00) doi:0.006/jcph.00.79 On the Construction, Comparison, and Local Characteristic Decomposition for High-Order Central WENO Schemes Jianxian Qiu, and Chi-Wang
More informationModule 7 VIDEO CODING AND MOTION ESTIMATION
Module 7 VIDEO CODING AND MOTION ESTIMATION Lesson 20 Basic Building Blocks & Temporal Redundancy Instructional Objectives At the end of this lesson, the students should be able to: 1. Name at least five
More informationAn Explicit Implicit Scheme for Cut Cells in Embedded Boundary Meshes
An Explicit Implicit Scheme for Cut Cells in Embedded Boundary Meshes S. May and M. Berger Research Report No. 2015-44 December 2015 Seminar für Angewandte Mathematik Eidgenössische Technische Hochschule
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
More informationLarge scale CFD computations at CEA
1 Large scale CFD computations at CEA G. Meurant a,h.jourdren a and B. Meltz a a CEA/DIF, PB 12, 91680 Bruyères le Châtel, France This paper describes some large scale Computational Fluid Dynamics computations
More information3 No-Wait Job Shops with Variable Processing Times
3 No-Wait Job Shops with Variable Processing Times In this chapter we assume that, on top of the classical no-wait job shop setting, we are given a set of processing times for each operation. We may select
More informationCentral-Upwind Schemes on Triangular Grids for Hyperbolic Systems of Conservation Laws
Central-Upwind Schemes on Triangular Grids for Hyperbolic Systems of Conservation Laws Alexander Kurganov, 1 Guergana Petrova 2 1 Department of Mathematics, Tulane University, New Orleans, Louisiana 70118
More informationA limiting strategy for the back and forth error compensation and correction method for solving advection equations
A limiting strategy for the back and forth error compensation and correction method for solving advection equations Lili Hu, Yao Li, Yingjie Liu April 22, 2014 Abstract We further study the properties
More informationOverview of Traditional Surface Tracking Methods
Liquid Simulation With Mesh-Based Surface Tracking Overview of Traditional Surface Tracking Methods Matthias Müller Introduction Research lead of NVIDIA PhysX team PhysX GPU acc. Game physics engine www.nvidia.com\physx
More informationIMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING
Image Anal Stereol 2001;20:143-148 Original Research Paper IMAGE ANALYSIS DEDICATED TO POLYMER INJECTION MOLDING DAVID GARCIA 1, GUY COURBEBAISSE 2 AND MICHEL JOURLIN 3 1 European Polymer Institute (PEP),
More informationReconstructing Volume Tracking 1
JOURNAL OF COMPUTATIONAL PHYSICS 141, 112 152 (1998) ARTICLE NO. CP985906 Reconstructing Volume Tracking 1 William J. Rider and Douglas B. Kothe Applied Theoretical and Computational Physics Division,
More informationFully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon
Fully discrete Finite Element Approximations of Semilinear Parabolic Equations in a Nonconvex Polygon Tamal Pramanick 1,a) 1 Department of Mathematics, Indian Institute of Technology Guwahati, Guwahati
More informationDiscontinuous Fluctuation Distribution
Discontinuous Fluctuation Distribution Matthew Hubbard School of Computing, University of Leeds, Leeds, LS2 9JT, UK. Abstract This paper describes a new numerical scheme for the approximation of steady
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 informationHigh Order Weighted Essentially Non-Oscillatory Schemes for Convection. Dominated Problems. Chi-Wang Shu 1
High Order Weighted Essentially Non-Oscillatory Schemes for Convection Dominated Problems Chi-Wang Shu Division of Applied Mathematics, Brown University, Providence, Rhode Island 09 ABSTRACT High order
More informationChapter 6. Petrov-Galerkin Formulations for Advection Diffusion Equation
Chapter 6 Petrov-Galerkin Formulations for Advection Diffusion Equation In this chapter we ll demonstrate the difficulties that arise when GFEM is used for advection (convection) dominated problems. Several
More informationStudy on the Numerical Accuracy for the CFD
Study on the Numerical Accuracy for the CFD T.Yamanashi 1, H.Uchida, and M.Morita 1 Department of Mathematics, Master s Research Course of Faculty of Science, Tokyo University of Science,1-3 Kagurazaka,
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
AIAA SciTech Forum 8 January 08, Kissimmee, Florida 08 AIAA Aerospace Sciences Meeting 0.54/6.08-0067 A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley and Dr.
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 informationThe Immersed Interface Method
The Immersed Interface Method Numerical Solutions of PDEs Involving Interfaces and Irregular Domains Zhiiin Li Kazufumi Ito North Carolina State University Raleigh, North Carolina Society for Industrial
More informationVon Neumann Analysis for Higher Order Methods
1. Introduction Von Neumann Analysis for Higher Order Methods Von Neumann analysis is a widely used method to study how an initial wave is propagated with certain numerical schemes for a linear wave equation
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 informationThe Level Set Method. Lecture Notes, MIT J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations
The Level Set Method Lecture Notes, MIT 16.920J / 2.097J / 6.339J Numerical Methods for Partial Differential Equations Per-Olof Persson persson@mit.edu March 7, 2005 1 Evolving Curves and Surfaces Evolving
More informationCONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS
CONSERVATIVE AND NON-CONSERVATIVE METHODS BASED ON HERMITE WEIGHTED ESSENTIALLY-NON-OSCILLATORY RECONSTRUCTION FOR VLASOV EQUATIONS CHANG YANG AND FRANCIS FILBET Abstract. We develop weighted essentially
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationcomputational field which is always rectangular by construction.
I. INTRODUCTION The numerical solution of partial differential equations requires some discretization of the field into a collection of points or elemental volumes (cells). The differential equations are
More informationCecil Jones Academy Mathematics Fundamental Map
Fundamentals Year 7 (Core) Knowledge Unit 1 Unit 2 Solve problems involving prime numbers Use highest common factors to solve problems Use lowest common multiples to solve problems Explore powers and roots
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More informationOn the Number of Tilings of a Square by Rectangles
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange University of Tennessee Honors Thesis Projects University of Tennessee Honors Program 5-2012 On the Number of Tilings
More informationJ. Vira, M. Sofiev SILAM winter school, February 2013, FMI
Numerical aspects of the advection-diffusion equation J. Vira, M. Sofiev SILAM winter school, February 2013, FMI Outline Intro Some common requirements for numerical transport schemes Lagrangian approach
More informationVOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, VOL. 24, 67 69 (997) VOLUME-TRACKING METHODS FOR INTERFACIAL FLOW CALCULATIONS MURRAY RUDMAN Commonwealth Scientific and Industrial Research Organisation
More informationLecture 15: The subspace topology, Closed sets
Lecture 15: The subspace topology, Closed sets 1 The Subspace Topology Definition 1.1. Let (X, T) be a topological space with topology T. subset of X, the collection If Y is a T Y = {Y U U T} is a topology
More informationComputational Astrophysics 5 Higher-order and AMR schemes
Computational Astrophysics 5 Higher-order and AMR schemes Oscar Agertz Outline - The Godunov Method - Second-order scheme with MUSCL - Slope limiters and TVD schemes - Characteristics tracing and 2D slopes.
More informationScientific Visualization Example exam questions with commented answers
Scientific Visualization Example exam questions with commented answers The theoretical part of this course is evaluated by means of a multiple- choice exam. The questions cover the material mentioned during
More informationRealistic Animation of Fluids
1 Realistic Animation of Fluids Nick Foster and Dimitris Metaxas Presented by Alex Liberman April 19, 2005 2 Previous Work Used non physics-based methods (mostly in 2D) Hard to simulate effects that rely
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 informationWater. Notes. Free surface. Boundary conditions. This week: extend our 3D flow solver to full 3D water We need to add two things:
Notes Added a 2D cross-section viewer for assignment 6 Not great, but an alternative if the full 3d viewer isn t working for you Warning about the formulas in Fedkiw, Stam, and Jensen - maybe not right
More informationInvestigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow
Investigating The Stability of The Balance-force Continuum Surface Force Model of Surface Tension In Interfacial Flow Vinh The Nguyen University of Massachusetts Dartmouth Computational Science Training
More informationNon-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws
Non-Oscillatory Third Order Fluctuation Splitting Schemes for Steady Scalar Conservation Laws Matthew Hubbard School of Computing, University of Leeds, Leeds, LS 9JT, UK. Abstract This paper addresses
More informationMERBEIN P-10 COLLEGE MATHS SCOPE & SEQUENCE
MERBEIN P-10 COLLEGE MATHS SCOPE & SEQUENCE Year Number & Algebra Measurement & Geometry Statistics & Probability P Numbers to 20 Location Data collection Counting and comparing groups Length Ordinal numbers
More informationarxiv: v1 [math.co] 25 Sep 2015
A BASIS FOR SLICING BIRKHOFF POLYTOPES TREVOR GLYNN arxiv:1509.07597v1 [math.co] 25 Sep 2015 Abstract. We present a change of basis that may allow more efficient calculation of the volumes of Birkhoff
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 informationImage Sampling and Quantisation
Image Sampling and Quantisation Introduction to Signal and Image Processing Prof. Dr. Philippe Cattin MIAC, University of Basel 1 of 46 22.02.2016 09:17 Contents Contents 1 Motivation 2 Sampling Introduction
More informationA step towards the Bermond-Thomassen conjecture about disjoint cycles in digraphs
A step towards the Bermond-Thomassen conjecture about disjoint cycles in digraphs Nicolas Lichiardopol Attila Pór Jean-Sébastien Sereni Abstract In 1981, Bermond and Thomassen conjectured that every digraph
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 informationFrom Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
Preprint accepted in Journal of Computational Physics. https://doi.org/10.1016/j.jcp.2018.06.019 From Hyperbolic Diffusion Scheme to Gradient Method: Implicit Green-Gauss Gradients for Unstructured Grids
More informationImage Sampling & Quantisation
Image Sampling & Quantisation Biomedical Image Analysis Prof. Dr. Philippe Cattin MIAC, University of Basel Contents 1 Motivation 2 Sampling Introduction and Motivation Sampling Example Quantisation Example
More information3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition
3D Finite Difference Time-Domain Modeling of Acoustic Wave Propagation based on Domain Decomposition UMR Géosciences Azur CNRS-IRD-UNSA-OCA Villefranche-sur-mer Supervised by: Dr. Stéphane Operto Jade
More informationLevel Set Methods and Fast Marching Methods
Level Set Methods and Fast Marching Methods I.Lyulina Scientific Computing Group May, 2002 Overview Existing Techniques for Tracking Interfaces Basic Ideas of Level Set Method and Fast Marching Method
More informationThe Proportional Colouring Problem: Optimizing Buffers in Radio Mesh Networks
The Proportional Colouring Problem: Optimizing Buffers in Radio Mesh Networks Florian Huc, Claudia Linhares Sales, Hervé Rivano To cite this version: Florian Huc, Claudia Linhares Sales, Hervé Rivano.
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationA more efficient algorithm for perfect sorting by reversals
A more efficient algorithm for perfect sorting by reversals Sèverine Bérard 1,2, Cedric Chauve 3,4, and Christophe Paul 5 1 Département de Mathématiques et d Informatique Appliquée, INRA, Toulouse, France.
More informationScientific Computing: Interpolation
Scientific Computing: Interpolation Aleksandar Donev Courant Institute, NYU donev@courant.nyu.edu Course MATH-GA.243 or CSCI-GA.22, Fall 25 October 22nd, 25 A. Donev (Courant Institute) Lecture VIII /22/25
More informationweighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces.
weighted minimal surface model for surface reconstruction from scattered points, curves, and/or pieces of surfaces. joint work with (S. Osher, R. Fedkiw and M. Kang) Desired properties for surface reconstruction:
More informationPreferred directions for resolving the non-uniqueness of Delaunay triangulations
Preferred directions for resolving the non-uniqueness of Delaunay triangulations Christopher Dyken and Michael S. Floater Abstract: This note proposes a simple rule to determine a unique triangulation
More informationAurélien Thinat Stéphane Cordier 1, François Cany
SimHydro 2012:New trends in simulation - Hydroinformatics and 3D modeling, 12-14 September 2012, Nice Aurélien Thinat, Stéphane Cordier, François Cany Application of OpenFOAM to the study of wave loads
More informationAdaptive-Mesh-Refinement Pattern
Adaptive-Mesh-Refinement Pattern I. Problem Data-parallelism is exposed on a geometric mesh structure (either irregular or regular), where each point iteratively communicates with nearby neighboring points
More informationFlow Structures Extracted from Visualization Images: Vector Fields and Topology
Flow Structures Extracted from Visualization Images: Vector Fields and Topology Tianshu Liu Department of Mechanical & Aerospace Engineering Western Michigan University, Kalamazoo, MI 49008, USA We live
More informationHigh Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study
Commun. Comput. Phys. doi:.48/cicp.375.6a Vol., No. 4, pp. 835-869 October 6 High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study Liang
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 6:
file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_1.htm 1 of 1 6/20/2012 12:24 PM The Lecture deals with: ADI Method file:///d:/chitra/nptel_phase2/mechanical/cfd/lecture6/6_2.htm 1 of 2 6/20/2012
More informationStage 7 Checklists Have you reached this Standard?
Stage 7 Checklists Have you reached this Standard? Main Criteria for the whole year. J K L Use positive integer powers and associated real roots Apply the four operations with decimal numbers Write a quantity
More information