du dθ We wish to approximate the initial value problem for the fractional ordinary differential equation of order β :
|
|
- Mervin Lawrence
- 5 years ago
- Views:
Transcription
1 Mixed Collocation for Fractional Differential Equations François Dubois and Stéphanie Mengué December 5, 3 This contribution is dedicaced to Claude Brézinski at the occasion of his sixtieth birthday. Abstract We present the mixed collocation method for numerical integration of fractional differential equations of the type D β u Φ(u, t). Given a regular mesh with constant discretization step, the unknown u(t) is considered as continuous and affine in each cell, and the dynamics Φ(u, t) as a constant. After a fractional integration, the equation is written strongly at the mesh vertices and the dynamics weakly in each cell. The Semidif software has been developed for the particular case of numerical integration of order. The validation for analytical results and published solutions is established and experimental convergence as the mesh size tends to zero is obtained. Good results are obtained for a nonlinear model with a strong singularity. Keywords : Mixed finite elements, numerical combustion. AMS(MOS) classification : 65N3. Applications Scientifiques du Calcul Intensif, bâtiment 56, B.P. 67, F-943 Orsay, Union Européenne. Mel : dubois@asci.fr. Conservatoire National des Arts et Métiers, Equipe de recherche associée n o 396, 5, rue Marat, F-78 Saint-Cyr-L Ecole, Union Européenne. Université de Marne-La-Vallée, Laboratoire Systèmes de Communication, bâtiment Copernic, F Marne-La-Vallée cedex, Union Européenne. Mel : smengue@univ-mlv.fr. Numerical Algorithms, volume 34, numéro, p. 33-3, dec 3, édition du mars 9.
2 François Dubois and Stéphanie Mengué ) Introduction Let β be a real number, < β <, Γ( ) the classical Euler function, and IR [, [ (u,t) Φ(u,t) IR a regular mapping. Following for example the now classical work of Caputo [Ca67], we define the fractional differential operator D β ( ) by () (D β t du dθ u)(t) Γ( β) dθ (t θ) β. We wish to approximate the initial value problem for the fractional ordinary differential equation of order β : { () D β (u u ) Φ(u (t),t), t > u u, t. We introduce a discretization step h, and denote by K h the interval of j+ the form K h j+ ]t j, t j+ [ ]jh, (j + )h[ for j IN. We use the discrete spaces P h and Q h defined as follows. The linear space P h is composed by functions that are continuous and affine in each mesh element K h. A j+ function v P h can be expanded on the basis ( ) ϕ h k of classical hat k IN functions : { ϕ h (θ) θ h, θ h, θ h ( θ t k ), t k θ t k+, ϕ h k (θ) h + ( θ t k ), t k θ t k, k, h, θ t k h and by using the nodal values v k v ( t k) at the grid vertices, we have : v k v k ϕ h k for v P h. The discrete space Q h is composed by functions that take a constant value in each element K h. The basis functions of the discrete space Q h j+ are null everywhere except inside the interval K h where it is equal to the j+ unity. Then, w Q h can be expanded on the χ h basis according to the k+ relation w k w k+ χh. k+ ) Mixed collocation scheme We integrate the differential equation () with the fractional integrator I β of order β defined by
3 Mixed Collocation for Fractional Differential Equations I β t (v ( ),t) Γ(β) (t θ)β v (θ) dθ and satisfying the functional equation I β D β Id when applied for a regular function u satisfying u(). The relation () can therefore be rewritten as (3) u (t) u I β (Φ(u ( ), t)), t. The P Q mixed collocation method consists in choosing a discrete state u h ( ) satisfying (4) u h u h j ϕ h j, u h P h, j and the so-called flux Φ(u ( ),t) by a discontinuous time function f h ( ) that is constant in each cell K h : j+ (5) f h, f h Q h. j f h j+ χ h j+ On one hand, the equation (3) is written strongly at the mesh vertices jh (j IN) and we can speak of a collocation method : (6) u h (jh) u I β ( f h ( ), jh ), j IN. On the other hand, it is a method inspired from mixed finite elements. We can not impose that for each time value t, we have f h (t) Φ(u h (t), t). Then the P Q mixed collocation method claims only that for the mean values, the approached flux f h ( ) and the exact flux Φ ( u ( ), ) h of the approached solution have the same mean value in each interval ]jh, (j +)h[. We impose : (j+)h jh f h (θ) dθ (j+)h jh Φ ( u h (θ),θ ) dθ, j IN. This projection step on the discrete space Q h can be written as : (7) f h j+ Φ ( u h j ( θ) + θu h j+, jh + θh ) dθ, j IN. Proposition. State-flux constraint for the P Q scheme With the choice (4) (5), the relation (6) can be written in the following way : (8) u h j+ h β Γ(β + ) fh j+ u + 3 h β Γ(β + ) j k α j k f h k+, j IN,
4 François Dubois and Stéphanie Mengué with (9) α k (k + ) β k β, k IN. Proof of proposition. We evaluate the right hand side of relation (6) for the discrete time t j+. We have : I β ( f h ( ),t j+) Γ(β) Γ(β) Γ(β) Γ(β) t j + j k j k j k Γ(β + ) h β Γ(β + ) ( t j+ θ ) β f h (θ)dθ (k+)h kh f h k+ f h k+ j k j k (t j+ θ) β f h (θ) dθ [(j + ) h (kh + hξ)] β h dξ h β ( β f h k+ and the relation (8) is a direct consequence of (6). ) [ (j + k ξ) β] ξ ξ h β [ (j + k) β (j k) β] f h k+ α j k taking into account (9) 3) First numerical tests We have compared in [DMk] several numerical schemes i.e. two and three point Grünwald-Letnikov scheme, Msallam scheme and P finite element scheme and we have also compared in [DM] thoses schemes with the mixed collocation method. These tests are done for a semi-derivation, i.e. β.5. The dynamics (u,t) Φ(u,t) is parameterized by the functions u f (u) and t g (t) and by the real numbers θ f and θ g : Φ(u,t) θ f f (u) + θ g g (t). The parameters θ f and θ g are chosen in order to obtain as solution the following particular functions : ( ) j () u j (t) t, j,,5. In a first approach, we consider θ f and θ g ; a simple semi-quadrature has to be done. The equation () corresponds to the following t g (t) choice : g (t) π, g (t) π t,, g3 (t) 3 4 π t, g4 (t) 8 3 π t t and g 5 (t) 5 6 π t. 4
5 Mixed Collocation for Fractional Differential Equations In a second approach, we take θ f and θ g. This dynamics is associated to the following functions : f (u) π, f (u) u, f3 (u) 3 π u 3, π 4 () f 4 (u) 8 3 π u 3 4, f5 (t) 5 4 π u 5. 6 In a third approach, whave also considered the following numerical test cases previously proposed by K. Diethelm [Di97], Diethelm and Ford [DF], Diethelm and Luchko [DL3], Blank [Bl96] and Lubich [Lu86] : () u, t, D / u u π t 3, t > : u (t) t ( (3) u, t, D / u u, t > : u (t) e t erf ) t + ( (4) u, t, D / u u, t > : u (t) e t erf ) t (5) u, t, D / u sin t, t > : u (t) ( ) tj t. π We remark that for the nonlinear cases proposed in (), the sufficient Lipschitz continuity condition for the function Φ(, ), revisited by Nagumo and presented e.g. the book [GV9] of Gorenflo and Vessela, that assumes that the nonlinear Abel-Volterra equation (3) has a unique continuous solution, is a priori not satisfied. Nevertheless, the functions presented at the relations () are clearly solution of problem (3) and our numerical results have captured this fact. The numerical simulation is done with our Fortran software Semidif. We have computed systematically the orders of convergence for mesh step sizes h of the type for 3 n 3. We have measured the L and L errors denotes respectively n by e n and e n and defined according to : e n max{ u(jh) u j, j,, n } and e n h u () u + n j ( ) j u u j + u () u n. The orders of convergence are summarized in the following table : n g (t) g (t) g 3 (t) g 4 (t) g 5 (t) () (3) (4) L L
6 François Dubois and Stéphanie Mengué Our results show a good agreement with the studies done by these authors. We refer to our report [DM] for the details. 4) Nonlinear model with a singularity We explain in this section how to consider the P Q method when we study the Joulin model [Jo85] for spherical flammes : { (6) D / (u) Φ(u (t),t), t > u, t (7) Φ(u (t),t) log u + E tγ ( t) H( t), u where E and γ.3 are positive constants and θ H(θ) the Heaviside fonction equal to if θ < end to if θ >. Proposition. Computation of the first point for Joulin equation For the equation (6), and the representation (8) of the unknown u h, the first unknown point u h is solution of the following equation whose unknown is denoted by x : h (8) x {log x + E ( hγ π x γ h )}. γ + Proof of proposition. We integration the equation (6) one half time on the interval [, h] : (9) u h h π (h θ) Φ ( u h (θ),θ ) dθ. Then we make the hypothesis that the flux Φ(, ) is equal to a constant f h on the interval [,h] : f h h h Φ( u h (θ),θ ) dθ. We inject this hypothesis at the right hand side of equation (9) : () u h h Φ ( u h (θ),θ ) dθ π h and we compute the term on the right of () by using the function Φ(, ) of the relation (7). We note the singularity in the vincinity of zero. The hypothesis of affine representation of u between and h, show the following expression : Φ ( u h (θ),θ ) ( ) log u h θ + E θ γ ( θ) h u h θ if h h log ( u h ) θ + log h + E h θγ ( θ) u h. 6
7 Mixed Collocation for Fractional Differential Equations Then after integration, we obtain h Φ ( u h (θ),θ ) dθ h log ( u h ) + log tdt + E u h h γ t γ ( ht) dt log ( u h ) + [t log t t] + E [ u h γ h ] γ + log ( ( ) u h E + u h γ h ). γ + and the proposition is established. The equation (8) is solved with the Newton algorithm. For the other grid points, some algebra and a numerical integration (see the details in [DM]) show that we have to solve the system composed by (8) and the following representation : f h j+ u h j+ uh j + E hγ [ u h j+ log u h j+ u h ] j+ ( j γ ( jh) u h j ) + (j + )γ ( (j + )h) u h j points 4 points 394 points 393 points 393 points 39 points 38 points 3 points 5 points Change of comportment between np393 and np393 mesh points 6 Simulations with E7.6 and t max Figure. Qualitative change of the solution according to various discretizations. We have followed in [DM] the working plan proposed by Audounet and Roquejoffre [AR98] for their simulations of the semidifferential equation (6) 7
8 François Dubois and Stéphanie Mengué (7) : γ.3 and 7.6 E 7.7. According to these authors, a bifurcation in the asymptotic comportment occurs for the solution of model (6)(7) with γ.3. For E 7.7, we have extinction of the flamme in a finite time whereas for E 7.8, its radius tends to infinity as time tends to infinity. We first observe that we have to control very precisely the discretization step h, as proposed in the following experiment. We consider the parameter value E 7.6 and t max 4. When we use less than N 393 mesh points, the comportment of the discrete solution is to tends to infinity and for more than N 393 mesh points, we observe an extinction in a finite time (see the figure )! points 5 points 5 points points 5 points 5 points points E Figure. Numerical solution of the model (6) (7), E In what follows, we systematically refine the mesh in order to make in evidence the independence of the numerical solution with the number of grid points. When we use sufficiently refined meshes, we observe (see [DM]) an extinction in finite time for E 7.6 and a diverging process for E 7.7. We have experimentally determined the value of bifurcation parameter E. We use a number of mesh points N according to the rule N, k,.5 k, 5 k, k+,. The same conclusions hold for the extinction parameter E 7.66 and we have a diverging flame for E 7.67 (Figure ), but the numerical proof needs meshes that use up to 5 points. The case E (Figure 3) shows an extinction of the flamme in a finite time, but the high cost of computer ressources (more than 8
9 Mixed Collocation for Fractional Differential Equations two millions mesh points to observe asymptotically numerical convergence with the time discretization!) makes in evidence the actual limitations of the method. Therefore, according to our numerical results, the bifurcation parameter E for the model (6) (7) satisfies the inequalities : < E < 7.67 with γ points 5 points 5 points points 5 points 5 points points 5 points 5 points points 5 points E Figure 3. Numerical solution of the model (6) (7), E ) Conclusion and aknowledgments In our Semidif software, the system of equations (7) (8) that characterizes the P Q mixed collocation method is solved with a Newton method. We have validated our approach with exact solutions and recent published numerical approximations about fractional differential equations. For a difficult test case issued from combustion modelling, we show that our mixed collocation method is relevant and conducts to significative results. We have noticed an accuracy of order for nonlinear test problems whereas the interpolation process is of order. The natural question is therefore to obtain second order accuracy in the future. Last but not least, we thank the unknown referee for helpful scientific comments on the first version of this article. 6) References. [AR98] J. Audounet and J.M. Roquejoffre. An asymptotic fractional differential model of spherical flame, ESAIM Proceedings, vol. 5, p. 5-7,
10 [Bl96] [Ca67] [DF] [Di97] [DL3] François Dubois and Stéphanie Mengué L. Blank. Numerical treatment of differential equations of fractional order, Numerical analysis report n o 87, Manchester University, March 996. M. Caputo. Linear models of dissipation whose Q is almost frequency independent, Part, Geophys. J.R. Astr. Soc., vol. 3, p , 967. K. Diethelm and N. Ford. Analysis of fractional differential equations, J. Math. Anal. Appl., vol. 65, p. 9-48,. K. Diethelm. An algorithm for the numerical solution of differential equations of fractional order, Electronic Transactions on Numerical Analysis, vol. 5, p. -6, 997. K. Diethelm and Y. Luchko. Numerical solution of linear multiterm differential equations of fractional order, J. Comput. Anal. Appl., 3, to appear. [DMk] F. Dubois and S. Mengué. Schémas numériques implicites pour les équations semi-différentielles, IAT-CNAM internal report n o 334, see also June. [DM] F. Dubois and S. Mengué. Collocation mixte pour les équations différentielles non-linéaires d ordre fractionnaire, IAT-CNAM internal report n o 348 and Semidif software, July. [GV9] R. Gorenflo and S. Vessella. Abel Integral Equations : Analysis and Applications, Lecture Notes in Mathematics, vol. 46, Springer, New York, 99. [Jo85] [Lu86] G. Joulin. Point source initiation of lean spherical flames of light reactants: an asymptotic theory, Combustion Science and Technology, vol. 43, p. 99-3, 985. C. Lubich. Discretized fractional calculus, SIAM Journal on Mathematical Analalysis, vol. 7, n o 3, p , 986.
An adaptation of the Gear scheme for fractional derivatives
An adaptation of the Gear scheme for fractional derivatives Ana Cristina Galucio, Jean-François Deü, Stéphanie Mengué, François Dubois To cite this version: Ana Cristina Galucio, Jean-François Deü, Stéphanie
More informationUnconditionally stable scheme for Riccati equation
ESAIM: Proceedings, Vol. 8, 2, 39-52 Contrôle des systèmes gouvernés par des équations aux dérivées partielles http://www.emath.fr/proc/vol.8/ Unconditionally stable scheme for Riccati equation François
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 informationHigher order numerical methods for solving fractional differential equations
Noname manuscript No. will be inserted by the editor Higher order numerical methods for solving fractional differential equations Yubin Yan Kamal Pal Neville J Ford Received: date / Accepted: date Abstract
More informationA Survey of Numerical Methods in Fractional Calculus
A Survey of Numerical Methods in Fractional Calculus Fractional Derivatives in Mechanics: State of the Art CNAM Paris, 17 November 2006 Kai Diethelm diethelm@gns-mbh.com, k.diethelm@tu-bs.de Gesellschaft
More informationHomographic scheme for Riccati equation
Homographic scheme for Riccati equation François Dubois and Abdelkader Saïdi Abstract In this paper we present a numerical scheme for the resolution of matrix Riccati equation, usualy used in control problems.
More informationBLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED
BLOWUP THEORY FOR THE CRITICAL NONLINEAR SCHRÖDINGER EQUATIONS REVISITED TAOUFIK HMIDI AND SAHBI KERAANI Abstract. In this note we prove a refined version of compactness lemma adapted to the blowup analysis
More informationTHE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS
JOURNAL OF INTEGRAL EQUATIONS AND APPLICATIONS Volume 4, Number 1, Winter 1992 THE METHOD OF LINES FOR PARABOLIC PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS J.-P. KAUTHEN ABSTRACT. We present a method of lines
More informationHigh-order ADI schemes for convection-diffusion equations with mixed derivative terms
High-order ADI schemes for convection-diffusion equations with mixed derivative terms B. Düring, M. Fournié and A. Rigal Abstract We consider new high-order Alternating Direction Implicit ADI) schemes
More informationSuperconvergence analysis of multistep collocation method for delay Volterra integral equations
Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 4, No. 3, 216, pp. 25-216 Superconvergence analysis of multistep collocation method for delay Volterra integral equations
More informationMEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
PROCEEDINGS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 40, Number 7, July 0, Pages 453 463 S 000-9939(0)8-X Article electronically published on November, 0 MEAN VALUE PROPERTY FOR p-harmonic FUNCTIONS
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 informationConvergence rate estimates for the gradient differential inclusion
Convergence rate estimates for the gradient differential inclusion Osman Güler November 23 Abstract Let f : H R { } be a proper, lower semi continuous, convex function in a Hilbert space H. The gradient
More informationFractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials
Applied Mathematical Sciences, Vol. 5, 211, no. 45, 227-2216 Fractional Calculus for Solving Abel s Integral Equations Using Chebyshev Polynomials Z. Avazzadeh, B. Shafiee and G. B. Loghmani Department
More informationarxiv: v1 [math.ap] 20 Nov 2007
Long range scattering for the Maxwell-Schrödinger system with arbitrarily large asymptotic data arxiv:0711.3100v1 [math.ap] 20 Nov 2007 J. Ginibre Laboratoire de Physique Théorique Université de Paris
More informationA novel difference schemes for analyzing the fractional Navier- Stokes equations
DOI: 0.55/auom-207-005 An. Şt. Univ. Ovidius Constanţa Vol. 25(),207, 95 206 A novel difference schemes for analyzing the fractional Navier- Stokes equations Khosro Sayevand, Dumitru Baleanu, Fatemeh Sahsavand
More informationSelf-similar solutions for the diffraction of weak shocks
Self-similar solutions for the diffraction of weak shocks Allen M. Tesdall John K. Hunter Abstract. We numerically solve a problem for the unsteady transonic small disturbance equations that describes
More informationMath 1B Final Exam, Solution. Prof. Mina Aganagic Lecture 2, Spring (6 points) Use substitution and integration by parts to find:
Math B Final Eam, Solution Prof. Mina Aganagic Lecture 2, Spring 20 The eam is closed book, apart from a sheet of notes 8. Calculators are not allowed. It is your responsibility to write your answers clearly..
More informationDifferentiability with respect to initial data for a scalar conservation law
Differentiability with respect to initial data for a scalar conservation law François BOUCHUT François JAMES Abstract We linearize a scalar conservation law around an entropy initial datum. The resulting
More informationCollocation and iterated collocation methods for a class of weakly singular Volterra integral equations
Journal of Computational and Applied Mathematics 229 (29) 363 372 Contents lists available at ScienceDirect Journal of Computational and Applied Mathematics journal homepage: www.elsevier.com/locate/cam
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 informationON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * 1. Introduction
Journal of Computational Mathematics, Vol.6, No.6, 008, 85 837. ON SPECTRAL METHODS FOR VOLTERRA INTEGRAL EQUATIONS AND THE CONVERGENCE ANALYSIS * Tao Tang Department of Mathematics, Hong Kong Baptist
More informationAN EXAMPLE OF FUNCTIONAL WHICH IS WEAKLY LOWER SEMICONTINUOUS ON W 1,p FOR EVERY p > 2 BUT NOT ON H0
AN EXAMPLE OF FUNCTIONAL WHICH IS WEAKLY LOWER SEMICONTINUOUS ON W,p FOR EVERY p > BUT NOT ON H FERNANDO FARRONI, RAFFAELLA GIOVA AND FRANÇOIS MURAT Abstract. In this note we give an example of functional
More informationarxiv: v1 [math.ca] 31 Dec 2018
arxiv:181.1173v1 [math.ca] 31 Dec 18 Some trigonometric integrals and the Fourier transform of a spherically symmetric exponential function Hideshi YAMANE Department of Mathematical Sciences, Kwansei Gakuin
More informationNONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS. Shouchuan Hu Nikolas S. Papageorgiou. 1. Introduction
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 34, 29, 327 338 NONTRIVIAL SOLUTIONS FOR SUPERQUADRATIC NONAUTONOMOUS PERIODIC SYSTEMS Shouchuan Hu Nikolas S. Papageorgiou
More informationAn Attempt of Characterization of Functions With Sharp Weakly Complete Epigraphs
Journal of Convex Analysis Volume 1 (1994), No.1, 101 105 An Attempt of Characterization of Functions With Sharp Weakly Complete Epigraphs Jean Saint-Pierre, Michel Valadier Département de Mathématiques,
More informationPERIODIC SOLUTIONS OF THE FORCED PENDULUM : CLASSICAL VS RELATIVISTIC
LE MATEMATICHE Vol. LXV 21) Fasc. II, pp. 97 17 doi: 1.4418/21.65.2.11 PERIODIC SOLUTIONS OF THE FORCED PENDULUM : CLASSICAL VS RELATIVISTIC JEAN MAWHIN The paper surveys and compares some results on the
More informationLayer structures for the solutions to the perturbed simple pendulum problems
Layer structures for the solutions to the perturbed simple pendulum problems Tetsutaro Shibata Applied Mathematics Research Group, Graduate School of Engineering, Hiroshima University, Higashi-Hiroshima,
More informationQualitative Properties of Numerical Approximations of the Heat Equation
Qualitative Properties of Numerical Approximations of the Heat Equation Liviu Ignat Universidad Autónoma de Madrid, Spain Santiago de Compostela, 21 July 2005 The Heat Equation { ut u = 0 x R, t > 0, u(0,
More informationExistence of Minimizers for Fractional Variational Problems Containing Caputo Derivatives
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 3 12 (2013) http://campus.mst.edu/adsa Existence of Minimizers for Fractional Variational Problems Containing Caputo
More informationInverse Brascamp-Lieb inequalities along the Heat equation
Inverse Brascamp-Lieb inequalities along the Heat equation Franck Barthe and Dario Cordero-Erausquin October 8, 003 Abstract Adapting Borell s proof of Ehrhard s inequality for general sets, we provide
More informationGaussian interval quadrature rule for exponential weights
Gaussian interval quadrature rule for exponential weights Aleksandar S. Cvetković, a, Gradimir V. Milovanović b a Department of Mathematics, Faculty of Mechanical Engineering, University of Belgrade, Kraljice
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 informationNONLINEAR DECAY AND SCATTERING OF SOLUTIONS TO A BRETHERTON EQUATION IN SEVERAL SPACE DIMENSIONS
Electronic Journal of Differential Equations, Vol. 5(5), No. 4, pp. 7. ISSN: 7-669. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu ftp ejde.math.txstate.edu (login: ftp) NONLINEAR DECAY
More informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationEXISTENCE AND REGULARITY RESULTS FOR SOME NONLINEAR PARABOLIC EQUATIONS
EXISTECE AD REGULARITY RESULTS FOR SOME OLIEAR PARABOLIC EUATIOS Lucio BOCCARDO 1 Andrea DALL AGLIO 2 Thierry GALLOUËT3 Luigi ORSIA 1 Abstract We prove summability results for the solutions of nonlinear
More informationSpectrum and Exact Controllability of a Hybrid System of Elasticity.
Spectrum and Exact Controllability of a Hybrid System of Elasticity. D. Mercier, January 16, 28 Abstract We consider the exact controllability of a hybrid system consisting of an elastic beam, clamped
More informationON THE INVERSE FUNCTION THEOREM
PACIFIC JOURNAL OF MATHEMATICS Vol. 64, No 1, 1976 ON THE INVERSE FUNCTION THEOREM F. H. CLARKE The classical inverse function theorem gives conditions under which a C r function admits (locally) a C Γ
More informationEquilibrium analysis for a mass-conserving model in presence of cavitation
Equilibrium analysis for a mass-conserving model in presence of cavitation Ionel S. Ciuperca CNRS-UMR 58 Université Lyon, MAPLY,, 696 Villeurbanne Cedex, France. e-mail: ciuperca@maply.univ-lyon.fr Mohammed
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationChapter 4: Constrained estimators and tests in the multiple linear regression model (Part III)
Chapter 4: Constrained estimators and tests in the multiple linear regression model (Part III) Florian Pelgrin HEC September-December 2010 Florian Pelgrin (HEC) Constrained estimators September-December
More informationarxiv: v1 [math-ph] 16 Dec 2008
Some Remarks on Effective Range Formula in Potential Scattering arxiv:812.347v1 [math-ph] 16 Dec 28 Khosrow Chadan Laboratoire de Physique Théorique Université de Paris XI, Bâtiment 21, 9145 Orsay Cedex,
More informationCONVERGENCE OF POLYNOMIAL ERGODIC AVERAGES
CONVERGENCE OF POLYNOMIAL ERGODIC AVERAGES BERNARD HOST AND BRYNA KRA Abstract. We prove the L 2 convergence for an ergodic average of a product of functions evaluated along polynomial times in a totally
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationCONVERGENCE AND STABILITY CONSTANT OF THE THETA-METHOD
Conference Applications of Mathematics 2013 in honor of the 70th birthday of Karel Segeth. Jan Brandts, Sergey Korotov, Michal Křížek, Jakub Šístek, and Tomáš Vejchodský (Eds.), Institute of Mathematics
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 informationarxiv: v1 [math.ap] 28 Mar 2014
GROUNDSTATES OF NONLINEAR CHOQUARD EQUATIONS: HARDY-LITTLEWOOD-SOBOLEV CRITICAL EXPONENT VITALY MOROZ AND JEAN VAN SCHAFTINGEN arxiv:1403.7414v1 [math.ap] 28 Mar 2014 Abstract. We consider nonlinear Choquard
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationRelaxation methods and finite element schemes for the equations of visco-elastodynamics. Chiara Simeoni
Relaxation methods and finite element schemes for the equations of visco-elastodynamics Chiara Simeoni Department of Information Engineering, Computer Science and Mathematics University of L Aquila (Italy)
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationNumerical methods for non-standard fractional operators in the simulation of dielectric materials
Numerical methods for non-standard fractional operators in the simulation of dielectric materials Roberto Garrappa Department of Mathematics - University of Bari - Italy roberto.garrappa@uniba.it Fractional
More informationTopology of the set of singularities of a solution of the Hamilton-Jacobi Equation
Topology of the set of singularities of a solution of the Hamilton-Jacobi Equation Albert Fathi IAS Princeton March 15, 2016 In this lecture, a singularity for a locally Lipschitz real valued function
More informationHandout on Newton s Method for Systems
Handout on Newton s Method for Systems The following summarizes the main points of our class discussion of Newton s method for approximately solving a system of nonlinear equations F (x) = 0, F : IR n
More informationThe Milne error estimator for stiff problems
13 R. Tshelametse / SAJPAM. Volume 4 (2009) 13-28 The Milne error estimator for stiff problems Ronald Tshelametse Department of Mathematics University of Botswana Private Bag 0022 Gaborone, Botswana. E-mail
More informationGlobal Solutions for a Nonlinear Wave Equation with the p-laplacian Operator
Global Solutions for a Nonlinear Wave Equation with the p-laplacian Operator Hongjun Gao Institute of Applied Physics and Computational Mathematics 188 Beijing, China To Fu Ma Departamento de Matemática
More informationON FRACTIONAL ORDER CANCER MODEL
Journal of Fractional Calculus and Applications, Vol.. July, No., pp. 6. ISSN: 9-5858. http://www.fcaj.webs.com/ ON FRACTIONAL ORDER CANCER MODEL E. AHMED, A.H. HASHIS, F.A. RIHAN Abstract. In this work
More informationEXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM. Saeid Shokooh and Ghasem A. Afrouzi. 1. Introduction
MATEMATIČKI VESNIK MATEMATIQKI VESNIK 69 4 (217 271 28 December 217 research paper originalni nauqni rad EXISTENCE OF THREE WEAK SOLUTIONS FOR A QUASILINEAR DIRICHLET PROBLEM Saeid Shokooh and Ghasem A.
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationNonplanar Periodic Solutions for Restricted N+1-body Problems
Non-planar Periodic Solutions for Spatial Restricted N+1-body Problems Xiaoxiao Zhao, Fengying Li and Shiqing Zhang Sichuan University New Perspectives on the N-body Problem, BIRS, 2013 Construction Content
More information1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel
1462. Jacobi pseudo-spectral Galerkin method for second kind Volterra integro-differential equations with a weakly singular kernel Xiaoyong Zhang 1, Junlin Li 2 1 Shanghai Maritime University, Shanghai,
More informationA Systematic Approach to Extremum Seeking Based on Parameter Estimation
49th IEEE Conference on Decision and Control December 15-17, 21 Hilton Atlanta Hotel, Atlanta, GA, USA A Systematic Approach to Extremum Seeking Based on Parameter Estimation Dragan Nešić, Alireza Mohammadi
More informationAsymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface
Asymptotic Analysis of the Approximate Control for Parabolic Equations with Periodic Interface Patrizia Donato Université de Rouen International Workshop on Calculus of Variations and its Applications
More informationNonlinear Schrödinger problems: symmetries of some variational solutions
Nonlinear Differ. Equ. Appl. (3), 5 5 c Springer Basel AG -97/3/35- published online April 3, DOI.7/s3--3- Nonlinear Differential Equations and Applications NoDEA Nonlinear Schrödinger problems: symmetries
More informationA LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE. 1.
A LOWER BOUND ON BLOWUP RATES FOR THE 3D INCOMPRESSIBLE EULER EQUATION AND A SINGLE EXPONENTIAL BEALE-KATO-MAJDA ESTIMATE THOMAS CHEN AND NATAŠA PAVLOVIĆ Abstract. We prove a Beale-Kato-Majda criterion
More informationThe Journal of Nonlinear Sciences and Applications
J. Nonlinear Sci. Appl. 2 (2009), no. 3, 195 203 The Journal of Nonlinear Sciences Applications http://www.tjnsa.com ON MULTIPOINT ITERATIVE PROCESSES OF EFFICIENCY INDEX HIGHER THAN NEWTON S METHOD IOANNIS
More informationDerivation of relativistic Burgers equation on de Sitter background
Derivation of relativistic Burgers equation on de Sitter background BAVER OKUTMUSTUR Middle East Technical University Department of Mathematics 06800 Ankara TURKEY baver@metu.edu.tr TUBA CEYLAN Middle
More informationA Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems
A Hamiltonian Numerical Scheme for Large Scale Geophysical Fluid Systems Bob Peeters Joint work with Onno Bokhove & Jason Frank TW, University of Twente, Enschede CWI, Amsterdam PhD-TW colloquium, 9th
More informationExistence of at least two periodic solutions of the forced relativistic pendulum
Existence of at least two periodic solutions of the forced relativistic pendulum Cristian Bereanu Institute of Mathematics Simion Stoilow, Romanian Academy 21, Calea Griviţei, RO-172-Bucharest, Sector
More informationNUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS
TWMS J Pure Appl Math V5, N2, 24, pp22-228 NUMERICAL SOLUTION OF DELAY DIFFERENTIAL EQUATIONS VIA HAAR WAVELETS S ASADI, AH BORZABADI Abstract In this paper, Haar wavelet benefits are applied to the delay
More informationA Product Property of Sobolev Spaces with Application to Elliptic Estimates
Rend. Sem. Mat. Univ. Padova Manoscritto in corso di stampa pervenuto il 23 luglio 2012 accettato l 1 ottobre 2012 A Product Property of Sobolev Spaces with Application to Elliptic Estimates by Henry C.
More informationA Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals
arxiv:1503.00037v2 [math.na] 19 Mar 2015 A Posteriori Error Estimator for a Non-Standard Finite Difference Scheme Applied to BVPs on Infinite Intervals Riccardo Fazio and Alessandra Jannelli Department
More informationEvolution equations with spectral methods: the case of the wave equation
Evolution equations with spectral methods: the case of the wave equation Jerome.Novak@obspm.fr Laboratoire de l Univers et de ses Théories (LUTH) CNRS / Observatoire de Paris, France in collaboration with
More informationPeriodic solutions of weakly coupled superlinear systems
Periodic solutions of weakly coupled superlinear systems Alessandro Fonda and Andrea Sfecci Abstract By the use of a higher dimensional version of the Poincaré Birkhoff theorem, we are able to generalize
More informationOn Ekeland s variational principle
J. Fixed Point Theory Appl. 10 (2011) 191 195 DOI 10.1007/s11784-011-0048-x Published online March 31, 2011 Springer Basel AG 2011 Journal of Fixed Point Theory and Applications On Ekeland s variational
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 informationGeneralised Summation-by-Parts Operators and Variable Coefficients
Institute Computational Mathematics Generalised Summation-by-Parts Operators and Variable Coefficients arxiv:1705.10541v [math.na] 16 Feb 018 Hendrik Ranocha 14th November 017 High-order methods for conservation
More informationLegendre modified moments for Euler s constant
Journal of Computational and Applied Mathematics 29 (28) 484 492 www.elsevier.com/locate/cam Legendre modified moments for Euler s constant Marc Prévost Laboratoire de Mathématiques Pures et Appliquées,
More informationMath 2250 Final Exam Practice Problem Solutions. f(x) = ln x x. 1 x. lim. lim. x x = lim. = lim 2
Math 5 Final Eam Practice Problem Solutions. What are the domain and range of the function f() = ln? Answer: is only defined for, and ln is only defined for >. Hence, the domain of the function is >. Notice
More informationThe effects of a discontinues weight for a problem with a critical nonlinearity
arxiv:1405.7734v1 [math.ap] 9 May 014 The effects of a discontinues weight for a problem with a critical nonlinearity Rejeb Hadiji and Habib Yazidi Abstract { We study the minimizing problem px) u dx,
More informationOn the Euler Lagrange Equation in Calculus of Variations
On the Euler Lagrange Equation in Calculus of Variations Ivar Ekeland Vietnam Journal of Mathematics ISSN 235-221X Vietnam J. Math. DOI 1.17/s113-18-285-z 1 23 Your article is protected by copyright and
More informationETNA Kent State University
Electronic Transactions on Numerical Analysis. Volume 4, pp. 373-38, 23. Copyright 23,. ISSN 68-963. ETNA A NOTE ON THE RELATION BETWEEN THE NEWTON HOMOTOPY METHOD AND THE DAMPED NEWTON METHOD XUPING ZHANG
More informationComparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations
Comparing Numerical Methods for Solving Nonlinear Fractional Order Differential Equations Farhad Farokhi, Mohammad Haeri, and Mohammad Saleh Tavazoei Abstract This paper is a result of comparison of some
More informationA FV Scheme for Maxwell s equations
A FV Scheme for Maxwell s equations Convergence Analysis on unstructured meshes Stephanie Lohrengel * Malika Remaki ** *Laboratoire J.A. Dieudonné (UMR CNRS 6621), Université de Nice Sophia Antipolis,
More informationThe collocation method for ODEs: an introduction
058065 - Collocation Methods for Volterra Integral Related Functional Differential The collocation method for ODEs: an introduction A collocation solution u h to a functional equation for example an ordinary
More informationOptimal stopping time formulation of adaptive image filtering
Optimal stopping time formulation of adaptive image filtering I. Capuzzo Dolcetta, R. Ferretti 19.04.2000 Abstract This paper presents an approach to image filtering based on an optimal stopping time problem
More informationThe Rademacher Cotype of Operators from l N
The Rademacher Cotype of Operators from l N SJ Montgomery-Smith Department of Mathematics, University of Missouri, Columbia, MO 65 M Talagrand Department of Mathematics, The Ohio State University, 3 W
More informationExact multiplicity of boundary blow-up solutions for a bistable problem
Computers and Mathematics with Applications 54 (2007) 1285 1292 www.elsevier.com/locate/camwa Exact multiplicity of boundary blow-up solutions for a bistable problem Junping Shi a,b,, Shin-Hwa Wang c a
More informationA Fixed Point Theorem for G-Monotone Multivalued Mapping with Application to Nonlinear Integral Equations
Filomat 31:7 (217), 245 252 DOI 1.2298/FIL17745K Published by Faculty of Sciences and Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat A Fixed Point Theorem for G-Monotone
More informationA collocation method for solving the fractional calculus of variation problems
Bol. Soc. Paran. Mat. (3s.) v. 35 1 (2017): 163 172. c SPM ISSN-2175-1188 on line ISSN-00378712 in press SPM: www.spm.uem.br/bspm doi:10.5269/bspm.v35i1.26333 A collocation method for solving the fractional
More informationError estimates for the Raviart-Thomas interpolation under the maximum angle condition
Error estimates for the Raviart-Thomas interpolation under the maximum angle condition Ricardo G. Durán and Ariel L. Lombardi Abstract. The classical error analysis for the Raviart-Thomas interpolation
More informationarxiv:math/ v1 [math.ap] 28 Oct 2005
arxiv:math/050643v [math.ap] 28 Oct 2005 A remark on asymptotic completeness for the critical nonlinear Klein-Gordon equation Hans Lindblad and Avy Soffer University of California at San Diego and Rutgers
More informationHe s Homotopy Perturbation Method for Nonlinear Ferdholm Integro-Differential Equations Of Fractional Order
H Saeedi, F Samimi / International Journal of Engineering Research and Applications (IJERA) ISSN: 2248-9622 wwwijeracom Vol 2, Issue 5, September- October 22, pp52-56 He s Homotopy Perturbation Method
More informationMAXIMALITY OF SUMS OF TWO MAXIMAL MONOTONE OPERATORS
MAXIMALITY OF SUMS OF TWO MAXIMAL MONOTONE OPERATORS JONATHAN M. BORWEIN, FRSC Abstract. We use methods from convex analysis convex, relying on an ingenious function of Simon Fitzpatrick, to prove maximality
More informationA Simple Compact Fourth-Order Poisson Solver on Polar Geometry
Journal of Computational Physics 182, 337 345 (2002) doi:10.1006/jcph.2002.7172 A Simple Compact Fourth-Order Poisson Solver on Polar Geometry Ming-Chih Lai Department of Applied Mathematics, National
More informationOptimized Schwarz methods in spherical geometry with an overset grid system
Optimized Schwarz methods in spherical geometry with an overset grid system J. Côté 1, M. J. Gander 2, L. Laayouni 3, and A. Qaddouri 4 1 Recherche en prévision numérique, Meteorological Service of Canada,
More informationIn this paper we study periodic solutions of a second order differential equation
Topological Methods in Nonlinear Analysis Journal of the Juliusz Schauder Center Volume 5, 1995, 385 396 A GRANAS TYPE APPROACH TO SOME CONTINUATION THEOREMS AND PERIODIC BOUNDARY VALUE PROBLEMS WITH IMPULSES
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More informationPERTURBATION OF THE FREE BOUNDARY IN ELLIPTIC PROBLEM WITH DISCONTINUITIES
Electronic Journal of Differential Equations, Vol. 16 (16), No. 13, pp. 1 14. IN: 17-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu PERTURBATION OF THE FREE BOUNDARY IN ELLIPTIC PROBLEM
More informationObservability for deterministic systems and high-gain observers
Observability for deterministic systems and high-gain observers design. Part 1. March 29, 2011 Introduction and problem description Definition of observability Consequences of instantaneous observability
More informationDate of Birth 18 July, CNIC Number Contact Number Pakistan
Salman AminMALIK Personal Information Date of Birth 18 July, 1979 Place of Birth Lahore, Pakistan Nationality Pakistani CNIC Number 35202-6719738-5 Gender Male Marital Status Married Contact Number Pakistan
More information