A Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory
|
|
- Henry Cummings
- 6 years ago
- Views:
Transcription
1 A Semi-Lagrangian Scheme for the Open Table Problem in Granular Matter Theory M. Falcone and S. Finzi Vita Abstract We introduce and analyze a new scheme for the approximation of the two-layer model proposed in [] for growing sandpiles over an open flat table. The method is based on a semi-lagrangian approximation of the nonlinear terms which allows to obtain more accurate results with respect to the standard finite difference approximation. We present several features of the scheme and give some hints on its implementation. Finally we show some tests where it is compared with a previously studied finite difference approach. Introduction Mathematical models for the dynamics of granular materials have recently been studied in several papers (see e.g. []), since a complete and realistic description of many phenomena in this field is far from being achieved. In this paper we restrict ourselves to the rather simple case of growing sandpiles over a flat bounded open table under the action of a given vertical source, neglecting wind effects and avalanches. For this phenomenon two main differential models have been proposed. The first one is the variational model of Prigozhin [], where the surface flow of sand is supposed to exist only at critical slope, that is the maximal admissible slope α for any stationary configuration of sand (here for simplicity it is assumed α = ). For that model there are a sufficiently developed theory and efficient numerical schemes which use duality arguments (see [2]). The second approach has been introduced by Hadeler and Kuttler [] as an extension of the known BCRE model [3] (see also [4]). In their model the pile is obtained summing two distinct Maurizio Falcone Dipartimento di Matematica, SAPIENZA - Università di Roma, falcone@mat.uniroma.it Stefano Finzi Vita Dipartimento di Matematica, SAPIENZA - Università di Roma, finzi@mat.uniroma.it 69
2 7 M. Falcone and S. Finzi Vita layers: the standing layer u and the rolling layer v. The two layers interact during the growth process, the exchange term between u and the v being described by the nonlinear term in the following system of partial differential equations: u t = ( Du )v in Ω (,T) v t = div(vdu) ( Du )v+ f in Ω (,T) () u(,) = v(,) = in Ω, u = on Γ (,T). In () f represents the source, the table is assumed initially empty, and no boundary condition is needed on v (the characteristic lines for the second equation are all directed towards the interior of Ω). Note that this model allows for the movement of rolling grains even at sub-critical slope. It also seems to be better suited for the description of fast processes or small details over the surface (see [2] for a detailed comparison between the two models). However, the mathematical theory for the two-layer model is still incomplete (even a general result of existence for its solutions is lacking). It is interesting to note that the two distinct models have essentially the same set of admissible stationary solutions (u, v) and that a characterization of such equilibria has been given, using the theory of viscosity solutions, in [5]: the maximal equilibrium u is given by the function d(x) measuring at any point x Ω its distance to the boundary, and an integral representation is proved for v along of the transport rays of sand. Anyway, as it was observed in [8], only in particular cases the asymptotic behavior is the same (and then u(x) = d(x)). In general, the effective stationary standing solution u for the model () is not explicitly known. It is also for such a reason that in [8] we proposed a finite difference scheme which preserves the properties of the model at the discrete level, and which is able to give a numerical characterization of stationary solutions for every general source support in a square open table. A disadvantage of this approach is that finite difference schemes on a regular structured grid restrict the numerical discrete flow essentially to the axis directions, breaking the homogeneous character of the real phenomenon (cones of sand are in general approximated by pyramids oriented along the grid axes). This numerical anisotropy is typical of finite difference (FD) schemes and trying to eliminate it is particularly important in view of the extension of schemes to the so-called partially-open table problem, that is when high vertical walls bound portions of the table boundary. In fact, in such a problem an infinite number of transport rays can meet at certains points on the boundary, yielding a discontinuous and even unbounded rolling layer v L (Ω) at the equilibrium, and strong numerical difficulties in its description (see [6]). This has motivated the development of a semi-lagrangian (SL) approach which will be discussed in this paper. By construction, the SL approach mimics the method of characteristics in order to follow in a more accurate way the real direction of sand flow. We introduce the scheme and its properties in Section 2, and in Section 3 we present first numerical tests in the regular case of the open table problem, showing a slight but clear improvement in the accuracy with respect to the standard FD approach. These results let us conjecture that in the partially-open table
3 A Semi-Lagrangian Scheme for the Open Table Problem 7 problem the SL approach could be of interest, despite of its time-consuming feature. Experiments in this direction are actually under consideration. A detailed presentation of SL methods for evolutive nonlinear partial differential equations of the first order can be found in [7]; see also [9] for a convergence result of a first order SL scheme to the viscosity solution. 2 The semi-lagrangian scheme and its properties Let us start from the first equation in () and observe that it can be equivalently written in one dimension as ( ) u t = v( u x ) = v + min au x, (2) a B(,) where B(,) denotes the closed unit ball of IR, that is the interval [,]. In this way the nonlinear term u x is written as the minimum over the directional derivatives of u and is actually computed comparing the values for a =, a = and a =. Note that in two dimensions we can use exactly the same formula to compute u. At the discrete level, if we introduce in Ω = (α,β) a uniform grid of points G = {x i : x i = α + i x,i =,..,N}, (2) can be written as u n+ i = u n i + tv n i ( u n (x i + a t) u n i + min a {,,} t ), (3) where, as usual, u n i denotes the approximate value for u(x i,n t). If t x and a, then the point z i (a) x i + a t is not a node of the grid G and the value of u n in it has to be computed, for example, by linear interpolation. For the equation in v, we have also second order terms. In fact, v t = v x u x + vu xx v( u x )+ f ; (4) let us use the standard second order central difference D 2 u n i to approximate u xx (x i ) whereas u x at any node is replaced by the maximal (in absolute value) finite difference in the left and right directions, Du n i. A simple trick is to replace the nonlinear term at the node x i by the previously computed difference (u n+ i u n i ). Finally, for v x we use the discrete directional derivative in the direction of u x. After a simplification, we finally obtain the following discretization of (4) v n+ i = v n (x i + tdu n i ) (un+ i u n i )+ t(vn i D2 u n i + f i). (5) Of course, also the term v n (x i + tdu n i ) requires a local reconstruction by interpolation. In order to complete the scheme we need to add initial and boundary conditions u i = v i = i ; u n = un N = n. (6)
4 72 M. Falcone and S. Finzi Vita The extension of this approach to IR 2 is straightforward.we assume for simplicity that Ω = (,) (,), that is the unit square table, and we introduce a uniform grid of points G = {x i, j = (i x, j x),i, j =,..,N}, where x = /N is again the space discretization step. Then the analogue of (3)-(5) and (6) is ( u n+ i, j = u n i, j + tv n i, j + ) min t a B(,) (un (x i, j + a t) u n i, j), v n+ i, j = v n (x i, j + tdu n i, j ) (un+ i, j u n i, j )+ t(vn i, j D2 u n i, j + f (7) i, j), u i, j = v i, j = i, j ; un i, j = n, if x i, j Ω, where u n i, j denotes the approximate value for u(x i, j,n t). In order to compute the value of u n (x i, j + a t) we use now a bilinear interpolation with respect to the four vertices of the cell containing x i, j + a t. The same local reconstruction is required for the term v n (x i, j + tdu n i, j ). The minimum term in (7) is approximated by comparing the values of u n on a finite set of directions (for example the eight directions θ k = kπ/4, k =,..,8, plus the origin, that is the node under consideration). D 2 u n i, j denotes the standard five-points second order difference which replaces u(x i, j ). 2. Properties of the SL scheme Here we present some features of the above scheme, showing that it preserves the physical properties of the continuous model. At the moment we are able to give a complete proof of that only in the one-dimensional case, but these properties are confirmed by all the 2D experiments. Further details and complete proofs will be presented in a forthcoming paper. Theorem. Let f in Ω and ( ) t x min 2,, (8) C f where C is a positive constant. Then, the sequences u n and v n defined in (3)-(5)-(6) satisfy the following properties for any n:. positivity and monotonicity in u: u n u n+ ; 2. positivity in v: v n ; 3. sub-critical slope in u: Du n. Proof. The structure of the proof is the same of that of Theorem 3. of [8]. The above properties are proved by induction. Here we only discuss some details which are characteristic of the semi-lagrangian scheme. In particular, for the nonlinear term in (3) we remark that min a {,,} u n (x i + a t) u n i t, n,i ; (9)
5 A Semi-Lagrangian Scheme for the Open Table Problem 73 moreover, since the intermediate values are computed in this case by linear interpolation and t x, we have for example that u n (x i + t) u n i t = un i+ un i x At the same time, if for example Du n i >, we can rewrite. () v n (x i + tdu n i ) = vn i + t x (vn i+ vn i )Dun i, () and the proof of positivity for v n follows the same arguments used for the FD scheme in [8], yielding the bound 2 t x. The stronger stability bound (8) has to be assumed in order to achieve the gradient constraint in every interval [x i,x i+ ]. An accurate discussion on admissible configurations allows to rewrite the different quotient in (9) always in terms of the slope of u n inside the interval under consideration. Then, considering the difference (u n i+ un i ) as defined by (3), one gets sub-criticality if ( t vn i ) x, which holds for the uniform boundedness of v n in terms of the source f. 3 Numerical results In this section we present some numerical tests in two dimensions to compare the SL method with the old FD method. For computation we used MATLAB 7. on a Processor Intel Pentium M74,.73GHz with 8Gb RAM. Test : constant source on all of Ω We assume f.5 in all the domain Ω = (,) (,). In this case the ridge set S (the set where d is singular) is given by the two diagonals of the square and, following [5], the stationary solutions are given by u(x) = d(x) and τ(x) v(x) = on S, v(x) = f(x+ tdu(x)) dt in Ω \ S, (2) where τ(x) denotes the distance of x from the set S along the transport ray through x. In Fig. we show the computed solutions u and v for x =.2, t =. (the stopping criterion for the equilibrium was: u n+ u n 6 ). We show in Table the results of the two schemes at the same final time when t/ x=.4, comparing the L and L norms of the error for u. It can be seen that the two methods are both approximately of first order. Nevertheless, the SL scheme is slightly more accurate in the L (+2%) and the L (+3%) norms respectively. Naturally, it is also much more expensive due to the minimum term computation at any iteration.
6 74 M. Falcone and S. Finzi Vita strato stazionario strato rotolante Fig. Test : Numerical results for u and v Table Errors for Test at CFL=.4 method x L error L error FD SL FD SL FD SL Test 2: constant source on a connected subdomain of Ω Assume now that the constant source is concentrated in a small square inside Ω (see Fig. 2), and it is zero elsewhere. An explicit formula for u is not known in this case (with some efforts v could be still computed by formula (2)), and Fig. 2 shows the stationary solution u as computed by the SL algorithm. In Fig. 3 the level lines found for u are compared with those found by the FD scheme u.6.8 Fig. 2 Test 2: The source support and the computed standing layer u
7 A Semi-Lagrangian Scheme for the Open Table Problem 75 curve di livello schema SL curve di livello schema DF Fig. 3 Test 2: Level curves for u with the SL and the FD schemes Test 3: constant source on a disconnected subdomain of Ω Assume now that the constant source is concentrated in two distinct small squares inside Ω, and it is zero elsewhere. Fig. 4 shows this source support together with the stationary solution u as computed by the SL algorithm. In Fig. 5 the level lines found for u by the two algorithms are compared u.6.8 Fig. 4 Test 3: The source support and the computed standing layer u Conclusions The above results seem to show the SL scheme, although enough expensive, can be more accurate than the corresponding FD scheme in the description of singularity regions for the solution. This can be seen by comparing the higher level sets of the two numerical solutions, those corresponding to the crests of the piles. This feature
8 76 M. Falcone and S. Finzi Vita curve di livello schema SL curve di livello schema DF Fig. 5 Test 3: Level curves for u with the SL and the FD schemes can be very important for more complicated problems and geometries, like for the wall and the silos problems, where hardest singularities can arise (see [6]). Acknowledgements We thank Alessia Pacella for the numerical experiments. References. Aranson, I.S., Tsimring, L.S.: Patterns and collective behavior in granular media: theoretical concepts. arxiv:cond-mat/5749 (25) 2. Barrett, J.W., Prigozhin, L.: Dual formulations in critical state problems. Interfaces and Free Boundaries 8, (26) 3. Bouchaud, J.-P., Cates, M.E., Ravi Prakash, J., Edwards, S.F.: A model for the dynamics of sandpile surfaces. J. Phys. I France 4, (994) 4. Boutreux, T., de Gennes, P.-G.: Surface flows of granular mixture, I. General principles and minimal model. J. Phys. I France 6, (996) 5. Cannarsa, P., Cardaliaguet, P.: Representation of equilibrium solutions to the table problem for growing sandpiles. J. Eur. Math. Soc. (JEMS) 6, (24) 6. Crasta, G., Finzi VIta, S.: An existence result for the sandpile problem on flat tables with walls. To appear; preprint available at (28) 7. Falcone, M., Ferretti, R.: Semi-Lagrangian schemes for Hamilton-Jacobi equations, discrete representation formulae and Godunov methods. J. Comp. Physics 75, (22) 8. Falcone, M., Finzi Vita, S.: A finite difference approximation of a two-layer system for growing sandpiles. SIAM J. Sci. Comput. 28, 2 32 (26) 9. Falcone M., Giorgi T.: An approximation scheme for evolutive Hamilton-Jacobi equations. In: W.M. McEneaney, G. Yin, Q. Zhang (eds.), Stochastic Analysis, Control, Optimization and Applications: A Volume in Honor of W.H. Fleming, Birkhäuser, (999). Hadeler, K.P., Kuttler, C.: Dynamical models for granular matter. Granular Matter 2, 9 8 (999). Prigozhin, L.: Variational model of sandpile growth. Euro. J. Appl. Math. 7, (996) 2. Prigozhin, L., Zaltzman, B.: Two continuous models for the dynamics of sandpiles surface. Physical Review E 63, 455 (2)
Differential Models for Sandpile Growth
Differential Models for Sandpile Growth Piermarco Cannarsa University of Rome Tor Vergata (Italy) http://www.mat.uniroma2.it LABORATOIRE JACQUES-LOUIS LIONS Universite Pierre et Marie Curie (Paris VI)
More informationSur le tas de Sable. Noureddine Igbida (LAMFA CNRS-UMR 6140, Université de Picardie Jules Sur le Verne, tas de SableAmiens)
Sur le tas de Sable Noureddine Igbida LAMFA CNRS-UMR 6140, Université de Picardie Jules Verne, 80000 Amiens Nancy, 4 Mars 2008 Nancy, 4 Mars 2008 1 / 26 Une source de sable qui tourne : résultat d une
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationFast Marching Methods for Front Propagation Lecture 1/3
Outline for Front Propagation Lecture 1/3 M. Falcone Dipartimento di Matematica SAPIENZA, Università di Roma Autumn School Introduction to Numerical Methods for Moving Boundaries ENSTA, Paris, November,
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationAn adaptive RBF-based Semi-Lagrangian scheme for HJB equations
An adaptive RBF-based Semi-Lagrangian scheme for HJB equations Roberto Ferretti Università degli Studi Roma Tre Dipartimento di Matematica e Fisica Numerical methods for Hamilton Jacobi equations in optimal
More informationNumerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle
Numerical Methods for Optimal Control Problems. Part I: Hamilton-Jacobi-Bellman Equations and Pontryagin Minimum Principle Ph.D. course in OPTIMAL CONTROL Emiliano Cristiani (IAC CNR) e.cristiani@iac.cnr.it
More informationConvergence of a first order scheme for a non local eikonal equation
Convergence of a first order scheme for a non local eikonal equation O. Alvarez, E. Carlini, R. Monneau, E. Rouy March 22, 2005 Abstract We prove the convergence of a first order finite difference scheme
More information1.The anisotropic plate model
Pré-Publicações do Departamento de Matemática Universidade de Coimbra Preprint Number 6 OPTIMAL CONTROL OF PIEZOELECTRIC ANISOTROPIC PLATES ISABEL NARRA FIGUEIREDO AND GEORG STADLER ABSTRACT: This paper
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationThe HJB-POD approach for infinite dimensional control problems
The HJB-POD approach for infinite dimensional control problems M. Falcone works in collaboration with A. Alla, D. Kalise and S. Volkwein Università di Roma La Sapienza OCERTO Workshop Cortona, June 22,
More informationSuboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids
Suboptimal feedback control of PDEs by solving Hamilton-Jacobi Bellman equations on sparse grids Jochen Garcke joint work with Axel Kröner, INRIA Saclay and CMAP, Ecole Polytechnique Ilja Kalmykov, Universität
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationFinite Volume Schemes: an introduction
Finite Volume Schemes: an introduction First lecture Annamaria Mazzia Dipartimento di Metodi e Modelli Matematici per le Scienze Applicate Università di Padova mazzia@dmsa.unipd.it Scuola di dottorato
More informationMEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 40, Number 7, July 0, Pages 453 463 S 000-9939(0)8-X Article electronically published on November, 0 MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationminimize x subject to (x 2)(x 4) u,
Math 6366/6367: Optimization and Variational Methods Sample Preliminary Exam Questions 1. Suppose that f : [, L] R is a C 2 -function with f () on (, L) and that you have explicit formulae for
More informationRiemann integral and volume are generalized to unbounded functions and sets. is an admissible set, and its volume is a Riemann integral, 1l E,
Tel Aviv University, 26 Analysis-III 9 9 Improper integral 9a Introduction....................... 9 9b Positive integrands................... 9c Special functions gamma and beta......... 4 9d Change of
More informationOn Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities
On Nonlinear Dirichlet Neumann Algorithms for Jumping Nonlinearities Heiko Berninger, Ralf Kornhuber, and Oliver Sander FU Berlin, FB Mathematik und Informatik (http://www.math.fu-berlin.de/rd/we-02/numerik/)
More informationA Semi-Lagrangian discretization of non linear fokker Planck equations
A Semi-Lagrangian discretization of non linear fokker Planck equations E. Carlini Università Sapienza di Roma joint works with F.J. Silva RICAM, Linz November 2-25, 206 / 3 Outline A Semi-Lagrangian scheme
More informationSYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS. Massimo Grosi Filomena Pacella S. L. Yadava. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 21, 2003, 211 226 SYMMETRY RESULTS FOR PERTURBED PROBLEMS AND RELATED QUESTIONS Massimo Grosi Filomena Pacella S.
More informationScalable Non-Linear Compact Schemes
Scalable Non-Linear Compact Schemes Debojyoti Ghosh Emil M. Constantinescu Jed Brown Mathematics Computer Science Argonne National Laboratory International Conference on Spectral and High Order Methods
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationMIT (Spring 2014)
18.311 MIT (Spring 014) Rodolfo R. Rosales May 6, 014. Problem Set # 08. Due: Last day of lectures. IMPORTANT: Turn in the regular and the special problems stapled in two SEPARATE packages. Print your
More informationNumerical approximation for optimal control problems via MPC and HJB. Giulia Fabrini
Numerical approximation for optimal control problems via MPC and HJB Giulia Fabrini Konstanz Women In Mathematics 15 May, 2018 G. Fabrini (University of Konstanz) Numerical approximation for OCP 1 / 33
More informationA Ginzburg-Landau approach to dislocations. Marcello Ponsiglione Sapienza Università di Roma
Marcello Ponsiglione Sapienza Università di Roma Description of a dislocation line A deformed crystal C can be described by The displacement function u : C R 3. The strain function β = u. A dislocation
More informationMean field games and related models
Mean field games and related models Fabio Camilli SBAI-Dipartimento di Scienze di Base e Applicate per l Ingegneria Facoltà di Ingegneria Civile ed Industriale Email: Camilli@dmmm.uniroma1.it Web page:
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationPrimal-Dual Interior-Point Methods for Linear Programming based on Newton s Method
Primal-Dual Interior-Point Methods for Linear Programming based on Newton s Method Robert M. Freund March, 2004 2004 Massachusetts Institute of Technology. The Problem The logarithmic barrier approach
More informationOn self-organised criticality in one dimension
On self-organised criticality in one dimension Kim Christensen Imperial College ondon Department of Physics Prince Consort Road SW7 2BW ondon United Kingdom Abstract In critical phenomena, many of the
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationFinite difference method for solving Advection-Diffusion Problem in 1D
Finite difference method for solving Advection-Diffusion Problem in 1D Author : Osei K. Tweneboah MATH 5370: Final Project Outline 1 Advection-Diffusion Problem Stationary Advection-Diffusion Problem in
More informationFinite Difference Methods for Boundary Value Problems
Finite Difference Methods for Boundary Value Problems October 2, 2013 () Finite Differences October 2, 2013 1 / 52 Goals Learn steps to approximate BVPs using the Finite Difference Method Start with two-point
More informationA convergence result for an Outer Approximation Scheme
A convergence result for an Outer Approximation Scheme R. S. Burachik Engenharia de Sistemas e Computação, COPPE-UFRJ, CP 68511, Rio de Janeiro, RJ, CEP 21941-972, Brazil regi@cos.ufrj.br J. O. Lopes Departamento
More informationA variational approach to the macroscopic. Electrodynamics of hard superconductors
A variational approach to the macroscopic electrodynamics of hard superconductors Dept. of Mathematics, Univ. of Rome La Sapienza Torino, July 4, 2006, joint meeting U.M.I. S.M.F. Joint work with Annalisa
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationNonlinear aspects of Calderón-Zygmund theory
Ancona, June 7 2011 Overture: The standard CZ theory Consider the model case u = f in R n Overture: The standard CZ theory Consider the model case u = f in R n Then f L q implies D 2 u L q 1 < q < with
More informationBlock-Structured Adaptive Mesh Refinement
Block-Structured Adaptive Mesh Refinement Lecture 2 Incompressible Navier-Stokes Equations Fractional Step Scheme 1-D AMR for classical PDE s hyperbolic elliptic parabolic Accuracy considerations Bell
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationOptimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations
Martino Bardi Italo Capuzzo-Dolcetta Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations Birkhauser Boston Basel Berlin Contents Preface Basic notations xi xv Chapter I. Outline
More informationOn an initial-value problem for second order partial differential equations with self-reference
Note di Matematica ISSN 113-536, e-issn 159-93 Note Mat. 35 15 no. 1, 75 93. doi:1.185/i15993v35n1p75 On an initial-value problem for second order partial differential equations with self-reference Nguyen
More informationMultiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions
Multiple integrals: Sufficient conditions for a local minimum, Jacobi and Weierstrass-type conditions March 6, 2013 Contents 1 Wea second variation 2 1.1 Formulas for variation........................
More informationGeneralized Hypothesis Testing and Maximizing the Success Probability in Financial Markets
Generalized Hypothesis Testing and Maximizing the Success Probability in Financial Markets Tim Leung 1, Qingshuo Song 2, and Jie Yang 3 1 Columbia University, New York, USA; leung@ieor.columbia.edu 2 City
More informationUnbounded Convex Semialgebraic Sets as Spectrahedral Shadows
Unbounded Convex Semialgebraic Sets as Spectrahedral Shadows Shaowei Lin 9 Dec 2010 Abstract Recently, Helton and Nie [3] showed that a compact convex semialgebraic set S is a spectrahedral shadow if the
More informationHomogenization and error estimates of free boundary velocities in periodic media
Homogenization and error estimates of free boundary velocities in periodic media Inwon C. Kim October 7, 2011 Abstract In this note I describe a recent result ([14]-[15]) on homogenization and error estimates
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationMean Field Games on networks
Mean Field Games on networks Claudio Marchi Università di Padova joint works with: S. Cacace (Rome) and F. Camilli (Rome) C. Marchi (Univ. of Padova) Mean Field Games on networks Roma, June 14 th, 2017
More informationChemical Equilibrium: A Convex Optimization Problem
Chemical Equilibrium: A Convex Optimization Problem Linyi Gao June 4, 2014 1 Introduction The equilibrium composition of a mixture of reacting molecules is essential to many physical and chemical systems,
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationFast convergent finite difference solvers for the elliptic Monge-Ampère equation
Fast convergent finite difference solvers for the elliptic Monge-Ampère equation Adam Oberman Simon Fraser University BIRS February 17, 211 Joint work [O.] 28. Convergent scheme in two dim. Explicit solver.
More informationFoliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary
Foliations of hyperbolic space by constant mean curvature surfaces sharing ideal boundary David Chopp and John A. Velling December 1, 2003 Abstract Let γ be a Jordan curve in S 2, considered as the ideal
More informationOn intermediate value theorem in ordered Banach spaces for noncompact and discontinuous mappings
Int. J. Nonlinear Anal. Appl. 7 (2016) No. 1, 295-300 ISSN: 2008-6822 (electronic) http://dx.doi.org/10.22075/ijnaa.2015.341 On intermediate value theorem in ordered Banach spaces for noncompact and discontinuous
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationCHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS. W. Erwin Diewert January 31, 2008.
1 ECONOMICS 594: LECTURE NOTES CHAPTER 2: CONVEX SETS AND CONCAVE FUNCTIONS W. Erwin Diewert January 31, 2008. 1. Introduction Many economic problems have the following structure: (i) a linear function
More information3.10 Lagrangian relaxation
3.10 Lagrangian relaxation Consider a generic ILP problem min {c t x : Ax b, Dx d, x Z n } with integer coefficients. Suppose Dx d are the complicating constraints. Often the linear relaxation and the
More informationSung-Ik Sohn and Jun Yong Shin
Commun. Korean Math. Soc. 17 (2002), No. 1, pp. 103 120 A SECOND ORDER UPWIND METHOD FOR LINEAR HYPERBOLIC SYSTEMS Sung-Ik Sohn and Jun Yong Shin Abstract. A second order upwind method for linear hyperbolic
More information1. Nonlinear Equations. This lecture note excerpted parts from Michael Heath and Max Gunzburger. f(x) = 0
Numerical Analysis 1 1. Nonlinear Equations This lecture note excerpted parts from Michael Heath and Max Gunzburger. Given function f, we seek value x for which where f : D R n R n is nonlinear. f(x) =
More informationConvergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates
Convergence of an immersed interface upwind scheme for linear advection equations with piecewise constant coefficients I: L 1 -error estimates Xin Wen and Shi Jin Abstract We study the L 1 -error estimates
More information256 Summary. D n f(x j ) = f j+n f j n 2n x. j n=1. α m n = 2( 1) n (m!) 2 (m n)!(m + n)!. PPW = 2π k x 2 N + 1. i=0?d i,j. N/2} N + 1-dim.
56 Summary High order FD Finite-order finite differences: Points per Wavelength: Number of passes: D n f(x j ) = f j+n f j n n x df xj = m α m dx n D n f j j n= α m n = ( ) n (m!) (m n)!(m + n)!. PPW =
More informationRational Chebyshev pseudospectral method for long-short wave equations
Journal of Physics: Conference Series PAPER OPE ACCESS Rational Chebyshev pseudospectral method for long-short wave equations To cite this article: Zeting Liu and Shujuan Lv 07 J. Phys.: Conf. Ser. 84
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationReceived 15 September 2006; accepted 30 November 2006 Available online 29 December 2006
Ann. I. H. Poincaré AN 5 008) 163 171 www.elsevier.com/locate/anihpc An explicit solution to a system of implicit differential equations Une solution explicite d un système implicite d équations différentielles
More informationAdvection Diffusion Problems with Pure Advection Approximation in Subregions
Advection Diffusion Problems with Pure Advection Approximation in Subregions M. J. Gander, L. Halpern 2, C. Japhet 2, and V. Martin 3 Université de Genève, 2-4 rue du Lièvre, CP 64, CH-2 Genève. Switzerland.
More informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More information26.2. Cauchy-Riemann Equations and Conformal Mapping. Introduction. Prerequisites. Learning Outcomes
Cauchy-Riemann Equations and Conformal Mapping 26.2 Introduction In this Section we consider two important features of complex functions. The Cauchy-Riemann equations provide a necessary and sufficient
More informationA numerical study of SSP time integration methods for hyperbolic conservation laws
MATHEMATICAL COMMUNICATIONS 613 Math. Commun., Vol. 15, No., pp. 613-633 (010) A numerical study of SSP time integration methods for hyperbolic conservation laws Nelida Črnjarić Žic1,, Bojan Crnković 1
More informationNoncooperative continuous-time Markov games
Morfismos, Vol. 9, No. 1, 2005, pp. 39 54 Noncooperative continuous-time Markov games Héctor Jasso-Fuentes Abstract This work concerns noncooperative continuous-time Markov games with Polish state and
More informationMath 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework
Math 7824 Spring 2010 Numerical solution of partial differential equations Classroom notes and homework Jan Mandel University of Colorado Denver May 12, 2010 1/20/09: Sec. 1.1, 1.2. Hw 1 due 1/27: problems
More informationLongitudinal waves of a finite amplitude in nonlinear elastic media
Longitudinal waves of a finite amplitude in nonlinear elastic media Ivan A. Molotkov 1), Ivan Pšenčík 2) 1) IZMIRAN, Troitsk, Moscow Region, 142190, Russia. E-mail molotkov@izmiran.rssi.ru 2) Geophysical
More informationA Computational Approach to Study a Logistic Equation
Communications in MathematicalAnalysis Volume 1, Number 2, pp. 75 84, 2006 ISSN 0973-3841 2006 Research India Publications A Computational Approach to Study a Logistic Equation G. A. Afrouzi and S. Khademloo
More informationOn the approximation of the principal eigenvalue for a class of nonlinear elliptic operators
On the approximation of the principal eigenvalue for a class of nonlinear elliptic operators Fabio Camilli ("Sapienza" Università di Roma) joint work with I.Birindelli ("Sapienza") I.Capuzzo Dolcetta ("Sapienza")
More informationSHARP BOUNDARY TRACE INEQUALITIES. 1. Introduction
SHARP BOUNDARY TRACE INEQUALITIES GILES AUCHMUTY Abstract. This paper describes sharp inequalities for the trace of Sobolev functions on the boundary of a bounded region R N. The inequalities bound (semi-)norms
More informationFast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma
Fast Algorithms for SDPs derived from the Kalman-Yakubovich-Popov Lemma Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 8 September 2003 European Union RTN Summer School on Multi-Agent
More informationNumerical Simulations on Two Nonlinear Biharmonic Evolution Equations
Numerical Simulations on Two Nonlinear Biharmonic Evolution Equations Ming-Jun Lai, Chun Liu, and Paul Wenston Abstract We numerically simulate the following two nonlinear evolution equations with a fourth
More informationConstruction of Lyapunov functions by validated computation
Construction of Lyapunov functions by validated computation Nobito Yamamoto 1, Kaname Matsue 2, and Tomohiro Hiwaki 1 1 The University of Electro-Communications, Tokyo, Japan yamamoto@im.uec.ac.jp 2 The
More informationOptimal shape and position of the support for the internal exact control of a string
Optimal shape and position of the support for the internal exact control of a string Francisco Periago Abstract In this paper, we consider the problem of optimizing the shape and position of the support
More informationNumerical explorations of a forward-backward diffusion equation
Numerical explorations of a forward-backward diffusion equation Newton Institute KIT Programme Pauline Lafitte 1 C. Mascia 2 1 SIMPAF - INRIA & U. Lille 1, France 2 Univ. La Sapienza, Roma, Italy September
More informationCONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS AND FREE BOUNDARY PROBLEMS
SIAM J. NUMER. ANAL. Vol. 44, No. 2, pp. 879 895 c 26 Society for Industrial and Applied Mathematics CONVERGENT DIFFERENCE SCHEMES FOR DEGENERATE ELLIPTIC AND PARABOLIC EQUATIONS: HAMILTON JACOBI EQUATIONS
More informationA VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION 1. INTRODUCTION
A VARIATIONAL METHOD FOR THE ANALYSIS OF A MONOTONE SCHEME FOR THE MONGE-AMPÈRE EQUATION GERARD AWANOU AND LEOPOLD MATAMBA MESSI ABSTRACT. We give a proof of existence of a solution to the discrete problem
More informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationStochastic Homogenization for Reaction-Diffusion Equations
Stochastic Homogenization for Reaction-Diffusion Equations Jessica Lin McGill University Joint Work with Andrej Zlatoš June 18, 2018 Motivation: Forest Fires ç ç ç ç ç ç ç ç ç ç Motivation: Forest Fires
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationWaves in a Shock Tube
Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid
More informationCubic Splines. Antony Jameson. Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305
Cubic Splines Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, California, 94305 1 References on splines 1. J. H. Ahlberg, E. N. Nilson, J. H. Walsh. Theory of
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
More informationS chauder Theory. x 2. = log( x 1 + x 2 ) + 1 ( x 1 + x 2 ) 2. ( 5) x 1 + x 2 x 1 + x 2. 2 = 2 x 1. x 1 x 2. 1 x 1.
Sep. 1 9 Intuitively, the solution u to the Poisson equation S chauder Theory u = f 1 should have better regularity than the right hand side f. In particular one expects u to be twice more differentiable
More informationOn the construction of discrete gradients
On the construction of discrete gradients Elizabeth L. Mansfield School of Mathematics, Statistics and Actuarial Science, University of Kent, Canterbury, CT2 7NF, U.K. G. Reinout W. Quispel Department
More informationWeighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods
Weighted Essentially Non-Oscillatory limiters for Runge-Kutta Discontinuous Galerkin Methods Jianxian Qiu School of Mathematical Science Xiamen University jxqiu@xmu.edu.cn http://ccam.xmu.edu.cn/teacher/jxqiu
More informationA finite element level set method for anisotropic mean curvature flow with space dependent weight
A finite element level set method for anisotropic mean curvature flow with space dependent weight Klaus Deckelnick and Gerhard Dziuk Centre for Mathematical Analysis and Its Applications, School of Mathematical
More informationA NONLINEAR DIFFERENTIAL EQUATION INVOLVING REFLECTION OF THE ARGUMENT
ARCHIVUM MATHEMATICUM (BRNO) Tomus 40 (2004), 63 68 A NONLINEAR DIFFERENTIAL EQUATION INVOLVING REFLECTION OF THE ARGUMENT T. F. MA, E. S. MIRANDA AND M. B. DE SOUZA CORTES Abstract. We study the nonlinear
More informationEXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN EQUATION WITH DEGENERATE MOBILITY
Electronic Journal of Differential Equations, Vol. 216 216), No. 329, pp. 1 22. ISSN: 172-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS TO THE CAHN-HILLIARD/ALLEN-CAHN
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationPropagation d interfaces avec termes non locaux
Propagation d interfaces avec termes non locaux P. Cardaliaguet Univ. Brest Janvier 2008 Joint works with G. Barles (Tours), O. Alvarez (Rouen), O. Ley (Tours), R. Monneau (CERMICS), A. Monteillet (Brest).
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationDefinition 3 (Continuity). A function f is continuous at c if lim x c f(x) = f(c).
Functions of Several Variables A function of several variables is just what it sounds like. It may be viewed in at least three different ways. We will use a function of two variables as an example. z =
More informationA Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations
2005 American Control Conference June 8-10, 2005. Portland, OR, USA WeB10.3 A Remark on IVP and TVP Non-Smooth Viscosity Solutions to Hamilton-Jacobi Equations Arik Melikyan, Andrei Akhmetzhanov and Naira
More information