Title of your thesis here
|
|
- Sylvia Sanders
- 5 years ago
- Views:
Transcription
1 Title of your thesis here Your name here A Thesis presented for the degree of Doctor of Philosophy Your Research Group Here Department of Mathematical Sciences University of Durham England Month and Year
2 Dedicated to Someone here
3 Title of your thesis here Your name here Submitted for the degree of Doctor of Philosophy month year here Abstract Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here Your abstract here
4 Declaration The work in this thesis is based on research carried out at the YOUR RESEARCH GROUP HERE, the Department of Mathematical Sciences, England. No part of this thesis has been submitted elsewhere for any other degree or qualification and it is all my own work unless referenced to the contrary in the text. [If your thesis based on joint research, you have to mention what part of it is your individual constribution, see Rules for the Submission of Work for Higher Degrees for detail. You will get one from the Graduate School.] Copyright c 2001 by YOUR NAME HERE. The copyright of this thesis rests with the author. No quotations from it should be published without the author s prior written consent and information derived from it should be acknowledged. iv
5 Acknowledgements Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template.thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. Thank to someone who prepared this template. v
6 Contents Abstract Declaration Acknowledgements iii iv v 1 Introduction 1 2 Evolutionary problem Notation The Existence and Uniqueness of the Continuous Problem A Semidiscrete approximation Notations A fully discrete approximation Scheme Existence and Uniqueness Numerical Experiments Practical Algorithms Iterative Method for Scheme Conclusions 11 Bibliography 12 Appendix 13 vi
7 Contents vii A Basic and Auxiliary Results 13 A.1 Basic Results October 10, 2014
8 List of Figures viii
9 List of Tables ix
10 Chapter 1 Introduction Let Ω be bounded domain in R d (d 3) with Lipschitz boundary Ω. We consider a coupled pair of Cahn-Hilliard Equations modelling a phase separation on a thin film of binary liquid mixture coating substrate, which is wet by one component denoted by A and the other by B: Find {u 1 (x, t), u 2 (x, t)} R R such that where u 1 t = w 1 in Ω, t > 0, (1.0.1a)[1a0001a] u 2 t = w 2 in Ω, t > 0, (1.0.1b)[1a0001b] w 1 = δf (u 1, u 2 ), δu 1 (1.0.1c)[1a0001c] w 2 = δf (u 1, u 2 ), δu 2 (1.0.1d)[1a0001d] F (u 1, u 2 ) = b 1 u 4 1 a 1 u c 1 u b 2 u 4 2 a 2 u c 2 u 2 2 ( ) 2 ( a1 + D u 1 + u 2 + 2b 1 a2 2b 2 ) 2.(1.0.1e)[1a0001e] Here δf (u 1, u 2 )/δu i, for i = 1, 2, indicates the functional derivative. The variable u 1 denotes a local concentration of A or B and u 2 indicates the presence of a liquid or a vapour phase. c i denote the surface tension of u i. The coefficient a i is proportional to T ci T, where T c1 corresponds to the critical temperature of the A-B phase 1
11 Chapter 1. Introduction 2 separation, and T c2 separation. represents the critical temperature of the liquid-vapour phase If a 1 > 0, a 2 > 0, there are two equilibrium phases for each field corresponding to u 1 = ± a1 2b 1 and u 2 = ± a2 2b 2, denote as u + 1, u 1, u + 2, and u 2. The coupling D energetically inhibits the existence of the phase denoted by the (u + 1, u + 2 ). Thus we have three-phase system: liquid A correspond to (u 1, u 2 ) regions, liquid B to (u + 1, u 2 ) regions and the vapour phase to (u 1, u + 2 ) regions. If D = 0, the problem reduces to two decoupled Cahn-Hilliard equations, which has been discussed at length in the mathematical literature. And for this type of problem, we do not have liquid-vapour interfaces. In Chapter 5 two practical algorithms (implicit an explicit methods) for solving the finite element at each time step are suggested. We discuss the convergence theory for the implicit scheme, which used to solve the system arising from Scheme 1. We also discuss in this chapter some computational results for one and two space dimensions. We only use the implicit scheme for all simulations. Before showing some computational results, we discuss linear stability solutions in one space dimension. October 10, 2014
12 Chapter 2 Evolutionary problem In this chapter a global existence and uniqueness theorem for a weak formulation possessing a Lyapunov function is proven. Regularity results are presented for the weak solution. 2.1 Notation Let Ω be a bounded domain in R d, d 3 with boundary Ω. For d = 2, 3 we assume that Ω is Lipschitz boundary. Through out this thesis we adopt the standard notation for Sobolev spaces, denoting the norm of W m,p (Ω) (m N, p [1, ]) by m,p and semi-norm by m,p. For p = 2, W m,p (Ω) will be denoted by H m with the associated norm and semi-norm written as m and m, respectively. In addition we denote the L 2 (Ω) inner product over Ω by (, ) and define the mean of integral η := 1 (η, 1) Ω η L1 (Ω). We also use the following notation, for 1 q <, { L q (0, T ; W m,p (Ω)) := η(x, t) : η(, t) W m,p (Ω), { L (0, T ; W m,p (Ω)) := T 0 } η(, t) q m,p dt <, } η(x, t) : η(, t) W m,p (Ω), ess sup η(, t) m,p < t (0,T ) For later purposes, we recall the Hölder inequality for u L p, v L q 1 < p <, ( uv dx Ω u p dx Ω ) 1 p ( Ω ) 1 v q q 1 dx, where 3 p + 1 q = 1,, and (2.1.1)[Holder]
13 2.2. The Existence and Uniqueness of the Continuous Problem 4 and the following well-known Sobolev interpolation results: Let p [1, ], m 1 ( ) and v W m,p (Ω). Then there are constants C and µ = d 1 1 such that the m p r inequality v 0,r C v 1 µ 0,p v µ m,p, holds for r [p, ] if m d > 0, p [p, ) if m d = 0, p [p, d ] if m d < 0. m d/p p (2.1.2)[gagliar] We also state the following lemma, which will prove useful in our subsequent analysis. Lemma [Lem201] Let u, v, η H 1 (Ω), f = u v, g = u m v n m, m, n = 0, 1, 2, and n m 0. Then for d = 1, 2, 3, fgηdx C u v 0 u m 1 Ω v n m 1 η 1. (2.1.3)[2le000] Proof: Note that using the Cauchy-Schwarz inequality we have u m (u) m v n m 0,2mp v (n m) 0,2(n m)p for n m 0, and m 0, 0,p u m 0,mp or v (n m) 0,(n m)p for m = 0, or n m = 0 respectively. Noting the generalise Hölder inequality and the result above we have fgηdx u v 0 u m v n m 0,3 η 0,6, Ω u 0,6 for m = 2, u v 0 η 0,6 u 0,6 v 0,6 for m = 1, v 2 0,6 for m = 0, C u v 0 u m 1 v n m 1 η 1, where we have noted (2.1.2) to obtain the last inequality. This ends the proof. 2.2 The Existence and Uniqueness of the Continuous Problem Given γ > 0 and u 0 i H 1 (Ω), for i = 1, 2, such that u u C. We consider the problem: October 10, 2014
14 2.2. The Existence and Uniqueness of the Continuous Problem 5 (P) Find {u i, w i } such that u i H 1 (0, T ; (H 1 (Ω)) ) L (0, T ; H 1 (Ω)) for a.e. t (0, T ), w i L 2 (0, T ; H 1 (Ω)) u1 t, η = ( w 1, η), (2.2.1a)[2B0000a] (w 1, η) = (φ(u 1 ), η)+ γ( u 1, η) + 2D(Ψ 1 (u 1, u 2 ), η),(2.2.1b)[2b0000b] u 1 (x, 0) = u 0 1(x), (2.2.1c)[2B0000c] and u2 t, η = ( w 2, η), (2.2.1d)[2B0000d] (w 2, η) = (φ(u 2 ), η)+ γ( u 2, η) + 2D(Ψ 2 (u 1, u 2 ), η),(2.2.1e)[2b0000e] u 2 (x, 0) = u 0 2(x), (2.2.1f)[2B0000f] for all η H 1 (Ω) for a.e. t (0, T ). October 10, 2014
15 Chapter 3 A Semidiscrete approximation In this chapter we introduce some notation which will be used in the current and following chapters. For completeness, we prove interpolation error estimates in the finite element space as these are necessary tools for analysis the current chapter and chapter 4. Then a semidiscrete finite element approximation is proposed where the existence is shown and the uniqueness is proven for one and two dimensions. An error bound between the semidiscrete and continuous solution is given is the final section. 3.1 Notations Let S h H 1 (Ω) be finite element space defined by S h := {χ C( Ω) : χ τ is linear τ T h } H 1 (Ω). Denote by {x i } J i=1 the set of nodes of T h and let {η i } J i=1 be basis for S h defined by η i (x j ) = δ ij, for i, j = 1,..., J. Let π h : C( Ω) S h be the interpolation operator such that π h χ(x i ) = χ(x i ), for i = 1,..., J and define a discrete inner product on C( Ω) as follows J (χ 1, χ 2 ) h := π h (χ 1 (x)χ 2 (x))dx m i χ 1 (x i )χ 2 (x i ), (3.1.1)[3S0000] Ω i=1 where m i = (η i, η i ) h. The induced norm h := [(, ) h ] 1 2 on S h is equivalent to 0 := [(, )] 1 2. Note that the integral (3.1.1) can easily be computed by means of vertex quadrature rule, which exact for piecewise linear functions. 6
16 Chapter 4 A fully discrete approximation In this chapter we introduce a numerical scheme (we call scheme 1) to solve the weak formulation form we mentioned in Chapter 1. We discussed the existence and uniqueness of the solutions for the scheme. We also discuss stability and convergence of the solution to the continuous problem in the weak formulation. We briefly mention a second scheme (we call scheme 2) and show existence, uniqueness, and the stability. We do not discuss the convergence of the scheme. 4.1 Scheme Existence and Uniqueness Given N, a positive integer, let t = T/N denote a fixed time step, and t k = k t where k = 0,..., N. We focus our attention on approximating (P) by the discrete scheme defined as follows: (P h, t 1 ) Given U 0 1, U 0 2, find {U n 1, U n 2, W n 1, W n 2 } S h S h S h S h, for n = 1,..., N, such that η S h ( U 1 n U1 n 1, η) h = ( W1 n, η), (4.1.1a)[4F0008] t (W n 1, η) h = (F 1 (U n 1, U n 2 ),η) h + γ( U n 1, η), (4.1.1b)[4F0009] U 0 1 = P h u 0 1, (4.1.1c)[4F0010] 7
17 4.1. Scheme 1 8 and ( U 2 n U2 n 1, η) h = ( W2 n, η), (4.1.1d)[4F0011] t (W n 2, η) h = (F 2 (U n 1, U n 2 ),η) h + γ( U n 2, η), (4.1.1e)[4F0012] U 0 2 = P h u 0 2, (4.1.1f)[4F0013] where F 1 (U1 n, U2 n ) = (U1 n ) 3 U1 n 1 + D(U1 n + U1 n 1 + 2)(U2 n 1 + 1) (4.1.1g)[4F0014] 2, F 2 (U n 1, U n 2 ) = (U n 2 ) 3 U n D(U n 2 + U n )(U n 1 + 1) 2.(4.1.1h)[4F0015] Note that (4.1.1a c) is independent of U n 2 and (4.1.1d f) is dependent on U n 1. October 10, 2014
18 Chapter 5 Numerical Experiments In this chapter we discuss two practical algorithms (implicit and explicit method) that are used to solve an algebraic system arising in the problem discussed in this thesis. We discuss the convergence theory for the implicit scheme as used to solve the system arising from the scheme 1. We also discuss some computational results for one and two dimensions. We use the implicit scheme in all simulations in this chapter. We have made a comparison with Scheme 2 and the results are compatible. Before showing some computational results, we discuss linear stability solution for the problem. 5.1 Practical Algorithms Iterative Method for Scheme 1 [implicit] Let us expand U i and W i, i = 1, 2, in terms of the standard nodal basis functions of the finite element space S h, that is, U n 1 = U n 2 = J U1,iη n i, W1 n = i=1 J U2,iη n i, W2 n = i=1 J W1,iη n i, i=1 J W2,iη n i, i=1 (5.1.1a)[5E0001a] (5.1.1b)[5E0001b] where J be the number of node points. 9
19 5.1. Practical Algorithms 10 Concluding Remarks To see a clear interaction between the solutions U 1 and U 2 in terms of their physical meaning, it is worthwhile doing computational experiment in three space dimensions. October 10, 2014
20 Chapter 6 Conclusions It was shown using a Faedo Galerkin approximation that there exist a unique solution for the coupled pair of Cahn-Hilliard equations modelling a phase separation on a thin film of binary liquid mixture coating substrate, which is wet by one component. This solution satisfies certain stability bounds. The regularity result was driven at the end of Chapter 2 as having major contribution in obtaining the error bound for the method proposed. 11
21 Bibliography [1] J. W. Barrett and J. F. Blowey (1995), An error bound for the finite element approximation of the Cahn-Hilliard equation with logarithmic free energy, Numerische Mathematics, 72, pp [2] J. W. Barrett and J. F. Blowey (1997), Finite element approximation of a model for phase separation of a multi-component alloy with non-smooth free energy, Numerische Mathematics, 77, pp [3] J. W. Barrett and J. F. Blowey (1999a), An improve error bound for finite element approximation of a model for phase separation of a multi-component alloy, IMA J. Numer. Anal. 19, pp [4] P. G. Ciarlet (1978), The Finite Element Method for Elliptic Problems, North-Holland. [5] J. L. Lions (1969), Quelques Móthodes de Résolution des Problémes aux Limites, Dunod. 12
22 Appendix A Basic and Auxiliary Results A.1 Basic Results 13
AN OPTIMAL ERROR BOUND FOR A FINITE ELEMENT APPROXIMATION OF A MODEL FOR PHASE SEPARATION OF A MULTI-COMPONENT ALLOY WITH NON-SMOOTH FREE ENERGY
Mathematical Modelling and Numerical Analysis MAN, Vol. 33, N o 5, 999, p. 97 987 Modélisation Mathématique et Analyse Numérique AN OPTIMA ERROR BOUND FOR A FINITE EEMENT APPROXIMATION OF A MODE FOR PHASE
More informationNumerical Approximation of Phase Field Models
Numerical Approximation of Phase Field Models Lecture 2: Allen Cahn and Cahn Hilliard Equations with Smooth Potentials Robert Nürnberg Department of Mathematics Imperial College London TUM Summer School
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 informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationAsymptotic Behavior for Semi-Linear Wave Equation with Weak Damping
Int. Journal of Math. Analysis, Vol. 7, 2013, no. 15, 713-718 HIKARI Ltd, www.m-hikari.com Asymptotic Behavior for Semi-Linear Wave Equation with Weak Damping Ducival Carvalho Pereira State University
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 information4 Sobolev spaces, trace theorem and normal derivative
4 Sobolev spaces, trace theorem and normal derivative Throughout, n will be a sufficiently smooth, bounded domain. We use the standard Sobolev spaces H 0 ( n ) := L 2 ( n ), H 0 () := L 2 (), H k ( n ),
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 informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationExistence of minimizers for the pure displacement problem in nonlinear elasticity
Existence of minimizers for the pure displacement problem in nonlinear elasticity Cristinel Mardare Université Pierre et Marie Curie - Paris 6, Laboratoire Jacques-Louis Lions, Paris, F-75005 France Abstract
More informationStationary mean-field games Diogo A. Gomes
Stationary mean-field games Diogo A. Gomes We consider is the periodic stationary MFG, { ɛ u + Du 2 2 + V (x) = g(m) + H ɛ m div(mdu) = 0, (1) where the unknowns are u : T d R, m : T d R, with m 0 and
More informationNumerical Methods for Differential Equations Mathematical and Computational Tools
Numerical Methods for Differential Equations Mathematical and Computational Tools Gustaf Söderlind Numerical Analysis, Lund University Contents V4.16 Part 1. Vector norms, matrix norms and logarithmic
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
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 informationANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE INTO CARTOON PLUS TEXTURE. C.M. Elliott and S.A.
COMMUNICATIONS ON Website: http://aimsciences.org PURE AND APPLIED ANALYSIS Volume 6, Number 4, December 27 pp. 917 936 ANALYSIS OF THE TV REGULARIZATION AND H 1 FIDELITY MODEL FOR DECOMPOSING AN IMAGE
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationMathematical Analysis Outline. William G. Faris
Mathematical Analysis Outline William G. Faris January 8, 2007 2 Chapter 1 Metric spaces and continuous maps 1.1 Metric spaces A metric space is a set X together with a real distance function (x, x ) d(x,
More informationA generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem
A generalised Ladyzhenskaya inequality and a coupled parabolic-elliptic problem Dave McCormick joint work with James Robinson and José Rodrigo Mathematics and Statistics Centre for Doctoral Training University
More informationOn m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry
On m-accretive Schrödinger operators in L p -spaces on manifolds of bounded geometry Ognjen Milatovic Department of Mathematics and Statistics University of North Florida Jacksonville, FL 32224 USA. Abstract
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationThe mortar element method for quasilinear elliptic boundary value problems
The mortar element method for quasilinear elliptic boundary value problems Leszek Marcinkowski 1 Abstract We consider a discretization of quasilinear elliptic boundary value problems by the mortar version
More informationInfinite-Dimensional Dynamical Systems in Mechanics and Physics
Roger Temam Infinite-Dimensional Dynamical Systems in Mechanics and Physics Second Edition With 13 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vii ix GENERAL
More informationON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS. Chun Liu and Jie Shen
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS Volume 7, Number2, April2001 pp. 307 318 ON LIQUID CRYSTAL FLOWS WITH FREE-SLIP BOUNDARY CONDITIONS Chun Liu and Jie Shen Department
More informationMULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
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 informationNonsmooth Schur Newton Methods and Applications
Nonsmooth Schur Newton Methods and Applications C. Gräser, R. Kornhuber, and U. Sack IMA Annual Program Year Workshop Numerical Solutions of Partial Differential Equations: Fast Solution Techniques November
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
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 informationMultigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids
Multigrid Methods for Elliptic Obstacle Problems on 2D Bisection Grids Long Chen 1, Ricardo H. Nochetto 2, and Chen-Song Zhang 3 1 Department of Mathematics, University of California at Irvine. chenlong@math.uci.edu
More informationFinite difference method for elliptic problems: I
Finite difference method for elliptic problems: I Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More informationSobolev spaces. May 18
Sobolev spaces May 18 2015 1 Weak derivatives The purpose of these notes is to give a very basic introduction to Sobolev spaces. More extensive treatments can e.g. be found in the classical references
More informationA PROPERTY OF SOBOLEV SPACES ON COMPLETE RIEMANNIAN MANIFOLDS
Electronic Journal of Differential Equations, Vol. 2005(2005), No.??, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) A PROPERTY
More informationWeak solutions for the Cahn-Hilliard equation with degenerate mobility
Archive for Rational Mechanics and Analysis manuscript No. (will be inserted by the editor) Shibin Dai Qiang Du Weak solutions for the Cahn-Hilliard equation with degenerate mobility Abstract In this paper,
More informationNumerical schemes for short wave long wave interaction equations
Numerical schemes for short wave long wave interaction equations Paulo Amorim Mário Figueira CMAF - Université de Lisbonne LJLL - Séminaire Fluides Compréssibles, 29 novembre 21 Paulo Amorim (CMAF - U.
More informationOn some nonlinear parabolic equation involving variable exponents
On some nonlinear parabolic equation involving variable exponents Goro Akagi (Kobe University, Japan) Based on a joint work with Giulio Schimperna (Pavia Univ., Italy) Workshop DIMO-2013 Diffuse Interface
More informationPart IV: Numerical schemes for the phase-filed model
Part IV: Numerical schemes for the phase-filed model Jie Shen Department of Mathematics Purdue University IMS, Singapore July 29-3, 29 The complete set of governing equations Find u, p, (φ, ξ) such that
More informationSufficient conditions on Liouville type theorems for the 3D steady Navier-Stokes equations
arxiv:1805.07v1 [math.ap] 6 May 018 Sufficient conditions on Liouville type theorems for the D steady Navier-Stokes euations G. Seregin, W. Wang May 8, 018 Abstract Our aim is to prove Liouville type theorems
More informationAppendix A Functional Analysis
Appendix A Functional Analysis A.1 Metric Spaces, Banach Spaces, and Hilbert Spaces Definition A.1. Metric space. Let X be a set. A map d : X X R is called metric on X if for all x,y,z X it is i) d(x,y)
More informationFoundations of Analysis. Joseph L. Taylor. University of Utah
Foundations of Analysis Joseph L. Taylor University of Utah Contents Preface vii Chapter 1. The Real Numbers 1 1.1. Sets and Functions 2 1.2. The Natural Numbers 8 1.3. Integers and Rational Numbers 16
More informationGoal-Oriented Error Estimation for Cahn Hilliard Models of Binary Phase Transition
Goal-Oriented Error Estimation for Cahn Hilliard Models of Binary Phase Transition Kristoffer G. van der Zee, J. Tinsley Oden, Serge Prudhomme, Andrea Hawkins-Daarud Institute for Computational Engineering
More informationExamination paper for TMA4215 Numerical Mathematics
Department of Mathematical Sciences Examination paper for TMA425 Numerical Mathematics Academic contact during examination: Trond Kvamsdal Phone: 93058702 Examination date: 6th of December 207 Examination
More informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationParallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain
Parallel Galerkin Domain Decomposition Procedures for Parabolic Equation on General Domain Keying Ma, 1 Tongjun Sun, 1 Danping Yang 1 School of Mathematics, Shandong University, Jinan 50100, People s Republic
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 informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationCoupled second order singular perturbations for phase transitions
Coupled second order singular perturbations for phase transitions CMU 06/09/11 Ana Cristina Barroso, Margarida Baía, Milena Chermisi, JM Introduction Let Ω R d with Lipschitz boundary ( container ) and
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationResearch Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations
Abstract and Applied Analysis Volume 212, Article ID 391918, 11 pages doi:1.1155/212/391918 Research Article A Two-Grid Method for Finite Element Solutions of Nonlinear Parabolic Equations Chuanjun Chen
More informationEXISTENCE AND REGULARITY OF SOLUTIONS FOR STOKES SYSTEMS WITH NON-SMOOTH BOUNDARY DATA IN A POLYHEDRON
Electronic Journal of Differential Equations, Vol. 2017 (2017), No. 147, pp. 1 10. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu EXISTENCE AND REGULARITY OF SOLUTIONS FOR
More informationON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS Abdelhafid Younsi To cite this version: Abdelhafid Younsi. ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS. 4 pages. 212. HAL Id:
More informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
More informationFinite Volume Discretization of a Class of Parabolic Equations with Discontinuous and Oscillating Coefficients
Applied Mathematical Sciences, Vol, 8, no 3, - 35 HIKARI Ltd, wwwm-hikaricom https://doiorg/988/ams876 Finite Volume Discretization of a Class of Parabolic Equations with Discontinuous and Oscillating
More informationInner Product Spaces
Inner Product Spaces Linear Algebra Josh Engwer TTU 28 October 2015 Josh Engwer (TTU) Inner Product Spaces 28 October 2015 1 / 15 Inner Product Space (Definition) An inner product is the notion of a dot
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 informationarxiv: v2 [math.ap] 28 Nov 2016
ONE-DIMENSIONAL SAIONARY MEAN-FIELD GAMES WIH LOCAL COUPLING DIOGO A. GOMES, LEVON NURBEKYAN, AND MARIANA PRAZERES arxiv:1611.8161v [math.ap] 8 Nov 16 Abstract. A standard assumption in mean-field game
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM PART I: ERROR ANALYSIS UNDER MINIMUM REGULARITIES XIAOBING FENG AND ANDREAS PROHL Abstract. In this first part
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 informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION. Xiaofeng Yang. (Communicated by Jie Shen)
DISCRETE AND CONTINUOUS doi:1.3934/dcdsb.29.11.157 DYNAMICAL SYSTEMS SERIES B Volume 11, Number 4, June 29 pp. 157 17 ERROR ANALYSIS OF STABILIZED SEMI-IMPLICIT METHOD OF ALLEN-CAHN EQUATION Xiaofeng Yang
More informationPreparatory Material for the European Intensive Program in Bydgoszcz 2011 Analytical and computer assisted methods in mathematical models
Preparatory Material for the European Intensive Program in Bydgoszcz 2011 Analytical and computer assisted methods in mathematical models September 4{18 Basics on the Lebesgue integral and the divergence
More information(1) Consider the space S consisting of all continuous real-valued functions on the closed interval [0, 1]. For f, g S, define
Homework, Real Analysis I, Fall, 2010. (1) Consider the space S consisting of all continuous real-valued functions on the closed interval [0, 1]. For f, g S, define ρ(f, g) = 1 0 f(x) g(x) dx. Show that
More informationSobolev Spaces. Chapter 10
Chapter 1 Sobolev Spaces We now define spaces H 1,p (R n ), known as Sobolev spaces. For u to belong to H 1,p (R n ), we require that u L p (R n ) and that u have weak derivatives of first order in L p
More informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
More informationCAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING CONTACT LINES
CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING CONTACT LINES C.G. Gal M Grasselli A Miranville To cite this version: C.G. Gal M Grasselli A Miranville. CAHN-HILLIARD-NAVIER-STOKES SYSTEMS WITH MOVING
More informationVariable Exponents Spaces and Their Applications to Fluid Dynamics
Variable Exponents Spaces and Their Applications to Fluid Dynamics Martin Rapp TU Darmstadt November 7, 213 Martin Rapp (TU Darmstadt) Variable Exponent Spaces November 7, 213 1 / 14 Overview 1 Variable
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationNUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE
NUMERICAL ANALYSIS OF THE CAHN-HILLIARD EQUATION AND APPROXIMATION FOR THE HELE-SHAW PROBLEM, PART II: ERROR ANALYSIS AND CONVERGENCE OF THE INTERFACE XIAOBING FENG AND ANDREAS PROHL Abstract. In this
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationMaximum-norm stability of the finite element Ritz projection with mixed boundary conditions
Noname manuscript No. (will be inserted by the editor) Maximum-norm stability of the finite element Ritz projection with mixed boundary conditions Dmitriy Leykekhman Buyang Li Received: date / Accepted:
More informationLocal flux mimetic finite difference methods
Local flux mimetic finite difference methods Konstantin Lipnikov Mikhail Shashkov Ivan Yotov November 4, 2005 Abstract We develop a local flux mimetic finite difference method for second order elliptic
More informationPREPRINT 2009:20. Finite Element Approximation of the Linearized Cahn-Hilliard-Cook Equation STIG LARSSON ALI MESFORUSH
PREPRINT 9: Finite Element Approximation of the Linearized Cahn-Hilliard-Coo Equation STIG LARSSON ALI MESFORUSH Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY
More informationOverview. A Posteriori Error Estimates for the Biharmonic Equation. Variational Formulation and Discretization. The Biharmonic Equation
Overview A Posteriori rror stimates for the Biharmonic quation R Verfürth Fakultät für Mathematik Ruhr-Universität Bochum wwwruhr-uni-bochumde/num1 Milan / February 11th, 013 The Biharmonic quation Summary
More informationNonexistence of solutions for quasilinear elliptic equations with p-growth in the gradient
Electronic Journal of Differential Equations, Vol. 2002(2002), No. 54, pp. 1 8. ISSN: 1072-6691. URL: http://ejde.math.swt.edu or http://ejde.math.unt.edu ftp ejde.math.swt.edu (login: ftp) Nonexistence
More informationLinear Analysis Lecture 5
Linear Analysis Lecture 5 Inner Products and V Let dim V < with inner product,. Choose a basis B and let v, w V have coordinates in F n given by x 1. x n and y 1. y n, respectively. Let A F n n be the
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 informationhal , version 6-26 Dec 2012
ON THE UNIQUENESS IN THE 3D NAVIER-STOKES EQUATIONS ABDEHAFID YOUNSI Abstract. In this paper, we give a new regularity criterion on the uniqueness results of weak solutions for the 3D Navier-Stokes equations
More informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
More informationTD M1 EDP 2018 no 2 Elliptic equations: regularity, maximum principle
TD M EDP 08 no Elliptic equations: regularity, maximum principle Estimates in the sup-norm I Let be an open bounded subset of R d of class C. Let A = (a ij ) be a symmetric matrix of functions of class
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 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 informationPartial Differential Equations, 2nd Edition, L.C.Evans Chapter 5 Sobolev Spaces
Partial Differential Equations, nd Edition, L.C.Evans Chapter 5 Sobolev Spaces Shih-Hsin Chen, Yung-Hsiang Huang 7.8.3 Abstract In these exercises always denote an open set of with smooth boundary. As
More informationCOSC 3361 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods
COSC 336 Numerical Analysis I Ordinary Differential Equations (II) - Multistep methods Fall 2005 Repetition from the last lecture (I) Initial value problems: dy = f ( t, y) dt y ( a) = y 0 a t b Goal:
More informationError estimates for a finite volume element method for parabolic equations in convex polygonal domains
Error estimates for a finite volume element method for parabolic equations in convex polygonal domains P Chatzipantelidis, R D Lazarov Department of Mathematics, Texas A&M University, College Station,
More informationRemarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable?
Remarks on the analysis of finite element methods on a Shishkin mesh: are Scott-Zhang interpolants applicable? Thomas Apel, Hans-G. Roos 22.7.2008 Abstract In the first part of the paper we discuss minimal
More informationA priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces
A priori error analysis of the BEM with graded meshes for the electric eld integral equation on polyhedral surfaces A. Bespalov S. Nicaise Abstract The Galerkin boundary element discretisations of the
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6.1 Mathematik Preprint Nr. 225 Estimates of the second-order derivatives for solutions to the two-phase parabolic problem
More information2 A Model, Harmonic Map, Problem
ELLIPTIC SYSTEMS JOHN E. HUTCHINSON Department of Mathematics School of Mathematical Sciences, A.N.U. 1 Introduction Elliptic equations model the behaviour of scalar quantities u, such as temperature or
More informationMass Lumping for Constant Density Acoustics
Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial
More informationLecture 17. Higher boundary regularity. April 15 th, We extend our results to include the boundary. Let u C 2 (Ω) C 0 ( Ω) be a solution of
Lecture 7 April 5 th, 004 Higher boundary regularity We extend our results to include the boundary. Higher a priori regularity upto the boundary. Let u C () C 0 ( ) be a solution of Lu = f on, u = ϕ on.
More informationAnalysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient
Analysis of an Adaptive Finite Element Method for Recovering the Robin Coefficient Yifeng Xu 1 Jun Zou 2 Abstract Based on a new a posteriori error estimator, an adaptive finite element method is proposed
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
More informationAnalysis of a high order trace finite element method for PDEs on level set surfaces
N O V E M B E R 2 0 1 6 P R E P R I N T 4 5 7 Analysis of a high order trace finite element method for PDEs on level set surfaces Jörg Grande *, Christoph Lehrenfeld and Arnold Reusken * Institut für Geometrie
More informationPartitioned Methods for Multifield Problems
C Partitioned Methods for Multifield Problems Joachim Rang, 6.7.2016 6.7.2016 Joachim Rang Partitioned Methods for Multifield Problems Seite 1 C One-dimensional piston problem fixed wall Fluid flexible
More informationGLOBAL EXISTENCE AND ENERGY DECAY OF SOLUTIONS TO A PETROVSKY EQUATION WITH GENERAL NONLINEAR DISSIPATION AND SOURCE TERM
Georgian Mathematical Journal Volume 3 (26), Number 3, 397 4 GLOBAL EXITENCE AND ENERGY DECAY OF OLUTION TO A PETROVKY EQUATION WITH GENERAL NONLINEAR DIIPATION AND OURCE TERM NOUR-EDDINE AMROUN AND ABBE
More informationWeak Solutions to Nonlinear Parabolic Problems with Variable Exponent
International Journal of Mathematical Analysis Vol. 1, 216, no. 12, 553-564 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.216.6223 Weak Solutions to Nonlinear Parabolic Problems with Variable
More informationModelling of interfaces and free boundaries
University of Regensburg Regensburg, March 2009 Outline 1 Introduction 2 Obstacle problems 3 Stefan problem 4 Shape optimization Introduction What is a free boundary problem? Solve a partial differential
More informationApplied/Numerical Analysis Qualifying Exam
Applied/Numerical Analysis Qualifying Exam August 9, 212 Cover Sheet Applied Analysis Part Policy on misprints: The qualifying exam committee tries to proofread exams as carefully as possible. Nevertheless,
More informationi=1 α i. Given an m-times continuously
1 Fundamentals 1.1 Classification and characteristics Let Ω R d, d N, d 2, be an open set and α = (α 1,, α d ) T N d 0, N 0 := N {0}, a multiindex with α := d i=1 α i. Given an m-times continuously differentiable
More information