Comparison of the ENATE approach and discontinuous Galerkin spectral element method in 1D nonlinear transport equations

In this paper a comparison of the performance of two ways of discretizing the nonlinear convection-diffusion equation in a one-dimensional (1D) domain is performed. The two approaches can be considered within the class of high-order methods. The first one is the discontinuous Galerkin method, which has been profusely used to solve general transport equations, either coupled as the Navier-Stokes equations, or on their own. On the other hand, the ENATE procedure (Enhanced Numerical Approximation of a Transport Equation), uses the exact solution to obtain an exact algebraic equation with integral coefficients that link nodal values with a three-point stencil. This paper is the first of thorough assessments of ENATE by comparing it with well-established high-order methods. Several test cases of the steady Burgers' equation with and without source have been chosen for comparison

Wavelet Methods for the Solutions of Partial and Fractional Differential Equations Arising in Physical Problems

The subject of fractional calculus has gained considerable popularity and importance during the past three decades or so, mainly due to its demonstrated applications in numerous seemingly diverse and widespread fields of science and engineering. It deals with derivatives and integrals of arbitrary orders. The fractional derivative has been occurring in many physical problems, such as frequency-dependent damping behavior of materials, motion of a large thin plate in a Newtonian fluid, creep and relaxation functions for viscoelastic materials, the PI D controller for the control of dynamical systems etc. Phenomena in electromagnetics, acoustics, viscoelasticity, electrochemistry, control theory, neutron point kinetic model, anomalous diffusion, Brownian motion, signal and image processing, fluid dynamics and material science are well described by differential equations of fractional order. Generally, nonlinear partial differential equations of fractional order are difficult to solve. So for the last few decades, a great deal of attention has been directed towards the solution (both exact and numerical) of these problems. The aim of this dissertation is to present an extensive study of different wavelet methods for obtaining numerical solutions of mathematical problems occurring in disciplines of science and engineering. This present work also provides a comprehensive foundation of different wavelet methods comprising Haar wavelet method, Legendre wavelet method, Legendre multi-wavelet methods, Chebyshev wavelet method, Hermite wavelet method and Petrov-Galerkin method. The intension is to examine the accuracy of various wavelet methods and their efficiency for solving nonlinear fractional differential equations. With the widespread applications of wavelet methods for solving difficult problems in diverse fields of science and engineering such as wave propagation, data compression, image processing, pattern recognition, computer graphics and in medical technology, these methods have been implemented to develop accurate and fast algorithms for solving integral, differential and integro-differential equations, especially those whose solutions are highly localized in position and scale. The main feature of wavelets is its ability to convert the given differential and integral equations to a system of linear or nonlinear algebraic equations, which can be solved by numerical methods. Therefore, our main focus in the present work is to analyze the application of wavelet based transform methods for solving the problem of fractional order partial differential equations. The introductory concept of wavelet, wavelet transform and multi-resolution analysis (MRA) have been discussed in the preliminary chapter. The basic idea of various analytical and numerical methods viz. Variational Iteration Method (VIM), Homotopy Perturbation Method (HPM), Homotopy Analysis Method (HAM), First Integral Method (FIM), Optimal Homotopy Asymptotic Method (OHAM), Haar Wavelet Method, Legendre Wavelet Method, Chebyshev Wavelet Method and Hermite Wavelet Method have been presented in chapter 1. In chapter 2, we have considered both analytical and numerical approach for solving some particular nonlinear partial differential equations like Burgers’ equation, modified Burgers’ equation, Huxley equation, Burgers-Huxley equation and modified KdV equation, which have a wide variety of applications in physical models. Variational Iteration Method and Haar wavelet Method are applied to obtain the analytical and numerical approximate solution of Huxley and Burgers-Huxley equations. Comparisons between analytical solution and numerical solution have been cited in tables and also graphically. The Haar wavelet method has also been applied to solve Burgers’, modified Burgers’, and modified KdV equations numerically. The results thus obtained are compared with exact solutions as well as solutions available in open literature. Error of collocation method has been presented in this chapter. Methods like Homotopy Perturbation Method (HPM) and Optimal Homotopy Asymptotic Method (OHAM) are very powerful and efficient techniques for solving nonlinear PDEs. Using these methods, many functional equations such as ordinary, partial differential equations and integral equations have been solved. We have implemented HPM and OHAM in chapter 3, in order to obtain the analytical approximate solutions of system of nonlinear partial differential equation viz. the Boussinesq-Burgers’ equations. Also, the Haar wavelet method has been applied to obtain the numerical solution of BoussinesqBurgers’ equations. Also, the convergence of HPM and OHAM has been discussed in this chapter. The mathematical modeling and simulation of systems and processes, based on the description of their properties in terms of fractional derivatives, naturally leads to differential equations of fractional order and the necessity to solve such equations. The mathematical preliminaries of fractional calculus, definitions and theorems have been presented in chapter 4. Next, in this chapter, the Haar wavelet method has been analyzed for solving fractional differential equations. The time-fractional Burgers-Fisher, generalized Fisher type equations, nonlinear time- and space-fractional Fokker-Planck equations have been solved by using two-dimensional Haar wavelet method. The obtained results are compared with the Optimal Homotopy Asymptotic Method (OHAM), the exact solutions and the results available in open literature. Comparison of obtained results with OHAM, Adomian Decomposition Method (ADM), VIM and Operational Tau Method (OTM) has been demonstrated in order to justify the accuracy and efficiency of the proposed schemes. The convergence of two-dimensional Haar wavelet technique has been provided at the end of this chapter. In chapter 5, the fractional differential equations such as KdV-Burger-Kuramoto (KBK) equation, seventh order KdV (sKdV) equation and Kaup-Kupershmidt (KK) equation have been solved by using two-dimensional Legendre wavelet and Legendre multi-wavelet methods. The main focus of this chapter is the application of two-dimensional Legendre wavelet technique for solving nonlinear fractional differential equations like timefractional KBK equation, time-fractional sKdV equation in order to demonstrate the efficiency and accuracy of the proposed wavelet method. Similarly in chapter 6, twodimensional Chebyshev wavelet method has been implemented to obtain the numerical solutions of the time-fractional Sawada-Kotera equation, fractional order Camassa-Holm equation and Riesz space-fractional sine-Gordon equations. The convergence analysis has been done for these wavelet methods. In chapter 7, the solitary wave solution of fractional modified Fornberg-Whitham equation has been attained by using first integral method and also the approximate solutions obtained by optimal homotopy asymptotic method (OHAM) are compared with the exact solutions acquired by first integral method. Also, the Hermite wavelet method has been implemented to obtain approximate solutions of fractional modified Fornberg-Whitham equation. The Hermite wavelet method is implemented to system of nonlinear fractional differential equations viz. the fractional Jaulent-Miodek equations. Convergence of this wavelet methods has been discussed in this chapter. Chapter 8 emphasizes on the application of Petrov-Galerkin method for solving the fractional differential equations such as the fractional KdV-Burgers’ (KdVB) equation and the fractional Sharma-TassoOlver equation with a view to exhibit the capabilities of this method in handling nonlinear equation. The main objective of this chapter is to establish the efficiency and accuracy of Petrov-Galerkin method in solving fractional differential equtaions numerically by implementing a linear hat function as the trial function and a quintic B-spline function as the test function. Various wavelet methods have been successfully employed to numerous partial and fractional differential equations in order to demonstrate the validity and accuracy of these procedures. Analyzing the numerical results, it can be concluded that the wavelet methods provide worthy numerical solutions for both classical and fractional order partial differential equations. Finally, it is worthwhile to mention that the proposed wavelet methods are promising and powerful methods for solving fractional differential equations in mathematical physics. This work also aimed at, to make this subject popular and acceptable to engineering and science community to appreciate the universe of wonderful mathematics, which is in between classical integer order differentiation and integration, which till now is not much acknowledged, and is hidden from scientists and engineers. Therefore, our goal is to encourage the reader to appreciate the beauty as well as the usefulness of these numerical wavelet based techniques in the study of nonlinear physical system

Instanton based importance sampling for rare events in stochastic PDEs

We present a new method for sampling rare and large fluctuations in a non-equilibrium system governed by a stochastic partial differential equation (SPDE) with additive forcing. To this end, we deploy the so-called instanton formalism that corresponds to a saddle-point approximation of the action in the path integral formulation of the underlying SPDE. The crucial step in our approach is the formulation of an alternative SPDE that incorporates knowledge of the instanton solution such that we are able to constrain the dynamical evolutions around extreme flow configurations only. Finally, a reweighting procedure based on the Girsanov theorem is applied to recover the full distribution function of the original system. The entire procedure is demonstrated on the example of the one-dimensional Burgers equation. Furthermore, we compare our method to conventional direct numerical simulations as well as to Hybrid Monte Carlo methods. It will be shown that the instanton-based sampling method outperforms both approaches and allows for an accurate quantification of the whole probability density function of velocity gradients from the core to the very far tails.Comment: 8 pages, 4 figure

Error Estimations, Error Computations, and Convergence Rates in FEM for BVPs

This paper presents derivation of a priori error estimates and convergence rates of finite element processes for boundary value problems (BVPs) described by self adjoint, non-self adjoint, and nonlinear differential operators. A posteriori error estimates are discussed in context with local approximations in higher order scalar product spaces. A posteriori error computational framework (without the knowledge of theoretical solution) is presented for all BVPs regardless of the method of approximation employed in constructing the integral form. This enables computations of local errors as well as the global errors in the computed finite element solutions. The two most significant and essential aspects of the research presented in this paper that enable all of the features described above are: 1) ensuring variational consistency of the integral form(s) resulting from the methods of approximation for self adjoint, non-self adjoint, and nonlinear differential operators and 2) choosing local approximations for the elements of a discretization in a subspace of a higher order scalar product space that is minimally conforming, hence ensuring desired global differentiability of the approximations over the discretizations. It is shown that when the theoretical solution of a BVP is analytic, the a priori error estimate (in the asymptotic range, discussed in a later section of the paper) is independent of the method of approximation or the nature of the differential operator provided the resulting integral form is variationally consistent. Thus, the finite element processes utilizing integral forms based on different methods of approximation but resulting in VC integral forms result in the same a priori error estimate and convergence rate. It is shown that a variationally consistent (VC) integral form has best approximation property in some norm, conversely an integral form with best approximation property in some norm is variationally consistent. That is best approximation property of the integral form and the VC of the integral form is equivalent, one cannot exist without the other, hence can be used interchangeably. Dimensional model problems consisting of diffusion equation, convection-diffusion equation, and Burgers equation described by self adjoint, non-self adjoint, and nonlinear differential operators are considered to present extensive numerical studies using Galerkin method with weak form (GM/WF) and least squares process (LSP) to determine computed convergence rates of various error norms and present comparisons with the theoretical convergence rates

The instanton method and its numerical implementation in fluid mechanics

A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional fluid dynamical problems. We illustrate these ideas using the two-dimensional Burgers equation and the three-dimensional Navier-Stokes equations

Wave Structures and Nonlinear Balances in a Family of 1+1 Evolutionary PDEs

We study the following family of evolutionary 1+1 PDEs that describe the balance between convection and stretching for small viscosity in the dynamics of 1D nonlinear waves in fluids: m_t + \underbrace{um_x \} _{(-2mm)\hbox{convection}(-2mm)} + \underbrace{b u_xm \} _{(-2mm)\hbox{stretching}(-2mm)} = \underbrace{\nu m_{xx}\ }_{(-2mm)\hbox{viscosity}}, \quad\hbox{with}\quad u=g*m . Here $u=g*m$ denotes $u(x)=\int_{-\infty}^\infty g(x-y)m(y) dy .$ We study exchanges of stability in the dynamics of solitons, peakons, ramps/cliffs, leftons, stationary solutions and other solitary wave solutions associated with this equation under changes in the nonlinear balance parameter $b$.Comment: 69 pages, 26 figure

Gr\"obner Bases and Generation of Difference Schemes for Partial Differential Equations

In this paper we present an algorithmic approach to the generation of fully conservative difference schemes for linear partial differential equations. The approach is based on enlargement of the equations in their integral conservation law form by extra integral relations between unknown functions and their derivatives, and on discretization of the obtained system. The structure of the discrete system depends on numerical approximation methods for the integrals occurring in the enlarged system. As a result of the discretization, a system of linear polynomial difference equations is derived for the unknown functions and their partial derivatives. A difference scheme is constructed by elimination of all the partial derivatives. The elimination can be achieved by selecting a proper elimination ranking and by computing a Gr\"obner basis of the linear difference ideal generated by the polynomials in the discrete system. For these purposes we use the difference form of Janet-like Gr\"obner bases and their implementation in Maple. As illustration of the described methods and algorithms, we construct a number of difference schemes for Burgers and Falkowich-Karman equations and discuss their numerical properties.Comment: Published in SIGMA (Symmetry, Integrability and Geometry: Methods and Applications) at http://www.emis.de/journals/SIGMA

The probability density function tail of the Kardar-Parisi-Zhang equation in the strongly non-linear regime

An analytical derivation of the probability density function (PDF) tail describing the strongly correlated interface growth governed by the nonlinear Kardar-Parisi-Zhang equation is provided. The PDF tail exactly coincides with a Tracy-Widom distribution i.e. a PDF tail proportional to $\exp( - c w_2^{3/2})$, where $w_2$ is the the width of the interface. The PDF tail is computed by the instanton method in the strongly non-linear regime within the Martin-Siggia-Rose framework using a careful treatment of the non-linear interactions. In addition, the effect of spatial dimensions on the PDF tail scaling is discussed. This gives a novel approach to understand the rightmost PDF tail of the interface width distribution and the analysis suggests that there is no upper critical dimension.Comment: 17 pages, 2 figure

Sharp bounds on enstrophy growth in the viscous Burgers equation

We use the Cole--Hopf transformation and the Laplace method for the heat equation to justify the numerical results on enstrophy growth in the viscous Burgers equation on the unit circle. We show that the maximum enstrophy achieved in the time evolution is scaled as $\mathcal{E}^{3/2}$, where $\mathcal{E}$ is the large initial enstrophy, whereas the time needed for reaching the maximal enstrophy is scaled as $\mathcal{E}^{-1/2}$. These bounds are sharp for sufficiently smooth initial conditions.Comment: 12 page