arxiv: v3 [math.ap] 7 May 2016
|
|
- Ernest Hall
- 6 years ago
- Views:
Transcription
1 Complex order fractional derivatives in viscoelasticity arxiv: v3 [math.ap] 7 May 16 Teodor M. Atanacković Sanja Konjik Stevan Pilipović Dušan Zorica Abstract We introduce complex order fractional derivatives in models that describe viscoelastic materials. This can not be carried out unrestrictedly, and therefore we derive, for the first time, real valued compatibility constraints, as well as physical constraints that lead to acceptable models. As a result, we introduce a new form of complex order fractional derivative. Also, we consider a fractional differential equation with complex derivatives, and study its solvability. Results obtained for stress relaxation and creep are illustrated by several numerical examples. Mathematics Subject Classification 1): Primary: 6A33; Secondary: 74D5 Keywords: real and complex order fractional derivatives, constitutive equations, the Laplace transform, the Fourier transform, thermodynamical restrictions 1 Introduction Fractional calculus is a powerful tool for modeling various phenomena in mechanics, physics, biology, chemistry, medicine, economy, etc. Last few decades have brought a rapid expansion of the non-integer order differential and integral calculus, from which both the theory and its applications benefit significantly. However, most of the work done in this field so far has been based on the use of real order fractional derivatives and integrals. It is worth to mention that there are several authors who also applied complex order fractional derivatives to model various phenomena, see the work of Machado or Makris, [, 3, 4]. In all of these papers, restrictions on constitutive parameters that follow from the Second Law of Thermodynamics were not examined. In the analysis that follows this issue will be addressed. The main goal of this paper is to motivate and explain basic concepts of fractional calculus with complex order fractional derivatives. Throughout the paper we will investigate constitutive equations in the dimensionless form, for all independent t and x) and dependent σ Faculty of Technical Sciences, Institute of Mechanics, University of Novi Sad, Trg D. Obradovića 6, 1 Novi Sad, Serbia, Electronic mail: atanackovic@uns.ac.rs Faculty of Sciences, Department of Mathematics and Informatics, University of Novi Sad, Trg D. Obradovića 4, 1 Novi Sad, Serbia, Electronic mail: sanja.konjik@dmi.uns.ac.rs Faculty of Sciences, Department of Mathematics and Informatics, University of Novi Sad, Trg D. Obradovića 4, 1 Novi Sad, Serbia, Electronic mail: pilipovic@dmi.uns.ac.rs Institute of Mathematics, Serbian Academy of Sciences and Arts, Kneza Mihaila 36, 11 Belgrade, Serbia, Electronic mail: dusan zorica@mi.sanu.ac.rs 1
2 and ε) variables. Thus, consider a constitutive equation, given by 1), connecting the stress σt, x) at the point x R and time t R +, with the strain εt, x): N n= a n D αn t σt, x) = M m= b m D βm t εt, x), 1) that contains fractional derivatives of complex order α 1,..., α N, β 1,..., β M. The precise definition of the operator D η t of fractional differentiation with respect to t is given below. In order to make a useful framework for the study of 1) we involve two types of conditions: 1. real valued compatibility constraints, and. thermodynamical constraints. Since this paper deals only with the well-posedness of constitutive equations of type 1) and their solvability for strain if stress is prescribed, we may, without loss of generality, assume that both σ and ε are functions only of t. Also, equation 1) can be seen as a generalization of different models considered in the literature so far see e.g. [5, 7, 8, 14, 16, 17,, 5]), since by taking all α n and β m to be real numbers, the problem is reduced to the real case studied in the mentioned papers. Our results will show that constitutive equations of the form 1) may lead to creep and stress relaxation curves that are not monotonic. We note that the conditions of complete monotonicity required in e.g. [1, 13, 1] turn out not to be necessary but rather sufficient to meet the thermodynamic restrictions. It seems that conditions of monotonicity are stronger than the conditions following from the Second Law of Thermodynamics in the form of [6] analogously as in [1] where it is shown that the asymptotic stability, for certain class of constitutive equations, requires more extensive conditions on the coefficients compared with the restrictions that the classical formulation of the Second Law of Thermodynamics imposes. Non-monotonic creep curves were observed experimentally, however such a behavior is attributed to inertia either of the rod itself or of the rheometer. The paper is organized as follows. In Section we investigate conditions leading to constitutive equations containing complex derivatives of stress and strain that can be used in viscoelastic models of the wave equation. More precisely, we derive restrictions on parameters in constitutive equation of the form 1) under which the Laplace and Fourier transforms, as well as their inverses, of real-valued functions will remain real-valued. Also, in order to impose the validity of the Second Law of Thermodynamics we follow the procedure presented in [6], and obtain additional restrictions on model parameters in constitutive equations. As a result, we shall introduce a new form of the fractional derivative of complex order. Then, in Section 3, we treat the fractional Kelvin-Voigt complex order constitutive equation for viscoelastic body, find thermodynamical restrictions and prove sufficient conditions for its invertibility. Several numerical examples are presented in Section 4, as an illustration of creep and stress relaxation in viscoelastic materials. To the end of this section we recall basic definitions and results that will be used in our work. Fractional operators of complex order are introduced as follows see [19, 7]): For η C with < Re η < 1, definition of the left Riemann-Liouville fractional integral of an absolutely continuous function on [, T ], T > y AC[, T ])) coincides with the case of real η, i.e., I η 1 t yτ) t yt) := Γη) dτ, t [, T ], where Γ is the Euler gamma function. t τ) 1 η If η = iθ, θ R, then the latter integral diverges, and hence one introduces the fractional integration of imaginary order as I iθ t yt) := d dt I 1+iθ t yt) = 1 d Γ1 + iθ) dt t τ) iθ yτ) dτ, t [, T ].
3 However, in both cases the left Riemann-Liouville fractional derivative of order η C with Re η < 1 is given by D η t yt) := d dt I 1 η t yt) = 1 d Γ1 η) dt yτ) dτ, t [, T ]. t τ) η The basic tool for our study will be the Laplace and Fourier transforms. In order to have a good framework we will perform these transforms in S R), the space of tempered distributions. It is the dual space for the Schwartz space of rapidly decreasing functions SR). In particular, we are interested in the space S +R) whose elements are of the form y = P D)Y, where Y is a locally integrable polynomial bounded function on R that vanishes on, ), and P D) denotes a partial differential operator. The Fourier transform of y L 1 R) or y L R)) is defined as Fyω) = ŷω) = e iωx yx) dx, ω R. ) In the distributional setting, one has Fy, ϕ = y, Fϕ, y S R), ϕ SR), where Fϕ is defined by ). For y L 1 R) with yt) =, t <, and yt) Ae at, a, A >, the Laplace transform is given by Lys) = ỹs) = e st yt) dt, Re s > a. If y S +R) then a = since y is bounded by a polynomial). Then Ly is a holomorphic function in the half plane Re s > see e.g. [8]). Let Y s), Re s >, be a holomorphic function bounded by a polynomial in that domain. Then, for a suitable polynomial P, Y s)/p s) is integrable along the line Γ = a i, a+i ), and the inverse Laplace transform of Y is a tempered distribution yt) = P d dt )Y t), where Y t) = L 1 [Y ]t) = 1 πi Let y S +. Recall: [ ] d n F dx n y Γ Y s) P s) est ds. ω) = iω) n Fyω) ω R), [ ] d n L dt n y s) = s n Lys) Re s > ), n N, F[ D α x y]ω) = iω) α Fyω) ω R), L[ D α t y]s) = s α Lys) Re s > ), α C. Linear fractional constitutive equations with complex derivatives In what follows, we shall denote by α and β the orders of fractional derivatives. α will be assumed to be a real number, while β will be an element of C which is not real, i.e., β = A+iB, with B. Also, we shall assume that < α, A < 1..1 Real valued compatibility constraints Similarly as in the real case, see [], it is quite difficult to begin with a study of the most general case of 1). Therefore, in order to try to recognize the essence of the problem and find possibilities for overcoming it, we shall first concentrate to simpler forms of constitutive 3
4 equations that contain complex derivatives. Hooke law in the complex setting: Consider the following generalization of the σt) = b D β t εt), 3) where β C and b R. In order to find restrictions on parameters b and β in 3) which yield a physically acceptable constitutive equation, we shall verify the next two conditions: For real strain ε, the stress σ has to be real valued function of t. We call this a real valued compatibility requirement. Thermodynamical restrictions will result from the Second Law of Thermodynamics, and will be studied in the next section. Note that in the case of constitutive equations with only real-valued fractional derivatives, the real valued compatibility requirement always holds true, while the thermodynamical restrictions had to be investigated cf. []). Theorem.1 Let ε AC[, T ]) be real-valued, for all T >, < A < 1 and b. Then function σ defined by 3) is real valued if and only if β R. Proof. It follows from 3), with β = A + ib, and 1/Γ1 β) = h + ir, that Im σt) = b d dt Denoting by r h := tg φ, for h, we obtain Im σt) = εt τ)τ A ) r cosb ln τ) h sinb ln τ) dτ, t. bh cos φ d dt εt τ)τ A sinφ B ln τ) dτ, t. εt τ)τ A cosb ln τ) dτ. In the case h =, the imaginary part of σ reduces to br d dt If B = then r = and φ =, hence Im σ = and σ is a real-valued function. Next, suppose that B. But then one can find a subinterval t 1, t ) of, T ) where sinφ B ln τ), respectively cosb ln τ), is positive or negative), and choose ε AC[, T ]) which is compactly supported in t 1, t ) and strictly positive. This leads to a contradiction with the assumption Im σ =. The previous theorem implies that equations of form 3) with β C\R can not be a constitutive equation for a viscoelastic body. Next, consider the equation σt) = b 1 D β 1 t εt) + b D β t εt), t, 4) where b 1, b R, and β 1, β C\R, i.e., β k = A k + ib k and B k k = 1, ). Suppose again that ε AC[, T ]) is a real-valued function, for every T >. Remark. Note that dimension [ D β 1 t ε] is T β 1, where T is time unit. Therefore, 4) makes sense if [b 1 ] = T β 1 and [b ] = T β. Theorem.3 Function σ given by 4) is real-valued for all real-valued positive ε AC[, T ]) if and only if b 1 = b and β = β 1. 4
5 Proof. We continue with the notation of Theorem.1. Let t. Then σt) = b 1 d Γ1 β 1 ) dt εt τ)τ β 1 b d dτ + Γ1 β ) dt εt τ)τ β dτ Denote by h k + ir k := 1/Γ1 β k ), k = 1,. Then the imaginary part of the right hand side reads: Im σt) = d dt + d dt ) εt τ)τ A 1 b 1 r 1 cosb 1 ln τ) b 1 h 1 sinb 1 ln τ) dτ ) εt τ)τ A b r cosb ln τ) b h sinb ln τ) dτ 5) Using the identity Γ z) = Γz), it is straight forward to check that b 1 = b and β = β 1 imply that Im σ =, and hence σ is a real-valued function on [, T ], for every T >. Conversely, we want to find conditions on parameters which yield a real-valued function σ. Thus, we look at 5), with the change of variables p = ln τ, τ, t), t T, and solutions of the equation Set r k h k d dt d dt ln t e p εt e p ) b 1 r 1 e A1p cosb 1 p) b 1 h 1 e A1p sinb 1 p) ) +b r e Ap cosb p) b h e Ap sinb p) dp =, t [, T ]. 6) := tg φ k, k = 1,, for h 1, h. For h k = set φ k := π, k = 1,.) Then 6) gives ln t e p εt e p b1 h 1 ) e A1p sinφ 1 B 1 p) + b ) h e Ap sinφ B p) dp =. cos φ 1 cos φ Assume first that B 1 B, say B 1 > B. Since the basic period of sinφ 1 B 1 p) T = π/ B 1 ) is smaller than for sinφ B p), it follows that for every k N, k > k, where k depends on ln t, the function sinφ 1 B 1 p) changes its sign at least three times in the interval kπ, kπ +π/ B ). Thus, on that interval, there exist at least two intervals where sinφ 1 B 1 p) and sinφ B p) have the same sign, and two intervals where they have opposite signs. We conclude that there exists an interval [a, b] kπ, kπ + π/ B ), ln t), so that b 1 h 1 e A1p sinφ 1 B 1 p) + b h e Ap sinφ B p) >, p [a, b]. 7) cos φ 1 cos φ Choose δ > and k N so that {t e p ; t T/ δ, T/ + δ), p [a, b]} = T/ δ e b, T/ + δ e a ) = I, T ). Now, we choose a non-negative function ε AC[, T ]) with the following properties: supp ε I, so that the function p εt e p ), p [a, b], is strictly positive on some [a 1, b 1 ] a, b). This implies that for t T/ δ, T/ + δ), b a e p εt e p b1 h 1 ) e A1p sinφ 1 B 1 p) + b ) h e Ap sinφ B p) dp cos φ 1 cos φ is not a constant function. This is in contradiction with 6). 5
6 Therefore, in order to have 6), one must have B 1 = B. Moreover, arguing as above, one concludes that φ 1 B 1 p = φ B p must hold, for all p, ln t), t. Then we examine the equation b 1 h 1 e A 1p b h e A p =, or b 1 h 1 e A 1p + b h e A p =, p [kπ, kπ + π/ B 1 ]. Now in both cases, b 1 b > or b 1 b <, it is easy to conclude that A 1 = A, b 1 = b and B 1 = B have to be satisfied. This proves the theorem. Remark.4 i) Theorem.3 states that a real valued compatibility constraint for constitutive equations of form 4) holds if they contain complex fractional derivatives of strain, whose orders have to be complex conjugated numbers. Therefore, we may assume in the sequel, without loss of generality, that B >. ii) According to the above analysis, one can take arbitrary linear combination of pairs of complex conjugated fractional derivatives of strain. Moreover, one can also allow the same type of fractional derivatives of stress. Thus, one can consider the most general stress-strain constitutive equation with fractional derivatives of complex order: σt) + N i=1 c i D γ i t ) + D γ i t σt) = εt) + M j=1 b j D β j t + D β ) j t εt), where c i, b j R and γ i, β j C, i = 1,..., N, j = 1,..., M. iii) As a consequence one has that stress-strain relations can also contain arbitrary real order fractional derivatives, without any additional restrictions. This fact has already been known from previous work. iv) The same conclusions can also be obtained using a different approach. One can apply the result from [11, p. 93, Satz ], which tells that a function F is real-valued almost everywhere), if its Laplace transform is real-valued, for all real s in the half-plane of convergence right from some real x, in order to show that an admissible fractional constitutive equation 1) may be of complex order only if it contains pairs of complex conjugated fractional derivatives of stress and strain.. Thermodynamical restrictions In the analysis that follows we shall consider the isothermal processes only. For such processes the Second Law of Thermodynamics, i.e., the entropy of the system increases, is equivalent to the dissipativity condition T σt) εt) dt. 8) Equation 8) is one-dimensional version of the Second Principle of Thermodynamics for simple materials under isothermal conditions see [1, p. 83]). Also, in writing 8) we observe the fact that for t, ) the system is in virginal state. Note that T = in [1, p. 113]. Since we consider arbitrary T, condition 8) is stronger than the condition in [1]. In [1] stress is assumed as σt) = G εt) + 6 G u)εt u) du, 9)
7 where G is instanteneous elastic modulus, and G is relaxation function, while stain is taken in the form εs) = ε cosωs) + ε sinωs), s t, ω >, where ε is the amplitude of strain. These assumptions along with 8) and T = nπ ω Note that 1) is equivalent to F s [G t)]ω) = Im lead to G u) sinωu) du <. 1) ) F[G t)]ω) >, where F s [ft)]ω) = ft) sinωt) dt. Since constitutive equations containing fractional derivatives that we shall treat in this work do not obey 9), we follow the approach of [6]. Namely, in [6] it was assumed that periodic strain ε results in periodic stress σ after the end of transient regime. Positivity of dissipation work 8) for a cycle and each time instant within the cycle leads to Re Êω) and Im Êω), ω >, where Ê = ˆσˆε is the complex modulus obtained from the constitutive equation by applying the Fourier transform cf. [6, Eq. ), 1)]). Note that Re Ê and Im Ê are referred as storage and loss modulus respectively. We start with the constitutive equation σt) = b Dβ t εt), D β t := 1 D β t + D β ) t, t, 11) where we assume that b > and β = A + ib, < A < 1, B >. Note that 11) generalizes the Hooke law in the complex fractional framework. In the case β R, this new complex fractional operator Dβ t coincides with the usual left Riemann-Liouville fractional derivative. We apply the Fourier transform to 11): ˆσω) = b iω) β + iω) β)ˆεω), ω R. Then, define the complex modulus Ê such that ˆσω) = Êω) ˆεω), ω R, i.e., Êω) := b iω) β + iω) β ) = b ω A e Bπ e i Aπ +B ln ω) + e Bπ e i Aπ B ln ω)), ω R. Thermodynamical restrictions involve, for ω R +, Re Êω) = b ωa e Bπ Aπ ) cos + B ln ω + e Bπ Aπ )) cos B ln ω, 1) Im Êω) = b ωa e Bπ Aπ ) sin + B ln ω + e Bπ Aπ )) sin B ln ω. 13) But this is in contradiction with B >, since for ω >, 1) and 13) imply B =. In order not to confront the real valued compatibility requirement and the Second Law of Thermodynamics for 11), one may require that 1) and 13) hold for ω in some bounded interval instead of in all of R. Alternatively, as we shall do in the sequel, one can modify 11) by adding additional terms, in order to preserve the Second Law of Thermodynamics. Thus, we proceed by proposing the following constitutive equation σt) = a Dt α εt) + b Dβ t εt), t, 14) 7
8 where we assume that a, b >, α R, < α < 1, and β = A + ib, < A < 1, B >. Again we follow the procedure described above for deriving thermodynamical restrictions: ˆσω) = [aiω) α + biω) β + iω) β)]ˆεω), ω R. Consider the complex module ω R) Êω) = a ω α e i απ + b ω A e Bπ e i Aπ +B ln ω) + e Bπ e i Aπ B ln ω)) ; 15) Re Êω) = a ωα cos απ + b ωa e Bπ Aπ ) cos + B ln ω + e Bπ Aπ )) cos B ln ω,16) Im Êω) = a ωα sin απ + b ωa e Bπ Aπ ) sin + B ln ω + e Bπ Aπ )) sin B ln ω.17) We will investigate conditions Re Ê and Im Ê on R +. The assumption α > A leads to a contradiction, since for ω the sign of the second term in 16) determines the sign of Re Ê, and it can be negative. Thus, we must have α A. If α < A. then for ω, the second term in 16) could be negative. This together yields that the only possibility is A = α. The same conclusion is obtained if one considers Im Ê, ω >.) Therefore, with A = α, 16) and 17) become: with Re Êω) = a ωα cos απ + b ωα fω), ω >, 18) Im Êω) = a ωα sin απ + b ωα gω), ω >, 19) fω) := cos απ cosln ωb ) cosh Bπ απ + sin sinln ωb ) sinh Bπ, ω >, ) gω) := sin απ cosln ωb ) cosh Bπ απ cos sinln ωb ) sinh Bπ, ω >. 1) We further have that Re Êω) aωα cos απ + bωα min ω R+ fω), ω >, and the similar estimate for Im Ê, hence we now look for the minimums of functions f and g on R +. Using the substitution x = ln ω B we find that f x) = and g x) = at x f and x g so that tg x f = tg απ tgh Bπ and tg x g = ctg απ tgh Bπ. Solutions x f1, x f, x g1 and x g satisfy: x f1, π ), x f π, 3π ), and x g1 π, π), x g 3π απ απ Bπ, π), since tg, ctg, tgh >, and fx f ) = ± cos απ Bπ cosh 1 + tg απ Bπ ), tgh gx g ) = ± sin απ Bπ cosh 1 + ctg απ Bπ ). tgh Therefore, we have min x R fx) = fx f ) and min x R gx) = gx g1 ), so that 18) and 19) can be estimated as Re Êω) ωα cos απ a b cosh Bπ 1 + tg απ Bπ ) ) tgh, ω >, Im Êω) ωα sin απ a b cosh Bπ 1 + ctg απ Bπ ) ) tgh, ω >. 8
9 We obtain the thermodynamical restrictions for 14) by requiring Re Êω) and Im Êω), for ω R + : a b cosh Bπ 1 + ctg απ Bπ ), tgh if α, 1 ], ) and a b cosh Bπ 1 + tg απ Bπ ), tgh if [ 1 ) α, 1. 3) Notice that both restrictions further imply a b. Remark.5 i) Fix a and α. Inequalities ) and 3) impliy that as B increases, the constant b has to decrease, i.e., the contribution of complex fractional derivative of strain in the constitutive equation 14) is smaller if its imaginary part is larger. ii) Also, inequalities ) and 3) lead to the same restrictions on parameters a, b, α and B, since for α, 1 ] one has 1 α [ 1, 1), and the values of ) and 3) coincide. iii) Under the same conditions, constitutive equation 14) can be extended to σt) = εt) + a Dt α εt) + b Dβ t εt), t, which will be investigated in the next section. 3 Complex order fractional Kelvin-Voigt model Consider the constitutive equation involving the complex order fractional derivative σt) = εt) + a Dt α εt) + b Dβ t εt), t. 4) We assume a, b, E >, < α < 1, B >, β = α + ib, and σ and ε are real-valued functions. Note that 4) is a generalization of the model proposed in [6]. Namely, 4) agrees with the latter when β is real and positive. The inverse relation, i.e., ε as a function of σ is given in Theorem 3.. The Laplace transform of 4) is σs) = E 1 + a s α + b s β + s β) ) εs), Re s >, and hence 1 εs) = σs), Re s >. 5) 1 + a s α + b s β + s β) In order to determine ε from 5) we need to analyze zeros of ψs) = 1 + a s α + b s β + s β ), s C. 6) Note that if we put s = iω, ω R +, in 6), it becomes the complex modulus: ψiω) = 1 + Êω) = 1 + a iω)α + b iω) β + iω) β ), ω >, 7) where Ê is given in 15). Let s = re iϕ, r >, ϕ [, π]. Then with β = α + ib) and ψs) = 1 + ar α e iαϕ + br α e Bϕ e iln rb +αϕ) + e Bϕ e iln rb αϕ) ), Re ψs) = 1+ar α cosαϕ)+br α cosln r B ) cosαϕ) coshbϕ)+sinln r B ) sinαϕ) sinhbϕ) ), 8) Im ψs) = ar α sinαϕ) + br α cosln r B ) sinαϕ) coshbϕ) sinln r B ) cosαϕ) sinhbϕ) ). 9) 9
10 3.1 Thermodynamical restrictions In the case of 7), using 8) and 9) we obtain: Re ψiω) = 1 + Re Êω) 1 + aωα cos απ + bωα min x R fx), x = ln ωb, ω >, Im ψiω) = Im Êω) aωα sin απ + bωα min x R gx), x = ln ωb, ω >, where f and g are as in ) and 1). This leads to the same thermodynamical restrictions ) and 3), as in Section.. Therefore, from now on, we shall assume ) and 3) to hold true. Now we shall examine the zeros of ψ. 3. Zeros of ψ and solutions of 4) Theorem 3.1 Let ψ be the function given by 6). Then i) ψ has no zeros in the right complex half-plane Re s. ii) ψ has no zeros in C if the coefficients a, b, α and B satisfy a b coshbπ) 1 + tgαπ) tghbπ) ) [ 1, for α 4, 3 ] { 1 \, 4 } a b coshbπ) 1 + ctgαπ) tghbπ) ). for α, 1 ) 4 { 1 } 3 4, 1 ). 3) Proof. First, we notice that if s is a solution to ψs) =, then s the complex conjugate of s ) is also a solution, since ψ s) = 1 + a s α + b s β + s β) = ψs). Thus, we restrict our attention to the upper complex half-plane Im s, i.e., ϕ [, π]. Using 8) and 9) we have, with s = re iϕ, r >, ϕ [, π], and x = ln r B, where Re ψs) 1 + ar α cosαϕ) + br α min fx), x R 31) Im ψs) ar α sinαϕ) + br α min gx), x R 3) fx) = cosx) cosαϕ) coshbϕ) + sinx) sinαϕ) sinhbϕ), x R, 33) gx) = cosx) sinαϕ) coshbϕ) sinx) cosαϕ) sinhbϕ), x R. 34) The critical points x f and x g of f and g, respectively, satisfy tg x f = tgαϕ) tghbϕ) and tg x g = ctgαϕ) tghbϕ), 35) The proof of i) and ii) will be given by the argument principle. i) Consider ψ in the case Re s, Im s >. Choose a contour Γ = γ R1 γ R γ R3 γ R4, as it is shown in Fig. 1. γ R1 is parametrized by s = x, x ε, R) with ε and R, so that 8) and 9) yield Re ψx) = 1 + x α a + b cosln x β )) 1 + x α a b), Im ψx) =, 1
11 Im s R3 L L1 R4 R L4 L3 R1 Re s Figure 1: Contour Γ. since both ) and 3) imply a b. Moreover, we have lim x ψx) = 1 and lim x ψx) =. Along γ R one has s = Re iϕ, ϕ [, π ], with R. By 35) we have tg x f, so that and therefore 31) becomes min fx) = cosαϕ) coshbϕ) 1 + tgαϕ) tghbϕ)), x R Re ψs) 1 + R α cosαϕ)a b coshbϕ) 1 + tgαϕ) tghbϕ)) ). The previous inequality holds true, since for ϕ [, π ] we have that pϕ) = coshbϕ) 1 + tgαϕ) tghbϕ)) cosh Bπ 1 + tg απ Bπ ), tgh because of the fact that the function p monotonically increases on [, π ]. Moreover, by 8) and 9), we have Re ψs) and Im ψs) =, for ϕ =, R, Re ψs), for ϕ = π, R. The next segment is γ R3, which is parametrized by s = iω, ω [R, ε], with ε and R. Then 8) and 9) yield Re ψiω) = 1 + Re Êω) and Im ψiω) = Im Êω), ω ε, R), due to the thermodynamical requirements. Moreover, by 8) and 9), we have Re ψω) 1 and Im ψω), as ω, Re ψω) and Im ψω), as ω. The last part of the contour Γ is the arc γ R4, with s = εe iϕ, ϕ [, π ], with ε. Using the same arguments as for the contour γ R, we have Re ψs) 1 + ε α cosαϕ) a b coshbϕ) 1 + tgαϕ) tghbϕ)) ) 1. 36) 11
12 Also, by 9) and 36), we have Re ψs) 1 and Im ψs), as ε. We conclude that arg ψs) = so that, by the argument principle, there are no zeroes of ψ in the right complex half-plane Re s. ii) In order to discuss the zeros of ψ in the left complex half-plane, we use the contour Γ L = γ L1 γ L γ L3 γ L4, shown in Fig. 1. The contour γ L1 has the same parametrization as the contour γ R3, so the same conclusions as for γ R3 hold true. The parametrization of the contour γ L is s = Re iϕ, ϕ [ π, π], with R. Let us distinguish two cases. a1. αϕ, π ) Then sinαϕ) >, cosαϕ) >, so the critical points of f and g, see 33) and 34), given by 35), satisfy tg x f > and tg x g <. For the minimums of f and g in 31) and 3) we have min fx) = cosαϕ) coshbϕ) 1 + tgαϕ) tghbϕ)), x R min gx) = sinαϕ) coshbϕ) 1 + ctgαϕ) tghbϕ)), x R respectively, so that 31) and 3) become Re ψs) 1 + R α cosαϕ) a b coshbϕ) 1 + tgαϕ) tghbϕ) ) ), 37) Im ψs) R α sinαϕ) a b coshbϕ) 1 + ctgαϕ) tghbϕ) ) ). 38) Function H f ϕ) = coshbϕ) 1 + tgαϕ) tghbϕ)) is monotonically increasing for ϕ [ π, π], since αϕ, π ), thus Re ψs) 1 + R α cosαϕ) a b coshbπ) 1 + ) tgαπ) tghbπ)), if 3) is satisfied. Note that Re ψs) and Im ψs), for ϕ = π, R. b1. αϕ [ π, π) Then sinαϕ) >, cosαϕ), so the critical points of g, see 34), given by 35), satisfy tg x g. For the minimum of g in 34) we have min gx) = sinαϕ) coshbϕ) 1 + ctgαϕ) tghbϕ) ), x R so that 3) becomes Im ψs) R α sinαϕ) a b coshbϕ) 1 + ctgαϕ) tghbϕ) ) ). Function H g ϕ) = coshbϕ) 1 + ctgαϕ) tghbϕ)) is monotonically increasing for ϕ [ π, π], since αϕ [ π, π), and 3) implies Im ψs) R α sinαϕ) a b coshbπ) 1 + ctgαπ) tghbπ) ) ). Note that Im ψs), for ϕ = π, R. 1
13 Now we discuss possible situations for α, 1) and ϕ [ π, π]. If α, 1 ) then a1. holds so that Re ψs). If α [ 1, 1) then we distinguish two cases. For ϕ [ π, π α ) case a1. holds, so Re ψs). For ϕ [ π α, π) case b1. holds, and Im ψs). Parametrization of the contour γ L3 is s = xe iπ, x ε, R), with ε and R. Again, we have two cases. a. α, 1 ) Then, for x ε, R), using the same argumentation as in case a1, we have, by 37) and 38) with R = x), Re ψs) 1 and Im ψs), due to 3). Thus, looking at 37) and 38) with R = x), we conclude Re ψs) and Im ψs), for x, Re ψs) 1 and Im ψs), for x. b. α [ 1, 1) Then, for x ε, R), using the same argumentation as in case b1, we have by 38)) that Im ψs), and Im ψs), for x and Im ψs), for x. The parametrization of the contour γ L4 is s = εe iϕ, ϕ [ π, π], with ε. From 8) and 9), for sufficiently small ε, we have [ π ] Re ψs) 1 and Im ψs), for ε, ϕ, π. Summing up all results from cases a1, a, b1 and b, we obtain the following: For α, 1 ), Re ψs), for s Γ L, which implies that arg ψs) =. Therefore, using the argument principle, we conclude that in this case ψ has no zeros in the left complex half-plane. If α [ 1, 1), then for s γ L1 and s {z γ L arg z π α }, we have Re ψs), while for s {z γ L arg z > π α } and s γ L3, we have Im ψs). For s γ L4 we again have Re ψs). Hence, we conclude that arg ψs) =, and therefore, using the argument principle, neither in case α [ 1, 1) function ψ has zeros in the left complex half-plane. This completes the proof. Rewrite 5) as εs) = Ks) σs), 1 Ks) =, Re s >. 39) 1 + as α + bs β + s β) Theorem 3. Let ε be given by 39). Then εt) = Kt) σt), t. 4) 13
14 Moreover, if 3) holds, then Kt) = K I t) = 1 πi e qt ] 1 + q α e [a iαπ + be i ln qb e Bπ + e i ln qb e Bπ ) ) e qt ] 1 + q α e [a iαπ + be i ln qb e Bπ + e i ln qb e Bπ ) dq. 41) If condition 3) is violated, then ψ has at most finite number of zeros in the left complex half-plane, and K = K I or K = K I + K R, where K R t) = ψs i )= i=1,,...,n Res Ks)e st, s i ) + Res Ks)e ) st, s i ), 4) with K given by 39). Proof. The first part is clear. We invert now K, given in 39), by the use of the Cauchy residues theorem Ks)e st ds = πi Res Ks)e st, s) 43) Γ ψ s)= and the contour Γ = Γ 1 Γ Γ r Γ 3 Γ 4 γ shown in Fig.. Figure : Contour Γ. If condition 3) is satisfied, then, by Theorem 3.1, the residues equal zero. One can show that the integrals over the contours Γ 1, Γ r and Γ 4 tend to zero when R and r. The 14
15 remaining integrals give: lim R, r lim R, r lim R, r Ks)e st ds = Γ 3 Ks)e st ds = γ Ks)e st ds = πik I t), e qt ] dq, 1 + q α e [a iαπ + be i ln qb e Bπ + e i ln qb e Bπ ) e qt ] dq, 1 + q α e [a iαπ + be i ln qb e Bπ + e i ln qb e Bπ ) which, by the Cauchy residues theorem 43), leads to 41). If condition 3) is violated, then, by Theorem 3.1, the denominator of K either has no zeros in the complex plane, and so K = K I, or it has zeros in the left complex halfplane, which comes in pairs with complex conjugates. We show now that ψ which is the dominator of K) can have at most finite number of zeros for Im s. Rewrite ψs) = as a + b s ib + s ib ) = s α. If s = re iϕ, ϕ [ π, π], then we have a + b r ib e ϕb + r ib e ϕb ) = r α e αϕi. When r the right hand side tends to zero, while the left hand side tends to a. As r we see that the left hand side is bounded, while the right hand side is not bounded. Thus, the zeros in the left half-plane of function ψ, if exist, have to be bounded both from above and below. In that case we have K = K I + K R, where K R is given by 4). Remark 3.3 In the cases of standard and fractional linear solid model creep compliance and relaxation modulus are exponential functions and derivative of a one-parameter Mittag-Leffler function. Also, these models have different relaxation and retardation times as well as the limiting values of the material functions. For more complicated models, as our is, functions representing creep compliance and relaxation modulus are possibly not so well known special functions. We note that other generalizations of standard linear viscoelastic solids, e.g. the one presented in [6], also do not tackle this question. 4 Numerical verifications Here we present several examples of the proposed constitutive equation. We shall treat stress relaxation, creep and periodic loading cases. 4.1 Stress relaxation experiment We take 4) with εt) = Ht), H is the Heaviside function, and regularize it as H ε t) = 1 exp t/k), k. In order to determine σ we calculate σt) = H ε t) + a D α t H ε t) + b Dβ t H εt), t, 44) with σ) =, using the expansion formula see [3, 4]), D γ t yt) yt) t γ AN, γ) N p=1 C p 1 γ) V p 1y)t) t p+γ, 45) 15
16 where and AN, γ) = Inserting 45) into 44) we obtain { [ a AN, α) σt) 1 + where ΓN γ) γγ1 γ)γγ)n!, C p 1γ) = Γp + γ) Γ1 γ)γγ)p 1)!, V 1) p 1 y)t) = tp 1 yt), V p 1 y)) =, p = 1,, 3,... 46) N p=1 t α { Cp 1 α) t p+α + AN, α + ib) + b t α+ib [ Cp 1 α + ib) t p+α+ib AN, α ib) + t α ib + C p 1α ib) t p+α ib )]} H ε t) ]} V p 1 H ε )t), 47) V 1) p 1 H ε)t) = t p 1 H ε t), p = 1,, 3,... 48) We will compare 47) with the stress σ obtained by 44) and by the definition of fractional derivative 11): σt) = H ε t) + a d 1 dt Γ1 α) +b d [ 1 dt Γ1 α ib) H ε τ) dτ t τ) α H ε τ) dτ t τ) α+ib + 1 Γ1 α + ib) ] H ε τ) dτ t τ) α ib. 49) In Fig. 3 we show results obtained by determining σ from 49) for small times and different values of B. In the same figure we show, by dots, the values of σ, in several points, obtained by using 47), 48) for k =.1, N = 1. As could be seen from Fig. 3, the results obtained from 49) and 47), 48) agree well. The stress relaxation curves are shown in Fig. 4, for the same set of parameters and for larger times. As could be seen, regardless of the values of B, we have lim t σt) = 1. Note that in all cases of B, the restriction which follows from the dissipation inequality is satisfied. Σ t 5 4 B.4 B.6 3 B.8 1 B t Figure 3: Stress relaxation curves for α =.4 and B {.4,.6,.8,.99}, a =.8, b =.1, t [,.5]. 16
17 Σ t B.4 B.6 B.8 B t Figure 4: Stress relaxation curves for α =.4 and B {.4,.6,.8,.99}, a =.8, b =.1, t [, 1]. 4. Creep experiment Suppose that σt) = Ht), i.e., By 45), we obtain or εt) εt) Ht) = 1 + a Dt α + b Dβ t )εt), t. 5) Ht) + { N Cp 1 α) p=1 t p+α 1 + a AN,α) t α Ht)t α + { N Cp 1 α) p=1 t p By using 51) in 46) we obtain + C p 1α+iB) t p+α+ib AN,α+iB) t α+ib + b t α + a AN, α) + b + C p 1α+iB) t p+ib AN,α+iB) t ib { V 1) Ht)tα + N Cp 1 α) p 1 ε)t) p=1 tp 1 t p t α + a AN, α) + b V p 1 ε)) =, p = 1,, 3,... + C p 1α ib) t p+α ib }V p 1 ε)t) + AN,α ib) t α ib ), + C p 1α ib) t p ib }V p 1 ε)t) + C p 1α+iB) t p+ib AN,α+iB) t ib Equation 5) may also be solved by contour integration + AN,α ib) t ib ). 51) + C p 1α ib) t p ib }V p 1 ε)t) + AN,α ib) t ib ), εt) = Kt) Ht), t, 5) where K is given by 41), see 4). Finally, the values of ε, at discrete points, could be determined directly from εs) = 1 Ks), s Re s >, see 39), by the use of Post inversion formula, see [9]. Thus, 1) n ] n+1 dn 1 εt) = lim [s n n! ds n s1 + as α + bs β + s β)), t. 53) s= n t 17
18 ε t t Figure 5: Creep curve for α =.4 and B =.4, a =.8, b =.1, t [, 1]. In Fig. 5, 6 and 7 we show ε for several values of parameters determined from 5). In Fig. 5, for specified values of t we present values of ε, determined from 51), with N = 7, denoted by dots, as well as the values of ε, determined by 53), with n = 5, denoted by squares. As could be seen the agreement of results determined by different methods is significant. From Fig. 6 and 7 one sees that, regardless of the value of B, creep curves tend to ε = 1. In Fig. 6 the creep curves are monotonically increasing, while in Fig. 7 creep curves has oscillatory character, characteristic for the case when the mass of the rod is not neglected. Note that in all cases of B, the restriction determined by the dissipation inequality is satisfied. ε t B B.4 B t Figure 6: Creep curves for α =.4 and B {.,.4,.6}, a =.8, b =.1, t [, ]. 5 Conclusion In this work, we proposed a new constitutive equation with fractional derivatives of complex order for viscoelastic body of the Kelvin-Voigt type. The use of fractional derivatives of complex order, together with restrictions following from the Second Law of Thermodynamics, represent the main novelty of our work. Our results can be summarized as follows. 1. In order to obtain real stress for given real strain, we used two fractional derivatives of complex order that are complex conjugated numbers, see Theorem.3. 18
19 ε t B.7 B.8 B.9 B t Figure 7: Creep curves for α =.4 and B {.7,.8,.9,.99}, a =.8, b =.1, t [, ].. The restrictions that follow from the Second Law of Thermodynamics for isothermal process implied that the constitutive equation must additionally contain a fractional derivative of real order. Thus, the simplest constitutive equation that gives real stress for real strain and satisfies the dissipativity condition is given by 4). 3. We provided a complete analysis of solvability of the complex order fractional Kelvin- Voigt model given by 4) see Theorems 3.1 and 3.). 4. We studied stress relaxation and creep problems through equation 4). An increase of B implied that the stress relaxation decreases more rapidly to the limiting value of stress, i.e., lim t σt) = We presented numerical experiments when the dissipation inequality is satisfied. The creep experiment showed that the increase in the imaginary part of the complex derivative B changes the character of creep curve from monotonic to oscillatory form, see Fig. 6 and 7. However, the creep curves never cross the value equal to 1. The creep curve resembles the form of creep curve when either the mass of the rod, see [5, p. 14], or inertia of the rheometer, see [18], is taken into account. Recently, [9] presented experimental results of creep curves for some biological materials that exhibit nonmonotonic creep curve. Thus, our model 4) is applicable for such a case. 6. The parameters in the proposed model 4) could be determined from experimental data using for instance the least square method. 7. Our further study will be directed to the problems of vibration and wave propagation of a rod with finite mass and constitutive equation 4), along the lines presented in [5], see also [15]. Acknowledgements We would like to thank Marko Janev for several helpful discussions on the subject. This work is supported by Projects 1745 and 1744 of the Serbian Ministry of Science, and of the Provincial Secretariat for Science. 19
20 References [1] Amendola, G., Fabrizio, M., Golden, J. M. Thermodynamics of Materials with Memory. Springer, New York, 1. [] Atanacković, T.M., Konjik, S., Oparnica, Lj., Zorica, D. Thermodynamical restrictions and wave propagation for a class of fractional order viscoelastic rods. Abstr. Appl. Anal., 11: pp), 11. [3] Atanacković, T.M., Stanković, B. An expansion formula for fractional derivatives and its applications. Fract. Calc. Appl. Anal., 73): , 4. [4] Atanacković, T.M., Stanković, B. On a numerical scheme for solving differential equations of fractional order. Mech. Res. Commun., 35:49 438, 8. [5] Atanacković, T. M., Pilipović, S., Stanković, B., Zorica, D. Fractional Calculus with Applications in Mechanics: Vibrations and Diffusion Processes. Wiley-ISTE, London, 14. [6] Bagley, R. L., Torvik, P. J. On the fractional calculus model of viscoelastic behavior. J. Rheology, 31): , [7] Caputo, M., Mainardi, F. Linear models of dissipation in anelastic solids. Riv. Nuovo Cimento Ser. II), 1: , [8] Caputo, M., Mainardi, F. A new dissipation model based on memory mechanism. Pure Appl. Geophys., 911): , [9] Cohen, A. M. Numerical Methods for Laplace Transform Inversion. Springer, New York, 1. [1] Day, W. A. The Thermodynamics of Simple Materials with Fading Memory. Springer- Verlag, Berlin, 197. [11] Doetsch, G. Handbuch der Laplace-Transformationen I. Birkhäuser, Basel, 195. [1] Fabrizio, M., Lazzari, B. Stability and Second Law of Thermodynamics in dual-phase-lag heat conduction. Int. J. Heat Mass Transfer, 74: , 1 4. [13] Fabrizio, M., Morro, A. Mathematical Problems in Linear Viscoelasticity. SIAM Studies in Applied Mathematics, Philadelphia, 199. [14] Gonsovski, V. L., Rossikhin, Yu. A. Stress waves in a viscoelastic medium with a singular hereditary kernel. J. Appl. Mech. Tech. Phys., 14: , [15] Hanyga, A. Long-range asymptotics of a step signal propagating in a hereditary viscoelastic medium. Q. J. Mechanics Appl. Math., 6):85 98, 7. [16] Hanyga, A. Multi-dimensional solutions of space-time-fractional diffusion equations. Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci., 45818):49 45,. [17] Hanyga, A. Multidimensional solutions of time-fractional diffusion-wave equations. Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci., 458): ,.
21 [18] Jaishankar, A., McKinley, G. H. Power-law rheology in the bulk and at the interface: quasi-properties and fractional constitutive equations. Proc. R. Soc. Lond., Ser. A, Math. Phys. Eng. Sci., 469:184, 1. [19] Love, E. R. Fractional derivatives of imaginary order. J. London Math. Soc., -3):41 59, [] Machado, J. A. T. Optimal controllers with complex order derivatives. J. Optim. Theory. Appl., 1561): 1, 13. [1] Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity. Imperial College Press, London, 1. [] Mainardi, F., Pagnini, G., Gorenflo, R. Some aspects of fractional diffusion equations of single and distributed order. Appl. Math. Comput., 1871):95 35, 7. [3] Makris, N. Complex-parameter Kelvin model for elastic foundations. Earthq. Eng. Struct. Dyn., 33):51 64, [4] Makris, N., Constantinou, M. Models of viscoelasticity with complex-order derivatives. J. Eng. Mech., 1197): , [5] Rossikhin, Yu.A., Shitikova, M.V. A new method for solving dynamic problems of fractional derivative viscoelasticity. Int. J. Eng. Sci., 39): , 1. [6] Rossikhin, Yu.A., Shitikova, M.V. Analysis of rheological equations involving more than one fractional parameters by the use of the simplest mechanical systems based on these equations. Mech. Time-Depend. Mater., 5): , 1. [7] Samko, S.G., Kilbas, A.A., Marichev, O.I. Fractional Integrals and Derivatives - Theory and Applications. Gordon and Breach Science Publishers, Amsterdam, [8] Vladimirov V. S. Equations of Mathematical Physics. Mir Publishers, Moscow, [9] Zingales, M. A mechanical description of anomalous time evolution: fractional hereditariness, heat transport and mass diffusion. Mechanics through Mathematical Modelling Book of abstracts, Novi Sad 15. 1
NEW RHEOLOGICAL PROBLEMS INVOLVING GENERAL FRACTIONAL DERIVATIVES WITH NONSINGULAR POWER-LAW KERNELS
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Volume 19, Number 1/218, pp. 45 52 NEW RHEOLOGICAL PROBLEMS INVOLVING GENERAL FRACTIONAL DERIVATIVES WITH NONSINGULAR
More informationNEW GENERAL FRACTIONAL-ORDER RHEOLOGICAL MODELS WITH KERNELS OF MITTAG-LEFFLER FUNCTIONS
Romanian Reports in Physics 69, 118 217 NEW GENERAL FRACTIONAL-ORDER RHEOLOGICAL MODELS WITH KERNELS OF MITTAG-LEFFLER FUNCTIONS XIAO-JUN YANG 1,2 1 State Key Laboratory for Geomechanics and Deep Underground
More informationFractional acoustic wave equations from mechanical and. Sverre Holm
Fractional acoustic wave equations from mechanical and thermal constitutive equations Sverre Holm Outline 2 Perspectives Non fractional wave equations Low and high frequency regions Fractional viscoelastic
More informationFRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY
FRACTIONAL CALCULUS AND WAVES IN LINEAR VISCOELASTICITY by Francesco Mainardi (University of Bologna, Italy) E-mail: Francesco.Mainardi@bo.infn.it Imperial College Press, London 2, pp. xx+ 347. ISBN: 978--8486-329-4-8486-329-
More informationMulti-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
International Journal of Difference Equations ISSN 0973-6069, Volume 0, Number, pp. 9 06 205 http://campus.mst.edu/ijde Multi-Term Linear Fractional Nabla Difference Equations with Constant Coefficients
More informationSettlement analysis of viscoelastic foundation under vertical line load using a fractional Kelvin-Voigt model
Geomechanics and Engineering, Vol. 4, No. 1 (01) 67-78 67 Technical Note Settlement analysis of viscoelastic foundation under vertical line load using a fractional Kelvin-Voigt model Hong-Hu Zhu* 1, Lin-Chao
More informationCONTROL OF THERMAL STRESSES IN AXISSYMMETRIC PROBLEMS OF FRACTIONAL THERMOELASTICITY FOR AN INFINITE CYLINDRICAL DOMAIN
THERMAL SCIENCE: Year 17, Vol. 1, No. 1A, pp. 19-8 19 CONTROL OF THERMAL STRESSES IN AXISSYMMETRIC PROBLEMS OF FRACTIONAL THERMOELASTICITY FOR AN INFINITE CYLINDRICAL DOMAIN by Yuriy POVSTENKO a, Derya
More informationFRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS
FRACTIONAL FOURIER TRANSFORM AND FRACTIONAL DIFFUSION-WAVE EQUATIONS L. Boyadjiev*, B. Al-Saqabi** Department of Mathematics, Faculty of Science, Kuwait University *E-mail: boyadjievl@yahoo.com **E-mail:
More informationChapter 2: Complex numbers
Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We
More informationNewtonian Mechanics. Chapter Classical space-time
Chapter 1 Newtonian Mechanics In these notes classical mechanics will be viewed as a mathematical model for the description of physical systems consisting of a certain (generally finite) number of particles
More informationMathematics 104 Fall Term 2006 Solutions to Final Exam. sin(ln t) dt = e x sin(x) dx.
Mathematics 14 Fall Term 26 Solutions to Final Exam 1. Evaluate sin(ln t) dt. Solution. We first make the substitution t = e x, for which dt = e x. This gives sin(ln t) dt = e x sin(x). To evaluate the
More informationTime fractional Schrödinger equation
Time fractional Schrödinger equation Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 The Schrödinger equation is considered with the first order time
More informationComplex Numbers. The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition,
Complex Numbers Complex Algebra The set of complex numbers can be defined as the set of pairs of real numbers, {(x, y)}, with two operations: (i) addition, and (ii) complex multiplication, (x 1, y 1 )
More informationPartial Differential Equations for Engineering Math 312, Fall 2012
Partial Differential Equations for Engineering Math 312, Fall 2012 Jens Lorenz July 17, 2012 Contents Department of Mathematics and Statistics, UNM, Albuquerque, NM 87131 1 Second Order ODEs with Constant
More informationTHE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY
THE FUNDAMENTAL SOLUTION OF THE DISTRIBUTED ORDER FRACTIONAL WAVE EQUATION IN ONE SPACE DIMENSION IS A PROBABILITY DENSITY RUDOLF GORENFLO Free University of Berlin Germany Email:gorenflo@mi.fu.berlin.de
More informationIMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES
Dynamic Systems and Applications ( 383-394 IMPROVEMENTS OF COMPOSITION RULE FOR THE CANAVATI FRACTIONAL DERIVATIVES AND APPLICATIONS TO OPIAL-TYPE INEQUALITIES M ANDRIĆ, J PEČARIĆ, AND I PERIĆ Faculty
More informationMahmoud M. El-Borai a, Abou-Zaid H. El-Banna b, Walid H. Ahmed c a Department of Mathematics, faculty of science, Alexandria university, Alexandria.
International Journal of Basic & Applied Sciences IJBAS-IJENS Vol:13 No:01 52 On Some Fractional-Integro Partial Differential Equations Mahmoud M. El-Borai a, Abou-Zaid H. El-Banna b, Walid H. Ahmed c
More informationarxiv: v1 [math.na] 5 Mar 2017
Analytical and numerical treatment of the heat conduction equation obtained via time-fractional distributed-order heat conduction law Velibor Želi, Dušan Zorica arxiv:173.3v1 [math.na] 5 Mar 17 September
More informationA truncation regularization method for a time fractional diffusion equation with an in-homogeneous source
ITM Web of Conferences, 7 18) ICM 18 https://doi.org/1.151/itmconf/187 A truncation regularization method for a time fractional diffusion equation with an in-homogeneous source Luu Vu Cam Hoan 1,,, Ho
More informationMath 417 Midterm Exam Solutions Friday, July 9, 2010
Math 417 Midterm Exam Solutions Friday, July 9, 010 Solve any 4 of Problems 1 6 and 1 of Problems 7 8. Write your solutions in the booklet provided. If you attempt more than 5 problems, you must clearly
More informationMTH3101 Spring 2017 HW Assignment 4: Sec. 26: #6,7; Sec. 33: #5,7; Sec. 38: #8; Sec. 40: #2 The due date for this assignment is 2/23/17.
MTH0 Spring 07 HW Assignment : Sec. 6: #6,7; Sec. : #5,7; Sec. 8: #8; Sec. 0: # The due date for this assignment is //7. Sec. 6: #6. Use results in Sec. to verify that the function g z = ln r + iθ r >
More informationDETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL DIFFUSION EQUATION
Journal of Fractional Calculus and Applications, Vol. 6(1) Jan. 2015, pp. 83-90. ISSN: 2090-5858. http://fcag-egypt.com/journals/jfca/ DETERMINATION OF AN UNKNOWN SOURCE TERM IN A SPACE-TIME FRACTIONAL
More informationHere are brief notes about topics covered in class on complex numbers, focusing on what is not covered in the textbook.
Phys374, Spring 2008, Prof. Ted Jacobson Department of Physics, University of Maryland Complex numbers version 5/21/08 Here are brief notes about topics covered in class on complex numbers, focusing on
More informationMAT389 Fall 2016, Problem Set 11
MAT389 Fall 216, Problem Set 11 Improper integrals 11.1 In each of the following cases, establish the convergence of the given integral and calculate its value. i) x 2 x 2 + 1) 2 ii) x x 2 + 1)x 2 + 2x
More informationA Note on the Differential Equations with Distributional Coefficients
MATEMATIKA, 24, Jilid 2, Bil. 2, hlm. 115 124 c Jabatan Matematik, UTM. A Note on the Differential Equations with Distributional Coefficients Adem Kilicman Department of Mathematics, Institute for Mathematical
More informationSIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE. Yuriy Povstenko. Abstract
SIGNALING PROBLEM FOR TIME-FRACTIONAL DIFFUSION-WAVE EQUATION IN A HALF-PLANE Yuriy Povstenko Abstract The time-fractional diffusion-wave equation is considered in a half-plane. The Caputo fractional derivative
More informationFORCED VIBRATIONS WITH FRACTIONAL TYPE OF DISSIPATION
FORCED VIBRATIONS WITH FRACTIONAL TYPE OF DISSIPATION D. T. SPASIC & N. C. CHARALAMBAKIS University of Novi Sad, Institute of Applied Mechanics, Yugoslavia Aristotle University of Thessaloniki, Department
More informationMATH 452. SAMPLE 3 SOLUTIONS May 3, (10 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic.
MATH 45 SAMPLE 3 SOLUTIONS May 3, 06. (0 pts) Let f(x + iy) = u(x, y) + iv(x, y) be an analytic function. Show that u(x, y) is harmonic. Because f is holomorphic, u and v satisfy the Cauchy-Riemann equations:
More informationTime domain analysis of viscoelastic models
Time domain analysis of viscoelastic models L. Rouleau 1, J.-F. Deü 1 1 LMSSC, Conservatoire national des arts et métiers, 2D6R10, 292 rue Saint-Martin, 75141 Paris cedex 03, France e-mail: lucie.rouleau@lecnam.net
More informationSubmitted Version to CAMWA, September 30, 2009 THE LAPLACE TRANSFORM ON ISOLATED TIME SCALES
Submitted Version to CAMWA, September 30, 2009 THE LAPLACE TRANSFORM ON ISOLATED TIME SCALES MARTIN BOHNER AND GUSEIN SH. GUSEINOV Missouri University of Science and Technology, Department of Mathematics
More informationON FRACTIONAL RELAXATION
Fractals, Vol. 11, Supplementary Issue (February 2003) 251 257 c World Scientific Publishing Company ON FRACTIONAL RELAXATION R. HILFER ICA-1, Universität Stuttgart Pfaffenwaldring 27, 70569 Stuttgart,
More informationMSE 383, Unit 3-3. Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept.
Dynamic Mechanical Behavior MSE 383, Unit 3-3 Joshua U. Otaigbe Iowa State University Materials Science & Engineering Dept. Scope Why DMA & TTS? DMA Dynamic Mechanical Behavior (DMA) Superposition Principles
More informationTHE UNIVERSITY OF WESTERN ONTARIO. Applied Mathematics 375a Instructor: Matt Davison. Final Examination December 14, :00 12:00 a.m.
THE UNIVERSITY OF WESTERN ONTARIO London Ontario Applied Mathematics 375a Instructor: Matt Davison Final Examination December 4, 22 9: 2: a.m. 3 HOURS Name: Stu. #: Notes: ) There are 8 question worth
More informationON THE C-LAGUERRE FUNCTIONS
ON THE C-LAGUERRE FUNCTIONS M. Ishteva, L. Boyadjiev 2 (Submitted by... on... ) MATHEMATIQUES Fonctions Specialles This announcement refers to a fractional extension of the classical Laguerre polynomials.
More informationFunctions of a Complex Variable (S1) Lecture 11. VII. Integral Transforms. Integral transforms from application of complex calculus
Functions of a Complex Variable (S1) Lecture 11 VII. Integral Transforms An introduction to Fourier and Laplace transformations Integral transforms from application of complex calculus Properties of Fourier
More informationA review: The Laplacian and the d Alembertian. j=1
Chapter One A review: The Laplacian and the d Alembertian 1.1 THE LAPLACIAN One of the main goals of this course is to understand well the solution of wave equation both in Euclidean space and on manifolds
More informationSMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS. Kai Diethelm. Abstract
SMOOTHNESS PROPERTIES OF SOLUTIONS OF CAPUTO- TYPE FRACTIONAL DIFFERENTIAL EQUATIONS Kai Diethelm Abstract Dedicated to Prof. Michele Caputo on the occasion of his 8th birthday We consider ordinary fractional
More informationAnswers to Problem Set Number MIT (Fall 2005).
Answers to Problem Set Number 6. 18.305 MIT (Fall 2005). D. Margetis and R. Rosales (MIT, Math. Dept., Cambridge, MA 02139). December 12, 2005. Course TA: Nikos Savva, MIT, Dept. of Mathematics, Cambridge,
More informationFractional Schrödinger Wave Equation and Fractional Uncertainty Principle
Int. J. Contemp. Math. Sciences, Vol., 007, no. 9, 943-950 Fractional Schrödinger Wave Equation and Fractional Uncertainty Principle Muhammad Bhatti Department of Physics and Geology University of Texas
More informationLinear Fractionally Damped Oscillator. Mark Naber a) Department of Mathematics. Monroe County Community College. Monroe, Michigan,
Linear Fractionally Damped Oscillator Mark Naber a) Department of Mathematics Monroe County Community College Monroe, Michigan, 48161-9746 In this paper the linearly damped oscillator equation is considered
More informationOn the local existence for an active scalar equation in critical regularity setting
On the local existence for an active scalar equation in critical regularity setting Walter Rusin Department of Mathematics, Oklahoma State University, Stillwater, OK 7478 Fei Wang Department of Mathematics,
More informationSecond Midterm Exam Name: Practice Problems March 10, 2015
Math 160 1. Treibergs Second Midterm Exam Name: Practice Problems March 10, 015 1. Determine the singular points of the function and state why the function is analytic everywhere else: z 1 fz) = z + 1)z
More information1 Res z k+1 (z c), 0 =
32. COMPLEX ANALYSIS FOR APPLICATIONS Mock Final examination. (Monday June 7..am 2.pm) You may consult your handwritten notes, the book by Gamelin, and the solutions and handouts provided during the Quarter.
More informationWe denote the space of distributions on Ω by D ( Ω) 2.
Sep. 1 0, 008 Distributions Distributions are generalized functions. Some familiarity with the theory of distributions helps understanding of various function spaces which play important roles in the study
More informationThe Internal Friction and the Relaxation Time Spectrum of Ferroelectric Ceramic PZT Type
Vol. 114 008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology The Internal Friction and the Relaxation Time Spectrum of Ferroelectric Ceramic PZT Type J. Bartkowska
More informationConformal maps. Lent 2019 COMPLEX METHODS G. Taylor. A star means optional and not necessarily harder.
Lent 29 COMPLEX METHODS G. Taylor A star means optional and not necessarily harder. Conformal maps. (i) Let f(z) = az + b, with ad bc. Where in C is f conformal? cz + d (ii) Let f(z) = z +. What are the
More informationIndia
italian journal of pure and applied mathematics n. 36 216 (819 826) 819 ANALYTIC SOLUTION FOR RLC CIRCUIT OF NON-INTGR ORDR Jignesh P. Chauhan Department of Applied Mathematics & Humanities S.V. National
More informationThe solutions of time and space conformable fractional heat equations with conformable Fourier transform
Acta Univ. Sapientiae, Mathematica, 7, 2 (25) 3 4 DOI:.55/ausm-25-9 The solutions of time and space conformable fractional heat equations with conformable Fourier transform Yücel Çenesiz Department of
More informationSOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER
Dynamic Systems and Applications 2 (2) 7-24 SOME RESULTS FOR BOUNDARY VALUE PROBLEM OF AN INTEGRO DIFFERENTIAL EQUATIONS WITH FRACTIONAL ORDER P. KARTHIKEYAN Department of Mathematics, KSR College of Arts
More informationA sufficient condition for the existence of the Fourier transform of f : R C is. f(t) dt <. f(t) = 0 otherwise. dt =
Fourier transform Definition.. Let f : R C. F [ft)] = ˆf : R C defined by The Fourier transform of f is the function F [ft)]ω) = ˆfω) := ft)e iωt dt. The inverse Fourier transform of f is the function
More informationOn boundary value problems for fractional integro-differential equations in Banach spaces
Malaya J. Mat. 3425 54 553 On boundary value problems for fractional integro-differential equations in Banach spaces Sabri T. M. Thabet a, and Machindra B. Dhakne b a,b Department of Mathematics, Dr. Babasaheb
More informationFunctions of a Complex Variable and Integral Transforms
Functions of a Complex Variable and Integral Transforms Department of Mathematics Zhou Lingjun Textbook Functions of Complex Analysis with Applications to Engineering and Science, 3rd Edition. A. D. Snider
More informationComplex Numbers and the Complex Exponential
Complex Numbers and the Complex Exponential φ (2+i) i 2 θ φ 2+i θ 1 2 1. Complex numbers The equation x 2 + 1 0 has no solutions, because for any real number x the square x 2 is nonnegative, and so x 2
More informationMath Final Exam.
Math 106 - Final Exam. This is a closed book exam. No calculators are allowed. The exam consists of 8 questions worth 100 points. Good luck! Name: Acknowledgment and acceptance of honor code: Signature:
More information1 Discussion on multi-valued functions
Week 3 notes, Math 7651 1 Discussion on multi-valued functions Log function : Note that if z is written in its polar representation: z = r e iθ, where r = z and θ = arg z, then log z log r + i θ + 2inπ
More informationA generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives
A generalized Gronwall inequality and its application to fractional differential equations with Hadamard derivatives Deliang Qian Ziqing Gong Changpin Li Department of Mathematics, Shanghai University,
More informationChapter 31. The Laplace Transform The Laplace Transform. The Laplace transform of the function f(t) is defined. e st f(t) dt, L[f(t)] =
Chapter 3 The Laplace Transform 3. The Laplace Transform The Laplace transform of the function f(t) is defined L[f(t)] = e st f(t) dt, for all values of s for which the integral exists. The Laplace transform
More informationlim when the limit on the right exists, the improper integral is said to converge to that limit.
hapter 7 Applications of esidues - evaluation of definite and improper integrals occurring in real analysis and applied math - finding inverse Laplace transform by the methods of summing residues. 6. Evaluation
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationDIfferential equations of fractional order have been the
Multistage Telescoping Decomposition Method for Solving Fractional Differential Equations Abdelkader Bouhassoun Abstract The application of telescoping decomposition method, developed for ordinary differential
More informationLinear viscoelastic behavior
Harvard-MIT Division of Health Sciences and Technology HST.523J: Cell-Matrix Mechanics Prof. Ioannis Yannas Linear viscoelastic behavior 1. The constitutive equation depends on load history. 2. Diagnostic
More informationSOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD
SOLUTION OF FRACTIONAL INTEGRO-DIFFERENTIAL EQUATIONS BY ADOMIAN DECOMPOSITION METHOD R. C. Mittal 1 and Ruchi Nigam 2 1 Department of Mathematics, I.I.T. Roorkee, Roorkee, India-247667. Email: rcmmmfma@iitr.ernet.in
More informationSynopsis of Complex Analysis. Ryan D. Reece
Synopsis of Complex Analysis Ryan D. Reece December 7, 2006 Chapter Complex Numbers. The Parts of a Complex Number A complex number, z, is an ordered pair of real numbers similar to the points in the real
More informationMath 185 Fall 2015, Sample Final Exam Solutions
Math 185 Fall 2015, Sample Final Exam Solutions Nikhil Srivastava December 12, 2015 1. True or false: (a) If f is analytic in the annulus A = {z : 1 < z < 2} then there exist functions g and h such that
More informationMath 341: Probability Seventeenth Lecture (11/10/09)
Math 341: Probability Seventeenth Lecture (11/10/09) Steven J Miller Williams College Steven.J.Miller@williams.edu http://www.williams.edu/go/math/sjmiller/ public html/341/ Bronfman Science Center Williams
More informationMathematical Methods for Physics and Engineering
Mathematical Methods for Physics and Engineering Lecture notes for PDEs Sergei V. Shabanov Department of Mathematics, University of Florida, Gainesville, FL 32611 USA CHAPTER 1 The integration theory
More informationMath 411, Complex Analysis Definitions, Formulas and Theorems Winter y = sinα
Math 411, Complex Analysis Definitions, Formulas and Theorems Winter 014 Trigonometric Functions of Special Angles α, degrees α, radians sin α cos α tan α 0 0 0 1 0 30 π 6 45 π 4 1 3 1 3 1 y = sinα π 90,
More informationNUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX
Palestine Journal of Mathematics Vol. 6(2) (217), 515 523 Palestine Polytechnic University-PPU 217 NUMERICAL SOLUTION OF FRACTIONAL ORDER DIFFERENTIAL EQUATIONS USING HAAR WAVELET OPERATIONAL MATRIX Raghvendra
More information(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS
(2m)-TH MEAN BEHAVIOR OF SOLUTIONS OF STOCHASTIC DIFFERENTIAL EQUATIONS UNDER PARAMETRIC PERTURBATIONS Svetlana Janković and Miljana Jovanović Faculty of Science, Department of Mathematics, University
More informationDIFFERENTIAL EQUATIONS
DIFFERENTIAL EQUATIONS Chapter 1 Introduction and Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationJASSON VINDAS AND RICARDO ESTRADA
A QUICK DISTRIBUTIONAL WAY TO THE PRIME NUMBER THEOREM JASSON VINDAS AND RICARDO ESTRADA Abstract. We use distribution theory (generalized functions) to show the prime number theorem. No tauberian results
More informationApplied Mathematics Letters
Applied Mathematics Letters 24 (211) 219 223 Contents lists available at ScienceDirect Applied Mathematics Letters journal homepage: www.elsevier.com/locate/aml Laplace transform and fractional differential
More informationare harmonic functions so by superposition
J. Rauch Applied Complex Analysis The Dirichlet Problem Abstract. We solve, by simple formula, the Dirichlet Problem in a half space with step function boundary data. Uniqueness is proved by complex variable
More informationin Bounded Domains Ariane Trescases CMLA, ENS Cachan
CMLA, ENS Cachan Joint work with Yan GUO, Chanwoo KIM and Daniela TONON International Conference on Nonlinear Analysis: Boundary Phenomena for Evolutionnary PDE Academia Sinica December 21, 214 Outline
More informationExercises for Part 1
MATH200 Complex Analysis. Exercises for Part Exercises for Part The following exercises are provided for you to revise complex numbers. Exercise. Write the following expressions in the form x + iy, x,y
More informationdx n a 1(x) dy
HIGHER ORDER DIFFERENTIAL EQUATIONS Theory of linear equations Initial-value and boundary-value problem nth-order initial value problem is Solve: a n (x) dn y dx n + a n 1(x) dn 1 y dx n 1 +... + a 1(x)
More informationThis work has been submitted to ChesterRep the University of Chester s online research repository.
This work has been submitted to ChesterRep the University of Chester s online research repository http://chesterrep.openrepository.com Author(s): Kai Diethelm; Neville J Ford Title: Volterra integral equations
More informationSPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS
ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 43, Number 3, 2013 SPECTRAL PROPERTIES OF THE SIMPLE LAYER POTENTIAL TYPE OPERATORS MILUTIN R. OSTANIĆ ABSTRACT. We establish the exact asymptotical behavior
More informationCHAPTER 3 ELEMENTARY FUNCTIONS 28. THE EXPONENTIAL FUNCTION. Definition: The exponential function: The exponential function e z by writing
CHAPTER 3 ELEMENTARY FUNCTIONS We consider here various elementary functions studied in calculus and define corresponding functions of a complex variable. To be specific, we define analytic functions of
More informationThe local fractional Hilbert transform in fractal space
The local fractional ilbert transform in fractal space Guang-Sheng Chen Department of Computer Engineering, Guangxi Modern Vocational Technology College, echi,guangxi, 547000, P.. China E-mail address:
More informationOn the Finite Caputo and Finite Riesz Derivatives
EJTP 3, No. 1 (006) 81 95 Electronic Journal of Theoretical Physics On the Finite Caputo and Finite Riesz Derivatives A. M. A. El-Sayed 1 and M. Gaber 1 Faculty of Science University of Alexandria, Egypt
More informationEXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1. Yong Zhou. Abstract
EXISTENCE AND UNIQUENESS OF SOLUTIONS FOR A SYSTEM OF FRACTIONAL DIFFERENTIAL EQUATIONS 1 Yong Zhou Abstract In this paper, the initial value problem is discussed for a system of fractional differential
More informationarxiv: v2 [math.ca] 8 Nov 2014
JOURNAL OF THE AMERICAN MATHEMATICAL SOCIETY Volume 00, Number 0, Pages 000 000 S 0894-0347(XX)0000-0 A NEW FRACTIONAL DERIVATIVE WITH CLASSICAL PROPERTIES arxiv:1410.6535v2 [math.ca] 8 Nov 2014 UDITA
More informationEXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM
1 EXISTENCE AND UNIQUENESS FOR A THREE DIMENSIONAL MODEL OF FERROMAGNETISM V. BERTI and M. FABRIZIO Dipartimento di Matematica, Università degli Studi di Bologna, P.zza di Porta S. Donato 5, I-4126, Bologna,
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationThe Laplace Transform
C H A P T E R 6 The Laplace Transform Many practical engineering problems involve mechanical or electrical systems acted on by discontinuous or impulsive forcing terms. For such problems the methods described
More informationApplication of fractional order theory of thermoelasticity to a 1D problem for a cylindrical cavity
Arch. Mech., 66, 4, pp. 257 267, Warszawa 2014 Application of fractional order theory of thermoelasticity to a 1D problem for a cylindrical cavity W. E. RASLAN Department of Mathematics and Engineering
More informationIn the present work the diffraction of plane electromagnetic waves by an impedance loaded parallel plate waveguide formed by a twopart
Progress In Electromagnetics Research, PIER 6, 93 31, 6 A HYBRID METHOD FOR THE SOLUTION OF PLANE WAVE DIFFRACTION BY AN IMPEDANCE LOADED PARALLEL PLATE WAVEGUIDE G. Çınar Gebze Institute of Technology
More informationOrdinary Differential Equation Theory
Part I Ordinary Differential Equation Theory 1 Introductory Theory An n th order ODE for y = y(t) has the form Usually it can be written F (t, y, y,.., y (n) ) = y (n) = f(t, y, y,.., y (n 1) ) (Implicit
More informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationTwo special equations: Bessel s and Legendre s equations. p Fourier-Bessel and Fourier-Legendre series. p
LECTURE 1 Table of Contents Two special equations: Bessel s and Legendre s equations. p. 259-268. Fourier-Bessel and Fourier-Legendre series. p. 453-460. Boundary value problems in other coordinate system.
More informationNotes for Expansions/Series and Differential Equations
Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated
More information(x 1, y 1 ) = (x 2, y 2 ) if and only if x 1 = x 2 and y 1 = y 2.
1. Complex numbers A complex number z is defined as an ordered pair z = (x, y), where x and y are a pair of real numbers. In usual notation, we write z = x + iy, where i is a symbol. The operations of
More informationTravelling waves. Chapter 8. 1 Introduction
Chapter 8 Travelling waves 1 Introduction One of the cornerstones in the study of both linear and nonlinear PDEs is the wave propagation. A wave is a recognizable signal which is transferred from one part
More informationStieltjes Transformation as the Iterated Laplace Transformation
International Journal of Mathematical Analysis Vol. 11, 2017, no. 17, 833-838 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ijma.2017.7796 Stieltjes Transformation as the Iterated Laplace Transformation
More informationBulletin T.CXXXIII de l Académie serbe des sciences et des arts 2006 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 31
Bulletin T.CXXXIII de l Académie serbe des sciences et des arts 2006 Classe des Sciences mathématiques et naturelles Sciences mathématiques, No 31 GENERALIZED SOLUTIONS TO SINGULAR INITIAL-BOUNDARY HYPERBOLIC
More informationSome Fun with Divergent Series
Some Fun with Divergent Series 1. Preliminary Results We begin by examining the (divergent) infinite series S 1 = 1 + 2 + 3 + 4 + 5 + 6 + = k=1 k S 2 = 1 2 + 2 2 + 3 2 + 4 2 + 5 2 + 6 2 + = k=1 k 2 (i)
More informationSolutions to practice problems for the final
Solutions to practice problems for the final Holomorphicity, Cauchy-Riemann equations, and Cauchy-Goursat theorem 1. (a) Show that there is a holomorphic function on Ω = {z z > 2} whose derivative is z
More informationIntroduction to Differential Equations
Chapter 1 Introduction to Differential Equations 1.1 Basic Terminology Most of the phenomena studied in the sciences and engineering involve processes that change with time. For example, it is well known
More informationThe first order quasi-linear PDEs
Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.
More information