Analytical and numerical methods for the time-fractional nonlinear diusion
|
|
- Maximilian Rogers
- 5 years ago
- Views:
Transcription
1 Analytical and numerical methods for the time-fractional nonlinear diusion Šukasz Pªociniczak Faculty of Pure and Applied Mathematics, Wroclaw University of Science and Technology BCAM, Bilbao, Spain 1 / 26
2 Introduction This talk is about a 1-D diusion of moisture in building materials (such as siliceous brick). By u(x, t) we denote the moisture concentration at x at time t. We consider the following initial-boundary conditions: u(0, t) = C, u(x, 0) = 0. Self-similarity - a characteristic feature of diusion in our experiment. Moisture concentration u(x, t) can be drawn on a single curve [1]: u(x, t) = U(η), η = x/ t, for U(0) = C i U( ) = 0. 2 / 26
3 Nature is tricky (and thus very interesting)! As it turns out, the diusion not always behaves as we are used to. In a number of experiments (ex. [2-4]) the so-called Boltzmann scaling η = x/t 1/2 is not observed. A more appropriate and accurate is the anomalous diusion scaling (Figure from [2]) u(x, t) = U(η), η = x/t α/2, 0 < α < 2. 3 / 26
4 How to describe this mathematically? The classical diusion equation does not work: for our initial-boundary conditions it possesses the x/ t scaling. In [2,4] the following modication of the constitutive equation has been proposed ( ) 1 u α 1 q = D(u). x Result: very complicated equations (doubly degenerate PME) and average tting accuracy. As it turned out, a more appropriate is to model this phenomenon by an equation with fractional derivative (see [5-7]) α u t α = ( D(u) u ). x x We obtain the sought scaling x/t α/2 with very small tting errors. 4 / 26
5 Model We consider a 1D time-fractional porous medium equation α u t α = ( u m u ), 0 < α < 1, m 1, x x with initial-boundary conditions u(0, t) = 1, u(x, 0) = 0. We seek for a self-similar slution u(x, t) = U(η), where η = x/t α/2. We obtain an ordinary integro-dierential equation [ (U m U ) = (1 α) α 2 η d ] F α U, F α := I 0,1 α, m 1, dη 2 α with U(0) = 1 and U( ) = 0, where the integral operator is of the Erdelyi-Kober type 5 / 26 I a,b c U(η) := 1 Γ(b) 1 0 (1 z) b 1 z a U(ηz 1 c )dz.
6 Physical derivation Time-fractional porous medium equation is not an ad-hoc choice. We assume Water can be trapped in some regions of porous medium for prolonged periods of time non-local conservation law. Consitutive relation Darcy's Law. Waiting times have a power-type weight emergence of Caputo fractional derivative. Zero initial condition equivalence of Riemann-Liouville and Caputo derivatives. The relevant paper with full derivation: Š.Pªociniczak, Analytical studies of a time-fractional porous medium equation. Derivation, approximation and applications, Communications in Nonlinear Science and Numerical Simulation 24 (13) (2015), / 26
7 Existence and uniqueness First, we will show that the equation [ (U m U ) = (1 α) α 2 η d dη ] F α U, m 1, with U(0) = 1 and U( ) = 0 possesses a unique compactly-supported solution. This is a free-boundary problem. The relevant paper: Š.Pªociniczak, M. witaªa, Existence and uniqueness results for a time-fractional nonlinear diusion equation, under review. 7 / 26
8 Existence and uniqueness We assume the compact support, i.e. there exists some η such that U(η) = 0 for η η (physically motivated). The key-idea is to use the transformation (borrowed from the works of Okrasi«ski [9]) U(η) = ( m(η ) 2) 1 m y(z), z = 1 η η, which turns our free-boundary problem into an initial-value one [ m(y m y ) = (1 α) + α 2 (1 z) d ] G α y, 0 < α < 1, 0 z 1, dz which y(0) = 0 (G α is related to F α ). What about the second condition? Theorem The solution y = y(z) of the above problem has to satisfy ( α ) 2 α Γ ( ) 2 α m y(z) 2 Γ ( z 2 α m ) 1 α + 2 α m (2 α) ( ) for z 0 +. m 1 8 / 26
9 Existence and uniqueness Next, transform the problem into an equivalent integral equation ( z y(z) = 0 ) 1 m+1 G(z, t)f α y(t)dt =: T (y)(z), for some kernel G. Nonlinear operator T is monotone and homogeneous of degree 1 m+1. Theorem ([10], P.J. Bushell, 1976) Suppose that T : K K is a monotone increasing mapping which is homogeneous of degree p with 0 < p < 1. If there exists ϕ K 0 such that γ 1 ϕ T (ϕ) γ 2 ϕ, for some positive constants γ 1,2 then T has a unique xed-point y K 0. 9 / 26
10 Existence and uniqueness Theorem Let ϕ(z) := z 2 α m, then γ 1 ϕ T (ϕ) γ 2 ϕ, with γ 1 = ( ( α 2 ( ( α 2 ) 1 α Γ( 2 α m ) Γ(2 α+ 2 α ) 2 α Γ(1+ 2 α m ) Γ(2 α+ 2 α m ) m ) 2 α+m(3 α) 1 2 α ) 1 m+1 1, 0 < α 1 1 m+1 ; ) 1 m+1, 1 1 m+1 < α 1, γ 2 = Γ(3 α) 1 m+1. Existence and uniqueness follows from Bushell's Theorem. 10 / 26
11 Approximate solutions Now, we will nd some approximate solutions of our problem. The strategy is to substitute the E-K operator (nonlocal) by some local approximation. This can be accurate for α 1. The relevant papers: Š.Pªociniczak, Analytical studies of a time-fractional porous medium equation. Derivation, approximation and applications, Communications in Nonlinear Science and Numerical Simulation 24 (13) (2015), Š.Pªociniczak, Approximation of the Erdelyi-Kober fractional operator with application to the time-fractional porous medium equation, SIAM Journal of Applied Mathematics 74(4) (2014), Š.Pªociniczak, H.Okrasi«ska, Approximate self-similar solutions to a nonlinear diusion equation with time-fractional derivative, Physica D 261 (2013), / 26
12 Approximate solutions Theorem For analytic U and a > 1, b > 0, c > 0 we have the following representation I a,b c U(η) = k=0 λ k U (k) (η) ηk k!, where λ k = k ( k ) j=0 j ( 1) k j Γ(a+ j +1) c Γ(a+b+ j +1). c Moreover, we have an asymptotic expansion when k λ k ( 1) k c ( ) c(a+1) Γ (c(a + 1)) 1 Γ (c(a + 1)). Γ(b) k The series converges very fast, especially for η close to 0. We can use it and obtain a series of useful approximations U i. 12 / 26
13 Numerical results Wetting front position Cumulative moisture t t Figure: On the left: front position η (t) (solid line) and its approximation η 3 (t) = η 3 t α/2 (dashed line). On the right: cumulative moisture (solid line) and approximations I 1 (dashed line) i I 2 (dot-dashed line). Here, α = 0.95 and m = / 26
14 Numerical results 3 m 3 D UHhL h Figure: Fitting U3 with experimental data from [2]. On the left: a self-similar pro le; an the right: time evolution. Here α = 0.855, C = 0.71 m3 /m3, m = 6.98, D0 = 5.36 mm/s /
15 Inverse problems In applications it is possible to measure U, but how to infer about the material properties of the medium? The relevant quantity is the diusivity D(U). Our (inverse) problem is to nd D(U) provided we know U. The relevant paper: Š.Pªociniczak, Diusivity identication in a nonlinear time-fractional diusion equation, Fractional Calculus and Applied Analysis 19(4) (2016), / 26
16 Inverse Problems Denition A mathematical problem (for example a dierential or algebraic equation) is called a well-posed problem if it has all of the following properties it has a solution, it has only one solution, it is continuous with respect to the perturbations in the initial data (stability). If a problem is not well-posed it is called ill-posed. Inverse problems usually belong to this class. 16 / 26
17 Diusivity identication Existence and uniqueness (up to a constant D s ) - easy integration D(U(η)) = 1 [ ( U D s U (0) + 1 α ) η F α U(z)dz α ] (η) 2 2 ηf αu(η). What about stability and computational cost? For the latter we have an approximation (based on our previous results) 0 D(U(η)) = 1 U (η) [D su (0) ( ) 1 + Γ(1 α) + α η ] α U(z)dz 2Γ(2 α) 2Γ(2 α) ηu(η). As we can see, this formula does not require double integration (much cheaper in applications) / 26
18 Main results Denition For η 0 > 0 the supremum norm of a function U : [0, ] R is dened by the formula U,η0 := sup U(η). η [0,η 0] Denition A function U : [0, ] [0, ] is called admissible if it is bounded along with its rst derivative, nonincreasing and vanishing at innity. Theorem Let U be admissible and x η 0 > 0 such that E η0 := sup η [0,η0] η/u (η) is nite. We then have D(U) D(U) E η0 (2A(α) U + B(α)η 0 U ),,η0 where A(α) and B(α) are O(1 α) as α 1 and explicitly known. 18 / 26
19 Main results. Regularization Dierentiation is not a stable operation, thus we have to regularize our formula for D. Let U δ be the moisture with a measurement noise, i.e. U U δ δ. Introduce a regularization strategy for U (ex. nite dierence) U (U δ ) h R(h, δ), for which there exists h 0 = h 0 (δ) such as R(h 0 (δ), δ) 0 as δ 0. Moreover, ( U δ) is a stable operator for each h > 0. h 19 / 26
20 Main results. Regularization Theorem Let D h be a family of regularization operators for the diusivity and x η 0 > 0. Then for U X such that there exists an ɛ > 0 with ɛ inf η [0,η0] U (η) < we have D(U) Dh (U δ ),η0 1 ( [D s 1 + R(h, δ) + U ɛ ɛ + η ( 0 Γ(2 α) δ + δ + U R(h, δ) ɛ ) R(h, δ) )]. Theorem Let the assumptions of the above Theorem be fullled. If the regularization strategy is such that R(h 0 (δ), δ) = O(δ µ ) for δ 0 then D(U) D h (U δ ),η0 = O(δ µ ) as δ / 26
21 Numerical methods This is a work in progress but we have some preliminary results. Discretization of the E-K operator n 1 Ic a,b U(η) = w n,i U(η i ) + R n (η). i=0 We have several schemes for choosing w n,i along with exact asymptotic behaviour of R n as n. Our results can be applied to any integro-dierential equations with E-K operator. The relevant work: Š.Pªociniczak, S. Sobieszek, Numerical schemes for integro-dierential equations with Erdélyi-Kober fractional operator, Numerical Algorithms 76(1) (2016), pp / 26
22 Numerical methods The second problem is a numerical solution of ( z ) 1 m+1 y(z) = G(z, t)f α y(t)dt, 0 which is not easy since the nonlinearity is not Lipschitzian. The above equation appears in the above integral formulation of time-fractional porous medium equation. Results: Convergence of a family of nite-dierence schemes. An explicit construction of a convergent method. Main assumption: kernel has to be positively bounded from below (now working how to remove this: for anomalous diusion we have G(z, u) C(z u) 1 α ). The relevant work: H. Okrasi«ska-Pªociniczak, Š.Pªociniczak, Finite dierence method for Volterra equation with a power-type nonlinearity, arxiv: [math.na] 22 / 26
23 Numerical methods The crucial step of the convergence proof is done with aid of the following Gronwall-Bellman's Lemma. Lemma Let {e n }, n = 1, 2,... be a sequence of positive numbers satisfying ( ) e n 1 n 1 A e i + B, n 2, (1) n i=1 for A > 0. If we dene M := max{e 1, B} 1 Aζ(1 A), 0 < A < 1 1, A 1, (2) then 23 / 26 e n M. (3) n1 A
24 Numerical methods: a sketch of the proof The method's error can be estimated as follows (m + 1)ξ m n e n (m + 1)ξ m n y(nh) y n = y(nh) m+1 yn m+1 hwd n 1 e i + δ(h) Some work is required to estimate ξ n in terms of y(nh). Then, a comparison theorem for integral equations helps us to obtain a recurrence inequality only for e i. Gronwall-Bellman's Lemma gives the assertion. In the above way we can prove the following estimate i=1 e n const. h 1 C m max{h p, δ(h)h 1 } as h 0 + with nh x n, where e 1 = O(h p ) as h 0 + with p > 0 and C is a (known) constant. 24 / 26
25 Numerical example Consider the following equation y(x) m+1 = which can be readily solved with x y(x) = 0 y(t)dt, x [0, 1], ( ) 1 m m + 1 x m. Our theorem gives us the order of convergence equal to at least 1 1 m (which is not large, though). To calculate the numerical order of convergence we use Aitken's method based on Richardson's extrapolation. As a quadrature we use the midpoint's method. 25 / 26 m Order Est. order
26 References 1. L. Pel, K. Kopinga, G. Bertram and G. Lang, Water absorption in a red-clay brick observed by NMR scanning, J. Phys. D.: Appl. Phys. 28 (1995) Abd El-Ghany, El Abd and J.J. Milczarek, Neutron radiography study of water absorption in porous building materials: anomalous diusion analysis, J. Phys. D.: Appl. Phys. 37 (2004) S.C. Taylor, W.D. Ho, M.A. Wilson and K.M. Green, Anomalous water transport properties of Portland and blended cement-based materials, J. Mater. Sci. Lett. 18 (1999) M.Kuntz and P.Lavallee, Experimental evidence and theoretical analysis of anomalous diusion during water inltration in porous building materials, J. Phys. D: Appl. Phys. 34 (2001), E. Gerolymatou, I.Vardoulakis and R.Hilfer, Modelling inltration by means of a nonlinear fractional diusion model, J. Phys. D: Appl. Phys. 39 (2006), D.N. Gerasimov, V.A. Kondratieva and O.A. Sinkevich, An anomalous non-self-similar inltration and fractional diusion equation, Physica D 239 (2010) Y. Pachepsky, D. Timlin and W. Rawls, Generalized Richard's equation to simulate water transport in unsaturated soils, J. Hydrol. (2003) J.R. King, Approximate solutions to a nonlinear diusion equation, Journal of Engineering Mathematics 22: (1988) 9. W. Okrasi«ski, On nontrivial solutions to some nonlinear ordinary dierential equations, Journal of mathematical analysis and applications 190(2) 1995, P.J. Bushell, On a class of Volterra and Fredholm non-linear integral equations, 26 / 26 Mathematical Proceedings of the Cambridge Philosophical Society 79(02) 1976,
Maximum principle for the fractional diusion equations and its applications
Maximum principle for the fractional diusion equations and its applications Yuri Luchko Department of Mathematics, Physics, and Chemistry Beuth Technical University of Applied Sciences Berlin Berlin, Germany
More informationARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (1996), 13 { 27. ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION
ARCHIVUM MATHEMATICUM (BRNO) Tomus 32 (996), 3 { 27 ON THE OSCILLATION OF AN mth ORDER PERTURBED NONLINEAR DIFFERENCE EQUATION P. J. Y. Wong and R. P. Agarwal Abstract. We oer sucient conditions for the
More informationPointwise convergence rate for nonlinear conservation. Eitan Tadmor and Tao Tang
Pointwise convergence rate for nonlinear conservation laws Eitan Tadmor and Tao Tang Abstract. We introduce a new method to obtain pointwise error estimates for vanishing viscosity and nite dierence approximations
More informationRearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton D
Rearrangements and polar factorisation of countably degenerate functions G.R. Burton, School of Mathematical Sciences, University of Bath, Claverton Down, Bath BA2 7AY, U.K. R.J. Douglas, Isaac Newton
More informationAlgorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations
dmitry.lyakhov@kaust.edu.sa Symbolic Computations May 4, 2018 1 / 25 Algorithmic Lie Symmetry Analysis and Group Classication for Ordinary Dierential Equations Dmitry A. Lyakhov 1 1 Computational Sciences
More informationNew trends in applied mathematics: fractional calculus
New trends in applied mathematics: fractional calculus Šukasz Pªociniczak 1 Faculty of Pure and Applied Mathematics, Wroclaw University of Science and Technology, Wroclaw, Poland 27-29.11.2018 USC 1 Support
More informationRecent result on porous medium equations with nonlocal pressure
Recent result on porous medium equations with nonlocal pressure Diana Stan Basque Center of Applied Mathematics joint work with Félix del Teso and Juan Luis Vázquez November 2016 4 th workshop on Fractional
More informationOn a WKB-theoretic approach to. Ablowitz-Segur's connection problem for. the second Painleve equation. Yoshitsugu TAKEI
On a WKB-theoretic approach to Ablowitz-Segur's connection problem for the second Painleve equation Yoshitsugu TAKEI Research Institute for Mathematical Sciences Kyoto University Kyoto, 606-8502, JAPAN
More informationTHE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOZZI AND M. PUGH.
THE LUBRICATION APPROXIMATION FOR THIN VISCOUS FILMS: REGULARITY AND LONG TIME BEHAVIOR OF WEAK SOLUTIONS A.L. BERTOI AND M. PUGH April 1994 Abstract. We consider the fourth order degenerate diusion equation
More informationVolterra Integral Equations of the First Kind with Jump Discontinuous Kernels
Volterra Integral Equations of the First Kind with Jump Discontinuous Kernels Denis Sidorov Energy Systems Institute, Russian Academy of Sciences e-mail: contact.dns@gmail.com INV Follow-up Meeting Isaac
More informationON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND STEFAN DIEHL
ON SCALAR CONSERVATION LAWS WITH POINT SOURCE AND DISCONTINUOUS FLUX FUNCTION STEFAN DIEHL Abstract. The conservation law studied is @u(x;t) + @ F u(x; t); x @t @x = s(t)(x), where u is a f(u); x > 0 concentration,
More informationFrom Fractional Brownian Motion to Multifractional Brownian Motion
From Fractional Brownian Motion to Multifractional Brownian Motion Antoine Ayache USTL (Lille) Antoine.Ayache@math.univ-lille1.fr Cassino December 2010 A.Ayache (USTL) From FBM to MBM Cassino December
More informationCongurations of periodic orbits for equations with delayed positive feedback
Congurations of periodic orbits for equations with delayed positive feedback Dedicated to Professor Tibor Krisztin on the occasion of his 60th birthday Gabriella Vas 1 MTA-SZTE Analysis and Stochastics
More informationSpurious Chaotic Solutions of Dierential. Equations. Sigitas Keras. September Department of Applied Mathematics and Theoretical Physics
UNIVERSITY OF CAMBRIDGE Numerical Analysis Reports Spurious Chaotic Solutions of Dierential Equations Sigitas Keras DAMTP 994/NA6 September 994 Department of Applied Mathematics and Theoretical Physics
More informationTools from Lebesgue integration
Tools from Lebesgue integration E.P. van den Ban Fall 2005 Introduction In these notes we describe some of the basic tools from the theory of Lebesgue integration. Definitions and results will be given
More informationThe Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations
Symmetry in Nonlinear Mathematical Physics 1997, V.1, 185 192. The Higher Dimensional Bateman Equation and Painlevé Analysis of Nonintegrable Wave Equations Norbert EULER, Ove LINDBLOM, Marianna EULER
More informationGeneration of undular bores and solitary wave trains in fully nonlinear shallow water theory
Generation of undular bores and solitary wave trains in fully nonlinear shallow water theory Gennady El 1, Roger Grimshaw 1 and Noel Smyth 2 1 Loughborough University, UK, 2 University of Edinburgh, UK
More informationApplication of Taylor-Padé technique for obtaining approximate solution for system of linear Fredholm integro-differential equations
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 1 (June 2017), pp. 392 404 Applications and Applied Mathematics: An International Journal (AAM) Application of Taylor-Padé
More informationfor all subintervals I J. If the same is true for the dyadic subintervals I D J only, we will write ϕ BMO d (J). In fact, the following is true
3 ohn Nirenberg inequality, Part I A function ϕ L () belongs to the space BMO() if sup ϕ(s) ϕ I I I < for all subintervals I If the same is true for the dyadic subintervals I D only, we will write ϕ BMO
More informationViscosity approximation methods for the implicit midpoint rule of asymptotically nonexpansive mappings in Hilbert spaces
Available online at www.tjnsa.com J. Nonlinear Sci. Appl. 9 016, 4478 4488 Research Article Viscosity approximation methods for the implicit midpoint rule of asymptotically nonexpansive mappings in Hilbert
More informationA slow transient diusion in a drifted stable potential
A slow transient diusion in a drifted stable potential Arvind Singh Université Paris VI Abstract We consider a diusion process X in a random potential V of the form V x = S x δx, where δ is a positive
More informationConverse Lyapunov Functions for Inclusions 2 Basic denitions Given a set A, A stands for the closure of A, A stands for the interior set of A, coa sta
A smooth Lyapunov function from a class-kl estimate involving two positive semidenite functions Andrew R. Teel y ECE Dept. University of California Santa Barbara, CA 93106 teel@ece.ucsb.edu Laurent Praly
More informationSYSTEMS OF ODES. mg sin ( (x)) dx 2 =
SYSTEMS OF ODES Consider the pendulum shown below. Assume the rod is of neglible mass, that the pendulum is of mass m, and that the rod is of length `. Assume the pendulum moves in the plane shown, and
More informationComputation Of Asymptotic Distribution. For Semiparametric GMM Estimators. Hidehiko Ichimura. Graduate School of Public Policy
Computation Of Asymptotic Distribution For Semiparametric GMM Estimators Hidehiko Ichimura Graduate School of Public Policy and Graduate School of Economics University of Tokyo A Conference in honor of
More informationEntrance Exam, Real Analysis September 1, 2017 Solve exactly 6 out of the 8 problems
September, 27 Solve exactly 6 out of the 8 problems. Prove by denition (in ɛ δ language) that f(x) = + x 2 is uniformly continuous in (, ). Is f(x) uniformly continuous in (, )? Prove your conclusion.
More informationAsymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics.
Asymptotic Distributions for the Nelson-Aalen and Kaplan-Meier estimators and for test statistics. Dragi Anevski Mathematical Sciences und University November 25, 21 1 Asymptotic distributions for statistical
More informationContinuous Functions on Metric Spaces
Continuous Functions on Metric Spaces Math 201A, Fall 2016 1 Continuous functions Definition 1. Let (X, d X ) and (Y, d Y ) be metric spaces. A function f : X Y is continuous at a X if for every ɛ > 0
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT. A. Khani
Acta Universitatis Apulensis ISSN: 1582-5329 No 38/214 pp 11-22 ADOMIAN-TAU OPERATIONAL METHOD FOR SOLVING NON-LINEAR FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS WITH PADE APPROXIMANT A Khani Abstract In this
More informationThe Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive Mappings in Hilbert Spaces
Applied Mathematical Sciences, Vol. 11, 2017, no. 12, 549-560 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ams.2017.718 The Generalized Viscosity Implicit Rules of Asymptotically Nonexpansive
More informationA MULTIGRID ALGORITHM FOR. Richard E. Ewing and Jian Shen. Institute for Scientic Computation. Texas A&M University. College Station, Texas SUMMARY
A MULTIGRID ALGORITHM FOR THE CELL-CENTERED FINITE DIFFERENCE SCHEME Richard E. Ewing and Jian Shen Institute for Scientic Computation Texas A&M University College Station, Texas SUMMARY In this article,
More informationSPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING. splitting, which now seems to be the main cause of the stochastic behavior in
SPLITTING OF SEPARATRICES FOR (FAST) QUASIPERIODIC FORCING A. DELSHAMS, V. GELFREICH, A. JORBA AND T.M. SEARA At the end of the last century, H. Poincare [7] discovered the phenomenon of separatrices splitting,
More informationOn the singularities of non-linear ODEs
On the singularities of non-linear ODEs Galina Filipuk Institute of Mathematics University of Warsaw G.Filipuk@mimuw.edu.pl Collaborators: R. Halburd (London), R. Vidunas (Tokyo), R. Kycia (Kraków) 1 Plan
More informationLois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique
Lois de conservations scalaires hyperboliques stochastiques : existence, unicité et approximation numérique de la solution entropique Université Aix-Marseille Travail en collaboration avec C.Bauzet, V.Castel
More informationarxiv:cs/ v1 [cs.na] 23 Aug 2004
arxiv:cs/0408053v1 cs.na 23 Aug 2004 WEIGHTED AVERAGE FINITE DIFFERENCE METHODS FOR FRACTIONAL DIFFUSION EQUATIONS Santos B. Yuste,1 Departamento de Física, Universidad de Extremadura, E-0071 Badaoz, Spain
More informationConsistency of the maximum likelihood estimator for general hidden Markov models
Consistency of the maximum likelihood estimator for general hidden Markov models Jimmy Olsson Centre for Mathematical Sciences Lund University Nordstat 2012 Umeå, Sweden Collaborators Hidden Markov models
More informationSemi-Linear Diusive Representations for Non-Linear Fractional Dierential Systems Jacques AUDOUNET 1, Denis MATIGNON 2?, and Gerard MONTSENY 3 1 MIP /
Semi-Linear Diusive Representations for Non-Linear Fractional Dierential Systems Jacques AUDOUNET 1, Denis MATIGNON?, and Gerard MONTSENY 3 1 MIP / CNRS, Universite Paul Sabatier. 118, route de Narbonne.
More informationAn Efficient Numerical Method for Solving. the Fractional Diffusion Equation
Journal of Applied Mathematics & Bioinformatics, vol.1, no.2, 2011, 1-12 ISSN: 1792-6602 (print), 1792-6939 (online) International Scientific Press, 2011 An Efficient Numerical Method for Solving the Fractional
More informationIntroduction Multiphase ows in porous media occur in various practical situations including oil recovery and unsaturated groundwater ow, see for insta
Analysis of a Darcy ow model with a dynamic pressure saturation relation Josephus Hulshof Mathematical Department of the Leiden University Niels Bohrweg, 333 CA Leiden, The Netherlands E-mail: hulshof@wi.leidenuniv.nl
More informationNumerical Treatment of Fractional Dierential Equations
Dipartimento di Matematica e Fisica Sezione di Matematica Dottorato di Ricerca in Matematica XXVI Ciclo Numerical Treatment of Fractional Dierential Equations Candidate: Moreno Concezzi Advisor: Prof.
More informationStatistical Learning Theory
Statistical Learning Theory Fundamentals Miguel A. Veganzones Grupo Inteligencia Computacional Universidad del País Vasco (Grupo Inteligencia Vapnik Computacional Universidad del País Vasco) UPV/EHU 1
More informationTRACE FORMULAS FOR PERTURBATIONS OF OPERATORS WITH HILBERT-SCHMIDT RESOLVENTS. Bishnu Prasad Sedai
Opuscula Math. 38, no. 08), 53 60 https://doi.org/0.7494/opmath.08.38..53 Opuscula Mathematica TRACE FORMULAS FOR PERTURBATIONS OF OPERATORS WITH HILBERT-SCHMIDT RESOLVENTS Bishnu Prasad Sedai Communicated
More informationBifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces
Nonlinear Analysis 42 (2000) 561 572 www.elsevier.nl/locate/na Bifurcation from the rst eigenvalue of some nonlinear elliptic operators in Banach spaces Pavel Drabek a;, Nikos M. Stavrakakis b a Department
More informationApproximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle
Malaya J. Mat. 4(1)(216) 8-18 Approximating solutions of nonlinear second order ordinary differential equations via Dhage iteration principle B. C. Dhage a,, S. B. Dhage a and S. K. Ntouyas b,c, a Kasubai,
More informationContractive metrics for scalar conservation laws
Contractive metrics for scalar conservation laws François Bolley 1, Yann Brenier 2, Grégoire Loeper 34 Abstract We consider nondecreasing entropy solutions to 1-d scalar conservation laws and show that
More informationNodal O(h 4 )-superconvergence of piecewise trilinear FE approximations
Preprint, Institute of Mathematics, AS CR, Prague. 2007-12-12 INSTITTE of MATHEMATICS Academy of Sciences Czech Republic Nodal O(h 4 )-superconvergence of piecewise trilinear FE approximations Antti Hannukainen
More informationAlignment processes on the sphere
Alignment processes on the sphere Amic Frouvelle CEREMADE Université Paris Dauphine Joint works with : Pierre Degond (Imperial College London) and Gaël Raoul (École Polytechnique) Jian-Guo Liu (Duke University)
More informationSome nonlinear elliptic equations in R N
Nonlinear Analysis 39 000) 837 860 www.elsevier.nl/locate/na Some nonlinear elliptic equations in Monica Musso, Donato Passaseo Dipartimento di Matematica, Universita di Pisa, Via Buonarroti,, 5617 Pisa,
More informationON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS. G. Makay
ON THE ASYMPTOTIC STABILITY IN TERMS OF TWO MEASURES FOR FUNCTIONAL DIFFERENTIAL EQUATIONS G. Makay Student in Mathematics, University of Szeged, Szeged, H-6726, Hungary Key words and phrases: Lyapunov
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
More informationOn the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals
Journal of Computational and Applied Mathematics 49 () 38 395 www.elsevier.com/locate/cam On the convergence of interpolatory-type quadrature rules for evaluating Cauchy integrals Philsu Kim a;, Beong
More informationA large deviation principle for a RWRC in a box
A large deviation principle for a RWRC in a box 7th Cornell Probability Summer School Michele Salvi TU Berlin July 12, 2011 Michele Salvi (TU Berlin) An LDP for a RWRC in a nite box July 12, 2011 1 / 15
More informationCRITICAL ATOM NUMBER OF A HARMONICALLY TRAPPED 87 Rb BOSE GAS AT DIFFERENT TEMPERATURES
CRITICAL ATOM UMBER OF A HARMOICALLY TRAPPED 87 Rb BOSE GAS AT DIFFERET TEMPERATURES AHMED S. HASSA, AZZA M. EL-BADRY Department of Physics, Faculty of Science, El Minia University, El Minia, Egypt E-mail:
More informationA Representation of Excessive Functions as Expected Suprema
A Representation of Excessive Functions as Expected Suprema Hans Föllmer & Thomas Knispel Humboldt-Universität zu Berlin Institut für Mathematik Unter den Linden 6 10099 Berlin, Germany E-mail: foellmer@math.hu-berlin.de,
More informationMonte Carlo Methods for Statistical Inference: Variance Reduction Techniques
Monte Carlo Methods for Statistical Inference: Variance Reduction Techniques Hung Chen hchen@math.ntu.edu.tw Department of Mathematics National Taiwan University 3rd March 2004 Meet at NS 104 On Wednesday
More informationKirsten Lackner Solberg. Dept. of Math. and Computer Science. Odense University, Denmark
Inference Systems for Binding Time Analysis Kirsten Lackner Solberg Dept. of Math. and Computer Science Odense University, Denmark e-mail: kls@imada.ou.dk June 21, 1993 Contents 1 Introduction 4 2 Review
More informationMapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing
NASA/CR--1631 ICASE Report No. -48 Mapping Closure Approximation to Conditional Dissipation Rate for Turbulent Scalar Mixing Guowei He and R. Rubinstein ICASE, Hampton, Virginia ICASE NASA Langley Research
More informationmaximally charged black holes and Hideki Ishihara Department ofphysics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo 152, Japan
Quasinormal modes of maximally charged black holes Hisashi Onozawa y,takashi Mishima z,takashi Okamura, and Hideki Ishihara Department ofphysics, Tokyo Institute of Technology, Oh-okayama, Meguro, Tokyo
More informationAn Operator Theoretical Approach to Nonlocal Differential Equations
An Operator Theoretical Approach to Nonlocal Differential Equations Joshua Lee Padgett Department of Mathematics and Statistics Texas Tech University Analysis Seminar November 27, 2017 Joshua Lee Padgett
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 informationThis method is introduced by the author in [4] in the case of the single obstacle problem (zero-obstacle). In that case it is enough to consider the v
Remarks on W 2;p -Solutions of Bilateral Obstacle Problems Srdjan Stojanovic Department of Mathematical Sciences University of Cincinnati Cincinnati, OH 4522-0025 November 30, 995 Abstract The well known
More informationAbstract. Previous characterizations of iss-stability are shown to generalize without change to the
On Characterizations of Input-to-State Stability with Respect to Compact Sets Eduardo D. Sontag and Yuan Wang Department of Mathematics, Rutgers University, New Brunswick, NJ 08903, USA Department of Mathematics,
More informationNumerical Solution of Linear Fredholm Integro-Differential Equations by Non-standard Finite Difference Method
Available at http://pvamu.edu/aam Appl. Appl. Math. ISS: 1932-9466 Vol. 10, Issue 2 (December 2015), pp. 1019-1026 Applications and Applied Mathematics: An International Journal (AAM) umerical Solution
More informationProceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005
Proceedings of the 5th International Conference on Inverse Problems in Engineering: Theory and Practice, Cambridge, UK, 11-15th July 2005 SOME INVERSE SCATTERING PROBLEMS FOR TWO-DIMENSIONAL SCHRÖDINGER
More informationPutzer s Algorithm. Norman Lebovitz. September 8, 2016
Putzer s Algorithm Norman Lebovitz September 8, 2016 1 Putzer s algorithm The differential equation dx = Ax, (1) dt where A is an n n matrix of constants, possesses the fundamental matrix solution exp(at),
More informationThe viscosity technique for the implicit midpoint rule of nonexpansive mappings in
Xu et al. Fixed Point Theory and Applications 05) 05:4 DOI 0.86/s3663-05-08-9 R E S E A R C H Open Access The viscosity technique for the implicit midpoint rule of nonexpansive mappings in Hilbert spaces
More informationA higher-order Boussinesq equation in locally non-linear theory of one-dimensional non-local elasticity
IMA Journal of Applied Mathematics 29 74, 97 16 doi:1.193/imamat/hxn2 Advance Access publication on July 27, 28 A higher-order Boussinesq equation in locally non-linear theory of one-dimensional non-local
More informationLinearized methods for ordinary di erential equations
Applied Mathematics and Computation 104 (1999) 109±19 www.elsevier.nl/locate/amc Linearized methods for ordinary di erential equations J.I. Ramos 1 Departamento de Lenguajes y Ciencias de la Computacion,
More informationNumerical Solution of Second Order Singularly Perturbed Differential Difference Equation with Negative Shift. 1 Introduction
ISSN 749-3889 print), 749-3897 online) International Journal of Nonlinear Science Vol.8204) No.3,pp.200-209 Numerical Solution of Second Order Singularly Perturbed Differential Difference Equation with
More informationFRACTIONAL BOUNDARY VALUE PROBLEMS ON THE HALF LINE. Assia Frioui, Assia Guezane-Lakoud, and Rabah Khaldi
Opuscula Math. 37, no. 2 27), 265 28 http://dx.doi.org/.7494/opmath.27.37.2.265 Opuscula Mathematica FRACTIONAL BOUNDARY VALUE PROBLEMS ON THE HALF LINE Assia Frioui, Assia Guezane-Lakoud, and Rabah Khaldi
More information5 Banach Algebras. 5.1 Invertibility and the Spectrum. Robert Oeckl FA NOTES 5 19/05/2010 1
Robert Oeckl FA NOTES 5 19/05/2010 1 5 Banach Algebras 5.1 Invertibility and the Spectrum Suppose X is a Banach space. Then we are often interested in (continuous) operators on this space, i.e, elements
More informationA brief introduction to ordinary differential equations
Chapter 1 A brief introduction to ordinary differential equations 1.1 Introduction An ordinary differential equation (ode) is an equation that relates a function of one variable, y(t), with its derivative(s)
More informationNumerical Methods. King Saud University
Numerical Methods King Saud University Aims In this lecture, we will... find the approximate solutions of derivative (first- and second-order) and antiderivative (definite integral only). Numerical Differentiation
More informationNew Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems
Systems & Control Letters 43 (21 39 319 www.elsevier.com/locate/sysconle New Lyapunov Krasovskii functionals for stability of linear retarded and neutral type systems E. Fridman Department of Electrical
More informationHopfLax-Type Formula for u t +H(u, Du)=0*
journal of differential equations 126, 4861 (1996) article no. 0043 HopfLax-Type Formula for u t +H(u, Du)=0* E. N. Barron, - R. Jensen, and W. Liu 9 Department of Mathematical Sciences, Loyola University,
More informationOn Symmetry Group Properties of the Benney Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2004, Vol. 50, Part 1, 214 218 On Symmetry Group Properties of the Benney Equations Teoman ÖZER Istanbul Technical University, Faculty of Civil
More informationAdvances in Difference Equations 2012, 2012:7
Advances in Difference Equations This Provisional PDF corresponds to the article as it appeared upon acceptance. Fully formatted PDF and full text (HTML) versions will be made available soon. Extremal
More informationScaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations
Scaling Limits of Waves in Convex Scalar Conservation Laws under Random Initial Perturbations Jan Wehr and Jack Xin Abstract We study waves in convex scalar conservation laws under noisy initial perturbations.
More informationNumerical analysis of the one-dimensional Wave equation subject to a boundary integral specification
Numerical analysis of the one-dimensional Wave equation subject to a boundary integral specification B. Soltanalizadeh a, Reza Abazari b, S. Abbasbandy c, A. Alsaedi d a Department of Mathematics, University
More informationStochastic Dynamic Programming. Jesus Fernandez-Villaverde University of Pennsylvania
Stochastic Dynamic Programming Jesus Fernande-Villaverde University of Pennsylvania 1 Introducing Uncertainty in Dynamic Programming Stochastic dynamic programming presents a very exible framework to handle
More informationBackward error analysis
Backward error analysis Brynjulf Owren July 28, 2015 Introduction. The main source for these notes is the monograph by Hairer, Lubich and Wanner [2]. Consider ODEs in R d of the form ẏ = f(y), y(0) = y
More informationApplications of the compensated compactness method on hyperbolic conservation systems
Applications of the compensated compactness method on hyperbolic conservation systems Yunguang Lu Department of Mathematics National University of Colombia e-mail:ylu@unal.edu.co ALAMMI 2009 In this talk,
More informationBoundary layers in a two-point boundary value problem with fractional derivatives
Boundary layers in a two-point boundary value problem with fractional derivatives J.L. Gracia and M. Stynes Institute of Mathematics and Applications (IUMA) and Department of Applied Mathematics, University
More informationBoundary value problems for fractional differential equations with three-point fractional integral boundary conditions
Sudsutad and Tariboon Advances in Difference Equations 212, 212:93 http://www.advancesindifferenceequations.com/content/212/1/93 R E S E A R C H Open Access Boundary value problems for fractional differential
More informationThe Uniformity Principle: A New Tool for. Probabilistic Robustness Analysis. B. R. Barmish and C. M. Lagoa. further discussion.
The Uniformity Principle A New Tool for Probabilistic Robustness Analysis B. R. Barmish and C. M. Lagoa Department of Electrical and Computer Engineering University of Wisconsin-Madison, Madison, WI 53706
More informationSPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS
SPECTRAL PROPERTIES AND NODAL SOLUTIONS FOR SECOND-ORDER, m-point, BOUNDARY VALUE PROBLEMS BRYAN P. RYNNE Abstract. We consider the m-point boundary value problem consisting of the equation u = f(u), on
More information3.1 Basic properties of real numbers - continuation Inmum and supremum of a set of real numbers
Chapter 3 Real numbers The notion of real number was introduced in section 1.3 where the axiomatic denition of the set of all real numbers was done and some basic properties of the set of all real numbers
More informationPathwise Construction of Stochastic Integrals
Pathwise Construction of Stochastic Integrals Marcel Nutz First version: August 14, 211. This version: June 12, 212. Abstract We propose a method to construct the stochastic integral simultaneously under
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationInstitute for Advanced Computer Studies. Department of Computer Science. On the Convergence of. Multipoint Iterations. G. W. Stewart y.
University of Maryland Institute for Advanced Computer Studies Department of Computer Science College Park TR{93{10 TR{3030 On the Convergence of Multipoint Iterations G. W. Stewart y February, 1993 Reviseed,
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 informationGarrett: `Bernstein's analytic continuation of complex powers' 2 Let f be a polynomial in x 1 ; : : : ; x n with real coecients. For complex s, let f
1 Bernstein's analytic continuation of complex powers c1995, Paul Garrett, garrettmath.umn.edu version January 27, 1998 Analytic continuation of distributions Statement of the theorems on analytic continuation
More informationExistence, Uniqueness and Stability of Hilfer Type Neutral Pantograph Differential Equations with Nonlocal Conditions
International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 6, Issue 8, 2018, PP 42-53 ISSN No. (Print) 2347-307X & ISSN No. (Online) 2347-3142 DOI: http://dx.doi.org/10.20431/2347-3142.0608004
More informationThe generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations.
The generalized Euler-Maclaurin formula for the numerical solution of Abel-type integral equations. Johannes Tausch Abstract An extension of the Euler-Maclaurin formula to singular integrals was introduced
More informationIntroduction to the Numerical Solution of IVP for ODE
Introduction to the Numerical Solution of IVP for ODE 45 Introduction to the Numerical Solution of IVP for ODE Consider the IVP: DE x = f(t, x), IC x(a) = x a. For simplicity, we will assume here that
More informationFUZZY SOLUTIONS FOR BOUNDARY VALUE PROBLEMS WITH INTEGRAL BOUNDARY CONDITIONS. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXV 1(26 pp. 119 126 119 FUZZY SOLUTIONS FOR BOUNDARY VALUE PROBLEMS WITH INTEGRAL BOUNDARY CONDITIONS A. ARARA and M. BENCHOHRA Abstract. The Banach fixed point theorem
More informationSTOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER. Department of Mathematics. University of Wisconsin.
STOCHASTIC DIFFERENTIAL EQUATIONS WITH EXTRA PROPERTIES H. JEROME KEISLER Department of Mathematics University of Wisconsin Madison WI 5376 keisler@math.wisc.edu 1. Introduction The Loeb measure construction
More informationAnalysis of Fractional Differential Equations. Kai Diethelm & Neville J. Ford
ISSN 136-1725 UMIST Analysis of Fractional Differential Equations Kai Diethelm & Neville J. Ford Numerical Analysis Report No. 377 A report in association with Chester College Manchester Centre for Computational
More informationA symmetry analysis of Richard s equation describing flow in porous media. Ron Wiltshire
BULLETIN OF THE GREEK MATHEMATICAL SOCIETY Volume 51, 005 93 103) A symmetry analysis of Richard s equation describing flow in porous media Ron Wiltshire Abstract A summary of the classical Lie method
More informationHigh Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional Differential Equation
International Symposium on Fractional PDEs: Theory, Numerics and Applications June 3-5, 013, Salve Regina University High Order Numerical Methods for the Riesz Derivatives and the Space Riesz Fractional
More information