Energy Splitting Theorems for Materials with Memory

Size: px
Start display at page:

Download "Energy Splitting Theorems for Materials with Memory"

Transcription

1 J Elast : DOI /s y Energy Splitting Theorems for Materials with Memory Antonino Favata Paolo Poio-Guiugli Giuseppe Tomassetti Receive: 29 July 2009 / Publishe online: 2 March 2010 Springer Science+Business Meia B.V Abstract We exten to materials with faing memory an materials with internal variables a result previously establishe for materials with instantaneous memory: the aitive ecomposability of the total energy into an internal an a kinetic part, an a representation of the latter an the inertial forces in terms of one an the same mass tensor. Keywors Internal energy Kinetic energy Simple materials Faing memory Internal variables Mathematics Subject Classification A20 74D99 1 Introuction The purpose of this paper is to exten to two classes of materials with memory a result establishe in [6] for materials that, as exemplifie by stanar thermoelastic materials, can only respon to the current values of their state variables. The result we aim to exten is calle in [6] the Energy Splitting Theorem: it is shown that the total energy an the inertia force have consistent representations, uner the assumptions that i the power expeniture of the inertia force be linear in the velocity; an that ii the inertial power plus the rate of change of the energy be translationally invariant. More precisely, it is shown that the energy can be split in two parts, internal an kinetic, with the internal energy inepenent of velocity an the kinetic energy a quaratic form in the A. Favata P. Poio-Guiugli G. Tomassetti Dipartimento i Ingegneria Civile, Università i Roma Tor Vergata, Via Politecnico 1, Rome, Italy tomassetti@ing.uniroma2.it A. Favata favata@ing.uniroma2.it P. Poio-Guiugli ppg@uniroma2.it

2 60 A. Favata et al. velocity, base on a time-inepenent mass tensor, the same that etermines also the workeffective part of the inertial force. The two material classes we here consier are: the class of simple materials in the sense of Truesell an Noll [8], whose mechanical response is etermine by the history of the eformation graient; an the class of materials with internal state variables, as consiere, e.g., by Coleman an Gurtin [1] an Lubliner [5], whose evolution is governe by a generally nonlinear ifferential equation that the Energy Splitting Theorem ha to be extenable to this material class was suggeste by M.E. Gurtin in 1994, on reaing a preprint of [6]. Since these two material classes have a nonempty intersection but o not overlap, we are oblige to prove the entry part of our generalize Energy Splitting Theorem twice; we give the reasons for this at the en of next section. Luckily, as we shall see, the rest of the proof is not as sensitive to the chosen class. Our paper is organize as follows. In Sect. 1, we introuce the quantities that are subject to a constitutive prescription, we stipulate their invariance properties, an we summarize the Energy Splitting Theorem which we aim to generalize. In Sect. 2, we provie a constructive proof of the Energy Splitting Theorem uner the assumption that the constitutive functionals be smooth relative to a norm having the faing-memory property. In Sect. 3, we sketch a proof of the Energy Splitting Theorem for materials with internal variables. Apart for some technicalities that we try to explain as they arise, the structure of the proofs we give is the same as the variant of the proof in [6] givenin[7]. 2 Setting the Stage We work in a referential setting. At any given boy point in this paper, we leave all space epenencies tacit, we introuce scalar volume ensities of the total energy, enote by τ, an of the inertial power: π in = in v, where v is the velocity vector an in is the inertia force vector. Next, we efine the internal power ensity to be α = τ + π in 1 a superpose ot enotes time ifferentiation. Both τ an in are constitutively prescribe at a later stage. For now, it suffices for us to stipulate that, in principle, they both epen on one an the same list of state variables, that we split as follows:, v, where the list inclues only translationally invariant variables. Precisely, a translational change in observer is a mapping leaving the time line unchange: t, x t, x + = t, x +, such that, at some fixe time t, the current shape of the boy uner stuy is pointwise preserve, while the velocity fiel varies by a uniform amount w: x x + = x + t tw, v v + = v + w. 2 Thus, as to state-variable pairs,, v, v + =, v + w.

3 Energy Splitting Theorems for Materials with Memory 61 In the next two sections, we generalize the following result in [6]. Energy Splitting Theorem Let the inertial power be linear in the velocity, in the sense that there is a mapping ˆ 0, referre to as the work-effective inertia force mapping, such that ˆ in, v v = ˆ 0 v, 3 for every state an velocity v. Moreover, let the internal power be invariant uner translational changes in observer: ˆα +, v + = α + = α =ˆα,v. 4 Finally, let the constitutive functions ˆ in, v an ˆτ,v be, respectively, continuous an twice-continuously ifferentiable. Then, the energy τ an the inertial force in have consistent representations, parameterize by i the mass tensor, a symmetric tensor M, inepenent of, v an obeying the mass conservation law Ṁ = 0; ii the internal energy, a scalar-value mapping ˆɛ efine over the state space. These representations are: with D in, v a skew-symmetric tensor. ˆτ,v =ˆɛ + 1 v Mv, 2 5 ˆ in, v = M v + D in, vv, 6 Note that 1, 3, an 2 imply that the invariance requirement 4 takes the form: ˆτ t, ˆvt + w ˆτ t, ˆvt + ˆ 0 t w = 0, 7 t for every constitutive process t t, ˆvt an for every vector w. With this in min, we are in a position to inicate why, in the last part of the Introuction, we state that the entry part of the proof of a theorem of this sort epens on the material class for which it is meant to hol: the first an crucial step in the proof is to achieve a preliminary aitive splitting of the energy into an internal part, that oes not epen on velocity, an a kinetic part. To take that step, it is necessary to compute the erivative of ˆτ with respect to its first argument. This is easy in the case consiere in [6]. It is not so when, as we here o, the constitutive epenence of energy an work-effective inertia force on the current value of is replace either by a functional epenence on the history of up to time t or by the current value of the pair, ˆβ, with ˆβ a solution of the generally nonlinear orinary ifferential equation governing the time evolution of a chosen list of internal variables β: β = f, β. Both replacements entail a rethinking of the structure of state space, as to the accessibility of its points an, more importantly, as to the possibility of giving any process an arbitrary short continuation in time. We will iscuss these technical issues at the appropriate stage of our evelopments.

4 62 A. Favata et al. 3 Materials with Faing Memory Our proof of an Energy Splitting Theorem for faing-memory materials is constructive, an is organize in five steps. Step 1. Translational Invariance of the Internal Power. We let the space of the translationally invariant state variables be a open set C of a finite-imensional inner prouct space L, an we let V enote the velocity space. By a state process an a velocity process ˆv we mean a smooth ifferentiable curve in C an V, respectively. Given a state process, its history up to time t is the mapping t :[0, + C, t s := t s; its past history up to time t is the restriction t r of t to the open half-line 0, + ; an its instantaneous value is t 0 = t; neeless to say, t an t, t r carry the same information. To escribe the energy an inertia-force response of materials with faing memory to state,velocity processes, we introuce three constitutive functionals, τ for the energy an in, 0 for the inertia force; we set: an we consistently rephrase assumption 3: τ v t := τ t, t r, ˆvt, in v t := in t, t r, ˆvt, 8 0 t := 0 t, r t ; in v t ˆvt = 0 t ˆvt. 9 Consequently, the invariance requirement 7 can now be written formally as follows: τv+w t τ v t + 0 t w = 0; 10 t for this to have a precise mathematical sense, we have to specify the regularity of the functionals τ an 0. Step 2. Faing-Memory Property an Chain-Rule Formula. For h : 0, + R + a nonnegative measurable function, chosen once an for all an such that + 0 hs 2 s<+, we enote by L r the Banach space of all measurable functions r : 0, + L, with the norm 1 2 r = hs r s s Two histories are close in the topology etermine by the norm if their values are close in the recent past no matter how far apart they are in the istant past. Thus, since the

5 Energy Splitting Theorems for Materials with Memory 63 constitutive mappings are smooth, changing the history of in the istant past oes not affect appreciably the instantaneous values of τ an in. In Coleman an Noll s terminology [3, 8], introucing the norm enows a material class with the Faing-Memory Property. We let the omain of the functional 0 be an open subset C of the Banach space L = L L r, an we let the common omain of functionals τ an in be C V. Moreover, we require that in be continuous an that τ be twice-continuously Fréchet ifferentiable. As pointe out in Remark 1 of [2], the continuous ifferentiability of τ an the smoothness of t t guarantee that the time erivatives in 10 are well efine, an that the following Chain-Rule Formula hols true: t τ vt = τ t, t r, vt t + δ r τ t, t r, vt[ t r ]+ v τ t, t r, vt vt. Here, t = t an vt = t t ˆvt; moreover, τ an v τ are the partial erivatives of τ with respect to its first an thir argument, respectively, an δ r τ is the unique boune linear functional on L r satisfying τ, t r + r, v = τ, t r, v + δ r τ, t r, v[ r]+o r, for all r in L r. By combining the invariance requirement 10 with the Chain-Rule Formula, we obtain an expression involving both the instantaneous value an the past history of the time erivative of the state process ˆ : τ t, tr, ˆvt + w τ t, tr, ˆvt t + δ r τ t, tr, ˆvt + w δ r τ t, tr, ˆvt [ t r ] + v τ t, tr, ˆvt + w v τ t, tr, ˆvt vt + 0 t, t r w = Step 3. Energy Splitting. As in [6] for instantaneous-memory materials, we achieve the esire result for materials with faing memory by showing that τ oes not epen on v. The argument use in [6] relies on the existence, given a state process an an arbitrary in L, of a state process Ɣ whose instantaneous value coincies with t, an whose rate Ɣt equals ; by replacing with Ɣ, an by invoking the arbitrariness of, one euces that the term multiplying t in the first line of 12 must vanish, an then conclues that τ oes not epen on v. The argument we use in the present proof is similar, an is base on a result ue to Coleman an Mizel [2, Remark 2]: for any given L an L, an for every positive ε, there is a state process ε Ɣ such that whose past history ε Ɣ t r satisfies t ε Ɣt <ε an ε Ɣt =, t r ε Ɣ t r <ε, t r ε Ɣ t r <ε. Replacing with ε Ɣ in 12, letting ε vanish, an using the smoothness of τ an in,we obtain that τ t, tr, ˆvt + w τ t, tr, ˆvt + =0

6 64 A. Favata et al. the ots stan for the remaining terms of 12, whence, by the arbitrariness of, it follows that, for all w in V, τ t, t r, ˆvt + w τ t, t r, ˆvt = 0, i.e., that τ oes not epen on velocity. We conclue that, for every triplet, t r, v of a state, a past history t r,anavelocityv, the constitutive mapping that elivers the total energy splits as follows: τ, t r, v = ɛ, t r + κ t r, v. 13 Step 4. Representation of the Kinetic Energy. In view of 13, relation 12 becomes: δ r κ tr, ˆvt + w δ r κ tr, ˆvt [ t r ] + v κ tr, ˆvt + w v κ tr, ˆvt vt + 0 t, t r w = 0, 14 for all w in V. Moreover, since both mappings δ r τ an v τ are continuously-ifferentiable by assumption, we have that v δ r τ, t r, v[ r] = δ r v τ, t r, v[ r], 15 for every triplet, t r, v of a state, a past history t r,anavelocityv, an for every r in L r. 1 On ifferentiating 14 with respect to w an using 15, we obtain: δ r v κ t r, ˆvt + w[ t r ]+ vv κ t r, ˆvt + w v + 0 t, t r = 0, 16 for all w in V ; upon choosing w = ˆvt in 16, we have that δ r v κ t r, 0[ t r ]+ vv κ t r, 0 vt + 0 t, t r = 0; 17 an, on subtracting 17 from 16 an selecting w = 0, we obtain: δ r v κ tr, ˆvt δ r v κ tr, 0 [ tr ]+ vv κ tr, ˆvt vv κ tr, 0 vt = Finally, since the choice of ˆvt an vt in 18 is arbitrary, we have that 18 itself is equivalent to δ r v κ tr, v δ r v κ tr, 0 [ t r ]=0 19 an vv κ t r, v = M t r, 20 for every past history t r an for every velocity v; note that M t r is a secon-orer symmetric tensor. An integration of 20 yiels: κ t r, v = 1 2 v M t r v + m t r v + μ t r This ientity generalizes the theorem on the inversion of the orer of partial ifferentiation for functions of real variables [4].

7 Energy Splitting Theorems for Materials with Memory 65 We incorporate the constant part μ t r of κ t r, in the internal energy ɛ, t r, without affecting the splitting 13; an we ispose of the linear part by requiring that the total energy be an even function of the velocity, an assumption that, with 13, entails that m t r 0. With these measures, we obtain the following representation for the kinetic energy: an, on combining 13an22, we arrive at κ t r, v = 1 2 v M t r v 22 τ, t r, v = ɛ, t r v M t r v, 23 an aitive energy splitting that generalizes the one in 5. Step 5. Mass Conservation Law an Representation of the Inertial Force. 22 into19, we obtain that, for all v V, By substituting δr M t r v [ t r ]=0, that is, in view also of the Chain Rule Formula, that mass is conserve in the following general sense: t M t r = 0, It remains for us to substitute 22 in17. This yiels: t M t r := δ M r t r [ t r ]. 24 0, t r = M t r v. 25 With 25 an assumption 9 of linearity in the velocity of the power expeniture of the inertia force, we obtain a representation of the inertia force that generalizes 6: in, t r, v = M t r v + D in, t r, vv, 26 where D in, t r, v is a skew-symmetric tensor. We are now in position to summarize our finings uner form of our Energy Splitting Theorem for Materials with Faing Memory Let the constitutive epenence of energy an inertia force be specifie by 8 1,2, with the mappings τ an in presume to be, respectively, twice-continuously ifferentiable an continuous an with the set of past histories enowe with the faing memory norm efine in 11. Assume that: i the inertial power satisfies 9; ii the internal power satisfies the invariance requirement 10; iii the energy mapping is an even function of its thir argument. Then, the mappings τ an in amit the representations 23 an 26, parameterize by a scalar internal-energy mapping ɛ an a symmetric-value mass-tensor mapping M, which obeys the conservation law 24 1.

8 66 A. Favata et al. Remark 1 Assumption iii which is not state in [6] an is not neee to arrive at 20 is instrumental in ruling out the linear term in 21. The same can be achieve by the strengthene invariance requirement that the internal power be invariant uner time-epenent translational changes of observer: if 2 2 is accoringly replace by v v + = v + wt, then, as a bonus, one can also show that the internal energy cannot epen on acceleration, a epenence that woul be incompatible with the Secon Law of Thermoynamics [6]. 4 Materials with Internal Variables For materials with internal variables, the total energy, the inertia force, an the workeffective part of the inertia force are given by ˇτ v t := ˇτ t, ˆβt, ˆvt, ď in v t := ďin t, ˆβt, ˆvt, 27 ď 0 t := ď 0 t, ˆβt, where ˆβt is a translationally-invariant vector of internal variables, whose evolution is governe by the orinary ifferential equation: t ˆβt = f t, ˆβt. 28 By 27 an28, the invariance statement 7 leas to ˇτ t, ˆβt, ˆvt + w ˇτ t, ˆβt, ˆvt t + β ˇτ t, ˆβt, ˆvt + w β ˇτ t, ˆβt, ˆvt f t, ˆβt + v ˇτ t, ˆβt, ˆvt + w v ˇτ t, ˆβt, ˆvt vt + ď 0 t, ˆβt w = By a free-continuation argument borrowe from [6], the time erivative t appearing in the first line of 29 can be replace by an arbitrary rate ; the resulting relation allows one to conclue that ˇτ t, ˆβt, ˆvt + w ˇτ t, ˆβt, ˆvt = 0, for all w in V, whence the energy splitting: ˇτ t, ˆβt, ˆvt =ˇɛ t, ˆβt +ˇκˆβt, ˆvt. From this point on, the proof procees along the steps liste in the previous section, with ˆβt in the place of t r. The conclusions are, mutatis mutanis, the same as those state in

9 Energy Splitting Theorems for Materials with Memory 67 the Energy Splitting Theorem for materials with faing memory: the total energy an the inertial force amit the representations: ˇτ,β, v =ˇɛ,β v ˇMβv, ď in, β, v = ˇMβ v + Ď in, β, vv, with ˇMβ symmetric an Ď in, β, v skew; an the mass-tensor mapping ˇM satisfies the conservation law t ˇM ˆβt = 0. Acknowlegements We thank an anonymous referee for pointing out to us that some justification was neee to rule out the linear term in 21. G.T. gratefully acknowleges the financial support of INAM- GNFM Research Project: Moellazione fisico-matematica i materiali e strutture intelligenti. References 1. Coleman, B.D., Gurtin, M.: Thermoynamics with internal state variables. J. Phys. Chem. 47, Coleman, B.D., Mizel, V.J.: A general theory of issipation in materials with memory. Arch. Ration. Mech. Anal. 27, Coleman, B.D., Noll, W.: An approximation theorem for functionals, with applications in continuum mechanics. Arch. Ration. Mech. Anal. 6, Graves, L.M.: Riemann integration an Taylor s theorem in general analysis. Trans. Am. Math. Soc. 29, Lubliner, J.: On faing memory in materials with evolutionary type. Acta Mech. 8, Poio-Guiugli, P.: Inertia an invariance. Ann. Mat. Pure Appl. CLXXII, Poio-Guiugli, P.: Sparse notes in thermoynamics 2009, forthcoming 8. Truesell, C., Noll, W.: The non-linear fiel theories of mechanics. In: Flugge, S. e. Encyclopeia of Physics, vol. III/3. Springer, Berlin 1965

Abstract A nonlinear partial differential equation of the following form is considered:

Abstract A nonlinear partial differential equation of the following form is considered: M P E J Mathematical Physics Electronic Journal ISSN 86-6655 Volume 2, 26 Paper 5 Receive: May 3, 25, Revise: Sep, 26, Accepte: Oct 6, 26 Eitor: C.E. Wayne A Nonlinear Heat Equation with Temperature-Depenent

More information

A note on the Mooney-Rivlin material model

A note on the Mooney-Rivlin material model A note on the Mooney-Rivlin material moel I-Shih Liu Instituto e Matemática Universiae Feeral o Rio e Janeiro 2945-97, Rio e Janeiro, Brasil Abstract In finite elasticity, the Mooney-Rivlin material moel

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA

TOWARDS THERMOELASTICITY OF FRACTAL MEDIA ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618

More information

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential

A Note on Exact Solutions to Linear Differential Equations by the Matrix Exponential Avances in Applie Mathematics an Mechanics Av. Appl. Math. Mech. Vol. 1 No. 4 pp. 573-580 DOI: 10.4208/aamm.09-m0946 August 2009 A Note on Exact Solutions to Linear Differential Equations by the Matrix

More information

Lie symmetry and Mei conservation law of continuum system

Lie symmetry and Mei conservation law of continuum system Chin. Phys. B Vol. 20, No. 2 20 020 Lie symmetry an Mei conservation law of continuum system Shi Shen-Yang an Fu Jing-Li Department of Physics, Zhejiang Sci-Tech University, Hangzhou 3008, China Receive

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Momentum and Energy. Chapter Conservation Principles

Momentum and Energy. Chapter Conservation Principles Chapter 2 Momentum an Energy In this chapter we present some funamental results of continuum mechanics. The formulation is base on the principles of conservation of mass, momentum, angular momentum, an

More information

Euler equations for multiple integrals

Euler equations for multiple integrals Euler equations for multiple integrals January 22, 2013 Contents 1 Reminer of multivariable calculus 2 1.1 Vector ifferentiation......................... 2 1.2 Matrix ifferentiation........................

More information

Chapter 2 Lagrangian Modeling

Chapter 2 Lagrangian Modeling Chapter 2 Lagrangian Moeling The basic laws of physics are use to moel every system whether it is electrical, mechanical, hyraulic, or any other energy omain. In mechanics, Newton s laws of motion provie

More information

Basic Thermoelasticity

Basic Thermoelasticity Basic hermoelasticity Biswajit Banerjee November 15, 2006 Contents 1 Governing Equations 1 1.1 Balance Laws.............................................. 2 1.2 he Clausius-Duhem Inequality....................................

More information

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1

Assignment 1. g i (x 1,..., x n ) dx i = 0. i=1 Assignment 1 Golstein 1.4 The equations of motion for the rolling isk are special cases of general linear ifferential equations of constraint of the form g i (x 1,..., x n x i = 0. i=1 A constraint conition

More information

Conservation Laws. Chapter Conservation of Energy

Conservation Laws. Chapter Conservation of Energy 20 Chapter 3 Conservation Laws In orer to check the physical consistency of the above set of equations governing Maxwell-Lorentz electroynamics [(2.10) an (2.12) or (1.65) an (1.68)], we examine the action

More information

Conservation laws a simple application to the telegraph equation

Conservation laws a simple application to the telegraph equation J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness

More information

Energy behaviour of the Boris method for charged-particle dynamics

Energy behaviour of the Boris method for charged-particle dynamics Version of 25 April 218 Energy behaviour of the Boris metho for charge-particle ynamics Ernst Hairer 1, Christian Lubich 2 Abstract The Boris algorithm is a wiely use numerical integrator for the motion

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

On the Inclined Curves in Galilean 4-Space

On the Inclined Curves in Galilean 4-Space Applie Mathematical Sciences Vol. 7 2013 no. 44 2193-2199 HIKARI Lt www.m-hikari.com On the Incline Curves in Galilean 4-Space Dae Won Yoon Department of Mathematics Eucation an RINS Gyeongsang National

More information

A vanishing viscosity approach to rate-independent modelling in metric spaces

A vanishing viscosity approach to rate-independent modelling in metric spaces RATE-INDEPENDENT EVOLUTIONS AND MATERIAL MODELING (Special Section of EQUADIFF 2007) Es.: T.Roubíček, U.Stefanelli. Pubblicazione IMATI-CNR, 29PV10/27/0, Pavia, 2010, pp.33 38, ISSN 1772-8964 A vanishing

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

arxiv: v1 [physics.flu-dyn] 8 May 2014

arxiv: v1 [physics.flu-dyn] 8 May 2014 Energetics of a flui uner the Boussinesq approximation arxiv:1405.1921v1 [physics.flu-yn] 8 May 2014 Kiyoshi Maruyama Department of Earth an Ocean Sciences, National Defense Acaemy, Yokosuka, Kanagawa

More information

ELEC3114 Control Systems 1

ELEC3114 Control Systems 1 ELEC34 Control Systems Linear Systems - Moelling - Some Issues Session 2, 2007 Introuction Linear systems may be represente in a number of ifferent ways. Figure shows the relationship between various representations.

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Lecture 2 Lagrangian formulation of classical mechanics Mechanics

Lecture 2 Lagrangian formulation of classical mechanics Mechanics Lecture Lagrangian formulation of classical mechanics 70.00 Mechanics Principle of stationary action MATH-GA To specify a motion uniquely in classical mechanics, it suffices to give, at some time t 0,

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

The Principle of Least Action

The Principle of Least Action Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of

More information

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure

Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure J. Math. Anal. Appl. 341 (28) 457 466 www.elsevier.com/locate/jmaa Exponential asymptotic property of a parallel repairable system with warm stanby uner common-cause failure Zifei Shen, Xiaoxiao Hu, Weifeng

More information

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that 25. More on bivariate functions: partial erivatives integrals Although we sai that the graph of photosynthesis versus temperature in Lecture 16 is like a hill, in the real worl hills are three-imensional

More information

Noether s theorem applied to classical electrodynamics

Noether s theorem applied to classical electrodynamics Noether s theorem applie to classical electroynamics Thomas B. Mieling Faculty of Physics, University of ienna Boltzmanngasse 5, 090 ienna, Austria (Date: November 8, 207) The consequences of gauge invariance

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Lagrangian and Hamiltonian Dynamics

Lagrangian and Hamiltonian Dynamics Lagrangian an Hamiltonian Dynamics Volker Perlick (Lancaster University) Lecture 1 The Passage from Newtonian to Lagrangian Dynamics (Cockcroft Institute, 22 February 2010) Subjects covere Lecture 2: Discussion

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES Communications on Stochastic Analysis Vol. 2, No. 2 (28) 289-36 Serials Publications www.serialspublications.com SINGULAR PERTURBATION AND STATIONARY SOLUTIONS OF PARABOLIC EQUATIONS IN GAUSS-SOBOLEV SPACES

More information

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS

ON ISENTROPIC APPROXIMATIONS FOR COMPRESSIBLE EULER EQUATIONS ON ISENTROPIC APPROXIMATIONS FOR COMPRESSILE EULER EQUATIONS JUNXIONG JIA AND RONGHUA PAN Abstract. In this paper, we first generalize the classical results on Cauchy problem for positive symmetric quasilinear

More information

Generalization of the persistent random walk to dimensions greater than 1

Generalization of the persistent random walk to dimensions greater than 1 PHYSICAL REVIEW E VOLUME 58, NUMBER 6 DECEMBER 1998 Generalization of the persistent ranom walk to imensions greater than 1 Marián Boguñá, Josep M. Porrà, an Jaume Masoliver Departament e Física Fonamental,

More information

Martin Luther Universität Halle Wittenberg Institut für Mathematik

Martin Luther Universität Halle Wittenberg Institut für Mathematik Martin Luther Universität alle Wittenberg Institut für Mathematik Weak solutions of abstract evolutionary integro-ifferential equations in ilbert spaces Rico Zacher Report No. 1 28 Eitors: Professors of

More information

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena

Chaos, Solitons and Fractals Nonlinear Science, and Nonequilibrium and Complex Phenomena Chaos, Solitons an Fractals (7 64 73 Contents lists available at ScienceDirect Chaos, Solitons an Fractals onlinear Science, an onequilibrium an Complex Phenomena journal homepage: www.elsevier.com/locate/chaos

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS

ON THE RIEMANN EXTENSION OF THE SCHWARZSCHILD METRICS ON THE RIEANN EXTENSION OF THE SCHWARZSCHILD ETRICS Valerii Dryuma arxiv:gr-qc/040415v1 30 Apr 004 Institute of athematics an Informatics, AS R, 5 Acaemiei Street, 08 Chisinau, olova, e-mail: valery@ryuma.com;

More information

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS

ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ON THE GEOMETRIC APPROACH TO THE MOTION OF INERTIAL MECHANICAL SYSTEMS ADRIAN CONSTANTIN AND BORIS KOLEV Abstract. Accoring to the principle of least action, the spatially perioic motions of one-imensional

More information

Diagonalization of Matrices Dr. E. Jacobs

Diagonalization of Matrices Dr. E. Jacobs Diagonalization of Matrices Dr. E. Jacobs One of the very interesting lessons in this course is how certain algebraic techniques can be use to solve ifferential equations. The purpose of these notes is

More information

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa

ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS. Gianluca Crippa Manuscript submitte to AIMS Journals Volume X, Number 0X, XX 200X Website: http://aimsciences.org pp. X XX ORDINARY DIFFERENTIAL EQUATIONS AND SINGULAR INTEGRALS Gianluca Crippa Departement Mathematik

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

Tutorial Test 5 2D welding robot

Tutorial Test 5 2D welding robot Tutorial Test 5 D weling robot Phys 70: Planar rigi boy ynamics The problem statement is appene at the en of the reference solution. June 19, 015 Begin: 10:00 am En: 11:30 am Duration: 90 min Solution.

More information

Euler Equations: derivation, basic invariants and formulae

Euler Equations: derivation, basic invariants and formulae Euler Equations: erivation, basic invariants an formulae Mat 529, Lesson 1. 1 Derivation The incompressible Euler equations are couple with t u + u u + p = 0, (1) u = 0. (2) The unknown variable is the

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM

ON THE OPTIMALITY SYSTEM FOR A 1 D EULER FLOW PROBLEM ON THE OPTIMALITY SYSTEM FOR A D EULER FLOW PROBLEM Eugene M. Cliff Matthias Heinkenschloss y Ajit R. Shenoy z Interisciplinary Center for Applie Mathematics Virginia Tech Blacksburg, Virginia 46 Abstract

More information

05 The Continuum Limit and the Wave Equation

05 The Continuum Limit and the Wave Equation Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,

More information

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS

TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS TRAJECTORY TRACKING FOR FULLY ACTUATED MECHANICAL SYSTEMS Francesco Bullo Richar M. Murray Control an Dynamical Systems California Institute of Technology Pasaena, CA 91125 Fax : + 1-818-796-8914 email

More information

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis

The Kepler Problem. 1 Features of the Ellipse: Geometry and Analysis The Kepler Problem For the Newtonian 1/r force law, a miracle occurs all of the solutions are perioic instea of just quasiperioic. To put it another way, the two-imensional tori are further ecompose into

More information

Stable and compact finite difference schemes

Stable and compact finite difference schemes Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long

More information

A. Incorrect! The letter t does not appear in the expression of the given integral

A. Incorrect! The letter t does not appear in the expression of the given integral AP Physics C - Problem Drill 1: The Funamental Theorem of Calculus Question No. 1 of 1 Instruction: (1) Rea the problem statement an answer choices carefully () Work the problems on paper as neee (3) Question

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

Systems & Control Letters

Systems & Control Letters Systems & ontrol Letters ( ) ontents lists available at ScienceDirect Systems & ontrol Letters journal homepage: www.elsevier.com/locate/sysconle A converse to the eterministic separation principle Jochen

More information

Calculus and optimization

Calculus and optimization Calculus an optimization These notes essentially correspon to mathematical appenix 2 in the text. 1 Functions of a single variable Now that we have e ne functions we turn our attention to calculus. A function

More information

On classical orthogonal polynomials and differential operators

On classical orthogonal polynomials and differential operators INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF PHYSICS A: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 38 2005) 6379 6383 oi:10.1088/0305-4470/38/28/010 On classical orthogonal polynomials an ifferential

More information

An extension of Alexandrov s theorem on second derivatives of convex functions

An extension of Alexandrov s theorem on second derivatives of convex functions Avances in Mathematics 228 (211 2258 2267 www.elsevier.com/locate/aim An extension of Alexanrov s theorem on secon erivatives of convex functions Joseph H.G. Fu 1 Department of Mathematics, University

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

More information

A Weak First Digit Law for a Class of Sequences

A Weak First Digit Law for a Class of Sequences International Mathematical Forum, Vol. 11, 2016, no. 15, 67-702 HIKARI Lt, www.m-hikari.com http://x.oi.org/10.1288/imf.2016.6562 A Weak First Digit Law for a Class of Sequences M. A. Nyblom School of

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Sturm-Liouville Theory

Sturm-Liouville Theory LECTURE 5 Sturm-Liouville Theory In the three preceing lectures I emonstrate the utility of Fourier series in solving PDE/BVPs. As we ll now see, Fourier series are just the tip of the iceberg of the theory

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold

CHAPTER 1 : DIFFERENTIABLE MANIFOLDS. 1.1 The definition of a differentiable manifold CHAPTER 1 : DIFFERENTIABLE MANIFOLDS 1.1 The efinition of a ifferentiable manifol Let M be a topological space. This means that we have a family Ω of open sets efine on M. These satisfy (1), M Ω (2) the

More information

2Algebraic ONLINE PAGE PROOFS. foundations

2Algebraic ONLINE PAGE PROOFS. foundations Algebraic founations. Kick off with CAS. Algebraic skills.3 Pascal s triangle an binomial expansions.4 The binomial theorem.5 Sets of real numbers.6 Surs.7 Review . Kick off with CAS Playing lotto Using

More information

arxiv: v1 [math-ph] 2 May 2016

arxiv: v1 [math-ph] 2 May 2016 NONLINEAR HEAT CONDUCTION EQUATIONS WITH MEMORY: PHYSICAL MEANING AND ANALYTICAL RESULTS PIETRO ARTALE HARRIS 1 AND ROBERTO GARRA arxiv:165.576v1 math-ph] May 16 Abstract. We stuy nonlinear heat conuction

More information

arxiv: v1 [math.dg] 1 Nov 2015

arxiv: v1 [math.dg] 1 Nov 2015 DARBOUX-WEINSTEIN THEOREM FOR LOCALLY CONFORMALLY SYMPLECTIC MANIFOLDS arxiv:1511.00227v1 [math.dg] 1 Nov 2015 ALEXANDRA OTIMAN AND MIRON STANCIU Abstract. A locally conformally symplectic (LCS) form is

More information

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES

WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL PRESSURE IN SOBOLEV SPACES Electronic Journal of Differential Equations, Vol. 017 (017), No. 38, pp. 1 7. ISSN: 107-6691. URL: http://eje.math.txstate.eu or http://eje.math.unt.eu WELL-POSEDNESS OF A POROUS MEDIUM FLOW WITH FRACTIONAL

More information

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische

More information

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21

'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21 Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

The canonical controllers and regular interconnection

The canonical controllers and regular interconnection Systems & Control Letters ( www.elsevier.com/locate/sysconle The canonical controllers an regular interconnection A.A. Julius a,, J.C. Willems b, M.N. Belur c, H.L. Trentelman a Department of Applie Mathematics,

More information

A STABILITY STUDY OF NON-NEWTONIAN FLUID FLOWS

A STABILITY STUDY OF NON-NEWTONIAN FLUID FLOWS U.P.B. Sci. Bull., Series A, Vol. 71, Iss. 4, 009 ISSN 13-707 A STABILITY STUDY OF NON-NEWTONIAN FLUID FLOWS Corina CIPU 1, Carmen PRICINĂ, Victor ŢIGOIU 3 Se stuiază problema e curgere a unui flui Olroy

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Chapter 9 Method of Weighted Residuals

Chapter 9 Method of Weighted Residuals Chapter 9 Metho of Weighte Resiuals 9- Introuction Metho of Weighte Resiuals (MWR) is an approimate technique for solving bounary value problems. It utilizes a trial functions satisfying the prescribe

More information

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes

11.7. Implicit Differentiation. Introduction. Prerequisites. Learning Outcomes Implicit Differentiation 11.7 Introuction This Section introuces implicit ifferentiation which is use to ifferentiate functions expresse in implicit form (where the variables are foun together). Examples

More information

A nonlinear inverse problem of the Korteweg-de Vries equation

A nonlinear inverse problem of the Korteweg-de Vries equation Bull. Math. Sci. https://oi.org/0.007/s3373-08-025- A nonlinear inverse problem of the Korteweg-e Vries equation Shengqi Lu Miaochao Chen 2 Qilin Liu 3 Receive: 0 March 207 / Revise: 30 April 208 / Accepte:

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp . ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay

More information

How the potentials in different gauges yield the same retarded electric and magnetic fields

How the potentials in different gauges yield the same retarded electric and magnetic fields How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department

More information

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy IPA Derivatives for Make-to-Stock Prouction-Inventory Systems With Backorers Uner the (Rr) Policy Yihong Fan a Benamin Melame b Yao Zhao c Yorai Wari Abstract This paper aresses Infinitesimal Perturbation

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 24 Winter 207 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Applications of the Wronskian to ordinary linear differential equations

Applications of the Wronskian to ordinary linear differential equations Physics 116C Fall 2011 Applications of the Wronskian to orinary linear ifferential equations Consier a of n continuous functions y i (x) [i = 1,2,3,...,n], each of which is ifferentiable at least n times.

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con

2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con A MIXED FINITE ELEMENT METHOD FOR A THIRD ORDER PARTIAL DIFFERENTIAL EQUATION G. Li 1 an F. A. Milner 2 A mixe finite element metho is escribe for a thir orer partial ifferential equation. The metho can

More information

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is

SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. where L is some constant, usually called the Lipschitz constant. An example is SYSTEMS OF DIFFERENTIAL EQUATIONS, EULER S FORMULA. Uniqueness for solutions of ifferential equations. We consier the system of ifferential equations given by x = v( x), () t with a given initial conition

More information

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS MICHAEL HOLST, EVELYN LUNASIN, AND GANTUMUR TSOGTGEREL ABSTRACT. We consier a general family of regularize Navier-Stokes an Magnetohyroynamics

More information

Polynomial Inclusion Functions

Polynomial Inclusion Functions Polynomial Inclusion Functions E. e Weert, E. van Kampen, Q. P. Chu, an J. A. Muler Delft University of Technology, Faculty of Aerospace Engineering, Control an Simulation Division E.eWeert@TUDelft.nl

More information

Math 180, Exam 2, Fall 2012 Problem 1 Solution. (a) The derivative is computed using the Chain Rule twice. 1 2 x x

Math 180, Exam 2, Fall 2012 Problem 1 Solution. (a) The derivative is computed using the Chain Rule twice. 1 2 x x . Fin erivatives of the following functions: (a) f() = tan ( 2 + ) ( ) 2 (b) f() = ln 2 + (c) f() = sin() Solution: Math 80, Eam 2, Fall 202 Problem Solution (a) The erivative is compute using the Chain

More information

Section 7.2. The Calculus of Complex Functions

Section 7.2. The Calculus of Complex Functions Section 7.2 The Calculus of Complex Functions In this section we will iscuss limits, continuity, ifferentiation, Taylor series in the context of functions which take on complex values. Moreover, we will

More information

Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"

Mark J. Machina CARDINAL PROPERTIES OF LOCAL UTILITY FUNCTIONS Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information