Uniform exponential decay of the free energy for Voronoi finite volume discretized reaction-diffusion systems
|
|
- Meryl Cross
- 6 years ago
- Views:
Transcription
1 Weierstrass Institute for Applied Analysis and Stochastics Special Session 36 Reaction Diffusion Systems 8th AIMS Conference on Dynamical Systems, Differential Equations & Applications Dresden University of Technology, May 25 28, 2010 Uniform exponential decay of the free energy for Voronoi finite volume discretized reaction-diffusion systems Annegret Glitzky Mohrenstrasse Berlin Germany Tel th AIMS Conference
2 Contents Model equations Energy estimates for the non-discretized problem Discretization scheme Results for finite-dimensional problems Remarks and generalizations Formulation in space of piecewise constant functions Uniform exponential decay of the free energy for classes of discretizations References Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 2 (20)
3 Model equations X i species, i = 1,...,m ū i reference densities v i chemical potentials u i = ū i e v i densities reversible reactions α 1 X α m X m β 1 X β m X m, (α,β) R Z m + Z m + reaction rates k αβ (e v α e v β ), v = (v 1,...,v m ) generation rate of species X i R i := k αβ (e v α e v β )(β i α i ) (α,β) R flux density j i = u i µ i v i, i = 1,...,m, continuity equations u i t + j i = R i in R + Ω, ν j i = 0 on R + Γ, u i (0) = U i in Ω, i = 1,...,m. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 3 (20)
4 Weak formulation v = (v 1,...,v m ) V = H 1 (Ω;R m ), u = (u 1,...,u m ) V Operators A, E : V L (Ω,R m ) V { m Av, v V := Ω ū i e v i µ i v i v i + i=1 m Ev, v V := ū i e vi v i dx, v V Ω i=1 (α, β) R k αβ (e α v e β v) } (α β) v dx Problem (P) u (t) + Av(t) = 0, u(t) = Ev(t) f.a.a. t R +, u(0) = U, u H 1 loc (R +;V ), v L 2 loc (R +;V ) L loc (R +;L (Ω) m ). } (P) Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 4 (20)
5 Energy functionals Dissipation rate m D(v) = ū i e v i µ i v i v i dx Ω i=1 + k αβ (e v α e v β )(α β) vdx 0 v V L (Ω,R m ) Ω (α,β) R Free energy F(u) = Ω m i=1 ( ui (ln u i ū i 1) + ū i ) dx Stoichiometric subspace S = span { α β : (α,β) R } Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 5 (20)
6 Results for the continuous problem (G./Hünlich 97, G./Gärtner 09) Invariants: (u, v) solution to (P) = { } u(t) U U := u V : ( u 1,1,..., u m,1 ) S t > 0. Thermodynamic equilibrium: There exists a unique solution (u,v ) to Av = 0, u = Ev, u U U. (S) Moreover, v = 0 and v S. Monotonous decay of the free energy: Let (u,v) be a solution to (P). Then F(u(t 2 )) F(u(t 1 )) F(U) for t 2 t 1 0. Exponential decay of the free energy: Let (u,v) be a solution to (P) and let (u,v ) be the thermodynamic equilibrium. Then there exists a constant λ > 0 such that F(u(t)) F(u ) e λt (F(U) F(u )) t 0. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 6 (20)
7 Voronoi finite volume discretization Ω R N, N 3, open, bounded, polyhedral domain M = (P,T,E ) set P of grid points in Ω, P is assumed to be boundary conforming Delaunay family T of control volumes, family E of faces of the Voronoi boxes E int E, E K E, M = #P xk K mσ L σ = K L K L x L dσ xl x K For x K P the Voronoi box K is given by K = {x Ω : x x K < x x L x L P, x L x K }, K T. m σ = (n 1)-dimensional measure of Voronoi face σ, d σ = x K x L for σ = K L E int Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 7 (20)
8 Discretization scheme u K i Mass of species X i in the box K u = (u K i ) K T,i=1,...,m R Mm v ik chemical potential of the species X i in x K v = (v ik ) K T,i=1,...,m R Mm Voronoi finite volume - implicit Euler - discretization u K i (t n) u K i (t n 1) t n t n 1 Yi σ σ E K E int Z σ i (t n)(v il (t n ) v ik (t n )) m σ d σ u K i (t n) = ū ik e v ik(t n ) mes(k), n 1, u K i (0) = U i dx, K T, i = 1,...,m, K = R K i (t n), (P M ) Zi σ = 2 1 (ev ik + e v il ), σ = K L, Yi σ suitable average of ū i µ i corresponding to face σ ū ik, k αβk average of ū i, k αβ on the box K, R K i (t n) = k αβk (e α v K(t n ) e β v K(t n ) )(β i α i )mes(k) (α,β) R Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 8 (20)
9 Results for each finite-dimensional problem (P M ) (Glitzky 08) Invariants: ( u, v) (R Mm ) 2 solution to the discretized problem (P M ) = u(t n ) U { u R Mm : ( u K 1,..., u K ) } m S K T K T n N Stationary solutions: There is exactly one stationary solution ( u, v ) to (P M ) fulfilling the above invariance property. ( u, v ) is a discrete thermodynamic equilibrium. The scheme is dissipative and the discrete free energy m ˆF( u) = i=1 K T ( ) u K i v ik u K i + ū ik mes(k) decays monotonously and exponentially to its equilibrium value. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 9 (20)
10 Remarks Results are also valid for (Gärtner/G. 09, G. 08) anisotropic diffusion S i positive definite, symmetric flux densities j i = u i S i v i Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 10 (20)
11 Remarks Results are also valid for (Gärtner/G. 09, G. 08) anisotropic diffusion S i positive definite, symmetric flux densities j i = u i S i v i more general statistical relations u i = ū i g i (v i v i ) flux densities j i = ū i g i (v i v i )S i v i Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 10 (20)
12 Remarks Results are also valid for (Gärtner/G. 09, G. 08) anisotropic diffusion S i positive definite, symmetric flux densities j i = u i S i v i more general statistical relations u i = ū i g i (v i v i ) flux densities j i = ū i g i (v i v i )S i v i charged species X i with charge numbers q i in addition Poisson equation for the electrostatic potential v 0 m (S 0 v 0 ) = f + q i u i on R + Ω + BCs i=1 flux densities j i = ū i g i (v i v i )S i (v i + q i v 0 ) Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 10 (20)
13 Remarks concerning anisotropies Use for each species X i, i = 1,...,m, their own anisotropic Voronoi boxes K (i) = {x Ω : (x x K ) T S 1 i (x x K ) (x x L ) T S 1 i (x x L ) x L P} fluxes through faces σ (i) additionally contain the factor S i ν σ (i) reaction terms have to be evaluated on all different intersection sets K (1) k 1 K (2) k 2 K (m) k m, k i {1,...,M} Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 11 (20)
14 Formulation of (P M ) by piecewise constant functions Aim: Uniform exponential decay of the free energy for a class of Voronoi finite volume discretizations M of Ω with the properties diam(σ) κ 1 d σ σ E int, r K,out κ 2 r K,in K T (Q) X(M ) = set of functions from Ω to R, being constant on each K T. v K = value of v X(M ) on K. Discrete H 1 -seminorm for v X(M ) v 2 1,M = D σ v 2 m σ, d σ E int σ where D σ v = v K v L for σ = K L. use functions v i u i ū i k αβ a i = e v i /2 X(M ) v ik u ik ū ik k αβk e v ik/2 on K where u ik := uk i mes(k) Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 12 (20)
15 Formulation of (P M ) by piecewise constant functions E M : X(M ) X(M ), E M (v) = (ū i e v i )i=1,...,m, u = E M (v) Free energy: F M : X(M ) R, m F M (u) = (u i v i u i + ū i )dx F(u) Ω i=1 Dissipation rate: m D M (v) = Yi σ i=1 σ E int m c i=1 Zi σ D σ v i 2 m σ + d σ a i 2 1,M + c Ω (α,β) R [ m i=1 [ k αβ e α v e β v] (α β) v dx Ω (α,β) R a α i i m i=1 a β ] i 2 i dx 0, ai = e v i /2. Invariants: u(t n ) U U n N For the thermodynamic equilibrium (u,v ) to (P M ) (resulting from ( u, v )) we get u i = ūi ū i u i, i = 1,...,m, v = v S. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 13 (20)
16 Formulation of (P M ) by piecewise constant functions (G. Prep. 1443) Theorem A. [Uniform estimate of the free energy by the dissipation rate] Let a class of Voronoi finite volume discretizations fulfill (Q), moreover let (u,v ) be the thermodynamic equilibrium to (P M ). Then, for all ρ > 0 there exists a constant c ρ > 0 such that F M (E M v) F M (u ) c ρ D M (v) { } for all v N ρm := v X(M ) : F M (E M v) F M (u ) ρ, E M v U U uniformly for all Voronoi finite volume discretizations M. Indirect proof: Suppose there is {M l } with size(m l ) 0 and corresponding v l N ρml such that construct a contradiction. F Ml (E Ml v l ) F Ml (u l ) = C ld Ml (v l ) with C l, Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 14 (20)
17 Formulation of (P M ) by piecewise constant functions Essential tool in this proof: Theorem B. [Discrete Sobolev-Poincaré inequality] (G./Griepentrog 10) u m Ω (u) L q (Ω) c q u 1,M u X(M ) for q [1, ) if N = 2, q [1, N 2 2N ) if N > 2 uniformly for all M with the property (Q). D Ml (v l ) 0 = a Ml i 1,Ml 0 = a Ml i m Ω (a Ml i) 0 in L q (Ω), m Ω (a Ml i) e v i /2 guarantees the limit passage l in parts of the dissipation rate resulting from the reaction terms. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 15 (20)
18 Uniform exponential decay of the free energy for solutions to (P M ) Theorem C. [Uniform exponential decay of the free energy] We consider Voronoi finite volume discretizations M fulfilling the property (Q). Let (u,v) be a solution to (P M ) and let (u,v ) be the thermodynamic equilibrium to (P M ). Then F M (u(t n )) F M (u ) e λ t n ( FM (U) F M (u ) ) n 1 with a unified constant λ > 0. Idea of the proof: Discrete Brézis formula: For n 2 > n 1 0 and λ 0 we obtain ) ( ) e λt n 2 (F M (u(t n2 )) F M (u ) e λt n 1 F M (u(t n1 )) F M (u ) n 2 { e λt r 1 (t r t r 1 ) e λ h ( ) } λ F M (u(t r )) F M (u ) D M (v(t r )). r=n 1 +1 Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 16 (20)
19 Idea of the proof For n 2 > n 1 0 and λ 0 we obtain ) ( ) e λt n 2 (F M (u(t n2 )) F M (u ) e λt n 1 F M (u(t n1 )) F M (u ) n 2 { e λt r 1 (t r t r 1 ) e λ h ( ) } λ F M (u(t r )) F M (u ) D M (v(t r )). r=n 1 +1 λ = 0: F M (u(t n2 )) + n 2 (t r t r 1 )D M (v(t r )) F M (u(t n1 )) F M (U) r=n 1 +1 F M decays monotonously, F M (U) F M (u ) ρ for all M with (Q), = v(t r ) N ρm for r 1 E M v(t r ) U U apply Theorem A, choose λ > 0 such that λ e λ h c ρ < 1 = result follows uniformly for all M with (Q). Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 17 (20)
20 References R. Eymard, T. Gallouët, and R. Herbin, The finite volume method, Handbook of numerical Analysis VII (2000), A. Glitzky, Exponential decay of the free energy for discretized electro-reaction-diffusion systems, Nonlinearity 21 (2008), A. Glitzky, Uniform exponential decay of the free energy for Voronoi finite volume discretized reaction- diffusion systems, WIAS-Preprint 1443 (2009). A. Glitzky and K. Gärtner, Energy estimates for continuous and discretized electroreaction-diffusion systems, Nonlin. Anal. 70 (2009), A. Glitzky and J. A. Giepentrog, Discrete Sobolev-Poincaré inequalities for Voronoi finite volume approximations, SINUM 48 (2010), A. Glitzky and R. Hünlich, Energetic estimates and asymptotics for reaction diffusion systems, Z. Angew. Math. Mech. 77 (1997), Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 18 (20)
21 Assumptions Ω R N, N = 2,3, open, bounded polyhedral domain; µ i, u i, U i L+(Ω), µ i, u i δ > 0, q i Z, i = 1,...,m; R Z m + Zm + finite subset, k αβ L+(Ω), k αβ δ > 0 for all (α,β) R. { } If N = 3 then max (α,β) R max m i=1 α i, m i=1 β i 3; (< 3 for uniform exponential decay of the free energy) m i=1 U i,λ i > 0 λ S + \ {0} (Slater condition) {b R m + : bα = b β (α,β) R, u U U, u i = ū i b i } R m + = /0; (no false equilibria in the sense of Prigogine) M Voronoi finite volume discretization, P boundary conforming Delaunay grid, h = sup n N (t n t n 1 ) <. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 19 (20)
22 Assumptions for uniform results Ω R N is star shaped with respect to some ball B(y 0,R). exp { } R 2 if y y ρ(y) = R 2 y y < R, 0 if y y 0 R. ρ M X(M ): ρk M (x) = min y K ρ(y) for x K. For all finite volume meshes M of Ω we suppose the following properties: Ω ρm (x)dx κ 0, E K E ext /0 implies x K Ω, The geometric weights fulfill property (Q). Ω = I I Ω I finite disjoint union of subdomains s.t. discontinuities of ū i, i = 1,...,m, coincide with subdomain boundaries. (N 1) dimensional measure of all internal subdomain boundaries is bounded. There exists some γ (0,1] such that ū i C 0,γ (Ω I ) := {w ΩI,w C 0,γ (R N )}, i = 1,...,m, I I. Uniform exponential decay of the free energy Dresden, May 26, 2010 Page 20 (20)
Weierstraß Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Discrete Sobolev-Poincaré inequalities for Voronoi finite volume approximations Annegret
More informationConvergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions
Convergence of Finite Volumes schemes for an elliptic-hyperbolic system with boundary conditions Marie Hélène Vignal UMPA, E.N.S. Lyon 46 Allée d Italie 69364 Lyon, Cedex 07, France abstract. We are interested
More informationOn p(x)-laplace thermistor models describing electrothermal feedback in organic semiconductor devices
Weierstrass Institute for Applied Analysis and Stochastics On p(x)-laplace thermistor models describing electrothermal feedback in organic semiconductor devices Matthias Liero joint work with Annegret
More informationSymmetry breaking for a problem in optimal insulation
Symmetry breaking for a problem in optimal insulation Giuseppe Buttazzo Dipartimento di Matematica Università di Pisa buttazzo@dm.unipi.it http://cvgmt.sns.it Geometric and Analytic Inequalities Banff,
More informationWeek 6 Notes, Math 865, Tanveer
Week 6 Notes, Math 865, Tanveer. Energy Methods for Euler and Navier-Stokes Equation We will consider this week basic energy estimates. These are estimates on the L 2 spatial norms of the solution u(x,
More informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
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 informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationNon-Constant Stable Solutions to Reaction- Diffusion Equations in Star-Shaped Domains
Rose-Hulman Undergraduate Mathematics Journal Volume 6 Issue Article 4 Non-Constant Stable Solutions to Reaction- iffusion Equations in Star-Shaped omains Greg rugan University of Texas at Austin, g_drugan@mail.utexas.edu
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationString-representation of hysteresis operators acting on vector-valued, left-continuous and piecewise monotaffine and continuous functions
Weierstrass Institute for Applied Analysis and Stochastics String-representation of hysteresis operators acting on vector-valued, left-continuous and piecewise monotaffine and continuous functions Olaf
More informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
More informationCharge cloud simulations, overview and very first results
W eierstraß -Institut für Angew andte Analysis und Stochastik K. Gärtner Charge cloud simulations, overview and very first results SRH generation 6M pairs, 7.5ns Berlin, July 28, 2008 charge cloud sim.
More informationNumerical Analysis and Methods for PDE I
Numerical Analysis and Methods for PDE I A. J. Meir Department of Mathematics and Statistics Auburn University US-Africa Advanced Study Institute on Analysis, Dynamical Systems, and Mathematical Modeling
More informationUniform-in-time convergence of numerical schemes for Richards and Stefan s models.
Uniform-in-time convergence of numerical schemes for Richards and Stefan s models. Jérôme Droniou, Robert Eymard and Cindy Guichard Abstract We prove that all Gradient Schemes which include Finite Element,
More informationPOINCARÉ S INEQUALITY AND DIFFUSIVE EVOLUTION EQUATIONS
POINCARÉ S INEQUALITY AND DIFFUSIVE EVOLUTION EQUATIONS CLAYTON BJORLAND AND MARIA E. SCHONBEK Abstract. This paper addresses the question of change of decay rate from exponential to algebraic for diffusive
More informationHeat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control
Outline Heat equations with singular potentials: Hardy & Carleman inequalities, well-posedness & control IMDEA-Matemáticas & Universidad Autónoma de Madrid Spain enrique.zuazua@uam.es Analysis and control
More informationStudy of a numerical scheme for miscible two-phase flow in porous media. R. Eymard, Université Paris-Est
Study of a numerical scheme for miscible two-phase flow in porous media R. Eymard, Université Paris-Est in collaboration with T. Gallouët and R. Herbin, Aix-Marseille Université C. Guichard, Université
More informationOn Multigrid for Phase Field
On Multigrid for Phase Field Carsten Gräser (FU Berlin), Ralf Kornhuber (FU Berlin), Rolf Krause (Uni Bonn), and Vanessa Styles (University of Sussex) Interphase 04 Rome, September, 13-16, 2004 Synopsis
More informationNumerical solution to the Complex 2D Helmholtz Equation based on Finite Volume Method with Impedance Boundary Conditions
Open Phys. 016; 14:436 443 Research Article Open Access Angela Handlovičová and Izabela Riečanová* Numerical solution to the Complex D Helmholtz Equation based on Finite Volume Method with Impedance Boundary
More informationHyperbolic Systems of Conservation Laws. in One Space Dimension. II - Solutions to the Cauchy problem. Alberto Bressan
Hyperbolic Systems of Conservation Laws in One Space Dimension II - Solutions to the Cauchy problem Alberto Bressan Department of Mathematics, Penn State University http://www.math.psu.edu/bressan/ 1 Global
More informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
More 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 informationOn the Numerics of 3D Kohn-Sham System
On the of 3D Kohn-Sham System Weierstrass Institute for Applied Analysis and Stochastics Workshop: Mathematical Models for Transport in Macroscopic and Mesoscopic Systems, Feb. 2008 On the of 3D Kohn-Sham
More informationPhase Transitions for Dilute Particle Systems with Lennard-Jones Potential
Weierstrass Institute for Applied Analysis and Stochastics Phase Transitions for Dilute Particle Systems with Lennard-Jones Potential Wolfgang König TU Berlin and WIAS joint with Andrea Collevecchio (Venice),
More informationON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR TERM
Internat. J. Math. & Math. Sci. Vol. 22, No. 3 (999 587 595 S 6-72 9922587-2 Electronic Publishing House ON WEAK SOLUTION OF A HYPERBOLIC DIFFERENTIAL INCLUSION WITH NONMONOTONE DISCONTINUOUS NONLINEAR
More informationConvergence of the MAC scheme for incompressible flows
Convergence of the MAC scheme for incompressible flows T. Galloue t?, R. Herbin?, J.-C. Latche??, K. Mallem????????? Aix-Marseille Universite I.R.S.N. Cadarache Universite d Annaba Calanque de Cortiou
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationThe parabolic Anderson model on Z d with time-dependent potential: Frank s works
Weierstrass Institute for Applied Analysis and Stochastics The parabolic Anderson model on Z d with time-dependent potential: Frank s works Based on Frank s works 2006 2016 jointly with Dirk Erhard (Warwick),
More informationHardy inequalities, heat kernels and wave propagation
Outline Hardy inequalities, heat kernels and wave propagation Basque Center for Applied Mathematics (BCAM) Bilbao, Basque Country, Spain zuazua@bcamath.org http://www.bcamath.org/zuazua/ Third Brazilian
More informationEntropy-dissipation methods I: Fokker-Planck equations
1 Entropy-dissipation methods I: Fokker-Planck equations Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Introduction Boltzmann equation Fokker-Planck equations Degenerate parabolic
More informationSome lecture notes for Math 6050E: PDEs, Fall 2016
Some lecture notes for Math 65E: PDEs, Fall 216 Tianling Jin December 1, 216 1 Variational methods We discuss an example of the use of variational methods in obtaining existence of solutions. Theorem 1.1.
More informationClassical solutions for the quasi-stationary Stefan problem with surface tension
Classical solutions for the quasi-stationary Stefan problem with surface tension Joachim Escher, Gieri Simonett We show that the quasi-stationary two-phase Stefan problem with surface tension has a unique
More informationA Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems
A Priori Error Analysis for Space-Time Finite Element Discretization of Parabolic Optimal Control Problems D. Meidner and B. Vexler Abstract In this article we discuss a priori error estimates for Galerkin
More informationarxiv: v3 [math.ds] 22 Feb 2012
Stability of interconnected impulsive systems with and without time-delays using Lyapunov methods arxiv:1011.2865v3 [math.ds] 22 Feb 2012 Sergey Dashkovskiy a, Michael Kosmykov b, Andrii Mironchenko b,
More informationComputation of homoclinic and heteroclinic orbits for flows
Computation of homoclinic and heteroclinic orbits for flows Jean-Philippe Lessard BCAM BCAM Mini-symposium on Computational Math February 1st, 2011 Rigorous Computations Connecting Orbits Compute a set
More informationLocally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem
56 Chapter 7 Locally convex spaces, the hyperplane separation theorem, and the Krein-Milman theorem Recall that C(X) is not a normed linear space when X is not compact. On the other hand we could use semi
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationDifferent types of phase transitions for a simple model of alignment of oriented particles
Different types of phase transitions for a simple model of alignment of oriented particles Amic Frouvelle Université Paris Dauphine Joint work with Jian-Guo Liu (Duke University, USA) and Pierre Degond
More informationDiscrete entropy methods for nonlinear diffusive evolution equations
Discrete entropy methods for nonlinear diffusive evolution equations 1 Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Joint work with E. Emmrich (TU Berlin), M. Bukal (Zagreb),
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
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 informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationNONLOCAL ELLIPTIC PROBLEMS
EVOLUTION EQUATIONS: EXISTENCE, REGULARITY AND SINGULARITIES BANACH CENTER PUBLICATIONS, VOLUME 52 INSTITUTE OF MATHEMATICS POLISH ACADEMY OF SCIENCES WARSZAWA 2 NONLOCAL ELLIPTIC PROBLEMS ANDRZE J KR
More informationBrownian intersection local times
Weierstrass Institute for Applied Analysis and Stochastics Brownian intersection local times Wolfgang König (TU Berlin and WIAS Berlin) [K.&M., Ann. Probab. 2002], [K.&M., Trans. Amer. Math. Soc. 2006]
More informationPIECEWISE POLYNOMIAL INTERPOLATION IN MUCKENHOUPT WEIGHTED SOBOLEV SPACES AND APPLICATIONS
PIECEWISE POLYNOMIAL INTERPOLATION IN MUCKENHOUPT WEIGHTED SOBOLEV SPACES AND APPLICATIONS RICARDO H. NOCHETTO, ENRIQUE OTÁROLA, AND ABNER J. SALGADO Abstract. We develop a constructive piecewise polynomial
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationOn a weighted total variation minimization problem
On a weighted total variation minimization problem Guillaume Carlier CEREMADE Université Paris Dauphine carlier@ceremade.dauphine.fr Myriam Comte Laboratoire Jacques-Louis Lions, Université Pierre et Marie
More informationGradient Schemes for an Obstacle Problem
Gradient Schemes for an Obstacle Problem Y. Alnashri and J. Droniou Abstract The aim of this work is to adapt the gradient schemes, discretisations of weak variational formulations using independent approximations
More informationFlore Nabet 1. s(c Γ ) εσ b n c, on (0,T) Γ;
ESAI: PROCEEDINGS AND SURVEYS, September 2014, Vol. 45, p. 502-511 J.-S. Dhersin, Editor FINIE-VOLUE EHOD FOR HE CAHN-HILLIARD EQUAION WIH DYNAIC BOUNDARY CONDIIONS Flore Nabet 1 Abstract. A numerical
More informationSt. Petersburg Math. J. 13 (2001),.1 Vol. 13 (2002), No. 4
St Petersburg Math J 13 (21), 1 Vol 13 (22), No 4 BROCKETT S PROBLEM IN THE THEORY OF STABILITY OF LINEAR DIFFERENTIAL EQUATIONS G A Leonov Abstract Algorithms for nonstationary linear stabilization are
More informationThe heat equation in time dependent domains with Neumann boundary conditions
The heat equation in time dependent domains with Neumann boundary conditions Chris Burdzy Zhen-Qing Chen John Sylvester Abstract We study the heat equation in domains in R n with insulated fast moving
More informationarxiv: v1 [math.ap] 24 Oct 2014
Multiple solutions for Kirchhoff equations under the partially sublinear case Xiaojing Feng School of Mathematical Sciences, Shanxi University, Taiyuan 030006, People s Republic of China arxiv:1410.7335v1
More informationTowards pressure-robust mixed methods for the incompressible Navier Stokes equations
Mohrenstrasse 39 10117 Berlin Germany Tel. +49 30 20372 0 www.wias-berlin.de 2016-05-05, Pittsburgh Weierstrass Institute for Applied Analysis and Stochastics Towards pressure-robust mixed methods for
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationOPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES
OPTIMAL CONVERGENCE RATES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS WITH POTENTIAL FORCES RENJUN DUAN Department of Mathematics, City University of Hong Kong 83 Tat Chee Avenue, Kowloon, Hong Kong,
More informationComparison Results for Semilinear Elliptic Equations in Equimeasurable Domains
Comparison Results for Semilinear Elliptic Equations in Equimeasurable Domains François HAMEL Aix-Marseille University & Institut Universitaire de France In collaboration with Emmanuel RUSS Workshop on
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationThreshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations
Threshold behavior and non-quasiconvergent solutions with localized initial data for bistable reaction-diffusion equations P. Poláčik School of Mathematics, University of Minnesota Minneapolis, MN 55455
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
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 informationKramers formula for chemical reactions in the context of Wasserstein gradient flows. Michael Herrmann. Mathematical Institute, University of Oxford
eport no. OxPDE-/8 Kramers formula for chemical reactions in the context of Wasserstein gradient flows by Michael Herrmann Mathematical Institute, University of Oxford & Barbara Niethammer Mathematical
More informationREAL VARIABLES: PROBLEM SET 1. = x limsup E k
REAL VARIABLES: PROBLEM SET 1 BEN ELDER 1. Problem 1.1a First let s prove that limsup E k consists of those points which belong to infinitely many E k. From equation 1.1: limsup E k = E k For limsup E
More informationThe inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method
The inviscid limit to a contact discontinuity for the compressible Navier-Stokes-Fourier system using the relative entropy method Alexis Vasseur, and Yi Wang Department of Mathematics, University of Texas
More informationConvergence and sharp thresholds for propagation in nonlinear diffusion problems
J. Eur. Math. Soc. 12, 279 312 c European Mathematical Society 2010 DOI 10.4171/JEMS/198 Yihong Du Hiroshi Matano Convergence and sharp thresholds for propagation in nonlinear diffusion problems Received
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationA class of non-convex polytopes that admit no orthonormal basis of exponentials
A class of non-convex polytopes that admit no orthonormal basis of exponentials Mihail N. Kolountzakis and Michael Papadimitrakis 1 Abstract A conjecture of Fuglede states that a bounded measurable set
More informationOn the Dependence of Euler Equations on Physical Parameters
On the Dependence of Euler Equations on Physical Parameters Cleopatra Christoforou Department of Mathematics, University of Houston Joint Work with: Gui-Qiang Chen, Northwestern University Yongqian Zhang,
More informationNumerical Simulation of Heat Transfer in Materials with Anisotropic Thermal Conductivity: A Finite Volume Scheme to Handle Complex Geometries
Numerical Simulation of Heat Transfer in Materials with Anisotropic Thermal Conductivity: A Finite Volume Scheme to Handle Complex Geometries Olaf Klein, Jürgen Geiser, Peter Philip 2 Weierstrass Institute
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationNon-radial solutions to a bi-harmonic equation with negative exponent
Non-radial solutions to a bi-harmonic equation with negative exponent Ali Hyder Department of Mathematics, University of British Columbia, Vancouver BC V6TZ2, Canada ali.hyder@math.ubc.ca Juncheng Wei
More informationNUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING TECHNIQUES SATISFYING THE DISCRETE MAXIMUM PRINCIPLE
Proceedings of the Czech Japanese Seminar in Applied Mathematics 2005 Kuju Training Center, Oita, Japan, September 15-18, 2005 pp. 69 76 NUMERICAL SOLUTION OF CONVECTION DIFFUSION EQUATIONS USING UPWINDING
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Proof of a Counterexample to the Finiteness Conjecture in the Spirit of the Theory of
More informationSobolevology. 1. Definitions and Notation. When α = 1 this seminorm is the same as the Lipschitz constant of the function f. 2.
Sobolevology 1. Definitions and Notation 1.1. The domain. is an open subset of R n. 1.2. Hölder seminorm. For α (, 1] the Hölder seminorm of exponent α of a function is given by f(x) f(y) [f] α = sup x
More informationA PROOF OF THE GAUSS-BONNET THEOREM. Contents. 1. Introduction. 2. Regular Surfaces
A PROOF OF THE GAUSS-BONNET THEOREM AARON HALPER Abstract. In this paper I will provide a proof of the Gauss-Bonnet Theorem. I will start by briefly explaining regular surfaces and move on to the first
More informationA REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0002-9939(XX)0000-0 A REMARK ON THE GLOBAL DYNAMICS OF COMPETITIVE SYSTEMS ON ORDERED BANACH SPACES KING-YEUNG LAM
More informationWeierstraß Institut. für Angewandte Analysis und Stochastik. im Forschungsverbund Berlin e.v. Preprint ISSN
Weierstraß Institut für Angewandte Analysis und Stochastik im Forschungsverbund Berlin e.v. Preprint ISSN 0946 8633 Existence of bounded steady state solutions to spin-polarized drift-diffusion systems
More informationTopics in Stochastic Geometry. Lecture 4 The Boolean model
Institut für Stochastik Karlsruher Institut für Technologie Topics in Stochastic Geometry Lecture 4 The Boolean model Lectures presented at the Department of Mathematical Sciences University of Bath May
More informationA mixed finite volume scheme for anisotropic diffusion problems on any grid
A mixed finite volume scheme for anisotropic diffusion problems on any grid Jérôme Droniou 1 and Robert Eymard 2 08/08/2006 Abstract We present a new finite volume scheme for anisotropic heterogeneous
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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
More informationHARDY INEQUALITIES WITH BOUNDARY TERMS. x 2 dx u 2 dx. (1.2) u 2 = u 2 dx.
Electronic Journal of Differential Equations, Vol. 003(003), No. 3, pp. 1 8. ISSN: 107-6691. UL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) HADY INEQUALITIES
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More informationfür Mathematik in den Naturwissenschaften Leipzig
Max-Planck-Institut für Mathematik in den Naturwissenschaften Leipzig A Note on the Poincaré Inequality for Convex Domains by Mario Bebendorf Preprint no.: 24 23 ANoteonthePoincaré Inequality for Convex
More informationINTRODUCTION TO TOPOLOGY, MATH 141, PRACTICE PROBLEMS
INTRODUCTION TO TOPOLOGY, MATH 141, PRACTICE PROBLEMS Problem 1. Give an example of a non-metrizable topological space. Explain. Problem 2. Introduce a topology on N by declaring that open sets are, N,
More informationA one-dimensional nonlinear degenerate elliptic equation
USA-Chile Workshop on Nonlinear Analysis, Electron. J. Diff. Eqns., Conf. 06, 001, pp. 89 99. http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu or ejde.math.unt.edu login: ftp)
More informationSymmetry and monotonicity of least energy solutions
Symmetry and monotonicity of least energy solutions Jaeyoung BYEO, Louis JEAJEA and Mihai MARIŞ Abstract We give a simple proof of the fact that for a large class of quasilinear elliptic equations and
More informationAbstract. 1. Introduction
Journal of Computational Mathematics Vol.28, No.2, 2010, 273 288. http://www.global-sci.org/jcm doi:10.4208/jcm.2009.10-m2870 UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS
More informationLYAPUNOV STABILITY OF CLOSED SETS IN IMPULSIVE SEMIDYNAMICAL SYSTEMS
Electronic Journal of Differential Equations, Vol. 2010(2010, No. 78, pp. 1 18. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu LYAPUNOV STABILITY
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationNon-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3
Non-linear Wave Propagation and Non-Equilibrium Thermodynamics - Part 3 Tommaso Ruggeri Department of Mathematics and Research Center of Applied Mathematics University of Bologna January 21, 2017 ommaso
More informationEXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS WITH UNBOUNDED POTENTIAL. 1. Introduction In this article, we consider the Kirchhoff type problem
Electronic Journal of Differential Equations, Vol. 207 (207), No. 84, pp. 2. ISSN: 072-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE OF SOLUTIONS FOR KIRCHHOFF TYPE EQUATIONS
More informationDiscrete entropy methods for nonlinear diffusive equations
Discrete entropy methods for nonlinear diffusive equations 1 Ansgar Jüngel Vienna University of Technology, Austria www.jungel.at.vu Joint with C. Chainais-Hillairet (Lille), J. A. Carrillo (Imperial),
More informationThe small ball property in Banach spaces (quantitative results)
The small ball property in Banach spaces (quantitative results) Ehrhard Behrends Abstract A metric space (M, d) is said to have the small ball property (sbp) if for every ε 0 > 0 there exists a sequence
More informationAsymptotic behavior of Ginzburg-Landau equations of superfluidity
Communications to SIMAI Congress, ISSN 1827-9015, Vol. 3 (2009) 200 (12pp) DOI: 10.1685/CSC09200 Asymptotic behavior of Ginzburg-Landau equations of superfluidity Alessia Berti 1, Valeria Berti 2, Ivana
More informationRegularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains
Regularity and nonexistence results for anisotropic quasilinear elliptic equations in convex domains Ilaria FRAGALÀ Filippo GAZZOLA Dipartimento di Matematica del Politecnico - Piazza L. da Vinci - 20133
More information