THE NONLINEAR DIFFUSION EQUATION

Transcription

1 THE NONLINEAR DIFFUSION EQUATION

2 THE NONLINEAR DIFFUSION EQUATION Asymptotic Solutions and Statistica! Problems by J. M. BURGERS Institute for Fluid Dynamics and Applied Mathematics, University of Maryland, College Park, Md., U.S.A. SPRINGER-SCIENCE+BUSINESS MEDIA, B.V.

3 First published under the title Statistica/ Problems Connected with Asymptotic Solutions of the One-Dimensional Non/inear Diffusion Equation in 1973 by the Institute for Fluid Dynamics and Applied Mathematics, University of Maryland, Lecture Series No. 52 Library ofcongress Catalog Card Number ISBN ISBN (ebook) DOI / Ali Rights Reserved Copyright 1974 by Springer Science+Business Media Dordrecht Originally published by D. Reidel Publishing Company, Dordrecht, Holland in Softcover reprint of the hardcover 1 st Edition 1974 No part of this book may be reproduced in any form, by print, photoprint, microfilm, or any other means, without written permission from the publisher

4 TABLE OF CONTENTS PREFACE IX INTRODUCTION CHAPTER I 1 THE HOPF-COLE SOLUTION OF THE NONLINEAR DIFFUSION EQUATION AND ITS GEOMETRICAL INTERPRETA TION FOR THE CASE OF SMALL DIFFUSIVITY 9 1. Basic Solution 9 2. Geometric Interpretation of the Solution. - Possibility of Multiple Roots Steep Fronts or 'Shock Waves'. General Form of the Solution for Large Values of t Alternative Derivation of the Approximate Solution of Equation (1.1) in the Neighborhood of a Steep Front. Additional Observations 18 CHAPTER II 1 DIGRESSION ON GENERALIZATIONS OF THE GEOMETRIC METHOD OF SOLUTION. - SOLUTIONS OF EQUATION (1.1) FOR THE DOMAIN x>o WITH A BOUNDARY CONDITIONAT x=o Application ofthe Geometric Method ofsolution to Slightly More General Equations Solutions of the Nonlinear Diffusion Equation with a Boundary Conditionat x=o. Two Preliminary Cases Generalization of Solution (6.6) Extension to a Continuous Function 31 CHAPTER III 1 STA TISTICAL PROBLEMS CONNECTED WITH THE SOL U- TIONS OF CHAPTER I, FOR V-> +0 AND t Ensemble of Initial Data-Curves. The Parameter J Statistics of a Chain of Parabolic Arcs Transformation of the Statistica! Problem Additional Observations 44 CHAPTER IV 1 SOLUTIONS OF THE LINEAR DIFFUSION EQUATION WITH A BOUNDARY CONDITION REFERRING TO A PARABOLA 13. Basic Functions 46 46

5 VI TABLE OF CONTENTS 14. Application of Green's Theorem Solution of Equation (13.1) Assuming Prescribed Values on the Boundary Change of the Form of the Boundary Transformation to a Coordinate System with Inclined x-axis Extension ofequation (16.2} Application of the Results of Sections 16 and 17 to a Parabolic Boundary Integro-Differential Equation for 8(z) Series Development for 8(z) The Function E(x) Application to a Parabolic Boundary Effect of a Change of the Boundary upon E(x) The Integral P The Normalization Condition (10.10) A Theorem Connected with the Integral P 1 for a Parabolic Boundary 70 CHAPTER V 1 DEVELOPMENT OF THE FUNCTIONS '1', E, F IN TERMS OF EXPONENTIALS MUL TIPLIED BY BESSEL FUNCTIONS The Transformation of the Differential Equation and its Solution Properties ofthe Functions q 0 (v), qn(v) The Green's Function Connected with (28.6} Summation Formulas Series for the Functions 1/J, '1', e, cp Expression for E(x) 82 CHAPTER VI 1 EVAL UATION OF INTEGRALS AND SUMS DEPENDING ON THE FUNCTIONS '1', E, F Basic Definitions Initial Integrations and Summations Series for the Functions Rm and /m More Summation Formulas Continuation from Section Results Obtained for the Integrals P 0, P 1,... P* Results Obtained for Some Additional Integrals 108 List oflntegrals 110 APPENDIX TO CHAPTER VI Additional Summation F ormulas 112

6 TABLE OF CONTENTS VII CHAPTER VII 1 MEAN VALUES CONNECTED WITH THE SA WTOOTH CURVE OF FIGURE Recapitulation ofprevious Results Mean Values Referring to a Single Shock Other Mean Va1ues Which can be Expressed in a Simple Way The Distribution Function for the Arc Lengths ~k 130 CHAPTER VIII 1 DISTRIBUTION FUNCTIONS REFERRING TO SETS OF TWO CONSECUTIVE ARCS Integrals with Two Successive 'l'-functions Weighting Function for the Wave1engths Ak Alternative Expressions for the Functions K* and L* Relations Invo1ving the Quantity clj The Value of clj Evaluation of< ~k)* for A-+ O Calcu1ation of< ~k~k+ 1 )* for A-+ O Calculation of <a~k+ 1)* for A-+ O Application of Expressions (48.5) and (48.6) for K* and L* 150 CHAPTER IX 1 CORRELA TION FUNCTIONS AND DISTRIBUTION FUNC- TIONS REFERRING TO SETS OF MORE THAN TWO CONSECUTIVE ARCS 55. Calculation of u 1 u Calculation of uiu Behavior for Small Values of z Behavior of the Sums Occurring in Equation (55.18) for Indefinitely Increasing Values of z Behavior of the Sums Occurring in Equation (56.10) or (56.12) for Indefinitely Increasing Va1ues of z Integrals Referring to a Series of Consecutive Arcs Continuation The Integra1s for the Functions W 0 ;(z),..., W 3 ;(z), lntroduced in Sections 55 and LIST OF PAPERS BY THE AUTHOR, CONNECTED WITH THE PRESENT MONOGRAPH 174

7 PREFACE Since the 'Introduction' to the main text gives an account of the way in which the problems treated in the following pages originated, this 'Preface' may be limited to an acknowledgement of the support the work has received. It started during the period when I was professor of aero- and hydrodynamics at the Technical University in Delft, Netherlands, and many discussions with colleagues ha ve in:fluenced its development. Oftheir names I mention here only that ofh. A. Kramers. Papers No ofthe list given at the end ofthe text were written during that period. Severa! ofthese were attempts to explore ideas which later had to be abandoned, but gradually a line of thought emerged which promised more definite results. This line began to come to the foreground in pa per No. 3 (1939}, while a preliminary formulation ofthe results was given in paper No. 12 (1954}. At that time, however, there still was missing a practica! method for manipulating a certain distribution function of central interest. A six months stay at the Hydrodynamics Laboratories ofthe California Institute of Technology, Pasadena, California ( }, was supported by a Contract with the Department of the Air F orce, N o. AF 33(038} A course of lectures was given during this period, which were published in typescript under the title 'On Turbulent Fluid Motion', as Report No. E-34.1, July 1951, of the Hydrodynamics Laboratory. Main progress carne after my return to Delft, when the Hopf-Cole transformation of the nonlinear equation into a linear one had opened a new way of attack. Since the end of 1955 until1971 the work was carried on at the Institute for Fluid Dynamics and Applied Mathematics ofthe University ofmaryland.lt was partially supported by the Air Force Oftice of Scientific Research, as a subject coming under Contracts AF-18(600}993, AF-49(638}401, F C-0023, and Grants AFOSR , AFOSR , AFOSR , and AFOSR For ali support and for the general interest of many colleagues, deep gratitude is expressed. At one stage ofthe work it was surmised that the integrals considered in Chapter VI below might be related to Wiener integrals. This was discussed in the beginning of 1963 with Dr F. Pollaczek, at that time a visitor at the Institute for Fluid Dynamics and Applied Mathematics; and also with Dr M. Kac of the Rockefeller University in New York. Dr Kac pointed out that an important property of Wiener integrals is that they lead to a partial differential equation, a subject which is treated in his book, Probability and Related Topics in Physical Sciences, Interscience Publishers, New York, 1959, Chapter IV, pp (in particular p. 171}.1t appeared that the corresponding differential equation for my case [Equation (28.3} below] could be obtained

8 X PREFACE by a direct transformation ofthe original equation, so that there was no need to give further attention to Wiener integrals. This is developed in Chapters V and VI of the monograph. It furnished the hasis for the evaluation of ali integrals connected with the central distribution function. The results of the calculations ha ve led to definite numbers, without adjustable constants. Numerica! calculations sometimes were desirable in order to obtain certain transcendental numbers. Mr Y.-K. Hsu and Mr S.-K. Oh gave help with the calcu Iation ofthe coefficients ofthe series (21.3) for the function <f; (z) and for the evaluation of E(O) according to Equation (23.7).- Mr David Wilson made calculations for the numbers m 1 (Equation (38.17a)) and m 3 (Equation (38.17b)). Special gratitude is due to the secretaries who have typed out, often repeatedly, the extensive manuscript notes with their long and complicated equations. I mention in particular Miss Helga Lincke for her work in ; and Mrs Irene Huvos for the still more extensive!ater work and her never ending patience, during the years To Professor Bruce Kellogg of the Institute I am indebted for much he1p and care given when the Institute decided to distribute the typewritten text as No. 52 of the Institute's 'Lecture Series', in November I consider it a great honor that D. Reidel Publishing Company in Dordrecht, The Netherlands, has decided to publish this monograph in the form of a book. I am well aware that the method followed here in treating asymptotic solutions of the non1inear diffusion equation and their statistics is a specialized one. Many approaches have been explored by a variety of authors starting from other points of view. It is impossib1e to extend the present monograph in such a way that an adequate digest of the existing literature can be given. This has therefore been left aside. A survey of published so1utions until 1972 has been presented by E. R. Benton and G. W. Platzmann; see the reference in footnote* at the beginning of Chapter 1. New papers are continually forthcoming. A final remark: whi1e great care has been taken to check the numerous equations by comparing them with the original material, a full guarantee against misprints cannot be given and the clemency of the reader is invoked. J.M.B.