DOCSLIB.ORG

Finite Part Integrals and Hypersingular Kernels

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Finite Part Integrals and Hypersingular Kernels DCDIS A Supplement, Advances in Dynamical Systems, Vol. 14(S2)264 --269 Copyright@2007 Watam Press Finite Part Integrals and Hypersingular Kernels Youn-Sha Chan1, Albert C. Fannjiang2, Glaucio H. Paulino3, and Bao-Feng Feng4 1Department of Computer and Mathematical Sciences, University of Houston–Downtown One Main Street, Houston, Texas 77002, U.S.A., Email: [email protected] 2Department of Mathematics, University of California, Davis, CA 95616, U.S.A. 3Department of Civil and Environmental Engineering, University of Illinois 2209 Newmark Laboratory, 205 North Mathews Avenue, Urbana, IL 61801, U.S.A. 4Department of Mathematics, The University of Texas–Pan American, 1201 W. University Drive, Edinburg, Texas 78541, U.S.A. y Abstract. Singular integral equation method is one of the most effective numerical methods solving a plane crack problem in frac- ture mechanics. Depending on the choice of the density function, very often a higher order of sigularity appears in the equation, and βx we need to give a proper meaning of the integration. In this article G(x) = G 0 e we address the Hadamard finite part integral and how it is used to solve the plane crack problems. Properties of the Hadamard finite part integral will be summarized and compared with other type c d x of integrals. Some numerical results for crack problems by using 2 a Hadamard finite part integral will be provided. Keywords. singular integral equations; Hadamard finite part in- tegrals; fracture mechanics; integral transform method. Figure 1: Antiplane shear problem for a nonhomogeneous AMS (MOS) subject classification: 35, 44, 45, and 73. βx material. Shear modulus G(x) = G0e ; c and d repre- sent the left and right crack tip, respectively; a is the half crack length. 1 Introduction ment vector w(x, y) is In general, the solution to the crack problems in the linear elastic fracture mechanics (LEFM) often leads to a system ∂w(x, y) 2w(x, y) + β = 0 (1) of Cauchy type singular integral equations [4, 5, 16] ∇ ∂x with the mixed boundary conditions a d φ (t) J d i i dt+ k (x, t)φ (t) dt+b φ (x) = p (x) , π − t x ij i i i i w(x, 0) = 0 , x / [c, d] , !c j=1 !c + ∈ (2) − " σyz(x, 0 ) = p(x) , x (c, d) , ∈ where c < x < d, a , b (i = 1, 2, , I) are real constants, where p(x) is the traction function along the crack sur- i i · · · and the kernels kij(x, t) are bounded in the closed domain faces (c, d). By a process of Fourier integral transform (x, t) [c, d] [c, d]. Each function p (x) is known and PDE (1) can be reduced to a hypersingular integral equa- ∈ × i given by the boundary condition(s). Functions φi(x) are tion [2]: the unknowns of the problems, also called by the density G eβx d 1 β functions which often are the derivatives of the displace- 0 = + +N(x, t) w(t, 0+) dt=p(x) , ments. However, if the unknown density function is cho- π (t x)2 2(t x) !c # − − $ sen to be the displacement, say wi(x), then the order of (3) singularity increases. Thus, a formulation of (a system where c < x < d, and N(x, t) takes an integral form of) hypersingular integral equations is made. 2 ∞ β √ξ cos[(t x)ξ] The choice of different unknown density function in the N(x, t)= − − dξ formulation leads to different order of singularity for the !0 2√ξ+ 2 A(ξ, β)+ξ ξ+ A(ξ, β) integral equation. For instance, consider a mode III crack % & () * problem in a nonhomogeneous elastic medium with the ' ' 3 2 shear modulus variation G(x) = G eβx (illustrated in ∞ β B(ξ, β) ξ sin[(t x)ξ]/2 0 − − dξ Figure 1), then the governing a partial differential equa- − & !0 B(ξ, β)'+ξ2 +√2 ξ 2ξ2 +β2 +2 B(ξ, β) tion (PDE) in terms of the z component of the displace- | | % ( &' ) ' * 264 with A(ξ, β) = ξ2 + β2 and B(ξ, β) = ξ4 + β2ξ2. Un- where x > c, and less a proper meaning of integration is given , the first d φ(t) d φ(t) integral in equations (3) is meaningless; the integral is dt := lim dt + φ(x) ln ' , regularized by “Hadamard finite part integral”, and we − t x " 0 t x !x − → +!x+" − , have used “double-bar integral” to denote it. Hadamard finite part (HFP) integral was first introduced by Jacques where x < d. Hadamard [8] to solve some linear PDE, which can be con- sidered as a generalization of the Cauchy principal value 2.2 HFP Integral (CPV) integral [6]. CPV integral does not work for a higher singularity. For instance, consider φ(t) = 1 and x = 0 in 2 HFP and CPV Integrals d φ(t) dt , c < x < d , (7) (t x)2 HFP integral is a generalization of CPV integral, thus let !c − us look at CPV integral first. that is, d dt , c < 0 < d . t2 2.1 CPV Integral !c The integral is not convergent, neither does the principal Equations that involve integrals of the type value exist, since d φ(t) dt 1 1 2 dt , x < 1 (4) = lim + 2 " 0 −c t x | | [c, d] ( ", ") t c − d ' ! − ! \ − → % ( in not integrable in the ordinary (Riemann or Lebesgue is not finite. Hadamard finite part integral is defined by integral) sense because of the kernel 1/(t x) is not in- disregarding the infinite part, 2/' , and keeping the finite tegrable over any interval that includes the− point t = x. part, i.e. d Thus, it is regularized by CPV integral [10, 13]: dt 1 1 = 2 = . (8) c t c − d d x " d ! φ(t) − φ(t) φ(t) dt := lim dt + dt , Definition 1 ( Hadamard finite part integral ) Let ' > − t x " 0 t x t x 0, and denote !c − → +!c − !x+" − , (5) where c < x < d. Notice that the '-neighborhood about F (', x) = f(t, x) dt , c < x < d , [c, d] (x ", x+") the singular point x = t must be symmetric, and it is how ! \ − CPV integral works out for canceling off the singularity. where the singularity appears at the point t = x. If F (', x) For the existence of the CPV integral, the function φ(x) is decomposed into in (5) needs to be at least H¨older continuous on (c, d), that is, φ(x) C0,α(c, d). This requirement of regularity can F (', x) = F0(', x) + F1(', x) , be easily∈checked by following manipulation: and " 0 d lim F0(', x) < , F1(', x) → , φ(t) " 0 ∞ → ∞ dt → − t x then the finite part integral is defined by keeping the “finite !c − φ(t) φ(x) dt part”, i.e. = lim − dt + φ(x) d " 0 t x " t x t x " t x → +!| − |≥ !| − |≥ , − − = f(t, x) dt = lim F0(', x) . d d φ(t) φ(x) dt c " 0 = − dt + φ(x) . (6) ! → c t x −c t x Notice that HFP integral can be considered as a gen- ! − ! − eralization of the CPV integral in the sense that if the Thus, for any φ C0,α, α > 0, the first integral on the ∈ principal value integral exists, then they give the same right side of (6) is an ordinary Riemann integral and the result [6]. We shall define the HFP integral for integrals second integral is with quadratic singularity as in (7). Denote by Cm,α(c, d) d the space of functions whose m-th derivatives are H¨older dt d x continuous on (c, d) with index 0 < α 1. = log − , c < x < d. ≤ −c t x x c ! − − Definition 2 If φ(x) C1,α(c, d), then ∈ Although Cauchy principal value integral is defined for d φ(t) an interior point in (c, d) above, it can be evaluated sep- = dt := (t x)2 arately on both sides of the end points: !c − x " d x x " − φ(t) φ(t) 2φ(x) φ(t) − φ(t) lim 2 dt + 2 dt .(9) dt := lim dt φ(x) ln ' , " 0 / c (t x) x+" (t x) − ' 0 − t x " 0 t x − → ! − ! − !c − → -!c − . 265 The condition φ(x) C1,α(c, d) is required for the exis- For φ Cn,α(c, d) L1+, the first integral on the right tence of the defined ∈HFP integral [12]. side of∈(13) is an ordinary∩ Riemann integral. Also, with Following observation may help to understand Defini- (13) in hand, integration by-parts formula holds for finite tion 2 for HFP. By a step of integration by-parts, the first part integrals [15]. integral under the limit ' 0 in (9) can be written as → Proposition 3 For φ Cn,α(c, d) L1+ x " x " ∈ ∩ − φ(t) φ(x ') c − φ%(t) d d φ%(t) φ(t) 2 dt = − + dt . c (t x) ' − c x c t x) = n dt = n = n+1 dt ! − − ! − c (t x) c (t x) ! − ! − Similarly, φ(d) φ(c) + n n , n 1 d d (d x) − (c x) ≥ φ(t) φ(x + ') d φ%(t) − − dt = + dt . (t x)2 ' − d x t x) and for φ Cα(c, d) L1+ !x+" − − !x+" − ∈ ∩ Thus, the term 2φ(x)/' in (9) will kill the singularity d − φ%(t) log t x dt [φ(x ') + φ(x + ')]/', and under the assumption that c | − | − 1,α ! φ(x) C (c, d) Definition 2 indeed takes the finite part d φ(t) of the∈integral according to Definition 1.
Load more

Recommended publications
  • The Hilbert Transform
    (February 14, 2017) The Hilbert transform Paul Garrett [email protected] http:=/www.math.umn.edu/egarrett/ [This document is http:=/www.math.umn.edu/egarrett/m/fun/notes 2016-17/hilbert transform.pdf] 1. The principal-value functional 2. Other descriptions of the principal-value functional 3. Uniqueness of odd/even homogeneous distributions 4. Hilbert transforms 5. Examples 6. ... 7. Appendix: uniqueness of equivariant functionals 8. Appendix: distributions supported at f0g 9. Appendix: the snake lemma Formulaically, the Cauchy principal-value functional η is Z 1 f(x) Z f(x) ηf = principal-value functional of f = P:V: dx = lim dx −∞ x "!0+ jxj>" x This is a somewhat fragile presentation. In particular, the apparent integral is not a literal integral! The uniqueness proven below helps prove plausible properties like the Sokhotski-Plemelj theorem from [Sokhotski 1871], [Plemelji 1908], with a possibly unexpected leading term: Z 1 f(x) Z 1 f(x) lim dx = −iπ f(0) + P:V: dx "!0+ −∞ x + i" −∞ x See also [Gelfand-Silov 1964]. Other plausible identities are best certified using uniqueness, such as Z f(x) − f(−x) ηf = 1 dx (for f 2 C1( ) \ L1( )) 2 x R R R f(x)−f(−x) using the canonical continuous extension of x at 0. Also, 2 Z f(x) − f(0) · e−x ηf = dx (for f 2 Co( ) \ L1( )) x R R R The Hilbert transform of a function f on R is awkwardly described as a principal-value integral 1 Z 1 f(t) 1 Z f(t) (Hf)(x) = P:V: dt = lim dt π −∞ x − t π "!0+ jt−xj>" x − t with the leading constant 1/π understandable with sufficient hindsight: we will see that this adjustment makes H extend to a unitary operator on L2(R).
    [Show full text]
  • An Extension of Hewitt's Inversion Formula and Its Application to Fluctuation Theory
    AN EXTENSION OF HEWITT’S INVERSION FORMULA AND ITS APPLICATION TO FLUCTUATION THEORY S¸erban E. B˘adil˘a1 Abstract. We analyze fluctuations of random walks with generally distributed increments. Integral representations for key performance measures are obtained by extending an inversion theorem of Hewitt [11] for Laplace-Stieltjes transforms. Another important step in the analysis involves introducing the so-called harmonic measures associated to the walk. It is also pointed out that such representations can be explicitly calculated, if one assumes a form of rational structure for the increment transform. Applications include, but are not restricted to, queueing and insurance risk problems. Keywords: fluctuations of random walks, harmonic analysis, singular integrals 2010 Mathematics Subject Classification: Primary 60G50 Secondary 30E20 1. Introduction A standard assumption in the theory of fluctuations of random walks as they appear in, e.g., queueing or insurance applications, is that the increment of the random walk can be represented as the difference of two independent random variables. In the context of queueing theory, random walks with increments which do not have this property appear as embedded at arrival epochs of customers in a single server queue in which the service requirement of the current customer is correlated with the time until the next arrival. In risk reserve processes that appear in insurance, the independence assumption is violated when the current claim size depends on the time elapsed since the previous arrival, and hence on the premium gained meanwhile. The purpose of the present paper is to show how much can be done for random walks which do not satisfy this independence assumption, regarding their maxima (as in the waiting time/maximum aggregate loss), their minima (idle periods/deficit at ruin) and their excur- arXiv:1508.00751v1 [math.PR] 4 Aug 2015 sions which are related to busy periods or the time to ruin.
    [Show full text]
  • 1 Introduction
    Tel Aviv University, 2015 Functions of real variables 1 1 Introduction 1a Integral buyer's guide: know what type is better 1 1b Dimension 2 is a different story . 3 1c Gauge: a bridge from Riemann to Lebesgue . 6 1d Small but dense . 8 1e Sandwich (Riemann), zigzag sandwich (Lebesgue) 9 1f Some achievements of Lebesgue theory . 10 1g Integral or measure? Functions or sets? . 11 This lightweight introduction is not a part of the theory of Lebesgue in- tegration. It is rather a \comparative anatomy" of integration theories1 in- tended to explain what distinguishes Lebesgue's theory from others, and why it is the best but not all-mighty. 1a Integral buyer's guide: know what type is better Several types of integral are well-known: the Riemann integral; ∗ the improper Riemann integral (two types: 1-dim and n-dim); ∗ the Lebesgue integral; ∗ the gauge integral; ∗ the Newton integral (that is, antiderivative). ∗ Probably, the most videly used in practice is the improper Riemann inte- gral; the most well-known to a wide audience is the antiderivative; the most teached are the Riemann and Lebesgue integrals; and the gauge integral is controversial, seldom used, seldom teached, but contains a bright idea able to illuminate the Lebesgue integral. This course is devoted to the Lebesgue integral, but relations to other integrals are also of some interest. For each type of integral we have the definition of integrable function; ∗ the definition of the integral of an integrable function; ∗ theorems on properties of the integral. ∗ 1This phrasing is borrowed from S. Berberian, see the quote on pp.
    [Show full text]
  • Problem Set 4. Generalized Functions
    Problem Set 4. GeneralizedProblem functions.set 3 MathMath 212a 212a OctoberSept. 30, 2, 2014,2008 due Due Oct. Oct. 15 1 The purpose of this exercise set, which is mainly computational, is to give you some experience with the ideas associated to generalized functions. Recall that in our study of the Fourier transform we introduced the space consisting of all functions on R which have derivatives of all order which vanishS rapidly at infinity, and that on this space we have a countable family of norms m n f m,n := x D f , ∞ where m and n are non-negative integers. Convergence in means convergence with respect to all of these norms. So a linear functionS on this space is continuous if (fk) 0 whenever fk m,n 0 for every m and n.Welet denote this dual| space,| → the space of continuous → linear functions on . S For example, consider the “Heaviside function” S 1 for x>0 H(x):= . 0 for x 0 ≤ This function does not belong to , but it defines a continuous linear function on by S S ∞ ∞ H, f = f(x)H(x)dx = f(x)dx. 0 −∞ Of course, any element of certainly defines a linear function on by the same procedure, the linear functionS associated to g is just the HilbertS space scalar product ∈ S ∞ f (f,g)= f(x)g(x)dx. → −∞ This also works for any g L2(R), but H L2(R). From the Riesz representa- tion theorem we thus know∈ that the linear function∈ given by H is not continuous with respect to the L2 norm, but it is continuous relative to the topology given above by the countable family of norms .
    [Show full text]
  • Cauchy Principal Value Contour Integral with Applications
    entropy Article Cauchy Principal Value Contour Integral with Applications Matilde Legua 1 and Luis M. Sánchez-Ruiz 2,* 1 EINA-Departamento de Matemática Aplicada, Universidad de Zaragoza, Zaragoza 50018, Spain; [email protected] 2 CITG-Departamento de Matemática Aplicada, Universitat Politècnica de València, Valencia 46022, Spain * Correspondence: [email protected]; Tel.: +34-963-877-665 Academic Editors: Carlo Cattani and Kevin Knuth Received: 28 March 2017; Accepted: 3 May 2017; Published: 10 May 2017 Abstract: Cauchy principal value is a standard method applied in mathematical applications by which an improper, and possibly divergent, integral is measured in a balanced way around singularities or at infinity. On the other hand, entropy prediction of systems behavior from a thermodynamic perspective commonly involves contour integrals. With the aim of facilitating the calculus of such integrals in this entropic scenario, we revisit the generalization of Cauchy principal value to complex contour integral, formalize its definition and—by using residue theory techniques—provide an useful way to evaluate them. Keywords: Cauchy principal value; contour integral; entropy as measurement; information extraction; thermodynamics; aerodynamics 1. Introduction Continuous intrinsic entropy treatment requires the assistance of measure theory. This involves adequate integrals intended to evaluate the unexpected uncertainty, indeed this is the case when differential entropy is studied [1–5] in which case improper integrals are commonly used. A number of analytical and numerical methods have been developed such as the Glauert and Goldstein methods developed to study improper integrals that appear in this setting, [6–11], which are relevant in aerodynamics. In fact, aerodynamic entropy study focuses on the overall performance neglecting singularities that enjoy some symmetric entropic phenomena such as turbulence on wing or rotor blades and for this Cauchy principal values are commonly used [12–14].
    [Show full text]
  • Advanced Mathematical Methods in Theoretical Physics
    Advanced Mathematical Methods in Theoretical Physics Wintersemester 2017/2018 Technische Universit¨atBerlin PD Dr. Gernot Schaller March 15, 2021 i This lecture aims at providing physics students and neighboring disciplines with a heuristic toolbox that can be used to tackle all kinds of problems. As time and resources are limited, we will explore only small fractions of different much wider fields, and students are invited to take the lecture only as a starting point. The lecture will try to be as self-contained as possible and aims at providing rather recipes than strict mathematical proofs. It is not an original contribution but an excerpt of many papers, books, other lectures, and own experiences. Only a few of these have been included as references (footnotes describing famous scientists originate from Wikipedia and solely serve for better identification). As successful learning requires practice, the lecture (Mondays) will be accompanied by exer- cises. Students should turn in individual solutions to the exercises at the beginning of Wednesday seminars (computer algebra may be used if applicable). Students may earn points for the exercise sheets, to be admitted to the final exam, they should have earned 50% of the points. After passing the final exam, the students earn in total six ECTS credit points. To get these credit points, students therefore need to get 50 % of the exercise points and 50 % to pass the final exam. The lecture script will be made available online at http://www1.itp.tu-berlin.de/schaller/. Corrections and suggestions for improvements should be addressed to [email protected].
    [Show full text]