scispace - formally typeset
Search or ask a question

Showing papers on "Partial differential equation published in 2012"


Journal ArticleDOI
TL;DR: First the freefem++ software deals with mesh adaptation for problems in two and three dimension, second, it solves numerically a problem with phase change and natural convection, and finally to show the possibilities for HPC the software solves a Laplace equation by a Schwarz domain decomposition problem on parallel computer.
Abstract: This is a short presentation of the freefem++ software. In Section 1, we recall most of the characteristics of the software, In Section 2, we recall how to to build the weak form of a partial differential equation (PDE) from the strong form. In the 3 last sections, we present different examples and tools to illustrated the power of the software. First we deal with mesh adaptation for problems in two and three dimension, second, we solve numerically a problem with phase change and natural convection, and the finally to show the possibilities for HPC we solve a Laplace equation by a Schwarz domain decomposition problem on parallel computer.

2,867 citations


Book
07 Feb 2012
TL;DR: In this paper, the Fokker-Planck Equation for N Variables (FPE) was extended to N = 1 variable and N = 2 variables, where N is the number of variables in the system.
Abstract: 1. Introduction.- 1.1 Brownian Motion.- 1.1.1 Deterministic Differential Equation.- 1.1.2 Stochastic Differential Equation.- 1.1.3 Equation of Motion for the Distribution Function.- 1.2 Fokker-Planck Equation.- 1.2.1 Fokker-Planck Equation for One Variable.- 1.2.2 Fokker-Planck Equation for N Variables.- 1.2.3 How Does a Fokker-Planck Equation Arise?.- 1.2.4 Purpose of the Fokker-Planck Equation.- 1.2.5 Solutions of the Fokker-Planck Equation.- 1.2.6 Kramers and Smoluchowski Equations.- 1.2.7 Generalizations of the Fokker-Planck Equation.- 1.3 Boltzmann Equation.- 1.4 Master Equation.- 2. Probability Theory.- 2.1 Random Variable and Probability Density.- 2.2 Characteristic Function and Cumulants.- 2.3 Generalization to Several Random Variables.- 2.3.1 Conditional Probability Density.- 2.3.2 Cross Correlation.- 2.3.3 Gaussian Distribution.- 2.4 Time-Dependent Random Variables.- 2.4.1 Classification of Stochastic Processes.- 2.4.2 Chapman-Kolmogorov Equation.- 2.4.3 Wiener-Khintchine Theorem.- 2.5 Several Time-Dependent Random Variables.- 3. Langevin Equations.- 3.1 Langevin Equation for Brownian Motion.- 3.1.1 Mean-Squared Displacement.- 3.1.2 Three-Dimensional Case.- 3.1.3 Calculation of the Stationary Velocity Distribution Function.- 3.2 Ornstein-Uhlenbeck Process.- 3.2.1 Calculation of Moments.- 3.2.2 Correlation Function.- 3.2.3 Solution by Fourier Transformation.- 3.3 Nonlinear Langevin Equation, One Variable.- 3.3.1 Example.- 3.3.2 Kramers-Moyal Expansion Coefficients.- 3.3.3 Ito's and Stratonovich's Definitions.- 3.4 Nonlinear Langevin Equations, Several Variables.- 3.4.1 Determination of the Langevin Equation from Drift and Diffusion Coefficients.- 3.4.2 Transformation of Variables.- 3.4.3 How to Obtain Drift and Diffusion Coefficients for Systems.- 3.5 Markov Property.- 3.6 Solutions of the Langevin Equation by Computer Simulation.- 4. Fokker-Planck Equation.- 4.1 Kramers-Moyal Forward Expansion.- 4.1.1 Formal Solution.- 4.2 Kramers-Moyal Backward Expansion.- 4.2.1 Formal Solution.- 4.2.2 Equivalence of the Solutions of the Forward and Backward Equations.- 4.3 Pawula Theorem.- 4.4 Fokker-Planck Equation for One Variable.- 4.4.1 Transition Probability Density for Small Times.- 4.4.2 Path Integral Solutions.- 4.5 Generation and Recombination Processes.- 4.6 Application of Truncated Kramers-Moyal Expansions.- 4.7 Fokker-Planck Equation for N Variables.- 4.7.1 Probability Current.- 4.7.2 Joint Probability Distribution.- 4.7.3 Transition Probability Density for Small Times.- 4.8 Examples for Fokker-Planck Equations with Several Variables.- 4.8.1 Three-Dimensional Brownian Motion without Position Variable.- 4.8.2 One-Dimensional Brownian Motion in a Potential.- 4.8.3 Three-Dimensional Brownian Motion in an External Force.- 4.8.4 Brownian Motion of Two Interacting Particles in an External Potential.- 4.9 Transformation of Variables.- 4.10 Covariant Form of the Fokker-Planck Equation.- 5. Fokker-Planck Equation for One Variable Methods of Solution.- 5.1 Normalization.- 5.2 Stationary Solution.- 5.3 Ornstein-Uhlenbeck Process.- 5.4 Eigenfunction Expansion.- 5.5 Examples.- 5.5.1 Parabolic Potential.- 5.5.2 Inverted Parabolic Potential.- 5.5.3 Infinite Square Well for the Schrudinger Potential.- 5.5.4 V-Shaped Potential for the Fokker-Planck Equation.- 5.6 Jump Conditions.- 5.7 A Bistable Model Potential.- 5.8 Eigenfunctions and Eigenvalues of Inverted Potentials.- 5.9 Approximate and Numerical Methods for Determining Eigenvalues and Eigenfunctions.- 5.9.1 Variational Method.- 5.9.2 Numerical Integration.- 5.9.3 Expansion into a Complete Set.- 5.10 Diffusion Over a Barrier.- 5.10.1 Kramers' Escape Rate.- 5.10.2 Bistable and Metastable Potential.- 6. Fokker-Planck Equation for Several Variables Methods of Solution.- 6.1 Approach of the Solutions to a Limit Solution.- 6.2 Expansion into a Biorthogonal Set.- 6.3 Transformation of the Fokker-Planck Operator, Eigenfunction Expansions.- 6.4 Detailed Balance.- 6.5 Ornstein-Uhlenbeck Process.- 6.6 Further Methods for Solving the Fokker-Planck Equation.- 6.6.1 Transformation of Variables.- 6.6.2 Variational Method.- 6.6.3 Reduction to an Hermitian Problem.- 6.6.4 Numerical Integration.- 6.6.5 Expansion into Complete Sets.- 6.6.6 Matrix Continued-Fraction Method.- 6.6.7 WKB Method.- 7. Linear Response and Correlation Functions.- 7.1 Linear Response Function.- 7.2 Correlation Functions.- 7.3 Susceptibility.- 8. Reduction of the Number of Variables.- 8.1 First-Passage Time Problems.- 8.2 Drift and Diffusion Coefficients Independent of Some Variables.- 8.2.1 Time Integrals of Markovian Variables.- 8.3 Adiabatic Elimination of Fast Variables.- 8.3.1 Linear Process with Respect to the Fast Variable.- 8.3.2 Connection to the Nakajima-Zwanzig Projector Formalism.- 9. Solutions of Tridiagonal Recurrence Relations, Application to Ordinary and Partial Differential Equations.- 9.1 Applications and Forms of Tridiagonal Recurrence Relations.- 9.1.1 Scalar Recurrence Relation.- 9.1.2 Vector Recurrence Relation.- 9.2 Solutions of Scalar Recurrence Relations.- 9.2.1 Stationary Solution.- 9.2.2 Initial Value Problem.- 9.2.3 Eigenvalue Problem.- 9.3 Solutions of Vector Recurrence Relations.- 9.3.1 Initial Value Problem.- 9.3.2 Eigenvalue Problem.- 9.4 Ordinary and Partial Differential Equations with Multiplicative Harmonic Time-Dependent Parameters.- 9.4.1 Ordinary Differential Equations.- 9.4.2 Partial Differential Equations.- 9.5 Methods for Calculating Continued Fractions.- 9.5.1 Ordinary Continued Fractions.- 9.5.2 Matrix Continued Fractions.- 10. Solutions of the Kramers Equation.- 10.1 Forms of the Kramers Equation.- 10.1.1 Normalization of Variables.- 10.1.2 Reversible and Irreversible Operators.- 10.1.3 Transformation of the Operators.- 10.1.4 Expansion into Hermite Functions.- 10.2 Solutions for a Linear Force.- 10.2.1 Transition Probability.- 10.2.2 Eigenvalues and Eigenfunctions.- 10.3 Matrix Continued-Fraction Solutions of the Kramers Equation.- 10.3.1 Initial Value Problem.- 10.3.2 Eigenvalue Problem.- 10.4 Inverse Friction Expansion.- 10.4.1 Inverse Friction Expansion for K0(t), G0,0(t) and L0(t).- 10.4.2 Determination of Eigenvalues and Eigenvectors.- 10.4.3 Expansion for the Green's Function Gn,m(t).- 10.4.4 Position-Dependent Friction.- 11. Brownian Motion in Periodic Potentials.- 11.1 Applications.- 11.1.1 Pendulum.- 11.1.2 Superionic Conductor.- 11.1.3 Josephson Tunneling Junction.- 11.1.4 Rotation of Dipoles in a Constant Field.- 11.1.5 Phase-Locked Loop.- 11.1.6 Connection to the Sine-Gordon Equation.- 11.2 Normalization of the Langevin and Fokker-Planck Equations.- 11.3 High-Friction Limit.- 11.3.1 Stationary Solution.- 11.3.2 Time-Dependent Solution.- 11.4 Low-Friction Limit.- 11.4.1 Transformation to E and x Variables.- 11.4.2 'Ansatz' for the Stationary Distribution Functions.- 11.4.3 x-Independent Functions.- 11.4.4 x-Dependent Functions.- 11.4.5 Corrected x-Independent Functions and Mobility.- 11.5 Stationary Solutions for Arbitrary Friction.- 11.5.1 Periodicity of the Stationary Distribution Function.- 11.5.2 Matrix Continued-Fraction Method.- 11.5.3 Calculation of the Stationary Distribution Function.- 11.5.4 Alternative Matrix Continued Fraction for the Cosine Potential.- 11.6 Bistability between Running and Locked Solution.- 11.6.1 Solutions Without Noise.- 11.6.2 Solutions With Noise.- 11.6.3 Low-Friction Mobility With Noise.- 11.7 Instationary Solutions.- 11.7.1 Diffusion Constant.- 11.7.2 Transition Probability for Large Times.- 11.8 Susceptibilities.- 11.8.1 Zero-Friction Limit.- 11.9 Eigenvalues and Eigenfunctions.- 11.9.1 Eigenvalues and Eigenfunctions in the Low-Friction Limit.- 12. Statistical Properties of Laser Light.- 12.1 Semiclassical Laser Equations.- 12.1.1 Equations Without Noise.- 12.1.2 Langevin Equation.- 12.1.3 Laser Fokker-Planck Equation.- 12.2 Stationary Solution and Its Expectation Values.- 12.3 Expansion in Eigenmodes.- 12.4 Expansion into a Complete Set Solution by Matrix Continued Fractions.- 12.4.1 Determination of Eigenvalues.- 12.5 Transient Solution.- 12.5.1 Eigenfunction Method.- 12.5.2 Expansion into a Complete Set.- 12.5.3 Solution for Large Pump Parameters.- 12.6 Photoelectron Counting Distribution.- 12.6.1 Counting Distribution for Short Intervals.- 12.6.2 Expectation Values for Arbitrary Intervals.- Appendices.- A1 Stochastic Differential Equations with Colored Gaussian Noise.- A2 Boltzmann Equation with BGK and SW Collision Operators.- A3 Evaluation of a Matrix Continued Fraction for the Harmonic Oscillator.- A4 Damped Quantum-Mechanical Harmonic Oscillator.- A5 Alternative Derivation of the Fokker-Planck Equation.- A6 Fluctuating Control Parameter.- S. Supplement to the Second Edition.- S.1 Solutions of the Fokker-Planck Equation by Computer Simulation (Sect. 3.6).- S.2 Kramers-Moyal Expansion (Sect. 4.6).- S.3 Example for the Covariant Form of the Fokker-Planck Equation (Sect. 4.10).- S.4 Connection to Supersymmetry and Exact Solutions of the One Variable Fokker-Planck Equation (Chap. 5).- S.5 Nondifferentiability of the Potential for the Weak Noise Expansion (Sects. 6.6 and 6.7).- S.6 Further Applications of Matrix Continued-Fractions (Chap. 9).- S.7 Brownian Motion in a Double-Well Potential (Chaps. 10 and 11).- S.8 Boundary Layer Theory (Sect. 11.4).- S.9 Calculation of Correlation Times (Sect. 7.12).- S.10 Colored Noise (Appendix A1).- S.11 Fokker-Planck Equation with a Non-Positive-Definite Diffusion Matrix and Fokker-Planck Equation with Additional Third-Order-Derivative Terms.- References.

2,582 citations


Book
01 Mar 2012
TL;DR: The relationship between symmetries of a linear second-order partial differential equation of mathematical physics, the coordinate systems in which the equation admits solutions via separation of variables, and the properties of the special functions that arise in this manner are discussed in this paper.
Abstract: Originally published in 1977, this volume is concerned with the relationship between symmetries of a linear second-order partial differential equation of mathematical physics, the coordinate systems in which the equation admits solutions via separation of variables, and the properties of the special functions that arise in this manner. Some group-theoretic twists in the ancient method of separation of variables that can be used to provide a foundation for much of special function theory are shown. In particular, it is shown explicitly that all special functions that arise via separation of variables in the equations of mathematical physics can be studied using group theory.

800 citations


Journal ArticleDOI
TL;DR: It is shown that expansion in DMD modes is unique under certain conditions, and an “optimized” DMD is introduced that computes an arbitrary number of dynamical modes from a data set and is superior at calculating physically relevant frequencies, and is less numerically sensitive.
Abstract: Dynamic mode decomposition (DMD) is an Arnoldi-like method based on the Koopman operator. It analyzes empirical data, typically generated by nonlinear dynamics, and computes eigenvalues and eigenmodes of an approximate linear model. Without explicit knowledge of the dynamical operator, it extracts frequencies, growth rates, and spatial structures for each mode. We show that expansion in DMD modes is unique under certain conditions. When constructing mode-based reduced-order models of partial differential equations, subtracting a mean from the data set is typically necessary to satisfy boundary conditions. Subtracting the mean of the data exactly reduces DMD to the temporal discrete Fourier transform (DFT); this is restrictive and generally undesirable. On the other hand, subtracting an equilibrium point generally preserves the DMD spectrum and modes. Next, we introduce an “optimized” DMD that computes an arbitrary number of dynamical modes from a data set. Compared to DMD, optimized DMD is superior at calculating physically relevant frequencies, and is less numerically sensitive. We test these decomposition methods on data from a two-dimensional cylinder fluid flow at a Reynolds number of 60. Time-varying modes computed from the DMD variants yield low projection errors.

672 citations


Book
06 Nov 2012
TL;DR: The Navier-Stokes equations as mentioned in this paper are the basic model for the flow of incompressible viscous fluids, and they are used in many other kinds of nonlinear evolution problems.
Abstract: The objective of this self-contained book is two-fold. First, the reader is introduced to the modelling and mathematical analysis used in fluid mechanics, especially concerning the Navier-Stokes equations which is the basic model for the flow of incompressible viscous fluids. Authors introduce mathematical tools so that the reader is able to use them for studying many other kinds of partial differential equations, in particular nonlinear evolution problems. The background needed are basic results in calculus, integration, and functional analysis. Some sections certainly contain more advanced topics than others. Nevertheless, the authors' aim is that graduate or PhD students, as well as researchers who are not specialized in nonlinear analysis or in mathematical fluid mechanics, can find a detailed introduction to this subject.

518 citations


Book
16 Jun 2012
TL;DR: In this article, the authors present a framework for stochastic ODEs driven by Brownian White Noise processes, and apply it to Stochastic Ordinary Differential Equations.
Abstract: Framework.- Applications to Stochastic Ordinary Differential Equations.- Stochastic Partial Differential Equations Driven by Brownian White Noise.- Stochastic Partial Differential Equations Driven by L#x00E9 vy Processes.

455 citations


Book ChapterDOI
01 Jan 2012
TL;DR: In this paper, a set of short tables of integral transforms of the functions that are either cited in the text or are in most common use in mathematical, physical, and engineering applications are provided.
Abstract: In this chapter, we provide a set of short tables of integral transforms of the functions that are either cited in the text or are in most common use in mathematical, physical, and engineering applications. For exhaustive lists of integral transforms, the reader is referred to Erdelyi et al. (Tables of Integral Transforms, Vols. 1 and 2, 1954), Campbell and Foster (Fourier Integrals for Practical Applications, 1948), Ditkin and Prudnikov (Integral Transforms and Operational Calculus, 1965), Doetsch (Introduction to the Theory and Applications of the Laplace Transformation, 1970), Marichev (1983), Debnath (1995), Debnath and Bhatta (Integral Transforms and Their Applications, 2nd edition, 2007), Oberhettinger (Tables of Bessel Transforms, 1972).

377 citations


Journal ArticleDOI
TL;DR: In this article, a simplified electrochemical and thermal model of LiFePO4-graphite based Li-ion batteries is developed based on the Pseudo Two-Dimensional mathematical structure.
Abstract: Derived from the Pseudo Two-Dimensional mathematical structure, a simplified electrochemical and thermal model of LiFePO4-graphite based Li-ion batteries is developed in this paper. Embedding the porous electrode theory, this model integrates the main design parameters of Li-ion systems and its partial differential equations mathematical structure makes it a promising candidate for battery management system (BMS) applications and comprehensive aging investigations. Based on a modified Single-Particle approach, the model is used to simulate and discuss capacity restitution in galvanostatic charges and discharges at various rates and temperatures. Constant high-rate solicitations similar to fast charge of plug-in electric vehicles or electric vehicles, are experimentally tested and simulated with the present model. Also, thermal issues occurring during these specific operating conditions are quantitatively pointed out. The concept of current-dependent spherical particle radius is used to obtain good agreement with experimental data related to galvanostatic charges and discharges. The capabilities and limits of this preliminary modeling work are discussed in detail and ways to extend the potentialities of this approach to BMS applications are proposed.

279 citations


Journal ArticleDOI
Bin Lu1
TL;DR: In this article, the fractional derivatives in the sense of modified Riemann-Liouville derivative and the first integral method are employed for constructing the exact solutions of nonlinear time-fractional partial differential equations.

276 citations


Journal ArticleDOI
TL;DR: This work gives a systematic method for discretizing Hamiltonian partial differential equations (PDEs) with constant symplectic structure, while preserving their energy exactly.

268 citations


Journal ArticleDOI
TL;DR: In this article, the phase field partial differential equations (PDEs) are solved simultaneously together with PDEs describing additional physics, such as solid mechanics and heat conduction, using the Jacobian-free Newton Krylov Method.

Journal ArticleDOI
TL;DR: In this article, the nonlinear Fourier transform (NFT) is used for solving integrable partial differential equations governing wave propagation in certain nonlinear media, such as optical fibers.
Abstract: The nonlinear Fourier transform (NFT), a powerful tool in soliton theory and exactly solvable models, is a method for solving integrable partial differential equations governing wave propagation in certain nonlinear media. The NFT decorrelates signal degrees-of-freedom in such models, in much the same way that the Fourier transform does for linear systems. In this three-part series of papers, this observation is exploited for data transmission over integrable channels such as optical fibers, where pulse propagation is governed by the nonlinear Schr\"odinger equation. In this transmission scheme, which can be viewed as a nonlinear analogue of orthogonal frequency-division multiplexing commonly used in linear channels, information is encoded in the nonlinear frequencies and their spectral amplitudes. Unlike most other fiber-optic transmission schemes, this technique deals with both dispersion and nonlinearity directly and unconditionally without the need for dispersion or nonlinearity compensation methods. This first paper explains the mathematical tools that underlie the method.

Journal ArticleDOI
TL;DR: In this paper, a survey of local and global Carleman estimates for elliptic and parabolic partial differential equations is presented, and the connexion of these optimality results to the local phase-space geometry after conjugation with the weight function is pointed out.
Abstract: Local and global Carleman estimates play a central role in the study of some partial differential equations regarding questions such as unique continuation and controllability. We survey and prove such estimates in the case of elliptic and parabolic operators by means of semi-classical microlocal techniques. Optimality results for these estimates and some of their consequences are presented. We point out the connexion of these optimality results to the local phase-space geometry after conjugation with the weight function. Firstly, we introduce local Carleman estimates for elliptic operators and deduce unique continuation properties as well as interpolation inequalities. These latter inequalities yield a remarkable spectral inequality and the null controllability of the heat equation. Secondly, we prove Carleman estimates for parabolic operators. We state them locally in space at first, and patch them together to obtain a global estimate. This second approach also yields the null controllability of the heat equation.

Journal ArticleDOI
TL;DR: In this paper, the improved fractional sub-equation method is proposed to construct analytical solutions of nonlinear evolution equations involving Jumarie's modified Riemann-Liouville derivative.

Journal ArticleDOI
TL;DR: In this article, the effects of magnetic field and nanoparticle on the Jeffery-Hamel flow are studied using a powerful analytical method called the Adomian decomposition method (ADM), which reduces the traditional Navier-Stokes equation of fluid mechanics and Maxwell's electromagnetism governing equations to nonlinear ordinary differential equations to model the problem.
Abstract: In this study, the effects of magnetic field and nanoparticle on the Jeffery-Hamel flow are studied using a powerful analytical method called the Adomian decomposition method (ADM). The traditional Navier-Stokes equation of fluid mechanics and Maxwell’s electromagnetism governing equations are reduced to nonlinear ordinary differential equations to model the problem. The obtained results are well agreed with that of the Runge-Kutta method. The present plots confirm that the method has high accuracy for different α, Ha, and Re numbers. The flow field inside the divergent channel is studied for various values of Hartmann number and angle of channel. The effect of nanoparticle volume fraction in the absence of magnetic field is investigated.

Posted Content
TL;DR: The Unified Form Language is presented, which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation and generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations.
Abstract: We present the Unified Form Language (UFL), which is a domain-specific language for representing weak formulations of partial differential equations with a view to numerical approximation. Features of UFL include support for variational forms and functionals, automatic differentiation of forms and expressions, arbitrary function space hierarchies for multi-field problems, general differential operators and flexible tensor algebra. With these features, UFL has been used to effortlessly express finite element methods for complex systems of partial differential equations in near-mathematical notation, resulting in compact, intuitive and readable programs. We present in this work the language and its construction. An implementation of UFL is freely available as an open-source software library. The library generates abstract syntax tree representations of variational problems, which are used by other software libraries to generate concrete low-level implementations. Some application examples are presented and libraries that support UFL are highlighted.

Journal ArticleDOI
TL;DR: A thermodynamically consistent four-species model of tumor growth on the basis of the continuum theory of mixtures, unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction-diffusion equation.
Abstract: In this paper, we develop a thermodynamically consistent four-species model of tumor growth on the basis of the continuum theory of mixtures. Unique to this model is the incorporation of nutrient within the mixture as opposed to being modeled with an auxiliary reaction-diffusion equation. The formulation involves systems of highly nonlinear partial differential equations of surface effects through diffuse-interface models. A mixed finite element spatial discretization is developed and implemented to provide numerical results demonstrating the range of solutions this model can produce. A time-stepping algorithm is then presented for this system, which is shown to be first order accurate and energy gradient stable. The results of an array of numerical experiments are presented, which demonstrate a wide range of solutions produced by various choices of model parameters.

Journal ArticleDOI
TL;DR: This work studies the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation, included in the energy equation, and variable wall temperature.
Abstract: In this work, we study the flow and heat transfer characteristics of a viscous nanofluid over a nonlinearly stretching sheet in the presence of thermal radiation, included in the energy equation, and variable wall temperature. A similarity transformation was used to transform the governing partial differential equations to a system of nonlinear ordinary differential equations. An efficient numerical shooting technique with a fourth-order Runge-Kutta scheme was used to obtain the solution of the boundary value problem. The variations of dimensionless surface temperature, as well as flow and heat-transfer characteristics with the governing dimensionless parameters of the problem, which include the nanoparticle volume fraction ϕ, the nonlinearly stretching sheet parameter n, the thermal radiation parameter NR, and the viscous dissipation parameter Ec, were graphed and tabulated. Excellent validation of the present numerical results has been achieved with the earlier nonlinearly stretching sheet problem of Cortell for local Nusselt number without taking the effect of nanoparticles.

Journal ArticleDOI
TL;DR: It is proved that a propagation of chaos phenomenon takes place, namely that in the mean-field limit, any finite number of neurons become independent and, within each population, have the same probability distribution.
Abstract: We derive the mean-field equations arising as the limit of a network of interacting spiking neurons, as the number of neurons goes to infinity. The neurons belong to a fixed number of populations and are represented either by the Hodgkin-Huxley model or by one of its simplified version, the FitzHugh-Nagumo model. The synapses between neurons are either electrical or chemical. The network is assumed to be fully connected. The maximum conductances vary randomly. Under the condition that all neurons' initial conditions are drawn independently from the same law that depends only on the population they belong to, we prove that a propagation of chaos phenomenon takes place, namely that in the mean-field limit, any finite number of neurons become independent and, within each population, have the same probability distribution. This probability distribution is a solution of a set of implicit equations, either nonlinear stochastic differential equations resembling the McKean-Vlasov equations or non-local partial differential equations resembling the McKean-Vlasov-Fokker-Planck equations. We prove the wellposedness of the McKean-Vlasov equations, i.e. the existence and uniqueness of a solution. We also show the results of some numerical experiments that indicate that the mean-field equations are a good representation of the mean activity of a finite size network, even for modest sizes. These experiments also indicate that the McKean-Vlasov-Fokker-Planck equations may be a good way to understand the mean-field dynamics through, e.g. a bifurcation analysis.

Journal ArticleDOI
TL;DR: The paper exposes some of the potential of a recently introduced backstepping transformation for linear uncertain time-delay systems to address the classic problems of equilibrium regulation under partial measurements, disturbance rejection, parameter or delay adaptation.

Journal ArticleDOI
TL;DR: The explicit analytic interaction solutions between solitary waves and cnoidal waves are obtained through the localization procedure of nonlocal symmetries which are related to Darboux transformation for the well-known KdV equation.
Abstract: In nonlinear science, it is very difficult to find exact interaction solutions among solitons and other kinds of complicated waves such as cnoidal waves and Painlevwaves. Actually, even if for the most well-known prototypical models such as the Kortewet-de Vries (KdV) equation and the Kadomtsev-Petviashvili (KP) equation, this kind of problem has not yet been solved. In this paper, the explicit analytic interaction solutions betweensolitarywavesandcnoidalwavesareobtainedthroughthelocalizationprocedureofnonlocalsymmetries which are related to Darboux transformation for the well-known KdV equation. The same approach also yields some other types of interaction solutions among different types of solutions such as solitary waves, rational solutions, Bessel function solutions, and/or general Painlev´ e II solutions.


Journal ArticleDOI
TL;DR: A bird’s eye view on the development of numerical methods for solving partial differential equations with a particular emphasis on nonlinear PDEs is provided.
Abstract: Numerical methods were first put into use as an effective tool for solving partial differential equations (PDEs) by John von Neumann in the mid1940s. In a 1949 letter von Neumann wrote “the entire computing machine is merely one component of a greater whole, namely, of the unity formed by the computing machine, the mathematical problems that go with it, and the type of planning which is called by both.” The “greater whole” is viewed today as scientific computation: over the past sixty years, scientific computation has emerged as the most versatile tool to complement theory and experiments, and numerical methods for solving PDEs are at the heart of many of today’s advanced scientific computations. Numerical solutions found their way from financial models on Wall Street to traffic models on Main Street. Here we provide a bird’s eye view on the development of these numerical methods with a particular emphasis on nonlinear PDEs.

Book
10 Jan 2012
TL;DR: In this article, the authors extend the major finite-dimensional state-space results to a large class of distributed parameter systems, which contain models for delay systems, as well as partial differential equations, and allow for unbounded inputs and outputs.
Abstract: The aim of this book is to extend the major finite-dimensional state-space results to a large class of distributed parameter systems. These distributed parameter systems contain models for delay systems, as well as partial differential equations, and allow for unbounded inputs and outputs.

Journal ArticleDOI
TL;DR: The steady boundary layer flow of nanofluid over an exponential stretching surface is investigated analytically and expressions for velocity, temperature and nanoparticle volume fraction are computed for some values of the parameters.
Abstract: The steady boundary layer flow of nanofluid over an exponential stretching surface is investigated analytically. The transport equations include the effects of Brownian motion parameter and thermophoresis parameter. The highly nonlinear coupled partial differential equations are simplified with the help of suitable similarity transformations. The reduced equations are then solved analytically with the help of homotopy analysis method (HAM). The convergence of HAM solutions are obtained by plotting h-curve. The expressions for velocity, temperature and nanoparticle volume fraction are computed for some values of the parameters namely, suction injection parameter α, Lewis number Le, the Brownian motion parameter Nb and thermophoresis parameter Nt.

Book
15 Feb 2012
TL;DR: In this paper, Tataru et al. showed that the wave map evolution of any smooth data exists globally as a smooth function, and that for large data, singularities may occur in finite time for $M=S^2$ as target.
Abstract: Wave maps are the simplest wave equations taking their values in a Riemannian manifold (M,g). Their Lagrangian is the same as for the scalar equation, the only difference being that lengths are measured with respect to the metric g. By Noether's theorem, symmetries of the Lagrangian imply conservation laws for wave maps, such as conservation of energy. In coordinates, wave maps are given by a system of semilinear wave equations. Over the past 20 years important methods have emerged which address the problem of local and global wellposedness of this system. Due to weak dispersive effects, wave maps defined on Minkowski spaces of low dimensions, such as $R_{t,x}^{2+1}$, present particular technical difficulties. This class of wave maps has the additional important feature of being energy critical, which refers to the fact that the energy scales exactly like the equation. Around 2000 Daniel Tataru and Terence Tao, building on earlier work of Klainerman–Machedon, proved that smooth data of small energy lead to global smooth solutions for wave maps from 2+1 dimensions into target manifolds satisfying some natural conditions. In contrast, for large data, singularities may occur in finite time for $M=S^2$ as target. This monograph establishes that for $H$ as target the wave map evolution of any smooth data exists globally as a smooth function. While we restrict ourselves to the hyperbolic plane as target the implementation of the concentration-compactness method, the most challenging piece of this exposition, yields more detailed information on the solution. This monograph will be of interest to experts in nonlinear dispersive equations, in particular to those working on geometric evolution equations.

Journal ArticleDOI
TL;DR: In this article, the authors derived analytical solutions of three types of the MT-TSCR-FADE with Dirichlet boundary conditions, which are called respectively the multi-term time-fractional diffusion terms, the multiscale time fractional wave terms, and the multiiscale mixed diffusion-wave terms.

Journal ArticleDOI
TL;DR: A family of time-varying hyperbolic systems of balance laws is considered and a strictification approach is taken to obtain a time-Varying strict Lyapunov function, which allows to establish asymptotic stability in the general case and a robustness property with respect to additive disturbances of input-to-state stability (ISS) type.
Abstract: A family of time-varying hyperbolic systems of balance laws is considered. The partial differential equations of this family can be stabilized by selecting suitable boundary conditions. For the stabilized systems, the classical technique of construction of Lyapunov functions provides a function which is a weak Lyapunov function in some cases, but is not in others. We transform this function through a strictification approach to obtain a time-varying strict Lyapunov function. It allows us to establish asymptotic stability in the general case and a robustness property with respect to additive disturbances of input-to-state stability (ISS) type. Two examples illustrate the results.

Journal ArticleDOI
TL;DR: Stable variational formulations are proposed for certain linear, unsymmetric operators with first order transport equations in bounded domains serving as the primary focus of this paper to adaptively resolve anisotropic solution features such as propagating singularities.
Abstract: We propose stable variational formulations for certain linear, unsymmetric operators with first order transport equations in bounded domains serving as the primary focus of this paper. The central objective is to develop for such classes adaptive solution concepts with provable error reduction. To adaptively resolve anisotropic solution features such as propagating singularities, the presently proposed variational formulations allow, in particular, the employment of trial spaces spanned by directional representation systems. Since such systems, typically given as frames, are known to be stable only in $L_2$, special emphasis is placed on $L_2$-stable formulations. The proposed stability concept is based on perturbations of certain “ideal” test spaces in Petrov--Galerkin formulations; see also [L. F. Demkowicz and J. Gopalakrishnan, Comput. Methods Appl. Mech. Engrg., 199 (2010), pp. 1558--1572], [L. Demkowicz and J. Gopalakrishnan, Numer. Methods Partial Differential Equations, 27 (2011), pp. 70--105], [J...

Journal ArticleDOI
TL;DR: In this paper, a unified gas-kinetic scheme based on the Shakhov model in two-dimensional space is presented, which can capture non-equilibrium flow physics in the transition and rarefied flow regimes.
Abstract: With discretized particle velocity space, a multi-scale unified gas-kinetic scheme for entire Knudsen number flows has been constructed based on the kinetic model in one-dimensional case (J. Comput. Phys., vol. 229 (2010), pp. 7747-7764). For the kinetic equation, to extend a one-dimensional scheme to multidimensional flow is not so straightforward. The major factor is that addition of one dimension in physical space causes the distribution function to become two-dimensional, rather than axially symmetric, in velocity space. In this paper, a unified gas-kinetic scheme based on the Shakhov model in two-dimensional space will be presented. Instead of particle-based modeling for the rarefied flow, such as the direct simulation Monte Carlo (DSMC) method, the philosophical principal underlying the current study is a partial-different- ial-equation (PDE)-based modeling. Since the valid scale of the kinetic equation and the scale of mesh size and time step may be significantly different, the gas evolu- tion in a discretized space is modeled with the help of kinetic equation, instead of directly solving the partial differential equation. Due to the use of both hydrody- namic and kinetic scales flow physics in a gas evolution model at the cell interface, the unified scheme can basically present accurate solution in all flow regimes from the free molecule to the Navier-Stokes solutions. In comparison with the DSMC and Navier-Stokes flow solvers, the current method is much more efficient than DSMC in low speed transition and continuum flow regimes, and it has better capability than NS solver in capturing of non-equilibrium flow physics in the transition and rarefied flow regimes. As a result, the current method can be useful in the flow simulation where both continuum and rarefied flow physics needs to be resolved in a single com- putation. This paper will extensively evaluate the performance of the unified scheme from free molecule to continuum NS solutions, and from low speed micro-flow to high speed non-equilibrium aerodynamics. The test cases clearly demonstrate that the uni- fied scheme is a reliable method for the rarefied flow computations, and the scheme provides an important tool in the study of non-equilibrium flow.