EXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES

Size: px
Start display at page:

Download "EXISTENCE OF SOLUTIONS FOR A CLASS OF VARIATIONAL INEQUALITIES"

Transcription

1 Journal of Matematics and Statistics 9 (4: 35-39, 3 ISSN: doi:.3844/jmssp Publised Online 9 (4 3 (ttp:// EXISTENCE OF SOLUTIONS FOR A CLASS OF ARIATIONAL INEQUALITIES Taalla Frek Department of Matematics, Faculty of Science, Badji Moktar University, Annaba, Algeria Received 3--6; Revised 3--9; Accepted 3--9 ABSTRACT In tis study we considered a deformed elastic solid wit a unilateral contact of a rigid body. We studied te existence, uniqueness and continuity of te deformation of tis solid wit respect to te data. We proved te existence of solutions for a class of variational inequalities. Keywords: ariational Inequalities, Unilateral Contact. INTRODUCTION Several problems in mecanics, pysics, control and tose dealing wit contacts, lead to te study of systems of variational inequalities. Tis model as been studied by Slimane et al. (4; Bernardi et al. (4; Brezis (983; Brezis and Stampaccia (968; Ciarlet (978; Grisvard (985; Haslinger et al. (996; Lions and Stampaccia (967. We consider a solid occupying an open bounded domain Ω of a sufficiently regular boundary Γ = Ω wit unilateral contact wit a rigid obstacle. Teorem. Let P L (ΩR 3 be te resulting of force density. Ten tere exists a unique solution for te variational problem: find u suc tat: Wit: a(u,v = Te bilinear form l(v = Te linear form a ( u,v = l( v, v { ( 3 = v H Ω,R To prove tis teorem we make use of te Lax-Milgram wic is based on proving te continuity and -ellipticity of te bilinear form a(u,v and te continuity of l(v. 35. FORMULATION OF THE CONTACT PROBLEM Here we consider a solid occupying an open bounded domain Ω of a sufficiently regular boundary Γ = Ω. Te solid is supposed to ave: A density on te volume, of force P in Ω Homogenous boundary conditions on Γ Unilateral contact wit a rigid obstacle of equation x 3 = on contact surface Ω c = Ω/Γ. Te displacement is given by: ( u ( x.e3, in Ω c Wit (e, e, e 3 Cartesian base we denote by η te reaction of te obstacle on te solid. Te relations leading to a unilateral contact (witout friction are given by: ( ( u x.e 3, in Ω c η.e 3, in Ω c u x η.e =, in Ω We use te space H ( 3,R 3 c Ω of functions in H (Ω, R 3 equals to zero on Г. Let us introduce te convex subspace K for te autorized displacements, to be defined as:

2 Taalla Frek / Journal of Matematics and Statistics 9 (4: 35-39, 3 { 3 3 c K = v H Ω,R, v.e, in Ω We consider te following variationnal formulation Find: Suc tat: Wit: ( u, η H ( Ω,R 3 H ( Ω ( P ( 3 e a u,v c η,v = l v, v H Ω,R c η,v = ηvdx Ωc And te reduced problem becomes: Find u K suc tat: Teorem. ( P I ( u l( v u a u,v For any solution (u, η of problem (P e, u is a solution of problem (P I. Let (u,η be a solution of problem (P e and u K, v K and we ave: Problem (P e leads to: We assume tat x = : η,v η,v χ η,u, χ K η,u η,u By replacing v by v-u in line one of problem (P e, we get: Were: a ( u,v u c ( η,v u = l( v u c η,v u = η,v u = η,v + η,u a u, v u l v u, v K 36 Let u be a solution of problem (P I ten (u,η is a solution of (P e : a ( u, v u l( v u, v K By using Green's formula, we get: a ( u, v u η, v u l( v u We assume tat v = i±ϕ, wit ϕ D(ΩR 3, (i.e., ϕ is of a compact support, ten te integral on te contour is zero: a ( u, φ = l ( φ, φ Te integral on a contact area leads to: By assuming tat: η,v u, v K v = v = u η,u = And wit te property of convexity of K, we get: Teorem. χ η,u = χ,u η,u = χ,u For any P H (Ω, R 3, te problem (P e as a unique solution (u,η H (ΩR 3 H (Ω Te existence of te solution u of problem is a direct application of Slimane et al. (4. Let us consider: Remark L( v = a ( u,v l( v In problem (P I, we ave: if v =, ten: l( u a u,u if v = u, ten:

3 Taalla Frek / Journal of Matematics and Statistics 9 (4: 35-39, 3 a ( u,u l( u L( u = Te Ker of te form (η, v is caracterized by: { ( 3 3 = v H Ω,R, u.e =, in Ω Let v, ten v and -v are in K from te problem (P I and L(u =, we ave: + ( λ + l( v l( u a ( u,v l( v a ( u,u l( u a ( u, v l( v L( u a u, v a u, u b v, l v l u + = We remplace v by -v in L(u to get: a ( u, v l( v L( u a u, v l v a u, v + l v = L is of a compact support in and from te following inf-sup condition: η,v sup β η v H We can prove tat tere exists η H (Ω. Ten (u, η satisfies line one of problem (P e. Te definition of K and L(u =, leads to: χ η,u = χ,u η,u = χ,u, χ K Tis proves te existence of te solution. Let U and U be two solutions of problem (P I. Wit U = u and U = u ten: ( ( a U,W U l W U, W K a U,W U l W U, W K By adding tat W = U and W = U we ave: ( ( ( + ( λ l( U U a ( U,U U l( U U ( ( a U,U U l U U a U,U U b U U, a U,U U l U U 37 ( a U U,U U a U U,U U U U U = U By te inf-sup condition of problem (P e gives us: v H Ω,R, η,v = η,v η = η 3 3. THE DISCRETE PROBLEM We introduce a discrete subspace of suc tat: { 3 = v C Ω,R, v P k v =,on Ω And dim <, terefore tere exists a basis: {ω i, I = to N, we can ten write: v N = β ω i i i= Now, let us construct a closed convex subset K of suc tat K sould be reduced to a finite number of constraints on te β i : K v, v.e ϕ 3 = at every vertex of eac triangle K Ten K K and K. We remark tat problem (P I is equivalent to find u K suc tat: P a u,v u l v u, v K We assume U = u and W = v. Teorem 3. Let U and U k be te solutions of problems (P I and (P, respectively. Let us denote by A L (, te map defined, by a(u, W = (AU, W, ten: U W + M U U = P AU ' U W + U W

4 Taalla Frek / Journal of Matematics and Statistics 9 (4: 35-39, 3 Wit P is te resultant of te volume force. By te definitions of U and W, we ave: a ( U,U W ( P,U W, W Ka ( U,U W ( P,U W, W K By adding tese inequalities and transposing terms, we obtain: + ( ( ( P,U W a ( U,W a ( U,W a U,U a U,U P,U W Fig.. Mes By subtracting a(u, U + a(u, U from bot sides and grouping terms and by using te continuity and te coercively of te bilinear form a(u, W, we deduce: P AU U W ' ' U U + P AU U W + M U U U W Since: M M U U U W U W Fig.. Isovalue of deplacement u x We obtain: U W + P AU ' M U U = U W + U W W K and W K 4. NUMERICAL RESULTS Consider elastic plate wit te undeformed rectangle sape (, (,. Te body force is te gravity force f and te boundary force g is zero on lower and upper side. On te two vertical sides of te beam are fixed (Fig Fig. 3. Isovalue of deplacement u y 5. CONCLUSION By starting wit te classical model for a deformed elastic solid wit a unilateral contact of a rigid body, we

5 Taalla Frek / Journal of Matematics and Statistics 9 (4: 35-39, 3 proved te existence of solutions for a class of variational inequalities. 6. REFERENCES Bernardi, C., Y. Maday and F. Rapetti, 4. Discrétisations ariationnelles De Problèmes Aux Limites Elliptiques. 3rd Edn., Matématiques et Applications, Springer, ISBN-: , pp: 3. Brezis, H. and G. Stampaccia, 968. Sur la régularité de la solution d'inéquations elliptiques. Bull. Soc. Mat. France, 96: Brezis, H., 983. Analyse fonctionnelle, Collection Matématiques Appliquées pour la Maîtrise. Ciarlet, P.G., 978. Te Finite Element Metod for Elliptic Problems. st Edn., Nort-Holland, Amsterdam, New-York, Oxford, ISBN-: , pp: 53. Grisvard, P., 985. Elliptic Problems in Nonsmoot Domains. st Edn., SIAM, Piladelpia, ISBN-: 69735, pp: 4. Haslinger, J., I. Hlavacek and J. Necas, 996. Numerical metods for unilateral problems in solid mecanics. Handbook Numer. Anal., 4: Lions, J.L. and G. Stampaccia, 967. ariational inequalities. Comm. Pure Appl. Mat., : DOI:./cpa.363 Slimane, L., A. Bendali and P. Laborde, 4. Mixed formulations for a class of variational inequalities. Model. Mat. Anal. Numer., 38: 77-. DOI:.5/man:49 39

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions

Overlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related

More information

Posteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell

Posteriori Analysis of a Finite Element Discretization for a Penalized Naghdi Shell International Journal of Difference Equations ISSN 0973-6069, Volume 8, Number 1, pp. 111 124 (2013) ttp://campus.mst.edu/ijde Posteriori Analysis of a Finite Element Discretization for a Penalized Nagdi

More information

arxiv: v1 [math.na] 11 May 2018

arxiv: v1 [math.na] 11 May 2018 Nitsce s metod for unilateral contact problems arxiv:1805.04283v1 [mat.na] 11 May 2018 Tom Gustafsson, Rolf Stenberg and Jua Videman May 14, 2018 Abstract We derive optimal a priori and a posteriori error

More information

Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method

Different Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles

More information

232 Calculus and Structures

232 Calculus and Structures 3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE

More information

A posteriori error estimation for unilateral contact with matching and non-matching meshes

A posteriori error estimation for unilateral contact with matching and non-matching meshes Comput. Metods Appl. Mec. Engrg. 186 (2000) 65±83 www.elsevier.com/locate/cma A posteriori error estimation for unilateral contact wit matcing and non-matcing meses Patrice Coorevits a, *, Patrick Hild

More information

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for

More information

Mass Lumping for Constant Density Acoustics

Mass 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 information

arxiv: v1 [math.na] 28 Apr 2017

arxiv: v1 [math.na] 28 Apr 2017 THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular

More information

Numerical analysis of a free piston problem

Numerical analysis of a free piston problem MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of

More information

LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS

LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO

More information

LAYER UNDULATIONS IN SMECTIC A LIQUID CRYSTALS

LAYER UNDULATIONS IN SMECTIC A LIQUID CRYSTALS LAYER UNDULATIONS IN SMECTIC A LIQUID CRYSTALS CARLOS J. GARCíA-CERVERA AND SOOKYUNG JOO Abstract. We investigate instabilities of smectic A liquid crystals wen a magnetic field is applied in te direction

More information

MANY scientific and engineering problems can be

MANY scientific and engineering problems can be A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial

More information

A domain decomposition convergence for elasticity equations

A domain decomposition convergence for elasticity equations A domain decomposition convergence for elasticity equations A. Cakib a, A. Ellabib b, A. Nacaoui c, a Département de Matématiques Appliquées et Informatique FST de Beni-Mellal, Université Sultan My Slimane

More information

arxiv: v1 [math.na] 12 Mar 2018

arxiv: v1 [math.na] 12 Mar 2018 ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure

More information

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801

Department of Mathematical Sciences University of South Carolina Aiken Aiken, SC 29801 RESEARCH SUMMARY AND PERSPECTIVES KOFFI B. FADIMBA Department of Matematical Sciences University of Sout Carolina Aiken Aiken, SC 29801 Email: KoffiF@usca.edu 1. Introduction My researc program as focused

More information

Preconditioning in H(div) and Applications

Preconditioning in H(div) and Applications 1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition

More information

c 2004 Society for Industrial and Applied Mathematics

c 2004 Society for Industrial and Applied Mathematics SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares

More information

1. Introduction. We consider the model problem: seeking an unknown function u satisfying

1. Introduction. We consider the model problem: seeking an unknown function u satisfying A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced

More information

arxiv: v1 [math.na] 17 Jul 2014

arxiv: v1 [math.na] 17 Jul 2014 Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order

More information

MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM

MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM MIXED AND STABILIZED FINITE ELEMENT METHODS FOR THE OBSTACLE PROBLEM TOM GUSTAFSSON, ROLF STENBERG, AND JUHA VIDEMAN Abstract. We discretize te Lagrange multiplier formulation of te obstacle problem by

More information

Nonlinear weakly curved rod by Γ -convergence

Nonlinear weakly curved rod by Γ -convergence Noname manuscript No. (will be inserted by te editor Nonlinear weakly curved rod by Γ -convergence Igor Velčić te date of receipt and acceptance sould be inserted later Abstract We present a nonlinear

More information

Poisson Equation in Sobolev Spaces

Poisson Equation in Sobolev Spaces Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on

More information

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim

Math Spring 2013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, (1/z) 2 (1/z 1) 2 = lim Mat 311 - Spring 013 Solutions to Assignment # 3 Completion Date: Wednesday May 15, 013 Question 1. [p 56, #10 (a)] 4z Use te teorem of Sec. 17 to sow tat z (z 1) = 4. We ave z 4z (z 1) = z 0 4 (1/z) (1/z

More information

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.

4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these. Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra

More information

A method for solving first order fuzzy differential equation

A method for solving first order fuzzy differential equation Available online at ttp://ijim.srbiau.ac.ir/ Int. J. Industrial Matematics (ISSN 2008-5621) Vol. 5, No. 3, 2013 Article ID IJIM-00250, 7 pages Researc Article A metod for solving first order fuzzy differential

More information

Differentiation in higher dimensions

Differentiation in higher dimensions Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends

More information

ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000

ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000 ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is

More information

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS

CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION SCHEME WITHOUT FICTITIOUS DOMAINS INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 7, Number 1, Pages 1 9 c 010 Institute for Scientific Computing and Information CELL CENTERED FINITE VOLUME METHODS USING TAYLOR SERIES EXPANSION

More information

ON A PROBLEM OF SHALLOW WATER TYPE

ON A PROBLEM OF SHALLOW WATER TYPE 24-Fez conference on Differential Equations and Mecanics Electronic Journal of Differential Equations, Conference 11, 24, pp. 19 116. ISSN: 172-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu

More information

arxiv: v1 [math.na] 9 Sep 2015

arxiv: v1 [math.na] 9 Sep 2015 arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on

More information

A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions

A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Numerisce Matematik manuscript No. will be inserted by te editor A New Class of Zienkiewicz-Type Nonconforming Element in Any Dimensions Wang Ming 1, Zong-ci Si 2, Jincao Xu1,3 1 LMAM, Scool of Matematical

More information

Math 161 (33) - Final exam

Math 161 (33) - Final exam Name: Id #: Mat 161 (33) - Final exam Fall Quarter 2015 Wednesday December 9, 2015-10:30am to 12:30am Instructions: Prob. Points Score possible 1 25 2 25 3 25 4 25 TOTAL 75 (BEST 3) Read eac problem carefully.

More information

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix

OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra

More information

A = h w (1) Error Analysis Physics 141

A = h w (1) Error Analysis Physics 141 Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.

More information

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders

Some Results on the Growth Analysis of Entire Functions Using their Maximum Terms and Relative L -orders Journal of Matematical Extension Vol. 10, No. 2, (2016), 59-73 ISSN: 1735-8299 URL: ttp://www.ijmex.com Some Results on te Growt Analysis of Entire Functions Using teir Maximum Terms and Relative L -orders

More information

Weak Solutions for Nonlinear Parabolic Equations with Variable Exponents

Weak Solutions for Nonlinear Parabolic Equations with Variable Exponents Communications in Matematics 25 217) 55 7 Copyrigt c 217 Te University of Ostrava 55 Weak Solutions for Nonlinear Parabolic Equations wit Variable Exponents Lingeswaran Sangerganes, Arumugam Gurusamy,

More information

STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES

STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES JEAN-LUC GUERMOND,, JOSEPH E PASCIAK Abstract Using a general approximation setting aving te generic properties of finite-elements, we

More information

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*

ERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER* EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T

More information

Numerical Differentiation

Numerical Differentiation Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function

More information

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness

Key words. Finite element method; convection-diffusion-reaction; nonnegativity; boundedness PRESERVING NONNEGATIVITY OF AN AFFINE FINITE ELEMENT APPROXIMATION FOR A CONVECTION-DIFFUSION-REACTION PROBLEM JAVIER RUIZ-RAMÍREZ Abstract. An affine finite element sceme approximation of a time dependent

More information

INJECTIVE AND PROJECTIVE PROPERTIES OF REPRESENTATIONS OF QUIVERS WITH n EDGES. Sangwon Park

INJECTIVE AND PROJECTIVE PROPERTIES OF REPRESENTATIONS OF QUIVERS WITH n EDGES. Sangwon Park Korean J. Mat. 16 (2008), No. 3, pp. 323 334 INJECTIVE AND PROJECTIVE PROPERTIES OF REPRESENTATIONS OF QUIVERS WITH n EDGES Sanwon Park Abstract. We define injective and projective representations of quivers

More information

A SHORT INTRODUCTION TO BANACH LATTICES AND

A SHORT INTRODUCTION TO BANACH LATTICES AND CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,

More information

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point

1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note

More information

INTRODUCTION AND MATHEMATICAL CONCEPTS

INTRODUCTION AND MATHEMATICAL CONCEPTS Capter 1 INTRODUCTION ND MTHEMTICL CONCEPTS PREVIEW Tis capter introduces you to te basic matematical tools for doing pysics. You will study units and converting between units, te trigonometric relationsips

More information

Convexity and Smoothness

Convexity and Smoothness Capter 4 Convexity and Smootness 4.1 Strict Convexity, Smootness, and Gateaux Differentiablity Definition 4.1.1. Let X be a Banac space wit a norm denoted by. A map f : X \{0} X \{0}, f f x is called a

More information

ANALYTIC SMOOTHING EFFECT FOR THE CUBIC HYPERBOLIC SCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS

ANALYTIC SMOOTHING EFFECT FOR THE CUBIC HYPERBOLIC SCHRÖDINGER EQUATION IN TWO SPACE DIMENSIONS Electronic Journal of Differential Equations, Vol. 2016 2016), No. 34, pp. 1 8. ISSN: 1072-6691. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu ANALYTIC SMOOTHING EFFECT

More information

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015

Math 31A Discussion Notes Week 4 October 20 and October 22, 2015 Mat 3A Discussion Notes Week 4 October 20 and October 22, 205 To prepare for te first midterm, we ll spend tis week working eamples resembling te various problems you ve seen so far tis term. In tese notes

More information

Integral Calculus, dealing with areas and volumes, and approximate areas under and between curves.

Integral Calculus, dealing with areas and volumes, and approximate areas under and between curves. Calculus can be divided into two ke areas: Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and minima problems Integral

More information

Uniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems

Uniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on

More information

Subdifferentials of convex functions

Subdifferentials of convex functions Subdifferentials of convex functions Jordan Bell jordan.bell@gmail.com Department of Matematics, University of Toronto April 21, 2014 Wenever we speak about a vector space in tis note we mean a vector

More information

The Mortar Finite Element Method for Contact Problems

The Mortar Finite Element Method for Contact Problems The Mortar Finite Element Method for Contact Problems F. Ben Belgacem, P. Hild, P. Laborde Mathématiques pour l Industrie et la Physique, Unité Mixte de Recherche CNRS UPS INSAT (UMR 5640), Université

More information

arxiv: v1 [math.dg] 4 Feb 2015

arxiv: v1 [math.dg] 4 Feb 2015 CENTROID OF TRIANGLES ASSOCIATED WITH A CURVE arxiv:1502.01205v1 [mat.dg] 4 Feb 2015 Dong-Soo Kim and Dong Seo Kim Abstract. Arcimedes sowed tat te area between a parabola and any cord AB on te parabola

More information

INTRODUCTION TO CALCULUS LIMITS

INTRODUCTION TO CALCULUS LIMITS Calculus can be divided into two ke areas: INTRODUCTION TO CALCULUS Differential Calculus dealing wit its, rates of cange, tangents and normals to curves, curve sketcing, and applications to maima and

More information

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations

A finite element approximation for the quasi-static Maxwell Landau Lifshitz Gilbert equations ANZIAM J. 54 (CTAC2012) pp.c681 C698, 2013 C681 A finite element approximation for te quasi-static Maxwell Landau Lifsitz Gilbert equations Kim-Ngan Le 1 T. Tran 2 (Received 31 October 2012; revised 29

More information

A SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS

A SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR

More information

Exam 1 Solutions. x(x 2) (x + 1)(x 2) = x

Exam 1 Solutions. x(x 2) (x + 1)(x 2) = x Eam Solutions Question (0%) Consider f() = 2 2 2 2. (a) By calculating relevant its, determine te equations of all vertical asymptotes of te grap of f(). If tere are none, say so. f() = ( 2) ( + )( 2)

More information

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems

5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems 5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we

More information

We name Functions f (x) or g(x) etc.

We name Functions f (x) or g(x) etc. Section 2 1B: Function Notation Bot of te equations y 2x +1 and y 3x 2 are functions. It is common to ave two or more functions in terms of x in te same problem. If I ask you wat is te value for y if x

More information

Decay of solutions of wave equations with memory

Decay of solutions of wave equations with memory Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,

More information

Function Composition and Chain Rules

Function Composition and Chain Rules Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function

More information

COMPOSITION IN FRACTIONAL SOBOLEV SPACES. 1. Introduction. A classical result about composition in Sobolev spaces asserts

COMPOSITION IN FRACTIONAL SOBOLEV SPACES. 1. Introduction. A classical result about composition in Sobolev spaces asserts DISCRETE AND CONTINUOUS Website: ttp://mat.smsu.edu/journal DYNAMICAL SYSTEMS Volume 7, Number, April 001 pp. 41 46 COMPOSITION IN FRACTIONAL SOBOLEV SPACES HAIM BREZIS (1)() AND PETRU MIRONESCU (3) 1.

More information

On convergence of the immersed boundary method for elliptic interface problems

On convergence of the immersed boundary method for elliptic interface problems On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many

More information

Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations

Smoothness of solutions with respect to multi-strip integral boundary conditions for nth order ordinary differential equations 396 Nonlinear Analysis: Modelling and Control, 2014, Vol. 19, No. 3, 396 412 ttp://dx.doi.org/10.15388/na.2014.3.6 Smootness of solutions wit respect to multi-strip integral boundary conditions for nt

More information

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.

Consider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx. Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions

More information

arxiv: v1 [math.na] 27 Jan 2016

arxiv: 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 information

arxiv: v1 [math.na] 19 Mar 2018

arxiv: v1 [math.na] 19 Mar 2018 A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne

More information

How to Find the Derivative of a Function: Calculus 1

How to Find the Derivative of a Function: Calculus 1 Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te

More information

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS

AN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES

More information

Continuity. Example 1

Continuity. Example 1 Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)

More information

Sharp Korn inequalities in thin domains: The rst and a half Korn inequality

Sharp Korn inequalities in thin domains: The rst and a half Korn inequality Sarp Korn inequalities in tin domains: Te rst and a alf Korn inequality Davit Harutyunyan (University of Uta) joint wit Yury Grabovsky (Temple University) SIAM, Analysis of Partial Dierential Equations,

More information

CLOSED CONVEX SHELLS Meunargia T. Mixed forms of stress-strain relations are given in the form. λ + 2µ θ + 1

CLOSED CONVEX SHELLS Meunargia T. Mixed forms of stress-strain relations are given in the form. λ + 2µ θ + 1 Seminar of I. Vekua Institute of Applied Matematics REPORTS, Vol. 43, 207 CLOSED CONVEX SHELLS Meunargia T. Abstract. If Ω is a closed convex sell, ten S : x 3 = 0 is an ovaloid. It is proved tat in tis

More information

The Verlet Algorithm for Molecular Dynamics Simulations

The Verlet Algorithm for Molecular Dynamics Simulations Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical

More information

EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD

EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD

More information

COMPOSITION IN FRACTIONAL SOBOLEV SPACES

COMPOSITION IN FRACTIONAL SOBOLEV SPACES COMPOSITION IN FRACTIONAL SOBOLEV SPACES HAIM BREZIS (1)() AND PETRU MIRONESCU (3) 1. Introduction A classical result about composition in Sobolev spaces asserts tat if u W k,p (Ω) L (Ω) and Φ C k (R),

More information

Continuity and Differentiability of the Trigonometric Functions

Continuity and Differentiability of the Trigonometric Functions [Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te

More information

Sample Problems for Exam II

Sample Problems for Exam II Sample Problems for Exam 1. Te saft below as lengt L, Torsional stiffness GJ and torque T is applied at point C, wic is at a distance of 0.6L from te left (point ). Use Castigliano teorem to Calculate

More information

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set

WYSE Academic Challenge 2004 Sectional Mathematics Solution Set WYSE Academic Callenge 00 Sectional Matematics Solution Set. Answer: B. Since te equation can be written in te form x + y, we ave a major 5 semi-axis of lengt 5 and minor semi-axis of lengt. Tis means

More information

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations

Influence of the Stepsize on Hyers Ulam Stability of First-Order Homogeneous Linear Difference Equations International Journal of Difference Equations ISSN 0973-6069, Volume 12, Number 2, pp. 281 302 (2017) ttp://campus.mst.edu/ijde Influence of te Stepsize on Hyers Ulam Stability of First-Order Homogeneous

More information

2.3 Product and Quotient Rules

2.3 Product and Quotient Rules .3. PRODUCT AND QUOTIENT RULES 75.3 Product and Quotient Rules.3.1 Product rule Suppose tat f and g are two di erentiable functions. Ten ( g (x)) 0 = f 0 (x) g (x) + g 0 (x) See.3.5 on page 77 for a proof.

More information

Math 1241 Calculus Test 1

Math 1241 Calculus Test 1 February 4, 2004 Name Te first nine problems count 6 points eac and te final seven count as marked. Tere are 120 points available on tis test. Multiple coice section. Circle te correct coice(s). You do

More information

Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method

Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using the differential transform method Meccanica 2006) 41:661 670 DOI 10.1007/s11012-006-9012-z Flapwise bending vibration analysis of double tapered rotating Euler Bernoulli beam by using te differential transform metod Ozge Ozdemir Ozgumus

More information

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS

FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XL 2002 FINITE DIFFERENCE APPROXIMATIONS FOR NONLINEAR FIRST ORDER PARTIAL DIFFERENTIAL EQUATIONS by Anna Baranowska Zdzis law Kamont Abstract. Classical

More information

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225

THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Math 225 THE IDEA OF DIFFERENTIABILITY FOR FUNCTIONS OF SEVERAL VARIABLES Mat 225 As we ave seen, te definition of derivative for a Mat 111 function g : R R and for acurveγ : R E n are te same, except for interpretation:

More information

Some Error Estimates for the Finite Volume Element Method for a Parabolic Problem

Some Error Estimates for the Finite Volume Element Method for a Parabolic Problem Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod

More information

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS

AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous

More information

Anisotropic, adaptive finite elements for the computation of a solutal dendrite

Anisotropic, adaptive finite elements for the computation of a solutal dendrite Interfaces and Free Boundaries 5 (23), 13 127 Anisotropic, adaptive finite elements for te computation of a solutal dendrite E. BURMAN AND M. PICASSO Institut de Matématiques, École Polytecnique Fédérale

More information

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations

Global Existence of Classical Solutions for a Class Nonlinear Parabolic Equations Global Journal of Science Frontier Researc Matematics and Decision Sciences Volume 12 Issue 8 Version 1.0 Type : Double Blind Peer Reviewed International Researc Journal Publiser: Global Journals Inc.

More information

The Column and Row Hilbert Operator Spaces

The Column and Row Hilbert Operator Spaces Te Column and Row Hilbert Operator Spaces Roy M Araiza Department of Matematics Purdue University Abstract Given a Hilbert space H we present te construction and some properties of te column and row Hilbert

More information

Computing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach

Computing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New

More information

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households

Volume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of

More information

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract

Superconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of

More information

3. Using your answers to the two previous questions, evaluate the Mratio

3. Using your answers to the two previous questions, evaluate the Mratio MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 0219 2.002 MECHANICS AND MATERIALS II HOMEWORK NO. 4 Distributed: Friday, April 2, 2004 Due: Friday,

More information

Continuity and Differentiability Worksheet

Continuity and Differentiability Worksheet Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;

More information

Chapter 1D - Rational Expressions

Chapter 1D - Rational Expressions - Capter 1D Capter 1D - Rational Expressions Definition of a Rational Expression A rational expression is te quotient of two polynomials. (Recall: A function px is a polynomial in x of degree n, if tere

More information

APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS

APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,

More information

4.2 - Richardson Extrapolation

4.2 - Richardson Extrapolation . - Ricardson Extrapolation. Small-O Notation: Recall tat te big-o notation used to define te rate of convergence in Section.: Definition Let x n n converge to a number x. Suppose tat n n is a sequence

More information

GELFAND S PROOF OF WIENER S THEOREM

GELFAND S PROOF OF WIENER S THEOREM GELFAND S PROOF OF WIENER S THEOREM S. H. KULKARNI 1. Introduction Te following teorem was proved by te famous matematician Norbert Wiener. Wiener s proof can be found in is book [5]. Teorem 1.1. (Wiener

More information

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS

ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS ERROR ESTIMATES FOR A FULLY DISCRETIZED SCHEME TO A CAHN-HILLIARD PHASE-FIELD MODEL FOR TWO-PHASE INCOMPRESSIBLE FLOWS YONGYONG CAI, AND JIE SHEN Abstract. We carry out in tis paper a rigorous error analysis

More information

1. Introduction. We study the numerical analysis of fully nonlinear secondorder elliptic Hamilton Jacobi Bellman (HJB) equations of the form

1. Introduction. We study the numerical analysis of fully nonlinear secondorder elliptic Hamilton Jacobi Bellman (HJB) equations of the form DISCONTINUOUS GALERIN FINITE ELEMENT APPROXIMATION OF HAMILTON JACOBI BELLMAN EQUATIONS WITH CORDÈS COEFFICIENTS IAIN SMEARS AND ENDRE SÜLI Abstract. We propose an p-version discontinuous Galerkin finite

More information